Foliation boudinage

In this thesis foliation boudinage and related structures have been studied based on field observations and numerical modeling. Foliation boudinage occurs in foliated rocks independent of lithology contrast. The developing structures are called ‘Foliation boudinage structures (FBSs)’ and show evidence for both ductile and brittle deformation. They are recognized in rocks by perturbations in monotonous foliation adjacent to a central discontinuity, mostly filled with vein material. Foliation boudinage structures have been studied in the Çine Massif in SW-Turkey and the Furka Pass-Urseren Zone in central Switzerland. Four common types have been distinguished in the field, named after vein geometries in their boudin necks in sections normal to the boudin axis: lozenge-, crescent-, X- and double crescent- type FBSs. Lozengetype FBSs are symmetric and characterized by lozenge-shaped veins in their boudin neck with two cusps facing opposite sides. A symmetrical pair of flanking folds occurs on the two sides of the vein. Crescent-type FBSs are asymmetric with a single smoothly curved vein in the boudin neck, with vein contacts facing to one side. X- and double crescent- type FBSs are asymmetric. The geometry of the neck veins resembles that of cuspate-lobate structures. The geometry of flanking structures is related to the shape of the veins. The veins are mostly filled with massive quartz in large single crystals, commonly associated with tourmaline, feldspar and biotite and in some cases with chlorite. The dominance of large facetted single quartz crystals and spherulitic chlorite in the veins suggest that the minerals grew into open fluidfilled space. FLAC experiments show that fracture propagation during ductile deformation strongly influences the geometry of developing veins. The cusps of the veins are better developed in the case of propagating fractures. The shape of the boudin neck veins in foliation boudinage depends on the initial orientation and shape of the fracture, the propagation behaviour of the fracture, the geometry of bulk flow, and the stage at which mineral filling takes place. A two dimensional discrete element model was used to study the progressive development of foliation boudinage structures and the behavior of visco-elastic material deformed under pure shear conditions. Discrete elements are defined by particles that are connected by visco-elastic springs. Springs can break. A number of simulations was Abstract vii performed to investigate the effect of material properties (Young’s modulus, viscosity and breaking strength) and anisotropy on the developing structures. The models show the development of boudinage in single layers, multilayers and in anisotropic materials with random mica distribution. During progressive deformation different types of fractures develop from mode I, mode II to the combination of both. Voids develop along extension fractures, at intersections of conjugate shear fractures and in small pull-apart structures along shear fractures. These patterns look similar to the natural examples. Fractures are more localized in the models where the elastic constants are low and the competence contrast is high between the layers. They propagate through layers where the constants are high and the competence contrast is relatively low. Flow localize around these fractures and voids. The patterns similar to symmetric boudinage structures and extensional neck veins (e.g. lozenge type) more commonly develop in the models with lower elastic constants and anisotropy. The patterns similar to asymmetric foliation boudinage structures (e.g. X-type) develop associated with shear fractures in the models where elastic constants and anisotropy of the materials are relatively high. In these models boudin neck veins form commonly at pull-aparts along the shear fractures and at the intersection of fractures.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Structural studies and gabbro mylonitization within the Barton Bay Deformation Zone, Geraldton, Ontario /

Structures related to ductile siMple shear parallel to the Bankf ield-Tonbill Fault, define a 5km wide zone, the Barton Bay Deformation Zone. Structures present within this zone Include; simple shear fabrics S, C and C , asymmetric Z shaped folds with rotated axes, boudinage and pinch and swell structures and a subhorlzontal extension llneation. The most highly deformed rock is a gabbro mylonite which occurs in the fault zone. The deformation of this gabbro has been traced in stages from a protomylonite to an ultramylonite In which feldspar and chlorite grainslze has been reduced from over 100 microns to as little as 5 microns. Evidence from the mylonite and the surrounding structure indicates that deformation within the Barton Bay Deformation Zone is related to a regional simple shear zone, the Bankf ield-Tombill Fault. Movement along this shear zone was in a south over north oblique strike slip fashion with a dextral sense of displacement.

Structural, stratigraphic and geochemical studies of the Horwood Peninsula - Gander Bay Area, Northeast Newfoundland /

