TanDEM-X High Resolution DEMs and Their Applications to Flow Modeling

Lava flow modeling can be a powerful tool in hazard assessments; however, the ability to produce accurate models is usually limited by a lack of high resolution, up-to-date Digital Elevation Models (DEMs). This is especially obvious in places such as Kilauea Volcano (Hawaii), where active lava flows frequently alter the terrain. In this study, we use a new technique to create high resolution DEMs on Kilauea using synthetic aperture radar (SAR) data from the TanDEM-X (TDX) satellite. We convert raw TDX SAR data into a geocoded DEM using GAMMA software [Werner et al., 2000]. This process can be completed in several hours and permits creation of updated DEMs as soon as new TDX data are available. To test the DEMs, we use the Harris and Rowland [2001] FLOWGO lava flow model combined with the Favalli et al. [2005] DOWNFLOW model to simulate the 3-15 August 2011 eruption on Kilauea's East Rift Zone. Results were compared with simulations using the older, lower resolution 2000 SRTM DEM of Hawaii. Effusion rates used in the model are derived from MODIS thermal infrared satellite imagery. FLOWGO simulations using the TDX DEM produced a single flow line that matched the August 2011 flow almost perfectly, but could not recreate the entire flow field due to the relatively high DEM noise level. The issues with short model flow lengths can be resolved by filtering noise from the DEM. Model simulations using the outdated SRTM DEM produced a flow field that followed a different trajectory to that observed. Numerous lava flows have been emplaced at Kilauea since the creation of the SRTM DEM, leading the model to project flow lines in areas that have since been covered by fresh lava flows. These results show that DEMs can quickly become outdated on active volcanoes, but our new technique offers the potential to produce accurate, updated DEMs for modeling lava flow hazards.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

(Tables 1,2) Tracking records and space use statistics of polar bears (Ursus maritimus) in Storfjorden, Svalbard

Arctic sea ice is declining rapidly, making it vital to understand the importance of different types of sea ice for ice-dependent species such as polar bears Ursus maritimus. In this study we used GPS telemetry (25 polar bear tracks obtained in Svalbard, Norway, during spring) and high-resolution synthetic aperture radar (SAR) sea-ice data to investigate fine-scale space use by female polar bears. Space use patterns differed according to reproductive state; females with cubs of the year (COYs) had smaller home ranges and used fast-ice areas more frequently than lone females. First-passage time (FPT) analysis revealed that females with COYs displayed significantly longer FPTs near (<10 km) glacier fronts than in other fast-ice areas; lone females also increased their FPTs in such areas, but they also frequently used drifting pack ice. These results clearly demonstrate the importance of fast-ice areas, in particular close to glacier fronts, especially for females with COYs. Access to abundant and predictable prey (ringed seal pups), energy conservation and reluctance to cross large open water areas are possible reasons for the observed patterns. However, glacier fronts are retracting in Svalbard, and declines in land-fast ice have been notable over the past decade. The eventual disappearance of these important habitats might become critical for the survival of polar bear cubs in Svalbard and other regions with similar habitat characteristics. Given the relatively small size of many fast-ice areas in Svalbard, the results observed in this study would not have been revealed using less accurate location data or lower-resolution sea-ice data.

Estimación de la Reflectividad del Multitrayecto en Entornos Urbanos en Radares de Apertura Sintética

Synthetic Aperture Radar (SAR) images a target region reflectivity function in the multi-dimensional spatial domain of range and cross-range. SAR synthesizes a large aperture radar in order to achieve finer azimuth resolution than the one provided by any on-board real antenna. Conventional SAR techniques assume a single reflection of transmitted waveforms from targets. Nevertheless, today¿s new scenes force SAR systems to work in urban environments. Consequently, multiple-bounce returns are added to direct-scatter echoes. We refer to these as ghost images, since they obscure true target image and lead to poor resolution. By analyzing the quadratic phase error (QPE), this paper demonstrates that Earth¿s curvature influences the defocusing degree of multipath returns. In addition to the QPE, other parameters such as integrated sidelobe ratio (ISLR), peak sidelobe ratio (PSLR), contrast and entropy provide us with the tools to identify direct-scatter echoes in images containing undesired returns coming from multipath.

