Simulating Granite Emplacement and Deformation with Centrifuge Analogue Modelling

The Greater Himalayan leucogranites are a discontinuous suite of intrusions emplaced in a thickened crust during the Miocene southward ductile extrusion of the Himalayan metamorphic core. Melt-induced weakening is thought to have played a critical role in strain localization that facilitated the extrusion. Recent advancements in centrifuge analogue modelling techniques allow for the replication of a broader range of crustal deformation behaviors, enhancing our understanding of large hot orogens. Polydimethylsiloxane (PDMS) is commonly used in centrifuge experiments to model weak melt zones. Difficulties in handling PDMS had, until now, limited its emplacement in models prior to any deformation. A new modelling technique has been developed where PDMS is emplaced into models that have been subjected to some shortening. This technique aims to better understand the effects of melt on strain localization and potential decoupling between structural levels within an evolving orogenic system. Models are subjected to an early stage of shortening, followed by the introduction of PDMS, and then a final stage of shortening. Theoretical percentages of partial melt and their effect on rock strength are considered when adding a specific percentage of PDMS in each model. Due to the limited size of the models, only PDMS sheets of 3 mm thickness were used, which varied in length and width. Within undeformed packages, minimal surface and internal deformation occurred when PDMS is emplaced in the lower layer of the model, showing a vertical volume increase of ~20% within the package; whereas the emplacement of PDMS into the middle layer showed internal dragging of the middle laminations into the lower layer and a vertical volume increase ~30%. Emplacement of PDMS results in ~7% shortening for undeformed and deformed models. Deformed models undergo ~20% additional shortening after two rounds of deformation. Strain localization and decoupling between units occur in deformed models where the degree of deformation changes based on the amount of partial melt present. Surface deformation visible by the formation of a bulge, mode 1 extension cracks and varying surface strain ellipses varies depending if PDMS is present. Better control during emplacement is exhibited when PDMS is added into cooler models, resulting in reduced internal deformation within the middle layer.

Simulating Granite Emplacement and Deformation with Centrifuge Analogue Modelling

The Greater Himalayan leucogranites are a discontinuous suite of intrusions emplaced in a thickened crust during the Miocene southward ductile extrusion of the Himalayan metamorphic core. Melt-induced weakening is thought to have played a critical role in strain localization that facilitated the extrusion. Recent advancements in centrifuge analogue modelling techniques allow for the replication of a broader range of crustal deformation behaviors, enhancing our understanding of large hot orogens. Polydimethylsiloxane (PDMS) is commonly used in centrifuge experiments to model weak melt zones. Difficulties in handling PDMS had, until now, limited its emplacement in models prior to any deformation. A new modelling technique has been developed where PDMS is emplaced into models that have been subjected to some shortening. This technique aims to better understand the effects of melt on strain localization and potential decoupling between structural levels within an evolving orogenic system. Models are subjected to an early stage of shortening, followed by the introduction of PDMS, and then a final stage of shortening. Theoretical percentages of partial melt and their effect on rock strength are considered when adding a specific percentage of PDMS in each model. Due to the limited size of the models, only PDMS sheets of 3 mm thickness were used, which varied in length and width. Within undeformed packages, minimal surface and internal deformation occurred when PDMS is emplaced in the lower layer of the model, showing a vertical volume increase of ~20% within the package; whereas the emplacement of PDMS into the middle layer showed internal dragging of the middle laminations into the lower layer and a vertical volume increase ~30%. Emplacement of PDMS results in ~7% shortening for undeformed and deformed models. Deformed models undergo ~20% additional shortening after two rounds of deformation. Strain localization and decoupling between units occur in deformed models where the degree of deformation changes based on the amount of partial melt present. Surface deformation visible by the formation of a bulge, mode 1 extension cracks and varying surface strain ellipses varies depending if PDMS is present. Better control during emplacement is exhibited when PDMS is added into cooler models, resulting in reduced internal deformation within the middle layer.

Properties of granular analogue model materials: A community wide survey

We report the material properties of 26 granular analogue materials used in 14 analogue modelling laboratories. We determined physical characteristics such as bulk density, grain size distribution, and grain shape, and performed ring shear tests to determine friction angles and cohesion, and uniaxial compression tests to evaluate the compaction behaviour. Mean grain size of the materials varied between c. 100 and 400 μm. Analysis of grain shape factors shows that the four different classes of granular materials (14 quartz sands, 5 dyed quartz sands, 4 heavy mineral sands and 3 size fractions of glass beads) can be broadly divided into two groups consisting of 12 angular and 14 rounded materials. Grain shape has an influence on friction angles, with most angular materials having higher internal friction angles (between c. 35° and 40°) than rounded materials, whereas well-rounded glass beads have the lowest internal friction angles (between c. 25° and 30°). We interpret this as an effect of intergranular sliding versus rolling. Most angular materials have also higher basal friction angles (tested for a specific foil) than more rounded materials, suggesting that angular grains scratch and wear the foil. Most materials have an internal cohesion in the order of 20–100 Pa except for well-rounded glass beads, which show a trend towards a quasi-cohesionless (C < 20 Pa) Coulomb-type material. The uniaxial confined compression tests reveal that rounded grains generally show less compaction than angular grains. We interpret this to be related to the initial packing density after sifting, which is higher for rounded grains than for angular grains. Ring-shear test data show that angular grains undergo a longer strain-hardening phase than more rounded materials. This might explain why analogue models consisting of angular grains accommodate deformation in a more distributed manner prior to strain localisation than models consisting of rounded grains.

