Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interferometric Synthetic Aperture Radars (InSAR)”

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Registration of interferometric SAR images

Interferometric synthetic aperture radar (INSAR) is a new way of performing topography mapping. Among the factors critical to mapping accuracy is the registration of the complex SAR images from repeated orbits. A new algorithm for registering interferometric SAR images is presented. A new figure of merit, the average fluctuation function of the phase difference image, is proposed to evaluate the fringe pattern quality. The process of adjusting the registration parameters according to the fringe pattern quality is optimized through a downhill simplex minimization algorithm. The results of applying the proposed algorithm to register two pairs of Seasat SAR images with a short baseline (75 m) and a long baseline (500 m) are shown. It is found that the average fluctuation function is a very stable measure of fringe pattern quality allowing very accurate registration.

Lin, Qian↗
NASA NTRS2006

Rapid, Repeat-sample Monitoring of Crustal Deformations and Environmental Phenomena with the Uninhabited Aerial Vehicle Synthetic Aperture Radar

The Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) is a precision repeat-pass Interferometric Synthetic Aperture Radar (InSAR) mission being developed by the Jet Propulsion Laboratory and the Dryden Flight Research Center in support of NASA s Science Mission Directorate. UAVSAR's unique ability to fly a repeatable flight path, along with an electronically steerable array, allows interferometric data to be obtained with accuracies measured in millimeters. Deploying the radar on an airborne platform will also allow for radar images to be collected and compared with images from the same area taken hours or even years later - providing for long-term trending and near real-time notification of changes and deformations. UAVSAR s data processing algorithms will provide for near-real time data reduction providing disaster planning and response teams with highly accurate data to aid in the prediction of, and response to, natural phenomena. UAVSAR data can be applied to increasing our understanding of the processes behind solid earth, cryosphere, carbon cycle and other areas of interest in earth science. Technologies developed for UAVSAR may also be applicable to a future earth-orbiting InSAR mission and possibly for missions to the Moon or Mars. The UAVSAR is expected to fly on a Gulfstream III aircraft this winter, followed by a flight test program lasting until the second half of 2007. Following radar calibration and data reduction activities, the platform will be ready for science users in the summer of 2008.

Smith, Robert C.↗
DOE OSTI2021

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Time–series Analysis of Volume Change at Brady Hot Springs, Nevada, USA using Geodetic Data from 2003 – 2018

Brady Hot Springs geothermal field has exhibited subsidence, as measured by interferometric synthetic aperture radar (InSAR). Previous studies have examined both the temporal evolution of the deformation from 2004 through 2016 and the spatial extent of the deformation, directly relating the observed subsidence to volumetric changes below the surface. We extend the modeling at Brady to analyze a data set of interferometric pairs spanning from the end of 2003 through 2018. We examine spatial and temporal trends in the observed deformation by time–series analysis of each of the 1656 cubic voxels in a parameterized elastic dislocation model to identify areas where the subsurface volume changes as a function of time. Joint time–series analysis of Global Positioning System and InSAR pairs confirm significant changes in rates of volume change during time intervals when well operations were varied. Here, the rate of subsidence increases with increased injection, consistent with the identification of thermal contraction of the rock matrix as the dominant driving mechanism. Conversely, the modeled volume increases when pumping ceases, suggesting thermal expansion of the rock matrix.

58 GEOSCIENCES↗
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↗
DOE OSTI2022

Shallow Aseismic Slip in the Delaware Basin Determined by Sentinel-1 InSAR

