Kinematic models of interseismic deformation from inversion of GPS and InSAR measurements to estimate fault parameters and coupling degree

We have used kinematic models in two Italian regions to reproduce surface interseismic velocities obtained from InSAR and GPS measurements. We have considered a Block modeling, BM, approach to evaluate which fault system is actively accommodating the occurring deformation in both considered areas. We have performed a study for the Umbria-Marche Apennines, obtaining that the tectonic extension observed by GPS measurements is explained by the active contribution of at least two fault systems, one of which is the Alto Tiberina fault, ATF. We have estimated also the interseismic coupling distribution for the ATF using a 3D surface and the result shows an interesting correlation between the microseismicity and the uncoupled fault portions. The second area analyzed concerns the Gargano promontory for which we have used jointly the available InSAR and GPS velocities. Firstly we have attached the two datasets to the same terrestrial reference frame and then using a simple dislocation approach, we have estimated the best fault parameters reproducing the available data, providing a solution corresponding to the Mattinata fault. Subsequently we have considered within a BM analysis both GPS and InSAR datasets in order to evaluate if the Mattinata fault may accommodate the deformation occurring in the central Adriatic due to the relative motion between the North-Adriatic and South-Adriatic plates. We obtain that the deformation occurring in that region should be accommodated by more that one fault system, that is however difficult to detect since the poor coverage of geodetic measurement offshore of the Gargano promontory. Finally we have performed also the estimate of the interseismic coupling distribution for the Mattinata fault, obtaining a shallow coupling pattern. Both of coupling distributions found using the BM approach have been tested by means of resolution checkerboard tests and they demonstrate that the coupling patterns depend on the geodetic data positions.

Towards the 3D attenuation imaging of active volcanoes: methods and tests on real and simulated data

The purpose of my PhD thesis has been to face the issue of retrieving a three dimensional attenuation model in volcanic areas. To this purpose, I first elaborated a robust strategy for the analysis of seismic data. This was done by performing several synthetic tests to assess the applicability of spectral ratio method to our purposes. The results of the tests allowed us to conclude that: 1) spectral ratio method gives reliable differential attenuation (dt*) measurements in smooth velocity models; 2) short signal time window has to be chosen to perform spectral analysis; 3) the frequency range over which to compute spectral ratios greatly affects dt* measurements. Furthermore, a refined approach for the application of spectral ratio method has been developed and tested. Through this procedure, the effects caused by heterogeneities of propagation medium on the seismic signals may be removed. The tested data analysis technique was applied to the real active seismic SERAPIS database. It provided a dataset of dt* measurements which was used to obtain a three dimensional attenuation model of the shallowest part of Campi Flegrei caldera. Then, a linearized, iterative, damped attenuation tomography technique has been tested and applied to the selected dataset. The tomography, with a resolution of 0.5 km in the horizontal directions and 0.25 km in the vertical direction, allowed to image important features in the off-shore part of Campi Flegrei caldera. High QP bodies are immersed in a high attenuation body (Qp=30). The latter is well correlated with low Vp and high Vp/Vs values and it is interpreted as a saturated marine and volcanic sediments layer. High Qp anomalies, instead, are interpreted as the effects either of cooled lava bodies or of a CO2 reservoir. A pseudo-circular high Qp anomaly was detected and interpreted as the buried rim of NYT caldera.

Statistical analysis of the error associated with the simplification of the stratigraphy in geotechnical models

The uncertainties in the determination of the stratigraphic profile of natural soils is one of the main problems in geotechnics, in particular for landslide characterization and modeling. The study deals with a new approach in geotechnical modeling which relays on a stochastic generation of different soil layers distributions, following a boolean logic – the method has been thus called BoSG (Boolean Stochastic Generation). In this way, it is possible to randomize the presence of a specific material interdigitated in a uniform matrix. In the building of a geotechnical model it is generally common to discard some stratigraphic data in order to simplify the model itself, assuming that the significance of the results of the modeling procedure would not be affected. With the proposed technique it is possible to quantify the error associated with this simplification. Moreover, it could be used to determine the most significant zones where eventual further investigations and surveys would be more effective to build the geotechnical model of the slope. The commercial software FLAC was used for the 2D and 3D geotechnical model. The distribution of the materials was randomized through a specifically coded MatLab program that automatically generates text files, each of them representing a specific soil configuration. Besides, a routine was designed to automate the computation of FLAC with the different data files in order to maximize the sample number. The methodology is applied with reference to a simplified slope in 2D, a simplified slope in 3D and an actual landslide, namely the Mortisa mudslide (Cortina d’Ampezzo, BL, Italy). However, it could be extended to numerous different cases, especially for hydrogeological analysis and landslide stability assessment, in different geological and geomorphological contexts.

