Update on CO2 and temperature profiles retrievals from NOMAD-SO on board ExoMars TGO

2021 ◽  
Author(s):  
Loïc Trompet ◽  
Ann Carine Vandaele ◽  
Shohei Aoki ◽  
Justin Erwin ◽  
Ian Thomas ◽  
...  
Keyword(s):  
Local Density ◽  
Regularization Parameter ◽  
Estimation Method ◽  
Diffraction Order ◽  
Tunable Filter ◽  
Fine Tuning ◽  
Temperature Profiles ◽  
Transfer Program ◽  
Solar Occultation ◽  
Retrieval Scheme

<ul> <li>The SO channel of the NOMAD instrument</li> </ul> <p>The NOMAD-SO channel [1] is an infrared spectrometer working in the 2.2 to 4.3 µm spectral range (2200-4500 cm<sup>-1</sup>) and started to perform solar occultation measurement on April 21, 2018. The instrument is composed of an echelle grating coupled to an Acousto-Optical Tunable Filter for the diffraction order selection. As TGO is on a quasi-circular orbit at around 400 km of altitude, it performs one orbit every two hours. During a solar occultation measurement, SO scans six diffraction orders each second. These diffraction orders are recorded on four bins leading to a vertical sampling lower than one km. The calibration of the SO channel is described in [2] and is still being fine-tuned.</p> <ul> <li>CO<sub>2</sub> density and temperature profiles retrievals</li> </ul> <p>Several diffraction orders probe different altitude ranges as they contain CO<sub>2</sub> lines with different intensities that appear and saturate at different altitudes. Correct temperature profiles are necessary for the retrieval of several species and the profiles have to be carefully retrieved as their inversion is very sensitive to noise. We use the following retrieval scheme:</p> <p>For each solar occultation measurement, we derive a slant column profile of CO<sub>2</sub> using ASIMUT-ALVL [3]. ASIMUT is a radiative transfer program developed at BIRA-IASB and uses the Optimal Estimation Method for regularization [4]. The GEM-Mars GCM provides the <em>a priori</em> profiles of CO<sub>2</sub> local density, pressure and temperature. We then apply an iterated Tikhonov regularization to derive a regularized local density profile using an improved version of the algorithm described in [5]. This method requires the selection of a regularization parameter to reduce as much as possible the presence of noise in the profile while keeping the real variations. This retrieval scheme allows a fine-tuning of the reguralization parameter. We finally apply the hydrostatic equilibrium equation to derive the temperature profiles [6]. We derived the NOMAD-SO CO<sub>2</sub> and temperature profiles for MY34 and 35.</p>

Ozone vertical profiles from TGO/NOMAD-UVIS: an inter-comparison of three retrieval schemes

2021 ◽  
Author(s):  
Arianna Piccialli ◽  
Ann Carine Vandaele ◽  
Yannick Willame ◽  
Shohei Aoki ◽  
Cedric Depiesse ◽  
...  
Keyword(s):  
Local Density ◽  
Estimation Method ◽  
Optimal Estimation ◽  
Space Science ◽  
Solar Occultation ◽  
Retrieval Scheme ◽  
Local Densities ◽  
The Vertical Distribution ◽  
Applied Optics ◽  
Infrared Channel

