Automated mapping of rock slope geometry,kinematics and stability with RSS-GIS

A GIS-implemented, deterministic approach for the automated spatial evaluation of geometrical and kinematical properties of rock slope terrains is presented. Based on spatially distributed directional information on planar geological fabrics and DEM-derived topographic attribute data, the internal geometry of rock slopes can be characterized on a grid cell basis. For such computations, different approaches for the analysis and regionalization of available structural directional information applicable in specific tectonic settings are demonstrated and implemented in a GIS environment. Simple kinematical testing procedures based on feasibility criteria can be conducted on a pixel basis to determine which failure mechanisms are likely to occur at particular terrain locations. In combination with hydraulic and strength data on geological discontinuities, scenario-based rock slope stability evaluations can be performed. For conceptual investigations on rock slope failure processes, a GIS-based specification tool for a 2-D distinct element code (UDEC) was designed to operate with the GIS-encoded spatially distributed rock slope data. The concepts of the proposed methodology for rock slope hazard assessments are demonstrated at three different test sites in Germany.     

Robust biased estimation and its applications in geodetic adjustments

In this paper, a number of robust biased estimators (e.g. ordinary robust ridge estimator, robust principal components estimator, robust combined principal components estimator, robust single-parametric principal components estimator, robust root-root estimator) are established by means of a unified expression of biased estimators and based on the principle of equivalent weight. The most attractive advantage of these new estimators is that they can not only overcome the ill-conditioning of the normal equation but also have the ability to resist outliers. A numerical example is used to illustrate that these new estimators are much better than the least-squares estimator and various biased estimators even when both ill-conditioning and outliers exist. Received: 14 November 1995/Accepted: 11 February 1998     

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.     

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.     

The determination of potential difference by the joint application of measured and synthetical gravity data: a case study in Hungary

In an elementary approach every geometrical height difference between the staff points of a levelling line should have a corresponding average g value for the determination of potential difference in the Earth’s gravity field. In practice this condition requires as many gravity data as the number of staff points if linear variation of g is assumed between them. Because of the expensive fieldwork, the necessary data should be supplied from different sources. This study proposes an alternative solution, which is proved at a test bed located in the Mecsek Mountains, Southwest Hungary, where a detailed gravity survey, as dense as the staff point density (~1 point/34 m), is available along a 4.3-km-long levelling line. In the first part of the paper the effect of point density of gravity data on the accuracy of potential difference is investigated. The average g value is simply derived from two neighbouring g measurements along the levelling line, which are incrementally decimated in the consecutive turns of processing. The results show that the error of the potential difference between the endpoints of the line exceeds 0.1 mm in terms of length unit if the sampling distance is greater than 2 km. Thereafter, a suitable method for the densification of the decimated g measurements is provided. It is based on forward gravity modelling utilising a high-resolution digital terrain model, the normal gravity and the complete Bouguer anomalies. The test shows that the error is only in the order of 10⁻³mm even if the sampling distance of g measurements is 4 km. As a component of the error sources of levelling, the ambiguity of the levelled height difference which is the Euclidean distance between the inclined equipotential surfaces is also investigated. Although its effect accumulated along the test line is almost zero, it reaches 0.15 mm in a 1-km-long intermediate section of the line.