VALIDATION OF A 'DISPLACEMENT TOMOGRAPHY' INVERSION METHOD FOR MODELING SHEET INTRUSIONS

The study of volcano deformation data can provide information on magma processes and help assess the potential for future eruptions. In employing inverse deformation modeling on these data, we attempt to characterize the geometry, location and volume/pressure change of a deformation source. Techniques currently used to model sheet intrusions (e.g., dikes and sills) often require significant a priori assumptions about source geometry and can require testing a large number of parameters. Moreover, surface deformations are a non-linear function of the source geometry and location. This requires the use of Monte Carlo inversion techniques which leads to long computation times. Recently, ‘displacement tomography’ models have been used to characterize magma reservoirs by inverting source deformation data for volume changes using a grid of point sources in the subsurface. The computations involved in these models are less intensive as no assumptions are made on the source geometry and location, and the relationship between the point sources and the surface deformation is linear. In this project, seeking a less computationally intensive technique for fracture sources, we tested if this displacement tomography method for reservoirs could be used for sheet intrusions. We began by simulating the opening of three synthetic dikes of known geometry and location using an established deformation model for fracture sources. We then sought to reproduce the displacements and volume changes undergone by the fractures using the sources employed in the tomography methodology. Results of this validation indicate the volumetric point sources are not appropriate for locating fracture sources, however they may provide useful qualitative information on volume changes occurring in the surrounding rock, and therefore indirectly indicate the source location.

Reassessing the Evolutionary History of the 17q21 Inversion Polymorphism

A polymorphic inversion that lies on chromosome 17q21 comprises two major haplotype families (H1 and H2) that not only differ in orientation but also in copy-number. Although the processes driving the spread of the inversion-associated lineage (H2) in humans remain unclear, a selective advantage has been proposed for one of its subtypes. Here, we genotyped a large panel of individuals from previously overlooked populations using a custom array with a unique panel of H2-specific single nucleotide polymorphisms and found a patchy distribution of H2 haplotypes in Africa, with North Africans displaying a higher frequency of inverted subtypes, when compared with Sub-Saharan groups. Interestingly, North African H2s were found to be closer to "non-African" chromosomes further supporting that these populations may have diverged more recently from groups outside Africa. Our results uncovered higher diversity within the H2 family than previously described, weakening the hypothesis of a strong selective sweep on all inverted chromosomes and suggesting a rather complex evolutionary history at this locus.

A qualitative impairment in face perception in Alzheimer’s disease : evidence from a reduced face inversion effect

Prevalent face recognition difficulties in Alzheimer’s disease (AD) have typically been attributed to the underlying episodic and semantic memory impairment. The aim of the current study was to determine if AD patients are also impaired at the perceptual level for faces, more specifically at extracting a visual representation of an individual face. To address this question, we investigated the matching of simultaneously presented individual faces and of other nonface familiar shapes (cars), at both upright and inverted orientation, in a group of mild AD patients and in a group of healthy older controls matched for age and education. AD patients showed a reduced inversion effect (i.e., larger performance for upright than inverted stimuli) for faces, but not for cars, both in terms of error rates and response times. While healthy participants showed a much larger decrease in performance for faces than for cars with inversion, the inversion effect did not differ significantly for faces and cars in AD. This abnormal inversion effect for faces was observed in a large subset of individual patients with AD. These results suggest that AD patients have deficits in higher-level visual processes, more specifically at perceiving individual faces, a function that relies on holistic representations specific to upright face stimuli. These deficits, combined with their memory impairment, may contribute to the difficulties in recognizing familiar people that are often reported in patients suffering from the disease and by their caregivers.

A Trans-dimensional inversion algorithm to model deformation sources with unconstrained shape in finite element domains

Ground deformation provides valuable insights on subsurface processes with pattens reflecting the characteristics of the source at depth. In active volcanic sites displacements can be observed in unrest phases; therefore, a correct interpretation is essential to assess the hazard potential. Inverse modeling is employed to obtain quantitative estimates of parameters describing the source. However, despite the robustness of the available approaches, a realistic imaging of these reservoirs is still challenging. While analytical models return quick but simplistic results, assuming an isotropic and elastic crust, more sophisticated numerical models, accounting for the effects of topographic loads, crust inelasticity and structural discontinuities, require much higher computational effort and information about the crust rheology may be challenging to infer. All these approaches are based on a-priori source shape constraints, influencing the solution reliability. In this thesis, we present a new approach aimed at overcoming the aforementioned limitations, modeling sources free of a-priori shape constraints with the advantages of FEM simulations, but with a cost-efficient procedure. The source is represented as an assembly of elementary units, consisting in cubic elements of a regular FE mesh loaded with a unitary stress tensors. The surface response due to each of the six stress tensor components is computed and linearly combined to obtain the total displacement field. In this way, the source can assume potentially any shape. Our tests prove the equivalence of the deformation fields due to our assembly and that of corresponding cavities with uniform boundary pressure. Our ability to simulate pressurized cavities in a continuum domain permits to pre-compute surface responses, avoiding remeshing. A Bayesian trans-dimensional inversion algorithm implementing this strategy is developed. 3D Voronoi cells are used to sample the model domain, selecting the elementary units contributing to the source solution and those remaining inactive as part of the crust.

