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.

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.

The influence of nano-scale second-phase particles on deformation of fine grained calcite mylonites

Grey and white carbonate mylonites were collected along thrust planes of the Helvetic Alps. They are characterised by very small grain sizes and non-random grain shape (SPO) and crystallographic preferred orientation (CPO). Presumably they deformed in the field of grain size sensitive flow by recrystallisation accommodated intracrystalline deformation in combination with granular flow. Both mylonites show a similar mean grain size, but in the grey mylonites the grain size range is larger, the grain shapes are more elongate and the dynamically recrystallised calcite grains are more often twinned. Grey mylonites have an oblique CPO, while the CPO in white mylonites is symmetric with respect to the shear plane. Combustion analysis and TEM investigations revealed that grey mylonites contain a higher amount of highly structured kerogens with particle sizes of a few tens of nanometers, which are finely dispersed at the grain boundaries. During deformation of the rock, nano-scale particles reduced the migration velocity of grain boundaries by Zener drag resulting in slower recrystallisation rates of the calcite aggregate. In the grey mylonites, more strain increments were accommodated by individual grains before they became refreshed by dynamic recrystallisation than in white mylonites, where grain boundary migration was less hindered and recrystallisation cycles were faster. Consequently, grey mylonites represent ‘deformation’ microfabrics while white mylonites are characterised by ‘recrystallisation’ microfabrics. Field geologists must utilise this different deformation behavior when applying the obliquity in CPO and SPO of the respective mylonites as reliable shear sense indicators.

Modelling the mechanical aspects of the curing process of magneto-sensitive elastomeric materials

In this paper, a phenomenologically motivated magneto-mechanically coupled finite strain elastic framework for simulating the curing process of polymers in the presence of a magnetic load is proposed. This approach is in line with previous works by Hossain and co-workers on finite strain curing modelling framework for the purely mechanical polymer curing (Hossain et al., 2009b). The proposed thermodynamically consistent approach is independent of any particular free energy function that may be used for the fully-cured magneto-sensitive polymer modelling, i.e. any phenomenological or micromechanical-inspired free energy can be inserted into the main modelling framework. For the fabrication of magneto-sensitive polymers, micron-size ferromagnetic particles are mixed with the liquid matrix material in the uncured stage. The particles align in a preferred direction with the application of a magnetic field during the curing process. The polymer curing process is a complex (visco) elastic process that transforms a fluid to a solid with time. Such transformation process is modelled by an appropriate constitutive relation which takes into account the temporal evolution of the material parameters appearing in a particular energy function. For demonstration in this work, a frequently used energy function is chosen, i.e. the classical Mooney-Rivlin free energy enhanced by coupling terms. Several representative numerical examples are demonstrated that prove the capability of our approach to correctly capture common features in polymers undergoing curing processes in the presence of a magneto-mechanical coupled load.

Quantitative analysis of crystal/grain sizes and their distributions in 2D and 3D

We review methods to estimate the average crystal (grain) size and the crystal (grain) size distribution in solid rocks. Average grain sizes often provide the base for stress estimates or rheological calculations requiring the quantification of grain sizes in a rock's microstructure. The primary data for grain size data are either 1D (i.e. line intercept methods), 2D (area analysis) or 3D (e.g., computed tomography, serial sectioning). These data have been used for different data treatments over the years, whereas several studies assume a certain probability function (e.g., logarithm, square root) to calculate statistical parameters as the mean, median, mode or the skewness of a crystal size distribution. The finally calculated average grain sizes have to be compatible between the different grain size estimation approaches in order to be properly applied, for example, in paleo-piezometers or grain size sensitive flow laws. Such compatibility is tested for different data treatments using one- and two-dimensional measurements. We propose an empirical conversion matrix for different datasets. These conversion factors provide the option to make different datasets compatible with each other, although the primary calculations were obtained in different ways. In order to present an average grain size, we propose to use the area-weighted and volume-weighted mean in the case of unimodal grain size distributions, respectively, for 2D and 3D measurements. The shape of the crystal size distribution is important for studies of nucleation and growth of minerals. The shape of the crystal size distribution of garnet populations is compared between different 2D and 3D measurements, which are serial sectioning and computed tomography. The comparison of different direct measured 3D data; stereological data and direct presented 20 data show the problems of the quality of the smallest grain sizes and the overestimation of small grain sizes in stereological tools, depending on the type of CSD. (C) 2011 Published by Elsevier Ltd.

Modelling the curing process in magneto-sensitive polymers: Rate-dependence and shrinkage

