Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 19 records

Leveraging Multi-Primary PS-InSAR Configurations for the Robust Estimation of Coastal Subsidence

InSAR is a key technique used to constrain contributions of diverse processes to coastal subsidence, also known as vertical land motion (VLM). However, coastal environments can pose major challenges for InSAR due to natural disturbances that degrade interferogram quality. We describe a new multi-primary pairing strategy for persistent scatterer InSAR (PS-InSAR) to estimate subsidence in challenging coastal environments. Our method retains only consistent PS candidates across multiprimary substacks and solves for redundant velocity observations using SVD-based inversion, similar to the conventional small baseline subset (SBAS) method. Through simulations and a case study comparing with single-primary PS-InSAR and conventional SBAS techniques, we show that our pairing strategy reduces temporal and spatial uncertainty in subsidence estimates in the presence of strong but temporary decorrelation loss, even with increased distance from the reference point. Moreover, our method visibly dampens time-series variation and decreases standard error in our time series fit by nearly 2x in our case study. Thus, we find that implementing a multi-primary PS-InSAR configuration is a simple method of increasing the robustness of VLM estimates in challenging coastal environments.

Persistent scatterer InSAR↗

Time Series Analysis of Insar Data: Methods and Trends

Time series analysis of InSAR data has emerged as an important tool for monitoring and measuring the displacement of the Earth's surface. Changes in the Earth's surface can result from a wide range of phenomena such as earthquakes, volcanoes, landslides, variations in ground water levels, and changes in wetland water levels. Time series analysis is applied to interferometric phase measurements, which wrap around when the observed motion is larger than one-half of the radar wavelength. Thus, the spatio-temporal ''unwrapping" of phase observations is necessary to obtain physically meaningful results. Several different algorithms have been developed for time series analysis of InSAR data to solve for this ambiguity. These algorithms may employ different models for time series analysis, but they all generate a first-order deformation rate, which can be compared to each other. However, there is no single algorithm that can provide optimal results in all cases. Since time series analyses of InSAR data are used in a variety of applications with different characteristics, each algorithm possesses inherently unique strengths and weaknesses. In this review article, following a brief overview of InSAR technology, we discuss several algorithms developed for time series analysis of InSAR data using an example set of results for measuring subsidence rates in Mexico City.

InSAR↗

Salt Kinematics and InSAR

As part of a long-term attempt to learn how the climatic and tectonic signal interact to shape a steady state mountain monitored displacement of a markers in SE termination and also near the summit of a small viscous salt fountain extruding onto the Central plateau of Iran. The marker displacements relate to the first InSAR interferograms of salt extrusion (980913 to 990620) calculated Earth tides, winds, air pressures and temperatures. In the first documented staking exercise, hammered wooden stakes vertically through the surgical marl (c. 1 Ocm deep) onto the top of crystalline salt. These stakes installed in an irregular array elongate E-W along the c.50 m high cliff marking the effective SE terminus of the glacier at Qum Kuh(Centra1 Iran) ,just to the E of a NE trending river cliff about 40 m high. We merely measured the distances between pairs of stakes with known azimuth about 2 m apart to calculate sub horizontal strain in a small part of Qum Kuh. Stakes moved and micro strains for up to 46 pairs of stakes (p strain= ((lengthl-length2)/1engthl) x 10-1) was calculated for each seven stake epochs and plotted against their azimuth on simplified array maps. The data fit well the sine curves cxpected of the maximum and minimum strain ellipses. The first documented stakes located on the SE where the InSAR image show -1 1 to 0 mm pink to purple, 0 to lOmm purple to blue, and show high activity of salt in low activity area of the InSAR image (980913 to 990620).Short term micro strains of stake tie lines record anisotropic expansions due to heating and contraction due to cooling. All epochs changed between 7 to 1 17 days (990928 to000 1 16), showed 200 to 400 micro strain lengthening and shortening. The contraction and extension existed in each epoch, but the final strain was extension in E-W in Epoch land 6, contraction in E-W direction during epochs 2-3-4-5 and 7. The second pair of stakes hammered about 20 cm deep into the deep soils(more than 1 m) , near summit, where the colors change between 19 to 29mm in InSAR image(9809 13 to 990620). Additional information is included in the original abstract.

Aftabi, Pedarm↗

Formation Flying for Distributed InSAR

