SEARCH · Search NASA

Results for “inversion and monitoring”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Development OF A Multi-Scale Framework for Mapping Global Evapotranspiration

As the worlds water resources come under increasing tension due to dual stressors of climate change and population growth, accurate knowledge of water consumption through evapotranspiration (ET) over a range in spatial scales will be critical in developing adaptation strategies. Remote sensing methods for monitoring consumptive water use (e.g, ET) are becoming increasingly important, especially in areas of significant water and food insecurity. One method to estimate ET from satellite-based methods, the Atmosphere Land Exchange Inverse (ALEXI) model uses the change in mid-morning land surface temperature to estimate the partitioning of sensible and latent heat fluxes which are then used to estimate daily ET. This presentation will outline several recent enhancements to the ALEXI modeling system, with a focus on global ET and drought monitoring.

Drought↗

NASA NTRSDate not supplied

Data Fusion and Mining Techniques to Map Water Use and Drought across Spatial and Temporal Scales

As the world’s water resources come under increasing tension due to dual stressors of climate change and population growth, accurate knowledge of water consumption through evapotranspiration (ET) over a range in spatial scales will be critical in developing adaptation strategies. Remote sensing methods for monitoring consumptive water use (e.g, ET) are becoming increasingly important, especially in areas of significant water and food insecurity. One method to estimate ET from satellite-based methods, the Atmosphere Land Exchange Inverse (ALEXI) model uses the change in mid-morning land surface temperature to estimate the partitioning of sensible and latent heat fluxes which are then used to estimate daily ET. This presentation will outline several recent enhancements to the ALEXI modeling system, with a focus on global ET and drought monitoring. Until recently, ALEXI has been limited to areas with high resolution temporal sampling of geostationary sensors. The use of geostationary sensors makes global mapping a complicated process, especially for real-time applications, as data from as many as five different sensors are required to be ingested and harmonized to create a global mosaic. However, our research team has developed a new and novel method of using twice-daily observations from polar-orbiting sensors such as MODIS and VIIRS to estimate the mid-morning rise in LST that is used to drive the energy balance estimations within ALEXI. This allows the method to be applied globally using a single sensor (in this case, initially MODIS with a planned transition to VIIRS) rather than a global compositing of all available geostationary data. Other advantages of this new method include the higher spatial resolution provided by MODIS and VIIRS and the increased sampling at high latitudes where oblique view angles limit the utility of geostationary sensors. This presentation will focus on global applications for mapping water use and drought using data mining and data fusion across spatial scales extending from 5-km to 30-m “field-scale” estimates.

Christopher Hain↗

NASA NTRS2016

Surface Pressure Dependencies in the Geos-Chem-Adjoint System and the Impact of the GEOS-5 Surface Pressure on CO2 Model Forecast

In the GEOS-Chem Adjoint (GCA) system, the total (wet) surface pressure of the GEOS meteorology is employed as dry surface pressure, ignoring the presence of water vapor. The Jet Propulsion Laboratory (JPL) Carbon Monitoring System (CMS) research team has been evaluating the impact of the above discrepancy on the CO2 model forecast and the CO2 flux inversion. The JPL CMS research utilizes a multi-mission assimilation framework developed by the Multi-Mission Observation Operator (M2O2) research team at JPL extending the GCA system. The GCA-M2O2 framework facilitates mission-generic 3D and 4D-variational assimilations streamlining the interfaces to the satellite data products and prior emission inventories. The GCA-M2O2 framework currently integrates the GCA system version 35h and provides a dry surface pressure setup to allow the CO2 model forecast to be performed with the GEOS-5 surface pressure directly or after converting it to dry surface pressure.

Carbon Monitoring System↗

NASA NTRS2013

Uncertainty Quantification in Remaining Useful Life of Aerospace Components using State Space Models and Inverse FORM

This paper investigates the use of the inverse first-order reliability method (inverse- FORM) to quantify the uncertainty in the remaining useful life (RUL) of aerospace components. The prediction of remaining useful life is an integral part of system health prognosis, and directly helps in online health monitoring and decision-making. However, the prediction of remaining useful life is affected by several sources of uncertainty, and therefore it is necessary to quantify the uncertainty in the remaining useful life prediction. While system parameter uncertainty and physical variability can be easily included in inverse-FORM, this paper extends the methodology to include: (1) future loading uncertainty, (2) process noise; and (3) uncertainty in the state estimate. The inverse-FORM method has been used in this paper to (1) quickly obtain probability bounds on the remaining useful life prediction; and (2) calculate the entire probability distribution of remaining useful life prediction, and the results are verified against Monte Carlo sampling. The proposed methodology is illustrated using a numerical example.