The Horwood Peninsula - Gander Bay area is located at NE Newfoundland in the Botwood Zone (Williams et a1., 1974) or in the Dunnage Zone (Williams, 1979) of the Central Mobile Belt of the Newfoundland Appalachians. The area is underlain by Middle Ordovician to possible Lower Silurian rocks of the Davidsville and Indian Islands Groups, respectively. Three conformable formations named informally : the Mafic Volcanic Formation, the Greywacke and Siltstone Formation and the Black Slate Formation, have been recognized in the Davidsville Group. The Greywacke and the Black Slate Formations pass locally into a Melange Formation. From consideration of regional structure and abundant locally-derived mafic volcanic olisto- 1iths in the melange, it is considered to have originated by gravity sliding rather than thrusting. Four formations have been recognized in the Indian Islands Group. They mainly contain silty slate and phyllite, grey cherty siltstone, green to red micaceous siltstone and limestone horizons. Repetition of lithological units by F1 folding are well-demonstrated in one of formations in this Group. The major structure in this Group on the Horwood Peninsula is interpreted to be a synclinal complex. The lithology of this Group is different from the Botwood Group to the west and is probably Late Ordovician and/or Early Silurian in age. The effects of soft-sediment deformation can be seen from the lower part of the Davidsville Group to the middle part of the Indian Islands Group indicating continuous and/or episodic slumping and sliding activities throughout the whole area. However, no siginificant depOSitional and tectonic break that could be assigned to the Taconian Orogeny has been recognized in this study. Three periods of tectonic deformation were produced by the Acadian Orogeny. Double boudinage in thin dikes indicates a southeast-northwest sub-horizontal compression and main northeast-southwest sub-horizontal extension during the D1 deformation. A penetrative, axial planar slaty cleavage (Sl) and tight to isocJ.ina1 F1 folds are products of this deformation. The D2 and D3 deformations formed S2 and S3 fabrics associated with crenulations and kink bands which are well-shown in the slates and phyllites of the Indian Islands Group. The D2 and D3 deformations are the products of vertical and northeast-southwest horizontal shortening respectively. The inferred fault between the Ordovician slates (Davidsville Group) and the siltstones (Indian Islands Group) suggested by Williams (1963, 1964b, 1972, 1978) is absent. Formations can be followed without displacement across this inferred fault. Chemically, the pillow lavas, mafic agglomerates, tuff beds and diabase dikes are subdivided into three rock suites : (a) basaltic komatiite (Beaver Cove Assemblage), (b) tholeiitic basalt (diabase dikes), (c) alkaline basalt (Shoal Bay Assemblage). The high Ti02 , MgO, Ni contents and bimodal characteristic of the basaltic komatiite in the area are comparable to the Svartenhuk Peninsula at Baffin Bay and are interpreted to be the result of an abortive volcano-tectonic rift-zone in a rear-arc basin. Modal and chemical analyses of greywackes and siltstones show the trend of maturity of these rocks increasing from poorly sorted Ordovician greywackes to fairly well-sorted Silurian siltstones. Rock fragments in greywackes indicate source areas consisting of plagiogranite, low grade metamorphic rocks and ultramafic rocks. Rare sedimentary structures in both Groups indicate a southeasterly provenance. Trace element analyses of greywackes also reveal a possible island-arc affinity.

Mesozoic fold structure of the Otago and Alpine Schist Belt and its implications on the tectonic evolution of South Island, New Zealand

Erneute Untersuchungen der mesozoischen Faltenstruktur des Otago Schiefergürtels, Südinsel, Neuseeland, zeigen, dass diese aus zwei aufeinander folgenden, ähnlichen, asymmetrischen, offenen bis mäßig engen Großfaltengenerationen im km- Größenbereich besteht anstatt aus den vorher angenommenen Decken- oder Halbfalten. Hauptproblem der Großfaltenstruktur sind Zonen von durchgreifender Boudinage, die in der Nähe der Großfaltenscharniere entstanden sind. Vorherige Bearbeiter deuteten diese Zonen als 'starke Verformungszonen' oder Überschiebungszonen. Diese Arbeit zeigt, dass in diesen Zonen nur durch die asymmetrische Faltung die unteren liegenden Schenkel der Großfalten boudiniert und somit häufig die ansonsten typischen Faltenstrukturen des liegenden Schenkels einer symmetrischen Faltung überprägt wurden. Ein weiteres Problem dieser mesozoischen Großfaltenstruktur ist die Überprägung einer Faltengeneration auf eine frühere. Weil die Verkürzungsrichtung der überprägenden Faltengeneration nicht subparallel zur älteren Faltenachse ist, sondern einen Winkel von rund 30 Grad einschließt, ist ein Wechsel von orthogonalen zu koaxialen Interferenzmustern der Kleinfalten beobachtbar. Folglich ist die Orientierung der Scheitellinie einer überprägenden und überprägten Kleinfalte nicht unbedingt subparallel zur Orientierung der Faltenachse der Großfalte trotz zylindrischer Faltung. Im letzten Teil dieser Arbeit wird die Überprägung der mesozoischen Großfaltenstruktur durch das känozoisch entstandene, transpressionale Alpine Störungssystem, das einen zweiseitigen Falten- und Überschiebungsgürtel im Otago und im Nordwesten anschließenden Alpinen Schiefergürtel bildet, beschrieben.