Intraobserver and interobserver variability of ascending aorta diameter as assessed with ECG-gated MDCT: automatic versus manual measurements

Purpose: Recently morphometric measurements of the ascending aorta have been done with ECG-gated MDCT to help the development of future endovascular therapies (TCT) [1]. However, the variability of these measurements remains unknown. It will be interesting to know the impact of CAD (computer aided diagnosis) with automated segmentation of the vessel and automatic measurements of diameter on the management of ascending aorta aneurysms. Methods and Materials: Thirty patients referred for ECG-gated CT thoracic angiography (64-row CT scanner) were evaluated. Measurements of the maximum and minimum ascending aorta diameters were obtained automatically with a commercially available CAD and semi-manually by two observers separately. The CAD algorithms segment the iv-enhanced lumen of the ascending aorta into perpendicular planes along the centreline. The CAD then determines the largest and the smallest diameters. Both observers repeated the automatic measurements and the semimanual measurements during a different session at least one month after the first measurements. The Bland and Altman method was used to study the inter/intraobserver variability. A Wilcoxon signed-rank test was also used to analyse differences between observers. Results: Interobserver variability for semi-manual measurements between the first and second observers was between 1.2 to 1.0 mm for maximal and minimal diameter, respectively. Intraobserver variability of each observer ranged from 0.8 to 1.2 mm, the lowest variability being produced by the more experienced observer. CAD variability could be as low as 0.3 mm, showing that it can perform better than human observers. However, when used in nonoptimal conditions (streak artefacts from contrast in the superior vena cava or weak lumen enhancement), CAD has a variability that can be as high as 0.9 mm, reaching variability of semi-manual measurements. Furthermore, there were significant differences between both observers for maximal and minimal diameter measurements (p<0.001). There was also a significant difference between the first observer and CAD for maximal diameter measurements with the former underestimating the diameter compared to the latter (p<0.001). As for minimal diameters, they were higher when measured by the second observer than when measured by CAD (p<0.001). Neither the difference of mean minimal diameter between the first observer and CAD nor the difference of mean maximal diameter between the second observer and CAD was significant (p=0.20 and 0.06, respectively). Conclusion: CAD algorithms can lessen the variability of diameter measurements in the follow-up of ascending aorta aneurysms. Nevertheless, in non-optimal conditions, it may be necessary to correct manually the measurements. Improvements of the algorithms will help to avoid such a situation.

Effect of Sodium Loading/Depletion on Renal Oxygenation in Young Normo-and Hypertensive Men Measured with BOLD-MRI

Objective: To measure renal tissue oxygenation in young normo-and hypertensive volunteers under conditions of salt loading and depletion using blood oxygen level dependent magnetic resonance imaging (BOLD-MRI). Design and Methods: Ten normotensive (NT) male volunteers (age 26.5_7.4 y) and eight non-treated, hypertensive (HT) male volunteers (age 28.8_5.7 y) were studied after one week on a high salt (HS) regimen (6g of salt/day added to their normal regimen) and again after one week of a low sodium diet (LS). On the 8th day, BOLD-MRI was performed under standard hydration conditions. Four coronal slices were selected in each kidney, and combination sequence was used to acquire T2* weighted images. The mean R2* (1/T2*) was measured to determine cortical and medullar oxygenation. Results: Baseline characteristics and their changes are shown in the table. The mean cortical R2* was not different under conditions of HS or LS (17.8_1.3 vs. 18.2_0.6 respectively in NT group, p_0.27; 17.4_0.6 vs 17.8_0.9 in HT group, p_0.16). However, the mean medullary R2* was significantly lower under LS conditions in both groups (31.3_0.6 vs 28.1_0.8 in NT group, p_0.05; 30.3_0.8 vs 27.9_1.5 in HT group, p_0.05), corresponding to higher medullary oxygenation as compared to HS conditions, without significant changes in hemoglobin or hematocrit values. The salt induced changes in medullary oxygenation were comparable in the two groups (ANOVA, p_0.1). Conclusion: Dietary sodium restriction leads to increased renal medullary oxygenation compared to high sodium intake in normo-and hypertensive subjects. This observation may in part explain the potential renal benefits of a low sodium intake.

Ambulatory measurement of ankle kinetics for clinical applications.

