Separation of global time-variable gravity signals into maximally independent components

The Gravity Recovery and Climate Experiment (GRACE) products provide valuable information about total water storage variations over the whole globe. Since GRACE detects mass variations integrated over vertical columns, it is desirable to separate its total water storage anomalies into their original sources. Among the statistical approaches, the principal component analysis (PCA) method and its extensions have been frequently proposed to decompose the GRACE products into space and time components. However, these methods only search for decorrelated components that on the one hand are not always interpretable and on the other hand often contain a superposition of independent source signals. In contrast, independent component analysis (ICA) represents a technique that separates components based on assumed statistical independence using higher-order statistical information. If one assumes that independent physical processes generate statistically independent signal components added up in the GRACE observations, separating them by ICA is a reliable strategy to identify these processes. In this paper, the performance of the conventional PCA, its rotated extension and ICA are investigated when applied to the GRACE-derived total water storage variations. These analyses have been tested on both a synthetic example and on the real GRACE level-2 monthly solutions derived from GeoForschungsZentrum Potsdam (GFZ RL04) and Bonn University (ITG2010). Within the synthetic example, we can show how imposing statistical independence in the framework of ICA improves the extraction of the ‘original’ signals from a GRACE-type super-position. We are therefore confident that also for the real case the ICA algorithm, without making prior assumptions about the long-term behaviour or on the frequencies contained in the signal, improves over the performance of PCA and its rotated extension in the separation of periodical and long-term components.     

Evaluating the impact of data quantity, distribution and algorithm selection on the accuracy of 3D subsurface models using synthetic grid models of varying complexity

Testing the accuracy of 3D modelling algorithms used for geological applications is extremely difficult as model results cannot be easily validated. This paper presents a new approach to evaluate the effectiveness of common interpolation algorithms used in 3D subsurface modelling, utilizing four synthetic grids to represent subsurface environments of varying geological complexity. The four grids are modelled with Inverse Distance Weighting and Ordinary Kriging, using data extracted from the synthetic grids in different spatial distribution patterns (regular, random, clustered and sparse), and with different numbers of data points (100, 256, 676 and 1,600). Utilizing synthetic grids for this evaluation allows quantitative statistical assessment of the accuracy of both interpolation algorithms in a variety of sampling conditions. Data distribution proved to be an important factor; as in many geological situations, relatively small numbers of randomly distributed data points can generate more accurate 3D models than larger amounts of clustered data. This study provides insight for optimizing the quantity and distribution of data required to accurately and cost-effectively interpolate subsurface units of varying complexity.     

Sub-daily alias and draconitic errors in the IGS orbits

Harmonic signals with a fundamental period near the GPS draconitic year (351.2 days) and overtones up to at least the sixth multiple have been observed in the power spectra of nearly all products of the International GNSS Service (IGS), including station position time series, apparent geocenter motions, orbit jumps between successive days, and midnight discontinuities in earth orientation parameter (EOP) rates. Two main mechanisms have been suggested for the harmonics: mismodeling of orbit dynamics and aliasing of near-sidereal local station multipath effects. Others have studied the propagation of local multipath errors into draconitic position variations, but orbit-related processes have been less examined. We elaborate our earlier analysis of GPS day-boundary orbit discontinuities where we observed some draconitic features as well as prominent spectral bands near 29-, 14-, 9-, and 7-day periods. Finer structures within the sub-seasonal bands fall close to the expected alias frequencies for 24-h sampling of sub-daily EOP tide lines but do not coincide precisely. While once-per-revolution empirical orbit parameters should strongly absorb any sub-daily EOP tide errors due to near-resonance of their respective periods, the observed differences require explanation. This has been done by simulating EOP tidal errors and checking their impact on a long series of estimated daily GPS orbits and EOPs. Indeed, simulated tidal aliases are found to be very similar to the observed IGS orbital features in the sub-seasonal bands. Moreover and unexpectedly, some low draconitic harmonics were also produced, potentially a source for the widespread errors in most IGS products. The results from this study are further evidence for the need of an improved sub-daily EOP tide model.     

