Microwave Radiometry at Frequencies From 500 to 1400 MHz: An Emerging Technology for Earth Observations.

Microwave radiometry has provided valuable spaceborne observations of Earth's geophysical properties for decades. The recent SMOS, Aquarius, and SMAP satellites have demonstrated the value of measurements at 1400 MHz for observing surface soil moisture, sea surface salinity, sea ice thickness, soil freeze/thaw state, and other geophysical variables. However, the information obtained is limited by penetration through the subsurface at 1400 MHz and by a reduced sensitivity to surface salinity in cold or wind-roughened waters. Recent airborne experiments have shown the potential of brightness temperature measurements from 500-1400 MHz to address these limitations by enabling sensing of soil moisture and sea ice thickness to greater depths, sensing of temperature deep within ice sheets, improved sensing of sea salinity in cold waters, and enhanced sensitivity to soil moisture under vegetation canopies. However, the absence of significant spectrum reserved for passive microwave measurements in the 500-1400 MHz band requires both an opportunistic sensing strategy and systems for reducing the impact of radio-frequency interference. Here, we summarize the potential advantages and applications of 500-1400 MHz microwave radiometry for Earth observation and review recent experiments and demonstrations of these concepts. We also describe the remaining questions and challenges to be addressed in advancing to future spaceborne operation of this technology along with recommendations for future research activities.

Change-point detection for recursive Bayesian geoacoustic inversions.

In order to carry out geoacoustic inversion in low signal-to-noise ratio (SNR) conditions, extended duration observations coupled with source and/or receiver motion may be necessary. As a result, change in the underlying model parameters due to time or space is anticipated. In this paper, an inversion method is proposed for cases when the model parameters change abruptly or slowly. A model parameter change-point detection method is developed to detect the change in the model parameters using the importance samples and corresponding weights that are already available from the recursive Bayesian inversion. If the model parameters change abruptly, a change-point will be detected and the inversion will restart with the pulse measurement after the change-point. If the model parameters change gradually, the inversion (based on constant model parameters) may proceed until the accumulated model parameter mismatch is significant and triggers the detection of a change-point. These change-point detections form the heuristics for controlling the coherent integration time in recursive Bayesian inversion. The method is demonstrated in simulation with parameters corresponding to the low SNR, 100-900 Hz linear frequency modulation pulses observed in the Shallow Water 2006 experiment [Tan, Gerstoft, Yardim, and Hodgkiss, J. Acoust. Soc. Am. 136, 1187-1198 (2014)].

Recursive Bayesian synthetic aperture geoacoustic inversion in the presence of motion dynamics.

A low signal to noise ratio (SNR), single source/receiver, broadband, frequency-coherent matched-field inversion procedure recently has been proposed. It exploits coherently repeated transmissions to improve estimation of the geoacoustic parameters. The long observation time improves the SNR and creates a synthetic aperture due to relative source-receiver motion. To model constant velocity source/receiver horizontal motion, waveguide Doppler theory for normal modes is necessary. However, the inversion performance degrades when source/receiver acceleration exists. Furthermore processing a train of pulses all at once does not take advantage of the natural incremental acquisition of data along with the ability to assess the temporal evolution of parameter uncertainty. Here a recursive Bayesian estimation approach is developed that coherently processes the data pulse by pulse and incrementally updates estimates of parameter uncertainty. It also approximates source/receiver acceleration by assuming piecewise constant but linearly changing source/receiver velocities. When the source/receiver acceleration exists, it is shown that modeling acceleration can reduce further the parameter estimation biases and uncertainties. The method is demonstrated in simulation and in the analysis of low SNR, 100-900 Hz linear frequency modulated (LFM) pulses from the Shallow Water 2006 experiment.

Compressive geoacoustic inversion using ambient noise.

Surface generated ambient noise can be used to infer sediment properties. Here, a passive geoacoustic inversion method that uses noise recorded by a drifting vertical array is adopted. The array is steered using beamforming to compute the noise arriving at the array from various directions. This information is used in two different ways: Coherently (cross-correlation of upward/downward propagating noise using a minimum variance distortionless response fathometer), and incoherently (bottom loss vs frequency and angle using a conventional beamformer) to obtain the bottom properties. Compressive sensing is used to invert for the number of sediment layer interfaces and their depths using coherent passive fathometry. Then the incoherent bottom loss estimate is used to refine the sediment thickness, sound speed, density, and attenuation values. Compressive sensing fathometry enables automatic determination of the number of interfaces. It also tightens the sediment thickness priors for the incoherent bottom loss inversion which reduces the search space. The method is demonstrated on drifting array data collected during the Boundary 2003 experiment.

Particle smoothers in sequential geoacoustic inversion.

Sequential Bayesian methods such as particle filters have been used to track a moving source in an unknown and space/time-evolving ocean environment. These methods treat both the source and the ocean parameters as non-stationary unknown random variables and track them via the multivariate posterior probability density function. Particle filters are numerical methods that can operate on nonlinear systems with non-Gaussian probability density functions. Particle smoothers are a natural extension to these filters. A smoother is appropriate in applications where data before and after the time of interest are readily available. Both past and "future" measurements are exploited in smoothers, whereas filters just use past measurements. Geoacoustic and source tracking is performed here using two smoother algorithms, the forward-backward smoother and the two-filter smoother. Smoothing is demonstrated on experimental data from both the SWellEx-96 and SW06 experiments where the parameter uncertainty is reduced relative to just filtering alone.

