Efficient total variation algorithm for fetal brain MRI reconstruction.

Fetal MRI reconstruction aims at finding a high-resolution image given a small set of low-resolution images. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has considered several regularization terms s.a. Dirichlet/Laplacian energy, Total Variation (TV)- based energies and more recently non-local means. Although TV energies are quite attractive because of their ability in edge preservation, standard explicit steepest gradient techniques have been applied to optimize fetal-based TV energies. The main contribution of this work lies in the introduction of a well-posed TV algorithm from the point of view of convex optimization. Specifically, our proposed TV optimization algorithm or fetal reconstruction is optimal w.r.t. the asymptotic and iterative convergence speeds O(1/n2) and O(1/√ε), while existing techniques are in O(1/n2) and O(1/√ε). We apply our algorithm to (1) clinical newborn data, considered as ground truth, and (2) clinical fetal acquisitions. Our algorithm compares favorably with the literature in terms of speed and accuracy.

Kernel-based Mapping of Meteorological Fields in Complex Orography

Due to the advances in sensor networks and remote sensing technologies, the acquisition and storage rates of meteorological and climatological data increases every day and ask for novel and efficient processing algorithms. A fundamental problem of data analysis and modeling is the spatial prediction of meteorological variables in complex orography, which serves among others to extended climatological analyses, for the assimilation of data into numerical weather prediction models, for preparing inputs to hydrological models and for real time monitoring and short-term forecasting of weather.In this thesis, a new framework for spatial estimation is proposed by taking advantage of a class of algorithms emerging from the statistical learning theory. Nonparametric kernel-based methods for nonlinear data classification, regression and target detection, known as support vector machines (SVM), are adapted for mapping of meteorological variables in complex orography.With the advent of high resolution digital elevation models, the field of spatial prediction met new horizons. In fact, by exploiting image processing tools along with physical heuristics, an incredible number of terrain features which account for the topographic conditions at multiple spatial scales can be extracted. Such features are highly relevant for the mapping of meteorological variables because they control a considerable part of the spatial variability of meteorological fields in the complex Alpine orography. For instance, patterns of orographic rainfall, wind speed and cold air pools are known to be correlated with particular terrain forms, e.g. convex/concave surfaces and upwind sides of mountain slopes.Kernel-based methods are employed to learn the nonlinear statistical dependence which links the multidimensional space of geographical and topographic explanatory variables to the variable of interest, that is the wind speed as measured at the weather stations or the occurrence of orographic rainfall patterns as extracted from sequences of radar images. Compared to low dimensional models integrating only the geographical coordinates, the proposed framework opens a way to regionalize meteorological variables which are multidimensional in nature and rarely show spatial auto-correlation in the original space making the use of classical geostatistics tangled.The challenges which are explored during the thesis are manifolds. First, the complexity of models is optimized to impose appropriate smoothness properties and reduce the impact of noisy measurements. Secondly, a multiple kernel extension of SVM is considered to select the multiscale features which explain most of the spatial variability of wind speed. Then, SVM target detection methods are implemented to describe the orographic conditions which cause persistent and stationary rainfall patterns. Finally, the optimal splitting of the data is studied to estimate realistic performances and confidence intervals characterizing the uncertainty of predictions.The resulting maps of average wind speeds find applications within renewable resources assessment and opens a route to decrease the temporal scale of analysis to meet hydrological requirements. Furthermore, the maps depicting the susceptibility to orographic rainfall enhancement can be used to improve current radar-based quantitative precipitation estimation and forecasting systems and to generate stochastic ensembles of precipitation fields conditioned upon the orography.

Enhanced compressed sensing recovery with level set normals.

We propose a compressive sensing algorithm that exploits geometric properties of images to recover images of high quality from few measurements. The image reconstruction is done by iterating the two following steps: 1) estimation of normal vectors of the image level curves, and 2) reconstruction of an image fitting the normal vectors, the compressed sensing measurements, and the sparsity constraint. The proposed technique can naturally extend to nonlocal operators and graphs to exploit the repetitive nature of textured images to recover fine detail structures. In both cases, the problem is reduced to a series of convex minimization problems that can be efficiently solved with a combination of variable splitting and augmented Lagrangian methods, leading to fast and easy-to-code algorithms. Extended experiments show a clear improvement over related state-of-the-art algorithms in the quality of the reconstructed images and the robustness of the proposed method to noise, different kind of images, and reduced measurements.

