Coeval folding and boudinage under plane strain with the axis of no change perpendicular to the layer

Plane-strain coaxial deformation of a competent plasticine layer embedded in an incompetent plasticine matrix was carried out to improve our understanding about the evolution of folds and boudins if the layer is oriented perpendicular to the Y-axis of the finite strain ellipsoid. The rock analogues used were Beck’s green plasticine (matrix) and Beck’s black plasticine (competent layer), both of which are strain-rate softening modelling materials with a stress exponent n=ca. 8. The effective viscosity η of the matrix plasticine was changed by adding different amounts of oil to the original plasticine. At a strain rate of 10⁻³ s⁻¹ and a finite strain e of 10%, the effective viscosity of the matrix ranges from 1.2×10⁶ to 7.2×10⁶ Pa s. The effective viscosity of the competent layer has been determined as 4.2×10⁷ Pa s. If the viscosity ratio is large (ca. 20) and the initial thickness of the competent layer is small, both folds and boudins develop simultaneously. Although the growth rate of the folds seems to be higher than the growth rate of the boudins, the wavelength of both structures is approximately the same as is suggested by analytical solutions. A further unexpected, but characteristic, aspect of the deformed competent layer is a significant increase in thickness, which can be used to distinguish plane-strain folds and boudins from constrictional folds and boudins.An erratum to this article can be found at     

地震型砾岩崩塌体运动学特征分析

汶川地震在工程区诱发了大量的崩塌灾害,提供了大量的"天然崩塌试验"现象。本文以汶川地震诱发的都江堰青城后山红岩度假区的大型砾岩崩塌体作为研究对象,主要分析了它的运动学特征。通过测量获得崩塌体最大位移块石运移边界角(reach of angle)。结合现场调查和数值模拟,得出坡体植被及人工构筑物的遮挡作用是导致研究对象运移边界角数值偏大的主要因素。统计该崩塌堆积体沿沟谷延伸方向的块石分布情况,并分析沟谷拐点处块石大小比例及各大小块石数占统计区相应大小块石总数的比例。得出大块石(>1m3)在滚落过程中,平均约有总数量的1/4被截留在了拐点处。作者认为,充分分析自然力作用下的典型崩塌体的分布及运移等特征对于崩塌的运动学研究及人类的防灾减灾工程具有重要的参考和指导。     

Direct dating of folding events by 40Ar/39Ar analysis of synkinematic muscovite from flexural-slip planes

Timing of folding is usually dated indirectly, with limited isotopic dating studies reported in the literature. The present study investigated the timing of intracontinental, multi-stage folding in Upper Proterozoic sandstone, limestone, and marble near Beijing, North China, and adjacent regions. Detailed field investigations with microstructural, backscattered electron (BSE) images and electron microprobe analyses indicate that authigenic muscovite and sericite crystallized parallel to stretching lineations/striations or along thin flexural-slip surfaces, both developed during the complex deformation history of the study area, involving repeated compressional, extensional and strike-slip episodes. Muscovite/sericite separates from interlayer-slip surfaces along the limbs and from dilatant sites in the hinges of folded sandstones yield muscovite ⁴⁰Ar/³⁹Ar plateau ages of ∼158–159 Ma, whereas those from folded marble and limestone samples yield ages of 156 ± 1 Ma. Muscovite from thin flexural-slip planes on fold limbs and hinges yields ages within analytical error of ∼155–165 Ma. Further muscovite samples collected from extensionally folded limestone and strike-slip drag folds yield younger ages of 128–125 Ma with well-defined plateaus. To assess the potential influence of the detrital mica component of the host rock on the age data, two additional muscovite samples were investigated, one from a folded upper Proterozoic–Cambrian sandstone outside the Western Hills of Beijing and one from a folded sandstone sampled 20 cm from folding-related slip planes. Muscovite separates from these samples yield significantly older ages of 575 ± 2 Ma and 587 ± 2 Ma, suggesting that the timing of folding can be directly determined using the ⁴⁰Ar/³⁹Ar method. This approach enables the identification and dating of distinct deformation events that occur during multi-stage regional folding. ⁴⁰Ar/³⁹Ar dating can be used to constrain the timing of muscovite and sericite growth at moderate to low temperatures (<400 °C) during folding, yielding well-defined plateau ages and thereby the age of deformation in the upper crust.     

Structure of the rift basins in the central Gulf of California: Kinematic implications for oblique rifting

