Tectonics of the Simplon massif and Lepontine gneiss dome: deformation structures due to collision between the underthrusting European plate and the Adriatic indenter

this study presents a review of published geological data, combined with original observations on the tectonics of the simplon massif and the Lepontine gneiss dome in the Western Alps. New observations concern the geometry of the Oligocene Vanzone back fold, formed under amphibolite facies conditions, and of its root between Domodossola and Locarno, which is cut at an acute angle by the Miocene, epi- to anchizonal, dextral centovalli strike-slip fault. the structures of the simplon massif result from collision over 50 Ma between two plate boundaries with a different geometry: the underthrusted European plate and the Adriatic indenter. Detailed mapping and analysis of a complex structural interference pattern, combined with observations on the metamorphic grade of the superimposed structures and radiometric data, allow a kinematic model to be developed for this zone of oblique continental collision. the following main Alpine tectonic phases and structures may be distinguished: 1. NW-directed nappe emplacement, starting in the Early Eocene (similar to 50 Ma); 2. W, SW and S- verging transverse folds; 3. transpressional movements on the dextral simplon ductile shear zone since similar to 32 Ma; 4. formation of the Bergell - Vanzone backfolds and of the southern steep belt during the Oligocene, emplacement of the mantle derived 31 - 29 Ma Bergell and Biella granodiorites and porphyritic andesites as well as intrusions of 29-25 Ma crustal aplites and pegmatites; 5. formation of the dextral discrete Rhone-Simplon line and the centovalli line during the Miocene, accompanied by the pull-apart development of the Lepontine gneiss dome - Dent blanche (Valpelline) depression. It is suggested that movements of shortening in fan shaped NW, W and sW directions accompanied the more regular NW- to WNW-directed displacement of the Adriatic indenter during continental collision.

Dense deformation field estimation for atlas-based segmentation of pathological MR brain images.

