Processing methodology for the ITS_LIVE Sentinel-1 ice velocity product

The NASA MEaSUREs Inter-mission Time Series of Land Ice Velocity and Elevation (ITS_LIVE) project seeks to accelerate understanding of critical glaciers and ice sheet processes by providing researchers with global, low-latency, comprehensive and state-of-the-art records of surface velocities and elevations as observed from space. Here we describe the image-pair ice velocity product and processing methodology for ESA Sentinel-1 radar data. We demonstrate improvements to the core processing algorithm for dense offset tracking, “autoRIFT”, that provides finer resolution and higher accuracy data products with improved computational efficiency when compared to earlier versions. A novel calibration is applied to the data to correct for Sentinel-1A/B subswath- and full swath-dependent geolocation errors caused by systematic issues with the instruments. Sentinel-1’s C-band images are affected by variations in the total electron content of the ionosphere that results in large velocity errors in the azimuth (along-track) direction. To reduce these effects slant-range (line-of-sight or LOS) velocities are used and accompanied by LOS parameters that support map coordinate (x/y) velocity inversion from ascending and descending slant-range offset measurements, as derived from 2 image-pairs. The described product and methods comprise the MEaSUREs ITS_LIVE Sentinel-1 Image-Pair Glacier and Ice Sheet Surface Velocities: Version 2 (https://its-live.jpl.nasa.gov).

Vector Current Measurement Using Doppler Scatterometry with Optimally Selected Observation Azimuths

The Doppler scatterometer is a new style of remote sensing tool that can provide current measurements over a wide swath for rapid global coverage. The existing current estimation method for Doppler scatterometry uses the maximum likelihood method to jointly derive the wind and current fields but shows high computational complexity. Moreover, the current radial speeds measured along two arbitrary observation azimuths are used to derive the vector current according to the parallelogram rule, which is not applicable for the case where two observation azimuths are not perpendicular. In this paper, a vector current velocity inversion method using an optimally selected observation azimuth combination—as well as a general current velocity calculation method—is proposed for Doppler scatterometry. Firstly, current radial speeds along several different observation azimuths are estimated using an interferometric phase difference matching method with low computational complexity. Then, two current radial components of each point are arbitrarily selected to estimate a preliminary current direction using the proposed vector current velocity derivation method. Finally, two observation azimuths that have the smallest intersection angles with the preliminarily estimated current direction are selected for vector current velocity determination. With the Ocean Surface Current Analyses Real-time (OSCAR) data as current input, vector current estimation experiments were conducted based on simulation analysis using an instrument conceptual design model for a pencil-beam scatterometer. The results show that the standard deviation of the estimated current velocity magnitude is 0.06 m/s. Compared with the reported results obtained by the existing method, the inversion accuracy of velocity magnitude is improved by 67%.

Parsimonious velocity inversion applied to the Los Angeles Basin, CA

The proliferation of dense arrays promises to improve our ability to image geological structures at the scales necessary for accurate assessment of seismic hazard. However, combining the resulting local high-resolution tomography with existing regional models presents an ongoing challenge. We developed a framework based on the level-set method that provides a means to infer where local data provides meaningful constraints beyond those found in regional models - e.g. the Community Velocity Models (CVMs) of southern California. This technique defines a volume within which updates are made to a reference CVM, with the boundary of the volume being part of the inversion rather than explicitly defined. By penalizing the complexity of the boundary, a minimal update that sufficiently explains the data is achieved. To test this framework, we use data from the Community Seismic Network, a dense permanent urban deployment. We inverted Love wave dispersion and amplification data, from the Mw 6.4 and 7.1 2019 Ridgecrest earthquakes. We invert for an update to CVM-S4.26 using the Tikhonov Ensemble Sampling scheme, a highly efficient derivative-free approximate Bayesian method. We find the data is best explained by a deepening of the Los Angeles Basin with its deepest part south of downtown Los Angeles, along with a steeper northeastern basin wall. This result offers new progress towards the parsimonious incorporation of detailed local basin models within regional reference models utilizing an objective framework and highlights the importance of accurate basin models when accounting for the amplification of surface waves in the high-rise building response band.

Scholte wave inversion and passive source imaging with ocean-bottom DAS

Geotechnical characterization of marine sediments remains an outstanding challenge for offshore energy development, including foundation design and site selection of wind turbines and offshore platforms. We demonstrate that passive distributed acoustic sensing (DAS) surveys offer a new solution for shallow offshore geotechnical investigation where seafloor power or communications cables with fiber-optic links are available. We analyze Scholte waves recorded by DAS on a 42 km power cable in the Belgian offshore area of the southern North Sea. Ambient noise crosscorrelations converge acceptably with just over one hour of data, permitting multimodal Scholte wave dispersion measurement and shear-wave velocity inversion along the cable. We identify anomalous off-axis Scholte wave arrivals in noise crosscorrelations at high frequencies. Using a simple passive source imaging approach, we associate these arrivals with individual wind turbines, which suggests they are generated by structural vibrations. While many technological barriers must be overcome before ocean-bottom DAS can be applied to global seismic monitoring in the deep oceans, high-frequency passive surveys for high-resolution geotechnical characterization and monitoring in coastal regions are easily achievable today.

Accurate Thickness Inversion of Corrosion Using A1 Lamb Wave

Corrosion is one of the most critical issues in the oil and gas industry, leading to severe environmental and economic problems. Due to the high cost and serious safety risk of corrosion, it is essential to improve current corrosion testing techniques to detect corrosion damages at an early stage. Guided wave tomography (GWT) demonstrates its great potential to inspect and quantify the corrosion damage. GWT is capable of determining the residual life of corrosion structures by quantifying the remaining wall thickness. In this paper, an accurate guided wave tomography technique incorporating full waveform inversion (FWI) and higher-order Lamb waves (A1 mode) is presented for plate-like structures, which is able to get high-resolution reconstruction results. The technique consists of three steps: forward modeling, velocity inversion and thickness reconstruction. The forward modeling is computed by solving the elastic full-wave equations in 2-D time domain by using the finite difference method. High-resolution phase velocity inversion can then be obtained by minimizing the waveform misfit function between simulated and recorded data using a second order optimization method, which updates the velocity model from low to high frequencies iteratively. Finally, the velocity variations can be transformed into depth profiles by using the dispersive characteristics of ultrasonic guided waves; therefore, the thickness reconstruction can be obtained. The numerical simulations are performed on an aluminum plate with a complicated corrosion defect. By comparing the thickness reconstruction maps using both A1 and A0 modes, the results demonstrate that FWI with A1 mode can achieve significantly better resolution of corrosion imaging than that with A0 mode.