Synthetic generation of amyloid PET images by non-linear dimensionality reduction inversion

Privacy issues and data scarcity in PET field call for efficient methods to expand datasets via synthetic generation of new data that cannot be traced back to real patients and that are also realistic. In this thesis, machine learning techniques were applied to 1001 amyloid-beta PET images, which had undergone a diagnosis of Alzheimer’s disease: the evaluations were 540 positive, 457 negative and 4 unknown. Isomap algorithm was used as a manifold learning method to reduce the dimensions of the PET dataset; a numerical scale-free interpolation method was applied to invert the dimensionality reduction map. The interpolant was tested on the PET images via LOOCV, where the removed images were compared with the reconstructed ones with the mean SSIM index (MSSIM = 0.76 ± 0.06). The effectiveness of this measure is questioned, since it indicated slightly higher performance for a method of comparison using PCA (MSSIM = 0.79 ± 0.06), which gave clearly poor quality reconstructed images with respect to those recovered by the numerical inverse mapping. Ten synthetic PET images were generated and, after having been mixed with ten originals, were sent to a team of clinicians for the visual assessment of their realism; no significant agreements were found either between clinicians and the true image labels or among the clinicians, meaning that original and synthetic images were indistinguishable. The future perspective of this thesis points to the improvement of the amyloid-beta PET research field by increasing available data, overcoming the constraints of data acquisition and privacy issues. Potential improvements can be achieved via refinements of the manifold learning and the inverse mapping stages during the PET image analysis, by exploring different combinations in the choice of algorithm parameters and by applying other non-linear dimensionality reduction algorithms. A final prospect of this work is the search for new methods to assess image reconstruction quality.

Observational evidences on the modulation of the South American Low Level Jet east of the Andes according the ENSO variability

The differences on the phase and wavelength of the quasi-stationary waves over the South America generated by El Nino (EN) and La Nina (LN) events seem to affect the daily evolution of the South American Low Level Jet east of the Andes (SALLJ). For the austral summer period of 1977 2004 the SALLJ episodes detected according to Bonner criterion 1 show normal to above-normal frequency in EN years, and in LN years the episodes show normal to below-normal frequency. During EN and LN years the SALLJ episodes were associated with positive rainfall anomalies over the La Plata Basin, but more intense during LN years. During EN years the increase in the SALLJ cases were associated to intensification of the Subtropical Jet (SJ) around 30 degrees S and positive Sea Level Pressure (SLP) anomalies over the western equatorial Atlantic and tropical South America, particularly over central Brazil. This favored the intensification of the northeasterly trade winds over the northern continent and it channeled by the Andes mountain to the La Plata Basin region where negative SLP are found. The SALLJ cases identified during the LN events were weaker and less frequent when compared to those for EN years. In this case the SJ was weaker than in EN years and the negative SLP anomalies over the tropical continent contributed to the inversion of the northeasterly trade winds. Also a southerly flow anomaly was generated by the geostrophic balance due to the anomalous blocking over southeast Pacific and the intense cyclonic transient over the southern tip of South America. As result the warm tropical air brought by the SALLJ encounters the cold extratropical air from the southerly winds over the La Plata basin. This configuration can increase the conditional instability over the La Plata basin and may explain the more intense positive rainfall anomalies in SALLJ cases during LN years than in EN years.

Depositional environments for strata cored in CRP-3 (Cape Roberts Project), Victoria Land Basin, Antarctica: palaeoglaciological and palaeoclimatological inferences