Impact of loading displacements on SLR-derived parameters and on the consistency between GNSS and SLR results

Displacements of the Earth’s surface caused by tidal and non-tidal loading forces are relevant in high-precision space geodesy. Some of the corrections are recommended by the international scientific community to be applied at the observation level, e.g., ocean tidal loading (OTL) and atmospheric tidal loading (ATL). Non-tidal displacement corrections are in general recommended not to be applied in the products of the International Earth Rotation and Reference Systems Service, in particular atmospheric non-tidal loading (ANTL), oceanic and hydrological non-tidal corrections. We assess and compare the impact of OTL, ATL and ANTL on SLR-derived parameters by reprocessing 12 years of SLR data considering and ignoring individual corrections. We show that loading displacements have an influence not only on station long-term stability, but also on geocenter coordinates, Earth Rotation Parameters, and satellite orbits. Applying the loading corrections reduces the amplitudes of annual signals in the time series of geocenter and station coordinates. The general improvement of the SLR station 3D coordinate repeatability when applying OTL, ATL and ANTL corrections are 19.5 %, 0.2 % and 3.3 % respectively, w.r.t. the solutions without loading corrections. ANTL corrections play a crucial role in the combination of optical (SLR) and microwave (GNSS, VLBI, DORIS) space geodetic observation techniques, because of the so-called Blue-Sky effect: SLR measurements can be carried out only under cloudless sky conditions—typically during high air pressure conditions, when the Earth’s crust is deformed, whereas microwave observations are weather-independent. Thus, applying the loading corrections at the observation level improves SLR-derived products as well as the consistency with microwave-based results. We assess the Blue-Sky effect on SLR stations and the consistency improvement between GNSS and SLR solutions when ANTL corrections are included. The omission of ANTL corrections may lead to inconsistencies between SLR and GNSS solutions of up to 2.5 mm for inland stations. As a result, the estimated GNSS–SLR coordinate differences correspond better to the local ties at the co-located stations when applying ANTL corrections.     

Understanding the Spatial Variability of Chlorophyll a and Total Suspended Matter Distribution Along the Southwest Bay of Bengal Using In-Situ and OCM-2 & MODIS-Aqua Measurements

Spatial and temporal distribution of chlorophyll a (chl a) and Total Suspended Matter (TSM) and inter comparison of Ocean Color Monitor-2 (OCM-2) and Moderate Resolution Imaging Spectro-radiometer (MODIS-Aqua) derived chlorophyll a and TSM was made along the southwest Bay of Bengal (BoB). The in-situ chl a and TSM concentration measured during different seasons were ranged from 0.09 to 10.63 μgl^?1 and 11.04–43.75 mgl^?1 respectively. OCM-2 and MODIS derived chl a showed the maximum (6–8 μgl^?1) at nearshore waters and the minimum (0–1 μgl^?1) along the offshore waters. OCM-2 derived TSM imageries showed the maximum (50–60 mgl^?1) along the nearshore waters of Palk Strait and the moderate concentration (2–5 mgl^?1) was observed in the offshore waters. MODIS derived minimum TSM concentration (13.244 mgl^?1) was recorded along the offshore waters, while the maximum concentration of 15.78 mgl^?1 was found along the Kodiakarai region. The inter-comparison of OCM-2 and MODIS chl a data (R ² ?=?0.549, n?=?49, p?<?0.001, SEE?=?±0.117) indicate that MODIS data overestimates chl a concentration in the nearshore waters of the southern BoB compared to the OCM-2. The correlation between OCM-2 and MODIS-Aqua TSM data (R ² ?=?0.508, N?=?53, P?<?0.001 and SEE?=?±0.024) confirms that variation in the range of values measured by OCM-2 (2–60 mgl^?1) and the MODIS (13–16 mgl^?1) derived TSM values. Despite problems in range of measurements, persistent cloud cover etc., the launch of satellites like OCM-2 with relatively high spatial resolutions makes job easier and possible to monitor chl a distribution and sediment discharges on day to day basis in the southwest BoB.     

