Magma injections and destabilization of basaltic volcanoes: a numerical study. Application to La Reunion (Indian ocean, France) and Stromboli (Tyrrhenian sea, Italy)

Most basaltic volcanoes are affected by recurrent lateral instabilities during their evolution. Numerous factors have been shown to be involved in the process of flank destabilization occurring over long periods of time or by instantaneous failures. However, the role of these factors on the mechanical behaviour and stability of volcanic edifices is poorly-constrained as lateral failure usually results from the combined effects of several parameters. Our study focuses on the morphological and structural comparison of two end-member basaltic systems, La Reunion (Indian ocean, France) and Stromboli (southern Tyrrhenian sea, Italy). We showed that despite major differences on their volumes and geodynamic settings, both systems present some similarities as they are characterized by an intense intrusive activity along well-developed rift zones and recurrent phenomena of flank collapse during their evolution. Among the factors of instability, the examples of la Reunion and Stromboli evidence the major contribution of intrusive complexes to volcano growth and destruction as attested by field observations and the monitoring of these active volcanoes. Classical models consider the relationship between vertical intrusions of magma and flank movements along a preexisting sliding surface. A set of published and new field data from Piton des Neiges volcano (La Reunion) allowed us to recognize the role of subhorizontal intrusions in the process of flank instability and to characterize the geometry of both subvertical and subhorizontal intrusions within basaltic edifices. This study compares the results of numerical modelling of the displacements associated with high-angle and low-angle intrusions within basaltic volcanoes. We use a Mixed Boundary Element Method to investigate the mechanical response of an edifice to the injection of magmatic intrusions in different stress fields. Our results indicate that the anisotropy of the stress field favours the slip along the intrusions due to cointrusive shear stress, generating flank-scale displacements of the edifice, especially in the case of subhorizontal intrusions, capable of triggering large-scale flank collapses on basaltic volcanoes. Applications of our theoretical results to real cases of flank displacements on basaltic volcanoes (such as the 2007 eruptive crisis at La Reunion and Stromboli) revealed that the previous model of subvertical intrusions-related collapse is a likely mechanism affecting small-scale steeply-sloping basaltic volcanoes like Stromboli. Furthermore, our field study combined to modelling results confirms the importance of shallow-dipping intrusions in the morpho-structural evolution of large gently-sloping basaltic volcanoes like Piton de la Fournaise, Etna and Kilauea, with particular regards to flank instability, which can cause catastrophic tsunamis.

Methods for Clearance Influence Analysis in Planar and Spatial Mechanisms

This doctoral dissertation presents a new method to asses the influence of clearancein the kinematic pairs on the configuration of planar and spatial mechanisms. The subject has been widely investigated in both past and present scientific literature, and is approached in different ways: a static/kinetostatic way, which looks for the clearance take-up due to the external loads on the mechanism; a probabilistic way, which expresses clearance-due displacements using probability density functions; a dynamic way, which evaluates dynamic effects like the actual forces in the pairs caused by impacts, or the consequent vibrations. This dissertation presents a new method to approach the problem of clearance. The problem is studied from a purely kinematic perspective. With reference to a given mechanism configuration, the pose (position and orientation) error of the mechanism link of interest is expressed as a vector function of the degrees of freedom introduced in each pair by clearance: the presence of clearance in a kinematic pair, in facts, causes the actual pair to have more degrees of freedom than the theoretical clearance-free one. The clearance-due degrees of freedom are bounded by the pair geometry. A proper modelling of clearance-affected pairs allows expressing such bounding through analytical functions. It is then possible to study the problem as a maximization problem, where a continuous function (the pose error of the link of interest) subject to some constraints (the analytical functions bounding clearance- due degrees of freedom) has to be maximize. Revolute, prismatic, cylindrical, and spherical clearance-affected pairs have been analytically modelled; with reference to mechanisms involving such pairs, the solution to the maximization problem has been obtained in a closed form.

Polymer-liquid crystal interface: a molecular dynamics study

