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.

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.

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.