Earth Observation Scientific Workflows in a Distributed Computing Environment

Geospatially Enabled Scientific Workflows offer a promising toolset to help researchers in the earth observation domain with many aspects of the scientific process. One such aspect is that of access to distributed earth observation data and computing resources. Earth observation research often utilizes large datasets requiring extensive CPU and memory resources in their processing. These resource intensive processes can be chained; the sequence of processes (and their provenance) makes up a scientific workflow. Despite the exponential growth in capacity of desktop computers, their resources are often insufficient for the scientific workflow processing tasks at hand. By integrating distributed computing capabilities into a geospatially enabled scientific workflow environment, it is possible to provide researchers with a mechanism to overcome the limitations of the desktop computer. Most of the effort on extending scientific workflows with distributed computing capabilities has focused on the web services approach, as exemplified by the OGC's Web Processing Service and by GRID computing. The approach to leveraging distributed computing resources described in this article uses instead remote objects via RPyC and the dynamic properties of the Python programming language. The Vistrails environment has been extended to allow for geospatial processing through the EO4Vistrails package ( http://code.google.com/p/eo4vistrails/ ). In order to allow these geospatial processes to be seamlessly executed on distributed resources such as cloud computing nodes, the Vistrails environment has been extended with both multi‐tasking capabilities and distributed processing capabilities. The multi‐tasking capabilities are required in order to allow Vistrails to run side‐by‐side processes, a capability it does not currently have. The distributed processing capabilities are achieved through the use of remote objects and mobile code through RPyC.     

Detecting anomalous CO2 flux using space borne spectroscopy

Over the time-scale, earth's atmospheric CO₂ concentration has varied and that is mostly determined by balance among the geochemical processes including burial of organic carbon in sediments, silicate rock weathering and volcanic activity. The best recorded atmospheric CO₂ variability is derived from Vostok ice core that records last four glacial/interglacial cycles. The present CO₂ concentration of earth's atmosphere has exceeded far that it was predicted from the ice core data. Other than rapid industrialization and urbanization since last century, geo-natural hazards such as volcanic activity, leakage from hydrocarbon reservoirs and spontaneous combustion of coal contribute a considerable amount of CO₂ to the atmosphere. Spontaneous combustion of coal is common occurrence in most coal producing countries and sometimes it could be in an enormous scale. Remote sensing has already proved to be a significant tool in coalfire identification and monitoring studies. However, coalfire related CO₂ quantification from remote sensing data has not endeavoured yet by scientific communities because of low spectral resolution of commercially available remote sensing data and relatively sparse CO₂ plume than other geological hazards like volcanic activity. The present research has attempted two methods to identify the CO₂ flux emitted from coalfires in a coalmining region in north China. Firstly, a band rationing method was used for column atmospheric retrieval of CO₂ and secondly atmospheric models were simulated in fast atmospheric signature code (FASCOD) to understand the local radiation transport and then the model was implemented with the inputs from hyperspectral remote sensing data. It was observed that retrieval of columnar abundance of CO₂ with the band rationing method is faster as less simulation required in FASCOD. Alternatively, the inversion model could retrieve CO₂ concentration from a (certain) source because it excludes the uncertainties in the higher altitude.     

Strategies to mitigate aliasing of loading signals while estimating GPS frame parameters

Although GNSS techniques are theoretically sensitive to the Earth center of mass, it is often preferable to remove intrinsic origin and scale information from the estimated station positions since they are known to be affected by systematic errors. This is usually done by estimating the parameters of a linearized similarity transformation which relates the quasi-instantaneous frames to a long-term frame such as the International Terrestrial Reference Frame (ITRF). It is well known that non-linear station motions can partially alias into these parameters. We discuss in this paper some procedures that may allow reducing these aliasing effects in the case of the GPS techniques. The options include the use of well-distributed sub-networks for the frame transformation estimation, the use of site loading corrections, a modification of the stochastic model by downweighting heights, or the joint estimation of the low degrees of the deformation field. We confirm that the standard approach consisting of estimating the transformation over the whole network is particularly harmful for the loading signals if the network is not well distributed. Downweighting the height component, using a uniform sub-network, or estimating the deformation field perform similarly in drastically reducing the amplitude of the aliasing effect. The application of these methods to reprocessed GPS terrestrial frames permits an assessment of the level of agreement between GPS and our loading model, which is found to be about 1.5 mm WRMS in height and 0.8 mm WRMS in the horizontal at the annual frequency. Aliased loading signals are not the main source of discrepancies between loading displacement models and GPS position time series.     