This study aimed to design and validate the measurement of ankle kinetics (force, moment, and power) during consecutive gait cycles and in the field using an ambulatory system. An ambulatory system consisting of plantar pressure insole and inertial sensors (3D gyroscopes and 3D accelerometers) on foot and shank was used. To test this system, 12 patients and 10 healthy elderly subjects wore shoes embedding this system and walked many times across a gait lab including a force-plate surrounded by seven cameras considered as the reference system. Then, the participants walked two 50-meter trials where only the ambulatory system was used. Ankle force components and sagittal moment of ankle measured by ambulatory system showed correlation coefficient (R) and normalized RMS error (NRMSE) of more than 0.94 and less than 13% in comparison with the references system for both patients and healthy subjects. Transverse moment of ankle and ankle power showed R>0.85 and NRMSE<23%. These parameters also showed high repeatability (CMC>0.7). In contrast, the ankle coronal moment of ankle demonstrated high error and lower repeatability. Except for ankle coronal moment, the kinetic features obtained by the ambulatory system could distinguish the patients with ankle osteoarthritis from healthy subjects when measured in 50-meter trials. The proposed ambulatory system can be easily accessible in most clinics and could assess main ankle kinetics quantities with acceptable error and repeatability for clinical evaluations. This system is therefore suggested for field measurement in clinical applications.

High-resolution selective three-dimensional magnetic resonance coronary angiography with navigator-echo technique: segment-by-segment evaluation of coronary artery stenosis.

PURPOSE: To investigate the feasibility of high-resolution selective three-dimensional (3D) magnetic resonance coronary angiography (MRCA) in the evaluation of coronary artery stenoses. MATERIALS AND METHODS: In 12 patients with coronary artery stenoses, MRCA of the coronary artery groups, including the coronary segments with stenoses of 50% or greater based on conventional x-ray coronary angiography (CAG), was performed with double-oblique imaging planes by orienting the 3D slab along the major axis of each right coronary artery-left circumflex artery (RCA-LCX) group and each left main trunk-left anterior descending artery (LMT-LAD) group. Ten RCA-LCX and five LMT-LAD MR angiograms were obtained, and the results were compared with those of conventional x-ray angiography. RESULTS: Among 70 coronary artery segments expected to be covered, a total of 49 (70%) segments were fully demonstrated in diagnostic quality. The identification of segmental location of stenoses showed as high an accuracy as 96%. The retrospective analysis for stenosis of 50% or greater yielded the sensitivity, specificity, and accuracy of 80%, 85%, and 84%, respectively. CONCLUSION: Selective 3D MRCA has the potential for segment-by-segment evaluation of major portions of the right and left coronary arteries with high accuracy.

Structural characterization of Turtle Mountain anticline (Alberta, Canada) and impact on rock slope failure

his paper proposes a structural investigation of the Turtle Mountain anticline (Alberta, Canada) to better understand the role of the different tectonic features on the development of both local and large scale rock slope instabilities occurring in Turtle Mountain. The study area is investigated by combining remote methods with detailed field surveys. In particular, the benefit of Terrestrial Laser Scanning for ductile and brittle tectonic structure interpretations is illustrated. The proposed tectonic interpretation allows the characterization of the fracturing pattern, the fold geometry and the role of these tectonic features in rock slope instability development. Ten discontinuity sets are identified in the study area, their local variations permitting the differentiation of the study zone into 20 homogenous structural domains. The anticline is described as an eastern verging fold that displays considerable geometry differences along its axis and developed by both flexural slip and tangential longitudinal strain folding mechanisms. Moreover, the origins of the discontinuity sets are determined according to the tectonic phases affecting the region (pre-folding, folding, post-folding). The localization and interpretation of kinematics of the different instabilities revealed the importance of considering the discrete brittle planes of weakness, which largely control the kinematic release of the local instabilities, and also the rock mass damage induced by large tectonic structures (fold hinge, thrust).

Influence of the ankle arthropathies on the 3D kinematics of the lower limbs during gait: analysis with an ambulatory system.

Introduction: Ankle arthropathy is associated with a decreased motion of the ankle-hindfoot during ambulation. Ankle arthrodesis was shown to result in degeneration of the neighbour joints of the foot. Inversely, total ankle arthroplasty conceptually preserves the adjacent joints because of the residual mobility of the ankle but this has not been demonstrated yet in vivo. It has also been reported that degenerative ankle diseases, and even arthrodesis, do not result in alteration of the knee and hip joints. We present the preliminary results of a new approach of this problem based on ambulatory gait analysis. Patients and Methods: Motion analysis of the lower limbs was performed using a Physilog® (BioAGM, CH) system consisting of three-dimensional (3D) accelerometer and gyroscope, coupled to a magnetic system (Liberty©, Polhemus, USA). Both systems have been validated. Three groups of two patients were included into this pilot study and compared to healthy subjects (controls) during level walking: patients with ankle osteoarthritis (group 1), patients treated by ankle arthrodesis (group 2), patients treated by total ankle prosthesis (group 3). Results: Motion patterns of all analyzed joints over more than 20 gait cycles in each subject were highly repeatable. Motion amplitude of the ankle-hindfoot in control patients was similar to recently reported results. Ankle arthrodesis limited the motion of the ankle-hindfoot in the sagittal and horizontal planes. The prosthetic ankle allowed a more physiologic movement in the sagittal plane only. Ankle arthritis and its treatments did not influence the range of motion of the knee and hip joint during stance phase, excepted for a slight decrease of the hip flexion in groups 1 and 2. Conclusion: The reliability of the system was shown by the repeatability of the consecutive measurements. The results of this preliminary study were similar to those obtained through laboratory gait analysis. However, our system has the advantage to allow ambulatory analysis of 3D kinematics of the lower limbs outside of a gait laboratory and in real life conditions. To our knowledge this is a new concept in the analysis of ankle arthropathy and its treatments. Therefore, there is a potential to address specific questions like the difficult comparison of the benefits of ankle arthroplasty versus arthrodesis. The encouraging results of this pilot study offer the perspective to analyze the consequences of ankle arthropathy and its treatments on the biomechanics of the lower limbs ambulatory, in vivo and in daily life conditions.

