Simple Approaches to Oil Spill Detection Using Sentinel Application Platform (SNAP)-Ocean Application Tools and Texture Analysis: A Comparative Study

Effective and efficient monitoring of oil spills that originate from ships, offshore platforms and any accidents are of immense importance from the viewpoint of public safety and environmental protection. Detection of spilled oil is also essential to estimate the potential spread and drift from the source to the nearby coastal areas. In this regard, utilization of SAR data for the detection and monitoring of oil spills has received considerable attention in recent times, due to their wide area coverage, day-night and all-weather capabilities. In this paper, two oil spills incidents along the coast of Mumbai, India are investigated; (1) The 2010 oil spill that occurred after the MV MSC Chitra and MV Khalijia-3 collided and (2) the oil spill caused due to sinking of MV RAK carrier in 2011. Two simple and relatively quick approaches for oil spill detection have been applied to VV polarized Radarsat-2 imagery of the incidents and a comparison is made of the results obtained. The first approach utilizes the oil spill detection tool of Sentinel Application Platform (SNAP) and the second explores texture analysis using Grey Level co-occurrence matrix (GLCM). The results of the study show that texture analysis proves to be an efficient method for oil spill detection as compared to the SNAP oil spill detection tool. Nevertheless, both the proposed methodologies are useful for detection of oil spills and for subsequent utilization of the results, timely and cost effectively, for the calibration and validation of numerical models that predict oil spill dispersion trajectories.     

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.     

On the consistency of the current conventional EOP series and the celestial and terrestrial reference frames

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.     

Analysis of decade-long time series of GPS-based polar motion estimates at 15-min temporal resolution

We present results from the generation of 10-year-long continuous time series of the Earth’s polar motion at 15-min temporal resolution using Global Positioning System ground data. From our results, we infer an overall noise level in our high-rate polar motion time series of 60 \(\upmu \hbox {as}\) (RMS). However, a spectral decomposition of our estimates indicates a noise floor of 4 \(\upmu \hbox {as}\) at periods shorter than 2 days, which enables recovery of diurnal and semidiurnal tidally induced polar motion. We deliberately place no constraints on retrograde diurnal polar motion despite its inherent ambiguity with long-period nutation. With this approach, we are able to resolve damped manifestations of the effects of the diurnal ocean tides on retrograde polar motion. As such, our approach is at least capable of discriminating between a historical background nutation model that excludes the effects of the diurnal ocean tides and modern models that include those effects. To assess the quality of our polar motion solution outside of the retrograde diurnal frequency band, we focus on its capability to recover tidally driven and non-tidal variations manifesting at the ultra-rapid (intra-daily) and rapid (characterized by periods ranging from 2 to 20 days) periods. We find that our best estimates of diurnal and semidiurnal tidally induced polar motion result from an approach that adopts, at the observation level, a reasonable background model of these effects. We also demonstrate that our high-rate polar motion estimates yield similar results to daily-resolved polar motion estimates, and therefore do not compromise the ability to resolve polar motion at periods of 2–20 days.     

Adaptive filtering of GOCE-derived gravity gradients of the disturbing potential in the context of the space-wise approach

Filtering and signal processing techniques have been widely used in the processing of satellite gravity observations to reduce measurement noise and correlation errors. The parameters and types of filters used depend on the statistical and spectral properties of the signal under investigation. Filtering is usually applied in a non-real-time environment. The present work focuses on the implementation of an adaptive filtering technique to process satellite gravity gradiometry data for gravity field modeling. Adaptive filtering algorithms are commonly used in communication systems, noise and echo cancellation, and biomedical applications. Two independent studies have been performed to introduce adaptive signal processing techniques and test the performance of the least mean-squared (LMS) adaptive algorithm for filtering satellite measurements obtained by the gravity field and steady-state ocean circulation explorer (GOCE) mission. In the first study, a Monte Carlo simulation is performed in order to gain insights about the implementation of the LMS algorithm on data with spectral behavior close to that of real GOCE data. In the second study, the LMS algorithm is implemented on real GOCE data. Experiments are also performed to determine suitable filtering parameters. Only the four accurate components of the full GOCE gravity gradient tensor of the disturbing potential are used. The characteristics of the filtered gravity gradients are examined in the time and spectral domain. The obtained filtered GOCE gravity gradients show an agreement of 63–84 mEötvös (depending on the gravity gradient component), in terms of RMS error, when compared to the gravity gradients derived from the EGM2008 geopotential model. Spectral-domain analysis of the filtered gradients shows that the adaptive filters slightly suppress frequencies in the bandwidth of approximately 10–30 mHz. The limitations of the adaptive LMS algorithm are also discussed. The tested filtering algorithm can be connected to and employed in the first computational steps of the space-wise approach, where a time-wise Wiener filter is applied at the first stage of GOCE gravity gradient filtering. The results of this work can be extended to using other adaptive filtering algorithms, such as the recursive least-squares and recursive least-squares lattice filters.     

