scholarly journals A Spectra Classification Methodology of Hyperspectral Infrared Images for Near Real-Time Estimation of the SO2 Emission Flux from Mount Etna with LARA Radiative Transfer Retrieval Model

2020 ◽  
Vol 12 (24) ◽  
pp. 4107
Author(s):  
Charlotte Segonne ◽  
Nathalie Huret ◽  
Sébastien Payan ◽  
Mathieu Gouhier ◽  
Valéry Catoire

Fast and accurate quantification of gas fluxes emitted by volcanoes is essential for the risk mitigation of explosive eruption, and for the fundamental understanding of shallow eruptive processes. Sulphur dioxide (SO2), in particular, is a reliable indicator to predict upcoming eruptions, and its systemic characterization allows the rapid assessment of sudden changes in eruptive dynamics. In this regard, infrared (IR) hyperspectral imaging is a promising new technology for accurately measure SO2 fluxes day and night at a frame rate down to 1 image per second. The thermal infrared region is not very sensitive to particle scattering, which is an asset for the study of volcanic plume. A ground based infrared hyperspectral imager was deployed during the IMAGETNA campaign in 2015 and provided high spectral resolution images of the Mount Etna (Sicily, Italy) plume from the North East Crater (NEC), mainly. The LongWave InfraRed (LWIR) hyperspectral imager, hereafter name Hyper-Cam, ranges between 850–1300 cm−1 (7.7–11.8 µm). The LATMOS (Laboratoire Atmosphères Milieux Observations Spatiales) Atmospheric Retrieval Algorithm (LARA), which is used to retrieve the slant column densities (SCD) of SO2, is a robust and a complete radiative transfer model, well adapted to the inversion of ground-based remote measurements. However, the calculation time to process the raw data and retrieve the infrared spectra, which is about seven days for the retrieval of one image of SO2 SCD, remains too high to infer near real-time (NRT) SO2 emission fluxes. A spectral image classification methodology based on two parameters extracting spectral features in the O3 and SO2 emission bands was developed to create a library. The relevance is evaluated in detail through tests. From data acquisition to the generation of SO2 SCD images, this method requires only ~40 s per image, which opens the possibility to infer NRT estimation of SO2 emission fluxes from IR hyperspectral imager measurements.

2021 ◽  
Author(s):  
Charlotte Segonne ◽  
Nathalie Huret ◽  
Sébastien Payan ◽  
Mathieu Gouhier

<p>Monitoring active volcanoes activity passes through the detection of fluctuations in degassing levels which may reflect changes in the magma supply rate and help inform a short-term forecast of on-going eruptions. Infrared hyperspectral imagers, which is an imaging technology still little used for volcanoes monitoring, have been deployed for various field campaigns on active volcanoes recently. For example, the Hyper-Cam LWIR (LongWave InfraRed) ranging between 850-1300 cm<sup>-1</sup> (7.7 - 11.8 µm) with a spectral resolution up to 0.25 cm<sup>-1</sup>, provided high spectral resolution images from ground-based measurements of the Mount Etna (Sicily, Italy) plume during IMAGETNA campaign in June 2015. Processing the raw data and retrieving the infrared spectra with the LATMOS (Laboratoire Atmosphères Milieux Observations Spatiales) Atmospheric Retrieval Algorithm (LARA), a robust and a complete radiative transfer model, require a calculation time of ~7 days per image.</p><p>One of the main ways of risk mitigation effects of explosive eruptions is to get a fast and accurate quantification of SO<sub>2</sub> fluxes emitted by volcanoes. In this context, using the dataset acquired during IMAGETNA campaign at Mount Etna, a spectra classification methodology has been developed to drastically decrease the calculation time and reach near real-time retrievals of SO<sub>2</sub> slant column densities. The methodology is based on a network built on two layers of information from the extraction of spectral features in the O<sub>3</sub> and SO<sub>2</sub> emission bands. A training dataset of five SO<sub>2</sub> slant column densities images retrieved with the time-consuming pixel-by-pixel retrieval method allowed the creation of a library. The spectra classification makes it possible to process each hyperspectral image in less than 40 seconds. It opens the possibility to infer near real-time estimation of SO<sub>2</sub> emission fluxes from IR hyperspectral imager measurements.</p>


2020 ◽  
Vol 13 (1) ◽  
pp. 116
Author(s):  
Lucie Leonarski ◽  
Laurent C.-Labonnote ◽  
Mathieu Compiègne ◽  
Jérôme Vidot ◽  
Anthony J. Baran ◽  
...  