Liquid crystals (LCs) are an interesting class of soft condensed matter systems characterized by an unusual combination of fluidity and long-range order, mainly known for their applications in displays (LCDs). However, the interest in LC continues to grow pushed by their application in new technologies in medicine, optical imaging, micro and nano technologies etc. In LCDs uniaxial alignment of LCs is mainly achieved by a rubbing process. During this treatment, the surfaces of polymer coated display substrates are rubbed in one direction by a rotating cylinder covered with a rubbing cloth. Basically, LC alignment involves two possible aligning directions: uniaxial planar (homogeneous) and vertical (homeotropic) to the display substrate. An interesting unresolved question concerning LCs regards the origin of their alignment on rubbed surfaces, and in particular on the polymeric ones used in the display industry. Most studies have shown that LCs on the surface of the rubbed polymer film layer are lying parallel to the rubbing direction. In these systems, micrometric grooves are generated on the film surface along the rubbing direction and also the polymer chains are stretched in this direction. Both the parallel aligned microgrooves and the polymer chains at the film surface may play a role in the LC alignment and it is not easy to quantify the effect of each contribution. The work described in this thesis is an attempt to find new microscopic evidences on the origin of LC alignment on polymeric surfaces through molecular dynamics (MD) simulations, which allow the investigation of the phenomenon with atomic detail. The importance of the arrangement of the polymeric chains in LCs alignment was studied by performing MD simulations of a thin film of a typical nematic LC, 4-cyano-4’-pentylbiphenyl (5CB), in contact with two different polymers: poly(methyl methacrylate)(PMMA) and polystyrene (PS). At least four factors are believed to influence the LC alignment: 1. the interactions of LCs with the backbone vinyl chains; 2. the interactions of LCs with the oriented side groups; 3. the anisotropic interactions of LCs with nanometric grooves; 4. the presence of static surface charges. Here we exclude the effect of microgrooves and of static surface charges from our virtual experiment, by using flat and neutral polymer surfaces, with the aim of isolating the chemical driving factors influencing the alignment of LC phases on polymeric surfaces.

Strategie di personalizzazione della terapia di resincronizzazione cardiaca

Solo il 60% dei candidati alla resincronizzazione cardiaca risponde in termini di rimodellamento ventricolare inverso che è il più forte predittore di riduzione della mortalità e delle ospedalizzazioni. Due cause possibili della mancata risposta sono la programmazione del dispositivo e i limiti dell’ approccio transvenoso. Nel corso degli anni di dottorato ho effettuato tre studi per ridurre il numero di non responder. Il primo studio valuta il ritardo interventricolare. Al fine di ottimizzare le risorse e fornire un reale beneficio per il paziente ho ricercato la presenza di predittori di ritardo interventricolare diverso dal simultaneo, impostato nella programmazione di base. L'unico predittore è risultato essere l’ intervallo QRS> 160 ms, quindi ho proposto una flow chart per ottimizzare solo i pazienti che avranno nella programmazione ottimale un intervallo interventricolare non simultaneo. Il secondo lavoro valuta la fissazione attiva del ventricolo sinistro con stent. I dislocamenti, la soglia alta di stimolazione del miocardio e la stimolazione del nervo frenico sono tre problematiche che limitano la stimolazione biventricolare. Abbiamo analizzato più di 200 angiografie per vedere le condizioni anatomiche predisponenti la dislocazione del catetere. Prospetticamente abbiamo deciso di utilizzare uno stent per fissare attivamente il catetere ventricolare sinistro in tutti i pazienti che presentavano le caratteristiche anatomiche favorenti la dislocazione. Non ci sono più state dislocazioni, c’è stata una migliore risposta in termini di rimodellamento ventricolare inverso e non ci sono state modifiche dei parametri elettrici del catetere. Il terzo lavoro ha valutato sicurezza ed efficacia della stimolazione endoventricolare sinistra. Abbiamo impiantato 26 pazienti giudicati non responder alla terapia di resincronizzazione cardiaca. La procedura è risultata sicura, il rischio di complicanze è simile alla stimolazione biventricolare classica, ed efficace nell’arrestare la disfunzione ventricolare sinistra e / o migliorare gli effetti clinici in un follow-up medio.

Conductive Polymer Composites