We use structural and seismostratigraphic interpretation of multichannel seismic reflection data to understand the structure and kinematic history of the central Gulf of California. Our analysis reveals that oblique strain in the central Gulf formed two tectono–sedimentary domains during distinct deformation stages. The eastern domain, offshore Sonora, is bounded by the East and West Pedro Nolasco faults that may constitute the southernmost segments of the Tiburón Fault System. Within this domain, the dip-slip Yaqui Fault controlled deposition of 3.9 km of sediments in the half-graben Yaqui Basin. The western domain, offshore Baja California, is bounded by the Guaymas Transform Fault, which controlled the accumulation of 1.45 km of sediments within a half-graben that formed the early Guaymas Basin. The tectono–sedimentary activity offshore Sonoran likely ranges from Late Miocene–Pliocene to Late Pliocene time, while activity in the Guaymas Basin commenced in Late Pliocene time. Extinction of the main faults offshore Sonora was nearly coeval to the initiation of the Guaymas Transform Fault. Our results suggest that oblique strain has been accommodated by strain partition since the onset of rifting in the central Gulf. The Guaymas Basin is now a nascent spreading center, but prior to this, it evolved as a half-graben controlled by the Guaymas Transform Fault; such drastic transition is not constrained, but likely occurred during the Pleistocene time and must be localized < 30 km north of the axial troughs. The faults within the central Gulf transpose the Miocene N–S oriented grabens of Basin and Range style preserved onshore in the conjugate rifted margins.     

Role of lateral thickness variations on the development of oblique structures at the Western end of the South Pyrenean Central Unit

The Montsec unit is one of the most important detached South-verging nappes within the South Pyrenean Central Unit (SPCU, Southern Pyrenees). A N–S cross-section of its Western sector, based on seismic reflection profiles, shows a hangingwall ramp geometry in Mesozoic strata, overlain by a syntectonic series of Lower Eocene sediments with growth geometry. The geometry of growth strata constrains the age of its movement between the Paleocene and the Middle Eocene. The geometry of the Western, oblique ramp of the South Pyrenean Central Unit is defined by a series of N–S folds, in some cases associated with underlying West-verging thrusts, as indicated by seismic reflection profiles and field data. In this paper, we propose that the geometry of the thrust wedge of Mesozoic units, progressively thinning from East to West, strongly contributed to constrain the location and geometry of the Western termination of the Montsec thrust. The hypothesis proposed is checked by a series of experimental wedges developed in a sandpack with lateral and three-dimensional thickness variations. Oblique structures form as thrusting progresses at the tip of the sand wedge.     

Intergranular diffusion rates from the analysis of garnet surfaces: Implications for metamorphic equilibration

Novel approaches to garnet analysis have been used to assess rates of intergranular diffusion between different matrix phases and garnet porphyroblasts in a regionally metamorphosed staurolite‐mica‐schist from the Barrovian‐type area in Scotland. X‐ray maps and chemical traverses of planar porphyroblast surfaces reveal chemical heterogeneity of the garnet grain boundary linked to the nature of the adjacent matrix phase. The garnet preserves evidence of low temperature retrograde exchange with matrix minerals and diffusion profiles documenting cation movement along the garnet boundaries. Garnet–quartz and garnet–plagioclase boundaries preserve evidence of sluggish Mg, Mn and Fe diffusion at comparable rates to volume diffusion in garnet, whereas diffusion along garnet–biotite interfaces is much more effective. Evidence of particularly slow Al transport, probably coupled to Fe³⁺ exchange, is locally preserved on garnet surfaces adjacent to Fe‐oxide phases. The Ca distribution on the garnet surface shows the most complex behaviour, with long‐wavelength heterogeneities apparently unrelated to the matrix grain boundaries. This implies that the Ca content of garnet is controlled by local availability and is thought likely to reflect disequilibrium established during garnet growth. Geochemical anomalies on the garnet surfaces are also linked to the location of triple junctions between the porphyroblasts and the matrix phases, and imply enhanced transport along these channels. The slow rates of intergranular diffusion and the characteristics of different boundary types may explain many features associated with the prograde growth of garnet porphyroblasts. Thus, minerals such as quartz, Fe‐oxides and plagioclase whose boundaries with garnet are characterized by slow intergranular diffusion rates appear to be preferentially trapped as inclusions within porphyroblasts. As such grain boundary diffusion rates may be a significant kinetic impediment to metamorphic equilibrium and garnet may struggle to maintain chemical and textural equilibrium during growth in pelites.     