Negotiating on location, timing, duration, and participant in agent-mediated joint activity-travel scheduling

Agent-based simulation has become an important modeling approach in activity-travel analysis. Social activities account for a large amount of travel and have an important effect on activity-travel scheduling. Participants in joint activities usually have various options regarding location, participants, and timing and take different approaches to make their decisions. In this context, joint activity participation requires negotiation among agents involved, so that conflicts among the agents can be addressed. Existing mechanisms do not fully provide a solution when utility functions of agents are nonlinear and non-monotonic. Considering activity-travel scheduling in time and space as an application, we propose a novel negotiation approach, which takes into account these properties, such as continuous and discrete issues, and nonlinear and non-monotonic utility functions, by defining a concession strategy and a search mechanism. The results of experiments show that agents having these properties can negotiate efficiently. Furthermore, the negotiation procedure affects individuals’ choices of location, timing, duration, and participants.     

Semantically Enriching an Open Source Sensor Observation Service Implementation for Accessing Heterogeneous Environmental Data Sources

Many kinds of environmental data are nowadays publicly available, but spread over the web. This article discusses using the Sensor Observation Service (SOS) standard of the Open Geospatial Consortium (OGC) as a common interface for providing data from heterogeneous sources which can be integrated to a user tailored environmental information system. In order to allow for providing user‐tailored and problem‐specific information the adjusted SOS is augmented by a semantic layer which maps the environmental information to ontology concepts. The necessary information fusion from different domains and data types lead to several specific requirements for the SOS. Addressing these requirements we have implemented a SOS which still conforms to the OGC SOS 1.0.0 standard specification. The developed SOS has been integrated in a publicly available demonstrator of our personalized environmental information system. Additionally this article discusses future consequences for the SOS, caused by the recently published SOS 2.0 specification.     

Integrating a Raster Geographical Information System with Multi‐Objective Land Allocation Optimization for Conservation Reserve Design

The multi‐objective land allocation problem is to optimize the selection of land for different uses based on a set of decision objectives. For most applications, a geographical information system (GIS) is either absent or loosely coupled through file exchange. In this article the evolutionary algorithm (EA), a heuristic solution method for optimization problems, is integrated with a raster GIS to form a spatial decision support system (SDSS) for multi‐objective conservation reserve design. The SDSS effectively combines the functions of a GIS for data management, analysis, and visualization, with the optimization capability of the EA; and provides a uniform way to solve conservation reserve design problems with different types of constraints and objectives. The SDSS is demonstrated through application to the creation of conservation reserves in Bolivia to protect 17 endemic mammals.     

3-D Radargrammetric Modeling of RADARSAT-2 Ultrafine Mode: Preliminary Results of the Geometric Calibration

The geometry and the accuracy of the 3-D cartographic localization of RADARSAT-2 images are being evaluated as part of the Canadian Space Agency's Science and Operational Applications Research program. In a first step, the Toutin's 3-D physical model, previously developed for RADARSAT-1, was adapted to RADARSAT-2 sensor and applied to two ultrafine mode images (U2 and U25) acquired over an area in Beauport, Quebec. Both the 3-D modeling computed with only 12 ground control points and its geometric localization were evaluated with different check data: 1) independent check points; 2) the two quasi-epipolar images; 3) the two orthoimages; and 4) 1-m accurate orthophotos. All four results and validations are in agreement and confirm that the 3-D geometric localization and restitution accuracy are 1 m in planimetry and 2 m in elevation. The checked data error being included in these evaluations and the relative error computed from the quasi-epipolar comparison provided a high level of confidence that the precision of Toutin's 3-D radargrammetric model is better than 0.25 m.