In recent years, nanotechnologies have led to the production of materials with new and sometimes unexpected qualities through the manipulation of nanoscale components. This research aimed primarily to the study of the correlation between hierarchical structures of hybrid organic-inorganic materials such as conductive polymer composites (CPCs). Using a bottom-up methodology, we could synthesize a wide range of inorganic nanometric materials with a high degree of homogeneity and purity, such as thiol capped metal nanoparticles, stoichiometric geomimetic chrysotile nanotubes and metal dioxide nanoparticles. It was also possible to produce inorganic systems formed from the interaction between the synthesized materials. These synthesized materials and others like multiwalled carbon nanotubes and grapheme oxide were used to produce conductive polymer composites. Electrospinning causes polymer fibers to become elongated using an electric field. This technique was used to produce fibers with a nanometric diameter of a polymer blend based on two different intrinsically conducting polymers polymers (ICPs): polyaniline (PANI) and poly(3-hexylthiophene) (P3HT). Using different materials as second phase in the initial electrospun polymer fibers caused significant changes to the material hierarchical structure, leading to the creation of CPCs with modified electrical properties. Further study of the properties of these new materials resulted in a better understanding of the electrical conductivity mechanisms in these electrospun materials.

Störungstheoretische Behandlung relativistischer Korrekturen zu Gleichgewichtsstrukturen im Rahmen der Coupled-Cluster-Theorie

Quantenchemische Untersuchungen von Atomen und Molekülen haben in den letzten Jahren durch die systematische Erweiterung der Methoden und Computerresourcen zunehmend für die Interpretation und Vorhersage experimenteller Ergebnisse an Bedeutung gewonnen. Relativistische Effekte in der Chemie werden zum Beispiel für die gelbe Farbe von Gold und den flüssigen Aggregatzustand von Quecksilber verantwortlich gemacht und müssen daher in quantenchemischen Rechnungen berücksichtigt werden. Relativistische Effekte sind bei leichten Elementen oft so klein, daß sie in vielen quantenchemischen Betrachtungen vernachlässigt werden. Dennoch sind es gerade diese Beiträge, die verbleibende Abweichungen von noch so genauen nichtrelativistischen Rechnungen von ebenso genauen experimentellen Ergebnissen ausmachen können. Relativistische Effekte können auf viele Arten in quantenchemischen Rechnungen berücksichtigt werden. Eine Möglichkeit ist die Störungstheorie. Ein derartiger Ansatz ist die Mass-velocity-Darwin-Näherung, ein anderer die Direkte Störungstheorie. Hier entspricht die relativistische Energiekorrektur erster Ordnung der ersten Ableitung der Energie nach einem relativistischen Störparameter. Für eine Bestimmung der Gleichgewichtsstruktur eines Moleküls müssen die Kräfte auf die Atomkerne bestimmt werden. Diese entsprechen einer ersten Ableitung der Gesamtenergie nach den Kernkoordinaten. Eine Einbeziehung der relativistischen Effekte auf diese Kräfte erfordert daher die gemischte zweite Ableitung der Energie nach dem relativistischen Störparameter und den Kernkoordinaten. Diese relativistischen Korrekturen wurden in dem quantenchemischen Programmpaket ACES2 implementiert. Ein Resultat dieser Arbeit ist, daß nun erstmalig eine Implementierung analytischer Gradienten für die Berechnung relativistischer Korrekturen zu Strukturparametern mit Hilfe der relativistischen Störungstheorie für den Coupled-Cluster-Ansatz bereit steht. Die Coupled-Cluster-Theorie eignet sich besonders gut für die hochgenaue Vorhersage von molekularen Eigenschaften, wie der Gleichgewichtsstruktur. Im Rahmen dieser Arbeit wurde die Basissatzabhängigkeit der relativistischen Beiträge zu Energien, Strukturparametern und harmonischen Schwingungsfrequenzen im Detail untersucht. Für die hier untersuchten Moleküle sind die relativistischen Effekte und Effekte aufgrund der Elektronenkorrelation nicht additiv, so verkürzt die Berücksichtigung relativistischer Effekte bei Hartree-Fock-Rechnungen die Bindung in den Hydrogenhalogeniden, während die Einbeziehung der Elektronenkorrelation durch CCSD(T)-Rechnungen zu einer verlängerten Bindung im Fluorwasserstoff und weniger stark ausgeprägten Korrekturen im Chlor- und Bromwasserstoff führt. Für die anderen hier untersuchten mehratomigen Moleküle findet sich kein einheitlicher Trend; dies unterstreicht die Notwendigkeit expliziter Rechnungen. Damit steht ein leistungsfähiges und vielseitiges Werkzeug für die Berechnung relativistischer Korrekturen auf verschiedenste molekulare Eigenschaften zur Verfügung, das mit modernen, systematisch verbesserbaren quantenchemischen Methoden verknüpft ist. Hiermit ist es möglich, hochgenaue Rechnungen zur Vorhersage und Interpretation von Experimenten durchzuführen.