New Technique for Static Target Detection in Dense-Multipath Urban Environments

Synthetic Aperture Radar (SAR) images a target region reflectivity function in the multi-dimensional spatial domain of range and cross-range with a finer azimuth resolution than the one provided by any on-board real antenna. Conventional SAR techniques assume a single reflection of transmitted waveforms from targets. Nevertheless, new uses of Unmanned Aerial Vehicles (UAVs) for civilian-security applications force SAR systems to work in much more complex scenes such as urban environments. Consequently, multiple-bounce returns are additionally superposed to direct-scatter echoes. They are known as ghost images, since they obscure true target image and lead to poor resolution. All this may involve a significant problem in applications related to surveillance and security. In this work, an innovative multipath mitigation technique is presented in which Time Reversal (TR) concept is applied to SAR images when the target is concealed in clutter, leading to TR-SAR technique. This way, the effect of multipath is considerably reduced ?or even removed?, recovering the lost resolution due to multipath propagation. Furthermore, some focusing indicators such as entropy (E), contrast (C) and Rényi entropy (RE) provide us with a good focusing criterion when using TR-SAR.

Polarimetric Response of Rice Fields at C-Band: Analysis and Phenology Retrieval

A set of ten RADARSAT-2 images acquired in fully polarimetric mode over a test site with rice fields in Seville, Spain, has been analyzed to extract the main features of the C-band radar backscatter as a function of rice phenology. After observing the evolutions versus phenology of different polarimetric observables and explaining their behavior in terms of scattering mechanisms present in the scene, a simple retrieval approach has been proposed. This algorithm is based on three polarimetric observables and provides estimates from a set of four relevant intervals of phenological stages. The validation against ground data, carried out at parcel level for a set of six stands and up to nine dates per stand, provides a 96% rate of coincidence. Moreover, an equivalent compact-pol retrieval algorithm has been also proposed and validated, providing the same performance at parcel level. In all cases, the inversion is carried out by exploiting a single satellite acquisition, without any other auxiliary information.

Crop Phenology Estimation Using a Multitemporal Model and a Kalman Filtering Strategy

In this letter, a new approach for crop phenology estimation with remote sensing is presented. The proposed methodology is aimed to exploit tools from a dynamical system context. From a temporal sequence of images, a geometrical model is derived, which allows us to translate this temporal domain into the estimation problem. The evolution model in state space is obtained through dimensional reduction by a principal component analysis, defining the state variables, of the observations. Then, estimation is achieved by combining the generated model with actual samples in an optimal way using a Kalman filter. As a proof of concept, an example with results obtained with this approach over rice fields by exploiting stacks of TerraSAR-X dual polarization images is shown.

Facilitating soil-and plant-water research through the use of TDR

The development of TDR for measurement of soil water content and electrical conductivity has resulted in a large shift in measurement methods for a breadth of soil and hydrological characterization efforts. TDR has also opened new possibilities for soil and plant research. Five examples show how TDR has enhanced our ability to conduct our soil- and plant-water research. (i) Oxygen is necessary for healthy root growth and plant development but quantitative evaluation of the factors controlling oxygen supply in soil depends on knowledge of the soil water content by TDR. With water content information we have modeled successfully some impact of tillage methods on oxygen supply to roots and their growth response. (ii) For field assessment of soil mechanical properties influencing crop growth, water content capability was added to two portable soil strength measuring devices; (a) A TDT (Time Domain Transmittivity)-equipped soil cone penetrometer was used to evaluate seasonal soil strengthwater content relationships. In conventional tillage systems the relationships are dynamic and achieve the more stable no-tillage relationships only relatively late in each growing season; (b) A small TDR transmission line was added to a modified sheargraph that allowed shear strength and water content to be measured simultaneously on the same sample. In addition, the conventional graphing procedure for data acquisition was converted to datalogging using strain gauges. Data acquisition rate was improved by more than a factor of three with improved data quality. (iii) How do drought tolerant plants maintain leaf water content? Non-destructive measurement of TDR water content using a flat serpentine triple wire transmission line replaces more lengthy procedures of measuring relative water content. Two challenges remain: drought-stressed leaves alter salt content, changing electrical conductivity, and drought induced changes in leaf morphology affect TDR measurements. (iv) Remote radar signals are reflected from within the first 2 cm of soil. Appropriate calibration of radar imaging for soil water content can be achieved by a parallel pair of blades separated by 8 cm, reaching 1.7 cm into soil and forming a 20 cm TDR transmission line. The correlation between apparent relative permittivity from TDR and synthetic aperture radar (SAR) backscatter coefficient was 0.57 from an airborne flyover. These five examples highlight the diversity in the application of TDR in soil and plant research.