<p>We will present the vertical distribution of <strong>ozone</strong> obtained from <strong>NOMAD-UVIS solar occultations</strong> and we will compare the results of three retrieval schemes.</p><p><strong>NOMAD</strong> (Nadir and Occultation for MArs Discovery) is a spectrometer composed of 3 channels: 1) a solar occultation channel (SO) operating in the infrared (2.3-4.3 μm); 2) a second infrared channel LNO (2.3-3.8 μm) capable of doing nadir, as well as solar occultation and limb; and 3) an ultraviolet/visible channel <strong>UVIS</strong> (200-650 nm) that can work in the three observation modes [1,2].</p><p>The UVIS channel has a spectral resolution <1.5 nm. In the solar occultation mode it is mainly devoted to study the climatology of <strong>ozone</strong> and <strong>aerosols</strong> content [3].</p><p>Since the beginning of operations, on 21 April 2018, NOMAD UVIS acquired more than 4000 solar occultations with an almost complete coverage of the planet.</p><p>NOMAD-UVIS spectra are simulated using three different retrieval scheme:</p><p>1) An onion peeling approach based on [4,5] deriving slant columns at the different altitudes sounded, from which local densities are obtained;</p><p>2) The line-by-line radiative transfer code ASIMUT-ALVL developed at IASB-BIRA [6] using the Optimal Estimation Method to derive the local density profile in one go (on all transmittances of one occultation observation);</p><p>3) A direct onion peeling method deriving sequentially from top to bottom the local densities in the different layers probed.</p><p>We will compare results obtained from the different retrieval methods as well as their uncertainties and we will discuss the advantages and difficulties of each method.</p><p><strong>References</strong></p><p>[1] Vandaele, A.C., et al., Planetary and Space Science, Vol. 119, pp. 233–249, 2015.</p><p>[2] Neefs, E., et al., Applied Optics, Vol. 54 (28), pp. 8494-8520, 2015.</p><p>[3] M.R. Patel et al., In: Appl. Opt. 56.10 (2017), pp. 2771–2782. DOI: 10.1364/AO.56.002771.</p><p>[4] Quémerais, E.,et al. J.Geophys. Res. (Planets)111, 9, 2006.</p><p>[5] Piccialli, A. et al., Planetary and Space Science, 113-114(2015) 321–335</p><p>[6] Vandaele, A.C., et al., JGR, 2008. 113 doi:10.1029/2008JE003140.</p>

scholarly journals Jewel: A Novel Method for Joint Estimation of Gaussian Graphical Models

Mathematics ◽  
2021 ◽  
Vol 9 (17) ◽  
pp. 2105
Author(s):  
Claudia Angelini ◽  
Daniela De De Canditiis ◽  
Anna Plaksienko
Keyword(s):  
Graphical Models ◽  
Regulatory Networks ◽  
Regularization Parameter ◽  
Estimation Method ◽  
Real Data ◽  
Penalized Regression ◽  
Joint Estimation ◽  
Gaussian Graphical Models ◽  
Lasso Penalty ◽  
Independence Structure

In this paper, we consider the problem of estimating multiple Gaussian Graphical Models from high-dimensional datasets. We assume that these datasets are sampled from different distributions with the same conditional independence structure, but not the same precision matrix. We propose jewel, a joint data estimation method that uses a node-wise penalized regression approach. In particular, jewel uses a group Lasso penalty to simultaneously guarantee the resulting adjacency matrix’s symmetry and the graphs’ joint learning. We solve the minimization problem using the group descend algorithm and propose two procedures for estimating the regularization parameter. Furthermore, we establish the estimator’s consistency property. Finally, we illustrate our estimator’s performance through simulated and real data examples on gene regulatory networks.

scholarly journals Solar occultation with SCIAMACHY: algorithm description and first validation

2005 ◽  
Vol 5 (1) ◽  
pp. 17-66 ◽  
Author(s):  
J. Meyer ◽  
A. Bracher ◽  
A. Rozanov ◽  
A. C. Schlesier ◽  
H. Bovensmann ◽  
...  
Keyword(s):  
Vertical Profile ◽  
Estimation Method ◽  
Optimal Estimation ◽  
Trace Gas ◽  
Transmission Spectra ◽  
Atmospheric Transmission ◽  
Gas Concentration ◽  
Solar Occultation ◽  
First Results ◽  
The Vertical Distribution

Abstract. This presentation concentrates on solar occultation measurements with the spaceborne spectrometer SCIAMACHY in the UV-Vis wavelength range. Solar occultation measurements provide unique information about the vertical distribution of atmospheric constituents. For retrieval of vertical trace gas concentration profiles, an algorithm has been developed based on the optimal estimation method. The forward model is capable to simulate the extinction signals of different species as they occur in atmospheric transmission spectra obtained from occultation measurements. Furthermore, correction algorithms have been implemented to address shortcomings of the tangent height pre-processing and inhomogeneities of measured solar spectra. First results of O3 and NO2 vertical profile retrievals have been validated with data from ozone sondes and satellite based occultation instruments. The validation shows very promising results for SCIAMACHY O3 and NO2 values between 15 to 35 km with errors in the order of 10% and 15%, respectively.

scholarly journals A Novel Improved Maximum Entropy Regularization Technique and Application to Identification of Dynamic Loads on the Coal Rock