Sankararaman, Shankar↗

NASA NTRS2013

Discriminating Phytoplankton Functional Types (PFTs) in the Coastal Ocean Using the Inversion Algorithm Phydotax and Airborne Imaging Spectrometer Data

There is a need in the Biological Oceanography community to discriminate among phytoplankton groups within the bulk chlorophyll pool to understand energy flow through ecosystems, to track the fate of carbon in the ocean, and to detect and monitor-for harmful algal blooms (HABs). The ocean color community has responded to this demand with the development of phytoplankton functional type (PFT) discrimination algorithms. These PFT algorithms fall into one of three categories depending on the science application: size-based, biogeochemical function, and taxonomy. The new PFT algorithm Phytoplankton Detection with Optics (PHYDOTax) is an inversion algorithm that discriminates taxon-specific biomass to differentiate among six taxa found in the California Current System: diatoms, dinoflagellates, haptophytes, chlorophytes, cryptophytes, and cyanophytes. PHYDOTax was developed and validated in Monterey Bay, CA for the high resolution imaging spectrometer, Spectroscopic Aerial Mapping System with On-board Navigation (SAMSON - 3.5 nm resolution). PHYDOTax exploits the high spectral resolution of an imaging spectrometer and the improved spatial resolution that airborne data provides for coastal areas. The objective of this study was to apply PHYDOTax to a relatively lower resolution imaging spectrometer to test the algorithm's sensitivity to atmospheric correction, to evaluate capability with other sensors, and to determine if down-sampling spectral resolution would degrade its ability to discriminate among phytoplankton taxa. This study is a part of the larger Hyperspectral Infrared Imager (HyspIRI) airborne simulation campaign which is collecting Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) imagery aboard NASA's ER-2 aircraft during three seasons in each of two years over terrestrial and marine targets in California. Our aquatic component seeks to develop and test algorithms to retrieve water quality properties (e.g. HABs and river plumes) in both marine and in-land water bodies. Results presented are from the 10 April 2013 overflight of the Monterey Bay region and focus primarily on the first objective - sensitivity to atmospheric correction. On-going and future work will continue to evaluate if PHYDOTax can be applied to historical (SeaWiFS and MERIS), existing (MODIS, VIIRS, and HICO), and future (PACE, GEO-CAPE, and HyspIRI) satellite sensors. Demonstration of cross-platform continuity may aid in calibration and validation efforts of these sensors.

Biological Oceanography↗

NASA NTRS2024

Harnessing Collaborative Learning Automata to Guide Multi-objective Optimization based Inverse Analysis for Structural Damage Identification