Dorsal fractures of the triquetrum: MRI findings with an emphasis on dorsal carpal ligament injuries.

OBJECTIVE: The objective of our study was to report the MRI findings in dorsal fractures of the triquetrum, with an emphasis on dorsal carpal ligament injuries. MATERIALS AND METHODS: A total of 21 patients (16 men, five women; mean age, 41.9 years) with acute or subacute (≤ 6 weeks) dorsal triquetral fractures on radiography and MRI were included in this two-center retrospective study. MRI of the wrist was performed on 3-T units with transverse T1-weighted, coronal or transverse (or both) fat-suppressed T2weighted, transverse gadolinium-enhanced fat-suppressed T1-weighted turbo spin-echo, and 3D gadolinium-enhanced fat-suppressed T1-weighted gradient-recalled echo sequences. Three musculoskeletal radiologists evaluated the ulnar styloid process index (USPI) on radiographs and the following MRI features: fracture pattern (types 1-6), bone fragment size and displacement, bone marrow edema distribution, and dorsal carpal ligament tears. RESULTS: Eight type 1, one type 2, six type 3, five type 4, and one type 5 fractures were identified. These fractures were associated with 14 (66.7%), 17 (81.0%), and 16 (76.2%) tears of the dorsal radiocarpal, ulnotriquetral, and intercarpal ligaments, respectively. There was no correlation between bone marrow edema distribution and dorsal carpal ligament injuries (all p > 0.05). The mean (± SD) bone fragment volume and displacement were 205 ± 157 mm(3) and 1.0 ± 1.1 mm, respectively. The mean USPI was 0.21 ± 0.10. CONCLUSION: Dorsal fractures of the triquetrum are frequently associated with dorsal carpal ligament injuries. Bone marrow edema distribution is not correlated with these ligament tears.

Formation of chlorite during thrust fault reactivation. Record of fluid origin and P-T conditions in the Monte Perdido thrust fault (southern Pyrenees)

The chemical and isotopic compositions of clay minerals such as illite and chlorite are commonly used to quantify diagenetic and low-grade metamorphic conditions, an approach that is also used in the present study of the Monte Perdido thrust fault from the South Pyrenean fold-and-thrust belt. The Monte Perdido thrust fault is a shallow thrust juxtaposing upper Cretaceous-Paleocene platform carbonates and Lower Eocene marls and turbidites from the Jaca basin. The core zone of the fault, about 6 m thick, consists of intensely deformed clay-bearing rocks bounded by major shear surfaces. Illite and chlorite are the main hydrous minerals in the fault zone. Illite is oriented along cleavage planes while chlorite formed along shear veins (< 50 mu m in thickness). Authigenic chlorite provides essential information about the origin of fluids and their temperature. delta O-18 and delta D values of newly formed chlorite support equilibration with sedimentary interstitial water, directly derived from the local hanging wall and footwall during deformation. Given the absence of large-scale fluid flow, the mineralization observed in the thrust faults records the P-T conditions of thrust activity. Temperatures of chlorite formation of about 240A degrees C are obtained via two independent methods: chlorite compositional thermometers and oxygen isotope fractionation between cogenetic chlorite and quartz. Burial depth conditions of 7 km are determined for the Monte Perdido thrust reactivation, coupling calculated temperature and fluid inclusion isochores. The present study demonstrates that both isotopic and thermodynamic methods applied to clay minerals formed in thrust fault are useful to help constrain diagenetic and low-grade metamorphic conditions.

A comparison between joint coordinate system and attitude vector for multi-segment foot kinematics.

