Adaptation of Range-Doppler Algorithm for Efficient Beamforming of Monostatic and Multistatic Ultrasound Signals.

Algorithmic changes that increase beamforming speed have become increasingly relevant to processing synthetic aperture (SA) ultrasound data. In particular, beamforming SA data in a spatio-temporal frequency domain using the F-k (Stolt) migration have been shown to reduce the beamforming time by up to two orders of magnitude compared with the conventional delay-and-sum (DAS) beamforming, and it has been used in applications where large amounts of raw data make real-time frame rates difficult to attain, such as multistatic SA imaging and plane-wave Doppler imaging with large ensemble lengths. However, beamforming signals in a spatio-temporal Fourier space can require loading large blocks of data at once, making it memory-intensive and less suited for parallel (i.e., multithreaded) processing. As an alternative, we propose beamforming in a range-Doppler (RD) frequency domain using the range-Doppler algorithm (RDA) that has originally been developed for SA radar (SAR) imaging. Through simulation and phantom experiments, we show that RDA achieves similar lateral resolution and contrast compared with DAS and F-k migration. At the same time, higher axial sidelobes in RDA images can be reduced via (temporal) frequency binning. Like the F-k migration, RDA significantly reduces the overall number of computations relative to DAS, and it achieves ten times lower processing time on a single CPU. Because RDA uses only a spatial Fourier transform (FT), it requires two times less memory than the F-k migration to process the simulated multistatic data and can be executed on as many as a thousand parallel threads (compared with eight parallel threads for the F-k migration), making it more suitable for implementation on modern graphics processing units (GPUs). While RDA is not as parallelizable as DAS, it is expected to hold a significant speed advantage on devices with moderate parallel processing capabilities (up to several thousand cores), such as point-of-care and low-cost ultrasound devices.

Permafrost Dynamics Observatory-Part I: Postprocessing and Calibration Methods of UAVSAR L-Band InSAR Data for Seasonal Subsidence Estimation.

Interferometric synthetic aperture radar (InSAR) has been used to quantify a range of surface and near surface physical properties in permafrost landscapes. Most previous InSAR studies have utilized spaceborne InSAR platforms, but InSAR datasets over permafrost landscapes collected from airborne platforms have been steadily growing in recent years. Most existing algorithms dedicated toward retrieval of permafrost physical properties were originally developed for spaceborne InSAR platforms. In this study, which is the first in a two part series, we introduce a series of calibration techniques developed to apply a novel joint retrieval algorithm for permafrost active layer thickness retrieval to an airborne InSAR dataset acquired in 2017 by NASA's Uninhabited Aerial Vehicle Synthetic Aperture Radar over Alaska and Western Canada. We demonstrate how InSAR measurement uncertainties are mitigated by these calibration methods and quantify remaining measurement uncertainties with a novel method of modeling interferometric phase uncertainty using a Gaussian mixture model. Finally, we discuss the impact of native SAR resolution on InSAR measurements, the limitation of using few interferograms per retrieval, and the implications of our findings for cross-comparison of airborne and spaceborne InSAR datasets acquired over Arctic regions underlain by permafrost.

Global drainage patterns and the origins of topographic relief on Earth, Mars, and Titan.

Rivers have eroded the topography of Mars, Titan, and Earth, creating diverse landscapes. However, the dominant processes that generated topography on Titan (and to some extent on early Mars) are not well known. We analyzed drainage patterns on all three bodies and found that large drainages, which record interactions between deformation and erosional modification, conform much better to long-wavelength topography on Titan and Mars than on Earth. We use a numerical landscape evolution model to demonstrate that short-wavelength deformation causes drainage directions to diverge from long-wavelength topography, as observed on Earth. We attribute the observed differences to ancient long-wavelength topography on Mars, recent or ongoing generation of long-wavelength relief on Titan, and the creation of short-wavelength relief by plate tectonics on Earth.