Atlas registration is a recognized paradigm for the automatic segmentation of normal MR brain images. Unfortunately, atlas-based segmentation has been of limited use in presence of large space-occupying lesions. In fact, brain deformations induced by such lesions are added to normal anatomical variability and they may dramatically shift and deform anatomically or functionally important brain structures. In this work, we chose to focus on the problem of inter-subject registration of MR images with large tumors, inducing a significant shift of surrounding anatomical structures. First, a brief survey of the existing methods that have been proposed to deal with this problem is presented. This introduces the discussion about the requirements and desirable properties that we consider necessary to be fulfilled by a registration method in this context: To have a dense and smooth deformation field and a model of lesion growth, to model different deformability for some structures, to introduce more prior knowledge, and to use voxel-based features with a similarity measure robust to intensity differences. In a second part of this work, we propose a new approach that overcomes some of the main limitations of the existing techniques while complying with most of the desired requirements above. Our algorithm combines the mathematical framework for computing a variational flow proposed by Hermosillo et al. [G. Hermosillo, C. Chefd'Hotel, O. Faugeras, A variational approach to multi-modal image matching, Tech. Rep., INRIA (February 2001).] with the radial lesion growth pattern presented by Bach et al. [M. Bach Cuadra, C. Pollo, A. Bardera, O. Cuisenaire, J.-G. Villemure, J.-Ph. Thiran, Atlas-based segmentation of pathological MR brain images using a model of lesion growth, IEEE Trans. Med. Imag. 23 (10) (2004) 1301-1314.]. Results on patients with a meningioma are visually assessed and compared to those obtained with the most similar method from the state-of-the-art.

Deformation properties of the Barcelona-Catania-Paris (BCP) energy density functional

We explore the deformation properties of the newly postulated Barcelona-Catania-Paris (BCP) energy density functional (EDF). The results obtained for three isotope chains of Mg, Dy, and Ra are compared to the available experimental data as well as to the results of the Gogny-D1S force. Results for the fission barrier of 240Pu are also discussed.

Correction of through-plane deformation artifacts in stimulated echo acquisition mode cardiac imaging.

Attempts to use a stimulated echo acquisition mode (STEAM) in cardiac imaging are impeded by imaging artifacts that result in signal attenuation and nulling of the cardiac tissue. In this work, we present a method to reduce this artifact by acquiring two sets of stimulated echo images with two different demodulations. The resulting two images are combined to recover the signal loss and weighted to compensate for possible deformation-dependent intensity variation. Numerical simulations were used to validate the theory. Also, the proposed correction method was applied to in vivo imaging of normal volunteers (n = 6) and animal models with induced infarction (n = 3). The results show the ability of the method to recover the lost myocardial signal and generate artifact-free black-blood cardiac images.

Hydrodynamic induced deformation of and orientation of a microscopic elastic filament

We describe simulations of an elastic filament immersed in a fluid and subjected to a body force. The coupling between the fluid flow and the friction that the filament experiences induces bending and alignment perpendicular to the force. With increasing force there are four shape regimes, ranging from slight distortion to an unsteady tumbling motion. We also find marginally stable structures. The instability of these shapes and the alignment are explained by induced bending and nonlocal hydrodynamic interactions. These effects are experimentally relevant for stiff microfilaments.

Nanomechanical properties of molecular organic thin films

Using atomic force microscopy we have studied the nanomechanical response to nanoindentations of surfaces of highly oriented molecular organic thin films (thickness¿1000¿nm). The Young¿s modulus E can be estimated from the elastic deformation using Hertzian mechanics. For the quasi-one-dimensional metal tetrathiafulvalene tetracyanoquinodimethane E~20¿GPa and for the ¿ phase of the p-nitrophenyl nitronyl nitroxide radical E~2GPa. Above a few GPa, the surfaces deform plastically as evidenced by discrete discontinuities in the indentation curves associated to molecular layers being expelled by the penetrating tip.

Detection of millimetric deformation using a terrestrial laser scanner: experiment and application to a rockfall event

Terrestrial laser scanning (TLS) is one of the most promising surveying techniques for rockslope characterization and monitoring. Landslide and rockfall movements can be detected by means of comparison of sequential scans. One of the most pressing challenges of natural hazards is combined temporal and spatial prediction of rockfall. An outdoor experiment was performed to ascertain whether the TLS instrumental error is small enough to enable detection of precursory displacements of millimetric magnitude. This consists of a known displacement of three objects relative to a stable surface. Results show that millimetric changes cannot be detected by the analysis of the unprocessed datasets. Displacement measurement are improved considerably by applying Nearest Neighbour (NN) averaging, which reduces the error (1¿) up to a factor of 6. This technique was applied to displacements prior to the April 2007 rockfall event at Castellfollit de la Roca, Spain. The maximum precursory displacement measured was 45 mm, approximately 2.5 times the standard deviation of the model comparison, hampering the distinction between actual displacement and instrumental error using conventional methodologies. Encouragingly, the precursory displacement was clearly detected by applying the NN averaging method. These results show that millimetric displacements prior to failure can be detected using TLS.

Beneficial Effects of Selected Additives on Asphalt Cement Mixes, HR-278, 1987

Effects of polyolefins, neoprene, styrene-butadiene-styrene (SBS) block copolymers, styrene-butadiene rubber (SBR) latex, and hydrated lime on two asphalt cements were evaluated. Physical and chemical tests were performed on a total of 16 binder blends. Asphalt concrete mixes were prepared and tested with these modified binders and two aggregates (crushed limestone and gravel), each at three asphalt content levels. Properties evaluated on the modified binders (original and thin-film oven aged) included: viscosity at 25 deg C, 60 deg C and 135 deg C with capillary tube and cone-plate viscometer, penetration at 5 deg C and 25 deg C, softening point, force ductility, and elastic recovery at 10 deg C, dropping ball test, tensile strength, and toughness and tenacity tests at 25 deg C. From these the penetration index, the viscosity-temperature susceptibility, the penetration-viscosity number, the critical low-temperature, long loading-time stiffness, and the cracking temperature were calculated. In addition, the binders were studied with x-ray diffraction, reflected fluorescence microscopy, and high-performance liquid chromatography techniques. Engineering properties evaluated on the 72 asphalt concrete mixes containing additives included: Marshall stability and flow, Marshall stiffness, voids properties, resilient modulus, indirect tensile strength, permanent deformation (creep), and effects of moisture by vacuum-saturation and Lottman treatments. Pavement sections of varied asphalt concrete thicknesses and containing different additives were compared to control mixes in terms of structural responses and pavement lives for different subgrades. Although all of the additives tested improved at least one aspect of the binder/mixture properties, no additive was found to improve all the relevant binder/mixture properties at the same time. On the basis of overall considerations, the optimum beneficial effects can be expected when the additives are used in conjunction with softer grade asphalts.

Comparison of Creep and Resilient Modulus Laboratory Results with Field Cores, HR-311, Part 3, 1993

This is Part 3 of a study of creep and resilient modulus testing of hot mix asphalt concrete. The creep and resilient modulus testing in Part 1 showed the improved load carrying characteristics of crushed particles. Cores from pavements drilled in Part 2 exhibited a poor correlation with rutting and creep/resilient modulus on pavement with a range of rut depths. The objective of Part 3 was to determine the relationship of creep and resilient modulus for 1) Marshall specimens from laboratory mixing for mix design; 2) Marshall specimens from construction plant mixing; and 3) cores drilled from the hot mixed asphalt pavement. The creep and resilient modulus data from these three sources exhibited substantial variations. No meaningful correlations of the results from these three sources were obtained.