Models and Algorihtm for the Optimization of Real-World Routing and Logistics Problems

Logistics involves planning, managing, and organizing the flows of goods from the point of origin to the point of destination in order to meet some requirements. Logistics and transportation aspects are very important and represent a relevant costs for producing and shipping companies, but also for public administration and private citizens. The optimization of resources and the improvement in the organization of operations is crucial for all branches of logistics, from the operation management to the transportation. As we will have the chance to see in this work, optimization techniques, models, and algorithms represent important methods to solve the always new and more complex problems arising in different segments of logistics. Many operation management and transportation problems are related to the optimization class of problems called Vehicle Routing Problems (VRPs). In this work, we consider several real-world deterministic and stochastic problems that are included in the wide class of the VRPs, and we solve them by means of exact and heuristic methods. We treat three classes of real-world routing and logistics problems. We deal with one of the most important tactical problems that arises in the managing of the bike sharing systems, that is the Bike sharing Rebalancing Problem (BRP). We propose models and algorithms for real-world earthwork optimization problems. We describe the 3DP process and we highlight several optimization issues in 3DP. Among those, we define the problem related to the tool path definition in the 3DP process, the 3D Routing Problem (3DRP), which is a generalization of the arc routing problem. We present an ILP model and several heuristic algorithms to solve the 3DRP.

Computational Modelling of Cardiac Electrophysiology: from Cell to Bedside

Heart diseases are the leading cause of death worldwide, both for men and women. However, the ionic mechanisms underlying many cardiac arrhythmias and genetic disorders are not completely understood, thus leading to a limited efficacy of the current available therapies and leaving many open questions for cardiac electrophysiologists. On the other hand, experimental data availability is still a great issue in this field: most of the experiments are performed in vitro and/or using animal models (e.g. rabbit, dog and mouse), even when the final aim is to better understand the electrical behaviour of in vivo human heart either in physiological or pathological conditions. Computational modelling constitutes a primary tool in cardiac electrophysiology: in silico simulations, based on the available experimental data, may help to understand the electrical properties of the heart and the ionic mechanisms underlying a specific phenomenon. Once validated, mathematical models can be used for making predictions and testing hypotheses, thus suggesting potential therapeutic targets. This PhD thesis aims to apply computational cardiac modelling of human single cell action potential (AP) to three clinical scenarios, in order to gain new insights into the ionic mechanisms involved in the electrophysiological changes observed in vitro and/or in vivo. The first context is blood electrolyte variations, which may occur in patients due to different pathologies and/or therapies. In particular, we focused on extracellular Ca2+ and its effect on the AP duration (APD). The second context is haemodialysis (HD) therapy: in addition to blood electrolyte variations, patients undergo a lot of other different changes during HD, e.g. heart rate, cell volume, pH, and sympatho-vagal balance. The third context is human hypertrophic cardiomyopathy (HCM), a genetic disorder characterised by an increased arrhythmic risk, and still lacking a specific pharmacological treatment.

Kinematics of local and high-z galaxies through 3D modeling of emission-line datacubes

The kinematics is a fundamental tool to infer the dynamical structure of galaxies and to understand their formation and evolution. Spectroscopic observations of gas emission lines are often used to derive rotation curves and velocity dispersions. It is however difficult to disentangle these two quantities in low spatial-resolution data because of beam smearing. In this thesis, we present 3D-Barolo, a new software to derive the gas kinematics of disk galaxies from emission-line data-cubes. The code builds tilted-ring models in the 3D observational space and compares them with the actual data-cubes. 3D-Barolo works with data at a wide range of spatial resolutions without being affected by instrumental biases. We use 3D-Barolo to derive rotation curves and velocity dispersions of several galaxies in both the local and the high-redshift Universe. We run our code on HI observations of nearby galaxies and we compare our results with 2D traditional approaches. We show that a 3D approach to the derivation of the gas kinematics has to be preferred to a 2D approach whenever a galaxy is resolved with less than about 20 elements across the disk. We moreover analyze a sample of galaxies at z~1, observed in the H-alpha line with the KMOS/VLT spectrograph. Our 3D modeling reveals that the kinematics of these high-z systems is comparable to that of local disk galaxies, with steeply-rising rotation curves followed by a flat part and H-alpha velocity dispersions of 15-40 km/s over the whole disks. This evidence suggests that disk galaxies were already fully settled about 7-8 billion years ago. In summary, 3D-Barolo is a powerful and robust tool to separate physical and instrumental effects and to derive a reliable kinematics. The analysis of large samples of galaxies at different redshifts with 3D-Barolo will provide new insights on how galaxies assemble and evolve throughout cosmic time.

