Benchmarking numerical models of brittle thrust wedges

We report quantitative results from three brittle thrust wedge experiments, comparing numerical results directly with each other and with corresponding analogue results. We first test whether the participating codes reproduce predictions from analytical critical taper theory. Eleven codes pass the stable wedge test, showing negligible internal deformation and maintaining the initial surface slope upon horizontal translation over a frictional interface. Eight codes participated in the unstable wedge test that examines the evolution of a wedge by thrust formation from a subcritical state to the critical taper geometry. The critical taper is recovered, but the models show two deformation modes characterised by either mainly forward dipping thrusts or a series of thrust pop-ups. We speculate that the two modes are caused by differences in effective basal boundary friction related to different algorithms for modelling boundary friction. The third experiment examines stacking of forward thrusts that are translated upward along a backward thrust. The results of the seven codes that run this experiment show variability in deformation style, number of thrusts, thrust dip angles and surface slope. Overall, our experiments show that numerical models run with different numerical techniques can successfully simulate laboratory brittle thrust wedge models at the cm-scale. In more detail, however, we find that it is challenging to reproduce sandbox-type setups numerically, because of frictional boundary conditions and velocity discontinuities. We recommend that future numerical-analogue comparisons use simple boundary conditions and that the numerical Earth Science community defines a plasticity test to resolve the variability in model shear zones.

Benchmarking analogue models of brittle thrust wedges

We performed a quantitative comparison of brittle thrust wedge experiments to evaluate the variability among analogue models and to appraise the reproducibility and limits of model interpretation. Fifteen analogue modeling laboratories participated in this benchmark initiative. Each laboratory received a shipment of the same type of quartz and corundum sand and all laboratories adhered to a stringent model building protocol and used the same type of foil to cover base and sidewalls of the sandbox. Sieve structure, sifting height, filling rate, and details on off-scraping of excess sand followed prescribed procedures. Our analogue benchmark shows that even for simple plane-strain experiments with prescribed stringent model construction techniques, quantitative model results show variability, most notably for surface slope, thrust spacing and number of forward and backthrusts. One of the sources of the variability in model results is related to slight variations in how sand is deposited in the sandbox. Small changes in sifting height, sifting rate, and scraping will result in slightly heterogeneous material bulk densities, which will affect the mechanical properties of the sand, and will result in lateral and vertical differences in peak and boundary friction angles, as well as cohesion values once the model is constructed. Initial variations in basal friction are inferred to play the most important role in causing model variability. Our comparison shows that the human factor plays a decisive role, and even when one modeler repeats the same experiment, quantitative model results still show variability. Our observations highlight the limits of up-scaling quantitative analogue model results to nature or for making comparisons with numerical models. The frictional behavior of sand is highly sensitive to small variations in material state or experimental set-up, and hence, it will remain difficult to scale quantitative results such as number of thrusts, thrust spacing, and pop-up width from model to nature.

Sandbox experiments on basin inversion: testing the influence of basin orientation and basin fill

Analogue modelling experiments using brittle materials are performed to study the inversion of extensional structures. Asymmetric grabens of two different orientations are first created during a phase of extension and progressively filled. They are subsequently shortened in the same direction. The aim of our experiments is to determine factors affecting the style of deformation during inversion. We specifically investigate variations in thickness and distribution of strong and weak layers constituting the graben fill and in initial basin orientation. The main advantage of our experimental set-up is that we have a complete control on graben location, width, infill and orientation before inversion. The experiments show that shortening results only in limited reactivation of pre-existing normal faults. In general, forward thrusts and backthrusts cut across normal faults into the footwall of the graben. The forward thrusts either propagate parallel to the enveloping surface of faulted blocks or they cut across basin-limiting normal faults at various angles. The graben fill is mechanically extruded by displacement along forward thrusts that accommodate most of the shortening. Both pre-existing faults and weak graben fill act as zones of weakness during inversion and determine the orientation and location of both backthrusts and forward thrusts. The results of our experiments conform well to natural examples of inverted graben structures.

Comparisons between analogue and numerical models of thrust wedge development

