Mass balances of different sectors in Antarctica in 1996, 2000 and 2006

Large uncertainties remain in the current and future contribution to sea level rise from Antarctica. Climate warming may increase snowfall in the continent's interior, but enhance glacier discharge at the coast where warmer air and ocean temperatures erode the buttressing ice shelves. Here, we use satellite interferometric synthetic-aperture radar observations from 1992 to 2006 covering 85% of Antarctica's coastline to estimate the total mass flux into the ocean. We compare the mass fluxes from large drainage basin units with interior snow accumulation calculated from a regional atmospheric climate model for 1980 to 2004. In East Antarctica, small glacier losses in Wilkes Land and glacier gains at the mouths of the Filchner and Ross ice shelves combine to a near-zero loss of 4 ± 61 Gt/yr. In West Antarctica, widespread losses along the Bellingshausen and Amundsen seas increased the ice sheet loss by 59% in 10 years to reach 132 ± 60 Gt/yr in 2006. In the Peninsula, losses increased by 140% to reach 60 ± 46 Gt/yr in 2006. Losses are concentrated along narrow channels occupied by outlet glaciers and are caused by ongoing and past glacier acceleration. Changes in glacier flow therefore have a significant, if not dominant impact on ice sheet mass balance.

Wetlands in central Siberia from ASAR WS 2003, link to shapefiles

An area in central Siberia (partial coverage of Turukhansky und Yeniseysky districts) was investigated using satellite data. It covers freshwater ecosystems of non-forested peatlands in boreal forests. The satellite data represent the growing seasons of 2003/2004. Microwave data were acquired by the Advanced Synthetic Aperture Radar (ASAR) instrument onboard ENVISAT. The multi-temporal capabilities and resolution (150mx150m in WS mode) of the ASAR wide swath mode enabled the detection of dynamic features >2ha over this vast area. Scatterometer (QuikScat) data could be employed to distinguish hydro-periods. Wetland types have been identified on the basis of seasonal changes in backscatter. Results for peatlands have been compared with Russian forest inventory data which contain information on wetland distribution.

Design and implementation of a non-destructive defect detection technique based on UWB-SAR imaging

In the last twenty years aerospace and automotive industries started working widely with composite materials, which are not easy to test using classic Non-Destructive Inspection (NDI) techniques. Pairwise, the development of safety regulations sets higher and higher standards for the qualification and certification of those materials. In this thesis a new concept of a Non-Destructive defect detection technique is proposed, based on Ultrawide-Band (UWB) Synthetic Aperture Radar (SAR) imaging. Similar SAR methods are yet applied either in minefield [22] and head stroke [14] detection. Moreover feasibility studies have already demonstrated the validity of defect detection by means of UWB radars [12, 13]. The system was designed using a cheap commercial off-the-shelf radar device by Novelda and several tests of the developed system have been performed both on metallic specimen (aluminum plate) and on composite coupon (carbon fiber). The obtained results confirm the feasibility of the method and highlight the good performance of the developed system considered the radar resolution. In particular, the system is capable of discerning healthy coupons from damaged ones, and correctly reconstruct the reflectivity image of the tested defects, namely a 8 x 8 mm square bulge and a 5 mm drilled holes on metal specimen and a 5 mm drilled hole on composite coupon.

Tile-Based Two-Dimensional Phase Unwrapping for Digital Holography Using a Modular Framework

A variety of physical and biomedical imaging techniques, such as digital holography, interferometric synthetic aperture radar (InSAR), or magnetic resonance imaging (MRI) enable measurement of the phase of a physical quantity additionally to its amplitude. However, the phase can commonly only be measured modulo 2π, as a so called wrapped phase map. Phase unwrapping is the process of obtaining the underlying physical phase map from the wrapped phase. Tile-based phase unwrapping algorithms operate by first tessellating the phase map, then unwrapping individual tiles, and finally merging them to a continuous phase map. They can be implemented computationally efficiently and are robust to noise. However, they are prone to failure in the presence of phase residues or erroneous unwraps of single tiles. We tried to overcome these shortcomings by creating novel tile unwrapping and merging algorithms as well as creating a framework that allows to combine them in modular fashion. To increase the robustness of the tile unwrapping step, we implemented a model-based algorithm that makes efficient use of linear algebra to unwrap individual tiles. Furthermore, we adapted an established pixel-based unwrapping algorithm to create a quality guided tile merger. These original algorithms as well as previously existing ones were implemented in a modular phase unwrapping C++ framework. By examining different combinations of unwrapping and merging algorithms we compared our method to existing approaches. We could show that the appropriate choice of unwrapping and merging algorithms can significantly improve the unwrapped result in the presence of phase residues and noise. Beyond that, our modular framework allows for efficient design and test of new tile-based phase unwrapping algorithms. The software developed in this study is freely available.