Development of Asphalt Dynamic Modulus Master Curve Using Falling Weight Deflectometer Measurements, TR-659

The asphalt concrete (AC) dynamic modulus (|E*|) is a key design parameter in mechanistic-based pavement design methodologies such as the American Association of State Highway and Transportation Officials (AASHTO) MEPDG/Pavement-ME Design. The objective of this feasibility study was to develop frameworks for predicting the AC |E*| master curve from falling weight deflectometer (FWD) deflection-time history data collected by the Iowa Department of Transportation (Iowa DOT). A neural networks (NN) methodology was developed based on a synthetically generated viscoelastic forward solutions database to predict AC relaxation modulus (E(t)) master curve coefficients from FWD deflection-time history data. According to the theory of viscoelasticity, if AC relaxation modulus, E(t), is known, |E*| can be calculated (and vice versa) through numerical inter-conversion procedures. Several case studies focusing on full-depth AC pavements were conducted to isolate potential backcalculation issues that are only related to the modulus master curve of the AC layer. For the proof-of-concept demonstration, a comprehensive full-depth AC analysis was carried out through 10,000 batch simulations using a viscoelastic forward analysis program. Anomalies were detected in the comprehensive raw synthetic database and were eliminated through imposition of certain constraints involving the sigmoid master curve coefficients. The surrogate forward modeling results showed that NNs are able to predict deflection-time histories from E(t) master curve coefficients and other layer properties very well. The NN inverse modeling results demonstrated the potential of NNs to backcalculate the E(t) master curve coefficients from single-drop FWD deflection-time history data, although the current prediction accuracies are not sufficient to recommend these models for practical implementation. Considering the complex nature of the problem investigated with many uncertainties involved, including the possible presence of dynamics during FWD testing (related to the presence and depth of stiff layer, inertial and wave propagation effects, etc.), the limitations of current FWD technology (integration errors, truncation issues, etc.), and the need for a rapid and simplified approach for routine implementation, future research recommendations have been provided making a strong case for an expanded research study.

Palaeozoic deformation and magmatism in the northern area of the Anatolide block (Konya), witness of the Palaeotethys active margin