2D and 3D numerical modelling of multilayer detachment folding and salt tectonics

Numerical modelling was performed to study the dynamics of multilayer detachment folding and salt tectonics. In the case of multilayer detachment folding, analytically derived diagrams show several folding modes, half of which are applicable to crustal scale folding. 3D numerical simulations are in agreement with 2D predictions, yet fold interactions result in complex fold patterns. Pre-existing salt diapirs change folding patterns as they localize the initial deformation. If diapir spacing is much smaller than the dominant folding wavelength, diapirs appear in fold synclines or limbs.rnNumerical models of 3D down-building diapirism show that sedimentation rate controls whether diapirs will form and influences the overall patterns of diapirism. Numerical codes were used to retrodeform modelled salt diapirs. Reverse modelling can retrieve the initial geometries of a 2D Rayleigh-Taylor instability with non-linear rheologies. Although intermediate geometries of down-built diapirs are retrieved, forward and reverse modelling solutions deviate. rnFinally, the dynamics of fold-and-thrusts belts formed over a tilted viscous detachment is studied and it is demonstrated that mechanical stratigraphy has an impact on the deformation style, switching from thrust- to folding-dominated. The basal angle of the detachment controls the deformation sequence of the fold-and-thrust belt and results are consistent with critical wedge theory.rn

3D-ultrastructure, functions and stress responses of rhogocytes from the freshwater gastropods Biomphalaria glabrata and Lymnaea stagnalis

Rhogocytes, also termed ‘pore cells’, exist free in the hemolymph or embedded in the connective tissue of different body parts of molluscs, notably gastropods. These unique cells can be round, elongated or irregularly shaped, and up to 30 μm in diameter. Their hallmark is the so-called slit apparatus: i.e. pocket-like invaginations of the plasma membrane creating extracellular lacunae, bridged by cytoplasmic bars. These bars form distinctive slits of ca. 20 nm width. A slit diaphragm composed of proteins establishes a molecular sieve with holes of 20 x 20 nm. Different functions have been assigned to this special molluscan cell type, notably biosynthesis of the hemolymph respiratory protein hemocyanin. It has further been proposed, but not proven, that in the case of red-blooded snail species rhogocytes might synthesize the hemoglobin. However, the secretion pathway of these hemolymph proteins, and the functional role of the enigmatic slit apparatus remained unclear. Additionally proposed functions of rhogocytes, such as heavy metal detoxification or hemolymph protein degradation, are also not well studied. This work provides more detailed electron microscopical, histological and immunobiochemical information on the structure and function of rhogocytes of the freshwater snails Biomphalaria glabrata and Lymnaea stagnalis. By in situ hybridization on mantle tissues, it proves that B. glabrata rhogocytes synthesize hemoglobin and L. stagnalis rhogocytes synthesize hemocyanin. Hemocyanin is present, in endoplasmic reticulum lacunae and in vesicles, as individual molecules or pseudo-crystalline arrays. The first 3D reconstructions of rhogocytes are provided by means of electron tomography and show unprecedented details of the slit apparatus. A highly dense material in the cytoplasmic bars close to the diaphragmatic slits was shown, by immunogold labeling, to contain actin. By immunofluorescence microscopy, the protein nephrin was localized at the periphery of rhogocytes. The presence of both proteins in the slit apparatus supports the previous hypothesis, hitherto solely based on similarities of the ultrastructure, that the molluscan rhogocytes are phylogenetically related to mammalian podocytes and insect nephrocytes. A possible secretion pathway of respiratory proteins that includes a transfer mechanism of vesicles through the diaphragmatic slits is proposed and discussed. We also studied, by electron microscopy, the reaction of rhogocytes in situ to two forms of animal stress: deprivation of food and cadmium contamination of the tank water. Significant cellular reactions to both stressors were observed and documented. Notably, the slit apparatus surface and the number of electron-dense cytoplasmic vesicles increased in response to cadmium stress. Food deprivation led to an increase in hemocyanin production. These observations are also discussed in the framework of using such animals as potential environmental biomarkers.

Visualization and validation of (Q)SAR models