Analyse des signaux radars polarimétriques en bandes C et L pour le suivi de l'humidité du sol de sites forestiers

Résumé : Dans les couverts forestiers, le suivi de l’humidité du sol permet de prévenir plusieurs désastres tels que la paludification, les incendies et les inondations. Comme ce paramètre est très dynamique dans l’espace et dans le temps, son estimation à grande échelle présente un grand défi, d’où le recours à la télédétection radar. Le capteur radar à synthèse d’ouverture (RSO) est couramment utilisé grâce à sa vaste couverture et sa résolution spatiale élevée. Contrairement aux sols nus et aux zones agricoles, le suivi de l’humidité du sol en zone forestière est très peu étudié à cause de la complexité des processus de diffusion dans ce type de milieu. En effet, la forte atténuation de la contribution du sol par la végétation et la forte contribution de volume issue de la végétation réduisent énormément la sensibilité du signal radar à l’humidité du sol. Des études portées sur des couverts forestiers ont montré que le signal radar en bande C provient principalement de la couche supérieure et sature vite avec la densité de la végétation. Cependant, très peu d’études ont exploré le potentiel des paramètres polarimétriques, dérivés d’un capteur polarimétrique comme RADARSAT-2, pour suivre l’humidité du sol sur les couverts forestiers. L’effet du couvert végétal est moins important avec la bande L en raison de son importante profondeur de pénétration qui permet de mieux informer sur l’humidité du sol. L’objectif principal de ce projet est de suivre l’humidité du sol à partir de données radar entièrement polarimétriques en bandes C et L sur des sites forestiers. Les données utilisées sont celles de la campagne terrain Soil Moisture Active Passive Validation EXperiment 2012 (SMAPVEX12) tenue du 6 juin au 17 juillet 2012 au Manitoba (Canada). Quatre sites forestiers de feuillus ont été échantillonnés. L’espèce majoritaire présente est le peuplier faux-tremble. Les données utilisées incluent des mesures de l’humidité du sol, de la rugosité de surface du sol, des caractéristiques des sites forestiers (arbres, sous-bois, litières…) et des données radar entièrement polarimétriques aéroportées et satellitaires acquises respectivement, en bande L (UAVSAR) à 30˚ et 40˚ et en bande C (RADARSAT-2) entre 20˚ et 30˚. Plusieurs paramètres polarimétriques ont été dérivés des données UAVSAR et RADARSAT-2 : les coefficients de corrélation (ρHHVV, φHHVV, etc); la hauteur du socle; l’entropie (H), l’anisotropie (A) et l’angle alpha extraits de la décomposition de Cloude-Pottier; les puissances de diffusion de surface (Ps), de double bond (Pd) extraites de la décomposition de Freeman-Durden, etc. Des relations entre les données radar (coefficients de rétrodiffusion multifréquences et multipolarisations (linéaires et circulaires) et les paramètres polarimétriques) et l’humidité du sol ont été développées et analysées. Les résultats ont montré que 1) En bande L, plusieurs paramètres optimaux permettent le suivi de l’humidité du sol en zone forestière avec un coefficient de corrélation significatif (p-value < 0,05): σ[indice supérieur 0] linéaire et σ[indice supérieur 0] circulaire (le coefficient de corrélation, r, varie entre 0,60 et 0,96), Ps (r entre 0,59 et 0,84), Pd (r entre 0,6 et 0,82), ρHHHV_30˚, ρVVHV_30˚, φHHHV_30˚ and φHHVV_30˚ (r entre 0,56 et 0,81) alors qu’en bande C, ils sont réduits à φHHHV, φVVHV et φHHVV (r est autour de 0,90). 2) En bande L, les paramètres polarimétriques n’ont pas montré de valeur ajoutée par rapport aux signaux conventionnels multipolarisés d’amplitude, pour le suivi de l’humidité du sol sur les sites forestiers. En revanche, en bande C, certains paramètres polarimétriques ont montré de meilleures relations significatives avec l’humidité du sol que les signaux conventionnels multipolarisés d’amplitude.

Tracking the attenuation and nonbreaking dissipation of swells using altimeters

A method for systematically tracking swells across oceanic basins is developed by taking advantage of high-quality data from space-borne altimeters and wave model output. The evolution of swells is observed over large distances based on 202 swell events with periods ranging from 12 to 18 s. An empirical attenuation rate of swell energy of about 4 × 10−7 m−1 is estimated using these observations, and the nonbreaking energy dissipation rates of swells far away from their generating areas are also estimated using a point source model. The resulting acceptance range of nonbreaking dissipation rates is −2.5 to 5.0 × 10−7 m−1, which corresponds to a dissipation e-folding scales of at least 2000 km for steep swells, to almost infinite for small-amplitude swells. These resulting rates are consistent with previous studies using in-situ and synthetic aperture radar (SAR) observations. The frequency dispersion and angular spreading effects during swell propagation are discussed by comparing the results with other studies, demonstrating that they are the two dominant processes for swell height attenuation, especially in the near field. The resulting dissipation rates from these observations can be used as a reference for ocean engineering and wave modeling, and for related studies such as air-sea and wind-wave-turbulence interactions.