Release of arsenic to deep groundwater in the Mekong Delta, Vietnam, linked to pumping-induced land subsidence.

Deep aquifers in South and Southeast Asia are increasingly exploited as presumed sources of pathogen- and arsenic-free water, although little is known of the processes that may compromise their long-term viability. We analyze a large area (>1,000 km(2)) of the Mekong Delta, Vietnam, in which arsenic is found pervasively in deep, Pliocene-Miocene-age aquifers, where nearly 900 wells at depths of 200-500 m are contaminated. There, intensive groundwater extraction is causing land subsidence of up to 3 cm/y as measured using satellite-based radar images from 2007 to 2010 and consistent with transient 3D aquifer simulations showing similar subsidence rates and total subsidence of up to 27 cm since 1988. We propose a previously unrecognized mechanism in which deep groundwater extraction is causing interbedded clays to compact and expel water containing dissolved arsenic or arsenic-mobilizing solutes (e.g., dissolved organic carbon and competing ions) to deep aquifers over decades. The implication for the broader Mekong Delta region, and potentially others like it across Asia, is that deep, untreated groundwater will not necessarily remain a safe source of drinking water.

Edgelist phase unwrapping algorithm for time series InSAR analysis.

We present here a new integer programming formulation for phase unwrapping of multidimensional data. Phase unwrapping is a key problem in many coherent imaging systems, including time series synthetic aperture radar interferometry (InSAR), with two spatial and one temporal data dimensions. The minimum cost flow (MCF) [IEEE Trans. Geosci. Remote Sens. 36, 813 (1998)] phase unwrapping algorithm describes a global cost minimization problem involving flow between phase residues computed over closed loops. Here we replace closed loops by reliable edges as the basic construct, thus leading to the name "edgelist." Our algorithm has several advantages over current methods-it simplifies the representation of multidimensional phase unwrapping, it incorporates data from external sources, such as GPS, where available to better constrain the unwrapped solution, and it treats regularly sampled or sparsely sampled data alike. It thus is particularly applicable to time series InSAR, where data are often irregularly spaced in time and individual interferograms can be corrupted with large decorrelated regions. We show that, similar to the MCF network problem, the edgelist formulation also exhibits total unimodularity, which enables us to solve the integer program problem by using efficient linear programming tools. We apply our method to a persistent scatterer-InSAR data set from the creeping section of the Central San Andreas Fault and find that the average creep rate of 22 mm/Yr is constant within 3 mm/Yr over 1992-2004 but varies systematically with ground location, with a slightly higher rate in 1992-1998 than in 1999-2003.

Size and shape of Saturn's moon Titan.

Cassini observations show that Saturn's moon Titan is slightly oblate. A fourth-order spherical harmonic expansion yields north polar, south polar, and mean equatorial radii of 2574.32 +/- 0.05 kilometers (km), 2574.36 +/- 0.03 km, and 2574.91 +/- 0.11 km, respectively; its mean radius is 2574.73 +/- 0.09 km. Titan's shape approximates a hydrostatic, synchronously rotating triaxial ellipsoid but is best fit by such a body orbiting closer to Saturn than Titan presently does. Titan's lack of high relief implies that most--but not all--of the surface features observed with the Cassini imaging subsystem and synthetic aperture radar are uncorrelated with topography and elevation. Titan's depressed polar radii suggest that a constant geopotential hydrocarbon table could explain the confinement of the hydrocarbon lakes to high latitudes.

Phase unwrapping in three dimensions with application to InSAR time series.

The problem of phase unwrapping in two dimensions has been studied extensively in the past two decades, but the three-dimensional (3D) problem has so far received relatively little attention. We develop here a theoretical framework for 3D phase unwrapping and also describe two algorithms for implementation, both of which can be applied to synthetic aperture radar interferometry (InSAR) time series. We test the algorithms on simulated data and find both give more accurate results than a two-dimensional algorithm. When applied to actual InSAR time series, we find good agreement both between the algorithms and with ground truth.