We consider two spacecraft flying in formation to create interferometric synthetic aperture radar (InSAR). Several candidate orbits for such in InSar formation have been previously determined based on radar performance and Keplerian orbital dynamics. However, with out active control, disturbance-induced drift can degrade radar performance and (in the worst case) cause a collision. This study evaluates the feasibility of operating the InSAR spacecraft as a formation, that is, with inner-spacecraft sensing and control. We describe the candidate InSAR orbits, design formation guidance and control architectures and algorithms, and report the (Delta)(nu) and control acceleration requirements for the candidate orbits for several tracking performance levels. As part of determining formation requirements, a formation guidance algorithm called Command Virtual Structure is introduced that can reduce the (Delta)(nu) requirements compared to standard Leader/Follower formation approaches.

control↗

Tropospheric Correction for InSAR Using Interpolated ECMWF Data and GPS Zenith Total Delay

To mitigate atmospheric errors caused by the troposphere, which is a limiting error source for spaceborne interferometric synthetic aperture radar (InSAR) imaging, a tropospheric correction method has been developed using data from the European Centre for Medium- Range Weather Forecasts (ECMWF) and the Global Positioning System (GPS). The ECMWF data was interpolated using a Stretched Boundary Layer Model (SBLM), and ground-based GPS estimates of the tropospheric delay from the Southern California Integrated GPS Network were interpolated using modified Gaussian and inverse distance weighted interpolations. The resulting Zenith Total Delay (ZTD) correction maps have been evaluated, both separately and using a combination of the two data sets, for three short-interval InSAR pairs from Envisat during 2006 on an area stretching from northeast from the Los Angeles basin towards Death Valley. Results show that the root mean square (rms) in the InSAR images was greatly reduced, meaning a significant reduction in the atmospheric noise of up to 32 percent. However, for some of the images, the rms increased and large errors remained after applying the tropospheric correction. The residuals showed a constant gradient over the area, suggesting that a remaining orbit error from Envisat was present. The orbit reprocessing in ROI_pac and the plane fitting both require that the only remaining error in the InSAR image be the orbit error. If this is not fulfilled, the correction can be made anyway, but it will be done using all remaining errors assuming them to be orbit errors. By correcting for tropospheric noise, the biggest error source is removed, and the orbit error becomes apparent and can be corrected for

Webb, Frank H.↗

InSAR-based detection method for mapping and monitoring slow-moving landslides in remote regions with steep and mountainous terrain: An application to Nepal

Mapping and monitoring landslides in remote areas with steep and mountainous terrain is logistically challenging, expensive, and time consuming. Yet, in order to mitigate hazards and prevent loss of life in these areas, and to better understand landslide processes, high-resolution measurements of landslide activity are necessary. Satellite-based synthetic aperture radar interferometry (InSAR) provides millimeter-scale measurements of ground surface deformation that can be used to identify and monitor landslides in remote areas where ground-based monitoring techniques are not feasible. Here we present a novel InSAR deformation detection approach, which uses double difference time-series with local and regional spatial filters and pixel clustering methods to identify and monitor slow-moving landslides without making a priori assumptions of the location of landslides. We apply our analysis to freely available Copernicus Sentinel-1 satellite data acquired between 2014 and 2017 centered on the Trishuli River drainage basin in Nepal. We found a minimum of 6 slow-moving landslides that all occur within the Ranimatta lithologic formation (phyllites, metasandstones, meta basics). These landslides have areas ranging from 0.39 to 1.66 sq.km and long-term dry-season displacement rates ranging from 2.1 to 8.8 cm/yr. Due to periods of low coherence during the monsoon season (June – September) each year, and following the 25 April 2015 Mw7.8 Gorkha earthquake, our time series analysis is limited to the 2014-2015 and 2016-2017 dry seasons (September - May). We found that each of the landslides displayed slightly higher rates during the 2014 period, likely as a result of higher cumulative rainfall that fell during the 2014 monsoon season. Although we do not have high quality InSAR data to show the landslide evolution directly following the Gorkha earthquake, the similar rates of movement before (2014-2015) and after (2016-2017) Gorkha suggest the earthquake had negligible long-term impact on these landslides. Our findings highlight the potential for region-wide mapping of slow-moving l and slides using freely available remote sensing data in remote areas such as Nepal and future work will benefit from expanding our methodology to other regions around the world.

Synthetic Aperture Radar↗

Simulation Studies of the Effect of Forest Spatial Structure on InSAR Signature

