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

Selection of the key earth observation sensors and platforms focusing on applications for Polar Regions in the scope of Copernicus system 2020-2030

An optimal payload selection conducted in the frame of the H2020 ONION project (id 687490) is presented based on the ability to cover the observation needs of the Copernicus system in the time period 2020–2030. Payload selection is constrained by the variables that can be measured, the power consumption, and weight of the instrument, and the required accuracy and spatial resolution (horizontal or vertical). It involved 20 measurements with observation gaps according to the user requirements that were detected in the top 10 use cases in the scope of Copernicus space infrastructure, 9 potential applied technologies, and 39 available commercial platforms. Additional Earth Observation (EO) infrastructures are proposed to reduce measurements gaps, based on a weighting system that assigned high relevance for measurements associated to Marine for Weather Forecast over Polar Regions. This study concludes with a rank and mapping of the potential technologies and the suitable commercial platforms to cover most of the requirements of the top ten use cases, analyzing the Marine for Weather Forecast, Sea Ice Monitoring, Fishing Pressure, and Agriculture and Forestry: Hydric stress as the priority use cases.Peer ReviewedPostprint (published version

Publications of the Jet Propulsion Laboratory 1989

This bibliography describes and indexes by primary author the externally distributed technical reporting, released during 1989, that resulted from scientific and engineering work performed, or managed, by JPL. Three classes of publications are included: JPL publications in which the information is complete for a specific accomplishment; articles from the quarterly Telecommunications and Data Acquisition (TDA) Progress Report; and articles published in the open literature

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications

Geosynchronous synthetic aperture radar for Earth continuous observation missions

This thesis belongs to the field of remote sensing, particularly Synthetic Aperture Radar (SAR) systems from the space. These systems acquire the signals along the orbital track of one or more satellites where the transmitter and receiver are mounted, and coherently process the echoes in order to form the synthetic aperture. So, high resolution images can be obtained without using large arrays of antennas. The study presented in this thesis is centred in a novel concept in SAR, which is known as Geosynchronous SAR or GEOSAR, where the transmitter and/or receiver are placed in a platform in a geostationary orbit. In this case, the small relative motions between the satellite and the Earth surface are taken to get the necessary motion to form the synthetic aperture and focus the image. The main advantage of these systems with respect to the current technology (where LEO satellites with lower height are considered) is the possibility of permanently acquire images from the same region thanks to the small motion of the platform. Therefore, the different possibilities in the orbital design that offer this novel technology as well as the geometric resolutions obtained in the final image have been firstly studied. However, the use of geosynchronous satellites as illuminators results in slant ranges between 35.000-38.000 Km, which are much higher than the typical values obtained in LEOSAR, under 1.000 Km. Fortunately, the slow motion of the satellite makes possible large integration of pulses during minutes or even hours, reaching Signal-to-Noise Ratio (SNR) levels in the order of LEO acquisitions without using high transmitted power or large antennas. Moreover, such large integration times, increases the length of the synthetic aperture to get the desired geometric resolutions of the image (in the order of a few meters or kilometres depending on the application). On the other hand, the use of long integration time presents some drawbacks such as the scene targets decorrelation, atmospheric artefacts due to the refraction index variations in the tropospheric layer, transmitter and receiver clock jitter, clutter decorrelation or orbital positioning errors; which will affect the correct focusing of the image. For this reason, a detailed theoretical study is presented in the thesis in order to characterize and model these artefacts. Several simulations have been performed in order to see their effects on the final images. Some techniques and algorithms to track and remove these errors from the focused image are presented and the improvement of the final focused image is analysed. Additionally, the real data from a GB-SAR (Ground-Based SAR) have been reused to simulate a long integration time acquisition and see the effects in the image focusing as well as to check the performance of compensation algorithms in the final image. Finally, a ground receiver to reuse signals of opportunity from a broadcasting satellite have been designed and manufactured. This hardware is expected to be an important tool for experimental testing in future GEOSAR analysis.Aquesta tesi s'emmarca dins de l'àmbit de la teledetecció, en particular, en els sistemes coneguts com a radar d'obertura sintètica (SAR en anglès) des de l'espai. Aquests sistemes adquireixen senyal al llarg de l'òrbita d'un o més satèl·lits on estan situats el transmissor i el receptor, i processa els ecos de forma coherent per a formar l'obertura sintètica. D'aquesta manera es poden aconseguir imatge d'alta resolució sense la necessitat d'emprar un array d'antenes molt gran. El treball realitzat en aquest estudi es centra en un nou concepte dins del món SAR que consisteix en l'ús de satèl·lits en òrbita geostacionària per a l'adquisició d'imatges, sistemes coneguts com a Geosynchronous SAR o GEOSAR. En aquest cas, els petits moviments relatius dels satèl·lits respecte de la superfície terrestre s'empren per a aconseguir el desplaçament necessari per a formar l'obertura sintètica i així obtenir la imatge. El principal avantatge d'aquests sistemes respecte a la tecnologia actual (on s'utilitzen satèl·lits en orbites més baixes LEO) és la possibilitat d'adquirir imatges d'una mateixa zona de forma permanent gràcies als petits desplaçaments del satèl·lit. Així doncs, en aquesta tesi s'estudien les diferents possibilitats en el disseny orbital que ofereixen aquests sistemes així com les resolucions d'imatge que s'obtindrien. Tot i així, l'ús de satèl·lits en òrbita geoestacionària, resulta en una distància entre el transmissor/receptor i l'escena entre 35000-38000 Km, molt més gran que les distàncies típiques en els sistemes LEO per sota dels 1000 Km. Tot i així, el moviment lent de les plataformes geostacionàries fa possible la integració de polsos durant minuts o hores, arribant a nivells acceptables de relació senyal a soroll (SNR) sense necessitat d'utilitzar potències transmeses i antenes massa grans. A més a més, aquesta llarga integració també permet assolir unes longituds d'obertura sintètica adients per a arribar a resolucions d'imatge desitjades (de l'ordre de pocs metres o kilòmetres segons l'aplicació). Malgrat això, l'ús de temps d'integració llargs té una sèrie d'inconvenients com poden ser la decorrelació dels blancs de l'escena, l'aparició d'artefactes atmosfèrics deguts als canvis d'índex de refracció en la troposfera, derives dels rellotges del transmissor i receptor, decorrelació del clutter o errors en el posicionament orbital, que poden afectar la correcta focalització de la imatge. Així doncs, en la tesi s'ha fet un detallat estudi teòric d'aquests problemes per tal de modelitzar-los i posteriorment s'han realitzat diverses simulacions per veure els seus efectes en una imatge. Diverses tècniques per a compensar aquests errors i millorar la qualitat de la imatge també s'han estudiat al llarg de la tesi. Per altra banda, dades reals d'un GB-SAR (SAR en una base terrestre) s'han reutilitzat per adaptar-les a una possible adquisició de llarga durada i veure així de forma experimental com afecta la llarga integració en les imatges i com millora l'enfocament després d'aplicar els algoritmes de compensació. Per últim, en la tesi es presenta un sistema receptor terrestre per a poder realitzar un anàlisi experimental del cas GEOSAR utilitzant un il·luminador d'oportunitat. Els primers passos en el disseny i la fabricació del hardware també es presenten en aquesta tes