Swell dissipation from 10 years of Envisat ASAR in wave mode

Swells are found in all oceans and strongly influence the wave climate and air-sea processes. The poorly known swell dissipation is the largest source of error in wave forecasts and hindcasts. We use synthetic aperture radar data to identify swell sources and trajectories, allowing a statistically significant estimation of swell dissipation. We mined the entire Envisat mission 2003–2012 to find suitable storms with swells (13 < T < 18 s) that are observed several times along their propagation. This database of swell events provides a comprehensive view of swell extending previous efforts. The analysis reveals that swell dissipation weakly correlates with the wave steepness, wind speed, orbital wave velocity, and the relative direction of wind and waves. Although several negative dissipation rates are found, there are uncertainties in the synthetic aperture radar-derived swell heights and dissipation rates. An acceptable range of the swell dissipation rate is −0.1 to 6 × 10−7 m−1 with a median of 1 × 10−7 m−1.

Extrapolating Satellite Winds to Turbine Operating Heights

Ocean wind retrievals from satellite sensors are typically performed for the standard level of 10 m. This restricts their full exploitation for wind energy planning, which requires wind information at much higher levels where wind turbines operate. A new method is presented for the vertical extrapolation of satellite-based wind maps. Winds near the sea surface are obtained from satellite data and used together with an adaptation of the Monin–Obukhov similarity theory to estimate the wind speed at higher levels. The thermal stratification of the atmosphere is taken into account through a long-term stability correction that is based on numerical weather prediction (NWP) model outputs. The effect of the long-term stability correction on the wind profile is significant. The method is applied to Envisat Advanced Synthetic Aperture Radar scenes acquired over the south Baltic Sea. This leads to maps of the long-term stability correction and wind speed at a height of 100 m with a spatial resolution of 0.02°. Calculations of the corresponding wind power density and Weibull parameters are shown. Comparisons with mast observations reveal that NWP model outputs can correct successfully for long-term stability effects and also, to some extent, for the limited number of satellite samples. The satellite-based and NWP-simulated wind profiles are almost equally accurate with respect to those from the mast. However, the satellite-based maps have a higher spatial resolution, which is particularly important in nearshore areas where most offshore wind farms are built.

Landscape Patterns of Permafrost Disturbance and Degradation in the Canadian High Arctic

Intensification of permafrost disturbances such as active layer detachments (ALDs) and retrogressive thaw slumps (RTS) have been observed across the circumpolar Arctic. These features are indicators of unstable conditions stemming from recent climate warming and permafrost degradation. In order to understand the processes interacting to give rise to these features, a multidisciplinary approach is required; i.e., interactions between geomorphology, hydrology, vegetation and ground thermal conditions. The goal of this research is to detect and map permafrost disturbance, predict landscape controls over disturbance and determine approaches for monitoring disturbance, all with the goal of contributing to the mitigation of permafrost hazards. Permafrost disturbance inventories were created by applying semi-automatic change detection techniques to IKONOS satellite imagery collected at the Cape Bounty Arctic Watershed Observatory (CBAWO). These methods provide a means to estimate the spatial distribution of permafrost disturbances for a given area for use as an input in susceptibility modelling. Permafrost disturbance susceptibility models were then developed using generalized additive and generalized linear models (GAM, GLM) fitted to disturbed and undisturbed locations and relevant GIS-derived predictor variables (slope, potential solar radiation, elevation). These models successfully delineated areas across the landscape that were susceptible to disturbances locally and regionally when transferred to an independent validation location. Permafrost disturbance susceptibility models are a first-order assessment of landscape susceptibility and are promising for designing land management strategies for remote permafrost regions. Additionally, geomorphic patterns associated with higher susceptibility provide important knowledge about processes associated with the initiation of disturbances. Permafrost degradation was analyzed at the CBAWO using differential interferometric synthetic aperture radar (DInSAR). Active-layer dynamics were interpreted using inter-seasonal and intra-seasonal displacement measurements and highlight the importance of hydroclimatic factors on active layer change. Collectively, these research approaches contribute to permafrost monitoring and the assessment of landscape-scale vulnerability in order to develop permafrost disturbance mitigation strategies.

Recent advances in radar technology