Cape Roberts Project drill core 3 (CRP-3) was obtained from Roberts ridge, a sea-floor high located at 77°S, 12 km offshore from Cape Roberts in western McMurdo Sound, Antarctica. The recovered core is about 939 m long and comprises strata dated as being early Oligocene (possibly latest Eocene) in age, resting unconformably on ∼ 116 m of basement rocks consisting of Palaeozoic Beacon Supergroup sediments. The core includes ten facies commonly occuring in five major associations that are repeated in particular sequences throughout the core and which are interpreted as representing different depositional environments through time. Depositional systems inferred to be represented in the succession include: outer shelf, inner shelf, nearshore to shoreface each under iceberg influence, deltaic and/or grounding-line fan, and ice proximal-ice marginal-subglacial (mass flow/rainout diamictite/subglacial till) singly or in combination. The record is taken to represent the initial talus/alluvial fan setting of a glaciated rift margin adjacent to the block-uplifted Transantarctic Mountains. Development of a deltaic succession upcore was probably associated with the formation of palaeo-Mackay valley with temperate glaciers in its headwaters. At that stage glaciation was intense enough to support glaciers ending in the sea elsewhere along the coast, but a local glacier was fluctuating down to the sea by the time the youngest part of CRP-3 was being deposited. Changes in palaeoenvironmental interpretations in this youngest part of the core are used to estimate relative glacial proximity to the drillsite through time. These inferred glacial fluctuations are compared with the global δ18O and Mg/Ca curves to evaluate the potential of glacial fluctuations on Antarctica for influencing these records of global change. Although the comparisons are tentative at present, the records do have similarities, but there are also some differences that require further evaluation.

Response of a multi-domain continental margin to compression: Study from seismic reflection-refraction and numerical modelling in the Tagus Abyssal Plain

The effects of the Miocene through Present compression in the Tagus Abyssal Plain are mapped using the most up to date available to scientific community multi-channel seismic reflection and refraction data. Correlation of the rift basin fault pattern with the deep crustal structure is presented along seismic line IAM-5. Four structural domains were recognized. In the oceanic realm mild deformation concentrates in Domain I adjacent to the Tore-Madeira Rise. Domain 2 is characterized by the absence of shortening structures, except near the ocean-continent transition (OCT), implying that Miocene deformation did not propagate into the Abyssal Plain, In Domain 3 we distinguish three sub-domains: Sub-domain 3A which coincides with the OCT, Sub-domain 3B which is a highly deformed adjacent continental segment, and Sub-domain 3C. The Miocene tectonic inversion is mainly accommodated in Domain 3 by oceanwards directed thrusting at the ocean-continent transition and continentwards on the continental slope. Domain 4 corresponds to the non-rifted continental margin where only minor extensional and shortening deformation structures are observed. Finite element numerical models address the response of the various domains to the Miocene compression, emphasizing the long-wavelength differential vertical movements and the role of possible rheologic contrasts. The concentration of the Miocene deformation in the transitional zone (TC), which is the addition of Sub-domain 3A and part of 3B, is a result of two main factors: (1) focusing of compression in an already stressed region due to plate curvature and sediment loading; and (2) theological weakening. We estimate that the frictional strength in the TC is reduced in 30% relative to the surrounding regions. A model of compressive deformation propagation by means of horizontal impingement of the middle continental crust rift wedge and horizontal shearing on serpentinized mantle in the oceanic realm is presented. This model is consistent with both the geological interpretation of seismic data and the results of numerical modelling.

A reconnaissance study of Upper Jurassic sediments of the Lusitanian Basin

The Upper Jurassic evolution of the Lusitanian Basin is shown to be linked to the rifting phase which preceded the separation of Iberia and the Grand Banks. Structural controls on sedimentation include both NNE-SSW trending faults in the Hercynian basement, and contemporaneous movement of salt diapirs. At the beginning of Upper Oxfordian times, the entire basin had been levelled to within a few metres of sea level, so that the freshwater algal marsh and marginal marine facies of the Cabaços and Vale Verde Beds rest on Triassic to Callovian strata. In the latter part of the Upper Oxfordian. carbonate sedimentation continued, with fluctuating salinity lagoons in the north (Pholodomya protei Beds) separated from shallow open marine carbonates in the south (Montejunto Beds) by the Caldas da Rainha diapir-barrier island complex. The commencement of rifting is recorded in the Kimmeridgian by the sudden influx of terrigenous clastics (developed in both fluviatile and deltaic/submarine fan environments) and accelerated depositional rates in excess of 10cm/10 k.yrs in association with contemporaneous faulting along the SE margin of the Arruda sub-basin. The Caldas-Santa Cruz chain of diapiric structures continued to influence the distribution of carbonate and clastic sediments. In the Portlandian, a simpler facies pattern occurs, with fluviatile clastics interfingering to the south with shallow low energy carbonates.