2019 ◽  
Vol 2019 ◽  
pp. 1-12
Author(s):  
Chunsheng Liu ◽  
Chunping Ren
Keyword(s):  
Maximum Entropy ◽  
Dynamic Load ◽  
Regularization Parameter ◽  
Estimation Method ◽  
Identification Problem ◽  
Iteration Algorithm ◽  
Load Identification ◽  
Signal Processing Algorithm ◽  
Rock Structure ◽  
Coal Rock

A new signal processing algorithm was proposed to identify the dynamic load acting on the coal-rock structure. First, the identification model for dynamic load is established through the relationship between the uncertain load vector, and the assembly matrix of the responses was measured by the machinery dynamic system. Then, the entropy item of maximum entropy regularization (MER) is redesigned using the robust estimation method, and the elongated penalty function according to the ill-posedness characteristics of load identification, which was named as a novel improved maximum entropy regularization (IMER) technique, was proposed to process the dynamic load signals. Finally, the load identification problem is transformed into an unconstrained optimization problem and an improved Newton iteration algorithm was proposed to solve the objective function. The result of IMER technique is compared with MER technique, and it is found that IMER technique is available for analyzing the dynamic load signals due to higher signal-noise ratio, lower restoration time, and fewer iterative steps. Experiments were performed to investigate the effect on the performance of dynamic load signals identification by different regularization parameters and calculation parameters, pi, respectively. Experimental results show that the identified dynamic load signals are closed to the actual load signals using IMER technique combined with the proposed PSO-L regularization parameter selection method. Selecting optimal calculated parameters pi is helpful to overcome the ill-condition of dynamic load signals identification and to obtain the stable and approximate solutions of inverse problems in practical engineering. Meanwhile, the proposed IMER technique can also play a guiding role for the coal-rock interface identification.

scholarly journals Improvements to a long-term Rayleigh-scatter lidar temperature climatology by using an optimal estimation method

2018 ◽  
Vol 11 (11) ◽  
pp. 6043-6058 ◽  
Author(s):  
Ali Jalali ◽  
Robert J. Sica ◽  
Alexander Haefele
Keyword(s):  
Traditional Method ◽  
Middle Atmosphere ◽  
Density Measurement ◽  
Estimation Method ◽  
Optimal Estimation ◽  
Temperature Profiles ◽  
Vertical Resolution ◽  
Order Of Magnitude ◽  
Lidar Equation ◽  
Model Temperature

Abstract. Hauchecorne and Chanin (1980) developed a robust method to calculate middle-atmosphere temperature profiles using measurements from Rayleigh-scatter lidars. This traditional method has been successfully used to greatly improve our understanding of middle-atmospheric dynamics, but the method has some shortcomings regarding the calculation of systematic uncertainties and the vertical resolution of the retrieval. Sica and Haefele (2015) have shown that the optimal estimation method (OEM) addresses these shortcomings and allows temperatures to be retrieved with confidence over a greater range of heights than the traditional method. We have calculated a temperature climatology from 519 nights of Purple Crow Lidar Rayleigh-scatter measurements using an OEM. Our OEM retrieval is a first-principle retrieval in which the forward model is the lidar equation and the measurements are the level-0 count returns. It includes a quantitative determination of the top altitude of the retrieved temperature profiles, the evaluation of nine systematic plus random uncertainties, and the vertical resolution of the retrieval on a profile-by-profile basis. Our OEM retrieval allows for the vertical resolution to vary with height, extending the retrieval in altitude 5 to 10 km higher than the traditional method. It also allows the comparison of the traditional method's sensitivity to two in-principle equivalent methods of specifying the seed pressure: using a model pressure seed versus using a model temperature combined with the lidar's density measurement to calculate the seed pressure. We found that the seed pressure method is superior to using a model temperature combined with the lidar-derived density. The increased altitude capability of our OEM retrievals allows for a comparison of the Rayleigh-scatter lidar temperatures throughout the entire altitude range of the sodium lidar temperature measurements. Our OEM-derived Rayleigh temperatures are shown to have improved agreement relative to our previous comparisons using the traditional method, and the agreement of the OEM-derived temperatures is the same as the agreement between existing sodium lidar temperature climatologies. This detailed study of the calculation of the new Purple Crow Lidar temperature climatology using the OEM establishes that it is both highly advantageous and practical to reprocess existing Rayleigh-scatter lidar measurements that cover long time periods, during which time the lidar may have undergone several significant equipment upgrades, while gaining an upper limit to useful temperature retrievals equivalent to an order of magnitude increase in power-aperture product due to the use of an OEM.