Shape fabric development during progressive deformation

Zusammenfassung:Das Ziel dieser Arbeit ist ein besseres Verständnis von der Art und Weise wie sich Formregelungsgefüge entwicklen. Auf dieser Basis wird der Nutzen von Formregelungsgefügen für die Geologie evaluiert. Untersuchungsmethoden sind Geländearbeit und -auswertung, numerische Simulationen und Analogexperimente. Untersuchungen an Formregelungsgefügen in Gesteinen zeigen, daß ein Formregelungsgefüge nur zu einem begrenzten Grad als Anzeiger für die Stärke der Verformung benutzt werden kann. Der angenommene Grund hierfür ist der Einfluß des Verhältnisses von ursprünglicher zu rekristallisierter Korngröße auf die Gefügeentwicklung und von der Art und Weise wie dynamische Rekristallisation ein Gefüge verändert. Um diese Beobachtung zu evaluieren, wurden verschiedene numerische Simulationen von dynamischer Rekristallisation durchgeführt. Ein neuer Deformationsapparat, mit dem generelle Fließregime modelliert werden können, wurde entwickelt. Die rheologischen Eigenschaften von Materialien, die für solche Experimente benutzt werden, wurden untersucht und diskutiert. Ergebnisse von Analogexperimenten zeigen, daß die Intensität eines Formregelungsgefüges positiv mit der Abnahme der 'kinematic vorticity number' und einem nicht-Newtonianischen, 'power law' Verhalten des Materixmaterials korreliert ist. Experimente, in denen die Formveränderung von viskosen Einschlüssen während der progressiven Verformung modelliert werden, zeigen, daß verschiedene Viskositätskontraste zwischen Matrix- und Einschlußmaterial in charakteristische Formgefüge resultieren.

TECHNOLOGICAL ADVANCEMENT IN REMOTE SENSING AND INNOVATIVE APPLICATION TO VOLCANOLOGY

The integration of remote monitoring techniques at different scales is of crucial importance for monitoring of volcanoes and assessment of the associated hazard. In this optic, technological advancement and collaboration between research groups also play a key role. Vhub is a community cyberinfrastructure platform designed for collaboration in volcanology research. Within the Vhub framework, this dissertation focuses on two research themes, both representing novel applications of remotely sensed data in volcanology: advancement in the acquisition of topographic data via active techniques and application of passive multi-spectral satellite data to monitoring of vegetated volcanoes. Measuring surface deformation is a critical issue in analogue modelling of Earth science phenomena. I present a novel application of the Microsoft Kinect sensor to measurement of vertical and horizontal displacements in analogue models. Specifically, I quantified vertical displacement in a scaled analogue model of Nisyros volcano, Greece, simulating magmatic deflation and inflation and related surface deformation, and included the horizontal component to reconstruct 3D models of pit crater formation. The detection of active faults around volcanoes is of importance for seismic and volcanic hazard assessment, but not a simple task to be achieved using analogue models. I present new evidence of neotectonic deformation along a north-south trending fault from the Mt Shasta debris avalanche deposit (DAD), northern California. The fault was identified on an airborne LiDAR campaign of part of the region interested by the DAD and then confirmed in the field. High resolution LiDAR can be utilized also for geomorphological assessment of DADs, and I describe a size-distance analysis to document geomorphological aspects of hummock in the Shasta DAD. Relating the remote observations of volcanic passive degassing to conditions and impacts on the ground provides an increased understanding of volcanic degassing and how satellite-based monitoring can be used to inform hazard management strategies in nearreal time. Combining a variety of satellite-based spectral time series I aim to perform the first space-based assessment of the impacts of sulfur dioxide emissions from Turrialba volcano, Costa Rica, on vegetation in the surrounding environment, and establish whether vegetation indices could be used more broadly to detect volcanic unrest.

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.

Influência da estratigrafia mecânica no desenvolvimento de falhas distensionais em bacias do tipo Rifte: contribuição da modelagem física