Hyperspectral remote sensing of evaporate minerals and associated sediments in Lake Magadi area, Kenya

Pleistocene to present evaporitic lacustrine sediments in Lake Magadi, East African Rift Valley, Kenya were studied and mapped using spectral remote sensing methods. This approach incorporated surface mineral mapping using space-borne hyperspectral Hyperion imagery together with laboratory analysis, including visible, near-infrared diffuse reflectance spectroscopy (VNIR) measurements and X-ray diffraction for selected rock and soil samples of the study area. The spectral signatures of Magadiite and Kenyaite, which have not been previously reported, were established and the spectral signatures of trona, chert series, volcanic tuff and the High Magadi bed were also analyzed.Image processing techniques, MNF (Minimum Noise Fraction) and MTMF (Mixture Tuned Matched Filtering) using a stratified approach (image analysis with and without the lake area), were used to enhance the mapping of evaporates. High Magadi beds, chert series and volcanic tuff were identified from the Hyperion image with an overall mapping accuracy of 84.3%. Even though, the spatial distribution of evaporites and sediments in Lake Magadi area change in response to climate variations, the mineralogy of this area has not been mapped recently. The results of this study shows the usefulness of the hypersspectral remote sensing to map the surface geology of this kind of environment and to locate promising sites for industrial open-pit trona mining in a qualitative and quantitative manner.     

UML‐Based Approach to Developing a CityGML Application Domain Extension

Recently a national 3D standard was established in the Netherlands as a CityGML Application Domain Extension (called IMGeo). In line with the Dutch practice of modeling geo‐information, the ADE is developed using a model driven approach. The classes are designed in UML and automatically mapped to GML schema. The current OGC CityGML specification does not provide rules or guidance on correctly modeling an ADE in UML. This article fills this gap by studying how CityGML can be extended for specific applications starting from the UML diagrams. Six alternatives for modeling ADEs in UML are introduced and compared. The optimal alternative is selected and applied to obtain the national 3D standard. The approach was extensively discussed with international experts, who were members of both SIG3D and other working groups. As a consequence the approach was adopted by the SIG3D, the Special Interest Group 3D which, among other things, work on the 3D standard CityGML in cooperation with OGC. Therefore the approach contains many issues that can be generalized and reused by future domain extensions of CityGML. To further support this, the article formulates a model‐driven framework to model CityGML ADEs. Open issues are described in the conclusions.     

Multi‐sensor Image Maps from SPOT,ERS and JERS

Abstract

Image mapping using data from visible and infrared sensors has, as a major drawback, the frequent cloud cover experienced in many countries. This is one of the main reasons why topographic maps at 1:100,000 scale and larger are often outdated. The results of a study which investigated the possibilities of fusing up‐to‐date spaceborne microwave data with existing images from optical sensors for topographic map updating at a scale of 1:100, 000 are presented in this paper. A key issue researched was the influence of geometric distortions and corrections of remote sensing data on the results of pixel based digital image fusion. After having terrain‐geocoded and radiometrically enhanced imagery from the Landsat, SPOT, ERS‐1 and JERS‐1 satellites, the data were fused applying a variety of colour transformation techniques as well as statistical or arithmetic methods. Initially, the image fusion was implemented using images covering a test site in the north of The Netherlands in order to calibrate specified combinations and techniques in a rather flat area. With the experience gained, the remote sensing data acquired over the research site were processed. The research test site is located in a typical Developing Country in the humid Tropics, on the mountainous south‐west coast of Sumatra in Indonesia. The results of the various applied techniques and image combinations were evaluated with reference to their capability to overcome the cloud cover problem. New combinations of techniques and images were developed as result of an optimisation process. The research produced two prototypes of annotated 1:100,000 scale image maps containing fused, cloud‐free optical/microwave imagery.     

GPS signal acquisition using the repeatability of successive code phase measurements

