Two-dimensional (2-D) phase unwrapping is a key step in the analysis of interferometric synthetic aperture radar (InSAR) data. While challenging even in the best of circumstances, this problem poses ...unique difficulties when the dimensions of the interferometric input data exceed the limits of one's computational capabilities. In order to deal with such cases, we propose a technique for applying the statistical-cost, network-flow phase-unwrapping algorithm (SNAPHU) of Chen and Zebker (2001) to large datasets. Specifically, we introduce a methodology whereby a large interferogram is partitioned into a set of several smaller tiles that are unwrapped individually and then divided further into independent, irregularly shaped reliable regions. These regions are subsequently assembled into a full unwrapped solution, with the phase offsets between regions computed in a secondary optimization problem whose objective is to maximize the a posteriori probability of the final solution. As this secondary problem assumes the same statistical models as employed in the initial tile-unwrapping stage, the technique results in a solution that approximates the solution that would have been obtained had the full-size interferogram been unwrapped as a single piece. The secondary problem is framed in terms of network-flow ideas, allowing the use of an existing nonlinear solver. Applying the algorithm to a large topographic interferogram acquired over central Alaska, we find that the technique is less prone to unwrapping artifacts than more simple tiling approaches.
In the San Joaquin Valley, California, recent droughts starting in 2007 have increased the pumping of groundwater, leading to widespread subsidence. In the southern portion of the San Joaquin Valley, ...vertical subsidence as high as 85 cm has been observed between June 2007 and December 2010 using Interferometric Synthetic Aperture Radar (InSAR). This study seeks to map regions where inelastic (not recoverable) deformation occurred during the study period, resulting in permanent compaction and loss of groundwater storage. We estimated the amount of permanent compaction by incorporating multiple data sets: the total deformation derived from InSAR, estimated skeletal‐specific storage and hydraulic parameters, geologic information, and measured water levels during our study period. We used two approaches, one that we consider to provide an estimate of the lowest possible amount of inelastic deformation, and one that provides a more reasonable estimate. These two approaches resulted in a spatial distribution of values for the percentage of the total deformation that was inelastic, with the former estimating a spatially averaged value of 54%, and the latter a spatially averaged value of 98%. The former corresponds to the permanent loss of
4.14 × 108 m3 of groundwater storage, or roughly 5% of the volume of groundwater used over the study time period; the latter corresponds to the loss of
7.48 × 108 m3 of groundwater storage, or roughly 9% of the volume of groundwater used. This study demonstrates that a data‐driven approach can be used effectively to estimate the permanent loss of groundwater storage.
Plain Language Summary
Subsidence due to groundwater pumping from 2007 to 2010 in the San Joaquin Valley, California was mapped using satellite data known as InSAR. These data were incorporated with additional datasets, including geological information, to estimate how much subsidence was permanent. This subsidence represents a permanent loss in groundwater storage. Using these methods, we estimated that a permanent loss of 7.48 × 108 m3 of groundwater storage occurred during our study period. This accounted for roughly 9% of groundwater pumping in our study area. While this is just a small fraction of the total water stored underground in this area, this 'water of compaction' is an important safeguard for times of drought that, once removed, cannot be replenished.
Key Points
Groundwater withdrawal during the 2007–2010 drought in the San Joaquin Valley, California caused significant subsidence
Geomechanical properties, geologic data, and water level data were integrated with InSAR to estimate inelastic deformation
Results indicate that the majority of deformation that occurred over this time period was inelastic
Rockslides have a high socioeconomic and environmental importance in many countries. Norway is particularly susceptible to large rockslides due to its many fjords and steep mountains. One of the most ...dangerous hazards related with rock slope failures are tsunamis that can lead to large loss of life. It is therefore very important to systematically identify potential unstable rock slopes. Traditional landslide monitoring techniques are expensive and time consuming. Differential satellite interferometric synthetic aperture radar (InSAR) is an invaluable tool for land displacement monitoring. Improved access to time series of satellite data has led to the development of several innovative multitemporal algorithms. Small baseline (SB) methods are based on combining and inverting a set of unwrapped interferograms that are computed with a small perpendicular baseline in order to reduce spatial phase decorrelation. Another well proven technique is the persistent scatterer interferometric method (PSI) that is based on analysis of persistent point targets. In this paper, we apply both approaches to study several rockslide sites in Troms County in the far north of Norway. Moreover, we take the opportunity to address the difference and similarities between the SB and the PSI multitemporal InSAR methods for displacement studies in rural terrain.
Radar interferograms provide precise (millimeter level) measurements of crustal deformation at fine resolution over wide areas. The dominant error source in many interferograms results from spatial ...heterogeneity of the wet component of atmospheric refractivity, resulting in excess path length of the radar signal propagating through the neutral atmosphere. In this report, we introduce a method to compensate for atmospheric phase artifacts in a radar interferogram using spatially interpolated zenith wet delay (ZWD) data obtained from a network of Global Positioning System (GPS) receivers in the region imaged by the radar. Our method, based on Taylor's “frozen‐flow” hypothesis, uses GPS data recorded prior to and after the individual radar acquisition times to infer denser spatial networks of control points for interpolation than is suggested by the sparse and irregular spatial distribution of GPS receivers. We test this approach by correcting phase fluctuations observed in a radar interferogram over an area covered by the Southern California Integrated GPS Network stations, using maps of atmospheric delay interpolated from the denser network of GPS ZWD measurements generated by our method. We find that the largest improvement results from removing an altitude‐dependent model derived from the GPS data: the uncorrected interferogram in our example has an RMS fluctuation of 16.1 mm, which falls to 8.6 mm after the altitude correction. A map derived from ZWD measurements recorded only at the instances of radar acquisition reduced overall fluctuations in the data from 8.6 to 8.1 mm RMS, while using the inferred denser network reduces the error further to 7.5 mm RMS. Most of this residual variation is due to decorrelation of the signal and does not result from atmospheric propagation, however. After smoothing to remove decorrelation noise, the residual drops from 4.7 to 3.9 mm (GPS stations only) and then to 2.7 mm RMS when the denser network algorithm is used.
Despite the abundance of rock glaciers in the Sierra Nevada of California, USA, few efforts have been made to measure their surface flow. Here we use the interferometric synthetic aperture radar ...(InSAR) technique to compile a benchmark inventory describing the kinematic state of 59 active rock glaciers in this region. In the late summer of 2007, these rock glaciers moved at speeds that range from 14 cm yr−1 to 87 cm yr−1, with a regional mean value of 53 cm yr−1. Our inventory reveals a spatial difference: rock glaciers in the southern Sierra Nevada moved faster than the ones in the central Sierra Nevada. In addition to the regional mapping, we also conduct a case study to measure the surface flow of the Mount Gibbs rock glacier in fine spatial and temporal detail. The InSAR measurements over this target reveal (1) that the spatial pattern of flow is correlated with surface geomorphic features and (2) a significant seasonal variation of flow speed whose peak value was 48 cm yr−1in the fall of 2007, more than twice the minimum value observed in the spring of 2008. The seasonal variation lagged air temperatures by three months. Our finding on the seasonal variation of surface speed reinforces the importance of a long time series with high temporal sampling rates to detect possible long-term changes of rock glacier kinematics in a warming climate.
Cassini RADAR altimetry data collected on the 49th flyby of Titan (2008 December 21) over Ontario Lacus in Titan's south polar region provides strong evidence for an extremely smooth surface, with ...less than 3 mm rms surface height variation over the 100m‐wide Fresnel zone. Histograms of the raw radar echoes imply a mirror‐like specular reflection of the transmitted signal. Such an echo is possible only if the surface is extremely flat relative to our 2.2‐cm wavelength. The 3 mm upper bound follows from analyzing the strength of the specular return, which declines exponentially with increasing surface height variance. In this experiment, the strength of the echo was larger than expected, severely saturating the receiver. We developed a method to partially correct the echoes for the distortion incurred. While the implied mm‐scale smoothness is not proof that the surface is liquid, it is unlikely that a solid surface is so smooth.
Drained thermokarst lake basins (DTLBs) are ubiquitous landforms on Arctic tundra lowland. Their dynamic states are seldom investigated, despite their importance for landscape stability, hydrology, ...nutrient fluxes, and carbon cycling. Here we report results based on high-resolution Interferometric Synthetic Aperture Radar (InSAR) measurements using space-borne data for a study area located on the North Slope of Alaska near Prudhoe Bay, where we focus on the seasonal thaw settlement within DTLBs, averaged between 2006 and 2010. The majority (14) of the 18 DTLBs in the study area exhibited seasonal thaw settlement of 3-4 cm. However, four of the DTLBs examined exceeded 4 cm of thaw settlement, with one basin experiencing up to 12 cm. Combining the InSAR observations with the in situ active layer thickness measured using ground penetrating radar and mechanical probing, we calculated thaw strain, an index of thaw settlement strength along a transect across the basin that underwent large thaw settlement. We found thaw strains of 10-35% at the basin center, suggesting the seasonal melting of ground ice as a possible mechanism for the large settlement. These findings emphasize the dynamic nature of permafrost landforms, demonstrate the capability of the InSAR technique to remotely monitor surface deformation of individual DTLBs, and illustrate the combination of ground-based and remote sensing observations to estimate thaw strain. Our study highlights the need for better description of the spatial heterogeneity of landscape-scale processes for regional assessment of surface dynamics on Arctic coastal lowlands.
A radar interferometric technique for topographic mapping of surfaces, implemented utilizing a single synthetic aperture radar (SAR) system in a nearly repeating orbit, is discussed. The authors ...characterize the various sources contributing to the echo correlation statistics, and isolate the term which most closely describes surficial change. They then examine the application of this approach to topographic mapping of vegetated surfaces which may be expected to possess varying backscatter over time. It is found that there is decorrelation increasing with time but that digital terrain model generation remains feasible. The authors present such a map of a forested area in Oregon which also includes some nearly unvegetated lava flows. Such a technique could provide a global digital terrain map.< >
Cassini RADAR images of Titan’s south polar region acquired during southern summer contain lake features which disappear between observations. These features show a tenfold increases in backscatter ...cross-section between images acquired one year apart, which is inconsistent with common scattering models without invoking temporal variability. The morphologic boundaries are transient, further supporting changes in lake level. These observations are consistent with the exposure of diffusely scattering lakebeds that were previously hidden by an attenuating liquid medium. We use a two-layer model to explain backscatter variations and estimate a drop in liquid depth of approximately 1-m-per-year. On larger scales, we observe shoreline recession between ISS and RADAR images of Ontario Lacus, the largest lake in Titan’s south polar region. The recession, occurring between June 2005 and July 2009, is inversely proportional to slopes estimated from altimetric profiles and the exponential decay of near-shore backscatter, consistent with a uniform reduction of 4
±
1.3
m in lake depth.
Of the potential explanations for observed surface changes, we favor evaporation and infiltration. The disappearance of dark features and the recession of Ontario’s shoreline represents volatile transport in an active methane-based hydrologic cycle. Observed loss rates are compared and shown to be consistent with available global circulation models. To date, no unambiguous changes in lake level have been observed between repeat images in the north polar region, although further investigation is warranted. These observations constrain volatile flux rates in Titan’s hydrologic system and demonstrate that the surface plays an active role in its evolution. Constraining these seasonal changes represents the first step toward our understanding of longer climate cycles that may determine liquid distribution on Titan over orbital time periods.
Interferogram images derived from repeat‐pass spaceborne synthetic aperture radar systems exhibit artifacts due to the time and space variations of atmospheric water vapor. Other tropospheric ...variations, such as pressure and temperature, also induce distortions, but the effects are smaller in magnitude and more evenly distributed throughout the interferogram than the wet troposphere term. Spatial and temporal changes of 20% in relative humidity lead to 10 cm errors in deformation products, and perhaps 100 m of error in derived topographic maps for those pass pairs with unfavorable baseline geometries. In wet regions such as Hawaii, these are by far the dominant errors in the Spaceborne Imaging Radar‐C and X Band Synthetic Aperature Radar (SIR‐C/X‐SAR) interferometric products. The unknown time delay from tropospheric distortion is independent of frequency, and thus multiwavelength measurements, such as those commonly used to correct radar altimeter and Global Positioning System (GPS) ionospheric biases, cannot be used to rectify the error. In the topographic case, the errors may be mitigated by choosing interferometric pairs with relatively long baselines, as the error amplitude is inversely proportional to the perpendicular component of the interferometer baseline. For the SIR‐C/X‐SAR Hawaii data we found that the best (longest) baseline pair produced a map supporting 100 m contouring, whereas the poorest baseline choice yielded an extremely noisy topographic map even at this coarse contour interval. In the case of deformation map errors the result is either independent of baseline parameters or else very nearly so. Here the only solution is averaging of independent interferograms, so in order to create accurate deformation products in wet regions many multiple passes may be required. Rules for designing optimal data acquisition and processing sequences for interferometric analyses in nondesert parts of the world are (1) to use the longest radar wavelengths possible, within ionospheric scintillation and Faraday rotation limits, (2) for topography, maximize interferometer baseline within decorrelation limits* and (3) for surface deformation, use multiple observations and average the derived products. Following the above recipe yields accuracies of 10 m for digital elevation models and 1 cm for deformation maps even in very wet regions, such as Hawaii.