The present study aims to quantify the potential of hyperspectral thermal infrared sounders such as the Infrared Atmospheric Sounding Interferometer (IASI) and the future IASI next generation (IASI-NG) for retrieving the ice cloud layer altitude and thickness together with the ice water path. We employed the radiative transfer model Radiative Transfer for TOVS (RTTOV) to simulate cloudy radiances using parameterized ice cloud optical properties. The radiances have been computed from an ice cloud profile database coming from global operational short-range forecasts at the European Center for Medium-range Weather Forecasts (ECMWF) which encloses the normal conditions, typical variability, and extremes of the atmospheric properties over one year (Eresmaa and McNally (2014)). We performed an information content analysis based on Shannon’s formalism to determine the amount and spectral distribution of the information about ice cloud properties. Based on this analysis, a retrieval algorithm has been developed and tested on the profile database. We considered the signal-to-noise ratio of each specific instrument and the non-retrieved atmospheric and surface parameter errors. This study brings evidence that the observing system provides information on the ice water path (IWP) as well as on the layer altitude and thickness with a convergence rate up to 95% and expected errors that decrease with cloud opacity until the signal saturation is reached (satisfying retrievals are achieved for clouds whose IWP is between about 1 and 300 g/m2).


Sensors ◽  
2019 ◽  
Vol 19 (7) ◽  
pp. 1532 ◽  
Author(s):  
Guido Masiello ◽  
Carmine Serio ◽  
Sara Venafra ◽  
Laurent Poutier ◽  
Frank-M. Göttsche

Timely processing of observations from multi-spectral imagers, such as SEVIRI (Spinning Enhanced Visible and Infrared Imager), largely depends on fast radiative transfer calculations. This paper mostly concerns the development and implementation of a new forward model for SEVIRI to be applied to real time processing of infrared radiances. The new radiative transfer model improves computational time by a factor of ≈7 compared to the previous versions and makes it possible to process SEVIRI data at nearly real time. The new forward model has been applied for the retrieval of surface parameters. Although the scheme can be applied for the simultaneous retrieval of temperature and emissivity, the paper mostly focuses on emissivity. The inverse scheme relies on a Kalman filter approach, which allows us to exploit a sequential processing of SEVIRI observations. Based on the new forward model, the paper also presents a validation retrieval performed with in situ observations acquired during a field experiment carried out in 2017 at Gobabeb (Namib desert) validation station. Furthermore, a comparison with IASI (Infrared Atmospheric Sounder Interferometer) emissivity retrievals has been performed as well. It has been found that the retrieved emissivities are in good agreement with each other and with in situ observations, i.e., average differences are generally well below 0.01.


2015 ◽  
Vol 12 (12) ◽  
pp. 13019-13067
Author(s):  
A. Barella-Ortiz ◽  
J. Polcher ◽  
P. de Rosnay ◽  
M. Piles ◽  
E. Gelati

Abstract. L-Band radiometry is considered to be one of the most suitable techniques to estimate surface soil moisture by means of remote sensing. Brightness temperatures are key in this process, as they are the main input in the retrieval algorithm. The work exposed compares brightness temperatures measured by the Soil Moisture and Ocean Salinity (SMOS) mission to two different sets of modelled ones, over the Iberian Peninsula from 2010 to 2012. The latter were estimated using a radiative transfer model and state variables from two land surface models: (i) ORganising Carbon and Hydrology In Dynamic EcosystEms (ORCHIDEE) and (ii) Hydrology – Tiled ECMWF Scheme for Surface Exchanges over Land (H-TESSEL). The radiative transfer model used is the Community Microwave Emission Model (CMEM). A good agreement in the temporal evolution of measured and modelled brightness temperatures is observed. However, their spatial structures are not consistent between them. An Empirical Orthogonal Function analysis of the brightness temperature's error identifies a dominant structure over the South-West of the Iberian Peninsula which evolves during the year and is maximum in Fall and Winter. Hypotheses concerning forcing induced biases and assumptions made in the radiative transfer model are analysed to explain this inconsistency, but no candidate is found to be responsible for it at the moment. Further hypotheses are proposed at the end of the paper.


2015 ◽  
Vol 15 (6) ◽  
pp. 3007-3020 ◽  
Author(s):  
R. Loughman ◽  
D. Flittner ◽  
E. Nyaku ◽  
P. K. Bhartia