This paper introduces a new method for GPS signal acquisition, which is based on the repeatability of successive code phase measurements and the M-of-N search algorithm. The performance of the proposed method in terms of probability of signal detection is similar to that of traditional methods, except that the calculation of the probability of detection does not rely on the noise distribution or the Carrier-to-Noise ratio (C/N₀). The code phase repeatability-based method is presented along with equations for probability of detection and probability of false detection. If the distribution of the noise is known, it also provides an estimate of the C/N₀. The proposed method is illustrated for coherent and non-coherent acquisition and C/N₀ estimation.     

Quality assessment of GPS reprocessed terrestrial reference frame

The International GNSS Service (IGS) contributes to the construction of the International Terrestrial Reference Frame (ITRF) by submitting time series of station positions and Earth Rotation Parameters (ERP). For the first time, its submission to the ITRF2008 construction is based on a combination of entirely reprocessed GPS solutions delivered by 11 Analysis Centers (ACs). We analyze the IGS submission and four of the individual AC contributions in terms of the GNSS frame origin and scale, station position repeatability and time series seasonal variations. We show here that the GPS Terrestrial Reference Frame (TRF) origin is consistent with Satellite laser Ranging (SLR) at the centimeter level with a drift lower than 1 mm/year. Although the scale drift compared to Very Long baseline Interferometry (VLBI) and SLR mean scale is smaller than 0.4 mm/year, we think that it would be premature to use that information in the ITRF scale definition due to its strong dependence on the GPS satellite and ground antenna phase center variations. The new position time series also show a better repeatability compared to past IGS combined products and their annual variations are shown to be more consistent with loading models. The comparison of GPS station positions and velocities to those of VLBI via local ties in co-located sites demonstrates that the IGS reprocessed solution submitted to the ITRF2008 is more reliable and precise than any of the past submissions. However, we show that some of the remaining inconsistencies between GPS and VLBI positioning may be caused by uncalibrated GNSS radomes.     

Soil Organic Carbon mapping of partially vegetated agricultural fields with imaging spectroscopy

Soil Organic Carbon (SOC) is one of the key soil properties, but the large spatial variation makes continuous mapping a complex task. Imaging spectroscopy has proven to be an useful technique for mapping of soil properties, but the applicability decreases rapidly when fields are partially covered with vegetation. In this paper we show that with only a few percent fractional maize cover the accuracy of a Partial Least Square Regression (PLSR) based SOC prediction model drops dramatically. However, this problem can be solved with the use of spectral unmixing techniques. First, the fractional maize cover is determined with linear spectral unmixing, taking the illumination and observation angles into account. In a next step the influence of maize is filtered out from the spectral signal by a new procedure termed Residual Spectral Unmixing (RSU). The residual soil spectra resulting from this procedure are used for mapping of SOC using PLSR, which could be done with accuracies comparable to studies performed on bare soil surfaces (Root Mean Standard Error of Calibration = 1.34 g/kg and Root Mean Standard Error of Prediction = 1.65 g/kg). With the presented RSU approach it is possible to filter out the influence of maize from the mixed spectra, and the residual soil spectra contain enough information for mapping of the SOC distribution within agricultural fields. This can improve the applicability of airborne imaging spectroscopy for soil studies in temperate climates, since the use of the RSU approach can extend the flight-window which is often constrained by the presence of vegetation.     

Remotely sensed estimation of forest canopy density: A comparison of the performance of four methods

In recent years, a number of alternative methods have been proposed to predict forest canopy density from remotely sensed data. To date, however, it remains difficult to decide which method to use, since their relative performance has never been evaluated. In this study the performance of: (1) an artificial neural network, (2) a multiple linear regression, (3) the forest canopy density mapper and (4) a maximum likelihood classification method was compared for prediction of forest canopy density using a Landsat ETM+ image. Comparison of confusion matrices revealed that the regression model performed significantly worse than the three other methods. These results were based on a z-test for comparison of weighted kappa statistics, which is an appropriate statistic for analysis of ranked categories. About 89% of the variance of the observed canopy density was explained by the artificial neural networks, which outperformed the other three methods in this respect. Moreover, the artificial neural networks gave an unbiased prediction, while other methods systematically under or over predicted forest canopy density. The choice of biased method could have a high impact on canopy density inventories.