The height of scattering phase retrieved from InSAR data is considered being correlated with the tree height and the spatial structure of the forest stand. Though some researchers have used simple backscattering models to estimate tree height from the height of scattering center, the effect of forest spatial structure on InSAR data is not well understood yet. A three-dimensional coherent radar backscattering model for forest canopies based on realistic three-dimensional scene was used to investigate the effect in this paper. The realistic spatial structure of forest canopies was established either by field measurements (stem map) or through use of forest growth model. Field measurements or a forest growth model parameterized using local environmental parameters provides information of forest species composition and tree sizes in certain growth phases. A fractal tree model (L-system) was used to simulate individual 3- D tree structure of different ages or heights. Trees were positioned in a stand in certain patterns resulting in a 3-D medium of discrete scatterers. The radar coherent backscatter model took the 3-D forest scene as input and simulates the coherent radar backscattering signature. Interferometric SAR images of 3D scenes were simulated and heights of scattering phase centers were estimated from the simulated InSAR data. The effects of tree height, crown cover, crown depth, and the spatial distribution patterns of trees on the scattering phase center were analyzed. The results will be presented in the paper.

Sun, Guoqing↗

Software for Generating Troposphere Corrections for InSAR Using GPS and Weather Model Data

Atmospheric errors due to the troposphere are a limiting error source for spaceborne interferometric synthetic aperture radar (InSAR) imaging. This software generates tropospheric delay maps that can be used to correct atmospheric artifacts in InSAR data. The software automatically acquires all needed GPS (Global Positioning System), weather, and Digital Elevation Map data, and generates a tropospheric correction map using a novel algorithm for combining GPS and weather information while accounting for terrain. Existing JPL software was prototypical in nature, required a MATLAB license, required additional steps to acquire and ingest needed GPS and weather data, and did not account for topography in interpolation. Previous software did not achieve a level of automation suitable for integration in a Web portal. This software overcomes these issues. GPS estimates of tropospheric delay are a source of corrections that can be used to form correction maps to be applied to InSAR data, but the spacing of GPS stations is insufficient to remove short-wavelength tropospheric artifacts. This software combines interpolated GPS delay with weather model precipitable water vapor (PWV) and a digital elevation model to account for terrain, increasing the spatial resolution of the tropospheric correction maps and thus removing short wavelength tropospheric artifacts to a greater extent. It will be integrated into a Web portal request system, allowing use in a future L-band SAR Earth radar mission data system. This will be a significant contribution to its technology readiness, building on existing investments in in situ space geodetic networks, and improving timeliness, quality, and science value of the collected data

Moore, Angelyn W.↗

The InSAR Scientific Computing Environment 3.0: A Flexible Framework for NISAR Operational and User-Led Science Processing

The InSAR Scientific Computing Environment (ISCE) was first developed under the NASA Advanced Information Systems Technology as a flexible, extensible object-oriented framework for Interferometric Synthetic Aperture Radar (InSAR) processing. The ISCE framework uses Python 3 at the workflow level, controlling modules of compiled code for functional processing, and managing inputs, outputs, and other flow control services. The currently released version, called ISCE 2.1, is distributed to the research community through the Western North America InSAR Consortium under a research license. The ISCE team is working on the next generation of the code in order to prepare for the NASAISRO SAR (NISAR) mission operational processing. Innovations in this code include augmentation or conversion of the custom Python framework elements in ISCE with the Pyre framework, new workflows for interferometric and polarimetric stack processing, a more intuitive and graphically based user interface, and flow control for hybrid computing environments including CPU/GPU clusters, logging and error tracking facilities, and new more efficient computational modules that exploit graphical processor units (GPUs) when available. The ISCE 3.0 framework is designed to work in an operational environment as well as on a single user’s laptop or compute cluster, with services to discover capabilities and scale computations accordingly.

Buckley, Sean M.↗

Tracking Coastal Change in American Samoa By Mapping Local Vertical Land Motion with PS-InSAR