The joint angles of multi-segment foot models have been primarily described using two mathematical methods: the joint coordinate system and the attitude vector. This study aimed to determine whether the angles obtained through these two descriptors are comparable, and whether these descriptors have similar sensitivity to experimental errors. Six subjects walked eight times on an instrumented walkway while the joint angles among shank, hindfoot, medial forefoot, and lateral forefoot were measured. The angles obtained using both descriptors and their sensitivity to experimental errors were compared. There was no overall significant difference between the ranges of motion obtained using both descriptors. However, median differences of more than 6° were noticed for the medial-lateral forefoot joint. For all joints and rotation planes, both descriptors provided highly similar angle patterns (median correlation coefficient: R>0.90), except for the medial-lateral forefoot angle in the transverse plane (median R=0.77). The joint coordinate system was significantly more sensitive to anatomical landmarks misplacement errors. However, the absolute differences of sensitivity were small relative to the joints ranges of motion. In conclusion, the angles obtained using these two descriptors were not identical, but were similar for at least the shank-hindfoot and hindfoot-medial forefoot joints. Therefore, the angle comparison across descriptors is possible for these two joints. Comparison should be done more carefully for the medial-lateral forefoot joint. Moreover, despite different sensitivities to experimental errors, the effects of the experimental errors on the angles were small for both descriptors suggesting that both descriptors can be considered for multi-segment foot models.

Breathhold three-dimensional coronary magnetic resonance angiography using real-time navigator technology.

The acquisition duration of most three-dimensional (3D) coronary magnetic resonance angiography (MRA) techniques is considerably prolonged, thereby precluding breathholding as a mechanism to suppress respiratory motion artifacts. Splitting the acquired 3D volume into multiple subvolumes or slabs serves to shorten individual breathhold duration. Still, problems associated with misregistration due to inconsistent depths of expiration and diaphragmatic drift during sustained respiration remain to be resolved. We propose the combination of an ultrafast 3D coronary MRA imaging sequence with prospective real-time navigator technology, which allows correction of the measured volume position. 3D volume splitting using prospective real-time navigator technology, was successfully applied for 3D coronary MRA in five healthy individuals. An ultrafast 3D interleaved hybrid gradient-echoplanar imaging sequence, including T2Prep for contrast enhancement, was used with the navigator localized at the basal anterior wall of the left ventricle. A 9-cm-thick volume, with in-plane spatial resolution of 1.1 x 2.2 mm, was acquired during five breathholds of 15-sec duration each. Consistently, no evidence of misregistration was observed in the images. Extensive contiguous segments of the left anterior descending coronary artery (48 +/- 18 mm) and the right coronary artery (75 +/- 5 mm) could be visualized. This technique has the potential for screening for anomalous coronary arteries, making it well suited as part of a larger clinical MR examination. In addition, this technique may also be applied as a scout scan, which allows an accurate definition of imaging planes for subsequent high-resolution coronary MRA.

Geological transect across the Northwestern Himalaya in eastern Ladakh and Lahul (A model for the continental collision of India and Asia)