The birth of the Rheic Ocean - Early Palaeozoic subsidence patterns and subsequent tectonic plate scenarios

New plate-tectonic reconstructions of the Gondwana margin suggest that the location of Gondwana-derived terranes should not only be guided by the models, but should also consider the possible detrital input from some Asian blocks (Hunia), supposed to have been located along the Cambrian Gondwana margin, and accreted in the Silurian to the North-Chinese block. Consequently, the Gondwana margin has to be subdivided into a more western domain, where the future Avalonian blocks will be separated from Gondwana by the opening Rheic Ocean, whereas in its eastern continuation, hosting the future basement areas of Central Europe, different periods of crustal extension should be distinguished. Instead of applying a rather cylindrical model, it is supposed that crustal extension follows a much more complex pattern, where local back-arcs or intra-continental rifts are involved. Guided by the age data of magmatic rocks and the pattern of subsidence curves, the following extensional events can be distinguished: During the early to middle Cambrian, a back-arc setting guided the evolution at the Gondwana margin. Contemporaneous intra-continental rift basins developed at other places related to a general post-PanAfrican extensional phase affecting Africa Upper Cambrian formation of oceanic crust is manifested in the Chamrousse area, and may have lateral cryptic relics preserved in other places. This is regarded as the oceanisation of some marginal basins in a context of back-arc rifting. These basins were closed in a mid-Ordovician tectonic phase, related to the subduction of buoyant material (mid-ocean ridge?) Since the Early Ordovician, a new phase of extension is observed, accompanied by a large-scale volcanic activity, erosion of the rift shoulders generated detritus (Armorican Quartzite) and the rift basins collected detrital zircons from a wide hinterland. This phase heralded the opening of Palaeotethys, but it failed due to the Silurian collision (Eo-Variscan phase) of an intra-oceanic arc with the Gondwana margin. During this time period, at the eastern wing of the Gondwana margin begins the drift of the future Hunia microcontinents, through the opening of an eastern prolongation of the already existing Rheic Ocean. The passive margin of the remaining Gondwana was composed of the Galatian superterranes, constituents of the future Variscan basement areas. Remaining under the influence of crustal extension, they will start their drift to Laurussia since the earliest Devonian during the opening of the Palaeotethys Ocean. (C) 2008 Elsevier B.V. All rights reserved.

Late Cretaceous to Eocene geology of the South Central American forearc area (southern Costa Rica and western Panama) : initiation and evolution of an intra-oceanic convergent margin