Structural damage identification based on physical models is often transformed into an optimization problem that minimizes the difference between measurement information of structure being monitored and the model prediction in the parametric space. However, the objective function in this context often exhibits multimodality, involving high-dimensional variables due to the reliance on finite element models for damage identification. These features pose challenges to optimization algorithms, where entrapment in local solutions can lead to false positives and false negatives in damage identification. In this research, we propose a reinforcement learning based multi-swarm optimizer to tackle such challenges in pursuit of a small yet diverse solution set that can capture the true damage scenario as one of the solutions. The proposed method leverages the flexibility of the particle swarm optimizer and incorporates novel strategies of metaheuristics to realize targeted improvement. To enable the particle swarm to adaptively select the appropriate search strategy based on the current environment, we adopt the learning automata technique, which sidesteps the need for reward strategy selection that is usually ad hoc at each step of the search. The integration harnesses the automatic learning and self-adaptation capabilities of learning automata, enabling the particles to navigate based on environmental signals. This leads to accumulated probabilities tied to advantageous movements, fostering an adaptive exploration of particles in the search space. The proposed approach is first validated through implementing into benchmark test cases with comparisons. It is then applied to structural damage identification with piezoelectric admittance experimental signals. `The results highlight the capability of the algorithm to identify a small solution set with high accuracy to match the actual damage scenario.

Yang Zhang↗

NASA NTRS1981

Clear air turbulence studies with microwave radiometers

The flight testing of two passive microwave radiometers for the detection of clear air turbulence (CAT) is discussed. A 55.3 GHz radiometer was used to measure altitude temperature profiles, and a 180 GHz radiometer was used for monitoring line-of-sight integrated water vapor content. The sequence of altitude temperature profiles, spaced 17 seconds apart (3.5 km), enable inversion layer and tropopause properties to be studied. On some occasions, the altitude distribution of CAT severity correlated well with inversion layer location. On other occasions, turbulence was not located within 55.3 GHz measured layers, possibly indicating cases of topography-generated CAT, where inversion layer would not necessarily be expected. Evidence was obtained supporting the hypothesis that CAT is generated within layers containing levels of wind shear that cannot be supported by the layer's lapse rate; i.e., that Kelvin-Helmholtz wave breakdown can generate CAT. The 180 GHz radiometer failed to warn of CAT events. It is suggested that the radiometer's sensitivity (1 K) was inadequate for detecting the small variations in line-of-sight water vapor content.

Gary, B. L.↗

NASA NTRS1989

Online submicron particle sizing by dynamic light scattering using autodilution

Efficient production of a wide range of commercial products based on submicron colloidal dispersions would benefit from instrumentation for online particle sizing, permitting real time monitoring and control of the particle size distribution. Recent advances in the technology of dynamic light scattering (DLS), especially improvements in algorithms for inversion of the intensity autocorrelation function, have made it ideally suited to the measurement of simple particle size distributions in the difficult submicron region. Crucial to the success of an online DSL based instrument is a simple mechanism for automatically sampling and diluting the starting concentrated sample suspension, yielding a final concentration which is optimal for the light scattering measurement. A proprietary method and apparatus was developed for performing this function, designed to be used with a DLS based particle sizing instrument. A PC/AT computer is used as a smart controller for the valves in the sampler diluter, as well as an input-output communicator, video display and data storage device. Quantitative results are presented for a latex suspension and an oil-in-water emulsion.

Nicoli, David F.↗

NASA NTRS1995

Surface ozone levels at Table Mountain during STOIC 1989

As a part of the routine operations of the Jet Propulsion Laboratory atmospheric measurements program at the Table Mountain Facility, the surface ozone concentration is continuously monitored using a Dasibi photometer. The influence of the Los Angeles basin to the southwest of the facility and the height of the inversion layer cause large fluctuations in the ozone concentration. Peaks as high as 200 parts per billion by volume (ppbv) were observed during the Stratospheric Ozone Intercomparison Campaign (STOIC) compared to a normal background level near 50 ppbv. These measurements, made during STOIC, were important in assessing the impact of the surface ozone concentration on the various instruments participating in the campaign.

McDermid, I. Stuart↗

NASA NTRS2006

Chlorophyll Fluorescence Emissions of Vegetation Canopies From High Resolution Field Reflectance Spectra

A two-year experiment was performed on corn (Zea mays L.) crops under nitrogen (N) fertilization regimes to examine the use of hyperspectral canopy reflectance information for estimating chlorophyll fluorescence (ChlF) and vegetation production. Fluorescence of foliage in the laboratory has proven more rigorous than reflectance for correlation to plant physiology. Especially useful are emissions produced from two stable red and far-red chlorophyll ChlF peaks centered at 685V10 nm and 735V5 nm. Methods have been developed elsewhere to extract steady state solar induced fluorescence (SF) from apparent reflectance of vegetation canopies/landscapes using the Fraunhofer Line Depth (FLD) principal. Our study utilized these methods in conjunction with field-acquired high spectral resolution canopy reflectance spectra obtained in 2004 and 2005 over corn crops, as part of an ongoing multi-year experiment at the USDA/Agriculture Research Service in Beltsville, MD. A spectroradiometer (ASD-FR Fieldspec Pro, Analytical Spectral Devices, Inc., Boulder, CO) was used to measure canopy radiances 1 m above plant canopies with a 22deg field of view and a 0deg nadir view zenith angle. Canopy and plant measurements were made at the R3 grain fill reproductive stage on 3-4 replicate N application plots provided seasonal inputs of 280, 140, 70, and 28 kg N/ha. Leaf level measurements were also made which included ChlF, photosynthesis, and leaf constituents (photosynthetic pigment, carbon (C), and N contents). Crop yields were determined at harvest. SIF intensities for ChlF were derived directly from canopy reflectance spectra in specific narrowband regions associated with atmospheric oxygen absorption features centered at 688 and 760 nm. The red/far-red S F ratio derived from these field reflectance spectra successfully discriminated foliar pigment levels (e.g., total chlorophyll, Chl) associated with N application rates in both corn crops. This canopy-level spectral ratio was also positively correlated to the foliar C/N ratio (r = 0.89, n = go), as was a leaf-level steady state fluorescence ratio (Fs/Chl, r = 0.92). The latter ratio was inversely correlated with crop grain yield (Kg 1 ha) (r = 0.9). This study has relevance to future passive satellite remote sensing approaches to monitoring C dynamics from space.

Middleton, E. M.↗

NASA NTRS2016

Surface Pressure Dependencies in the GEOS-Chem-Adjoint System and the Impact of the GEOS-5 Surface Pressure on CO2 Model Forecast

Lee, Meemong↗

NASA NTRS2018

An Efficient Inverse Finite Element Method for Shape and Stress Sensing of Laminated Composite and Sandwich Plates and Shells

The inverse Finite Element Method (iFEM) is a revolutionary methodology for real-time reconstruction of full-field structural displacements and stresses in structures that are instrumented with strain sensors. This inverse problem is commonly referred to as shape and stress sensing, which is well-recognized as an enabling technology for structural health monitoring systems. In this study, an improved iFEM formulation is proposed for shape and stress sensing of laminated composite and sandwich plates and shells. The formulation includes the kinematics of a shear deformation plate theory known as Refined Zigzag Theory (RZT) as its baseline. The present iFEM formulation is based upon the minimization of a weighted-least-squares functional that uses the complete set of section strains of RZT. The improved iFEM methodology is applicable for shape and stress sensing of thin and moderately thick plate and shell structures involving a relatively small number of strain gauges. The main advantage of the current formulation is that highly accurate through-the-thickness distributions of displacements, strains, and stresses are attainable using an element based on simple C0-continuous displacement interpolation functions. A three-node inverseshell element, named i3-RZT, is developed. Two example problems are examined in detail: (1) a simply supported rectangular laminated composite plate and (2) a wedge structure with a hole near one of the clamped ends. For both problems, the experimental strain data are generated numerically by the direct finite element analysis using high-fidelity discretizations. These strains are then regarded as the experimental strains obtained from surface mounted strain gauges or embedded fiber Bragg grating (FBG) sensors. The numerical results demonstrate the superior capability and potential applicability of the i3- RZT/iFEM methodology for performing accurate shape and stress sensing of complex composite structures.

Kefal, Adnan↗

NASA NTRS2020

Non-linear Least Square Fitting Technique for the Determination of Field Line Resonance Frequency in Ground Magnetometer Data: Application to Remote Sensing of Plasmaspheric Mass Density

The accurate determination of the Field Line Resonance (FLR) frequency of a resonating geomagnetic field line is necessary to remotely monitor the plasmaspheric mass density during geomagnetic storms and quiet times alike. Under certain assumptions the plasmaspheric mass density at the equator is inversely proportional to the square of the FLR frequency. The most common techniques to determine the FLR frequency from ground magnetometer measurements are the amplitude ratio and phase difference techniques, both based on geomagnetic field observations at two latitudinally separated ground stations along the same magnetic meridian. Previously developed automated techniques have used statistical methods to pinpoint the FLR frequency using the amplitude ratio and phase difference calculations. We now introduce a physics-based automated technique, using non-linear least square fitting of the ground magnetometer data to the analytical resonant wave equations, that reproduces the wave characteristics on the ground, and from those determine the FLR frequency. One of the advantages of the new technique is the estimation of physics-based errors of the FLR frequency, and as a result of the equatorial plasmaspheric mass density. We present analytical results of the new technique, and test it using data from the Inner-Magnetospheric Array for Geospace Science (iMAGS) ground magnetometer chain along the coast of Chile and the east coast of the United States. We compare the results with the results of previously published statistical automated techniques.

A. Boudouridis↗

NASA NTRS2022

Correlated Spatio-temporal Evolution of Extreme-Ultraviolet Ribbons and Hard X-Rays in a Solar Flare

We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet images from the Interface Region Imaging Spectrograph and the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA), magnetic data from the SDO/Helioseismic and Magnetic Imager, hard X-ray (HXR) images from the Reuven Ramaty High Energy Solar Spectroscopic Imager, and light curves from the Fermi/Gamma-ray Burst Monitor, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. The width of the newly brightened front along the extension of the ribbon is highly intermittent in both space and time, presumably reflecting nonuniformities in the structure and/or dynamics of the flare current sheet. Furthermore, the ribbon width grows most rapidly in regions exhibiting concentrated nonthermal HXR emission, with sharp increases slightly preceding the HXR bursts. The light curve of the ultraviolet emission matches the HXR light curve at photon energies above 25 keV. In other regions the ribbon-width evolution and light curves do not temporally correlate with the HXR emission. This indicates that the production of nonthermal electrons is highly nonuniform within the flare current sheet. Our results suggest a strong connection between the production of nonthermal electrons and the locally enhanced perpendicular extent of flare ribbon fronts, which in turn reflects the inhomogeneous structure and/or reconnection dynamics of the current sheet. Despite this variability, the ribbon fronts remain nearly continuous, quasi-one-dimensional features. Thus, although the reconnecting coronal current sheets are highly structured, they remain quasi-two-dimensional and the magnetic energy release occurs systematically, rather than stochastically, through the volume of the reconnecting magnetic flux.

S. J. Naus↗

NASA NTRS1985

The investigation of advanced remote sensing, radiative transfer and inversion techniques for the measurement of atmospheric constituents

The research program is documented for developing space and ground-based remote sensing techniques performed during the period from December 15, 1977 to March 15, 1985. The program involved the application of sophisticated radiative transfer codes and inversion methods to various advanced remote sensing concepts for determining atmospheric constituents, particularly aerosols. It covers detailed discussions of the solar aureole technique for monitoring columnar aerosol size distribution, and the multispectral limb scattered radiance and limb attenuated radiance (solar occultation) techniques, as well as the upwelling scattered solar radiance method for determining the aerosol and gaseous characteristics. In addition, analytical models of aerosol size distribution and simulation studies of the limb solar aureole radiance technique and the variability of ozone at high altitudes during satellite sunrise/sunset events are also described in detail.

Deepak, Adarsh↗

NASA NTRS2011

SO2 Emissions and Lifetimes: Estimates from Inverse Modeling Using In Situ and Global, Space-Based (SCIAMACHY and OMI) Observations

Top-down constraints on global sulfur dioxide (SO2) emissions are inferred through inverse modeling using SO2 column observations from two satellite instruments (SCIAMACHY and OMI). We first evaluated the S02 column observations with surface SO2 measurements by applying local scaling factors from a global chemical transport model (GEOS-Chem) to SO2 columns retrieved from the satellite instruments. The resulting annual mean surface SO2 mixing ratios for 2006 exhibit a significant spatial correlation (r=0.86, slope=0.91 for SCIAMACHY and r=0.80, slope = 0.79 for OMI) with coincident in situ measurements from monitoring networks throughout the United States and Canada. We evaluate the GEOS-Chem simulation of the SO2 lifetime with that inferred from in situ measurements to verity the applicability of GEOS-Chem for inversion of SO2 columns to emissions. The seasonal mean SO2 lifetime calculated with the GEOS-Chem model over the eastern United States is 13 h in summer and 48 h in winter, compared to lifetimes inferred from in situ measurements of 19 +/- 7 h in summer and 58 +/- 20 h in winter. We apply SO2 columns from SCIAMACHY and OMI to derive a top-down anthropogenic SO2 emission inventory over land by using the local GEOS-Chem relationship between SO2 columns and emissions. There is little seasonal variation in the top-down emissions (<15%) over most major industrial regions providing some confidence in the method. Our global estimate for annual land surface anthropogenic SO2 emissions (52.4 Tg S/yr from SCIAMACHY and 49.9 Tg S / yr from OMI) closely agrees with the bottom-up emissions (54.6 Tg S/yr) in the GEOS-Chem model and exhibits consistency in global distributions with the bottom-up emissions (r = 0.78 for SCIAMACHY, and r = 0.77 for OMI). However, there are significant regional differences.

Lee, Chulkyu↗

NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