Graph theory for analyzing pair-wise data: application to geophysical model parameters estimated from interferometric synthetic aperture radar data at Okmok volcano,Alaska

Graph theory is useful for analyzing time-dependent model parameters estimated from interferometric synthetic aperture radar (InSAR) data in the temporal domain. Plotting acquisition dates (epochs) as vertices and pair-wise interferometric combinations as edges defines an incidence graph. The edge-vertex incidence matrix and the normalized edge Laplacian matrix are factors in the covariance matrix for the pair-wise data. Using empirical measures of residual scatter in the pair-wise observations, we estimate the relative variance at each epoch by inverting the covariance of the pair-wise data. We evaluate the rank deficiency of the corresponding least-squares problem via the edge-vertex incidence matrix. We implement our method in a MATLAB software package called GraphTreeTA available on GitHub (https://github.com/feigl/gipht). We apply temporal adjustment to the data set described in Lu et al. (Geophys Res Solid Earth 110, 2005) at Okmok volcano, Alaska, which erupted most recently in 1997 and 2008. The data set contains 44 differential volumetric changes and uncertainties estimated from interferograms between 1997 and 2004. Estimates show that approximately half of the magma volume lost during the 1997 eruption was recovered by the summer of 2003. Between June 2002 and September 2003, the estimated rate of volumetric increase is \((6.2 \, \pm \, 0.6) \times 10^6~\mathrm{m}^3/\mathrm{year} \). Our preferred model provides a reasonable fit that is compatible with viscoelastic relaxation in the five years following the 1997 eruption. Although we demonstrate the approach using volumetric rates of change, our formulation in terms of incidence graphs applies to any quantity derived from pair-wise differences, such as range change, range gradient, or atmospheric delay.     

A novel mathematical morphology based algorithm for shoreline extraction from satellite images

Abstract

Shoreline extraction is fundamental and inevitable for several studies. Ascertaining the precise spatial location of the shoreline is crucial. Recently, the need for using remote sensing data to accomplish the complex task of automatic extraction of features, such as shoreline, has considerably increased. Automated feature extraction can drastically minimize the time and cost of data acquisition and database updating. Effective and fast approaches are essential to monitor coastline retreat and update shoreline maps. Here, we present a flexible mathematical morphology-driven approach for shoreline extraction algorithm from satellite imageries. The salient features of this work are the preservation of actual size and shape of the shorelines, run-time structuring element definition, semi-automation, faster processing, and single band adaptability. The proposed approach is tested with various sensor-driven images with low to high resolutions. Accuracy of the developed methodology has been assessed with manually prepared ground truths of the study area and compared with an existing shoreline classification approach. The proposed approach is found successful in shoreline extraction from the wide variety of satellite images based on the results drawn from visual and quantitative assessments.     

Modeling the environmental susceptibility of landfill sites in California

??The proper management of solid waste (SW) is a global environmental challenge. A major issue is the proper disposal of SW while balancing a wide range of criteria and working with different spatial data. In this study, we used geographic information system as a tool to perform multi-criteria decision analysis with an analytical hierarchy process to develop an environmental impact susceptibility model (EISM) for landfills. The model was applied to the state of California, USA and results are presented herein. In particular, the EISM considers factors such as geology, pedology, geomorphology, water resources, and climate as represented by 13 associated environmental indicators. The results of the EISM indicate that more than 75% of California’s territory is situated in areas with very low, low, and medium environmental impact susceptibility categories. However, in the remaining 25% of the state’s land, 61 landfills are located in the high and very high categories. These results are alarming because during the period from 2000 to 2015, these 61 landfills received approximately 308 million tons of SW, which corresponds to more than 57% of all SW disposed in California. The model results can be used toward mitigating the environmental impacts of these facilities.     

The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS

A new method for modeling the ionospheric delay using global positioning system (GPS) data is proposed, called the ionospheric eclipse factor method (IEFM). It is based on establishing a concept referred to as the ionospheric eclipse factor (IEF) λ of the ionospheric pierce point (IPP) and the IEF’s influence factor (IFF) . The IEF can be used to make a relatively precise distinction between ionospheric daytime and nighttime, whereas the IFF is advantageous for describing the IEF’s variations with day, month, season and year, associated with seasonal variations of total electron content (TEC) of the ionosphere. By combining λ and with the local time t of IPP, the IEFM has the ability to precisely distinguish between ionospheric daytime and nighttime, as well as efficiently combine them during different seasons or months over a year at the IPP. The IEFM-based ionospheric delay estimates are validated by combining an absolute positioning mode with several ionospheric delay correction models or algorithms, using GPS data at an international Global Navigation Satellite System (GNSS) service (IGS) station (WTZR). Our results indicate that the IEFM may further improve ionospheric delay modeling using GPS data.     

Atmospheric Artifact Compensation in Ground-Based DInSAR Applications

In this letter, a coherence-based technique for atmospheric artifact removal in ground-based (GB) zero-baseline synthetic aperture radar (SAR) acquisitions is proposed. For this purpose, polarimetric measurements acquired using the GB-SAR sensor developed at the Universitat Politecnica de Catalunya are employed. The heterogeneous environment of Collserola Park in the outskirts of Barcelona, Spain, was selected as the test area. Data sets were acquired at X-band during one week in June 2005. The effects of the atmosphere variations between successive zero-baseline SAR polarimetric acquisitions are treated here in detail. The need to compensate for the resulting phase-difference errors when retrieving interferometric information is put forward. A compensation technique is then proposed and evaluated using the control points placed inside the observed scene.     

Seasonal Variation of MODIS Vegetation Indexes and Their Statistical Relationship With Climate Over the Subtropic Evergreen Forest in Zhejiang, China

Moderate Resolution Imaging Spectroradiometer (MODIS) 16-day 1-km vegetation index products, daily temperature, photosynthetically active radiation (PAR), and precipitation from 2001 to 2004 were utilized to analyze the temporal variations of the MODIS normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI), as well as their correlations with climate over the evergreen forested sites in Zhejiang-a humid subtropical region in the southeast of China. The results showed that both NDVI and EVI could discern the seasonal variation of the evergreen forests. Attributed to the sufficient precipitation in the study area, the growth of vegetation is mainly controlled by energy; as a result, NDVI, and especially EVI, is more correlated with temperature and PAR than precipitation. Compared with NDVI, EVI is more sensitive to climate condition and is a better indicator to study vegetation variations in the study region     