Removal of the Solar Component in AVHRR 3.7-µm Radiances for the Retrieval of Cirrus Cloud Parameters

1995 ◽  
Vol 34 (2) ◽  
pp. 482-499 ◽  
Author(s):  
N. X. Rao ◽  
S. C. Ou ◽  
K. N. Liou
Keyword(s):  
Lidar Data ◽  
Cirrus Cloud ◽  
Transfer Program ◽  
Ice Crystal ◽  
Threshold Values ◽  
Cloud Parameters ◽  
Look Up Table ◽  
Retrieval Scheme ◽  
The Look ◽  
Very High

Abstract A numerical scheme has been developed to remove the solar component in the Advanced Very High Resolution Radiometer (AVHRR) 3.7-µm channel for the retrieval of cirrus parameters during daytime. This method uses a number of prescribed threshold values for AVHRR channels 1 (0.63 µm), 2 (0.8 µm), 3 (3.7 µm), 4 (10.9 µm), and 5 (12 µm) to separate clear and cloudy pixels. A look-up table relating channels 1 and 3 solar reflectances is subsequently constructed based on the prescribed mean effective ice crystal sizes and satellite geometric parameters. An adding–doubling radiative transfer program has been used to generate numerical values in the construction of the look-up table. Removal of the channel 3 solar component is accomplished by using the look-up table and the measured channel 1 reflectance. The cloud retrieval scheme described in Ou et al. has been modified in connection with the removal program. The authors have applied the removal–retrieval scheme to the AVHRR global area coverage daytime data, collected during the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment cirrus intensive field observation (FIRE IFO) at 2100 UTC 28 October 1986 over the Wisconsin area. Distributions of the retrieved cloud heights and optical depths are comparable to those determined from Geostationary Operational Environmental Satellite visible and IR channels data reported by Minnis et al. Morwver, verifications of the retrieved cirrus temperature and height against lidar data have been carried out using results reported from three FIRE IFO nations. The retrieved cloud heights are within 0.5 km of the measured lidar values.

THE PROBLEM OF THE BAND GAP IN LDA CALCULATIONS

2007 ◽  
Vol 14 (03) ◽  
pp. 481-487 ◽  
Author(s):  
I. N. YAKOVKIN ◽  
P. A. DOWBEN
Keyword(s):  
Density Functional Theory ◽  
Band Structure ◽  
Band Gap ◽  
Local Density Approximation ◽  
Density Functional ◽  
Local Density ◽  
Fine Tuning ◽  
Density Approximation ◽  
Functional Theory ◽  
Basic Physics

In calculating band structure, the local density approximation and density functional theory are widely popular and do reproduce a lot of the basic physics. Regrettably, without some fine tuning, the local density approximation and density functional theory do not generally get the details of the experimental band structure correct, in particular the band gap in semiconductors and insulators is generally found to be too small when compared with experiment. For experimentalists using commercial packages to calculate the electronic structure of materials, some caution is indicated, as some long-standing problems exist with the local density approximation and density functional theory.

scholarly journals Validating EEG, MEG and combined MEG and EEG beamforming for an estimation of the epileptogenic zone in focal cortical dysplasia

2021 ◽  
Author(s):  
Frank Neugebauer ◽  
Marios Antonakakis ◽  
Kanjana Unnwongse ◽  
Yaroslav Parpaley ◽  
Jörg Wellmer ◽  
...  
Keyword(s):  
Regularization Parameter ◽  
Focal Cortical Dysplasia ◽  
Estimation Method ◽  
Cortical Dysplasia ◽  
Source Analysis ◽  
Seizure Freedom ◽  
Epileptogenic Zone ◽  
Scanning Method ◽  
Eeg Source Analysis ◽  
Skull Conductivity