Analogue and finite element numerical models with frictional and viscous properties are used to model thrust wedge development. Comparison between model types yields valuable information about analogue model evolution, scaling laws and the relative strengths and limitations of the techniques. Both model types show a marked contrast in structural style between ‘frictional-viscous domains’ underlain by a thin viscous layer and purely ‘frictional domains’. Closely spaced thrusts form a narrow and highly asymmetric fold-and-thrust belt in the frictional domain, characterized by in-sequence propagation of forward thrusts. In contrast, the frictional-viscous domain shows a wide and low taper wedge and a thrust belt with a more symmetrical vergence, with both forward and back thrusts. The frictional-viscous domain numerical models show that the viscous layer initially simple shears as deformation propagates along it, while localized deformation resulting in the formation of a pop-up structure occurs in the overlying frictional layers. In both domains, thrust shear zones in the numerical model are generally steeper than the equivalent faults in the analogue model, because the finite element code uses a non-associated plasticity flow law. Nevertheless, the qualitative agreement between analogue and numerical models is encouraging. It shows that the continuum approximation used in numerical models can be used to model frictional materials, such as sand, provided caution is taken to properly scale the experiments, and some of the limitations are taken into account.

Fracture studies in the Zuni and Lucero uplifts, New Mexico.

"Contract No. AT(30-1)-1390."

On the new Red Sandstone Series in the Basin of the Eden, and north-western coasts of Cumberland and Lancashire /

Reprint. Originally published: London : Geological Society of London. (Geol. trans. 2nd ser. ; v. 4).

The nature and timing of tectonism in the western facies terrane of Nevada and California : an outline of evidence and interpretations derived from geologic maps of key areas /

Mode of access: Internet.

Die mutmasslichen geologischen Profile des neuen Hauenstein- und des Grenchenbergtunnels im Schweizer Jura, von A. Buxtorf.

Separatabdruck aus den Verhandlungen der Naturforschenden Gesellschaft in Basel. Bd. 24.

Seismically-Induced Rock-Slope Failure: Numerical Investigation using the Bonded Particle Model

Thesis (Ph.D.)--University of Washington, 2016-06

Partitioning of deformation within a subduction channel during exhumation of high-pressure rocks: a case study from the Western Alps

The metamorphic belt of the Western Alps was subjected to widespread extensional tectonism at the end of the Eocene (ca. 45-35 Ma). Extension was accommodated by hinterland-directed movements along gently inclined extensional shear zones, which facilitated rapid exhumation of high-pressure and ultra-high-pressure rocks. This deformation resulted in a normal metamorphic sequence. Extension in the inner parts of the Western Alps was coeval with shortening at the front of the belt (foreland-directed thrusts), which took place during decompression, and emplaced higher grade metamorphic units over lower grade metamorphic rocks, thus forming an inverse metamorphic sequence. Two mechanisms for this extensional episode are discussed: (1) collapse of an overthickened lithosphere, and (2) internal readjustments within the orogenic wedge due to subduction channel dynamics. We favour the latter mechanism because it can account for the development of the observed inverse and normal metamorphic sequences along foreland-directed thrusts and hinterland-directed detachments, respectively. This hypothesis is supported by published structural, metamorphic and geochronological data from four geological transects through the Western Alps. This study also emphasizes the importance of post-shearing deformation (e.g. horizontal buckling versus vertical flattening), which can modify the distribution of hinterland- and foreland-directed shear zones in orogenic belts. (c) 2006 Elsevier Ltd. All rights reserved.

Cycles of passive versus active diapirism recorded along an exposed salt wall

Acknowledgements GIA is grateful for funding from the Carnegie Trust for the Universities of Scotland that enabled fieldwork for this project. RW was supported by the Israel Science Foundation (ISF grant No. 1245/11). SM was supported by the Israel Science Foundation (ISF grant No. 1436/14). We would like to thank Chris Talbot and Yohann Poprawski for careful and constructive reviews. The authors appreciate the help of Nicolas Waldmann in precisely locating the positons of dated unconformities.

Ground penetrating radar (GPR) and electrical resistivity tomography (ERT) survey on the neotectonic activity of the Padul-Nigüelas Fault Zone (PNFZ), Southern Spain

The Padul-Nigüelas Fault Zone (PNFZ) is situated at the south-western mountain front of the Sierra Nevada (Spain) in an extensive regime and belongs to the internal zone of the Betic Cordilleras. The aim of this study is a collection of new evidence for neotectonic activity of the fault zone with classical geological field work and modern geophysical methods, such as ground penetrating radar (GPR). Among an apparently existing bed rock fault scarp with triangular facets, other evidences, such as deeply incised valleys and faults in the colluvial wedges, are present in the PNFZ. The preliminary results of our recent field work have shown that the synsedimentary faults within the colluvial sediments seem to propagate basinwards and the bed rock fault is only exhumed due to erosion for the studied segment (west of Marchena). We will use further GPR data and geomorphologic indices to gather further evidences of neotectonic activity of the PNFZ.