The detailed geological mapping and structural study of a complete transect across the northwestern Himalaya allow to describe the tectonic evolution of the north Indian continental margin during the Tethys ocean opening and the Himalayan Orogeny. The Late Paleozoic Tethys rifting is associated with several tectonomagmatic events. In Upper Lahul and SE Zanskar, this extensional phase is recorded by Lower Carboniferous synsedimentary transtensional faults, a Lower Permian stratigraphic unconformity, a Lower Permian granitic intrusion and middle Permian basaltic extrusions (Panjal Traps). In eastern Ladakh, a Permian listric normal fault is also related to this phase. The scarcity of synsedimentary faults and the gradual increase of the Permian syn-rift sediment thickness towards the NE suggest a flexural type margin. The collision of India and Asia is characterized by a succession of contrasting orogenic phases. South of the Suture Zone, the initiation of the SW vergent Nyimaling-Tsarap Nappe corresponds to an early phase of continental underthrusting. To the S, in Lahul, an opposite underthrusting within the Indian plate is recorded by the NE vergent Tandi Syncline. This structure is associated with the newly defined Shikar Beh Nappe, now partly eroded, which is responsible for the high grade (amphibolite facies) regional metamorphism of South Lahul. The main thrusting of the Nyimaling-Tsarap Nappe followed the formation of the Shikar Beh Nappe. The Nyimaling-Tsarap Nappe developed by ductile shear of the upper part of the subducted Indian continental margin and is responsible for the progressive regional metamorphism of SE Zanskar, reaching amphibolite facies below the frontal part of the nappe, near Sarchu. In Upper Lahul, the frontal parts of the Nyimaling-Tsarap and Shikar Beh nappes are separated by a zone of low grade metamorphic rocks (pumpellyite-actinolite facies to lower greenschist facies). At high structural level, the Nyimaling-Tsarap Nappe is characterized by imbricate structures, which grade into a large ductile shear zone with depth. The related crustal shortening is about 87 km. The root zone and the frontal part of this nappe have been subsequently affected by two zones of dextral transpression and underthrusting: the Nyimaling Shear Zone and the Sarchu Shear Zone. These shear zones are interpreted as consequences of the counterclockwise rotation of the continental underthrusting direction of India relative to Asia, which occurred some 45 and 36 Ma ago, according to plate tectonic models. Later, a phase of NE vergent `'backfolding'' developed on these two zones of dextral transpression, creating isoclinal folds in SE Zanskar and more open folds in the Nyimaling Dome and in the Indus Molasse sediments. During a late stage of the Himalayan Orogeny, the frontal part of the Nyimaling-Tsarap Nappe underwent an extension of about 15 km. This phase is represented by two types of structures, responsible for the tectonic unroofing of the amphibolite facies rocks of the Sarchu area: the Sarchu high angle Normal Fault, cutting a first set of low angle normal faults, which have been created by reactivation of older thrust planes related to the Nyimaling-Tsarap Nappe.

Effect of sodium loading/depletion on renal oxygenation in young normotensive and hypertensive men.

The goal of this study was to investigate the effect of sodium intake on renal tissue oxygenation in humans. To this purpose, we measured renal hemodynamics, renal sodium handling, and renal oxygenation in normotensive (NT) and hypertensive (HT) subjects after 1 week of a high-sodium and 1 week of a low-sodium diet. Renal oxygenation was measured using blood oxygen level-dependent magnetic resonance. Tissue oxygenation was determined by the measurement of R2* maps on 4 coronal slices covering both kidneys. The mean R2* values in the medulla and cortex were calculated, with a low R2* indicating a high tissue oxygenation. Ten male NT (mean age: 26.5+/-7.4 years) and 8 matched HT subjects (mean age: 28.8+/-5.7 years) were studied. Cortical R2* was not different under the 2 conditions of salt intake. Medullary R2* was significantly lower under low sodium than high sodium in both NT and HT subjects (28.1+/-0.8 versus 31.3+/-0.6 s(-1); P<0.05 in NT; and 27.9+/-1.5 versus 30.3+/-0.8 s(-1); P<0.05, in HT), indicating higher medullary oxygenation under low-sodium conditions. In NT subjects, medullary oxygenation was positively correlated with proximal reabsorption of sodium and negatively with absolute distal sodium reabsorption, but not with renal plasma flow. In HT subjects, medullary oxygenation correlated with the 24-hour sodium excretion but not with proximal or with the distal handling of sodium. These data demonstrate that dietary sodium intake influences renal tissue oxygenation, low sodium intake leading to an increased renal medullary oxygenation both in normotensive and young hypertensive subjects.

Physical considerations on discrepancies in target volume delineation.

BACKGROUND AND PURPOSE: To compare the delineations and interpretations of target volumes by physicians in different radio-oncology centers. MATERIALS AND METHODS: Eleven Swiss radio-oncology centers delineated volumes according to ICRU 50 recommendations for one prostate and one head and neck case. In order to evaluate the consistency of the volume delineations, the following parameters were determined: 1) the target volumes (GTV, CTV and manually expanded PTV) and their extensions in the three main axes and 2) the correlation of the volume delineated by each pair of centers using the ratio of the intersection to the union (called proximity index). RESULTS: The delineated prostate volume was 105+/-55cm(3) for the CTV and 218+/-44cm(3) for the PTV. The delineated head and neck volume was 46+/-15cm(3) for the GTV, 327+/-154cm(3) for the CTV and 528+/-106cm(3) for the PTV. The mean proximity index for the prostate case was 0.50+/-0.13 for the CTV and 0.57+/-0.11 for the PTV. The proximity index for the head and neck case was 0.45+/-0.09 for the GTV, 0.42+/-0.13 for the CTV and 0.59+/-0.06 for the PTV. CONCLUSIONS: Large discrepancies between all the delineated target volumes were observed. There was an inverse relationship between the CTV volume and the margin between CTV and PTV, leading to less discrepancies in the PTV than is the CTV delineations. There was more spread in the sagittal and frontal planes due to CT pixel anisotropy, which suggests that radiation oncologists should delineate the target volumes not only in the transverse plane, but also in the sagittal and frontal planes to improve the delineation by allowing a consistency check.

3D seismic facies characterization and geological patterns recognition (Australian North West Shelf)

EXECUTIVE SUMMARY This PhD research, funded by the Swiss Sciences Foundation, is principally devoted to enhance the recognition, the visualisation and the characterization of geobodies through innovative 3D seismic approaches. A series of case studies from the Australian North West Shelf ensures the development of reproducible integrated 3D workflows and gives new insight into local and regional stratigraphic as well as structural issues. This project was initiated in year 2000 at the Geology and Palaeontology Institute of the University of Lausanne (Switzerland). Several collaborations ensured the improvement of technical approaches as well as the assessment of geological models. - Investigations into the Timor Sea structural style were carried out at the Tectonics Special Research Centre of the University of Western Australia and in collaboration with Woodside Energy in Perth. - Seismic analysis and attributes classification approach were initiated with Schlumberger Oilfield Australia in Perth; assessments and enhancements of the integrated seismic approaches benefited from collaborations with scientists from Schlumberger Stavanger Research (Norway). Adapting and refining from "linear" exploration techniques, a conceptual "helical" 3D seismic approach has been developed. In order to investigate specific geological issues this approach, integrating seismic attributes and visualisation tools, has been refined and adjusted leading to the development of two specific workflows: - A stratigraphic workflow focused on the recognition of geobodies and the characterization of depositional systems. Additionally, it can support the modelling of the subsidence and incidentally the constraint of the hydrocarbon maturity of a given area. - A structural workflow used to quickly and accurately define major and secondary fault systems. The integration of the 3D structural interpretation results ensures the analysis of the fault networks kinematics which can affect hydrocarbon trapping mechanisms. The application of these integrated workflows brings new insight into two complex settings on the Australian North West Shelf and ensures the definition of astonishing stratigraphic and structural outcomes. The stratigraphic workflow ensures the 3D characterization of the Late Palaeozoic glacial depositional system on the Mermaid Nose (Dampier Subbasin, Northern Carnarvon Basin) that presents similarities with the glacial facies along the Neotethys margin up to Oman (chapter 3.1). A subsidence model reveals the Phanerozoic geodynamic evolution of this area (chapter 3.2) and emphasizes two distinct mode of regional extension for the Palaeozoic (Neotethys opening) and Mesozoic (abyssal plains opening). The structural workflow is used for the definition of the structural evolution of the Laminaria High area (Bonaparte Basin). Following a regional structural characterization of the Timor Sea (chapter 4.1), a thorough analysis of the Mesozoic fault architecture reveals a local rotation of the stress field and the development of reverse structures (flower structures) in extensional setting, that form potential hydrocarbon traps (chapter 4.2). The definition of the complex Neogene structural architecture associated with the fault kinematic analysis and a plate flexure model (chapter 4.3) suggest that the Miocene to Pleistocene reactivation phases recorded at the Laminaria High most probably result from the oblique normal reactivation of the underlying Mesozoic fault planes. This episode is associated with the deformation of the subducting Australian plate. Based on these results three papers were published in international journals and two additional publications will be submitted. Additionally this research led to several communications in international conferences. Although the different workflows presented in this research have been primarily developed and used for the analysis of specific stratigraphic and structural geobodies on the Australian North West Shelf, similar integrated 3D seismic approaches will have applications to hydrocarbon exploration and production phases; for instance increasing the recognition of potential source rocks, secondary migration pathways, additional traps or reservoir breaching mechanisms. The new elements brought by this research further highlight that 3D seismic data contains a tremendous amount of hidden geological information waiting to be revealed and that will undoubtedly bring new insight into depositional systems, structural evolution and geohistory of the areas reputed being explored and constrained and other yet to be constrained. The further development of 3D texture attributes highlighting specific features of the seismic signal, the integration of quantitative analysis for stratigraphic and structural processes, the automation of the interpretation workflow as well as the formal definition of "seismo-morphologic" characteristics of a wide range of geobodies from various environments would represent challenging examples of continuation of this present research. The 21st century will most probably represent a transition period between fossil and other alternative energies. The next generation of seismic interpreters prospecting for hydrocarbon will undoubtedly face new challenges mostly due to the shortage of obvious and easy targets. They will probably have to keep on integrating techniques and geological processes in order to further capitalise the seismic data for new potentials definition. Imagination and creativity will most certainly be among the most important quality required from such geoscientists.

Partial melting of basic dykes in the contact aureoles of two shallow-level, gabbro-pyroxenite intrusions, Fuerteventura

THESIS ABSTRACT Low-pressure anatexis of basic dykes gave rise to unusual, zebra-like migmatites, in the contact metamorphic aureoles of two layered gabbro-pyroxenite intrusions, PXl and PX2, in the root zone of an ocean island, Fuerteventura Basal Complex (Canary Islands). This thesis focuses on the understanding of processes attributing to the partial melting and formation of these migmatites, characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations with perfectly preserved igneous textures. The presence of fluids are required to decrease the solidus of basic igneous lithologies, to allow partial melting in such aloes-pressure (1-2 kb) environment. An oxygen isotope study was thus carried out on dykes inside and beyond the PX2 aureole, in order to decipher the nature and origin of such fluids. Low or negative δ18O values were obtained for whole rocks and mineral-separates, decreasing towards the contact, with the intrusion itself retaining fairly high values. This trend has been attributed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation. A comparison of whole rock and mineral major- and trace- element data allowed the redistribution of elements .between different mineral phases and generations, during contact metamorphism and partial melting to be assessed. Certain trace-elements, e.g. Zr, Hf, Y, and REEs, were internally redistributed during contact metamorphic recrystallisation, causing- the enrichment of neocrystallised diopsides compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post- solidus reequilibration on cooling of the migmatite terrain. The mineralogical, textural and geochemical evolution of dykes in a contact metamorphic aureole, is recorded around an apophysis of the PX1 intrusion, where there is evidence of incipient partial melting. Hydrothermal mineral pseudomorphs in the outer parts of the aureole are progressively replaced by dry mineral assemblages, with increasingly recrystallised diopside and evidence of partial melting -the extent of which varies from one lithology to another. The appearance of more mafic lithologies towards the intrusion, with lower whole rock SiO2 and mobile element abundances, e.g. Rb, Cs, K, has been explained by the migration and accumulation of feldspathic material into leucosomes outside the samples. A micro-structural study of leucosomes and leucocratic pods, with the aid of high-resolution X-ray computed micro-tomography (HRXµCT), allowing the visualization and quantification of shapes and orientations, was carried out in order to better understand the processes of melt segregation in the PX1 aureole. Leucocratic pods, representing former amygdales, are considered as natural strain ellipsoids. Their short axes are oriented perpendicular to leucosome planes, which sub-parallel the intrusive contact. Leucosomes thus effectively represent foliation planes. This implies that the direction of maximum shortening, during migmatisation, was perpendicular to the orientation of leucosomes, contradicting earlier models that suggest leucosomes represent tension veins. RESUME DE LA THESE Un phénomène rare de fusion partielle de filons basiques à basse pression a été étudié dans les auréoles de contact de deux intrusions litées de gabbro-pyroxénite, PX1 et PX2, localisées dans le soubassement de l'île volcanique de Fuerteventura aux Canaries. Cette anatexie a engendré des migmatites finement zébrées d'aspect très inhabituel, dont les processus de formation ont été étudiés dans le présent travail. Ces roches sont caractérisées par un réseau dense de veinules leucocrates d'épaisseur millimétrique, dont les textures ignées sont parfaitement préservées. La fusion partielle de roches basiques à basse pression (1-2 kbar) requiert la présence d'eau afin d'abaisser le solidus du système à des températures géologiquement réalistes. Une étude comparative des isotopes de l'oxygène a ainsi été menée sur des filons respectivement affectés et non affectés par le métamorphisme de contact, afin de confirmer la présence de ces fluides, de déterminer l'importance de leur interaction avec les roches et leur origine. Des valeurs de δ180 basses ou négatives ont été mesurées sur roche totale et minéraux séparés, décroissantes en direction du contact, alors que l'intrusion elle-même a conservé des valeurs élevées. Ce gradient a été attribué à l'advection d'eau météorique durant la mise en place du magma, les températures les plus élevées favorisant d'autant plus la circulation des fluides et les échanges isotopiques durant la recristallisation des roches. Cette recristallisation engendré une redistribution chimique complète des éléments entre les différentes générations de minéraux résultant du métamorphisme de contact et de l'anatexie, mise en évidence par microanalyse. Certains éléments traces comme Zr, Hf, Y et les REE ont été concentrés dans le diopside néoformé consécutivement à la déstabilisation de minéraux primaires riches en ces éléments comme la kaersutite ou le sphène. Les compositions en éléments majeurs et traces des minéraux des mélanosomes et leucosomes des migmatites sont pratiquement identiques, indiquant une rééquilibration syn- à postsolidus lors du refroidissement de l'auréole de contact. La transformation progressive des filons basiques au niveau de leur minéralogie, textures et composition chimique a pu être observée en détail à l'approche du contact d'une apophyse de l'intrusion PX1. La paragenèse magmatique initiale n'est jamais préservée, les faciès les plus distants du contact étant constitués d'un assemblage pseudomorphique hydrothermal. Ce dernier est progressivement remplacé par des assemblages anhydres incluant du diopside néoformé, puis apparaissent les premiers signes de fusion partielle, dont l'importance varie fortement d'une lithologie à l'autre. L'apparition de faciès plus basiques en direction du contact, avec des teneurs réduites en SiO2 et en éléments incompatibles tels Rb, Cs, K, a été attribuée à l'échappement de leucosomes feldspathiques hors du système. Une étude microstructurale de la distribution spatiale du matériel leucocrate au sein des migmatites par microtomographie X de haute résolution (HRXµCT) a été menée pour mieux comprendre les processus de ségrégation des liquides dans l'auréole de PX1. De petites entités ovoïdes, représentant d'anciennes structures amygdalaires au sein des filons, ont été considérées comme des ellipsoïdes marqueurs de la déformation finie. Leur petit axe est orienté perpendiculairement aux plans définis par les leucosomes, eux-mêmes subparallèles au contact intrusif. Les leucosomes matérialisent donc des plans de clivage. Ainsi, la direction de raccourcissement maximum durant la fusion partielle était perpendiculaire à l'orientation des leucosomes, contrairement à ce qui a été dit dans de précédentes publications, qui suggéraient que les leucosomes représentaient des veines de tension. RESUME DE LA THESE (POUR LE GRAND PUBLIC) L'observation directe du soubassement d'une île volcanique est une occasion rare, accessible dans le «complexe de base » de l'île canarienne de Fuerteventura. Ce dernier a enregistré divers phénomènes magmatiques, métamorphiques et de fusion partielle induits par l'intrusion répétée de magmas alimentant des appareils volcaniques sus jacents, sous forme de petits plutons, essaims de filons et complexes annulaires de gabbros alcalins, pyroxénites, syénites et carbonatites. Dans ce contexte de flux de chaleur élevé, des filons basiques ont subi une fusion partielle au contact de deux intrusions de gabbro-pyroxénite, un phénomène extrêmement rare à une profondeur aussi réduite, estimée à quelque 3-6 km. Les produits de cette fusion partielle sont des liquides très riches en feldspath, concentrés en un réseau dense de veinules blanches (leucosomes) au sein du matériau résiduel sombre non fondu (mélanosome) pour former ce qu'on appelle des migmatites. Outre les aspects pétrologiques liés à la formation de ces migmatites, l'intérêt majeur du phénomène réside dans le fait qu'il puisse représenter la source des magmas évolués parfois observés sur les îles océaniques. A des pressions aussi faibles que dans le soubassement de Fuerteventura (1-2 Kbar), la présence de fluides abondants est nécessaire pour abaisser la température de début de fusion des roches (solidus) à des valeurs géologiquement réalistes. Des expériences ont montré que même en présence de plusieurs %-poids d'eau, une température de 1000°C était encore nécessaire pour obtenir une proportion de liquide équivalente à celle observée sur le terrain, soit 25%. Or les magmas alcalins des îles océaniques, bien qu'hydratés, n'en contiennent de loin pas autant, ce qui implique une source d'eau externe. Une étude isotopique de l'oxygène a été entreprise afin de tester cette hypothèse. Les valeurs obtenues en 5180 sont basses ou négatives et indiquent l'influence d'eau d'origine météorique. Cette eau de pluie se serait infiltrée le long des filons depuis la surface du volcan et les aurait complètement hydrothermalisés en profondeur (situation encore visible à l'extérieur de l'auréole de contact), leur permettant ainsi de stocker l'eau nécessaire à leur fusion partielle ultérieure. L'interaction entre eau de pluie et filons a été d'autant plus importante que ces derniers étaient proches du contact avec l'intrusion, ce qui suggère que la circulation de ces eaux et leur interaction avec les roches a été favorisée par la chaleur fournie par l'intrusion elle-même. Un autre aspect de ce travail s'est focalisé sur la redistribution des éléments traces au sein des minéraux des filons basiques durant le métamorphisme de contact et la fusion partielle. Ainsi, le pyroxène de seconde génération est-il sensiblement enrichi en traces telles Zr, Hf, Y et les terres rares, par rapport au pyroxène magmatique originel, en relation avec la déstabilisation de minéraux primaires riches en ces éléments tels le sphène et la kaersutite. Cependant, les compositions en éléments majeurs et traces des minéraux recristallisés des migmatites sont pratiquement les mêmes dans les leucosomes et les mélanosomes, suggérant une rééquilibration chimique complète durant le refroidissement de ces lithologies. Si certaines migmatites se sont comportées en système chimiquement fermé (hormis l'eau météorique), d'autres filons ont manifestement perdu une partie de leurs leucosomes, ainsi qu'en témoigne leur composition progressivement appauvrie en silice et autres éléments incompatibles mobiles, tels K et Rb à l'approche du contact de l'intrusion. Parallèlement à cette évolution chimique, les paragenèses hydrothermales distantes du contact sont progressivement remplacées par des paragenèses anhydres, puis par l'apparition des premiers leucosomes, tandis que les textures magmatiques initiales sont complètement effacées au profit d'une combinaison de textures magmatique dans les leucosomes et en mosaïque dans les mélanosomes. Enfin, la distribution spatiale des liquides de fusion partielle a été étudiée par microtomographie X de haute résolution, dans des filons contenant des entités ovoïdes leucocrates, sans doute d'anciennes amygdales à zéolites. Ces dernières ont été considérées comme des ellipsoïdes de la déformation finie. L'orientation de leur petit axe, perpendiculaire au plan défini par les veinules de leucosomes, indique que ces derniers représentent des plans de clivage perpendiculaires à la direction de raccourcissement maximum. Ainsi, la ségrégation des liquides de fusion partielle se serait faite dans les plans de compression et non dans des plans de dilatation, contrairement à ce que laisserait penser le sens commun.