Crystal phases and glassy dynamics in monodisperse hard ellipsoids

We have performed Monte Carlo and molecular dynamics simulations of suspensions of monodisperse, hard ellipsoids of revolution. Hard-particle models play a key role in statistical mechanics. They are conceptually and computationally simple, and they offer insight into systems in which particle shape is important, including atomic, molecular, colloidal, and granular systems. In the high density phase diagram of prolate hard ellipsoids we have found a new crystal, which is more stable than the stretched FCC structure proposed previously . The new phase, SM2, has a simple monoclinic unit cell containing a basis of two ellipsoids with unequal orientations. The angle of inclination is very soft for length-to-width (aspect) ratio l/w=3, while the other angles are not. A symmetric state of the unit cell exists, related to the densest-known packings of ellipsoids; it is not always the stable one. Our results remove the stretched FCC structure for aspect ratio l/w=3 from the phase diagram of hard, uni-axial ellipsoids. We provide evidence that this holds between aspect ratios 3 and 6, and possibly beyond. Finally, ellipsoids in SM2 at l/w=1.55 exhibit end-over-end flipping, warranting studies of the cross-over to where this dynamics is not possible. Secondly, we studied the dynamics of nearly spherical ellipsoids. In equilibrium, they show a first-order transition from an isotropic phase to a rotator phase, where positions are crystalline but orientations are free. When over-compressing the isotropic phase into the rotator regime, we observed super-Arrhenius slowing down of diffusion and relaxation, and signatures of the cage effect. These features of glassy dynamics are sufficiently strong that asymptotic scaling laws of the Mode-Coupling Theory of the glass transition (MCT) could be tested, and were found to apply. We found strong coupling of positional and orientational degrees of freedom, leading to a common value for the MCT glass-transition volume fraction. Flipping modes were not slowed down significantly. We demonstrated that the results are independent of simulation method, as predicted by MCT. Further, we determined that even intra-cage motion is cooperative. We confirmed the presence of dynamical heterogeneities associated with the cage effect. The transit between cages was seen to occur on short time scales, compared to the time spent in cages; but the transit was shown not to involve displacements distinguishable in character from intra-cage motion. The presence of glassy dynamics was predicted by molecular MCT (MMCT). However, as MMCT disregards crystallization, a test by simulation was required. Glassy dynamics is unusual in monodisperse systems. Crystallization typically intervenes unless polydispersity, network-forming bonds or other asymmetries are introduced. We argue that particle anisometry acts as a sufficient source of disorder to prevent crystallization. This sheds new light on the question of which ingredients are required for glass formation.

Multicomponent nanodevices based on molecular and nanocrystal moieties

Nanoscience is an emerging and fast-growing field of science with the aim of manipulating nanometric objects with dimension below 100 nm. Top down approach is currently used to build these type of architectures (e.g microchips). The miniaturization process cannot proceed indefinitely due to physical and technical limitations. Those limits are focusing the interest on the bottom-up approach and construction of nano-objects starting from “nano-bricks” like atoms, molecules or nanocrystals. Unlike atoms, molecules can be “fully programmable” and represent the best choice to build up nanostructures. In the past twenty years many examples of functional nano-devices able to perform simple actions have been reported. Nanocrystals which are often considered simply nanostructured materials, can be active part in the development of those nano-devices, in combination with functional molecules. The object of this dissertation is the photophysical and photochemical investigation of nano-objects bearing molecules and semiconductor nanocrystals (QDs) as components. The first part focuses on the characterization of a bistable rotaxane. This study, in collaboration with the group of Prof. J.F. Stoddart (Northwestern University, Evanston, Illinois, USA) who made the synthesis of the compounds, shows the ability of this artificial machine to operate as bistable molecular-level memory under kinetic control. The second part concerns the study of the surface properties of luminescent semiconductor nanocrystals (QDs) and in particular the effect of acid and base on the spectroscopical properties of those nanoparticles. In this section is also reported the work carried out in the laboratory of Prof H. Mattoussi (Florida State University, Tallahassee, Florida, USA), where I developed a novel method for the surface decoration of QDs with lipoic acid-based ligands involving the photoreduction of the di-thiolane moiety.

Valutazione della risposta sismica di edifici mediante modelli anelastici equivalenti