Intelligent Computations for Flood Monitoring

Floods represent the most devastating natural hazards in the world, affecting more people and causing more property damage than any other natural phenomena. One of the important problems associated with flood monitoring is flood extent extraction from satellite imagery, since it is impractical to acquire the flood area through field observations. This paper presents a method to flood extent extraction from synthetic-aperture radar (SAR) images that is based on intelligent computations. In particular, we apply artificial neural networks, self-organizing Kohonen’s maps (SOMs), for SAR image segmentation and classification. We tested our approach to process data from three different satellite sensors: ERS-2/SAR (during flooding on Tisza river, Ukraine and Hungary, 2001), ENVISAT/ASAR WSM (Wide Swath Mode) and RADARSAT-1 (during flooding on Huaihe river, China, 2007). Obtained results showed the efficiency of our approach.

Evaluation of TerraSAR-X Observations for Wetland InSAR Application

This paper assesses the potential of using spaceborne X-band synthetic aperture radar (SAR) data for monitoring water-level changes over wetlands. Our analysis is based on three sets of TerraSAR-X (TSX) observations acquired over South Florida's Everglades wetlands during an eight-month period in 2008. The first set was acquired in single HH polarization stripmap mode over our northern study area, consisting of managed wetlands and urban environments. The second set was acquired in dual-polarization stripmap mode over the western half of the same area, consisting mostly of managed wetlands. The third set was also acquired with dual-polarization stripmap mode over our southern study area, consisting of natural flow freshand salt-water wetlands in the southern Everglades. The first data set was used for a proof-of-concept study to verify that X-band data can generate coherent interferograms in wetland areas. Interferometric processing of this data set shows a high level of coherence (> 0.35) over both wetland and urban regions, maintaining interferometric phase in all three interferograms spanning 11 days. Surprisingly, phase is maintained over some of the wetlands even for interferograms spanning 33 days. The other two data sets were used to evaluate interferometric coherence of all four polarization modes and to determine dominant scattering mechanism in each wetland environment. Our results show high coherence values (> 0.4) in all polarization modes, with highest values in HH, then VV, and lowest in HV or VH. Interferograms calculated from multipolarization data show very similar fringe patterns regardless of the polarization type, suggesting that the phase information in all polarization data reflects water-level changes in wetlands and that volume scattering may be less important than commonly believed. We also used the two multipolarization data sets to conduct the Pauli decomposition, finding a strong dependence of scattering mechanism on vegetation t- - ype. The high interferometric coherence level of all polarization data suggests that a significant part of the X-band scattered signal interacts with lower sections of the vegetation (trunks and branches), because scattering from wind-affected canopies cannot support such a high coherence level. The high spatial resolution of TSX, combined with its 11-day repeat orbit, makes this X-band sensor surprisingly suitable for wetland interferometric SAR applications.

Helicopter-borne sea ice thickness measurements in the Storfjord, Spitsbergen from flight HEM_STJ06_01