The study area. located north of Konva (Central Turkey), is composed of Silurian to Cretaceous metamorphosed rocks. The lower unit of the oldest formation (Silurian-Early Permian) is mostly made up of Silurian-Early Carboniferous metacarbonates. These rocks pass laterally and vertically to Devonian-Early Permian series having continental margin, shallow water and pelagic characteristics. They are intruded or juxtaposed to different kinds of metamagmatic rocks. which show MORB. continental arc and within plate characteristics. The Palaeozoic units are covered unconformably by Triassic-Cretaceous metasedimentary units. All these rocks are overthrusted by Mesozoic ophiolites. The Palaeozoic sequence can be seen as a northern Palaeotethys passive, then active margin. The northward subduction of the Palaeotethys ocean during the Carboniferous-Triassic times, induced the development of a magmatic arc and fore-arc sequence (Carboniferous-Permian). Before the Early Triassic (?Late Permian) time. the fore-arc sequence was uplifted above sea level and eroded. The Triassic sequences are regarded as marking the onset of back-arc opening and detachment of the Anatolian Konya block from the active Eurasian margin. Finally. a suture zone formed during the Carman between the Konya region and the Menderes-Tauride Cimmerian block due to the closing of Palaeotethvs. This geodynamic evolution can be correlated with the evolution of the Karaburun sequence in western Turkey.

