1 / 22

Elisa Carboni

Spectrally High resolution Infrared measurements for the characterisation of Volcanic Ash (SHIVA): a new way to study volcanic processes. Elisa Carboni. Objective. Identify ash composition from information contained in infrared spectra to study their change during an eruption,

latona
Download Presentation

Elisa Carboni

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Spectrally High resolution Infrared measurements for the characterisation of Volcanic Ash (SHIVA): a new way to study volcanic processes Elisa Carboni

  2. Objective Identify ash composition from information contained in infrared spectra to study their change during an eruption, to better understand the volcanic process that control eruptive activity. How we do it: Developing algorithms using new ash optical properties to retrieve ash characteristics from data collected with three different hyperspectral observing systems: ground-based absorption spectra FTIRand satellite-based emission spectra (from the limb MIPASand from nadir IASI). Validating the results, with geochemical and petrological analysis of ash samples and with reference to laboratory refractive index measurements of ash samples.

  3. Scientific Background Component and morphological analyses of the erupted ash, and comparison of these data with those from other monitoring techniques, demonstrates a clear relationship between ash features and styles of explosive activity. . Basaltic Andesitic Rhyolitic Increasing SiO2% Increasing explosivityof volcanic activity Extinction coefficients obtained for ice, volcanic ash (Peters, 2012), pumice, andesite, H2SO4, quartz, obsidian.

  4. Scientific Background IASI observations of three dust storms and seven volcanic ash plumes (Kasatochi 8/2008, Eyjafjallajoküll 4/2010, Sarychev 6/2009, Montserrat 2/2010, Kliuchevskoi 6/2010 and Chaitén 4/2008) from Clarisse et al. (2010).

  5. Method Optimal estimation approach (Rodgers, 2000) to retrieve ash composition and possibly size from infrared spectral measurements. We will study ash formed from different magmas and at different stages of evolution within a volcanic plume. one-to-one correspondence between refractive index spectra, compositions and remote sensing measurements. FTIR MIPAS IASI ground and satellite retrievalsash type/composition (with different size distributions) to be compared with the analysis done on ash sample for the same volcanic plume. + Ash sample analysis

  6. Ash optical properties – WP1 literature survey + Dan meas. + TEM, SEM meas. ash refractive index, morphology and size Mie(spherical particles) or T-matrix (non-spherical) spectral extinction coefficient, spectral single scattering albedo and phase function Opt. prop. simulations of spectrometer measurements Spectral extinction coefficient for different effective radius, obtained using the same volcanic ash refractive index (Peters, 2013). We can consider: - Log-normal size distributions, 2 or more components - External or internal mixing - Coating of the ash particles with water, ice or sulphuric acid one to one Deliverable: Database of spectral ash and volcanic aerosol optical properties available through the project web page.

  7. Ground transmittance analysis (FTIR) – WP2 Fourier transform infrared spectroscopy (FTIR) transmittance measurements from 800 to 5000 cm-1. Project partners and measurements available. Deliverable: (1) Forward model to reproduce the FTIR measurements and retrieval code for analysis. (2) A publication describing the ash/aerosol retrieval using FTIRmeasurements and including the results of the analysis. Forward Model Options: - Thin plume -> a single scattering model - Thicker plume -> DISORT (include multiple scattering and thermal emission) - others? Problems: transmittance? Radiometric drift? Different sources? Retrieved parameters: ash optical depth AODand effective radius Reff and possibly other (mode radius, spread and components) depending on the information content

  8. IASI – WP3 We have a forward model that simulate IASI spectra with ash and SO2 and can include two separate layers (e.g. ash above cloud). The RTTOV output for a clean atmosphere (containing gas but not cloud or aerosol/ash) will be combined with an the ash/cloud layer using the ORAC scheme Retrieved parameters: ash optical depth (at a reference wavelength tbd), ash effective radius, ash type and layer altitude. We will assess the ability of the retrieval to distinguish the type of layer over a range of ash-aerosol-cloud mixtures • Comparison with: • (1) independent IASI retrieval (project partner: University of Bristol), • (2) ASTER ash retrieval (project partner: University of Bristol) • (3) ORACash and cloud retrievals from the SEVIRI and AATSR. • (4) independent ash retrieval from SEVIRI and MODIS (project partner INGV) • (5) CALIPSOmeasurements. Deliverables:(1) retrieval code for IASI ash; (2) a publication.

  9. IASI FAST FORWARD MODEL Tac Iac Iac atm. above We need to divide the contributions between atmosphere and aerosol layer components The suffix ‘c’ refers to aerosol layer Rc, ec, Tc ,Bc aerosol layer Ibc ,Tbc atm. belowe surface Rs, es, Bs b a c d LUTs for aerosol layer parameters Other atmospheric parameters (radiances above/below aerosol layer going up/down) are computed with RTTOV using ECMWF atmospheric profiles.

  10. Ash and cloud IASI forward model IASI forward model include aerosol and cloud. Obtained with ‘Aso’ refractive index measured by Daniel Peters. Ash Re=2 [mm], h=3 [km] Water Cloud Re=10 [mm], h=3 [km] Dust Re=2 [mm], h=3 [km]

  11. MIPAS – WP4 The spatial distribution of measurements from MIPAS is more sparse than IASI, but provides accurate information on the ash or cloud top layer altitude (Hurley et al. 2011, Spang et al. ACP 2012). Two approaches for ash retrieval will be tested: (1) Fitting the measurements using an ash MIPAS forward model (to be developed) and the optical properties from WP1. (2) Using a singular vector decomposition (SVD) technique, as successfully applied for cloud retrieval (Hurley et al., 2011). Deliverables: (1) retrieval code for MIPAS ash; (2) A publication describing the ash/aerosol retrieval using MIPAS measurements and including the results of the analysis.

  12. Validation/comparison and interpretation– WP5 - Test consistency of retrieved composition from different sensors (IASI, MIPAS, FTIR) - Compare remotely sensed composition (ash type) with samples taken on ground. - Case studies to observe how ash/aerosol size distribution and composition changes as a function of time (or others parameters) to obtain information on the volcanic processes. - Improve prediction of ash characteristcs by studying relationships between ash concentration and size with other parameters: - measures of volcanic activity (plume height/seismicity etc) - the involvement of external water sources., - the gas (e.g. SO2) flux etc.

  13. Wish list: output of the day ? How we can define a finite number of ash/aerosol model (3-10?) that correspond to ash composition/optical properties to be used in the retrievals? (possibly distinguishable ash model) How to link satellite and ground base measurements? Can we relate FTIR measurements with satellite? With ash sample measurements? What is more interesting to retrieve/study? What can tell us something about volcanic processes? (aerosol type/ash composition, amount, dimension, variation in time?) Define tasks, deadline and who will work on them. (tasks should become papers!?!) Define how we will be in contact after the workshop: mail list? Teleconference? Others?

  14. ash cloud can be discriminated from other dust or aerosol clouds? : Aod=[0.2,0.5,1,1.5,2]; Reff=[0,2,0.5,1, 1.5,2,3,4,6,10] micron; H=[0,2,6,9,12,15] km. All simulations of IASI ash spectra have been used to compute the covariance matrix Sash and the mean spectra mashof ash. c2 = (ydust-mash)TSash-1 (ydust-mash). c2 between simulated IASI spectra with dust and ash, as function of different dust AOD and altitude Values higher then 1700 (dark blu) have probability less the 0.04 to be obtained (with the number of channel considered e.g. 1601), it means that the dust spectra is not fitted by the ash model, so ash and dust (with the assumed refractive index and optical properties) are distinguishable.

  15. For the same state vector the ash and dust case have different spectral shape of the jacobian. (IASI FM)

  16. RETRIEVAL METHOD OPTIMAL ESTIMATION [Rodgers 2000] Initial state estimate: x0 A priori: xa Run forward model: f(xi) Compare to J = [y - f(xi)]Se-1[y - f(xi)] + measurements (y): [xi - xa]Sa-1[xi - xa] Update state: xi → xi+1(Levenburg-Marquardt) Stop when: J is small , or when i is large. NB Optimal estimation method provides quality control and error estimate

  17. OE detection theory [Rodger 2000] The optimal estimate of x taking into account total measurement error may be computed as: Sytot is computed considering an appropriate ensemble of N measured spectra to construct an estimate of total measurement error variance-covariance Syobs [ Walker, Dudhia, Carboni, Atmos. Meas. Tech. Discuss., 2010 ]

  18. Thermal infrared spectra ash SO2 absorption bands Brightness Temperature [k] 1 3 Wavenumber [cm-1]

  19. AEROSOL OPTICAL PROPERTIES Microphysical properties Every component is characterized by: Spectral refractive index m(l) + i k(l) Mode radius rm and spread s Aerosol class with more components Changing the mixing ratio between component we obtain the optical properties corresponding to different effective radius Kext() () < 1 P(, ) Volcanic ash: 1 component only Changing the mode radius rm (= tranlation of the log normal)

More Related