Trypanites ichnofacies: palaeoenvironmental and tectonic implications. A case study from the miocene disconformity at Foz da Fonte (lower Tagus basin, Portugal)

Palaeogeography, Palaeoclimatology, Palaeoecology 292, 35–43

Paleoalterações e cimentações nos depósitos continentais terciários de Portugal central: importância na interpretação de processos antigos

This paper describes the palaeoweathering, cementation, clay minerals association and other closely related characteristics of central Portugal allostratigraphic Tertiary units (SLD's), that can be used for palaeoclimatic interpretation and palaeoenvironmental reconstruction. Lateral and vertical changes in palaeosols are of value for improving our understanding of the autocyclic and allocyclic controls on sediment acumulation in an alluvial basin, but they can also have stratigraphic importance. In some cases it is concluded that the geomorphological setting may have been more decisive than climatic conditions to the production of the palaeoweathering. During late Palaeogene (SLD7-8), surface and near-surface silicification were developed on tectonically stable land surfaces of minimal local relief under a semi-arid climate; groundwater flow was responsible for some eodiagenesis calcareous accumulations, with the neoformation of palygorskite. Conditions during the Miocene (SLD9-11) were favourable for the smectization of the metamorphic basement and arenization of granites. Intense rubefaction associated with basement conversion into clay (illite and kaolinite), is ascribed to internal drainage during late Messinian-Zanclean (SLD12). During Piacenzian (SLD13) intense kaolinization and hydromorphism are typical, reflecting a more humid and hot temperate climate and important Atlantic fluvial drainage. Later on (Gelasian-early Pleistocene ?; SLD14). more cold and dry conditicns are interpreted, at the beginning of the fluvial incision sage. Silica cementation is identified in the upper Eocence-Oligocene ? (SLD18; the major period of silicification), middle to upper Miocene (SLD10)and upper Tortonian-Messinian (SLD11); these occurrences are compatible with either arid or semi-arid conditions and the establishment of a flat landscape upon which a silcrete was developed.

Tectónica da região de Sintra

Photo-interpretation of aerial stereopairs of the Sintra region on the approx. 1/32 000 scale together with field work allowed the production of the present Tectonic Map of the Sintra region. It is now possible to separate structures which resulted from two different tectonic events: one, corresponding to the intrusion of the Late Cretaceous Sintra igneous diapir, and the other the Miocene compressive event, the most important tectonic inversion phase of the Lusitanian Basin. The former are present to the south, southeast and east of the intrusion and within the intrusion itself, affecting the peripheral granites and their contacts with the gabbro-syenite core. These structures comprehend: i) faults and conical fractures striking parallel to the massif boundary, which were intruded by dykes, ii) vertical faults and fractures of two conjugate sets, dextral NNW-SSE and sinistral NNE-SSW. These faults are certainly associated with the E-W striking massif's northwards directed thrust and indicate a N-S oriented horizontal maximum compressive stress. The Miocene compressive event reactivated most of the inherited structures as follows. The NNWSSE faults located on the Sintra southern platform were reactivated as dextral strike slip faults and the E-W thrust along the northern boundary of the massif was also reactivated. This thrust propagated to the east. It also enhanced the asymmetry of the rim-syncline, uplifted the massif and reactivated the NNE-SSW faults as sinistral lateral ramps, which also accommodated vertical throw. The present Tectonic Map of Sintra together with the available geophysical data (MOREIRA, 1984, KULLBERG et al., 1991, SILVA & MIRANDA, 1994) allowed reassessment of the models proposed for the emplacement of the Sintra, Sines and Monchique igneous massifs, which intruded during Late Cretaceous times along the deep dextral NNW-SSE oriented strike slip fault (RIBEIRO et al., 1979; TERRINHA, 1998; TERRINHA & KULLBERG, 1998).

Reinterpretation de l'Histoire Stratigraphique et Tectonique de la Structure de Serra de El-Rei (Portugal)

