Three essays on major trends in a slow clockspeed industry : the case of industrial automation

The motivation for this research initiated from the abrupt rise and fall of minicomputers which were initially used both for industrial automation and business applications due to their significantly lower cost than their predecessors, the mainframes. Later industrial automation developed its own vertically integrated hardware and software to address the application needs of uninterrupted operations, real-time control and resilience to harsh environmental conditions. This has led to the creation of an independent industry, namely industrial automation used in PLC, DCS, SCADA and robot control systems. This industry employs today over 200'000 people in a profitable slow clockspeed context in contrast to the two mainstream computing industries of information technology (IT) focused on business applications and telecommunications focused on communications networks and hand-held devices. Already in 1990s it was foreseen that IT and communication would merge into one Information and communication industry (ICT). The fundamental question of the thesis is: Could industrial automation leverage a common technology platform with the newly formed ICT industry? Computer systems dominated by complex instruction set computers (CISC) were challenged during 1990s with higher performance reduced instruction set computers (RISC). RISC started to evolve parallel to the constant advancement of Moore's law. These developments created the high performance and low energy consumption System-on-Chip architecture (SoC). Unlike to the CISC processors RISC processor architecture is a separate industry from the RISC chip manufacturing industry. It also has several hardware independent software platforms consisting of integrated operating system, development environment, user interface and application market which enables customers to have more choices due to hardware independent real time capable software applications. An architecture disruption merged and the smartphone and tablet market were formed with new rules and new key players in the ICT industry. Today there are more RISC computer systems running Linux (or other Unix variants) than any other computer system. The astonishing rise of SoC based technologies and related software platforms in smartphones created in unit terms the largest installed base ever seen in the history of computers and is now being further extended by tablets. An underlying additional element of this transition is the increasing role of open source technologies both in software and hardware. This has driven the microprocessor based personal computer industry with few dominating closed operating system platforms into a steep decline. A significant factor in this process has been the separation of processor architecture and processor chip production and operating systems and application development platforms merger into integrated software platforms with proprietary application markets. Furthermore the pay-by-click marketing has changed the way applications development is compensated: Three essays on major trends in a slow clockspeed industry: The case of industrial automation 2014 freeware, ad based or licensed - all at a lower price and used by a wider customer base than ever before. Moreover, the concept of software maintenance contract is very remote in the app world. However, as a slow clockspeed industry, industrial automation has remained intact during the disruptions based on SoC and related software platforms in the ICT industries. Industrial automation incumbents continue to supply systems based on vertically integrated systems consisting of proprietary software and proprietary mainly microprocessor based hardware. They enjoy admirable profitability levels on a very narrow customer base due to strong technology-enabled customer lock-in and customers' high risk leverage as their production is dependent on fault-free operation of the industrial automation systems. When will this balance of power be disrupted? The thesis suggests how industrial automation could join the mainstream ICT industry and create an information, communication and automation (ICAT) industry. Lately the Internet of Things (loT) and weightless networks, a new standard leveraging frequency channels earlier occupied by TV broadcasting, have gradually started to change the rigid world of Machine to Machine (M2M) interaction. It is foreseeable that enough momentum will be created that the industrial automation market will in due course face an architecture disruption empowered by these new trends. This thesis examines the current state of industrial automation subject to the competition between the incumbents firstly through a research on cost competitiveness efforts in captive outsourcing of engineering, research and development and secondly researching process re- engineering in the case of complex system global software support. Thirdly we investigate the industry actors', namely customers, incumbents and newcomers, views on the future direction of industrial automation and conclude with our assessments of the possible routes industrial automation could advance taking into account the looming rise of the Internet of Things (loT) and weightless networks. Industrial automation is an industry dominated by a handful of global players each of them focusing on maintaining their own proprietary solutions. The rise of de facto standards like IBM PC, Unix and Linux and SoC leveraged by IBM, Compaq, Dell, HP, ARM, Apple, Google, Samsung and others have created new markets of personal computers, smartphone and tablets and will eventually also impact industrial automation through game changing commoditization and related control point and business model changes. This trend will inevitably continue, but the transition to a commoditized industrial automation will not happen in the near future.

The handaxe and the microscope: individual and social learning in a multidimensional model of adaptation

When individuals learn by trial-and-error, they perform randomly chosen actions and then reinforce those actions that led to a high payoff. However, individuals do not always have to physically perform an action in order to evaluate its consequences. Rather, they may be able to mentally simulate actions and their consequences without actually performing them. Such fictitious learners can select actions with high payoffs without making long chains of trial-and-error learning. Here, we analyze the evolution of an n-dimensional cultural trait (or artifact) by learning, in a payoff landscape with a single optimum. We derive the stochastic learning dynamics of the distance to the optimum in trait space when choice between alternative artifacts follows the standard logit choice rule. We show that for both trial-and-error and fictitious learners, the learning dynamics stabilize at an approximate distance of root n/(2 lambda(e)) away from the optimum, where lambda(e) is an effective learning performance parameter depending on the learning rule under scrutiny. Individual learners are thus unlikely to reach the optimum when traits are complex (n large), and so face a barrier to further improvement of the artifact. We show, however, that this barrier can be significantly reduced in a large population of learners performing payoff-biased social learning, in which case lambda(e) becomes proportional to population size. Overall, our results illustrate the effects of errors in learning, levels of cognition, and population size for the evolution of complex cultural traits. (C) 2013 Elsevier Inc. All rights reserved.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Quantify this! Report on a round table discussion on quantitative mass spectrometry in proteomics.