Strain analysis and vorticity of flow in the Northern Sardinian Variscan Belt: Recognition of a partitioned oblique deformation event

A field example of strain partitioning has been analysed along the Nurra–Asinara transect of the NW Sardinian Variscan chain (Italy). The section in the Nurra–Asinara area is in a continuous sequence of tectono-metamorphic complexes made of low- to high-grade metamorphic rocks affected by a polyphase tectonic history. The principal fabric of the area is controlled by a D2 progressive deformation phase in which the strain is partitioned into folds and shear zone domains. The D2 stretching lineation and shear sense show a clear change from south to north. The principal meso- and micro-structures, vorticity gauges and a quantitative kinematic analysis of local strain suggest that the D2 kinematic history could be envisaged as an oblique heterogeneous deformation similar to the transpressive systems described in ancient and modern settings elsewhere. Using a simple kinematic model we also propose that both a transpressive system followed by “thrusting” or a partitioned transpressive system could be responsible for the fabric distribution and strain accumulation described in the study transect.     

An oblique factor analysis solution for the analysis of mixtures

The problems involved in the factor analysis of data consisting of measurements of material quantity or proportion are discussed, and the inability of existing factor analysis solutions to properly model such data is pointed out. A new factor analysis solution using the linear programming technique is presented which properly analyzes mixture data. A numerical example is presented in which a body of suspended sediment grain-size distribution data is analyzed using the new factor analysis solution. FORTRAN-IV subroutines for computing this solution are included in an Appendix.     

An attempt to use LA-ICP-SMS to quantify enrichment of trace elements on pyrite surfaces in oxidizing mine tailings

Metals released from oxidation and weathering of sulphide minerals in mine tailings are to a high degree retained at deeper levels within the tailings themselves. To be able to predict what could happen in the future with these secondarily retained metals, it is important to understand the retention mechanisms. In this study an attempt to use laser ablation high-resolution ICP-MS (LA-ICP-SMS) to quantify enrichment of trace elements on pyrite surfaces in mine tailings was performed. Pyrite grains were collected from a profile through the pyrite-rich tailings at the Kristineberg mine in northern Sweden. At each spot hit by the laser, the surface layer was analyzed in the first shot, and a second shot on the same spot gave the chemical composition of the pyrite immediately below. The crater diameter for a laser shot was known, and by estimating the crater depth and total pyrite surface, the total enrichment on pyrite grains was calculated. Results are presented for As, Cd, Co, Cu, Ni and Zn. The results clearly show that there was an enrichment of As, Cd, Cu and Zn on the pyrite surfaces below the oxidation front in the tailings, but not of Co and Ni. Arsenic was also enriched on the pyrite grains that survived in the oxidized zone. Copper has been enriched on pyrite surfaces in unoxidized tailings in the largest amount, followed by Zn and As. However, only 1.4 to 3.1% of the Cd and Zn released by sulphide oxidation in the oxidized zone have been enriched on the pyrite surfaces in the unoxidized tailings, but for As and Cu corresponding figures are about 64 and 43%, respectively. There were many uncertainties in these calculations, and the results shall not be taken too literally but allowed the conclusion that enrichment on pyrite surfaces is an important process for retention of As and Cu below the oxidation front in pyrite rich tailings. Laser ablation is not a surface analysis technique, but more of a thin layer method, and gives no information on the type of processes resulting in enrichment on the pyrite surfaces. Although only pyrite grains that appeared to be fresh and without surface coatings were used in this study, the possibility that a thin layer of Fe-hydroxides occurred must be considered. Both adsorption to the pyrite directly or to Fe-oxyhydroxides may explain the enrichment of As, Cd, Cu and Zn on the pyrite surfaces, and, in the case of Cu, also the replacement of Fe(II) by Cu(II) in pyrite.     

青藏高原东缘鲜水河断裂带磁组构特征及构造意义