Abstract. The Gauss–Seidel limb scattering (GSLS) radiative transfer (RT) model simulates the transfer of solar radiation through the atmosphere and is imbedded in the retrieval algorithm used to process data from the Ozone Mapping and Profiler Suite (OMPS) limb profiler (LP), which was launched on the Suomi NPP satellite in October 2011. A previous version of this model has been compared with several other limb scattering RT models in previous studies, including Siro, MCC++, CDIPI, LIMBTRAN, SASKTRAN, VECTOR, and McSCIA. To address deficiencies in the GSLS radiance calculations revealed in earlier comparisons, several recent changes have been added that improve the accuracy and flexibility of the GSLS model, including 1. improved treatment of the variation of the extinction coefficient with altitude, both within atmospheric layers and above the nominal top of the atmosphere; 2. addition of multiple-scattering source function calculations at multiple solar zenith angles along the line of sight (LOS); 3. introduction of variable surface properties along the limb LOS, with minimal effort required to add variable atmospheric properties along the LOS as well; 4. addition of the ability to model multiple aerosol types within the model atmosphere. The model improvements 1 and 2 are verified by comparison to previously published results (using standard radiance tables whenever possible), demonstrating significant improvement in cases for which previous versions of the GSLS model performed poorly. The single-scattered radiance errors that were as high as 4% in earlier studies are now generally reduced to 0.3%, while total radiance errors generally decline from 10% to 1–3%. In all cases, the tangent height dependence of the GSLS radiance error is greatly reduced.


2021 ◽  
Author(s):  
Ilaria Petracca ◽  
Davide De Santis ◽  
Stefano Corradini ◽  
Lorenzo Guerrieri ◽  
Matteo Picchiani ◽  
...  