Three-dimensional hydrodynamic numerical modelling of fold nappe formation, basement-cover deformation and slab detachment with applications to the helvetic nappe system (W Switzerland)

Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth's crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These struc-tures occur from the centimeter to the kilometer scale and the related deformation processes con-trol the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D defor-mation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as fold and cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. Especially, we are interested in fold nappes which are recumbent folds with amplitudes usually exceeding 10 km and they have been presumably formed by ductile shearing. They often exhibit a constant sense of shearing and a non-linear increase of shear strain towards their overturned limb. The fold axes of the Morcles fold nappe in western Switzerland plunges to the ENE whereas the fold axes in the more eastern Doldenhorn nappe plunges to the WSW. These opposite plunge direc-tions characterize the Rawil depression (Wildstrubel depression). The Morcles nappe is mainly the result of layer parallel contraction and shearing. During the compression the massive lime-stones were more competent than the surrounding marls and shales, which led to the buckling characteristics of the Morcles nappe, especially in the north-dipping normal limb. The Dolden-horn nappe exhibits only a minor overturned fold limb. There are still no 3-D numerical studies which investigate the fundamental dynamics of the formation of the large-scale 3-D structure including the Morcles and Doldenhorn nappes and the related Rawil depression. We study the 3-D evolution of geometrical instabilities and fold nappe formation with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FE algorithm combines a nu-merical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FE mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D de-formation for a visco-elastic rheology. The viscous rheology is described by a power-law flow law. The code is used to study the 3-D fold nappe formation, the lateral propagation of folding and also the lateral propagation of cusps due to initial half graben geometry. Thereby, the small initial geometrical perturbations for folding and necking are exactly followed by the FE mesh, whereas the initial large perturbation describing a half graben is defined by a contour line inter-secting the finite elements. Further, the 3-D algorithm is applied to 3-D viscous nacking during slab detachment. The results from various simulations are compared with 2-D resulats and a 1-D analytical solution. -- On retrouve beaucoup de structures en 3 dimensions (3-D) dans les roches qui ont pour origines une déformation de la lithosphère terrestre. Ces structures sont par exemple des plis, des boudins (pinch-and-swell) ou des mullions (cuspate-lobate) et sont présentés de l'échelle centimétrique à kilométrique. Mécaniquement, ces structures peuvent être expliquées par une différence de résistance entre les différentes unités de roches et sont généralement le fruit d'une instabilité géométrique. Ces différences mécaniques entre les unités contrôlent non seulement les types de structures rencontrées, mais également le type de déformation (thick skin, thin skin) et le style tectonique (bassin d'avant pays, chaîne d'avant pays). Les processus de la déformation en deux dimensions (2-D) formant ces structures sont relativement bien compris. Cependant, lorsque l'on ajoute la troisiéme dimension, plusieurs processus ne sont pas complètement compris lors de la déformation à large échelle. L'un de ces processus est la propagation latérale des structures, par exemple la propagation de plis ou de mullions dans la direction perpendiculaire à l'axe de com-pression, ou la propagation des zones d'amincissement des boudins perpendiculairement à la direction d'extension. Nous sommes particulièrement intéressés les nappes de plis qui sont des nappes de charriage en forme de plis couché d'une amplitude plurikilométrique et étant formées par cisaillement ductile. La plupart du temps, elles exposent un sens de cisaillement constant et une augmentation non linéaire de la déformation vers la base du flanc inverse. Un exemple connu de nappes de plis est le domaine Helvétique dans les Alpes de l'ouest. Une de ces nap-pes est la Nappe de Morcles dont l'axe de pli plonge E-NE tandis que de l'autre côté de la dépression du Rawil (ou dépression du Wildstrubel), la nappe du Doldenhorn (équivalent de la nappe de Morcles) possède un axe de pli plongeant O-SO. La forme particulière de ces nappes est due à l'alternance de couches calcaires mécaniquement résistantes et de couches mécanique-ment faibles constituées de schistes et de marnes. Ces différences mécaniques dans les couches permettent d'expliquer les plissements internes à la nappe, particulièrement dans le flanc inver-se de la nappe de Morcles. Il faut également noter que le développement du flanc inverse des nappes n'est pas le même des deux côtés de la dépression de Rawil. Ainsi la nappe de Morcles possède un important flanc inverse alors que la nappe du Doldenhorn en est presque dépour-vue. A l'heure actuelle, aucune étude numérique en 3-D n'a été menée afin de comprendre la dynamique fondamentale de la formation des nappes de Morcles et du Doldenhorn ainsi que la formation de la dépression de Rawil. Ce travail propose la première analyse de l'évolution 3-D des instabilités géométriques et de la formation des nappes de plis en utilisant des simulations numériques. Notre modèle est basé sur la méthode des éléments finis (FEM) qui permet de ré-soudre avec précision les interfaces entre deux matériaux ayant des propriétés mécaniques très différentes (par exemple entre les couches calcaires et les couches marneuses). De plus nous utilisons un maillage lagrangien déformable avec une fonction de re-meshing (production d'un nouveau maillage). Grâce à cette méthode combinée il nous est possible de suivre avec précisi-on les interfaces matérielles et de résoudre avec précision les instabilités géométriques lors de la déformation de matériaux visco-élastiques décrit par une rhéologie non linéaire (n>1). Nous uti-lisons cet algorithme afin de comprendre la formation des nappes de plis, la propagation latérale du plissement ainsi que la propagation latérale des structures de type mullions causé par une va-riation latérale de la géométrie (p.ex graben). De plus l'algorithme est utilisé pour comprendre la dynamique 3-D de l'amincissement visqueux et de la rupture de la plaque descendante en zone de subduction. Les résultats obtenus sont comparés à des modèles 2-D et à la solution analytique 1-D. -- Viele drei dimensionale (3-D) Strukturen, die in Gesteinen vorkommen und durch die Verfor-mung der Erdkruste und Litosphäre entstanden sind werden von den unterschiedlichen mechani-schen Eigenschaften der Gesteinseinheiten kontrolliert und sind häufig das Resulat von geome-trischen Istabilitäten. Zu diesen strukturen zählen zum Beispiel Falten, Pich-and-swell Struktu-ren oder sogenannte Cusbate-Lobate Strukturen (auch Mullions). Diese Strukturen kommen in verschiedenen Grössenordungen vor und können Masse von einigen Zentimeter bis zu einigen Kilometer aufweisen. Die mit der Entstehung dieser Strukturen verbundenen Prozesse kontrol-lieren die Entstehung von Gerbirgen und Sediment-Becken sowie die Verformung des Kontaktes zwischen Grundgebirge und Stedimenten. Die zwei dimensionalen (2-D) Verformungs-Prozesse die zu den genannten Strukturen führen sind bereits sehr gut untersucht. Einige Prozesse wäh-rend starker 3-D Verformung sind hingegen noch unvollständig verstanden. Einer dieser 3-D Prozesse ist die seitliche Fortpflanzung der beschriebenen Strukturen, so wie die seitliche Fort-pflanzung von Falten und Cusbate-Lobate Strukturen senkrecht zur Verkürzungsrichtung und die seitliche Fortpflanzung von Pinch-and-Swell Strukturen othogonal zur Streckungsrichtung. Insbesondere interessieren wir uns für Faltendecken, liegende Falten mit Amplituden von mehr als 10 km. Faltendecken entstehen vermutlich durch duktile Verscherung. Sie zeigen oft einen konstanten Scherungssinn und eine nicht-lineare zunahme der Scherverformung am überkipp-ten Schenkel. Die Faltenachsen der Morcles Decke in der Westschweiz fallen Richtung ONO während die Faltenachsen der östicher gelegenen Doldenhorn Decke gegen WSW einfallen. Diese entgegengesetzten Einfallrichtungen charakterisieren die Rawil Depression (Wildstrubel Depression). Die Morcles Decke ist überwiegend das Resultat von Verkürzung und Scherung parallel zu den Sedimentlagen. Während der Verkürzung verhielt sich der massive Kalkstein kompetenter als der Umliegende Mergel und Schiefer, was zur Verfaltetung Morcles Decke führ-te, vorallem in gegen Norden eifallenden überkippten Schenkel. Die Doldenhorn Decke weist dagegen einen viel kleineren überkippten Schenkel und eine stärkere Lokalisierung der Verfor-mung auf. Bis heute gibt es keine 3-D numerischen Studien, die die fundamentale Dynamik der Entstehung von grossen stark verformten 3-D Strukturen wie den Morcles und Doldenhorn Decken sowie der damit verbudenen Rawil Depression untersuchen. Wir betrachten die 3-D Ent-wicklung von geometrischen Instabilitäten sowie die Entstehung fon Faltendecken mit Hilfe von numerischen Simulationen basiert auf der Finite Elemente Methode (FEM). Die Simulation von geometrischen Instabilitäten, die aufgrund von Änderungen der Materialeigenschaften zwischen verschiedenen Gesteinseinheiten entstehen, erfortert einen numerischen Algorithmus, der in der Lage ist die Materialgrenzen mit starkem Kontrast der Materialeigenschaften (zum Beispiel zwi-schen Kalksteineinheiten und Mergel) für starke Verfomung genau aufzulösen. Um dem gerecht zu werden kombiniert unser FE Algorithmus eine numerische Contour-Linien-Technik und ein deformierbares Lagranges Netz mit Re-meshing. Mit dieser kombinierten Methode ist es mög-lich den anfänglichen Materialgrenzen mit dem FE Netz genau zu folgen und die geometrischen Instabilitäten genügend aufzulösen. Der Algorithmus ist in der Lage visko-elastische 3-D Ver-formung zu rechnen, wobei die viskose Rheologie mit Hilfe eines power-law Fliessgesetzes beschrieben wird. Mit dem numerischen Algorithmus untersuchen wir die Entstehung von 3-D Faltendecken, die seitliche Fortpflanzung der Faltung sowie der Cusbate-Lobate Strukturen die sich durch die Verkürzung eines mit Sediment gefüllten Halbgraben bilden. Dabei werden die anfänglichen geometrischen Instabilitäten der Faltung exakt mit dem FE Netz aufgelöst wäh-rend die Materialgranzen des Halbgrabens die Finiten Elemente durchschneidet. Desweiteren wird der 3-D Algorithmus auf die Einschnürung während der 3-D viskosen Plattenablösung und Subduktion angewandt. Die 3-D Resultate werden mit 2-D Ergebnissen und einer 1-D analyti-schen Lösung verglichen.