In the last decades, analogue modelling has been used in geology to improve the knowledge of how geological structures are nucleated, how they grow and what are the main important points in such processes. The use of this tool in the oil industry, to help seismic interpretations and mainly to search for structural traps contributed to disseminate the use of this tool in the literature. Nowadays, physical modelling has a large field of applications, since landslide to granite emplacement along shear zones. In this work, we deal with physical modelling to study the influence of mechanical stratifications in the nucleation and development of faults and fractures in a context of orthogonal and conjugated oblique basins. To simulate a mechanical stratigraphy we used different materials, with distinct physical proprieties, such as gypsum powder, glass beads, dry clay and quartz sand. Some experiments were run along with a PIV (Particle Image Velocimetry), an instrument that shows the movement of the particles to each deformation moment. Two series of experiments were studied: i) Series MO: We tested the development of normal faults in a context of an orthogonal (to the extension direction) basin. Experiments were run taking into account the change of materials and strata thickness. Some experiments were done with sintectonic sedimentation. We registered differences in the nucleation and growth of faults in layers with different rheological behavior. The gypsum powder layer behaves in a more competent mode, which generates a great number of high angle fractures. These fractures evolve to faults that exhibit a higher dip than when they cross less competent layers, like the one of quartz sand. This competent layer exhibits faulted blocks arranged in a typical domino-style. Cataclastic breccias developed along the faults affecting the competent layers and showed different evolutional history, depending on the deforming stratigraphic sequence; ii) Series MOS2: Normal faults were analyzed in conjugated sub-basins (oblique to the extension direction) developed in a sequence with and without rheological contrast. In experiments with rheological contrast, two important grabens developed along the faulted margins differing from the subbasins with mechanical stratigraphy. Both experiments developed oblique fault systems and, in the area of sub-basins intersection, faults traces became very curved.

Crystal and molecular structure of the thioether-analogue of PEEK from X-ray powder data and diffraction-modelling

Analysis of X-ray powder data for the melt-crystallisable aromatic poly(thioether thioether ketone) [-S-Ar-S-Ar-CO-Ar](n), ('PTTK', Ar= 1,4-phenylene), reveals that it adopts a crystal structure very different from that established for its ether-analogue PEEK. Molecular modelling and diffraction-simulation studies of PTTK show that the structure of this polymer is analogous to that of melt-crystallised poly(thioetherketone) [-SAr-CO-Ar](n) in which the carbonyl linkages in symmetry-related chains are aligned anti-parallel to one another. and that these bridging units are crystallographically interchangeable. The final model for the crystal structure of PTTK is thus disordered, in the monoclinic space group 121a (two chains per unit cell), with cell dimensions a = 7.83, b = 6.06, c = 10.35 angstrom, beta = 93.47 degrees. (c) 2005 Elsevier Ltd. All rights reserved.

Analog Modelling of N-S anticlines of the External Sierras (Jaca basin, Southern Pyrenees)

Report for the scientific sojourn carried out at the Uppsala Universitet, Sweden, from April to July the 2007. Two series of analogue models are used to explore ductile-frictional contrasts of the basal décollement in the development of oblique and transverse structures simultaneously to thin-skinned shortening. These models simulate the evolution of the Central External Sierras (Southern Pyrenees, Spain), which constitute the frontal emerging part of the southernmost Pyrenean thrust sheet. They are characterized by the presence of transverse N-S to NW-SE anticlines, which are perpendicular to the Pyrenean structural trend and developed in the hangingwall of the Santo Domingo thrust system, detaching on an unevenly distributed Triassic materials (evaporitic-dolomitic interfingerings). Model setup performs a décollement made by three patches of silicone neighbouring pure brittle sand. Model series A test the thickness ratio between overburden and décollement. Model series B test the width of frictional detachment areas. Model results show how deformation reaches further in areas detached on ductile layer whereas frictional décollement areas assimilate the strain by means of an additional uplift. This replicates the structural style of Central External Sierras: higher structural relief of N-S anticlines with regard to orogen-parallel structures, absence of a representative ductile décollement in the core, tilting towards the orogen and foreland-side closure not thrusted by the frontal emerging South-Pyrenean thrust.

Modelling ecological niches with support vector machines

