Proposed support design, Kaletepe tunnel, Turkey

In this paper, preliminary support design of Kaletepe tunnel, located on Bilecik-Istanbul highway, Turkey, was analyzed by empirical and numerical methods. The rock mass rating (RMR) and rock mass quality (Q) systems were employed for empirical rock mass quality determination. Numerical analysis for the stress–strain distribution of the tunnel excavation and support systems was also carried out. The applied support performance was investigated at different sections of the tunnel route. It was seen that empirical and numerical approaches showed similar results. This indicates that when the empirical method is supported by numerical method, the preliminary support design will be more reliable.     

Three-dimensional forward and backward modelling of diapirism: numerical approach and its applicability to the evolution of salt structures in the Pricaspian basin

Numerical studies of ductile deformations induced by salt movements have, until now, been restricted to two-dimensional (2D) modelling of diapirism. This paper suggests a numerical approach to model the evolution of three-dimensional (3D) salt structures toward increasing maturity. This approach is also used here to restore the evolution of salt structures through successive earlier stages. The numerical methodology is applied to study several model examples. We analyse a model of salt diapirs that develop from an initial random perturbation of the interface between salt and its overburden and restore the evolved salt diapirs to their initial stages. We show that the average restoration errors are less than 1%. An evolutionary model of a 2D salt wall loaded by a 2D pile of sediments predicts a decomposition of the salt wall into 3D diapiric structures when the overburden of salt is supplied by 3D synkinematic wedge of sediments. We model salt extrusion feeding a gravity current over the depositional surface and estimate an average rate of extrusion and horizontal velocity of salt spreading. Faulting of the overburden to salt overhangs initiates new secondary diapirs, and we analyse the growth of these secondary diapirs. We also study how lateral flow effects the evolution of salt diapirs. The shape of a salt diapir can be very different if the rate of horizontal flow is much greater than the initial rate of diapiric growth solely due to gravity. We discuss the applicability of the results of the models to the evolution of Late Permian salt structures in the Pricaspian basin (Russia and Kazakhstan). These structures are distinguishable into a variety of styles representing different stages of growth: salt pillows, diapirs, giant salt massifs, 2D diapiric walls and 3D stocks complicated by large overhangs. The different sizes, shapes and maturities of salt structures in different parts of the Pricaspian basin reflect areal differences in salt thickness and loading history. Our results suggest that the numerical methodology can be employed to analyse the evolution of all salt structures that have upbuilt through younger ductile overburdens.     

Crustal thickening in an extensional regime: application to the mid-Norwegian Vøring margin

The structure of the mid-Norwegian volcanic Vøring margin at the onset of the Maastrichtian–Paleocene extension phase reflects the cumulative effect of earlier consecutive rifting events. Lateral structural differences present on the margin at that time are a consequence of migration of the location of maximum extension in time between Norway and Greenland. The most important imprints (Moho depth, thermal structure) of these events on the lithosphere are incorporated in a numerical simulation of the final extension phase. We focus on a possible mechanism of formation of the Vøring Marginal High and address the relationship between spatial and temporal evolution of crustal thinning and thickening, uplift of the surface and strength of the lithosphere.It is found that the Vøring Basin formed the strongest part of the margin which explains why the Maastrichtian–Paleocene rift axis was not located here but instead jumped westward with respect to the earlier rift axes locations. The modeling study predicts that local crustal thickening during extension can be expected when large lateral thermal variations are present in the lithosphere at the onset of extension. Negative buoyancy induced by lateral temperature differences increases downwelling adjacent to the rifting zone; convergence of material at the particular part of the margin is mainly taken up by the lower crust. The model shows that during the final phase of extension, the crust in the Vøring Marginal High area was thickened and the surface uplifted. It is likely that this dynamic process and the effects of magmatic intrusions both acted in concert to form the Marginal High.     

黄土震陷灾害典型震例的综合研究

阐述了1995年甘肃永登5.8级地震时造成的黄土震陷灾害形成条件、震害形态特征、黄土的基本性质及湿陷性和静、动土力学参数与震陷灾害的成因关系。并用震陷量预测计算方法和三维有限元数值模拟方法检验了黄土震陷灾害预测的理论、方法并解释了永登5.8 级地震黄土震陷现象。     

On the variability of the flow along the Meso-American Barrier Reef system: a numerical model study of the influence of the Caribbean current and eddies