Estimation of Forest Biomass and Faraday Rotation using Ultra High Frequency Synthetic Aperture Radar

Synthetic Aperture Radar (SAR) data in the Ultra High Frequency (UHF; 300 MHz – 3 GHz)) band have been shown to be strongly dependent of forest biomass, which is a poorly estimated variable in the global carbon cycle. In this thesis UHF-band SAR data from the fairly flat hemiboreal test site Remningstorp in southern Sweden were analysed. The data were collected on several occasions with different moisture conditions during the spring of 2007. Regression models for biomass estimation on stand level (0.5-9 ha) were developed for each date on which SAR data were acquired. For L-band (centre frequency 1.3 GHz) the best estimation model was based on HV-polarized backscatter, giving a root mean squared error (rmse) between 31% and 46% of the mean biomass. For P-band (centre frequency 340 MHz), regression models including HH, HV or HH and HV backscatter gave an rmse between 18% and 27%. Little or no saturation effects were observed up to 290 t/ha for P-band. A model based on physical-optics has been developed and was used to predict HH-polarized SAR data with frequencies from 20 MHz to 500 MHz from a set of vertical trunks standing on an undulating ground surface. The model shows that ground topography is a critical issue in SAR imaging for these frequencies. A regression model for biomass estimation which includes a correction for ground slope was developed using multi-polarized P-band SAR data from Remningstorp as well as from the boreal test site Krycklan in northern Sweden. The latter test site has pronounced topographic variability. It was shown that the model was able to partly compensate for moisture variability, and that the model gave an rmse of 22-33% when trained using data from Krycklan and evaluated using data from Remningstorp. Regression modelling based on P-band backscatter was also used to estimate biomass change using data acquired in Remningstorp during the spring 2007 and during the fall 2010. The results show that biomass change can be measured with an rmse of about 15% or 20 tons/ha. This suggests that not only deforestation, but also forest growth and degradation (e.g. thinning) can be measured using P-band SAR data. The thesis also includes result on Faraday rotation, which is an ionospheric effect which can have a significant impact on spaceborne UHF-band SAR images. Faraday rotation angles are estimated in spaceborne L-band SAR data. Estimates based on distributed targets and calibration targets with high signal to clutter ratios are found to be in very good agreement. Moreover, a strong correlation with independent measurements of Total Electron Content is found, further validating the estimates

Satellite remote sensing of surface winds, waves, and currents: Where are we now?

This review paper reports on the state-of-the-art concerning observations of surface winds, waves, and currents from space and their use for scientific research and subsequent applications. The development of observations of sea state parameters from space dates back to the 1970s, with a significant increase in the number and diversity of space missions since the 1990s. Sensors used to monitor the sea-state parameters from space are mainly based on microwave techniques. They are either specifically designed to monitor surface parameters or are used for their abilities to provide opportunistic measurements complementary to their primary purpose. The principles on which is based on the estimation of the sea surface parameters are first described, including the performance and limitations of each method. Numerous examples and references on the use of these observations for scientific and operational applications are then given. The richness and diversity of these applications are linked to the importance of knowledge of the sea state in many fields. Firstly, surface wind, waves, and currents are significant factors influencing exchanges at the air/sea interface, impacting oceanic and atmospheric boundary layers, contributing to sea level rise at the coasts, and interacting with the sea-ice formation or destruction in the polar zones. Secondly, ocean surface currents combined with wind- and wave- induced drift contribute to the transport of heat, salt, and pollutants. Waves and surface currents also impact sediment transport and erosion in coastal areas. For operational applications, observations of surface parameters are necessary on the one hand to constrain the numerical solutions of predictive models (numerical wave, oceanic, or atmospheric models), and on the other hand to validate their results. In turn, these predictive models are used to guarantee safe, efficient, and successful offshore operations, including the commercial shipping and energy sector, as well as tourism and coastal activities. Long-time series of global sea-state observations are also becoming increasingly important to analyze the impact of climate change on our environment. All these aspects are recalled in the article, relating to both historical and contemporary activities in these fields