<p>When an eruption event occurs it is necessary to accurately and rapidly determine the position and evolution during time of the volcanic cloud and its parameters (such as Aerosol Optical Depth-AOD, effective radius-Re and mass-Ma of the ash particles), in order to ensure the aviation security and the prompt management of the emergencies.</p><p>Here we present different procedures for volcanic ash cloud detection and retrieval using S3 SLSTR (Sentinel-3 Sea and Land Surface Temperature Radiometer) data collected the 22 June at 00:07 UTC by the Sentinel-3A platform during the Raikoke (Kuril Islands) 2019 eruption.</p><p>The volcanic ash detection is realized by applying an innovative machine learning based algorithm, which uses a MultiLayer Perceptron Neural Network (NN) to classify a SLSTR image in eight different surfaces/objects, distinguishing volcanic and weather clouds, and the underlying surfaces. The results obtained with the NN procedure have been compared with two consolidated approaches based on an RGB channels combination in the visible (VIS) spectral range and the Brightness Temperature Difference (BTD) procedure that exploits the thermal infrared (TIR) channels centred at 11 and 12 microns (S8 and S9 SLSTR channels respectively). The ash volcanic cloud is correctly identified by all the models and the results indicate a good agreement between the NN classification approach, the VIS-RGB and BTD procedures.</p><p>The ash retrieval parameters (AOD, Re and Ma) are obtained by applying three different algorithms, all exploiting the volcanic cloud “mask” obtained from the NN detection approach. The first method is the Look Up Table (LUT<sub>p</sub>) procedure, which uses a Radiative Transfer Model (RTM) to simulate the Top Of Atmosphere (TOA) radiances in the SLSTR thermal infrared channels (S8, S9), by varying the aerosol optical depth and the effective radius. The second algorithm is the Volcanic Plume Retrieval (VPR), based on a linearization of the radiative transfer equation capable to retrieve, from multispectral satellite images, the abovementioned parameters. The third approach is a NN model, which is built on a training set composed by the inputs-outputs pairs TOA radiances vs. ash parameters. The results of the three retrieval methods have been compared, considering as reference the LUT<sub>p</sub> procedure, since that it is the most consolidated approach. The comparison shown promising agreement between the different methods, leading to the development of an integrated approach for the monitoring of volcanic ash clouds using SLSTR.</p><p>The results presented in this work have been obtained in the sphere of the VISTA (Volcanic monItoring using SenTinel sensors by an integrated Approach) project, funded by ESA and developed within the EO Science for Society framework [https://eo4society.esa.int/projects/vista/].</p>


2013 ◽  
Vol 52 (3) ◽  
pp. 710-726 ◽  
Author(s):  
Chenxi Wang ◽  
Ping Yang ◽  
Steven Platnick ◽  
Andrew K. Heidinger ◽  
Bryan A. Baum ◽  
...  

AbstractA computationally efficient high-spectral-resolution cloudy-sky radiative transfer model (HRTM) in the thermal infrared region (700–1300 cm−1, 0.1 cm−1 spectral resolution) is advanced for simulating the upwelling radiance at the top of atmosphere and for retrieving cloud properties. A precomputed transmittance database is generated for simulating the absorption contributed by up to seven major atmospheric absorptive gases (H2O, CO2, O3, O2, CH4, CO, and N2O) by using a rigorous line-by-line radiative transfer model (LBLRTM). Both the line absorption of individual gases and continuum absorption are included in the database. A high-spectral-resolution ice particle bulk scattering properties database is employed to simulate the radiation transfer within a vertically nonisothermal ice cloud layer. Inherent to HRTM are sensor spectral response functions that couple with high-spectral-resolution measurements in the thermal infrared regions from instruments such as the Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer. When compared with the LBLRTM and the discrete ordinates radiative transfer model (DISORT), the root-mean-square error of HRTM-simulated single-layer cloud brightness temperatures in the thermal infrared window region is generally smaller than 0.2 K. An ice cloud optical property retrieval scheme is developed using collocated AIRS and Moderate Resolution Imaging Spectroradiometer (MODIS) data. A retrieval method is proposed to take advantage of the high-spectral-resolution instrument. On the basis of the forward model and retrieval method, a case study is presented for the simultaneous retrieval of ice cloud optical thickness τ and effective particle size Deff that includes a cloud-top-altitude self-adjustment approach to improve consistency with simulations.


2021 ◽  
Vol 14 (10) ◽  
pp. 6483-6507
Author(s):  
Zhao-Cheng Zeng ◽  
Vijay Natraj ◽  
Feng Xu ◽  
Sihe Chen ◽  
Fang-Ying Gong ◽  
...  

Abstract. Remote sensing of greenhouse gases (GHGs) in cities, where high GHG emissions are typically associated with heavy aerosol loading, is challenging due to retrieval uncertainties caused by the imperfect characterization of scattering by aerosols. We investigate this problem by developing GFIT3, a full physics algorithm to retrieve GHGs (CO2 and CH4) by accounting for aerosol scattering effects in polluted urban atmospheres. In particular, the algorithm includes coarse- (including sea salt and dust) and fine- (including organic carbon, black carbon, and sulfate) mode aerosols in the radiative transfer model. The performance of GFIT3 is assessed using high-spectral-resolution observations over the Los Angeles (LA) megacity made by the California Laboratory for Atmospheric Remote Sensing Fourier transform spectrometer (CLARS-FTS). CLARS-FTS is located on Mt. Wilson, California, at 1.67 km a.s.l. overlooking the LA Basin, and it makes observations of reflected sunlight in the near-infrared spectral range. The first set of evaluations are performed by conducting retrieval experiments using synthetic spectra. We find that errors in the retrievals of column-averaged dry air mole fractions of CO2 (XCO2) and CH4 (XCH4) due to uncertainties in the aerosol optical properties and atmospheric a priori profiles are less than 1 % on average. This indicates that atmospheric scattering does not induce a large bias in the retrievals when the aerosols are properly characterized. The methodology is then further evaluated by comparing GHG retrievals using GFIT3 with those obtained from the CLARS-GFIT algorithm (used for currently operational CLARS retrievals) that does not account for aerosol scattering. We find a significant correlation between retrieval bias and aerosol optical depth (AOD). A comparison of GFIT3 AOD retrievals with collocated ground-based observations from AErosol RObotic NETwork (AERONET) shows that the developed algorithm produces very accurate results, with biases in AOD estimates of about 0.02. Finally, we assess the uncertainty in the widely used tracer–tracer ratio method to obtain CH4 emissions based on CO2 emissions and find that using the CH4/CO2 ratio effectively cancels out biases due to aerosol scattering. Overall, this study of applying GFIT3 to CLARS-FTS observations improves our understanding of the impact of aerosol scattering on the remote sensing of GHGs in polluted urban atmospheric environments. GHG retrievals from CLARS-FTS are potentially complementary to existing ground-based and spaceborne observations to monitor anthropogenic GHG fluxes in megacities.


2016 ◽  
Author(s):  
Christopher E. Sioris ◽  
Landon A. Rieger ◽  
Nicholas D. Lloyd ◽  
Adam E. Bourassa ◽  
Chris Z. Roth ◽  
...  

Abstract. A new retrieval algorithm for OSIRIS (Optical Spectrograph and Infrared Imager System) nitrogen dioxide (NO2) profiles is described and validated. The algorithm relies on spectral fitting to obtain line-of-sight (LOS) column densities of NO2 followed by inversion using an algebraic reconstruction technique and the SaskTran spherical radiative transfer model to obtain vertical profiles of local number density. The validation covers different latitudes (tropical to polar), years (2002–2012), all seasons (winter, spring, summer, and autumn), different concentrations of nitrogen dioxide (from deNOxified polar vortex to polar summer), a large range of solar zenith angles (68.6 to 90.5°) and altitudes between 10.5 and 39 km, thereby covering the full retrieval range of a typical OSIRIS NO2 profile. The use of a larger spectral fitting window than used in previous retrievals reduces retrieval uncertainties and the scatter in the retrieved profiles due to noisy radiances. Improvements are also demonstrated through the validation in terms of bias reduction at 15–17 km relative to the OSIRIS operational v3.0 algorithm. By accounting for the diurnal variation along the LOS in the two-dimensional radiative transfer model, the scatter of the differences relative to the correlative balloon NO2 profile data is reduced.


Sign in / Sign up

Export Citation Format

Share Document