Nella presente tesi è proposta una metodologia per lo studio e la valutazione del comportamento sismico di edifici a telaio. Il metodo prevede la realizzazione di analisi non-lineari su modelli equivalenti MDOF tipo stick, in accordo alla classificazione data nel report FEMA 440. Gli step per l’applicazione del metodo sono descritti nella tesi. Per la validazione della metodologia si sono utilizzati confronti con analisi time-history condotte su modelli tridimensionali dettagliati delle strutture studiate (detailed model). I parametri ingegneristici considerati nel confronto, nell’ottica di utilizzare il metodo proposto in un approccio del tipo Displacement-Based Design sono lo spostamento globale in sommità, gli spostamenti di interpiano, le forze di piano e la forza totale alla base. I risultati delle analisi condotte sui modelli stick equivalenti, mostrano una buona corrispondenza, ottima in certi casi, con quelli delle analisi condotte sui modelli tridimensionali dettagliati. Le time-history realizzate sugli stick model permettono però, un consistente risparmio in termini di onere computazionale e di tempo per il post-processing dei risultati ottenuti.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

The APSEL4D Monolithic Active Pixel Sensor and its Usage in a Single Electron Interference Experiment

We have realized a Data Acquisition chain for the use and characterization of APSEL4D, a 32 x 128 Monolithic Active Pixel Sensor, developed as a prototype for frontier experiments in high energy particle physics. In particular a transition board was realized for the conversion between the chip and the FPGA voltage levels and for the signal quality enhancing. A Xilinx Spartan-3 FPGA was used for real time data processing, for the chip control and the communication with a Personal Computer through a 2.0 USB port. For this purpose a firmware code, developed in VHDL language, was written. Finally a Graphical User Interface for the online system monitoring, hit display and chip control, based on windows and widgets, was realized developing a C++ code and using Qt and Qwt dedicated libraries. APSEL4D and the full acquisition chain were characterized for the first time with the electron beam of the transmission electron microscope and with 55Fe and 90Sr radioactive sources. In addition, a beam test was performed at the T9 station of the CERN PS, where hadrons of momentum of 12 GeV/c are available. The very high time resolution of APSEL4D (up to 2.5 Mfps, but used at 6 kfps) was fundamental in realizing a single electron Young experiment using nanometric double slits obtained by a FIB technique. On high statistical samples, it was possible to observe the interference and diffractions of single isolated electrons traveling inside a transmission electron microscope. For the first time, the information on the distribution of the arrival time of the single electrons has been extracted.

In vivo evaluation of the translations of the gleno-humeral joint using Magnetic resonance imaging

Gleno-humeral joint (GHJ) is the most mobile joint of the human body. This is related to theincongr uence between the large humeral head articulating with the much smaller glenoid (ratio 3:1). The GHJ laxity is the ability of the humeral head to be passively translated on the glenoid fossa and, when physiological, it guarantees the normal range of motion of the joint. Three-dimensional GHJ linear displacements have been measured, both in vivo and in vitro by means of different instrumental techniques. In vivo gleno-humeral displacements have been assessed by means of stereophotogrammetry, electromagnetic tracking sensors, and bio-imaging techniques. Both stereophotogrammetric systems and electromagnetic tracking devices, due to the deformation of the soft tissues surrounding the bones, are not capable to accurately assess small displacements, such as gleno-humeral joint translations. The bio-imaging techniques can ensure for an accurate joint kinematic (linear and angular displacement) description, but, due to the radiation exposure, most of these techniques, such as computer tomography or fluoroscopy, are invasive for patients. Among the bioimaging techniques, an alternative which could provide an acceptable level of accuracy and that is innocuous for patients is represented by magnetic resonance imaging (MRI). Unfortunately, only few studies have been conducted for three-dimensional analysis and very limited data is available in situations where preset loads are being applied. The general aim of this doctoral thesis is to develop a non-invasive methodology based on open-MRI for in-vivo evaluation of the gleno-humeral translation components in healthy subjects under the application of external loads.

Application of Submerged Grouted Anchors in Sheet Pile Quay Walls