Résumé : L'arc volcanique du sud de l'Amérique Centrale se situe sur la marge SW de la Plaque Caraïbe, au-dessus des plaques subduites de Cocos et Nazca. Il s'agit de l'un des arcs intra-océaniques les plus étudiés au monde, qui est généralement considéré comme s'étant développé à la fin du Crétacé le long d'un plateau océanique (le Plateau Caraïbe ou CLIP) et se trouvant actuellement dans un régime de subduction érosive. Au cours des dernières décennies, des efforts particuliers ont été faits pour comprendre les processus liés à la subduction sur la base d'études géophysiques et géochimiques. Au sud du Costa Rica et à l'ouest du Panama, des complexes d'accrétions et structures à la base de l'arc volcanique ont été exposés grâce à la subduction de rides asismiques et de failles transformantes. Des affleurements, situés jusqu'à seulement 15 km de la fosse, offrent une possibilité unique de mieux comprendre quelques uns des processus ayant lieu le long de la zone de subduction. Nous présentons de nouvelles contraintes sur l'origine de ces affleurements en alliant une étude de terrain poussée, de nouvelles données géochimiques, sédimentaires et paléontologiques, ainsi que des observations structurales effectuées en télédétection. Une nouvelle stratigraphie tectonique entre le Campanien et l'Éocène est définie pour la région d'avant-arc située entre la Péninsule d'Osa (Costa Rica) et la Péninsule d'Azuero (Panama). Nos résultats montrent que la partie externe de la marge est composée d'un arrangement complexe de roches ignées et de séquences sédimentaires de recouvrement qui comprennent principalement le socle de l'arc, des roches d'arc primitif, des fragments de monts sous-marins accrétés et des mélanges d'accrétion. Des preuves sont données pour le développement de l'arc volcanique du sud de l'Amérique Centrale sur un plateau océanique. Le début de la subduction le long de la marge SW de la Plaque Caraïbe a eu lieu au Campanien et a généré des roches d'arc primitif caractérisées par des affinités géochimiques particulières, globalement intermédiaires entre des affinités de plateau et d'arc insulaire. L'arc était mature au Maastrichtien et formait un isthme essentiellement continu entre l'Amérique du Nord et l'Amérique du Sud. Ceci a permis la migration de faunes terrestres entre les Amériques et pourrait avoir contribué à la crise fin Crétacé -Tertiaire en réduisant les courants océaniques subéquatoriaux entre le Pacifique et l'Atlantique. Plusieurs unités composées de fragments de monts sous-marins accrétés sont définies. La nature et l'arrangement structural de ces unités définissent de nouvelles contraintes sur les modes d'accrétion des monts sous-marins/îles océaniques et sur l'évolution de la marge depuis la formation de la zone de subduction. Entre la fin du Crétacé et l'Éocène moyen, la marge a enregistré plusieurs épisodes ponctuels d'accrétion de monts sous-marins alternant avec de la subduction érosive. A l'Éocène moyen, un événement tectonique régional pourrait avoir causé un fort couplage entre les plaques supérieure et inférieure, menant à des taux plus important d'accrétion de monts sous-marins. Durant cette période, la situation le long de la marge était très semblable à la situation actuelle et caractérisée par la présence de monts sous-marins subductants et l'absence d'accrétion de sédiments. L'enregistrement géologique montre qu'il n'est pas possible d'attribuer une nature érosive ou accrétionnaire à la marge dans le passé ou -par analogie- aujourd'hui, parce que (1) les processus d'accrétion et érosifs varient fortement spatialement et temporellement et (2) il est impossible d'évaluer la quantité exacte de matériel tectoniquement enlevé à la marge depuis le début de la subduction. Au sud du Costa Rica, certains fragments de monts sous-marins accrétés sont représentatifs d'une interaction entre une ride et un point chaud dans le Pacifique au Crétacé terminal/Paléocène. L'existence de ces fragments de monts sous-marins et la morphologie du fond de l'Océan Pacifique indiquent que la formation de la ride de Cocos-Nazca s'est formée au moins ~40 Ma avant l'âge proposé par les modèles tectoniques actuels. Au Panama, nous avons identifié une île océanique d'âge début Éocène qui a été accrétée à l'Éocène moyen. L'accrétion a eu lieu à très faible profondeur par détachement de l'île dans la fosse, et a mené à une exceptionnelle préservation des structures volcaniques. Des affleurement comprenant aussi bien des parties basses et hautes de l'édifice volcanique on été étudiées, depuis la phase sous-marine bouclier jusqu'à la phase subaérienne post-bouclier. La stratigraphie nous a permis de différencier les laves de la phase sous-marine de celles de la phase subaérienne. La composition des laves indique une diminution progressive de l'intensité de la fusion partielle de la source et une diminution de la température des laves produites durant les derniers stades de l'activité volcanique. Nous interprétons ces changements comme étant liés à l'éloignement progressif de l'île océanique de la zone de fusion ou point chaud. Abstract The southern Central American volcanic front lies on the SW edge of the Caribbean Plate, inboard of the subducting Cocos and Nazca Plates. It is one of the most studied intra-oceanic convergent margins around the world, which is generally interpreted to have developed in the late Cretaceous along an oceanic plateau (the Caribbean Large Igneous Province or CLIP) and to be currently undergoing a regime of subduction erosion. In the last decades a particular effort has been made to understand subduction-related processes on the basis of geophysical and geochemical studies. In southern Costa Rica and western Panama accretionary complexes and structures at the base of the volcanic front have been exposed in response to subduction of aseismic ridges and transforms. Onland exposures are located as close as to 15 km from the trench and provide a unique opportunity to better understand some of the processes occurring along the subduction zone. We provide new constraints on the origins of these exposures by integrating a comprehensive field work, new geochemical, sedimentary and paleontological data, as well as structural observations based on remote imaging. A new Campanian to Eocene tectonostratigraphy is defined for the forearc area located between the Osa Peninsula (Costa Rica) and the Azuero Peninsula (Panama). Our results show that the outer margin is composed of a complicated arrangement of igneous complexes and overlapping sedimentary sequences that essentially comprise an arc basement, primitive island-arc rocks, accreted seamount fragments and accretionary mélanges. Evidences are provided for the development of the southern Central American arc on the top an oceanic plateau. The subduction initiation along the SW edge of the Caribbean Plate occurred in the Campanian and led to formation of primitive island-arc rocks characterized by unusual geochemical affinities broadly intermediate between plateau and arc affinities. The arc was mature in the Maastrichtian and was forming a predominantly continuous landbridge between the North and South Americas. This allowed migration of terrestrial fauna between the Americas and may have contributed to the Cretaceous-Tertiary crisis by limiting trans-equatorial oceanic currents between the Pacific and the Atlantic. Several units composed of accreted seamount fragments are defined. The nature of the units and their structural arrangement provide new constraints on the modes of accretion of seamounts/oceanic islands and on the evolution of the margin since subduction initiation. Between the late Cretaceous and the middle Eocene, the margin recorded several local episodes of seamount accretion alternating with tectonic erosion. In the middle Eocene a regional tectonic event may have triggered strong coupling between the overriding and subducting plates, leading to higher rates of seamount accretion. During this period the situation along the margin was very similar to the present and characterized by subducting seamounts and absence of sediment accretion. The geological record shows that it is not possible to ascribe an overall erosive or accretionary nature to the margin in the past and, by analogy, today, because (1) accretionary and erosive processes exhibit significant lateral and temporal variations and (2) it is impossible to estimate the exact amount of material tectonically eroded from the margin since subduction initiation. In southern Costa Rica, accreted seamount fragments point toward a plume-ridge interaction in the Pacific in the late Cretaceous/Paleocene. This occurrence of accreted seamount fragments and morphology of the Pacific Ocean floor is indicative of the formation of the Cocos-Nazca spreading system at least ~40 Ma prior to the age proposed in current tectonic models. In Panama, we identified a remarkably-well preserved early Eocene oceanic island that accreted in the middle Eocene. The accretion probably occurred at very shallow depth by detachment of the island in the trench and led to an exceptional preservation of the volcanic structures. Exposures of both deep and superficial parts of the volcanic edifice have been studied, from the submarine-shield to subaerial-postshield stages. The stratigraphy allowed us to distinguish lavas produced during the submarine and subaerial stages. The lava compositions likely define a progressive diminution of source melting and a decrease in the temperature of erupted melts in the latest stages of volcanic activity. We interpret these changes to primarily reflect the progressive migration of the oceanic island out of the melting region or hotspot.