Analyzing and modeling relationships between the structure of chemical compounds, their physico-chemical properties, and biological or toxic effects in chemical datasets is a challenging task for scientific researchers in the field of cheminformatics. Therefore, (Q)SAR model validation is essential to ensure future model predictivity on unseen compounds. Proper validation is also one of the requirements of regulatory authorities in order to approve its use in real-world scenarios as an alternative testing method. However, at the same time, the question of how to validate a (Q)SAR model is still under discussion. In this work, we empirically compare a k-fold cross-validation with external test set validation. The introduced workflow allows to apply the built and validated models to large amounts of unseen data, and to compare the performance of the different validation approaches. Our experimental results indicate that cross-validation produces (Q)SAR models with higher predictivity than external test set validation and reduces the variance of the results. Statistical validation is important to evaluate the performance of (Q)SAR models, but does not support the user in better understanding the properties of the model or the underlying correlations. We present the 3D molecular viewer CheS-Mapper (Chemical Space Mapper) that arranges compounds in 3D space, such that their spatial proximity reflects their similarity. The user can indirectly determine similarity, by selecting which features to employ in the process. The tool can use and calculate different kinds of features, like structural fragments as well as quantitative chemical descriptors. Comprehensive functionalities including clustering, alignment of compounds according to their 3D structure, and feature highlighting aid the chemist to better understand patterns and regularities and relate the observations to established scientific knowledge. Even though visualization tools for analyzing (Q)SAR information in small molecule datasets exist, integrated visualization methods that allows for the investigation of model validation results are still lacking. We propose visual validation, as an approach for the graphical inspection of (Q)SAR model validation results. New functionalities in CheS-Mapper 2.0 facilitate the analysis of (Q)SAR information and allow the visual validation of (Q)SAR models. The tool enables the comparison of model predictions to the actual activity in feature space. Our approach reveals if the endpoint is modeled too specific or too generic and highlights common properties of misclassified compounds. Moreover, the researcher can use CheS-Mapper to inspect how the (Q)SAR model predicts activity cliffs. The CheS-Mapper software is freely available at http://ches-mapper.org.

A comparison between 2D and 3D cephalometry on CBCT scans of human skulls

The purpose of this study was to evaluate whether measurements on conventional cephalometric radiographs are comparable with 3D measurements on 3D models of human skulls, derived from cone beam CT (CBCT) data. A CBCT scan and a conventional cephalometric radiograph were made of 40 dry skulls. Standard cephalometric software was used to identify landmarks on both the 2D images and the 3D models. The same operator identified 17 landmarks on the cephalometric radiographs and on the 3D models. All images and 3D models were traced five times with a time-interval of 1 week and the mean value of repeated measurements was used for further statistical analysis. Distances and angles were calculated. Intra-observer reliability was good for all measurements. The reproducibility of the measurements on the conventional cephalometric radiographs was higher compared with the reproducibility of measurements on the 3D models. For a few measurements a clinically relevant difference between measurements on conventional cephalometric radiographs and 3D models was found. Measurements on conventional cephalometric radiographs can differ significantly from measurements on 3D models of the same skull. The authors recommend that 3D tracings for longitudinal research are not used in cases were there are only 2D records from the past.

Development of replication-defective lymphocytic choriomeningitis virus vectors for the induction of potent CD8+ T cell immunity

Lymphocytic choriomeningitis virus (LCMV) exhibits natural tropism for dendritic cells and represents the prototypic infection that elicits protective CD8(+) T cell (cytotoxic T lymphocyte (CTL)) immunity. Here we have harnessed the immunobiology of this arenavirus for vaccine delivery. By using producer cells constitutively synthesizing the viral glycoprotein (GP), it was possible to replace the gene encoding LCMV GP with vaccine antigens to create replication-defective vaccine vectors. These rLCMV vaccines elicited CTL responses that were equivalent to or greater than those elicited by recombinant adenovirus 5 or recombinant vaccinia virus in their magnitude and cytokine profiles, and they exhibited more effective protection in several models. In contrast to recombinant adenovirus 5, rLCMV failed to elicit vector-specific antibody immunity, which facilitated re-administration of the same vector for booster vaccination. In addition, rLCMV elicited T helper type 1 CD4+ T cell responses and protective neutralizing antibodies to vaccine antigens. These features, together with low seroprevalence in humans, suggest that rLCMV may show utility as a vaccine platform against infectious diseases and cancer.

3D bloodstain pattern analysis: ballistic reconstruction of the trajectories of blood drops and determination of the centres of origin of the bloodstains

