Sources of recent tectonic stress in the Pannonian region:inferences from finite element modelling

We present the results of finite element modelling of the recent stress field in the Pannonian basin and surrounding Alpine orogenic belt. Our results show that the recent, predominantly compressive, stress regime in the Alpine–Pannonian–Carpathian–Dinaric system is governed by distinct tectonic factors. Of great importance is the deformation of crustal blocks with different geometries and rigidities in an overall convergent setting associated with the Africa–Europe collision. The most important stress source appears to be the counterclockwise rotation of the Adriatic microplate at the southwest boundary of the Pannonian basin. This plate tectonic unit has been interpreted as moving independently of both the European plate and the African plate. Additional boundary conditions—active shortening and compression in the Vrancea zone and the Bohemian Massif, and the effect of the immobile Moesian Platform—also significantly influence the modelling results. The incorporation of additional stress sources such as crustal thickness variation and the presence of two main fault zones separating the primary tectonic units in the study area have only locally important effects but improve the fit between the calculated results and the observed stress pattern.     

A pseudo-coupled numerical approach for stability analysis of frozen soil slopes based on finite element limit analysis method

To simplify the stability analysis of frozen soil slope, a pseudo-coupled numerical approach is developed. In this approach, the coupled heat transfer and water flow in frozen soils are simulated first, and based on the computed thermal-hydro field, the stability of frozen soil slope is evaluated. Although the shear strength for frozen soil is very complicated and is usually represented by a nonlinear MC failure criterion, a simple linear MC yield criterion is utilized. In this method, the internal friction angle is expressed as a function of volumetric ice content and the cohesion is fitted as a simple bilinear expression of T and volumetric water content. To assess slope stability, the limit analysis is employed in conjunction with the recently developed α-section search algorithm. A frozen soil slope example is used to examine the proposed pseudo-coupled numerical approach, and numerical studies validate its effectiveness. Based on numerical results, it is seen that slope stability may be remarkably influenced by warming air (or ground surface) temperature. With increasing ground surface temperature, slope stability indicated by FOS may reduce to 1.0, implying that warming air temperature could be a trigger of frozen soil slope failure.     

A finite element analysis of tidal deformation of the entire earth with a discontinuous outer layer

Tidal deformation of the Earth is normally calculated using the analytical solution with some simplified assumptions, such as the Earth is a perfect sphere of continuous media. This paper proposes an alternative way, in which the Earth crust is discontinuous along its boundaries, to calculate the tidal deformation using a finite element method. An in-house finite element code is firstly introduced in brief and then extended here to calculate the tidal deformation. The tidal deformation of the Earth due to the Moon was calculated for an geophysical earth model with the discontinuous outer layer and compared with the continuous case. The preliminary results indicate that the discontinuity could have different effects on the tidal deformation in the local zone around the fault, but almost no effects on both the locations far from the fault and the global deformation amplitude of the Earth. The localized deformation amplitude seems to depend much on the relative orientation between the fault strike direction and the loading direction (i.e. the location of the Moon) and the physical property of the fault.     

Probabilistic analysis of embankment slope stability in frozen ground regions based on random finite element method

Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafrost embankment designin frozen ground regions. The prediction may be further complicated due to the inherent uncertainties of materialproperties. Hence, stochastic analyses should be conducted. Firstly, Karhunen-Loeve expansion is applied to attain the random fields for hydraulic and thermal conductions. Next, the mixed-form modified Richards equation for mass transfer (i.e., mass equation) and the heat transport equation for heat transient flow in a variably saturated frozen soil are combined into one equation with temperature unknown. Furthermore, the finite element formulation for the coupled thermal-hydraulic fields is derived. Based on the random fields, the stochastic finite element analyses on stability of embankment are carried out. Numerical results show that stochastic analyses of embankment stability may provide a more rational picture for the distribution of factors of safety (FOS), which is definitely useful forembankment design in frozen ground regions.     

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method.
  Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately.
  Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.     

Modelling the temporal element in land information systems

Abstract

This paper considers the use of a modern programming language (Modula-2) to develop a data model for a lot-based land information system. The emphasis is on the importance of maintaining the history of the lots, and a data model is developed which incorporates the history of each lot.     

Surface joint systems, tectonic stresses and geomorphology: a reconciliation of conflicting observations

This note seeks to reconcile the widespread small-scale fractographic observations indicating that surface joints are extension fractures with the frequently observed occurrence of such joints as what appear to be nearly orthogonal conjugate sets, their strikes fitting the shear lines of the neotectonic stress field. It is suggested that the local appearance of joints as extension fractures may have nothing to do with their orientation in the large-scale neotectonic stress field, inasmuch as extrapolations from a local scale to plate-tectonic dimensions are quite speculative. Thus, the local extension-characteristics of the joint surfaces may be acquired at the latest stage of their genesis as a result of the corresponding rock faces becoming exposed, but the geometrically orientational attributes may be conditioned by the shear in the surrounding large-scale neotectonic stress field; a new possible mechanism for reconciling the conflicting local and large-scale observations is suggested.