The Delaware Basin, Texas is currently a hot-spot of induced seismicity and ground deformation due to fluid extraction and injection associated with horizontal drilling techniques; however, the driving mechanism behind the seismicity and deformation remains under debate. Here, using vertical and east-west horizontal surface deformation measurements derived from Sentinel-1 interferometric synthetic aperture radar (InSAR), we show that the subsurface responds differently to oil and gas activity in the northern and southeastern portions of the basin. In the north, where there is little seismicity, deformation patterns display long-wavelengths and equidimensional patterns. In contrast, the southeast region hosts most of the seismicity and displays spatial deformation patterns with narrow linear features that strike parallel to the maximum principal horizontal stress and to trends in seismicity, suggesting movement along normal faults. We model a linear deformation feature using edge dislocations and show that the InSAR observations can be reproduced by slip on normal faults contained within the Delaware Mountain Group (DMG), the formation that hosts local wastewater injection and the majority of earthquakes. Our model consists of three parallel, high-angle normal faults, with two dipping toward one another in a graben structure. Slip magnitudes reach up to 25 cm and are spatially correlated with injection wells. Measured seismicity can only explain ~2% of the fault motion predicted by our fault model, suggesting that slip leading to the deformation is predominantly aseismic. We conclude that seismic and aseismic fault motion in the southeastern Delaware Basin is likely driven by wastewater injection near critically-stressed normal faults within the DMG.

58 GEOSCIENCES↗
NASA NTRS2007

Radar Interferometer for Topographic Mapping of Glaciers and Ice Sheets

A report discusses Ka-band (35-GHz) radar for mapping the surface topography of glaciers and ice sheets at high spatial resolution and high vertical accuracy, independent of cloud cover, with a swath-width of 70 km. The system is a single- pass, single-platform interferometric synthetic aperture radar (InSAR) with an 8-mm wavelength, which minimizes snow penetration while remaining relatively impervious to atmospheric attenuation. As exhibited by the lower frequency SRTM (Shuttle Radar Topography Mission) AirSAR and GeoSAR systems, an InSAR measures topography using two antennas separated by a baseline in the cross-track direction, to view the same region on the ground. The interferometric combination of data received allows the system to resolve the pathlength difference from the illuminated area to the antennas to a fraction of a wavelength. From the interferometric phase, the height of the target area can be estimated. This means an InSAR system is capable of providing not only the position of each image point in along-track and slant range as with a traditional SAR but also the height of that point through interferometry. Although the evolution of InSAR to a millimeter-wave center frequency maximizes the interferometric accuracy from a given baseline length, the high frequency also creates a fundamental problem of swath coverage versus signal-to-noise ratio. While the length of SAR antennas is typically fixed by mass and stowage or deployment constraints, the width is constrained by the desired illuminated swath width. As the across-track beam width which sets the swath size is proportional to the wavelength, a fixed swath size equates to a smaller antenna as the frequency is increased. This loss of antenna size reduces the two-way antenna gain to the second power, drastically reducing the signal-to-noise ratio of the SAR system. This fundamental constraint of high-frequency SAR systems is addressed by applying digital beam-forming (DBF) techniques to synthesize multiple simultaneous receive beams in elevation while maintaining a broad transmit illumination. Through this technique, a high antenna gain on receive is preserved, thereby reducing the required transmit power and thus enabling high-frequency SARs and high-precision InSAR from a single spacecraft.

Moller, Delwyn K.↗
NASA NTRS2018

Multibaseline TanDEM-X Mangrove Height Estimation: The Selection of the Vertical Wavenumber

We generated a large-scale mangrove forest height map using multiple TanDEM-X (TDX) interferometric synthetic aperture radar (InSAR) acquisitions with various spatial baselines in order to improve the height estimation accuracy across a wide range of forest heights. The forest height inversion using InSAR data is strongly dependent upon the vertical wavenumber (i.e., perpendicular baseline). First, we investigated the role of the vertical wavenumber in forest height inversion from InSAR data using the sensitivity of the interferometric (volume) coherence to forest height. We used corrected but lower resolution and accuracy Shuttle Radar Topography Mission (SRTM) mangrove height maps as a priori information over Akanda and Pongara National Parks in Gabon to estimate lower and upper boundaries of the vertical wavenumber over test sites from the measured coherence-to-height sensitivity. Only TDX acquisitions within the boundaries of the vertical wavenumber were selected and combined for multibaseline mangrove height inversion. Mangrove forest height was obtained with multibaseline TDX acquisitions and was validated against the reference height derived from field measurement data providing improvements in multibaseline inversion over existing height estimates (i.e., SRTM height) and single-baseline inversions (multibaseline inversion: r(exp 2) = 0.98, root mean square error (RMSE) of 2.73 m; SRTM height: r(exp 2) = 0.86, RMSE = 7.21 m; single-baseline inversions: r(exp 2) = 0.08-0.97, RMSE = 3.86-11.10 m). As a result, to accurately estimate forest heights over a wide range (3-60 m), multibaseline InSAR acquisitions (at least three different baselines) are needed to exclude biases associated with the vertical wavenumber in forest height inversion.