为了探究鲜水河断裂带的几何学、运动学特征,在野外构造、显微构造分析基础上,研究了鲜水河断裂带296块构造岩定向样品的磁组构特征和热磁特征,结果显示样品的平均磁化率km值总体较小,属微弱磁性到弱磁性;热磁实验及其显微构造表明顺磁性的页硅酸盐(如黑云母)等矿物对糜棱岩类样品磁化率贡献较大;磁化率各向异性度PJ总体较大,表明鲜水河断裂带构造变形强烈;磁化率椭球体形状参数T总体大于0,扁率E总体在1附近分布,说明鲜水河断裂带磁化率椭球体以扁圆形为主,整体上磁面理较磁线理发育,进一步显示出鲜水河断裂带构造变形样式以剪切、压扁为主,伴有拉伸的构造变形,同时也反映出鲜水河断裂带多次变形的综合特征;最小磁化率主轴Kmin方位表明鲜水河断裂带北段和南段分别受近EW向和NE-SW向主压应力控制;同时Kmin方位及其倾伏角特征显示鲜水河断裂带总体以左行走滑剪切为主,北段两侧块体在鲜水河断裂带两次不同的构造活动时期,各自有一定的相对抬升,但抬升幅度均不大;南段则是SW块体相对NE块体抬升,抬升幅度较大;整个断裂带特征显示出鲜水河断裂带在向南发展逐渐转化为挤压构造,这可能与青藏高原物质的向东逃逸受阻以及鲜水河断裂带与龙门山断裂带在此交接的地质背景不无关系.     

The transition from upright to recumbent folding in the Variscan fold belt of southwest England: a model based on the kinematics of simple shear

A transition from upright folds, at high structural levels, to recumbent folds at depth is described from the Variscan fold belt in southwest England. The folds tighten and cleavage intensifies progressively as the axial plane dip decreases. A simple shear model is developed in which the shortening of a multilayer and its folding produces initially upright open folds which tighten as they rotate during increasing shear strain. The model predicts the observed relationship between interlimb angle and axial plane dip and is used to discuss the development of the structure of north Cornwall.     

破碎岩体压实特性的三维离散元数值计算方法研究

研究破碎岩体的压实特性是矿井地下工程的基础工作之一,由于破碎岩体所处环境的隐蔽性与危险性,常采用实验室测试和数值计算的研究方法。提出了一种三维破碎岩体模型构建方法,即在3D Voronoi建立完整岩体数值模型的基础上,通过预定孔隙率,随机删除完整岩体中的块体反演破碎岩体结构,测定破碎岩体的压实特性。该方法可较真实地反映破碎岩体的块度特征、碎胀特性与压实特性,与现有研究方法有较高的吻合度,为矿山地下工程的安全控制提供了新的有效的研究方法。     

Preservation of cross-strata due to the migration of subaqueous dunes: an experimental investigation

This experimental investigation examined the controls on the geometry of cross‐sets formed by subaqueous dunes. A range of steady, unidirectional flow conditions spanning the field of dune existence was investigated, and aggradation rate ranged from 0 mm s^?1 to 0·014 mm s^?1. Data from an ultrasonic depth profiler consist of high‐resolution temporal and spatial series of bed profiles from which dune height and length, migration rate and the depth of trough scour were measured. Cross‐set thickness and length were measured from sediment peels. The size and shape of dunes from an equilibrium assemblage change continuously. Individual dunes commonly increase in height by trough scouring and, occasionally, by being caught‐up by the upstream dune. Both types of behaviour occur suddenly and irregularly in time and, hence, do not appear to depend on dunes further upstream. However, dune climbing or flattening is a typical response of dunes that disappear under the influence of the upstream dune. All types of behaviour occur at any flow velocity or aggradation rate. Successive dune‐trough trajectories, defined by dunes showing various behaviours, affect the geometry of the preserved cross‐sets. Mean cross‐set thickness/mean dune height averages 0·33 (±0·7), and mean cross‐set length/mean dune length averages 0·49 (±0·08), and both show no systematic variation with aggradation rate or flow velocity. Mean cross‐set thickness/mean cross‐set length tends to decrease with increasing flow velocity and Froude number, therefore allowing a qualitative estimation of flow conditions. Quantitative analysis of the temporal changes in the geometry and migration rate of individual dunes allows the development of a two‐dimensional stochastic model of dune migration and formation of cross‐sets. Computer realizations produced stacks of cross‐sets of comparable shape and thickness to laboratory flume observations, indicating a good empirical understanding of the variability of dune‐trough trajectories. However, interactions among dunes and aggradation rates of the order of 10^?2 mm s^?1 should be considered in future improved models.     

From discrete to continuum modelling of boundary value problems in geomechanics: An integrated FEM-DEM approach