Following the success of the first round table in 2001, the Swiss Proteomic Society has organized two additional specific events during its last two meetings: a proteomic application exercise in 2002 and a round table in 2003. Such events have as their main objective to bring together, around a challenging topic in mass spectrometry, two groups of specialists, those who develop and commercialize mass spectrometry equipment and software, and expert MS users for peptidomics and proteomics studies. The first round table (Geneva, 2001) entitled "Challenges in Mass Spectrometry" was supported by brief oral presentations that stressed critical questions in the field of MS development or applications (Stöcklin and Binz, Proteomics 2002, 2, 825-827). Topics such as (i) direct analysis of complex biological samples, (ii) status and perspectives for MS investigations of noncovalent peptide-ligant interactions; (iii) is it more appropriate to have complementary instruments rather than a universal equipment, (iv) standardization and improvement of the MS signals for protein identification, (v) what would be the new generation of equipment and finally (vi) how to keep hardware and software adapted to MS up-to-date and accessible to all. For the SPS'02 meeting (Lausanne, 2002), a full session alternative event "Proteomic Application Exercise" was proposed. Two different samples were prepared and sent to the different participants: 100 micro g of snake venom (a complex mixture of peptides and proteins) and 10-20 micro g of almost pure recombinant polypeptide derived from the shrimp Penaeus vannamei carrying an heterogeneous post-translational modification (PTM). Among the 15 participants that received the samples blind, eight returned results and most of them were asked to present their results emphasizing the strategy, the manpower and the instrumentation used during the congress (Binz et. al., Proteomics 2003, 3, 1562-1566). It appeared that for the snake venom extract, the quality of the results was not particularly dependant on the strategy used, as all approaches allowed Lication of identification of a certain number of protein families. The genus of the snake was identified in most cases, but the species was ambiguous. Surprisingly, the precise identification of the recombinant almost pure polypeptides appeared to be much more complicated than expected as only one group reported the full sequence. Finally the SPS'03 meeting reported here included a round table on the difficult and challenging task of "Quantification by Mass Spectrometry", a discussion sustained by four selected oral presentations on the use of stable isotopes, electrospray ionization versus matrix-assisted laser desorption/ionization approaches to quantify peptides and proteins in biological fluids, the handling of differential two-dimensional liquid chromatography tandem mass spectrometry data resulting from high throughput experiments, and the quantitative analysis of PTMs. During these three events at the SPS meetings, the impressive quality and quantity of exchanges between the developers and providers of mass spectrometry equipment and software, expert users and the audience, were a key element for the success of these fruitful events and will have definitively paved the way for future round tables and challenging exercises at SPS meetings.

Anévrismes complexes de l'aorte abdominale: la chirurgie ouverte est-elle obsolète [The complex aortic abdominal aneurysm: is open surgery old fashion?].

Open surgery is still the main treatment of complex abdominal aortic aneurysm. Nevertheless, this approach is associated with major complications and high mortality rate. Therefore the fenestrated endograft has been used to treat the juxtarenal aneurysms. Unfortunately, no randomised controlled study is available to assess the efficacy of such devices. Moreover, the costs are still prohibitive to generalise this approach. Alternative treatments such as chimney or sandwich technique are being evaluated in order to avoid theses disadvantages. The aim of this paper is to present the endovascular approach to treat juxtarenal aneurysm and to emphasize that this option should be used only by highly specialized vascular centres.

The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.

The geometric structure of the brain fiber pathways.

The structure of the brain as a product of morphogenesis is difficult to reconcile with the observed complexity of cerebral connectivity. We therefore analyzed relationships of adjacency and crossing between cerebral fiber pathways in four nonhuman primate species and in humans by using diffusion magnetic resonance imaging. The cerebral fiber pathways formed a rectilinear three-dimensional grid continuous with the three principal axes of development. Cortico-cortical pathways formed parallel sheets of interwoven paths in the longitudinal and medio-lateral axes, in which major pathways were local condensations. Cross-species homology was strong and showed emergence of complex gyral connectivity by continuous elaboration of this grid structure. This architecture naturally supports functional spatio-temporal coherence, developmental path-finding, and incremental rewiring with correlated adaptation of structure and function in cerebral plasticity and evolution.

Structural macro factors and the affine term structure of interest rates

This work consists of three essays investigating the ability of structural macroeconomic models to price zero coupon U.S. government bonds. 1. A small scale 3 factor DSGE model implying constant term premium is able to provide reasonable a fit for the term structure only at the expense of the persistence parameters of the structural shocks. The test of the structural model against one that has constant but unrestricted prices of risk parameters shows that the exogenous prices of risk-model is only weakly preferred. We provide an MLE based variance-covariance matrix of the Metropolis Proposal Density that improves convergence speeds in MCMC chains. 2. Affine in observable macro-variables, prices of risk specification is excessively flexible and provides term-structure fit without significantly altering the structural parameters. The exogenous component of the SDF is separating the macro part of the model from the term structure and the good term structure fit has as a driving force an extremely volatile SDF and an implied average short rate that is inexplicable. We conclude that the no arbitrage restrictions do not suffice to temper the SDF, thus there is need for more restrictions. We introduce a penalty-function methodology that proves useful in showing that affine prices of risk specifications are able to reconcile stable macro-dynamics with good term structure fit and a plausible SDF. 3. The level factor is reproduced most importantly by the preference shock to which it is strongly and positively related but technology and monetary shocks, with negative loadings, are also contributing to its replication. The slope factor is only related to the monetary policy shocks and it is poorly explained. We find that there are gains in in- and out-of-sample forecast of consumption and inflation if term structure information is used in a time varying hybrid prices of risk setting. In-sample yield forecast are better in models with non-stationary shocks for the period 1982-1988. After this period, time varying market price of risk models provide better in-sample forecasts. For the period 2005-2008, out of sample forecast of consumption and inflation are better if term structure information is incorporated in the DSGE model but yields are better forecasted by a pure macro DSGE model.

Supergenes and complex phenotypes.