vertical wavenumber.↗
NASA NTRSDate not supplied

Estimating Forest Stand Height in Savannakhet, Lao PDR Using InSAR and Backscatter Methods with L-Band SAR Data

Forest stand height (FSH), or average canopy height, serves as an important indicator for forest monitoring. The information provided about above ground biomass for greenhouse gas emissions reporting, and estimating carbon storage is relevant for reporting for Reducing Emissions from Deforestation and Forest Degradation (REDD+). A novel forest height estimation method utilizing a fusion of backscatter and Interferometric Synthetic Aperture Radar (InSAR) data from JAXA’s Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar (ALOS PALSAR) are applied to a use case in Savannakhet, Lao. Compared with LiDAR, the estimated height from the fusion method had an RMSE of 4.90 m and an R2 of 0.26. These results are comparable to previous studies using SAR estimation techniques. Despite limitations of data quality and quantity, the Savannakhet, Lao use case demonstrates the applicability of these techniques utilizing L-band SAR data for estimating FSH in tropical forests and can be used as a springboard for use of L-band data from the future NASA-ISRO SAR (NISAR) mission.

remote sensing↗
DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Calibration of the Geosar Dual Frequency Interferometric SAR

GeoSAR is an airborne, interferometric Synthetic Aperture Radar (INSAR) system for terrain mapping, currently under development by a consortium including NASA's Jet Propulsion Laboratory (JPL), Calgis, Inc., and the California Department of Conservation (CalDOC) with funding provided by the Topographic Engineering Center (TEC) of the U.S. Army Corps of Engineers and the Defense Advanced Research Projects Agency (DARPA). The radar simultaneously maps swaths on both sides of the aircraft at two frequencies, X-Band and P-Band. For the P-Band system, data is collected for two across track interferometric baselines and at the crossed polarization. The aircraft position and attitude are measured using two Honeywell Embedded GPS Inertial Navigation Units (EGI) and an Ashtech Z12 GPS receiver. The mechanical orientation and position of the antennas are actively measured using a Laser Baseline Metrology System (LBMS). In the GeoSAR motion measurement software, these data are optimally combined with data from a nearby ground station using Ashtech PNAV software to produce the position, orientation, and baseline information are used to process the dual frequency radar data. Proper calibration of the GeoSAR system is essential to obtaining digital elevation models (DEMS) with the required sub-meter level planimetric and vertical accuracies. Calibration begins with the determination of the yaw and pitch biases for the two EGI units. Common range delays are determined for each mode, along with differential time and phase delays between channels. Because the antennas are measured by the LBMS, baseline calibration consists primarily of measuring a constant offset between mechanical center and the electrical phase center of the antennas. A phase screen, an offset to the interferometric phase difference which is a function of absolute phase, is applied to the interferometric data to compensate for multipath and leakage. Calibration parameters are calculated for each of the ten processing modes, each of the operational bandwidths (80 and 160 MHZ), and each aircraft altitude. In this talk we will discuss the layout calibration sites, the synthesis of data from multiple flights to improve the calibration, methods for determining time and phase delays, and techniques for determining radiometric and polarimetric quantities. We will describe how calibration quantities are incorporated into the processor and pre-processor. We will demonstrate our techniques applied to GeoSar data and assess the stability and accuracy of the calibration. This will be compared to the modeled performance determined from calibrating the output of a point target simulator. The details of baseline determination and phase screen calculation are covered in related talks.