AbstractMEG and EEG source analysis is frequently used for the presurgical evaluation of pharma-coresistant epilepsy patients. The source localization of the epileptogenic zone depends, among other aspects, on the selected inverse and forward approaches and their respective parameter choices. In this validation study, we compare for the inverse problem the standard dipole scanning method with two beamformer approaches and we investigate the influence of the covariance estimation method and the strength of regularization on the localization performance for EEG, MEG and combined EEG and MEG. For forward modeling, we investigate the difference between calibrated six-compartment and standard three-compartment head modeling. In a retrospective study of two patients with focal epilepsy due to focal cortical dysplasia type IIb and seizure-freedom following lesionectomy or radiofrequency-guided thermocoagulation, we used the distance of the localization of interictal epileptic spikes to the resection cavity resp. rediofrequency lesion as reference for good localization. We found that beamformer localization can be sensitive to the choice of the regularization parameter, which has to be individually optimized. Estimation of the covariance matrix with averaged spike data yielded more robust results across the modalities. MEG was the dominant modality and provided a good localization in one case, while it was EEG for the other. When combining the modalities, the good results of the dominant modality were mostly not spoiled by the weaker modality. For appropriate regularization parameter choices, the beamformer localized better than the standard dipole scan. Compared to the importance of an appropriate regularization, the sensitivity of the localization to the head modeling was smaller, due to similar skull conductivity modeling and the fixed source space without orientation constraint.

scholarly journals Covariance Matrix Reconstruction for Direction Finding with Nested Arrays Using Iterative Reweighted Nuclear Norm Minimization

2019 ◽  
Vol 2019 ◽  
pp. 1-13 ◽  
Author(s):  
Weijie Tan ◽  
Xi’an Feng
Keyword(s):  
Covariance Matrix ◽  
Regularization Parameter ◽  
Estimation Method ◽  
Estimation Algorithm ◽  
Doa Estimation ◽  
Direction Finding ◽  
Nuclear Norm ◽  
Low Rank ◽  
Matrix Reconstruction ◽  
Linear Transform

In this paper, we address the direction finding problem in the background of unknown nonuniform noise with nested array. A novel gridless direction finding method is proposed via the low-rank covariance matrix approximation, which is based on a reweighted nuclear norm optimization. In the proposed method, we first eliminate the noise variance variable by linear transform and utilize the covariance fitting criteria to determine the regularization parameter for insuring robustness. And then we reconstruct the low-rank covariance matrix by iteratively reweighted nuclear norm optimization that imposes the nonconvex penalty. Finally, we exploit the search-free DoA estimation method to perform the parameter estimation. Numerical simulations are carried out to verify the effectiveness of the proposed method. Moreover, results indicate that the proposed method has more accurate DoA estimation in the nonuniform noise and off-grid cases compared with the state-of-the-art DoA estimation algorithm.

scholarly journals CLoNe: Automated clustering based on local density neighborhoods for application to biomolecular structural ensembles

2020 ◽  
Author(s):  
Träger Sylvain ◽  
Tamò Giorgio ◽  
Aydin Deniz ◽  
Fonti Giulia ◽  
Audagnotto Martina ◽  
...  
Keyword(s):  
Cluster Analysis ◽  
State Of The Art ◽  
Local Density ◽  
Binding Pocket ◽  
Fine Tuning ◽  
Supplementary Information ◽  
Original Algorithm ◽  
Density Peaks ◽  
Cluster Shape ◽  
Structural Ensembles

Abstract Motivation Proteins are intrinsically dynamic entities. Flexibility sampling methods, such as molecular dynamics or those arising from integrative modeling strategies are now commonplace and enable the study of molecular conformational landscapes in many contexts. Resulting structural ensembles increase in size as technological and algorithmic advancements take place, making their analysis increasingly demanding. In this regard, cluster analysis remains a go-to approach for their classification. However, many state-of-the-art algorithms are restricted to specific cluster properties. Combined with tedious parameter fine-tuning, cluster analysis of protein structural ensembles suffers from the lack of a generally applicable and easy to use clustering scheme. Results We present CLoNe, an original Python-based clustering scheme that builds on the Density Peaks algorithm of Rodriguez and Laio. CLoNe relies on a probabilistic analysis of local density distributions derived from nearest neighbors to find relevant clusters regardless of cluster shape, size, distribution and amount. We show its capabilities on many toy datasets with properties otherwise dividing state-of-the-art approaches and improves on the original algorithm in key aspects. Applied to structural ensembles, CLoNe was able to extract meaningful conformations from membrane binding events and ligand-binding pocket opening as well as identify dominant dimerization motifs or inter-domain organization. CLoNe additionally saves clusters as individual trajectories for further analysis and provides scripts for automated use with molecular visualization software. Availability www.epfl.ch/labs/lbm/resources, github.com/LBM-EPFL/CLoNe Supplementary information Supplementary data are available at Bioinformatics online.

Sign in / Sign up

Export Citation Format

Share Document