Understanding the molecular underpinnings of evolutionary adaptations is a central focus of modern evolutionary biology. Recent studies have uncovered a panoply of complex phenotypes, including locally adapted ecotypes and cryptic morphs, divergent social behaviours in birds and insects, as well as alternative metabolic pathways in plants and fungi, that are regulated by clusters of tightly linked loci. These 'supergenes' segregate as stable polymorphisms within or between natural populations and influence ecologically relevant traits. Some supergenes may span entire chromosomes, because selection for reduced recombination between a supergene and a nearby locus providing additional benefits can lead to locus expansions with dynamics similar to those known for sex chromosomes. In addition to allowing for the co-segregation of adaptive variation within species, supergenes may facilitate the spread of complex phenotypes across species boundaries. Application of new genomic methods is likely to lead to the discovery of many additional supergenes in a broad range of organisms and reveal similar genetic architectures for convergently evolved phenotypes.

How siblings adjust sib-sib communication and begging signals to each other

Parents allocate food resources to their offspring in proportion to the intensity of begging behaviour. Begging encompasses several activities including vocalizations that should honestly signal need and jostling for the position in the nest where parents predictably deliver food items. Although siblings are known to adjust begging level to each other, the underlying mechanism remains unknown. We examined this issue in experimental two-chick broods of the barn owl, Tyto alba, a species in which siblings communicate vocally with each other in the prolonged absence of parents. The function of sib-sib vocal communication, so-called sibling negotiation, is to resolve conflicts over which individual will have priority of access to the next delivered indivisible food item. We found that when a nestling produced longer negotiation calls and stood closer to the nestbox entrance in the absence of parents, its sibling vocally negotiated at a lower rate. Additionally, when an individual produced more negotiation calls in the absence of parents, its sibling begged less intensely at the parent's return, with begging being the key factor that determined which nestling obtained a food item. We conclude that position in the nest and the duration of negotiation calls produced in the absence of parents influence the rate of producing negotiation calls, which in turn influences the rate at which siblings beg for food from their parents. Adjusting begging behaviour could therefore depend on complex sib-sib interactions taking place in the prolonged absence of parents.

A wearable system to time gate crossing during alpine skiing slalom

INTRODUCTION: In alpine skiing, chronometry analysis is currently the most common tool to assess performance. It is widely used to rank competitors during races, as well as to manage athletes training and to evaluate material. Usually, this measurement is accurately realized using timing cells. Nevertheless, these devices are too complex and expensive to allow chronometry of every gates crossing. On the other side, differential GPS can be used for measuring gate crossing time (Waegli et al). However, this is complex (e.g. recording gate position with GPS) and mainly used in research applications. The aim of the study was to propose a wearable system to time gates crossing during alpine skiing slalom (SL), which is suitable for routine uses. METHODS: The proposed system was composed of a 3D accelerometer (ADXL320®, Analog Device, USA) placed at the sacrum of the athlete, a matrix of force sensors (Flexiforce®, Tekscan, USA) fixed on the right shin guard and a data logger (Physilog®, BioAGM, Switzerland). The sensors were sampled at 500 Hz. The crossing time were calculated in two phases. First, the accelerometer was used to detect the curves by considering the maximum of the mediolateral peak acceleration. Then, the force sensors were used to detect the impacts with the gates by considering maximum force variation. In case of non impact, the detection was realized based on the acceleration and features measured at the other gates. In order to assess the efficiency of the system, two different SL were monitored twice for two world cup level skiers, a male SL expert and a female downhill expert. RESULTS AND DISCUSSION: The combination of the accelerometer and force sensors allowed to clearly identify the gate crossing times. When comparing the runs of the SL expert and the downhill expert, we noticed that the SL expert was faster. For example for the first SL, the overall difference between the best run of each athlete was of 5.47s. At each gate, the SL expert increased the time difference slower at the beginning (0.27s/gate) than at the end (0.34s/gate). Furthermore, when comparing the runs of the SL expert, a maximum time difference of 20ms at each gate was noticed. This showed high repeatability skills of the SL expert. In opposite, the downhill expert with a maximum difference time of 1s at each gate was clearly less repeatable. Both skiers were not disturbed by the system. CONCLUSION: This study proposed a new wearable system to automatically time gates crossing during alpine skiing slalom combining force and accelerometer sensors. The system was evaluated with two professional world cup skiers and showed a high potential. This system could be extended to time other parameters. REFERENCES Waegli A, Skaloud J (2007). Inside GNSS, Spring, 24-34.

Local optima networks of hard combinatorial landscapes

Combinatorial optimization involves finding an optimal solution in a finite set of options; many everyday life problems are of this kind. However, the number of options grows exponentially with the size of the problem, such that an exhaustive search for the best solution is practically infeasible beyond a certain problem size. When efficient algorithms are not available, a practical approach to obtain an approximate solution to the problem at hand, is to start with an educated guess and gradually refine it until we have a good-enough solution. Roughly speaking, this is how local search heuristics work. These stochastic algorithms navigate the problem search space by iteratively turning the current solution into new candidate solutions, guiding the search towards better solutions. The search performance, therefore, depends on structural aspects of the search space, which in turn depend on the move operator being used to modify solutions. A common way to characterize the search space of a problem is through the study of its fitness landscape, a mathematical object comprising the space of all possible solutions, their value with respect to the optimization objective, and a relationship of neighborhood defined by the move operator. The landscape metaphor is used to explain the search dynamics as a sort of potential function. The concept is indeed similar to that of potential energy surfaces in physical chemistry. Borrowing ideas from that field, we propose to extend to combinatorial landscapes the notion of the inherent network formed by energy minima in energy landscapes. In our case, energy minima are the local optima of the combinatorial problem, and we explore several definitions for the network edges. At first, we perform an exhaustive sampling of local optima basins of attraction, and define weighted transitions between basins by accounting for all the possible ways of crossing the basins frontier via one random move. Then, we reduce the computational burden by only counting the chances of escaping a given basin via random kick moves that start at the local optimum. Finally, we approximate network edges from the search trajectory of simple search heuristics, mining the frequency and inter-arrival time with which the heuristic visits local optima. Through these methodologies, we build a weighted directed graph that provides a synthetic view of the whole landscape, and that we can characterize using the tools of complex networks science. We argue that the network characterization can advance our understanding of the structural and dynamical properties of hard combinatorial landscapes. We apply our approach to prototypical problems such as the Quadratic Assignment Problem, the NK model of rugged landscapes, and the Permutation Flow-shop Scheduling Problem. We show that some network metrics can differentiate problem classes, correlate with problem non-linearity, and predict problem hardness as measured from the performances of trajectory-based local search heuristics.