1. The ecological niche is a fundamental biological concept. Modelling species' niches is central to numerous ecological applications, including predicting species invasions, identifying reservoirs for disease, nature reserve design and forecasting the effects of anthropogenic and natural climate change on species' ranges. 2. A computational analogue of Hutchinson's ecological niche concept (the multidimensional hyperspace of species' environmental requirements) is the support of the distribution of environments in which the species persist. Recently developed machine-learning algorithms can estimate the support of such high-dimensional distributions. We show how support vector machines can be used to map ecological niches using only observations of species presence to train distribution models for 106 species of woody plants and trees in a montane environment using up to nine environmental covariates. 3. We compared the accuracy of three methods that differ in their approaches to reducing model complexity. We tested models with independent observations of both species presence and species absence. We found that the simplest procedure, which uses all available variables and no pre-processing to reduce correlation, was best overall. Ecological niche models based on support vector machines are theoretically superior to models that rely on simulating pseudo-absence data and are comparable in empirical tests. 4. Synthesis and applications. Accurate species distribution models are crucial for effective environmental planning, management and conservation, and for unravelling the role of the environment in human health and welfare. Models based on distribution estimation rather than classification overcome theoretical and practical obstacles that pervade species distribution modelling. In particular, ecological niche models based on machine-learning algorithms for estimating the support of a statistical distribution provide a promising new approach to identifying species' potential distributions and to project changes in these distributions as a result of climate change, land use and landscape alteration.

Diclofenac plasma protein binding: PK-PD modelling in cardiac patients submitted to cardiopulmonary bypass

Twenty-four surgical patients of both sexes without cardiac, hepatic, renal or endocrine dysfunctions were divided into two groups: 10 cardiac surgical patients submitted to myocardial revascularization and cardiopulmonary bypass (CPB), 3 females and 7 males aged 65 ± 11 years, 74 ± 16 kg body weight, 166 ± 9 cm height and 1.80 ± 0.21 m2 body surface area (BSA), and control, 14 surgical patients not submitted to CPB, 11 female and 3 males aged 41 ± 14 years, 66 ± 14 kg body weight, 159 ± 9 cm height and 1.65 ± 0.16 m2 BSA (mean ± SD). Sodium diclofenac (1 mg/kg, im Voltaren 75® twice a day) was administered to patients in the Recovery Unit 48 h after surgery. Venous blood samples were collected during a period of 0-12 h and analgesia was measured by the visual analogue scale (VAS) during the same period. Plasma diclofenac levels were measured by high performance liquid chromatography. A two-compartment open model was applied to obtain the plasma decay curve and to estimate kinetic parameters. Plasma diclofenac protein binding decreased whereas free plasma diclofenac levels were increased five-fold in CPB patients. Data obtained for analgesia reported as the maximum effect (EMAX) were: 25% VAS (CPB) vs 10% VAS (control), P<0.05, median measured by the visual analogue scale where 100% is equivalent to the highest level of pain. To correlate the effect versus plasma diclofenac levels, the EMAX sigmoid model was applied. A prolongation of the mean residence time for maximum effect (MRTEMAX) was observed without any change in lag-time in CPB in spite of the reduced analgesia reported for these patients, during the time-dose interval. In conclusion, the extent of plasma diclofenac protein binding was influenced by CPB with clinically relevant kinetic-dynamic consequences

Formal modelling for digital media distribution

Human beings have always strived to preserve their memories and spread their ideas. In the beginning this was always done through human interpretations, such as telling stories and creating sculptures. Later, technological progress made it possible to create a recording of a phenomenon; first as an analogue recording onto a physical object, and later digitally, as a sequence of bits to be interpreted by a computer. By the end of the 20th century technological advances had made it feasible to distribute media content over a computer network instead of on physical objects, thus enabling the concept of digital media distribution. Many digital media distribution systems already exist, and their continued, and in many cases increasing, usage is an indicator for the high interest in their future enhancements and enriching. By looking at these digital media distribution systems, we have identified three main areas of possible improvement: network structure and coordination, transport of content over the network, and the encoding used for the content. In this thesis, our aim is to show that improvements in performance, efficiency and availability can be done in conjunction with improvements in software quality and reliability through the use of formal methods: mathematical approaches to reasoning about software so that we can prove its correctness, together with the desirable properties. We envision a complete media distribution system based on a distributed architecture, such as peer-to-peer networking, in which different parts of the system have been formally modelled and verified. Starting with the network itself, we show how it can be formally constructed and modularised in the Event-B formalism, such that we can separate the modelling of one node from the modelling of the network itself. We also show how the piece selection algorithm in the BitTorrent peer-to-peer transfer protocol can be adapted for on-demand media streaming, and how this can be modelled in Event-B. Furthermore, we show how modelling one peer in Event-B can give results similar to simulating an entire network of peers. Going further, we introduce a formal specification language for content transfer algorithms, and show that having such a language can make these algorithms easier to understand. We also show how generating Event-B code from this language can result in less complexity compared to creating the models from written specifications. We also consider the decoding part of a media distribution system by showing how video decoding can be done in parallel. This is based on formally defined dependencies between frames and blocks in a video sequence; we have shown that also this step can be performed in a way that is mathematically proven correct. Our modelling and proving in this thesis is, in its majority, tool-based. This provides a demonstration of the advance of formal methods as well as their increased reliability, and thus, advocates for their more wide-spread usage in the future.