Chapine, Elaine↗
NASA NTRS2000

Estimating Forest Vertical Structure from Multialtitude, Fixed-Baseline Radar Interferometric and Polarimetric Data

Parameters describing the vertical structure of forests, for example tree height, height-to-base-of-live-crown, underlying topography, and leaf area density, bear on land-surface, biogeochemical, and climate modeling efforts. Single, fixed-baseline interferometric synthetic aperture radar (INSAR) normalized cross-correlations constitute two observations from which to estimate forest vertical structure parameters: Cross-correlation amplitude and phase. Multialtitude INSAR observations increase the effective number of baselines potentially enabling the estimation of a larger set of vertical-structure parameters. Polarimetry and polarimetric interferometry can further extend the observation set. This paper describes the first acquisition of multialtitude INSAR for the purpose of estimating the parameters describing a vegetated land surface. These data were collected over ponderosa pine in central Oregon near longitude and latitude -121 37 25 and 44 29 56. The JPL interferometric TOPSAR system was flown at the standard 8-km altitude, and also at 4-km and 2-km altitudes, in a race track. A reference line including the above coordinates was maintained at 35 deg for both the north-east heading and the return southwest heading, at all altitudes. In addition to the three altitudes for interferometry, one line was flown with full zero-baseline polarimetry at the 8-km altitude. A preliminary analysis of part of the data collected suggests that they are consistent with one of two physical models describing the vegetation: 1) a single-layer, randomly oriented forest volume with a very strong ground return or 2) a multilayered randomly oriented volume; a homogeneous, single-layer model with no ground return cannot account for the multialtitude correlation amplitudes. Below the inconsistency of the data with a single-layer model is followed by analysis scenarios which include either the ground or a layered structure. The ground returns suggested by this preliminary analysis seem too strong to be plausible, but parameters describing a two-layer compare reasonably well to a field-measured probability distribution of tree heights in the area.

Treuhaft, Robert N.↗
NASA NTRS2000

Shuttle Radar Topography Mission (SRTM) Flight System Design and Operations Overview

This paper provides an overview of the Shuttle Radar Topography Mission (SRTM), with emphasis on flight system implementation and mission operations from systems engineering perspective. Successfully flown in February, 2000, the SRTM's primary payload consists of several subsystems to form the first spaceborne dual-frequency (C-band and X-band) fixed baseline interferometric synthetic aperture radar (InSAR) system, with the mission objective to acquire data sets over 80% of Earth's landmass for height reconstruction. The paper provides system architecture, unique design features, engineering budgets, design verification, in-flight checkout and data acquisition of the SRTM payload, in particular for the C-band system. Mission operation and post-mission data processing activities are also presented. The complexity of the SRTM as a system, the ambitious mission objective, the demanding requirements and the high interdependency between multi-disciplined subsystems posed many challenges. The engineering experience and the insight thus gained have important implications for future spaceborne interferometric SAR mission design and implementation.

Shen, Yuhsyen↗
DOE OSTI2023

A Fusion of Geothermal and InSAR Data with Machine Learning for Enhanced Deformation Forecasting at the Geysers