Broadband synthetic aperture geoacoustic inversion.

A typical geoacoustic inversion procedure involves powerful source transmissions received on a large-aperture receiver array. A more practical approach is to use a single moving source and/or receiver in a low signal to noise ratio (SNR) setting. This paper uses single-receiver, broadband, frequency coherent matched-field inversion and exploits coherently repeated transmissions to improve estimation of the geoacoustic parameters. The long observation time creates a synthetic aperture due to relative source-receiver motion. This approach is illustrated by studying the transmission of multiple linear frequency modulated (LFM) pulses which results in a multi-tonal comb spectrum that is Doppler sensitive. To correlate well with the measured field across a receiver trajectory and to incorporate transmission from a source trajectory, waveguide Doppler and normal mode theory is applied. The method is demonstrated with low SNR, 100-900 Hz LFM pulse data from the Shallow Water 2006 experiment.

Using Fisher information to quantify uncertainty in environmental parameters estimated from correlated ambient noise.

Efforts to characterize environmental parameters from ambient noise must contend with uncertainty introduced by stochastic fluctuations of the noise itself. This Letter calculates the Fisher information and Cramer-Rao bound of an unbiased correlated ambient noise parameter estimate. As an illustration, lower bounds on the error covariance of medium speed and attenuation parameters are obtained for a two-dimensional isotropic ambient noise scenario. The results demonstrate that an optimal sensor separation exists for obtaining the minimum error and the predictions are validated using simulated parameter inversions. The influences of record length, bandwidth, signal-to-noise, and spatial resolution are discussed.

Sequential geoacoustic inversion at the continental shelfbreak.

Environmental parameters can have large spatial and temporal variability in shelfbreak regions. The capability of sequential Bayesian filters in tracking this variation is investigated. Particle filtering (PF) is used to extract the environmental parameters and their uncertainties. The method tracks the environment with fewer particles relative to conventional geoacoustic inversion methods using successive independent inversions. As an example, data from the Shallow Water 2006 Experiment are processed. The PF approach first is used to track the source and the environment with little spatial variation just northwest of the shelfbreak. Then the strongly range-dependent shelfbreak region is analyzed and the PF results are compared to previous geoacoustic inversion studies from the region.

Particle filtering for passive fathometer tracking.

Seabed interface depths and fathometer amplitudes are tracked for an unknown and changing number of sub-bottom reflectors. This is achieved by incorporating conventional and adaptive fathometer processors into sequential Monte Carlo methods for a moving vertical line array. Sediment layering information and time-varying fathometer response amplitudes are tracked by using a multiple model particle filter with an uncertain number of reflectors. Results are compared to a classical particle filter where the number of reflectors is considered to be known. Reflector tracking is demonstrated for both conventional and adaptive processing applied to the drifting array data from the Boundary 2003 experiment. The layering information is successfully tracked by the multiple model particle filter even for noisy fathometer outputs.

Range aliasing in frequency coherent geoacoustic inversion.

This paper discusses the effects of frequency selection on source localization and geoacoustic inversion methods that use frequency coherent objective functions. Matched-field processors based on frequency-coherent objective functions often have rapidly fluctuating range ambiguity surfaces. Insufficient sampling in frequency domain results in range aliasing terms that affect geoacoustic inversion. Range aliasing and its effects on source localization and environmental parameter inversion are demonstrated on data collected during the MAPEX2000 experiment. Guidance for frequency selection to avoid range aliasing is provided.

Geoacoustic and source tracking using particle filtering: experimental results.

A particle filtering (PF) approach is presented for performing sequential geoacoustic inversion of a complex ocean acoustic environment using a moving acoustic source. This approach treats both the environmental parameters [e.g., water column sound speed profile (SSP), water depth, sediment and bottom parameters] at the source location and the source parameters (e.g., source depth, range and speed) as unknown random variables that evolve as the source moves. This allows real-time updating of the environment and accurate tracking of the moving source. As a sequential Monte Carlo technique that operates on nonlinear systems with non-Gaussian probability densities, the PF is an ideal algorithm to perform tracking of environmental and source parameters, and their uncertainties via the evolving posterior probability densities. The approach is demonstrated on both simulated data in a shallow water environment with a sloping bottom and experimental data collected during the SWellEx-96 experiment.

Tracking of geoacoustic parameters using Kalman and particle filters.

This paper incorporates tracking techniques such as the extended Kalman, unscented Kalman, and particle (PF) filters into geoacoustic inversion problems. This enables spatial and temporal tracking of environmental parameters and their underlying probability densities, making geoacoustic tracking a natural extension to geoacoustic inversion techniques. Water column and seabed properties are tracked in simulation for both vertical (VLA) and horizontal (HLA) line arrays using the three tracking filters. Filter performances are compared in terms of filter efficiencies using the posterior Cramer-Rao lower bound. Tracking capabilities of the geoacoustic filters under slowly and quickly changing environments are studied in terms of divergence statistics. Geoacoustic tracking can provide continuously environmental estimates and their uncertainties using only a fraction of the computational power of classical geoacoustic inversion schemes. Interfilter comparison show that while a high-particle-number PF outperforms the Kalman filters, there are many cases where all three filters perform equally well depending on the inversion configuration (such as the HLA versus VLA and frequency) and the tracked parameters.