Growth mechanisms and timing of emplacement of a vertically layered mafic intrusion in the feeder zone of an ocean island volcano (Fuerteventura, Canary islands)

Résumé pour le grand public L'île de Fuerteventura (Canaries) offre l'occasion rare d'observer les racines d'un volcan océanique édifié il y a 25 à 30 millions d'années et complètement érodé. On y voit de nombreux petits plutons de forme et composition variées, témoignant d'autant d'épisodes de l'activité magmatique. L'un de ces plutons, appelé PX1, présente une structure inhabituelle formée d'une alternance de bandes verticales d'épaisseur métrique à hectométrique de roches sombres de composition pyroxénilique ou gabbroïque. Les pyroxénites résultent clairement de l'accumulation de cristaux de pyroxènes et non de la simple solidification d'un magma? Se pose dès lors la question de la nature du processus qui a conduit à l'accumulation verticale de niveaux concentrés en pyroxènes. En effet, les litages pyroxénitiques classiques sont subhorizontaux, car ils résultent de l'accumulation gravitaire des cristaux séparés du magma dont ils cristalli¬sent par sédimentation. Cette étude vise à identifier et comprendre les mécanismes qui ont engendré ce Iitage minéralogique vertical et l'im¬portant volume de ces faciès cumulatifs. Nous nous sommes également intéressés aux conditions de pression et de température régnant au moment de la mise en place du pluton, ainsi qu'à sa durée de vie et à sa vitesse de refroidis¬sement. Enfin une approche géochimique nous a permis de préciser la nature de la source mantellique des magmas liés à cette activité magmatique. PX1 est en réalité un complexe filonien formé à des conditions de pression et de température de 1-2 kbar et 1050- 1100°C; sa construction a nécessité au moins 150 km3 de magma. L'alternance d'horizons gabbroïques et pyroxéniti¬ques représente des injections successives de magma sous la forme de filons verticaux, mis en place dans un contexte régional en extension. L'étude des orientations des minéraux dans ces faciès révèle que les horizons gabbroïques enregistrent l'extension régionale, alors que les pyroxénites sont générées par une compaction au sein du pluton. Ceci suggère que le régime des contraintes, qui était extensif lors de l'initiation de la mise en place de PX1, est pério¬diquement devenu compressif au sein même du pluton. Cette compression serait liée à des cycles de mise en place où la vitesse de croissance du pluton dépassait celle de l'extension régionale. La différenciation observée au sein de chaque horizon, depuis des pyroxénites riches en olivine jusqu'à des pyroxé¬nites à plagioclase interstitiel et des gabbros, ainsi que la composition géochimique des minéraux qui les constituent suggèrent que chaque filon vertical s'est mis en place à partir d'un magma de composition identique, puis a évolué indépendamment des autres en fonction du régime thermique et du régime des contraintes local. Lorsque le magma en train de cristalliser s'est trouvé en compression, le liquide résiduel a été séparé des cristaux déjà formés et extrait du système, laissant derrière lui une accumulation de cristaux dont la nature et les proportions dépendaient du stade de cristallisation atteint par le magma au moment de l'extraction. Ainsi, les niveaux de pyroxénites à olivine (premier minéral à cristalliser) ont été formés lorsque le magma correspondant était encore peu cristallisé; à l'inverse, les py¬roxénites riches en plagioclase (minéral plus tardif dans la séquence de cristallisation) et certains gabbros à caractère cumulatif résultent d'une compression tardive dans le processus de cristallisation du filon concerné. Les liquides résiduels extraits des niveaux pyroxénitiques sont rarement observés dans PX1, certaines poches et filonets de com¬position anorthositique pourraient en être les témoins. L'essentiel de ces liquides a probablement gagné des niveaux supérieurs du pluton, voire la surface du volcan. L'origine du régime compressif périodique affectant les filons en voie de cristallisation est attribuée aux injections suivantes de magma au sein du pluton, qui se sont succédées à un rythme plus rapide que la vitesse de consolidation des filons. Des datations U/Pb de haute précision sur des cristaux de zircon et de baddeleyite ainsi que40Ar/39Ar sur des cris¬taux d'amphibole révèlent une initiation de la mise en place de PX1 il y a 22.1 ± 0,7 Ma; celle-ci a duré quelque 0,48 ± 0,22 à 0,52 ± 0,29 Ma. Ce laps de temps est compatible avec celui nécessaire à la cristallisation des filons individuels, qui va de moins d'une année lors de l'initiation du magmatisme à 5 ans lors du maximum d'activité de PX1. La présence de cristaux résorbés enregistrant une cristallisation complexe suggère l'existence d'une chambre mag¬matique convective sous-jacente à PX1 et périodiquement rechargée. Les compositions isotopiques des roches étu¬diées révèlent une source mantellique profonde de type point chaud avec une contribution du manteau lithosphéri- que métasomatisé présent sous les îles Canaries. Résumé L'intrusion mafique Miocène PX1 fait partie du soubassement superficiel (0.15-0.2 GPa, 1100 °Q d'un volcan d'île océanique. La particularité de ce pluton est l'existence d'alternances d'unités de gabbros et de pyroxénites qui met¬tent en évidence un litage magmatique vertical (NNE-SSW). Les horizons gabbroiques et pyroxénitiques sont constitués d'unités de différenciation métriques qui suggèrent tine mise en place par injections périodiques de filons verticaux de magma formant un complexe filonien. Chaque filon vertical a subi une différenciation parallèle à un front de solidification sub-vertical parallèle aux bords du filon. Les pyroxénites résultent du fractionnement et de l'accumulation d'olivine ± clinopyroxene ± plagioclase à partir d'un magma basaltique faiblement alcalin et sont interprétées comme étant des imités de différenciation tronquées dont le liquide interstitiel a été extrait par compaction. L'orientation préférentielle des clinopyroxènes dans ces pyroxe- nites (obtenues par analyse EBSD et micro-tomographique) révèle une composante de cisaillement simple dans la genèse de ces roches, ce qui confirme cette interprétation. La compaction des pyroxénites est probablement causée par a mise en place de filons de magma suivants. Le liquide interstitiel expulsé est probablement par ces derniers. Les clinopyroxènes des gabbros, montrent une composante de cisaillement pure suggérant qu'ils sont affectés par une déformation syn-magmatique parallèle aux zones de cisaillement NNE-SSW observées autour de PX1 et liées au contexte tectonique Miocène d'extension régionale. Ceci suggère que les gabbros sont liés à des taux de mise en place faibles à la fin de cycles d'activité magmatique et sont peu ou pas affectés par la compaction. L'initiation et la géométrie de PX1 sont donc contrôlées par le contexte tectonique régional d'extension alors que les taux et les volumes de magma dépendent de facteurs liés à la source. Des taux d'injection élevés résultent probable¬ment en une croissance du pluton supérieure à la place crée par cette extension. Dans ce cas de figure, la propagation des nouveaux dykes et l'inaptitude du magma à circuler à travers les anciens dykes cristallisés pourrait causer une augmentation de la pression non-lithostatique sur ces derniers, exprimée par un cisaillement simple et l'expulsion du liquide interstitiel qu'ils contiennent (documenté par les zones de collecte anorthositiques). Les compositions en éléments majeurs et traces des gabbros et pyroxenites de PX1 sont globalement homogènes et dépendent de la nature cumulative des échantillons. Cependant, de petites variations des concentrations en éléments traces ainsi que les teneurs en éléments traces des bordures de clinopyroxenes suggèrent que ces derniers ont subi un processus de rééquilibrage et de cristallisation in situ. L'homogénéité des compositions chimiques des échantillons, ainsi que la présence de grains de clinopyroxene résorbés suggère que le complexe filonien PX1 s'est mis en place au dessus d'une chambre magmatique périodiquement rechargée dans laquelle la convection est efficace. Chaque filon est donc issu d'un même magma, mais a subi une différenciation par cristallisation in situ (jusqu'à 70% de fraction¬nement) indépendamment des autres. Dans ces filons cristallisés, les minéraux cumulatifs subissent un rééquilibrage partiel avec les liquide interstitiel avant que ce dernier ne soit expulsé lors de la compaction (mettant ainsi un terme à la différenciation). Ce modèle de mise en place signifie qu'un minimum de 150Km3 de magma est nécessaire à la genèse de PX1, une partie de ce volume ayant été émis par le 'Central Volcanic Complex' de Fuerteventura. Les rapports isotopiques radiogéniques mesurés révèlent la contribution de trois pôles mantelliques dans la genèse du magma formant PX1. Le mélange de ces pôles HIMU, DMM et EM1 refléterai l'interaction du point chaud Cana¬rien avec un manteau lithosphérique hétérogène métasomatisé. Les petites variations de ces rapports et des teneurs en éléments traces au sein des faciès pourrait refléter des taux de fusion partielle variable de la source, résultant en un échantillonnage variable du manteau lithosphérique métasomatisé lors de son interaction avec le point chaud. Des datations U/Pb de haute précision (TIMS) sur des cristaux de zircon et de baddeleyite extraits de gabbros de PX1 révèlent que l'initiation de la cristallisation du magma a eu lieu il y a 22.10±0.07 Ma et que l'activité magmatique a duré un minimum de 0.48 à 0.52 Ma. Des âges 40Ar/39Ar obtenus sur amphibole sont de 21.9 ± 0.6 à 21.8 ± 0.3 Ma, identiques aux âges U/Pb. La combinaison de ces méthodes de datations, suggère que le temps maximum nécessaire à PX1 pour se refroidir en dessous de la température de fermeture de l'amphibole est de 0.8Ma. Ceci signifie que la durée de vie de PX1 est de 520 000 à 800 000 ans. La coexistence de cristaux de baddeleyite et de zircon dans un gabbro est attribuée à son interaction avec un fluide riche en C02 relâché par les carbonatites encaissantes lors du métamorphisme de contact généré par la mise en place de PX1 environ 160 000 ans après le début de sa mise en place. Les durées de vie obtenue sont en accord avec le modèle de mise en place suggérant une durée de cristallisation poux chaque filon allant de 1 an à 5 ans. Abstract The Miocene PX1 gabbro-pyroxenite intrusion (Fuerteventura, Canary Islands), is interpreted as the shallow-level feeder-zone (0.15-0.2 GPa and 1100-1120°C), to an ocean island volcano. The particularity of PX1 is that it displays a NNE-SSW trending vertical magmatic banding expressed by alternating gabbro and pyroxeriite sequences. The gabbro and pyroxenite sequences consist of metre-thick differentiation units, which suggest emplacement by pe¬riodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke com¬plex. Individual dykes underwent internal differentiation following a solidification front (favoured by a significant lateral/horizontal thermal gradient) parallel to the dyke edges. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene ± olivine ± plagioclase crystals from a mildly alkaline basaltic liquid and are interpre¬ted as truncated differentiation sequences, from which residual melts were extracted by compaction. Clinopyroxene mineral orientation in pyroxenites (evidenced by EBSD and micro X-ray tomography analysis) display a marked pure shear component, supporting this interpretation. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The resulting expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Gabbro sequences represent crystallised coalesced magma batches, emplaced at lower rates at the end of eruptive cycles, and underwent minor melt extraction as evi¬denced by clinopyroxene orientations that record a simple shear component suggesting syn-magmatic deformation parallel to observed NNF.-SSW trending shear-zones induced by the regional tensional Miocene stress-field. The initiation and geometry of PX1 is controlled by the regional extensional tectonic regime whereas rates and vo¬lumes of magma depend on source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke tip geometry and the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; generating strong pure-shear compaction and interstitial melt expulsion within the feeder-zone as recorded by the cumulitic pyroxenite bands and anorthositic collection zones. The whole-rock major and trace-element chemistry of PX1 gabbros and pyroxenites is globally homogeneous and controlled by the cumulate nature of the samples (i.e. on the modal proportions of olivine, pyroxene, plagioclase and oxides). However, small variations of whole-rock trace-element contents as well as trace-element contents of clinopyroxene rims suggest that in-situ re-equilibration and crystallisation has occurred. Additionally, the global homogeneity and presence of complex zoning of rare resorbed clinopyroxene crystals suggest that the PX1 feeder- zone overlies a periodically replenished and efficiently mixed magma chamber. Each individual dyke of magma thus originated from a compositionally constant mildly alkaline magma and differentiated independently from the others reaching up to 70% fractionation. Following dyke arrest these are affected by interaction with the trapped interstitial liquid prior to its compaction-linked expulsion (thus stopping the differentiation process). This emplacement model implies that minimum amount of approximately 150 km3 of magma is needed to generate PX1, part of it having been erupted through the overlying Central Volcanic Complex of Fuerteventura. The radiogenic isotope ratios of PX1 samples reveal the contribution on three end-members during magma genesis. This mixing of the H1MU, EMI and DMM end-members could reflect the interaction of the deep-seated Canarian mantle plume with a heterogeneous metasomatic and sepentininsed lithospheric mantle. Additionally, the observed trace-element and isotopic variations within the same fades groups could reflect varying degrees of partial melting of the source region, thus tapping more or less large areas of the metasomatised lithospheric mantle during interac¬tion with the plume. High precision ID-TIMS U/Pb zircon and baddeleyite ages from the PX1 gabbro samples, indicate initiation of magma crystallisation at 22.10 ± 0.07 Ma. The magmatic activity lasted a minimum of 0.48 to 0.52 Ma. 40Ar/39Ar amphibole ages are of 21.9 ± 0.6 to 21.8 ± 0.3, identical within errors to the U/Pb ages. The combination of the 40Ar/39Ar and U/Pb datasets imply that the maximum amount of time PX1 took to cool below amphibole Tc is 0.8 Ma, suggesting PX1 lifetime of 520 000 to 800 000 years. On top of this, the coexistence of baddeleyite and zircon in a single sample is ascribed to the interaction of PX1 with C02-rich carbonatite-derived fluids released from the host-rock carbonatites during contact metamorphism 160 000 years after PX1 initiation. These ages are in agreement with the emplacement model, implying a crystallisation time of less than 1 to 5 years for individual dykes.

Contribution of isotope ratio mass spectrometry to the investigation of improvised explosives: isotopic study of black powders ans ammonium nitrates

The establishment of legislative rules about explosives in the eighties has reduced the illicit use of military and civilian explosives. However, bomb-makers have rapidly taken advantage of substances easily accessible and intended for licit uses to produce their own explosives. This change in strategy has given rise to an increase of improvised explosive charges, which is moreover assisted by the ease of implementation of the recipes, widely available through open sources. While the nature of the explosive charges has evolved, instrumental methods currently used in routine, although more sensitive than before, have a limited power of discrimination and allow mostly the determination of the chemical nature of the substance. Isotope ratio mass spectrometry (IRMS) has been applied to a wide range of forensic materials. Conclusions drawn from the majority of the studies stress its high power of discrimination. Preliminary studies conducted so far on the isotopic analysis of intact explosives (pre-blast) have shown that samples with the same chemical composition and coming from different sources could be differentiated. The measurement of stable isotope ratios appears therefore as a new and remarkable analytical tool for the discrimination or the identification of a substance with a definite source. However, much research is still needed to assess the validity of the results in order to use them either in an operational prospect or in court. Through the isotopic study of black powders and ammonium nitrates, this research aims at evaluating the contribution of isotope ratio mass spectrometry to the investigation of explosives, both from a pre-blast and from a post-blast approach. More specifically, the goal of the research is to provide additional elements necessary to a valid interpretation of the results, when used in explosives investigation. This work includes a fundamental study on the variability of the isotopic profile of black powder and ammonium nitrate in both space and time. On one hand, the inter-variability between manufacturers and, particularly, the intra-variability within a manufacturer has been studied. On the other hand, the stability of the isotopic profile over time has been evaluated through the aging of these substances exposed to different environmental conditions. The second part of this project considers the applicability of this high-precision technology to traces and residues of explosives, taking account of the characteristics specific to the field, including their sampling, a probable isotopic fractionation during the explosion, and the interferences with the matrix of the site.

Geology, correlations, and geodynamic evolution of the Biga Peninsula (NW Turkey)

L?objectif de ce travail de recherche était de décrypter l?évolution géodynamique de la Péninsule de Biga (Turquie du N-O), à travers l?analyse de deux régions géologiques peu connues, le mélange de Çetmi et la zone d?Ezine (i.e. le Groupe d?Ezine et l?ophiolite de Denizgören). Une étude complète et détaillée de terrain (cartographie et échantillonnage) ainsi qu?une approche multidisciplinaire (sédimentologie de faciès, pétrographie sédimentaire et magmatique, micropaléontologie, datations absolues, géochimie sur roche totale, cristallinité de l?illite) ont permis d?obtenir de nouveaux éléments d?information sur la région considérée. ? Le mélange de Çetmi, de type mélange d?accrétion, affleure au nord et au sud de la Péninsule de Biga ; les principaux résultats de son étude peuvent se résumer comme suit: - Son aspect structural actuel (nature des contacts, organisation tectonique) est principalement dû au régime extensif Tertiaire présent dans la région. - Il est constitué de blocs de différentes natures : rares calcaires Scythien-Ladinien dans le faciès Han Bulog, blocs hectométriques de calcaires d?âge Norien-Rhaetien de rampe carbonatée, nombreux blocs décamétriques de radiolarites rouges d?âge Bajocien- Aptien, blocs/écailles de roches magmatiques de type spilites (basaltes à andésite), ayant des signatures géochimiques d?arcs ou intra-plaques. - La matrice du mélange est constituée d?une association greywacke-argilites dont l?âge Albien inférieur à moyen a été déterminé par palynologie. - L?activité du mélange s?est terminée avant le Cénomanien (discordance Cénomanienne au sommet du mélange, pas de bloc plus jeune que la matrice). - Du point de vue de ses corrélations latérales, le mélange de Çetmi partage plus de traits communs avec les mélanges se trouvant dans les nappes allochtones du Rhodope (nord de la Grèce et sud-ouest de la Bulgarie) qu?avec ceux de la suture Izmir-Ankara (Turquie); il apparaît finalement que sa mise en place s?est faite dans une logique balkanique (chevauchements vers le nord d?âge anté-Cénomanien). ? Le Groupe d?Ezine et l?ophiolite sus-jacente de Denizgören affleurent dans la partie ouest de la Péninsule de Biga. Le Groupe d?Ezine est une épaisse séquence sédimentaire continue (3000 m), subdivisée en trois formations, caractérisée chacune par un type de sédimentation spécifique, relatif à un environnement de dépôt particulier. De par ses caractéristiques (grande épaisseur, variations latérales de faciès et d?épaisseur dans les formations, érosion de matériel provenant de l?amont du bassin), le groupe d?Ezine est interprétée comme un dépôt syn-rift d?âge Permien moyen-Trias inférieur. Il pourrait représenter une partie de la future marge passive sud Rhodopienne à la suite de l?ouverture de l?océan Maliac/Méliata. L?ophiolite de Denizgören sus-jacente repose sur le Groupe d?Ezine par l?intermédiaire d?une semelle métamorphique à gradient inverse, du faciès amphibolite à schiste vert. L?âge du faciès amphibolite suggère une initiation de l?obduction au Barrémien (125 Ma, âge Ar/Ar); cet âge est unique dans le domaine égéen, mais il peut là aussi être relié à une logique balkanique, sur la base de comparaison avec le domaine Rhodopien. ? Toutes les unités précédentes (mélange de Çetmi, Groupe d?Ezine et ophiolite de Denizgören) ont passivement subi trois phases extensives pendant le Tertiaire. Dans la région d?Ezine et du mélange nord, les micaschistes HP sous-jacents ont été exhumés avant l?Eocène moyen. Dans le cas du mélange sud, cette exhumation Eocene est en partie enregistrée dans les mylonites séparant le mélange du dôme métamorphique sous-jacent du Kazda?. Le mélange sud est dans tous les cas fortement érodé à la suite de la double surrection du dôme du Kazda?, près de la lim ite Oligocène/Miocene et pendant le Plio- Quaternaire. Dans le premier cas, ce soulèvement est caractérisé par le développement d?une faille de détachement à faible pendage, qui contrôle à la fois l?exhumation du massif, et la formation d?un bassin sédimentaire syntectonique, de type bassin supradétachement; quant à la phase extensive la plus récente, elle est contrôlée par le jeu de failles normales à forts pendages qui remanient l?ensemble des structures héritées, et dictent la géomorphologie actuelle de la région. ? Il est possible de proposer un scénario pour l?évolution géodynamique de la Péninsule de Biga, basé sur l?ensemble des résultats précédents et sur les données de la géologie régionale ; ses points principaux sont: - La Péninsule de Biga fait partie de la marge Rhodopienne. - Le Groupe d?Ezine est un témoin de la marge passive nord Maliac/Méliata. - L?ophiolite de Denizgören et le mélange de Çetmi ont été mis en place tous deux vers le nord sur la marge précédente, respectivement au Barrémien et à l?Albien terminal- Cénomanien inférieur. - Une forte composante décrochante durant l?emplacement est suggérée par la préservation de fragments de la marge passive et l?absence de métamorphisme dans la plaque inférieure. - Tous les évènements précédents ont été largement affectés par le régime d?extension Tertiaire.<br/><br/>The purpose of this study is to unravel the geodynamic evolution of the Biga Peninsula (NW Turkey) through the detailed study of two poorly known areas, the Çetmi mélange and the Ezine zone (i.e. the Ezine Group and the Denizgören ophiolite). The methodology was based on a detailed field work and a multidisciplinary approach. ? The accretion-related Çetmi mélange is mainly cropping out north and south of the Biga Peninsula; the main results of its study can be summarized as follows: -Its present-day structural aspect (type of contacts, tectonic organisation) is largely inherited from the Tertiary extensional regime in the region. -It is made of blocks of various natures: Han Bulog limestones with a Scythian to Ladinian age, common carbonate ramp Norian-Rhaetian limestones (biggest blocks of the mélange), red radolarite with a Bajocian to Aptian age; the most common lithology of the mélange is made by block/slices of spilitic magmatic rocks (basalt to andesite); they have volcanic arc or within plate basalt geochemical signatures. -The matrix of the mélange is made of a greywacke-shale association of Early-Middle Albian age. - The mélange stopped its activity before the Cenomanian (no younger blocks than the matrix, and Cenomanian unconformity). - If compared to the regional geology, the Çetmi mélange shares some characteristics with the Izmir-Ankara mélanges (less), and with the mélanges from allochthonous nappes found in eastern Rhodope (more); it appears finally that its emplacement is related to a Balkanic logic (ante-Cenomanian northward thrusting). ? The Ezine Group and the overlying Denizgören ophiolite are cropping out in the western part of the Biga Peninsula. The Ezine Group is a thick sedimentary sequence interpreted as a syn-rift deposit of Middle Permian-Early Triassic age. It represents a part of the south Rhodopian passive margin, following the opening of the Maliac/Meliata oceanic domain. The Denizgören ophiolite has been emplaced northward on the Ezine Group in the Barremian (125 Ma, age of the amphibolitic sole); this age is unique in the Aegean domain, but here again, it may be related to a Balkan logic. ? All the previous units (Çetmi mélange, Ezine Group and Denizgören ophiolite) have passively suffered two extensional regimes during the Tertiary. In the Ezine and northern Çetmi mélange area, the underlying HP Çamlýca micaschists were exhumed before the Middle Eocene. As for the southern mélange, it was strongly eroded following the Late Oligocene to Quaternary uplift of the underlying Kazda? Massif. This uplift was characterized by the development of a low-angle detachment fault controlling a part of the exhumation, as well as the development of a supra-detachment basin. ? Based on the previous results, and on the data from the regional geology, one can propose a scenario for the geodynamic evolution of the Biga Peninsula. Its key points are:- The Biga Peninsula is belonging to the Rhodope margin. - The Ezine Group is a remnant of the northern Maliac/Meliata passive margin. - Both the Denizgören ophiolite and the Çetmi mélange have been emplaced northward on the previous margin, respectively in the Barremian and in the Late Albian-Early Cenomanian times. - The preservation of the remnants of the Rhodope margin, as well as the absence of metamorphism in the lower plate suggest a strong strike-slip component during the emplacements. - All the previous events are (at least) partly obliterated by the Tertiary extensional regime.<br/><br/>Le géologue est comme un «historien» de la Terre, qui porte un intérêt particulier à l?étude du passé de notre planète; ce dernier, très ancien, se mesure en dizaines ou centaines de millions d?années (Ma). Or le visage de la terre a constamment évolué au cours des ces millions d?années écoulés, car les plaques (continentales et océaniques) qui composent son enveloppe superficielle ne restent pas immobiles, mais se déplacent continuellement à sa surface, à une vitesse de l?ordre du cm/an (théorie de la tectonique des plaques); c?est ainsi, par exemple, que des océans naissent, grandissent, puis finissent par se refermer. On appelle sutures océaniques, les zones, aujourd?hui sur la terre ferme, où l?on retrouve les restes d?océans disparus. Ces sutures sont caractérisées par deux associations distinctes de roches, que l?on appelle les mélanges et les ophiolites; ces mélanges et ophiolites sont donc les témoins de l?activité passée d?un océan aujourd?hui refermé. L?équipe de recherche dans laquelle ce travail à été réalisé s?intéresse à un vaste domaine océanique fossile: l?océan Néotéthys. Cet océan, de plusieurs milliers de kilomètres de large, séparait alors l?Europe et l?Asie au nord, de l?Afrique, l?Inde et l?Australie au sud. De cet océan, il n?en subsiste aujourd?hui qu?une infime partie, qui se confond avec notre mer Méditerranée actuelle. Or, tout comme l?océan Pacifique est bordé de mers plus étroites (Mer de Chine, du Japon, etc?), l?océan Néotéthys était bordé au nord de mers marginales. C?est dans ce cadre que s?est inscrit mon travail de thèse, puisqu?il a consisté en l?étude d?une suture océanique (mélange plus ophiolite), témoin d?une des mers qui bordait l?océan Néotéthys sur sa marge nord. L?objectif était de préciser de quelle suture il s?agissait, puis de déterminer quand et comment elle avait fonctionné (i.e son évolution géologique). Les roches qui composent cette suture affleurent aujourd?hui en Turquie nord occidentale dans la Péninsule de Biga. Au nord et au sud de la péninsule se trouvent les zones géologique du mélange de Çetmi, et à l?ouest, le Groupe d?Ezine et l?ophiolite susjacente, dite ophiolite de Denizgören. Une étude complète et détaillée de terrain (cartographie, échantillonnage), suivie de diverses analyses en laboratoire (détermination de leur âge, de leur condition de formation, etc?), ont permis d?aboutir aux principaux résultats suivants : - Mise en évidence dans le mélange de Çetmi des témoins (1) de l?océan Lycien disparu (ancienne mer marginale de la Néotéthys), et (2) de la marge continentale qui le bordait au nord. - Fin de l?activité du mélange de Çetmi il y a environ 105 Ma (Albien). - Le mélange de Çetmi est difficilement corrélable dans le temps avec les unités semblables affleurant dans la région d?étude (unicité du mélange), ce qui implique des conditions particulière de formation. - L?ophiolite de Denizgören est un morceau d?océan Lycien posé sur un reste préservé de sa marge continentale nord. - Cette dernière est représentée sur le terrain par une succession de roches caractéristiques, le Groupe d?Ezine. Celui-ci est lui-même un témoin de l?ouverture d?un océan marginal de la Néotethys antérieur au Lycien, l?océan Maliac, qui s?est ouvert il y a 245 Ma (Permien-Trias). - La mise en place de l?ophiolite de Denizgören sur le Groupe d?Ezine (125 Ma, Barrémien) est antérieure à la mise en place du mélange de Çetmi. - Il apparaît que ces deux mises en place sont contemporaines de la formation de la chaîne des Balkans, terminée avant le Cénomanien (100 Ma). - L?évolution dans le temps des objets précédents (océans, marges continentales) montre de grands mouvements latéraux est-ouest entre ces objets (translation). Ce qui implique que les roches que l?on retrouve aujourd?hui sur un transect nord-sud ne l?étaient pas nécessairement auparavant. - Enfin, il s?avère que le mélange de Çetmi, l?ophiolite de Denizgören, et le Groupe d?Ezine ont subi par la suite des déformations extensives importantes qui ont considérablement perturbé le schéma post-mise en place.