Reduced-Dimensional Processing for Ground Moving Target Detection in Distributed Space-Based Radar

With multisatellite radar systems, several additional features are achieved: multistatic observation, interferometry, ground moving target indication (GMTI). In this letter, a new reduced-dimensional method based on joint pixels sum-difference (Sigma-Delta) data for clutter rejection and GMTI is proposed. The reduced-dimensional joint pixels Sigma-Delta data are obtained by the orthogonal projection of the joint pixels data of different synthetic aperture radar (SAR) images generated by a multisatellite radar system. In the sense of statistic expectation, the joint pixels Sigma-Delta data contain the common and different information among SAR images. Then, the objective of clutter cancellation and GMTI can be achieved by adaptive processing. Simulation results demonstrate the effectiveness and robustness of the proposed method even with clutter fluctuation and image coregistration errors     

On the link between GPS pseudorange noise and day-boundary discontinuities in geodetic time transfer solutions

When neglecting calibration issues, the accuracy of GPS-based time and frequency transfer using a combined analysis of code and carrier phase measurements highly depends on the noise of the GPS codes. In particular, the pseudorange noise is responsible for day-boundary discontinuities which can reach more than 1 ns in the time transfer results obtained from geodetic analysis. These discontinuities are caused by the fact that the data are analyzed in daily data batches where the absolute clock offset is determined by the mean code value during the daily data batch. This pseudorange noise is not a white noise, in particular due to multipath and variations of instrumental delays. In this paper, the pseudorange noise behavior is characterized in order to improve the understanding of the origin of the large day-boundary discontinuities in the geodetic time transfer results. In a first step, the effect of short-term noise and multipath is estimated, and shown to be responsible for only a maximum of 150 ps (picoseconds) of the day-boundary jumps, with only one exception at NRC1 where the correction provides a jump reduction of 300 ps. In a second step, a combination of time transfer results obtained with pseudoranges only and geodetic time transfer results is used to characterize the long-term evolution of pseudorange errors. It demonstrates that the day-boundary jumps, especially those of large amplitude, can be explained by an instrumental effect imposing a common behavior on all the satellite pseudoranges. Using known influences as temperature variations at ALGO or cable damages at HOB2, it is shown that the approach developed in this study can be used to look for the origin of the day-boundary discontinuities in other stations.     

Mapping impervious surfaces with the integrated use of Landsat Thematic Mapper and radar data: A case study in an urban–rural landscape in the Brazilian Amazon

This research explored the integrated use of Landsat Thematic Mapper (TM) and radar (i.e., ALOS PALSAR L-band and RADARSAT-2 C-band) data for mapping impervious surface distribution to examine the roles of radar data with different spatial resolutions and wavelengths. The wavelet-merging technique was used to merge TM and radar data to generate a new dataset. A constrained least-squares solution was used to unmix TM multispectral data and multisensor fusion images to four fraction images (high-albedo, low-albedo, vegetation, and soil). The impervious surface image was then extracted from the high-albedo and low-albedo fraction images. QuickBird imagery was used to develop an impervious surface image for use as reference data to evaluate the results from TM and fusion images. This research indicated that increasing spatial resolution by multisensor fusion improved spatial patterns of impervious surface distribution, but cannot significantly improve the statistical area accuracy. This research also indicated that the fusion image with 10-m spatial resolution was suitable for mapping impervious surface spatial distribution, but TM multispectral image with 30 m was too coarse in a complex urban–rural landscape. On the other hand, this research showed that no significant difference in improving impervious surface mapping performance by using either PALSAR L-band or RADARSAT C-band data with the same spatial resolution when they were used for multi-sensor fusion with the wavelet-based method.