Double-scale numerical methods constitute an effective tool for simultaneously representing the complex nature of geomaterials and treating real-scale engineering problems such as a tunnel excavation or a pressuremetre at a reasonable numerical cost. This paper presents an approach coupling discrete elements (DEM) at the microscale with finite elements (FEM) at the macroscale. In this approach, a DEM-based numerical constitutive law is embedded into a standard FEM formulation. In this regard, an exhaustive discussion is presented on how a 2D/3D granular assembly can be used to generate, step by step along the overall computation process, a consistent Numerically Homogenised Law. The paper also focuses on some recent developments including a comprehensive discussion of the efficiency of Newton-like operators, the introduction of a regularisation technique at the macroscale by means of a second gradient framework, and the development of parallelisation techniques to alleviate the computational cost of the proposed approach. Some real-scale problems taking into account the material spatial variability are illustrated, proving the numerical efficiency of the proposed approach and the benefit of a particle-based strategy.     

Optical ages on loess derived from outwash surfaces constrain the advance of the Laurentide Ice Sheet out of the Lake Superior Basin,USA

We present textural and thickness data on loess from 125 upland sites in west-central Wisconsin, which confirm that most of this loess was derived from the sandy outwash surfaces of the Chippewa River and its tributaries, which drained the Chippewa Lobe of the Laurentide front during the Wisconsin glaciation (MIS 2). On bedrock uplands southeast of the widest outwash surfaces in the Chippewa River valley, this loess attains thicknesses > 5 m. OSL ages on this loess constrain the advance of the Laurentide ice from the Lake Superior basin and into west-central Wisconsin, at which time its meltwater started flowing down the Chippewa drainage. The oldest MAR OSL age, 23.8 ka, from basal loess on bedrock, agrees with the established, but otherwise weakly constrained, regional glacial chronology. Basal ages from four other sites range from 13.2 to 18.5 ka, pointing to the likelihood that these sites remained geomorphically unstable and did not accumulate loess until considerably later in the loess depositional interval. Other OSL ages from this loess, taken higher in the stratigraphic column but below the depth of pedoturbation, range to nearly 13 ka, suggesting that the Chippewa River valley may have remained a loess source for several millennia.     

Structural analysis of a newly emerged accretionary prism along the Jinlun-Taimali coast,southeastern Taiwan: From subduction to arc-continent collision

In southern Taiwan the initial collision of the Luzon volcanic arc with the passive continental margin of China results in the emergence of an accretionary prism of, predominantly, turbidites in composition, thus providing an appropriate place to study the temporal and spatial variation of deformation during the transition of subduction to arc-continent collision. Field surveys have recently been carried out in slightly metamorphosed rocks along the well-exposed Jinlun-Taimali coast in southeastern Taiwan. Three folding phases are identified in the area. The first phase is characterized by gently dipping but widely distributed phyllitic cleavage (S₁). The second phase is represented by sparsely distributed crenulation cleavage (S₂) that folded the phyllitic cleavage. The third phase is characterized by E–W trending antiforms (F₃) that involved both types of pre-existing cleavages. Restoration of such an antiform in the north using a method proposed in this paper reveals that phyllitic cleavage in the overturned beds dips gently towards the southeast or east-southeast before the antiform, in relation to the first-phase thrusting or folding under regional ESE-WNW compression. From the first to third phase, the maximum horizontal compression underwent an about 90° anticlockwise rotation from ESE-WNW to E–W or NE–SW to N–S, and the deformation depth seems to decrease drastically, in terms of the decreasing proportion of pervasive deformation. All these variations are attributed to the oblique arc-continent collision that exhumed the whole accretionary prism and induced a local stress perturbation in southeastern Taiwan.     

河流相储层综合预测方法探索

本文以储层横向预测的研究内容为线索,从储层追踪、储层特征参数的提取与结构重建等方面入手,针对河流相储层薄、变化快且泥质干扰严重等特点,在立足于综合研究的基础上,特别注重波阻抗技术、地质统计学、多元统计等在河流相储层预测中的作用。目前,利用波阻抗剖面,采用钻井标定色标的办法,进行储层追踪是描述储层宏观特性的有利手段。虽测井信息与地震信息就其采样的储集空间来说处于两个极端,但应用统计学方法,如多元统计、地质统计学等,可以把这两种处于两个极端的信息综合起来。其具体实现是依据岩心、测井,参照井间地震建立储层变化模型,并借助于这一模型用合适的统计方法将这些不同分辨尺度的信息有机地融合在一起。     