The Geysers geothermal field in California is experiencing land subsidence due to the seismic and geothermal activities taking place. This poses a risk not only to the underlying infrastructure but also to the groundwater level which would reduce the water availability for the local community. Because of this, it is crucial to monitor and assess the surface deformation occurring and adjust geothermal operations accordingly. In this study, we examine the correlation between the geothermal injection and production rates as well as the seismic activity in the area, and we show the high correlation between the injection rate and the number of earthquakes. This motivates the use of this data in a machine learning model that would predict future deformation maps. First, we build a model that uses interferometric synthetic aperture radar (InSAR) images that have been processed and turned into a deformation time series using LiCSBAS, an open-source InSAR time series package, and evaluate the performance against a linear baseline model. The model includes both convolutional neural network (CNN) layers as well as long short-term memory (LSTM) layers and is able to improve upon the baseline model based on a mean squared error metric. Then, after getting preprocessed, we incorporate the geothermal data by adding them as additional inputs to the model. This new model was able to outperform both the baseline and the previous version of the model that uses only InSAR data, motivating the use of machine learning models as well as geothermal data in assessing and predicting future deformation at The Geysers as part of hazard mitigation models which would then be used as fundamental tools for informed decision making when it comes to adjusting geothermal operations.

Yazbeck, Joe (ORCID:0000000302235260)↗
NASA NTRS2009

UAVSAR Phased Array Aperture

This paper describes the development of a patch antenna array for an L-band repeat-pass interferometric synthetic aperture radar (InSAR) instrument that is to be flown on an unmanned aerial vehicle (UAV). The antenna operates at a center frequency of 1.2575 GHz and with a bandwidth of 80 MHz, consistent with a number of radar instruments that JPL has previously flown. The antenna is designed to radiate orthogonal linear polarizations in order to facilitate fully-polarimetric measurements. Beam-pointing requirements for repeat-pass SAR interferometry necessitate electronic scanning in azimuth over a range of -20degrees in order to compensate for aircraft yaw. Beam-steering is accomplished by transmit/receive (T/R) modules and a beamforming network implemented in a stripline circuit board. This paper, while providing an overview of phased array architecture, focuses on the electromagnetic design of the antenna tiles and associated interconnects. An important aspect of the design of this antenna is that it has an amplitude taper of 10dB in the elevation direction. This is to reduce multipath reflections from the wing that would otherwise be detrimental to interferometric radar measurements. This taper is provided by coupling networks in the interconnect circuits as opposed to attenuating the output of the T/R modules. Details are given of material choices and fabrication techniques that meet the demanding environmental conditions that the antenna must operate in. Predicted array performance is reported in terms of co-polarized and crosspolarized far-field antenna patterns, and also in terms of active reflection coefficient.

Unmanned Aerial Vehicle (UAV)↗
DOE OSTI2020

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗
NASA NTRS2020

Characterization of the 2008 Phreatomagmatic Eruption of Okmok From ArcticDEM and InSAR: Deposition, Erosion, and Deformation

Okmok volcano, one of the most active volcanoes in the Aleutian Islands, erupted between 12 July and 19 August 2008 and produced a new tephra cone (Ahmanilix) up to 253 m high. A coeruptive deposit thickness map and the posteruptive elevation change rate map are generated for this phreatomagmatic eruption using the new, high-resolution (2 m) digital elevation model (DEM) time series provided by ArcticDEM. The total dense rock equivalent (DRE) volume over the proximal deposit field classified from the elevation change map is estimated as 0.147 ± 0.001 km 3 DRE (bulk volume of 0.453 ± 0.003 km 3 ). Combined with the distal deposit volume derived from field measurements, the total DRE is 0.25 km 3 (bulk volume of 0.76 km 3 ). The estimated erosion rate is up to −15 ± 3 m/year along the flanks of Ahmanilix, and the redeposition rate is about 1 ± 0.3 m/year at the base and 6 ± 1 m/year within the center of Ahmanilix. The surface elevation change also reveals a large ground surface collapse that disrupted the old Cone D lava bench, with a total collapse area of 0.15 km2 and a volume of (7.1 ± 0.4) × 10−3 km 3 . We also produce a coeruptive deformation map from satellite interferometric synthetic aperture radar (InSAR) using topography from both ArcticDEM and the Shuttle Radar Topography Mission (SRTM) and compare the results. Finally, based on historical eruption volumes, the magma supply rate at Okmok is estimated as (1.77 ± 0.1) × 10−3 km 3 /year.

Chunli Dai↗