Sheet pile walls are one of the oldest earth retention systems utilized in civil engineering projects. They are used for various purposes; such as excavation support system, cofferdams, cut-off walls under dams, slope stabilization, waterfront structures, and flood walls. Sheet pile walls are one of the most common types of quay walls used in port construction. The worldwide increases in utilization of large ships for transportation have created an urgent need of deepening the seabed within port areas and consequently the rehabilitation of its wharfs. Several methods can be used to increase the load-carrying capacity of sheet-piling walls. The use of additional anchored tie rods grouted into the backfill soil and arranged along the exposed wall height is one of the most practical and appropriate solutions adopted for stabilization and rehabilitation of the existing quay wall. The Ravenna Port Authority initiated a project to deepen the harbor bottom at selected wharves. An extensive parametric study through the finite element program, PLAXIS 2D, version 2012 was carried out to investigate the enhancement of using submerged grouted anchors technique on the load response of sheet-piling quay wall. The influence of grout-ties area, length of grouted body, anchor inclination and anchor location were considered and evaluated due to the effect of different system parameters. Also a comparative study was conducted by Plaxis 2D and 3D program to investigate the behavior of these sheet pile quay walls in terms of horizontal displacements induced along the sheet pile wall and ground surface settlements as well as the anchor force and calculated factor of safety. Finally, a comprehensive study was carried out by using different constitutive models to simulate the mechanical behavior of the soil to investigate the effect of these two models (Mohr-Coulomb and Hardening Soil) on the behavior of these sheet pile quay walls.

The reconstruction of skeletal movement: the soft tissue artefact issue

In 3D human movement analysis performed using stereophotogrammetric systems and skin markers, bone pose can only be estimated in an indirect fashion. During a movement, soft tissue deformations make the markers move with respect to the underlying bone generating soft tissue artefact (STA). STA has devastating effects on bone pose estimation and its compensation remains an open question. The aim of this PhD thesis was to contribute to the solution of this crucial issue. Modelling STA using measurable trial-specific variables is a fundamental prerequisite for its removal from marker trajectories. Two STA model architectures are proposed. Initially, a thigh marker-level artefact model is presented. STA was modelled as a linear combination of joint angles involved in the movement. This model was calibrated using ex-vivo and in-vivo STA invasive measures. The considerable number of model parameters led to defining STA approximations. Three definitions were proposed to represent STA as a series of modes: individual marker displacements, marker-cluster geometrical transformations (MCGT), and skin envelope shape variations. Modes were selected using two criteria: one based on modal energy and another on the selection of modes chosen a priori. The MCGT allows to select either rigid or non-rigid STA components. It was also empirically demonstrated that only the rigid component affects joint kinematics, regardless of the non-rigid amplitude. Therefore, a model of thigh and shank STA rigid component at cluster-level was then defined. An acceptable trade-off between STA compensation effectiveness and number of parameters can be obtained, improving joint kinematics accuracy. The obtained results lead to two main potential applications: the proposed models can generate realistic STAs for simulation purposes to compare different skeletal kinematics estimators; and, more importantly, focusing only on the STA rigid component, the model attains a satisfactory STA reconstruction with less parameters, facilitating its incorporation in an pose estimator.

Spatial and temporal characterisation of ground deformation recorded by geodetic techniques

A critical point in the analysis of ground displacements time series is the development of data driven methods that allow the different sources that generate the observed displacements to be discerned and characterised. A widely used multivariate statistical technique is the Principal Component Analysis (PCA), which allows reducing the dimensionality of the data space maintaining most of the variance of the dataset explained. Anyway, PCA does not perform well in finding the solution to the so-called Blind Source Separation (BSS) problem, i.e. in recovering and separating the original sources that generated the observed data. This is mainly due to the assumptions on which PCA relies: it looks for a new Euclidean space where the projected data are uncorrelated. The Independent Component Analysis (ICA) is a popular technique adopted to approach this problem. However, the independence condition is not easy to impose, and it is often necessary to introduce some approximations. To work around this problem, I use a variational bayesian ICA (vbICA) method, which models the probability density function (pdf) of each source signal using a mix of Gaussian distributions. This technique allows for more flexibility in the description of the pdf of the sources, giving a more reliable estimate of them. Here I present the application of the vbICA technique to GPS position time series. First, I use vbICA on synthetic data that simulate a seismic cycle (interseismic + coseismic + postseismic + seasonal + noise) and a volcanic source, and I study the ability of the algorithm to recover the original (known) sources of deformation. Secondly, I apply vbICA to different tectonically active scenarios, such as the 2009 L'Aquila (central Italy) earthquake, the 2012 Emilia (northern Italy) seismic sequence, and the 2006 Guerrero (Mexico) Slow Slip Event (SSE).