Microfabrics of calcite ultramylonites as records of coaxial and non-coaxial deformation kinematics: Examples from the Rocher de l'Yret shear zone (Western Alps)

Microfabrics were analysed in calcite mylonites from the rim of the Pelvoux massif (Western Alps, France). WNW-directed emplacement of the internal Penninic units onto the Dauphinois domain produced intense deformation of an Eocene-age nummulitic limestone under lower anchizone metamorphic conditions (slightly below 300 °C). Two types of microfabrics developed primarily by dislocation creep accompanied by dynamic recrystallisation in the absence of twinning. Coaxial kinematics are inferred for samples exhibiting grain shape fabrics and textures with orthorhombic symmetry. Their texture (crystallographic preferred orientation, CPO) is characterised by two c-axis maxima, symmetrically oriented at 15° from the normal to the macroscopic foliation. Non-coaxial deformation is evident in samples with monoclinic shape fabrics and textures characterised by a single oblique c-axis maximum tilted with the sense of shear by about 15°. From the analysis of suitably oriented slip systems for the main texture components under given kinematics it is inferred that the orthorhombic textures, which developed in coaxial kinematics, favour activity of <10–11> and <02–21> slip along the f and r planes, respectively, with minor contributions of basal-<a> slip. In contrast, the monoclinic textures, which developed during simple shear, are most suited for duplex <a> slip along the basal plane. The transition between the dominating slip systems for the orthorhombic and monoclinic microfabrics is interpreted to be due to the effects of dynamic recrystallisation upon texture development. Since oblique c-axis maxima documented in the literature are most often rotated not with but against the shear sense, calcite textures alone should not be used as unequivocal shear sense indicators, but need to be complemented by microstructural criteria such as shape preferred orientations, grain size estimates and amount of twinning.     

Structural geometry and kinematics of the Ailao Shan shear zone: insights from integrated structural,microstructural, and fabric studies of the Yao Shan complex,Yunnan, Southwest China

ABSTRACT

The Yao Shan complex, a massif near the southern segment of the Ailao Shan–Red River (ASRR) shear zone, bears important information on the structural framework of the massif and the kinematics of ductile shearing along the ASRR shear zone. In this contribution, structural, microstructural, quartz c-axis fabric, magnetic fabric, and geochronologic data are used to determine the structural framework of the Yao Shan massif and its tectonic implications for the ASRR shear zone. The Yao Shan complex is characterized by an overall linear A-type antiform that contains a core of high-grade metamorphic rocks with Palaeoproterozoic to Mesozoic protoliths and a mantle of Permo-Triassic low-grade rocks. Both the high-grade metamorphic core and low-grade Permo-Triassic rocks have experienced progressive ductile shearing. Anisotropy of magnetic susceptibility (AMS) results from 17 samples collected along the Xinjie–Pingbian section across the complex show that magnetic lineation (K_max) and foliation (K_max–K_int) are generally subparallel to the corresponding structural elements in the sheared rocks. The shape parameter E values of the magnetic ellipsoids are indicative of dominantly oblate and plane strain, but vary with protolith type and degree of strain among the various rock types. In agreement with the field and microstructural observations, the corrected degree of anisotropy (Pj) values reflect high shear strain in the core rocks and relatively low shear strain in the low-grade strata. A kinematic analysis based on structural and magnetic fabric data shows that both left- and right-lateral shear occurred during the deformation of the Yao Shan complex. Therefore, instead of being an element of the ASRR shear zone, the Yao Shan complex constitutes a crustal-scale inharmonic A-type fold with a fold axis parallel to the stretching lineation. Geochronologic data reveal that the folding occurred coevally with ductile shearing of the middle to lower crust between ca. 30 and 21 Ma.     

A rapid approach to estimate the rockfall energies and distances at the base of rock cliffs

This paper presents a new methodology for estimating the expected energies and first impact distances at the base of a rock cliff, subject to the geometry and properties of the cliff and the representative block being known. The method is based on a sensitivity analysis, conducted by means of kinematic simulations and carried out for a large range of input parameters and their combinations, taking into account the uncertainty associated with their estimate. The proposed approach is validated by comparing predictions to experimental data and shows great potential for a quick qualitative hazard assessment.