Results from electromagnetic induction surveys of sea-ice thickness in Storfjorden, Svalbard, reveal large interannual ice-thickness variations in a region which is typically characterized by a reoccurring polynya. The surveys were performed in March 2003, May 2006 and March 2007 with helicopter- and ship-based sensors. The thickness distributions are influenced by sea-ice and atmospheric boundary conditions 2 months prior to the surveys, which are assessed with synthetic aperture radar (SAR) images, regional QuikSCAT backscatter maps and wind information from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis dataset. Locally formed thin ice from the Storfjorden polynya was frequently observed in 2003 and 2007 (mean thickness 0.55 and 0.37 m, respectively) because these years were characterized by prevailing northeasterly winds. In contrast, the entire fjord was covered with thick external sea ice in 2006 (mean thickness 2.21 m), when ice from the Barents Sea was driven into the fjord by predominantly southerly winds. The modal thickness of this external ice in 2006 increased from 1.2 m in the northern fjord to 2.4 m in the southern fjord, indicating stronger deformation in the southern part. This dynamically thickened ice was even thicker than multi-year ice advected from the central Arctic Ocean in 2003 (mean thickness 1.83 m). The thermodynamic ice thickness of fast ice as boundary condition is investigated with a one-dimensional sea-ice growth model (1DICE) forced with meteorological data from the weather station at the island of Hopen, southeast of Storfjorden. The model results are in good agreement with the modal thicknesses of fast-ice measurements in all years.

Circumarctic wetland dataset based on ENVISAT ASAR Global Monitoring Mode with links to geotiff image and NetCDF files (2005-01 to 2011-12)

A circumpolar representative and consistent wetland map is required for a range of applications ranging from upscaling of carbon fluxes and pools to climate modelling and wildlife habitat assessment. Currently available data sets lack sufficient accuracy and/or thematic detail in many regions of the Arctic. Synthetic aperture radar (SAR) data from satellites have already been shown to be suitable for wetland mapping. Envisat Advanced SAR (ASAR) provides global medium-resolution data which are examined with particular focus on spatial wetness patterns in this study. It was found that winter minimum backscatter values as well as their differences to summer minimum values reflect vegetation physiognomy units of certain wetness regimes. Low winter backscatter values are mostly found in areas vegetated by plant communities typically for wet regions in the tundra biome, due to low roughness and low volume scattering caused by the predominant vegetation. Summer to winter difference backscatter values, which in contrast to the winter values depend almost solely on soil moisture content, show expected higher values for wet regions. While the approach using difference values would seem more reasonable in order to delineate wetness patterns considering its direct link to soil moisture, it was found that a classification of winter minimum backscatter values is more applicable in tundra regions due to its better separability into wetness classes. Previous approaches for wetland detection have investigated the impact of liquid water in the soil on backscatter conditions. In this study the absence of liquid water is utilized. Owing to a lack of comparable regional to circumpolar data with respect to thematic detail, a potential wetland map cannot directly be validated; however, one might claim the validity of such a product by comparison with vegetation maps, which hold some information on the wetness status of certain classes. It was shown that the Envisat ASAR-derived classes are related to wetland classes of conventional vegetation maps, indicating its applicability; 30% of the land area north of the treeline was identified as wetland while conventional maps recorded 1-7%.

Circumpolar dataset of Soil Organic Carbon north of treeline derived from ENVISAT ASAR GM, link to GeoTIFF

A new approach for the estimation of soil organic carbon (SOC) pools north of the tree line has been developed based on synthetic aperture radar (SAR; ENVISAT Advanced SAR Global Monitoring mode) data. SOC values are directly determined from backscatter values instead of upscaling using land cover or soil classes. The multi-mode capability of SAR allows application across scales. It can be shown that measurements in C band under frozen conditions represent vegetation and surface structure properties which relate to soil properties, specifically SOC. It is estimated that at least 29 Pg C is stored in the upper 30 cm of soils north of the tree line. This is approximately 25 % less than stocks derived from the soil-map-based Northern Circumpolar Soil Carbon Database (NCSCD). The total stored carbon is underestimated since the established empirical relationship is not valid for peatlands or strongly cryoturbated soils. The approach does, however, provide the first spatially consistent account of soil organic carbon across the Arctic. Furthermore, it could be shown that values obtained from 1 km resolution SAR correspond to accounts based on a high spatial resolution (2 m) land cover map over a study area of about 7 × 7 km in NE Siberia. The approach can be also potentially transferred to medium-resolution C-band SAR data such as ENVISAT ASAR Wide Swath with ~120 m resolution but it is in general limited to regions without woody vegetation. Global Monitoring-mode-derived SOC increases with unfrozen period length. This indicates the importance of this parameter for modelling of the spatial distribution of soil organic carbon storage.

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.