Ocean Remote Sensing with Synthetic Aperture Radar

The ocean covers approximately 71% of the Earth’s surface, 90% of the biosphere and contains 97% of Earth’s water. The Synthetic Aperture Radar (SAR) can image the ocean surface in all weather conditions and day or night. SAR remote sensing on ocean and coastal monitoring has become a research hotspot in geoscience and remote sensing. This book—Progress in SAR Oceanography—provides an update of the current state of the science on ocean remote sensing with SAR. Overall, the book presents a variety of marine applications, such as, oceanic surface and internal waves, wind, bathymetry, oil spill, coastline and intertidal zone classification, ship and other man-made objects’ detection, as well as remotely sensed data assimilation. The book is aimed at a wide audience, ranging from graduate students, university teachers and working scientists to policy makers and managers. Efforts have been made to highlight general principles as well as the state-of-the-art technologies in the field of SAR Oceanography

Design Options For Low Cost, Low Power Microsatellite Based SAR.

This research aims at providing a system design that reduces the mass and cost of spaceborne Synthetic Aperture Radar (SAR) missions by a factor of two compared to current (TecSAR - 300 kg, ~ £ 127 M) or planned (NovaSAR-S — 400 kg, ~ £ 50 M) mission. This would enable the cost of a SAR constellation to approach that of the current optical constellation such as Disaster Monitoring Constellation (DMC). This research has identified that the mission cost can be reduced significantly by: focusing on a narrow range of applications (forestry and disasters monitoring); ensuring the final design has a compact stowage volume, which facilitates a shared launch; and building the payload around available platforms, rather than the platform around the payload. The central idea of the research has been to operate the SAR at a low instantaneous power level—a practical proposition for a micro-satellite based SAR. The use of a simple parabolic reflector with a single horn at L-band means that a single, reliable and efficient Solid State Power Amplifier (SSPA) can be used to lower the overall system cost, and to minimise the impact on the spacecraft power system. A detailed analysis of basic pulsed (~ 5 - 10 % duty cycle) and Continuous Wave (CW) SAR (100 % duty cycle) payloads has shown their inability to fit directly into existing microsatellite buses without involving major changes, or employing more than one platform. To circumvent the problems of pulsed and CW techniques, two approaches have been formulated. The first shows that a CW SAR can be implemented in a mono-static way with a single antenna on a single platform. In this technique, the SAR works in an Interrupted CW (ICW) mode, but these interruptions introduce periodic gaps in the raw data. On processing, these gapped data result in artefacts in the reconstructed images. By applying data based statistical estimation techniques to “fill in the gaps” in the simulated raw SAR data, this research has shown the possibility of minimising the effects of these artefacts. However, once the same techniques are applied to the real SAR data (in this case derived from RADARSAT-1), the artefacts are shown to be problematic. Because of this the ICW SAR design technique it is—set aside. The second shows that an extended chirp mode pulsed (ECMP) SAR (~ 20 - 54 % duty cycle) can be designed with a lowered peak power level which enables a single SSPA to feed a parabolic Cassegrain antenna. The detailed analysis shows the feasibility of developing a microsatellite based SAR design at a comparable price to those of optical missions

Arrayed synthetic aperture radar

In this thesis, the use of array processing techniques applied to Single Input Multiple Output (SIMO) SAR systems with enhanced capabilities is investigated. In Single Input Single Output (SISO) SAR systems there is a high resolution, wide swath contradiction, whereby it is not possible to increase both cross-range resolution and the imaged swath width simultaneously. To overcome this, a novel beamformer for SAR systems in the cross-range direction is proposed. In particular, this beamformer is a superresolution beamformer capable of forming wide nulls using subspace based approaches. SIMO SAR systems also give rise to additional sets of received data, which includes geometrical information about the SAR and target environment, and can be used for enhanced target parameter estimation. In particular, this thesis looks at round trip delay, joint azimuth and elevation angle, and relative target power estimation. For round trip delay estimation, the use of the traditional matched filter with subspace partitioning is proposed. Then by using a joint 2D Multiple Signal Classification (MUSIC) algorithm, joint Direction of Arrival (DOA) estimation can be achieved. Both the use of range lines of raw SAR data and the use of a Region of Interest (ROI) of a SAR image are investigated. However in terms of imaging, MUSIC is not well-suited for SAR, due to its target response not corresponding to the target's true power return. Therefore a joint DOA and target power estimation algorithm is proposed to overcome this limitation. These algorithms provide the framework for the development of three processing techniques. These allow sidelobe suppression in the slant range direction, along with the reconstruction of undersampled data and region enhancement using MUSIC with power preservation.Open Acces

Processing of multiple-receiver spaceborne arrays for wide-area SAR