Enriching the semantics of the directed polyline–polygon topological relationships: the $$\mathcal {}$$-intersection matrix

In this paper, we define an intersection matrix for enriching the semantics of the topological relationships between a directed polyline and a polygon. In particular, we propose the \(\mathcal {DLP}\)-intersection matrix which enables us to model the origin and destination points, as well as the right- and left-hand sides of the directed polyline. This matrix overcomes the limitation of the well-known DE-9IM, because it allows the representation of the different dimensions of the intersection results at the same time. Accordingly, the geo-operators have been revised and extended in order to address the notions of right- and left-hand sides of a directed polyline, as well as additional notions related to the orientation of the polyline. The \(\mathcal {DLP}\)-intersection matrix has been implemented by extending the Java Topology Suite methods in order to address the new geo-operators based on the notion of orientation.     

A New Combined Assessment of Mixed Uncertainty in Spatial Models: Conceptualization and Implementation

Uncertainty quantification is not often performed in spatial modeling applications, especially when there is a mixture of probabilistic and non‐probabilistic uncertainties. Furthermore, the effect of positional uncertainty is often not assessed, despite its relevance to geographical applications. Although there has been much work in investigating the aforementioned types of uncertainty in isolation, combined approaches have not been much researched. This has resulted in a lack of tools for conducting mixed uncertainty analyses that include positional uncertainty. This research addresses the issue by first presenting a new, flexible, simulation‐oriented conceptualization of positional uncertainty in geographic objects called F‐Objects. F‐Objects accommodates various representations of uncertainty, while remaining conceptually simple. Second, a new Python‐based framework is introduced, termed Wiggly and capable of conducting mixed uncertainty propagation using fuzzy Monte Carlo simulation (FMCS). FMCS combines both traditional Monte Carlo with fuzzy analysis in a so‐called hybrid approach. F‐Objects is implemented within the Wiggly framework, resulting in a tool capable of considering any combination of: (1) probabilistic variables; (2) fuzzy variables; and (3) positional uncertainty of objects (probabilistic/fuzzy). Finally, a realistic GIS‐based groundwater contamination problem demonstrates how F‐Objects and Wiggly can be used to assess the effect of positional uncertainty.     

Application of Supervised Classification and CrostaTechnique for Lithological Discrimination in Parts of South Khetri Belt,Sikar District,Rajasthan

Supervised classification and Crosta technique is widely used for lithological and alteration zone mapping respectively. Landsat ETM+ digital imagery has been used to generate Crosta and Supervised classification image using digital image processing (DIP) technique. In the study, The Crosta technique has been used for the first time in lithological mapping. These techniques were applied to distinguish the litho-contacts between Alwar and Ajabgarh quartzite in the south-western part of North Delhi Fold Belt (NDFB) on the basis of argillaceous and arenaceous nature of the rocks. Litho contacts between Alwar and Ajabgarh groups are varying from quartzite to biotite schist in the Khetri basin whereas in the study area the litho-contacts are Alwar quartzite to Ajabgarh quartzite which are very difficult to classify on naked eye. Therefore, this area has been selected for this study. Several authors have studied the geological setting of Khetri sub basin and classified on the basis of argillaceous, arenaceous nature and primary sedimentary structural features like ripple marks and cross bedding. This is the first attempt in the Khetri sub basin for lithological classification based on remote sensing and digital image processing technique. This study revealed that Crosta image has significant spatial correlation with the lithology discriminated using Supervised classification technique. Lithological variations were clearly demarcated using these technique around Rahunathgarh and Golyana areas. The litho units between these areas were marked under the Group of Ajabgarh quartzite by the GSI people during 1994. However, the present study classified those area under the Group of Alwar quartzite. Similar type of studies can also be carried out where these type of lithological problems arises.     

Removing Shadows Using RGB Color Space in Pairs of Optical Satellite Images

One of the greatest challenges for optical satellite images applications is the presence of shadows. In stereo correspondence of images for example, shadows obstruct the correct extraction of objects that degrades the quality of stereo matching results. The aim of this research is to present a new, simple and efficient shadow detection and removal approach. The proposed approach first detects shadows by operating directly in red, green and blue color space using a new method including spectral and spatial properties of shadow. Secondly, shadows are removed by supplying more light to the shadow’s region using an energy minimization concept. The edges of shadows are removed or attenuated using some filters. The experimental results show that the proposed shadow detection and removal approach can generate accurate and efficient recovered pairs of satellite images. Furthermore, we demonstrate its reliability on the application of a Hopfield neural matching by comparing the correspondence of images before and after shadow removal.