Characterizingdiverse contributionsto coastal land change is a key step to mitigating the effects ofrising sea levels that threaten coastal communities. This taskis particularly critical for small island communities in tectonically active regions, which are highlyvulnerableto the effects of sea level rise. We highlight here a case study to extract regional estimates of vertical land motion (VLM) over American Samoa, which in recent years has observedincreased nuisance flooding.Weusedpersistent scatterer Interferometric Synthetic Aperture Radar (PS-InSAR)to derive high-resolutionestimates of crustaldeformationinpopulated regions overTutuila Island from 2015 to 2021. While the area is small and highly vegetatedand poses challenges for InSAR, we wereable to construct a regional map of the estimated deformation ratein these areasandvalidate the time-series with alocal GPS station. Our preliminary results suggest that PS-InSAR has the ability to capture the local and regional deformation patterns, although further workis needed to compensate formore complexatmospheric effects, estimate error margins,and integrate ourresults intomodels of sea-level change that can be used local stakeholders

Stacey A Huang↗

Joint retrieval of soil moisture and permafrost active layer thickness using L-band InSAR and P-band PolSAR

Seasonal subsidence measured by repeat-pass interferometric synthetic aperture radar (InSAR) can be used to infer the active layer thickness (ALT) in permafrost regions. The differential volume of soil water undergoing the phase change over the thaw season is one of the factors impacting the seasonal subsidence and is a function of both soil moisture profile and thaw depth. Without the information about soil moisture, this InSAR approach can have large biases in the ALT estimates when soil moisture profile is below saturation. Soil moisture and ALT can also be estimated from polarimetric synthetic aperture radar (PolSAR) backscatter observations but the sensing depth of the PolSAR approach is limited when deep ALT is present. In this paper, we integrated these two approaches and applied a joint retrieval method to estimate the soil moisture profiles and ALT from the L-band InSAR and P-band PolSAR data acquired over the Arctic-boreal region during the 2017 Arctic-Boreal Vulnerability Experiment (ABoVE) airborne campaign.

Schaefer, Kevin M.↗

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↗

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↗

Mitigation of tropospheric InSAR phase artifacts through differential multisquint processing

We propose a technique for mitigating tropospheric phase errors in repeat-pass interferometric synthetic aperture radar (InSAR). The mitigation technique is based upon the acquisition of multisquint InSAR data. On each satellite pass over a target area, the radar instrument will acquire images from multiple squint (azimuth) angles, from which multiple interferograms can be formed. The diversity of viewing angles associated with the multisquint acquisition can be used to solve for two components of the 3-D surface displacement vector as well as for the differential tropospheric phase. We describe a model for the performance of the multisquint technique, and we present an assessment of the performance expected.

troposphere↗

InSAR Scientific Computing Environment

This computing environment is the next generation of geodetic image processing technology for repeat-pass Interferometric Synthetic Aperture (InSAR) sensors, identified by the community as a needed capability to provide flexibility and extensibility in reducing measurements from radar satellites and aircraft to new geophysical products. This software allows users of interferometric radar data the flexibility to process from Level 0 to Level 4 products using a variety of algorithms and for a range of available sensors. There are many radar satellites in orbit today delivering to the science community data of unprecedented quantity and quality, making possible large-scale studies in climate research, natural hazards, and the Earth's ecosystem. The proposed DESDynI mission, now under consideration by NASA for launch later in this decade, would provide time series and multiimage measurements that permit 4D models of Earth surface processes so that, for example, climate-induced changes over time would become apparent and quantifiable. This advanced data processing technology, applied to a global data set such as from the proposed DESDynI mission, enables a new class of analyses at time and spatial scales unavailable using current approaches. This software implements an accurate, extensible, and modular processing system designed to realize the full potential of InSAR data from future missions such as the proposed DESDynI, existing radar satellite data, as well as data from the NASA UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar), and other airborne platforms. The processing approach has been re-thought in order to enable multi-scene analysis by adding new algorithms and data interfaces, to permit user-reconfigurable operation and extensibility, and to capitalize on codes already developed by NASA and the science community. The framework incorporates modern programming methods based on recent research, including object-oriented scripts controlling legacy and new codes, abstraction and generalization of the data model for efficient manipulation of objects among modules, and well-designed module interfaces suitable for command- line execution or GUI-programming. The framework is designed to allow users contributions to promote maximum utility and sophistication of the code, creating an open-source community that could extend the framework into the indefinite future.

Rosen, Paul A.↗

Simplex GPS and InSAR Inversion Software