The technical developments and advances that have taken place thus far are reviewed in those areas impacting future phased array active aperture radar systems. The areas covered are printed circuit antennas and antenna arrays, GaAs MMIC design and fabrication leading to affordable transmitter-receiver (T-R) modules, and novel hardware and software developments. The use of fiber-optic distribution networks to interconnect the monolithically integrated optical components with the T-R modules is discussed. Beamforming and sidelobe control techniques for active phased array systems are also examined.

合成孔径激光成像雷达(Ⅰ):离焦和相位偏置望远镜接收天线

把合成孔径激光成像雷达的目标衍射区分为三个区域,提出采用离焦或者附加空间相位调制板的光学接收望远镜补偿回波像差。当目标处于菲涅耳衍射区时可采用离焦或偏置望远镜补偿回波二次项离焦像差并产生用于孔径合成的二次项相位历程;目标处于夫琅和费衍射区时可以采用离焦或偏置望远镜补偿回波二次项离焦像差但不产生相位历程;目标处于瑞利-索末菲衍射区域时不可能补偿回波高阶像差。

合成孔径激光成像雷达(Ⅱ):空间相位偏置发射望远镜

报道一种可以进行空间相位偏置的光学望远镜,用作合成孔径激光成像雷达中的光学发射天线。在望远镜内放置相位调制平板,控制望远镜的离焦量和位相调制平板的相位函数,能够在激光望远镜的照明区产生可控制的附加空间相位二次项,灵活改变激光照明波前,以在目标回波接收信号中产生雷达运动方向上的所需的二次项相位历程,因此能够实现特定的方位向成像分辨率。

合成孔径激光成像雷达(Ⅲ):双向环路发射接收望远镜

提出了一种用于合成孔径激光成像雷达的双向环路结构的发射接收望远镜,双向环路包括发射4-f转像系统、接收4-f转像系统和独立的望远镜。发射通道中设置离焦和相位调制平板偏置,接收通道中设置离焦和相位平板偏置。控制发射离焦量,发射相位调制函数,接收离焦量,接收相位调制函数,用同一个望远镜可以同时实现空间二次项相位附加偏置的激光发射和消除目标点散射回波接收波面像差的离焦光学接收,并产生雷达运动方向上合适的和可控制的相位二次项历程,从而实现孔径合成成像。详细介绍了系统设计,给出了从发射到光电外差接收的全过程传输方程。

Study of the effects of moving targets in SAR images for their detection.

Synthetic Aperture Radar’s (SAR) are systems designed in the early 50’s that are capable of obtaining images of the ground using electromagnetic signals. Thus, its activity is not interrupted by adverse meteorological conditions or during the night, as it occurs in optical systems. The name of the system comes from the creation of a synthetic aperture, larger than the real one, by moving the platform that carries the radar (typically a plane or a satellite). It provides the same resolution as a static radar equipped with a larger antenna. As it moves, the radar keeps emitting pulses every 1/PRF seconds —the PRF is the pulse repetition frequency—, whose echoes are stored and processed to obtain the image of the ground. To carry out this process, the algorithm needs to make the assumption that the targets in the illuminated scene are not moving. If that is the case, the algorithm is able to extract a focused image from the signal. However, if the targets are moving, they get unfocused and/or shifted from their position in the final image. There are applications in which it is especially useful to have information about moving targets (military, rescue tasks,studyoftheflowsofwater,surveillanceofmaritimeroutes...).Thisfeatureiscalled Ground Moving Target Indicator (GMTI). That is why the study and the development of techniques capable of detecting these targets and placing them correctly in the scene is convenient. In this document, some of the principal GMTI algorithms used in SAR systems are detailed. A simulator has been created to test the features of each implemented algorithm on a general situation with moving targets. Finally Monte Carlo tests have been performed, allowing us to extract conclusions and statistics of each algorithm.

Effective electromagnetic noise cancellation with beamformers and synthetic gradiometry in shielded and partly shielded environments

The major challenge of MEG, the inverse problem, is to estimate the very weak primary neuronal currents from the measurements of extracranial magnetic fields. The non-uniqueness of this inverse solution is compounded by the fact that MEG signals contain large environmental and physiological noise that further complicates the problem. In this paper, we evaluate the effectiveness of magnetic noise cancellation by synthetic gradiometers and the beamformer analysis method of synthetic aperture magnetometry (SAM) for source localisation in the presence of large stimulus-generated noise. We demonstrate that activation of primary somatosensory cortex can be accurately identified using SAM despite the presence of significant stimulus-related magnetic interference. This interference was generated by a contact heat evoked potential stimulator (CHEPS), recently developed for thermal pain research, but which to date has not been used in a MEG environment. We also show that in a reduced shielding environment the use of higher order synthetic gradiometry is sufficient to obtain signal-to-noise ratios (SNRs) that allow for accurate localisation of cortical sensory function.