This paper deals with a phenomenologically motivated magneto-viscoelastic coupled finite strain framework for simulating the curing process of polymers under the application of a coupled magneto-mechanical road. Magneto-sensitive polymers are prepared by mixing micron-sized ferromagnetic particles in uncured polymers. Application of a magnetic field during the curing process causes the particles to align and form chain-like structures lending an overall anisotropy to the material. The polymer curing is a viscoelastic complex process where a transformation from fluid. to solid occurs in the course of time. During curing, volume shrinkage also occurs due to the packing of polymer chains by chemical reactions. Such reactions impart a continuous change of magneto-mechanical properties that can be modelled by an appropriate constitutive relation where the temporal evolution of material parameters is considered. To model the shrinkage during curing, a magnetic-induction-dependent approach is proposed which is based on a multiplicative decomposition of the deformation gradient into a mechanical and a magnetic-induction-dependent volume shrinkage part. The proposed model obeys the relevant laws of thermodynamics. Numerical examples, based on a generalised Mooney-Rivlin energy function, are presented to demonstrate the model capacity in the case of a magneto-viscoelastically coupled load.

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.

Condicionamento estrutural das zonasa de cisalhamento da região de Forquilha, domínio Ceará Central: uma abordagem integrada de sensoriamento remoto e geologia estrutural

The studied region, named Forquilha and localized in northwestern Central Ceará domain (northern portion of Borborema Province), presents a lithostratigraphic framework constituted by paleoproterozoic metaplutonics, metasedimentary sequences and neoproterozoic granitoids. The metasedimentary rocks of Ceará group occupy most part of the area. This group is subdivided in two distinct units: Canindé and Independência. Canindé unit is represented basically by biotite paragneisses and muscovite paragneisses, with minor metabasic rocks (amphibolite lens). Independência sequence is composed by garnetiferous paragneisses, sillimanite-garnet-quartz-muscovite schists and quartz-muscovite schists, pure or muscovite quartzites and rare marbles. At least three ductile deformation events were recognized in both units of Ceará group, named D1, D2 and D3. The former one is interpreted as related to a low angle tangential tectonics which mass transport is southward. D2 event is marked by the development of close/isoclinal folds with a N-S oriented axis. Refolding patterns generated by F1 and F2 superposition are found in several places. The latest event (D3) corresponds to a transcurrent tectonics, which led to development of mega-folds and several shear zones, under a transpressional regime. The mapped shear zones are Humberto Monte (ZCHM), Poço Cercado (ZCPC) and Forquilha (ZCF). Digital image processing of enhanced Landsat 7-ETM+ satellite images, combined with field data, demonstrate that these penetrative structures are associated with positive and negative geomorphologic patterns, distributed in linear and curvilinear arrangements with tonal banding, corresponding to the ductile fabric and to crests. Diverse color composites were tested and RGB-531 and RGB-752 provided the best results for lineament analysis of the most prominent shear zones. Spatial filtering techniques (3x3 and 5x5 filters) were also used and the application of Prewitt filters generated the best products. The integrated analysis of morphological and textural aspects from filtered images, variation of tonalities related to the distribution of geologic units in color composites and the superposition over a digital elevation model, contributed to a characterization of the structural framework of the study area. Kinematic compatibility of ZCHM, ZCPC, ZCF shear zones, as well as Sobral-Pedro II (ZCSPII) shear zone, situated to the west of the study area, was one of the goal of this work. Two of these shear zones (ZCHM, ZCPC) display sinistral movements, while the others (ZCSPII, ZCF) exhibit dextral kinematics. 40Ar/39Ar ages obtained in this thesis for ZCSPII and ZCPC, associated with other 40Ar/39Ar data of adjacent areas, indicate that all these shear zones are related to Brasiliano orogeny. The trend of the structures, the opposite shear senses and the similar metamorphic conditions are fitted in a model based on the development of conjugate shear zones in an unconfined transpression area. A WNW-ESE bulk shortening direction is infered. The geometry and kinematic of the studied structures suggest that shortening was largely accommodated by lateral extrusion, with only minor amounts of vertical stretch

Estrutura, geocronologia e alojamento dos batólitos de Quixadá, Quixeramobim e Senandor Pompeu - Ceará central

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

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.

Structural verification of the ESEO spacecraft: from launch loads analysis to composite panels verification

The present thesis work was performed in the frame of ESEO (European Student Earth Orbiter) project. The activities that are described in this document were carried out in the Microsatellites and Space Micro systems Lab led by Professor Paolo Tortora and in ALMASpace company facilities. The thesis deals with ESEO structural analysis, at system and unit level, and verification: after determining the design limit loads to be applied to the spacecraft as an envelope of different launchers load profiles, a finite element structural analysis was performed on the model of the satellite in order to ensure the capability to withstand the loads encountered during the launch; all the analyses were performed according to ESA standards and using the software MSC NASTRAN SIMXPERT. Amplification factors were derived and used to determine loads to be considered at unit level. In particular structural analyses were carried out on the GPS unit, the payload developed for ESEO by students of University of Bologna and results were used in the preparation of GPS payload design definition file. As for the verification phase a study on the panels and inserts to be used in the spacecraft was performed: different designs were created exploiting methods to optimize weight and mechanical behavior. The configurations have been analyzed and results compared to select the final design.

Uncertainty in Structural Interpretation : Lessons to be learnt

Peer reviewed

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.

Guidebooks to the geology of Maryland.

Each part has special t.p. and separate paging.