The instantaneous area illuminated by a single-aperture synthetic aperture radar (SAR) is fundamentally limited by the minimum SAR antenna area constraint. This limitation is due to the fact that the number of illuminated resolution cells cannot exceed the number of collected data samples. However, if spatial sampling is added through the use of multiple-receiver arrays, then the maximum unambiguous illumination area is increased because multiple beams can be formed to reject range-Doppler ambiguities. Furthermore, the maximum unambiguous illumination area increases with the number of receivers in the array. One spaceborne implementation of multiple-aperture SAR that has been proposed is a constellation of formation-flying satellites. In this implementation, several satellites fly in a cluster and work together as a single coherent system. There are many advantages to the constellation implementation including cost benefits, graceful performance degradation, and the possibility of performing in multiple modes. The disadvantage is that the spatial samples provided by such a constellation will be sparse and irregularly spaced; consequently, traditional matched filtering produces unsatisfactory results. We investigate SAR performance and processing of sparse, multiple-aperture arrays. Three filters are evaluated: the matched filter, maximum-likelihood filter, and minimum mean-squared error filter. It is shown that the maximum-likelihood and minimum mean-squared error filters can provide quality SAR images when operating on data obtained from sparse satellite constellations. We also investigate the performance of the three filters versus system parameters such as SNR, the number of receivers in the constellation, and satellite positioning error

Learning to Estimate Sea Ice Concentration from SAR Imagery

Through the growing interest in the Arctic for shipping, mining and climate research, large-scale high quality ice concentration is of great interest. Due to the unavailability of suitable ice concentration estimation algorithms, ice concentration maps are interpreted from synthetic aperture radar (SAR) images manually by ice experts for operational uses. An automatic ice concentration estimation algorithm is required for accurate large-scale ice mapping. In this thesis, a set of algorithms are developed aiming at operational ice concentration estimation from SAR images. The major difficulty in designing a robust algorithm for ice concentration estimation from SAR images is the constantly changing SAR image features of ice and water in time and location. This difficulty is addressed by learning features instead of designing features from SAR images. A set of convolutional neural network based ice concentration estima- tion algorithms are developed to learn multi-scale SAR image features and simultaneously regress ice concentration from the learned image features. We first demonstrated the capa- bility of CNNs in ice concentration estimation from SAR images when trained using image analysis charts as ground truth. Then the model is further improved by taking into account the errors in the image analysis charts. Ice concentration estimates with improved robust- ness to training samples errors, accuracy and scale of details are obtained. The robustness of the developed methods are further demonstrated in the melt season of the Beaufort Sea, where reasonable ice concentration estimates are acquired. In order to reduce the model training time, it is desired to reuse existing models. The model transferability is evaluated and suggestions on using existing models to accelerate the training process are given, which is shown to reduce the training time by over 10 times in our case

Compaction of C-band synthetic aperture radar based sea ice information for navigation in the Baltic Sea

In this work operational sea ice synthetic aperture radar (SAR) data products were improved and developed. A SAR instrument is transmitting electromagnetic radiation at certain wavelengths and measures the radiation which is scattered back towards the instrument from the target, in our case sea and sea ice. The measured backscattering is converted to an image describing the target area through complex signal processing. The images, however, differ from optical images, i.e. photographs, and their visual interpretation is not straightforward. The main idea in this work has been to deliver the essential SAR-based sea ice information to end-users (typically on ships) in a compact and user-friendly format. The operational systems at Finnish Institute of Marine Research (FIMR) are currently based on the data received from a Canadian SAR-satellite, Radarsat-1. The operational sea ice classification, developed by the author with colleagues, has been further developed. One problem with the SAR data is typically that the backscattering varies depending on the incidence angle. The incidence angle is the angle in which the transmitted electromagnetic wave meets the target surface and it varies within each SAR image and between different SAR images depending on the measuring geometry. To improve this situation, an incidence angle correction algorithm to normalize the backscattering over the SAR incidence angle range for Baltic Sea ice has been developed as part of this work. The algorithm is based on SAR backscattering statistics over the Baltic Sea. To locate different sea ice areas in SAR images, a SAR segmentation algorithm based on pulse-coupled neural networks has been developed and tested. The parameters have been tuned suitable for the operational data in use at FIMR. The sea ice classification is based on this segmentation and the classification is segment-wise rather than pixel-wise. To improve SAR-based distinguishing between sea ice and open water an open water detection algorithm based on segmentation and local autocorrelation has been developed. Also ice type classification based on higher-order statistics and independent component analysis have been studied to get an improved SAR-based ice type classification. A compression algorithm for compressing sea ice SAR data for visual use has been developed. This algorithm is based on the wavelet decomposition, zero-tree structure and arithmetic coding. Also some properties of the human visual system were utilized. This algorithm was developed to produce smaller compressed SAR images, with a reasonable visual quality. The transmission of the compressed images to ships with low-speed data connections in reasonable time is then possible. One of the navigationally most important sea ice parameters is the ice thickness. SAR-based ice thickness estimation has been developed and evaluated as part of this work. This ice thickness estimation method uses the ice thickness history derived from digitized ice charts, made daily at the Finnish Ice Service, as its input, and updates this chart based on the novel SAR data. The result is an ice thickness chart representing the ice situation at the SAR acquisition time in higher resolution than in the manually made ice thickness charts. For the evaluation of the results a helicopter-borne ice thickness measuring instrument, based on electromagnetic induction and laser altimeter, was used.reviewe