Mesozoic oceanic terranes of southern Central America : geology, geochemistry and gedynamics

Abstract The Northwestern edge of the modern Caribbean Plate, located in central Middle America (S-Guatemala to N-Costa Rica), is characterized by a puzzle of oceanic and continental terranes that belonged originally to the Pacific façade of North America. South of the Motagua Fault Zone, the actual northern strike slip boundary of the Caribbean Plate, three continental slivers (Copán, Chortis s. str. and Patuca) are sandwiched between two complex suture zones that contain HP/LT mafic and ultramafic oceanic rocks: The Motagua Mélanges to the North, extensively studied in the last ten years and the' newly defined Mesquito Composite Oceanic Terrane (MCOT) to the South. No modem geological data were available for the oceanic terrane located in the southern part of the so called continental "Chortis Block". Classically, the southern limit of this block with the Caribbean Large Igneous Province (CLIP) was placed at a hypothetical fault line connecting the main E-W fault in the Santa Elena Peninsula (N-Costa Rica) with the Hess Escarpment. However, our study in eastern Nicaragua and northwestern Costa Rica evidences an extensive assemblage of oceanic upper mantle and crustal rocks outcropping between the Chortis/Patuca continental blocks and the CLIP. They comprise collided and accreted exotic terranes of Pacific origin recording a polyphased tectonic history. We distinguish: 1- The MCOT that comprises a Late Triassic to Early Cretaceous puzzle of oceanic crust and arc-derived rocks set in a serpentinite matrix, and 2- The Manzanillo and Nicoya Terranes that are made of Cretaceous plateau-like rocks associated with oceanic sediments older than the CLIP. This study has been focused on the rocks of the MCOT. The MCOT comprises the southern half of the former "Chortis Block" and is defined by 4 comer localities characterized by ultramafic and mafic oceanic rocks of Late Triassic, Jurassic and Early Cretaceous age: 1- The Siuna Serpentinite Mélange (NE-Nicaragua), 2- The El Castillo Mélange (Nicaragua/Costa Rica border), 3- DSDP Legs 67 and 84 (Guatemala fore-arc basin), and 4- The Santa Elena Peridiotite (NW-Costa Rica). The Siuna Serpentinite Mélange (SSM) is a HP/LT subduction zone mélange set in a serpentinite matrix that contains oceanic crust and arc-related greenschist to blueschist/eclogite facies metamafic and metasedimentary blocks. Middle Jurassic (Bajocian-Bathonian) radiolarites are found in original sedimentary contact with arc-derived greenstones. Late Jurassic black detrital chert possibly formed in a marginal (fore-arc?) basin shortly before subduction. A phengite 40Ar/39Ar -cooling age dates the exhumation of the high pressure rocks as 139 Ma. The El Castillo Mélange (ECM) is composed of serpentinite matrix with OIB metabasalts and Late Triassic (Rhaetian) red and green radiolarite blocks. Recent studies of the DSDP Legs 67/84 show that the Guatemala/Nicaragua fore-arc basin is composed of a pile of ultramafic, mafic (OIB-like) and arc related rocks with ages ranging from Late Triassic to Campanian. Finally, the Santa Elena peridiotites that mark the limit of the MCOT with the Manzanillo/Nicoya Terranes and correspond to an association of ultramafic rocks that comprise peridiotites, dunites and chromites of abyssal and fore-arc origin. The SSM is the result of a collision between a Middle Jurassic island arc and the Patuca Terrane, a fragment of the Western N-American active continental margin. The Siuna Mélange (SSM) and the South Montagna Mélange share common characteristics with the Pacific N-American suture zone (E-Franciscan and Vizcaino mélanges), in particular, the Mesozoic ages of HP/LT metamorphic and the arc-derived blocks. For us, these mélanges imply an originally continuous, but slightly diachronous suture that affected the entire W-American active margin. It may imply the arrival and collision of an exotic intraoceanic arc (Guerrero-Phoenix) related to the origin of the Pacific Plate that initiated as a back arc basin of this arc. The present disposition of the fragments of this suture zone is the result of a northward shift of the active left-lateral strike slip motion between the N-American and the Caribbean Plates. Résumé Le coin nord-ouest de la Plaque Caraïbe moderne se trouve en Amérique Centrale, entre le sud du Guatemala et le nord du Costa Rica. Cette région est composée d'un puzzle de terrains océaniques et continentaux dont les origines se situent sur la façade pacifique de l'Amérique du Nord. Au sud de la faille de Motagua, la limite septentrionale actuelle, décrochante, de la Plaque Caraïbe, se trouvent 3 copeaux continentaux (Copàn, Chortis s. str. et Patuca) coincés entre deux zones de suture complexes à roches mafiques et ultramafiques qui ont subi un métamorphisme de haute pression/basse température (HP/LT). Il s'agit des Mélanges de Motagua au nord, largement étudiés ces dernières années, et du Mesquito Composite Oceanic Terrane (MCOT), récemment défini par nous, au sud. En vue de l'absence de données géologiques modernes concernant les terrains océaniques qui se trouvent dans la partie sud du "Chortis Block" considérée comme continentale, nous avons dédié cette étude à cette région. Classiquement, la limite méridionale entre le "Chortis Block" et la "Caribbean Large Igneous Province" (CLIP) a été associée à une faille hypothétique reliant la faille E-W de Santa Elena (nord du Costa Rica) à l'Escarpement de Hess. Notre étude au Nicaragua oriental et au Costa Rica nord-occidental a révélé l'existence de larges terrains composés d'assemblages de roches mantéliques et océaniques qui se placent entre les blocs continentaux Chortis/Patuca et le CLIP. Ces assemblages révèlent des terrains collisionnés et accrétés d'origine pacifique enregistrant une histoire tectonique polyphasée. Nous distinguons: 1- Le MCOT, un puzzle de roches océaniques d'arc d'âge Triassique supérieur au Crétacée inférieur, 2- Les terrains de Manzanillo et de Nicoya, des morceaux de plateaux océaniques associés à des sédiments océaniques plus âgés que le CLIP. Cette étude se focalisera sur les roches du MCOT. Le MCOT occupe la moitié sud de l'ancien "Chortis Block" et peut se définir par 4 localités de référence qui montrent des roches mafiques et ultramafiques océaniques d'âges compris entre le Trias supérieur et le Crétacée inférieur. 1- Le Siuna Serpentinite Mélange (NE-Nicaragua), 2- Le El Castillo Mélange (Nicaragua/Costa Rica border), 3- Le DSDP Legs 67/84 (Guatemala fore-arc basin) et 4- La Santa Elena Peridiotite (nord-ouest du Costa Rica). Le Siuna Serpentinite Mélange (SSM) est un mélange de subduction HP/BT dans une matrice de serpentinite. On y trouve des éléments de croûte océanique et d'arc insulaire en faciès de schistes verts et schistes bleus. Des radiolarites du Jurassique moyen se trouvent en contact sédimentaire sur des roches vertes d'arc. En revanche, des cherts noirs détritiques datent du Jurassique supérieur et sont probablement issus d'un bassin marginal (fore-arc ?) peu avant leur subduction, car un âge 40Ar/39Ar de refroidissement des phengites date l'exhumation des roches de haute pression à 139 Ma. Le Mélange d'El Castillo (ECM) est constitué d'une matrice serpentinitique et contient des blocs de metabasaltes OIB et des blocs de radiolarites du Trias terminal. Des études récentes ont repris les roches forées lors des DSDP Legs 67 et 84 et montrent que le soubassement du bassin d'avant-arc du Guatemala-Nicaragua est composé de roches ultramafiques et mafiques (OIB et arc), dont les âges isotopiques vont du Trias au Crétacé supérieur. Finalement, les péridiotites de Santa Elena forment la limite sud du MCOT par rapport aux terrains de Manzanillo et Nicoya. Elles contiennent des serpentinites et localement des dunites et chromites à affinité abyssale et de fore-arc. Le SSM témoigne d'une collision entre un arc insulaire d'âge Jurassique moyen et le Patuca Terrane, un fragment de la marge active nord-américaine. Le SSM et le South Motagua Mélange ont des caractéristiques en commun avec les zones de suture de la façade pacifique de l'Amérique du nord (E-Franciscan et Vizcaino mélanges), notamment les âges Mésozoïques du métamorphisme HP/BT et les blocs de roches d'arc. Ce fait nous conduit à penser qu'il s'agit d'une grande zone de suture qui était à l'origine continue sur toute la marge ouest-américaine, mais légèrement diachrone. Cette suture implique l'arrivée et la collision d'un arc intraocéanique exotique (Guerrero-Phoenix) qui est à l'origine de la Plaque Pacifique qui s'ouvrait en back arc par rapport à celui-ci. La disposition actuelle des fragments de cette suture est due à la migration vers le nord du décrochement actif senestre entre la Plaque nord-américaine et la Plaque Caraïbe. K. Flores, 2009 Mesozoic oceanic terranes of southern central America Résumé Grand Public La présente thèse est le résultat de travaux de terrain effectués de 2005 à 2008 au nord-est et au sud du Nicaragua et au nord du Costa Rica, en Amérique Centrale, des analyses pétrologiques et géochimiques en laboratoire ainsi que de la modélisation de l'évolution géodynamique. La région étudiée se situe en bordure nord - ouest de la Plaque Caraïbe moderne. Dans la majorité des publications récentes cette région est représentée comme un vaste bloc continental (le "Bloc Chortis") qui serait limité, (i) au nord, par la faille décrochante de Motagua, la limite actuelle entre la Plaque Nord-Américaine et la Plaque Caraïbe, et (ii) au sud, par une suture hypothétique qui se trouverait aux confins entre le Nicaragua et le Costa Rica. La région du Costa Rica a été considérée presque entièrement comme une partie du Plateau Caraïbe ("Caribbean Large Igneous Province" (CLIP)). L'étude détaillée des affleurements nous a permis de mettre en évidence : - Au nord-est du Nicaragua (Siuna) : Des roches océaniques datées du Jurassique moyen, grâce aux faunes à radiolaires qui ont été extraites des radiolarites rouges. Ces roches ont subi un métamorphisme de haute pression typique des zones de collision. L'étude radio-isotopique Ar/Ar a permis de dater la collision du Crétacé basal (139 Ma). - Au sud du Nicaragua : Des roches océaniques d'âge Trias terminal (200 millions d'années), également datées à l'aide de faunes à radiolaires. Il s'agit actuellement des roches océaniques les plus anciennes connues de l'Amérique Centrale. - L'étude géochimique et les âges des fossiles démontrent que le tiers septentrional du Costa Rica possède un soubassement construit d'au moins deux terrains (Nicoya et Manzanillo), qui ont des caractéristiques de Plateau océanique (Nicoya) et d'arc volcanique du Crétacé moyen (Manzanillo). Ces deux terrains sont plus anciens que le CLIP. En conclusion, nous constatons que la région étudiée est constituée d'un puzzle de 3 blocs continentaux et d'un vaste terrain océanique composite que nous appelons Mesquito Composite Oceanic Terrane (MCOT). En plus, nous définissons les terrains de Nicoya et de Manzanillo comme plus âgés et distincts du CLIP. Le MCOT est caractérisé par la présence de roches du manteau supérieur (les serpentinites) et de la croûte océanique, ainsi que des morceaux d'arcs, d'âge allant du Trias supérieur au Crétacé. Ce terrain est comparable à d'autres zones de suture de la façade pacifique de l'Amérique du nord, notamment en ce qui concerne les âges Mésozoïques, le métamorphisme de haute pression et l'association de roches mantéliques et crustales océaniques. Ce fait nous conduit à penser qu'il s'agit d'une grande zone de suture qui était à l'origine continue sur toute la marge ouest-américaine. Cette suture implique l'arrivée et la collision d'un arc infra-océanique exotique qui serait à l'origine de la Plaque Pacifique qui se serait ouverte en bassin d'arrière arc par rapport à celui-ci. La disposition actuelle des fragments de cette suture est due à la migration vers le nord du décrochement actif senestre entre la Plaque nord-américaine et la Plaque Caraïbe.

Aggregation invariance in general clustering approaches

General clustering deals with weighted objects and fuzzy memberships. We investigate the group- or object-aggregation-invariance properties possessed by the relevant functionals (effective number of groups or objects, centroids, dispersion, mutual object-group information, etc.). The classical squared Euclidean case can be generalized to non-Euclidean distances, as well as to non-linear transformations of the memberships, yielding the c-means clustering algorithm as well as two presumably new procedures, the convex and pairwise convex clustering. Cluster stability and aggregation-invariance of the optimal memberships associated to the various clustering schemes are examined as well.

A case for hornblende dominated fractionation of arc magmas: the Chelan Complex (Washington Cascades)

Amphibole fractionation in the deep roots of subduction-related magmatic arcs is a fundamental process for the generation of the continental crust. Field relations and geochemical data of exposed lower crustal igneous rocks can be used to better constrain these processes. The Chelan Complex in the western U. S. forms the lowest level of a 40-km thick exposed crustal section of the North Cascades and is composed of olivine websterite, pyroxenite, hornblendite, and dominantly by hornblende gabbro and tonalite. Magmatic breccias, comb layers and intrusive contacts suggest that the Chelan Complex was build by igneous processes. Phase equilibria, textural observations and mineral chemistry yield emplacement pressures of similar to 1.0 GPa followed by isobaric cooling to 700 degrees C. The widespread occurrence of idiomorphic hornblende and interstitial plagioclase together with the lack of Eu anomalies in bulk rock compositions indicate that the differentiation is largely dominated by amphibole. Major and trace element modeling constrained by field observations and bulk chemistry demonstrate that peraluminous tonalite could be derived by removing successively 3% of olivine websterite, 12% of pyroxene hornblendite, 33% of pyroxene hornblendite, 19% of gabbros, 15% of diorite and 2% tonalite. Peraluminous tonalite with high Sr/Y that are worldwide associated with active margin settings can be derived from a parental basaltic melt by crystal fractionation at high pressure provided that amphibole dominates the fractionation process. Crustal assimilation during fractionation is thus not required to generate peraluminous tonalite.

Ruin problems under IBNR Dynamics

In this paper, we consider a discrete-time risk process allowing for delay in claim settlement, which introduces a certain type of dependence in the process. From martingale theory, an expression for the ultimate ruin probability is obtained, and Lundberg-type inequalities are derived. The impact of delay in claim settlement is then investigated. To this end, a convex order comparison of the aggregate claim amounts is performed with the corresponding non-delayed risk model, and numerical simulations are carried out with Belgian market data.

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.

An efficient total variation algorithm for super-resolution in fetal brain MRI with adaptive regularization.

Although fetal anatomy can be adequately viewed in new multi-slice MR images, many critical limitations remain for quantitative data analysis. To this end, several research groups have recently developed advanced image processing methods, often denoted by super-resolution (SR) techniques, to reconstruct from a set of clinical low-resolution (LR) images, a high-resolution (HR) motion-free volume. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has been quite attracted by Total Variation energies because of their ability in edge preserving but only standard explicit steepest gradient techniques have been applied for optimization. In a preliminary work, it has been shown that novel fast convex optimization techniques could be successfully applied to design an efficient Total Variation optimization algorithm for the super-resolution problem. In this work, two major contributions are presented. Firstly, we will briefly review the Bayesian and Variational dual formulations of current state-of-the-art methods dedicated to fetal MRI reconstruction. Secondly, we present an extensive quantitative evaluation of our SR algorithm previously introduced on both simulated fetal and real clinical data (with both normal and pathological subjects). Specifically, we study the robustness of regularization terms in front of residual registration errors and we also present a novel strategy for automatically select the weight of the regularization as regards the data fidelity term. Our results show that our TV implementation is highly robust in front of motion artifacts and that it offers the best trade-off between speed and accuracy for fetal MRI recovery as in comparison with state-of-the art methods.

Efficient total variation algorithm for fetal MRI reconstruction

Fetal MRI reconstruction aims at finding a high-resolution image given a small set of low-resolution images. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has considered several regularization terms s.a. Dirichlet/Laplacian energy [1], Total Variation (TV)based energies [2,3] and more recently non-local means [4]. Although TV energies are quite attractive because of their ability in edge preservation, standard explicit steepest gradient techniques have been applied to optimize fetal-based TV energies. The main contribution of this work lies in the introduction of a well-posed TV algorithm from the point of view of convex optimization. Specifically, our proposed TV optimization algorithm for fetal reconstruction is optimal w.r.t. the asymptotic and iterative convergence speeds O(1/n(2)) and O(1/root epsilon), while existing techniques are in O(1/n) and O(1/epsilon). We apply our algorithm to (1) clinical newborn data, considered as ground truth, and (2) clinical fetal acquisitions. Our algorithm compares favorably with the literature in terms of speed and accuracy.

Structured sparsity through reweighting and application to diffusion MRI

We consider the problem of multiple correlated sparse signals reconstruction and propose a new implementation of structured sparsity through a reweighting scheme. We present a particular application for diffusion Magnetic Resonance Imaging data and show how this procedure can be used for fibre orientation reconstruction in the white matter of the brain. In that framework, our structured sparsity prior can be used to exploit the fundamental coherence between fibre directions in neighbour voxels. Our method approaches the ℓ0 minimisation through a reweighted ℓ1-minimisation scheme. The weights are here defined in such a way to promote correlated sparsity between neighbour signals.

On evaluating the accuracy and biological plausibility of diffusion MRI tractograms

In diffusion MRI, traditional tractography algorithms do not recover truly quantitative tractograms and the structural connectivity has to be estimated indirectly by counting the number of fiber tracts or averaging scalar maps along them. Recently, global and efficient methods have emerged to estimate more quantitative tractograms by combining tractography with local models for the diffusion signal, like the Convex Optimization Modeling for Microstructure Informed Tractography (COMMIT) framework. In this abstract, we show the importance of using both (i) proper multi-compartment diffusion models and (ii) adequate multi-shell acquisitions, in order to evaluate the accuracy and the biological plausibility of the tractograms.