Les travaux faits sur le terrain et les interprétations cartographiques de ces dernières années, montrant la pertinence des observations de Paul Choffat au siècle dernier, permettent de préciser l'échelle stratigraphique et de mieux distinguer les phases de l'histoire tectonique de la structure de Serra de EI-Rei. En attendant la publication d'une cartographie très détaillée mettant en évidence le détail des déformations tectoniques, il faut souligner la dissymétrie des deux flancs de cette structure ainsi que le rôle de certaines failles transverses. Sur son flanc nord, la série liasique, allant des dolomies du Lias inférieur au Toarcien inférieur ou moyen, est recouverte stratigraphiquement et en légère discordance par le "Lusitanien" qui remanie à sa base des éléments du Toarcien inférieur ou moyen. Sur son flanc sud, dans la partie centrale, la série stratigraphique liasique est comparable à celle de l'autre flanc, mais les dépôts du Toarcien, plus complets, sont affectés au cours de la zone à Meneghinii par des perturbations synsédimentaires avec fossiles brisés et/ou remaniés et des réductions d'épaisseur. Plus à l'Est, le Toarcien est très réduit et (tectoniquement?) en contact avec le Dogger. Au Sud-Ouest en revanche, prés de Bolhos, la série toarcienne est plus complète et couronnée par l'Aalénien inférieur, mais tronquée par un contact oblique N-S, repris par la tectonique récente. Sur ce flanc, le Bajocien inférieur est représenté mais peu épais; il est suivi par le Bajocien supérieur, le Bathonien et le Callovien très développés qui sont recouverts par un "Lusitanien" épais, très différent de celui du flanc nord. Cette série du Jurassique moyen et supérieur a été décrite dans la thèse de Christiane Ruget (1963). La structure de Serra de EI-Rei constitue une sorte de raccord E-W entre les importantes structures (diapiriques ?) méridiennes de Caldas da Rainha au NE et de Bolhos au SW. II s'agit d'une structure de tectonique distensive, en demi-graben, d'âge jurassique, qui est responsable de la dissymétrie des séries lithostratigraphiques sur ses deux flancs; au Miocène, pendant la compression bétique, cette structure a subi une inversion en constriction. Un modèle de l'évolution sédimentaire du secteur est présenté.

Time-variations of equivalent water heights'from Grace Mission and in-situ river stages in the Amazon basin

Gravity Recovery and Climate Experiment (GRACE) mission is dedicated to measuring temporal variations of the Earth's gravity field. In this study, the Stokes coefficients made available by Groupe de Recherche en Géodésie Spatiale (GRGS) at a 10-day interval were converted into equivalent water height (EWH) for a ~4-year period in the Amazon basin (from July-2002 to May-2006). The seasonal amplitudes of EWH signal are the largest on the surface of Earth and reach ~ 1250mm at that basin's center. Error budget represents ~130 mm of EWH, including formal errors on Stokes coefficient, leakage errors (12 ~ 21 mm) and spectrum truncation (10 ~ 15 mm). Comparison between in situ river level time series measured at 233 ground-based hydrometric stations (HS) in the Amazon basin and vertically-integrated EWH derived from GRACE is carried out in this paper. Although EWH and HS measure different water bodies, in most of the cases a high correlation (up to ~80%) is detected between the HS series and EWH series at the same site. This correlation allows adjusting linear relationships between in situ and GRACE-based series for the major tributaries of the Amazon river. The regression coefficients decrease from up to down stream along the rivers reaching the theoretical value 1 at the Amazon's mouth in the Atlantic Ocean. The variation of the regression coefficients versus the distance from estuary is analysed for the largest rivers in the basin. In a second step, a classification of the proportionality between in situ and GRACE time-series is proposed.

Analog Modelling of N-S anticlines of the External Sierras (Jaca basin, Southern Pyrenees)

Report for the scientific sojourn carried out at the Uppsala Universitet, Sweden, from April to July the 2007. Two series of analogue models are used to explore ductile-frictional contrasts of the basal décollement in the development of oblique and transverse structures simultaneously to thin-skinned shortening. These models simulate the evolution of the Central External Sierras (Southern Pyrenees, Spain), which constitute the frontal emerging part of the southernmost Pyrenean thrust sheet. They are characterized by the presence of transverse N-S to NW-SE anticlines, which are perpendicular to the Pyrenean structural trend and developed in the hangingwall of the Santo Domingo thrust system, detaching on an unevenly distributed Triassic materials (evaporitic-dolomitic interfingerings). Model setup performs a décollement made by three patches of silicone neighbouring pure brittle sand. Model series A test the thickness ratio between overburden and décollement. Model series B test the width of frictional detachment areas. Model results show how deformation reaches further in areas detached on ductile layer whereas frictional décollement areas assimilate the strain by means of an additional uplift. This replicates the structural style of Central External Sierras: higher structural relief of N-S anticlines with regard to orogen-parallel structures, absence of a representative ductile décollement in the core, tilting towards the orogen and foreland-side closure not thrusted by the frontal emerging South-Pyrenean thrust.