SAR imaging of moving targets by subaperture based low-rank and sparse decomposition

Synthetic aperture radar (SAR) has gained significance as an indispensable instrument of remote sensing and airborne surveillance. Its applications extend to 3D terrain mapping, oil spill detection, crop yield estimation and disaster evaluation. SAR utilizes platform motion to synthesize a large antenna thus rendering a very fine spatial resolution. Nevertheless, imaging of moving targets with SAR is a challenging problem. In this thesis, we propose a moving target imaging approach for SAR which exploits the low-rank and sparse decomposition (LRSD) of the subaperture data. As a first step, multiple subapertures are constructed from the raw data using frequency domain filtering. In contrast to the stationary points, moving targets in the SAR scene shift their position in the various subapertures. This enables a successful low-rank and sparse decomposition of the subaperture data where the sparse component captures the moving targets’ phase histories and reflectivity profiles. On the other hand, the low-rank component consists of the static background due to fewer spatial variations in multiple subapertures. This framework allows the reconstruction of full-resolution sparse and low-rank images by combining the spectral information of the decomposed subapertures. Furthermore, it enhances the applicability of sparsity-driven moving target imaging frameworks to very low signal to clutter ratio (SCR) scenarios by offering a considerable SCR performance improvement. We manifest the effectiveness of our approach through experiments with synthetic as well as real SAR data. Our real SAR experiments were based on MiniSAR and EMISAR data

Applications of Satellite Earth Observations section - NEODAAS: Providing satellite data for efficient research

The NERC Earth Observation Data Acquisition and Analysis Service (NEODAAS) provides a central point of Earth Observation (EO) satellite data access and expertise for UK researchers. The service is tailored to individual users’ requirements to ensure that researchers can focus effort on their science, rather than struggling with correct use of unfamiliar satellite data

Satellite monitoring of harmful algal blooms (HABs) to protect the aquaculture industry

Harmful algal blooms (HABs) can cause sudden and considerable losses to fish farms, for example 500,000 salmon during one bloom in Shetland, and also present a threat to human health. Early warning allows the industry to take protective measures. PML's satellite monitoring of HABs is now funded by the Scottish aquaculture industry. The service involves processing EO ocean colour data from NASA and ESA in near-real time, and applying novel techniques for discriminating certain harmful blooms from harmless algae. Within the AQUA-USERS project we are extending this capability to further HAB species within several European countries

Retrieval of Ocean Surface Currents and Winds Using Satellite SAR backscatter and Doppler frequency shift

Ocean surface winds and currents play an important role for weather, climate, marine life, ship navigation, oil spill drift and search and rescue. In-situ observations of the ocean are sparse and costly. Satellites provide a useful complement to these observations. Synthetic aperture radar (SAR) is particularly attractive due to its high spatial resolution and its capability to extract both sea surface winds and currents day and night and almost independent of weather.The work in this thesis involves processing of along-track interferometric SAR (ATI-SAR) data, analysis of the backscatter and Doppler frequency shift, and development of wind and current retrieval algorithms. Analysis of the Doppler frequency shift showed a systematic bias. A calibration method was proposed and implemented to correct for this bias. Doppler analysis also showed that the wave contribution to the SAR Doppler centroid often dominates over the current contribution. This wave contribution is estimated using existing theoretical and empirical Doppler models. For wind and current retrieval, two methods were developed and implemented.The first method, called the direct method, consists of retrieval of the wind speed from SAR backscatter using an empirical backscatter model. In order to retrieve the radial current, the retrieved wind speed is used to correct for the wave contribution. The current retrieval was assessed using two different (theoretical and empirical) Doppler models and wind inputs (model and SAR-derived). It was found that the results obtained by combining the Doppler empirical model with the SAR-derived wind speed were more consistent with ocean models.The second method, called Bayesian method, consists of blending the SAR observables (backscatter and Doppler shift) with an atmospheric and an oceanic model to retrieve the total wind and current vector fields. It was shown that this method yields more accurate estimates, i.e. reduces the models biases against in-situ measurements. Moreover, the method introduces small scale features, e.g. fronts and meandering, which are weakly resolved by the models.The correlation between the surface wind vectors and the SAR Doppler shift was demonstrated empirically using the Doppler shift estimated from over 300 TanDEM-X interferograms and ECMWF reanalysis wind vectors. Analysis of polarimetric data showed that theoretical models such as Bragg and composite surface models over-estimate the backscatter polarization ratio and Doppler shift polarization difference. A combination of a theoretical Doppler model and an empirical modulation transfer function was proposed. It was found that this model is more consistent with the analyzed data than the pure theoretical models.The results of this thesis will be useful for integrating SAR retrievals in ocean current products and assimilating SAR observables in the atmospheric, oceanic or coupled models. The results are also relevant for preparation studies of future satellite missions