A high resolution (3–8 km grid), 3D numerical ocean model of the West Caribbean Sea (WCS) is used to investigate the variability and the forcing of flows near the Meso-American Barrier Reef System (MBRS) which runs along the coasts of Mexico, Belize, Guatemala and Honduras. Mesoscale variations in velocity and temperature along the reef were found in seasonal model simulations and in observations; these variations are associated with meandering of the Caribbean current (CC) and the propagation of Caribbean eddies. Diagnostic calculations and a simple assimilation technique are combined to infer the dynamically adjusted flow associated with particular eddies. The results demonstrate that when a cyclonic eddy (negative sea surface height anomaly (SSHA)) is found near the MBRS the CC shifts offshore, the cyclonic circulation in the Gulf of Honduras (GOH) intensifies, and a strong southward flow results along the reef. However, when an anticyclonic eddy (positive SSHA) is found near the reef, the CC moves onshore and the flow is predominantly westward across the reef. The model results help to explain how drifters are able to propagate in a direction opposite to the mean circulation when eddies cause a reversal of the coastal circulation. The effect of including the Meso-American Lagoon west of the Belize Reef in the model topography was also investigated, to show the importance of having accurate coastal topography in determining the variations of transports across the MBRS. The variations found in transports across the MBRS (on seasonal and mesoscale time scales) may have important consequences for biological activities along the reef such as spawning aggregations; better understanding the nature of these variations will help ongoing efforts in coral reef conservation and maintaining the health of the ecosystem in the region.     

The role of the seismic trigger in the Calitri landslide (Italy): historical reconstruction and dynamic analysis

The study covered by this paper was focused on the historical case of the Calitri landslide, which was repeatedly reactivated by earthquakes, as reported since 1694. The town of Calitri (Southern Italy) is located on a ridge whose southern slope, from its top to the Ofanto river valley floor, has been historically affected by major landsliding. The last record of recurrence of the Calitri landslide leads back to the 1980 Irpinia earthquake, which caused significant damage to the town and had pervasive and visible ground effects. Based on a detailed historical reconstruction of landsliding and seismicity at Calitri, the study analysed the current static and dynamic stability of the landsliding slope by means of a finite-difference numerical analysis, taking into account the various factors of landslide initiation and kinematics. The results of the numerical analysis: (i) were consistent with the roto-translational mechanism observed upon the latest reactivation of the landslide; (ii) demonstrated that excess pore pressure redistribution caused a lag between the seismic trigger and the initiation of landsliding; and (iii) showed the impact of seismic input frequency on propagation and depth of slope instability.     

Western Alpine back-thrusting as subsidence mechanism in the Tertiary Piedmont Basin (Western Po Plain, NW Italy)

Basin formation dynamics of the Tertiary Piedmont Basin (TPB) are here investigated by means of cross-section numerical modelling. Previous works hypothesised that basin subsidence occurred due first to extension (Oligocene) and then to subsequent loading due to back-thrusting (Miocene). However, structural evidence shows that the TPB was mainly under contraction from Oligocene until post Pliocene time while extension played a minor role. Furthermore, thermal indicators strongly call for a cold (flexure-induced) mechanism but are strictly inconsistent with a hot (thermally induced) mechanism. Our new modelling shows that the TPB stratigraphic features can be reproduced by flexure of a visco-elastic plate loaded by back-thrusts active in the Western Alps in Oligo-Miocene times. Far-field compression contributed to the TPB subsidence and controlled the basin infill geometry by enhancing basin tilting, forebulge uplift and erosion of the southern margin of the basin. These results suggest that the TPB subsidence is the result of a combination of mechanisms including thrust loading and far-field compressional stresses.     

From the geological to the numerical model in the analysis of gravity-induced slope deformations: An example from the Central Apennines (Italy)

This paper presents the findings from a study on gravity-induced slope deformations along the northern slope of Mt. Nuria (Rieti-Italy). The slope extends from the village of Pendenza to the San Vittorino plain and hosts the Peschiera River springs, i.e. the most important springs of the Central Apennines (average discharge: about 18 m³/s).

Detailed geological-geomorphological and geomechanical surveys, supported by a site stress-strain monitoring system and laboratory tests, led us to define the main evolutionary features of the studied phenomena. Based on the collected data, a “geological-evolutionary model” was developed with a view to identifying a spatio-temporal correlation between relief forms, jointing of the rock mass and its stress conditions. The geological-evolutionary model was expected to improve numerical simulations and to test our assumptions.

The numerical model also allowed us to simulate changes in the stress-strain conditions of the rock mass and correlate them with jointing, seepage, as well as with site-detected and site-monitored forms and deformations. In particular, significant relations between seepage, tensile stresses within the rock mass, karst solution and collapse of cavities were identified.     

Universal hydrofracturing algorithm for shear-thinning fluids: Particle velocity based simulation

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.     

Flexible membranes anchored to the ground for slope stabilisation: Numerical modelling of soil slopes using SPH

An alternative modelling for flexible membranes anchored to the ground for soil slope stabilisation is presented using Smoothed-Particle Hydrodynamics to model the unstable ground mass in a soil slope, employing a dynamic solve engine. A regression model of pressure normal to the ground, q_sim, and also membrane deflection, f_sim, have been developed using Design of Experiment. Finally, a comparison between the pressure obtained from numerical simulation and from a limit equilibrium analysis considering infinite slope has been carried out, showing differences in the results, mainly due to the membrane stiffness.