Changes in the shape of the Earth's surface can be routinely measured with precisions better than centimeters. Processes below the surface often drive these changes and as a result, investigators require models with inversion methods to characterize the sources. Simplex inverts any combination of GPS (global positioning system), UAVSAR (uninhabited aerial vehicle synthetic aperture radar), and InSAR (interferometric synthetic aperture radar) data simultaneously for elastic response from fault and fluid motions. It can be used to solve for multiple faults and parameters, all of which can be specified or allowed to vary. The software can be used to study long-term tectonic motions and the faults responsible for those motions, or can be used to invert for co-seismic slip from earthquakes. Solutions involving estimation of fault motion and changes in fluid reservoirs such as magma or water are possible. Any arbitrary number of faults or parameters can be considered. Simplex specifically solves for any of location, geometry, fault slip, and expansion/contraction of a single or multiple faults. It inverts GPS and InSAR data for elastic dislocations in a half-space. Slip parameters include strike slip, dip slip, and tensile dislocations. It includes a map interface for both setting up the models and viewing the results. Results, including faults, and observed, computed, and residual displacements, are output in text format, a map interface, and can be exported to KML. The software interfaces with the QuakeTables database allowing a user to select existing fault parameters or data. Simplex can be accessed through the QuakeSim portal graphical user interface or run from a UNIX command line.

Donnellan, Andrea↗

Central Valley Water Resources: Improving California Groundwater Assessments using GRACE and InSAR Datasets for Water Resource Management

California’s Central Valley is one of the most productive agricultural areas in the world, producing approximately $20 billion in crops annually. The recent California droughts of 2007-2010 and 2012-2019 resulted in increased groundwater pumping in the Central Valley to adequately irrigate farmland. Overdrafting of the Central Valley aquifer results in groundwater depletion, land subsidence, and permanent loss of groundwater storage. In 2014, depletion of groundwater led the state of California to enact the Sustainable Groundwater Management Act (SGMA),requiring critically overdrafted, high, and medium priority sub-basins to reach sustainable levels of groundwater pumping and recharge by 2042. SGMA allows local Groundwater Sustainability Agencies (GSAs) the authority to create Groundwater Sustainability Plans (GSPs) at the sub-basin level. To assist California’s Department of Water Resources (DWR), this project quantified groundwater change and land subsidence in Central Valley sub-basins with sparse or unreliable well and Geographic Positioning Systems (GPS) data. This was done using NASA’s Gravity Recovery and Climate Experiment (GRACE), GRACE FollowOn (GRACE-FO), and interferograms derived from Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Advanced Land Observing Satellite 2 (ALOS-2)Phased Array L-band Synthetic Aperture Radar 2 (PALSAR-2). Time series of the GRACE and InSAR data were compared with well and GPS data in data-dense sub-basins to determine the feasibility of these datasets for groundwater storage and subsidence monitoring. We found that GRACE and InSAR data are effective tools for determining groundwater change and land subsidence and can be used on their own to monitor sub-basins in the absence of well and GPS data

Water Resources↗

Central Valley Water Resources: Improving California Groundwater Assessments using GRACE and InSAR Datasets for Water Resource Management

California’s Central Valley is one of the most productive agricultural areas in the world, producing approximately $20 billion in crops annually. The recent California droughts of 2007-2010 and 2011-2017 resulted in increased groundwater pumping in the Central Valley to adequately irrigate farmland. Overdrafting of the Central Valley aquifer results in groundwater depletion, land subsidence, and permanent loss of groundwater storage. In 2014, depletion of groundwater led the state of California to enact the Sustainable Groundwater Management Act (SGMA), requiring critically overdrafted, high, and medium priority sub-basins to reach sustainable levels of groundwater pumping and recharge by 2042. SGMA allows local Groundwater Sustainability Agencies the authority to create Groundwater Sustainability Plans at the sub-basin level. To assist California’s Department of Water Resources, this project quantified groundwater change and land subsidence in Central Valley sub-basins with sparse or unreliable well and GPS data. This was done using NASA’s Gravity Recovery and Climate Experiment (GRACE), GRACE Follow-On (GRACE-FO), and interferograms derived from Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Advanced Land Observing Satellite 2 (ALOS-2) Phased Array L-band Synthetic Aperture Radar 2 (PALSAR-2). Time series of the GRACE and InSAR data were compared with well and GPS data in data-dense sub-basins to determine the feasibility of these datasets for groundwater storage and subsidence monitoring. We found thatGRACE and InSAR data are effective tools for determining groundwater change and land subsidence and can be used on their own to monitor sub-basins in the absence of well and GPS data.

Water Resources↗