For crime scene investigation in cases of homicide, the pattern of bloodstains at the incident site is of critical importance. The morphology of the bloodstain pattern serves to determine the approximate blood source locations, the minimum number of blows and the positioning of the victim. In the present work, the benefits of the three-dimensional bloodstain pattern analysis, including the ballistic approximation of the trajectories of the blood drops, will be demonstrated using two illustrative cases. The crime scenes were documented in 3D, using the non-contact methods digital photogrammetry, tachymetry and laser scanning. Accurate, true-to-scale 3D models of the crime scenes, including the bloodstain pattern and the traces, were created. For the determination of the areas of origin of the bloodstain pattern, the trajectories of up to 200 well-defined bloodstains were analysed in CAD and photogrammetry software. The ballistic determination of the trajectories was performed using ballistics software. The advantages of this method are the short preparation time on site, the non-contact measurement of the bloodstains and the high accuracy of the bloodstain analysis. It should be expected that this method delivers accurate results regarding the number and position of the areas of origin of bloodstains, in particular the vertical component is determined more precisely than using conventional methods. In both cases relevant forensic conclusions regarding the course of events were enabled by the ballistic bloodstain pattern analysis.

Getting in touch-3D printing in Forensic Imaging

With the increasing use of medical imaging in forensics, as well as the technological advances in rapid prototyping, we suggest combining these techniques to generate displays of forensic findings. We used computed tomography (CT), CT angiography, magnetic resonance imaging (MRI) and surface scanning with photogrammetry in conjunction with segmentation techniques to generate 3D polygon meshes. Based on these data sets, a 3D printer created colored models of the anatomical structures. Using this technique, we could create models of bone fractures, vessels, cardiac infarctions, ruptured organs as well as bitemark wounds. The final models are anatomically accurate, fully colored representations of bones, vessels and soft tissue, and they demonstrate radiologically visible pathologies. The models are more easily understood by laypersons than volume rendering or 2D reconstructions. Therefore, they are suitable for presentations in courtrooms and for educational purposes.

Personalized X-Ray Reconstruction of the Proximal Femur via Intensity-Based Non-rigid 2D-3D Registration

This paper presents a new approach for reconstructing a patient-specific shape model and internal relative intensity distribution of the proximal femur from a limited number (e.g., 2) of calibrated C-arm images or X-ray radiographs. Our approach uses independent shape and appearance models that are learned from a set of training data to encode the a priori information about the proximal femur. An intensity-based non-rigid 2D-3D registration algorithm is then proposed to deformably fit the learned models to the input images. The fitting is conducted iteratively by minimizing the dissimilarity between the input images and the associated digitally reconstructed radiographs of the learned models together with regularization terms encoding the strain energy of the forward deformation and the smoothness of the inverse deformation. Comprehensive experiments conducted on images of cadaveric femurs and on clinical datasets demonstrate the efficacy of the present approach.

On the traces of XPD: cell cycle matters - untangling the genotype-phenotype relationship of XPD mutations

Abstract Mutations in the human gene coding for XPD lead to segmental progeria - the premature appearance of some of the phenotypes normally associated with aging - which may or may not be accompanied by increased cancer incidence. XPD is required for at least three different critical cellular functions: in addition to participating in the process of nucleotide excision repair (NER), which removes bulky DNA lesions, XPD also regulates transcription as part of the general transcription factor IIH (TFIIH) and controls cell cycle progression through its interaction with CAK, a pivotal activator of cyclin dependent kinases (CDKs). The study of inherited XPD disorders offers the opportunity to gain insights into the coordination of important cellular events and may shed light on the mechanisms that regulate the delicate equilibrium between cell proliferation and functional senescence, which is notably altered during physiological aging and in cancer. The phenotypic manifestations in the different XPD disorders are the sum of disturbances in the vital processes carried out by TFIIH and CAK. In addition, further TFIIH- and CAK-independent cellular activities of XPD may also play a role. This, added to the complex feedback networks that are in place to guarantee the coordination between cell cycle, DNA repair and transcription, complicates the interpretation of clinical observations. While results obtained from patient cell isolates as well as from murine models have been elementary in revealing such complexity, the Drosophila embryo has proven useful to analyze the role of XPD as a cell cycle regulator independently from its other cellular functions. Together with data from the biochemical and structural analysis of XPD and of the TFIIH complex these results combine into a new picture of the XPD activities that provides ground for a better understanding of the patophysiology of XPD diseases and for future development of diagnostic and therapeutic tools.