Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

Spectral-Element and Adjoint 3D Full-Wave Tomography for the Lithosphere of Central Italy

The primary objective of this thesis is to obtain a better understanding of the 3D velocity structure of the lithosphere in central Italy. To this end, I adopted the Spectral-Element Method to perform accurate numerical simulations of the complex wavefields generated by the 2009 Mw 6.3 L’Aquila event and by its foreshocks and aftershocks together with some additional events within our target region. For the mainshock, the source was represented by a finite fault and different models for central Italy, both 1D and 3D, were tested. Surface topography, attenuation and Moho discontinuity were also accounted for. Three-component synthetic waveforms were compared to the corresponding recorded data. The results of these analyses show that 3D models, including all the known structural heterogeneities in the region, are essential to accurately reproduce waveform propagation. They allow to capture features of the seismograms, mainly related to topography or to low wavespeed areas, and, combined with a finite fault model, result into a favorable match between data and synthetics for frequencies up to ~0.5 Hz. We also obtained peak ground velocity maps, that provide valuable information for seismic hazard assessment. The remaining differences between data and synthetics led us to take advantage of SEM combined with an adjoint method to iteratively improve the available 3D structure model for central Italy. A total of 63 events and 52 stations in the region were considered. We performed five iterations of the tomographic inversion, by calculating the misfit function gradient - necessary for the model update - from adjoint sensitivity kernels, constructed using only two simulations for each event. Our last updated model features a reduced traveltime misfit function and improved agreement between data and synthetics, although further iterations, as well as refined source solutions, are necessary to obtain a new reference 3D model for central Italy tomography.

Implications of regional fault distribution and kinematics for the uplift of rift flanks around the Rwenzori Mountains, Southwestern Uganda

Das Ziel dieser Arbeit bestand in der Untersuchung der Störungsverteilung und der Störungskinematik im Zusammenhang mit der Hebung der Riftschultern des Rwenzori Gebirges.rnDas Rwenzori Gebirge befindet sich im NNE-SSWbis N-S verlaufenden Albertine Rift, des nördlichsten Segments des westlichen Armes des Ostafrikanischen Grabensystems. Das Albertine Rift besteht aus Becken unterschiedlicher Höhe, die den Lake Albert, Lake Edward, Lake George und Lake Kivu enthalten. Der Rwenzori horst trennt die Becken des Lake Albert und des Lake Edward. Es erstreckt sich 120km in N-S Richtung, sowie 40-50km in E-W Richtung, der h¨ochste Punkt befindet sich 5111 ü. NN. Diese Studie untersucht einen Abschnitt des Rifts zwischen etwa 1°N und 0°30'S Breite sowie 29°30' und 30°30' östlicher Länge ersteckt. Auch die Feldarbeit konzentrierte sich auf dieses Gebiet.rnrnHauptzweck dieser Studie bestand darin, die folgende These auf ihre Richtigkeit zu überprüfen: ’Wenn es im Verlauf der Zeit tatsächlich zu wesentlichen Änderungen in der Störungskinematik kam, dann ist die starke Hebung der Riftflanken im Bereich der Rwenzoris nicht einfach durch Bewegung entlang der Graben-Hauptst¨orungen zu erklären. Vielmehr ist sie ein Resultat des Zusammenspiels mehrerer tektonische Prozesse, die das Spannungsfeld beeinflussen und dadurch Änderungen in der Kinematik hervorrufen.’ Dadurch konzentrierte sich die Studie in erster Linie auf die Störungsanalyse.rnrnDie Kenntnis regionaler Änderungen der Extensionsrichtung ist entscheidend für das Verständnis komplexer Riftsysteme wie dem Ostafrikanischen Graben. Daher bestand der Kern der Untersuchung in der Kartierung von Störungen und der Untersuchung der Störungskinematik. Die Aufnahme strukturgeologischer Daten konzentrierte sich auf die Ugandische Seite des Rifts, und Pal¨aospannungen wurden mit Hilfe von St¨orungsdaten durch Spannungsinversion rekonstruiert.rnDie unterschiedliche Orientierung spr¨oder Strukturen im Gelände, die geometrische Analyse der geologischen Strukturen sowie die Ergebnisse von Mikrostrukturen im Dünnschliff (Kapitel 4) weisen auf verschiedene Spannungsfelder hin, die auf mögliche Änderungen der Extensionsrichtung hinweisen. Die Resultate der Spannungsinversion sprechen für Ab-, Über- und Blattverschiebungen sowie für Schrägüberschiebungen (Kapitel 5). Aus der Orientierung der Abschiebungen gehen zwei verschiedene Extensionsrichtungen hervor: im Wesentlichen NW-SE Extension in fast allen Gebieten, sowie NNE-SSW Extension im östlichen Zentralbereich.rnAus der Analyse von Blattverschiebungen ergaben sich drei unterschiedliche Spannungszustände. Zum Einen NNW-SSE bis N-S Kompression in Verbindung mit ENE-WSW bzw E-W Extension wurde für die nördlichen und die zentralen Ruwenzoris ausgemacht. Ein zweiter Spannungszustand mit WNW-ESE Kompression/NNE-SSW Extension betraf die Zentralen Rwenzoris. Ein dritter Spannungszustand mit NNW-SSE Extension betraf den östlichen Zentralteil der Rwenzoris. Schrägüberschiebungen sind durch dazu schräge Achsen charakterisiert, die für N-S bis NNW-SSE Kompression sprechen und ausschließlich im östlichen Zentralabschnitt auftreten. Überschiebungen, die hauptsächlich in den zentralen und den östlichen Rwenzoris auftreten, sprechen für NE-SW orientierten σ2-Achsen und NW-SE Extension.rnrnEs konnten drei unterschiedliche Spannungseinflüsse identifiziert werden: auf die kollisionsbedingte Bildung eines Überschiebungssystem folgte intra-kratonische Kompression und schließlich extensionskontrollierte Riftbildung. Der Übergang zwischen den beiden letztgenannten Spannungszuständen erfolgte Schrittweise und erzeugte vermutlich lokal begrenzte Transpression und Transtension. Gegenw¨artig wird die Störungskinematik der Region durch ein tensiles Spannungsregime in NW-SE bis N-S Richtung bestimmt.rnrnLokale Spannungsvariationen werden dabei hauptsächlich durch die Interferenzrndes regionalen Spannungsfeldes mit lokalen Hauptst¨orungen verursacht. Weitere Faktoren die zu lokalen Veränderungen des Spannungsfeldes führen können sind unterschiedliche Hebungsgeschwindigkeiten, Blockrotation oder die Interaktion von Riftsegmenten. Um den Einfluß präexistenter Strukturen und anderer Bedingungen auf die Hebung der Rwenzoris zu ermitteln, wurde der Riftprozeß mit Hilfe eines analogen ’Sandbox’-Modells rekonstruiert (Kapitel 6). Da sich die Moho-Diskontinuität im Bereich des Arbeitsgebietes in einer Tiefe von 25 km befindet, aktive Störungen aber nur bis zu einer Tiefe von etwa 20 km beobachtet werden können (Koehn et al. 2008), wurden nur die oberen 25 km im Modell nachbebildet. Untersucht und mit Geländebeobachtungen verglichen wurden sowohl die Reihenfolge, in der Riftsegmente entstehen, als auch die Muster, die sich im Verlauf der Nukleierung und des Wachstums dieser Riftsegmente ausbilden. Das Hauptaugenmerk wurde auf die Entwicklung der beiden Subsegmente gelegt auf denen sich der Lake Albert bzw. der Lake Edward und der Lake George befinden, sowie auf das dazwischenliegende Rwenzori Gebirge. Das Ziel der Untersuchung bestand darin herauszufinden, in welcher Weise das südwärts propagierende Lake Albert-Subsegment mit dem sinistral versetzten nordwärts propagierenden Lake Edward/Lake George-Subsegment interagiert.rnrnVon besonderem Interesse war es, in welcherWeise die Strukturen innerhalb und außerhalb der Rwenzoris durch die Interaktion dieser Riftsegmente beeinflußt wurden. rnrnDrei verschiedene Versuchsreihen mit unterschiedlichen Randbedingungen wurden miteinander verglichen. Abhängig vom vorherrschenden Deformationstyp der Transferzone wurden die Reihen als ’Scherungs-dominiert’, ’Extensions-dominiert’ und als ’Rotations-dominiert’ charakterisiert. Die Beobachtung der 3-dimensionalen strukturellen Entwicklung der Riftsegmente wurde durch die Kombination von Modell-Aufsichten mit Profilschnitten ermöglicht. Von den drei genannten Versuchsreihen entwickelte die ’Rotationsdominierten’ Reihe einen rautenförmiger Block im Tranferbereich der beiden Riftsegmente, der sich um 5−20° im Uhrzeigersinn drehte. DieserWinkel liegt im Bereich des vermuteten Rotationswinkel des Rwenzori-Blocks (5°). Zusammengefasst untersuchen die Sandbox-Versuche den Einfluss präexistenter Strukturen und der Überlappung bzw. Überschneidung zweier interagierender Riftsegmente auf die Entwicklung des Riftsystems. Sie befassen sich darüber hinaus mit der Frage, welchen Einfluss Blockbildung und -rotation auf das lokale Stressfeld haben.

Studio di fattibilità di una nuova tipologia di argano per il sollevamento di carichi nel settore off-shore

Studio di una nuova tipologia di argano per il sollevamento di carichi nel settore off-shore. In particolare, un argano funzionalmente identico ad uno in fase di realizzazione in azienda è stato modellato al CAD 3D e verificato strutturalmente utilizzando Ansys Workbench, software adatto per analisi F.E.M.. Ottenuti i risultati, ne è stato fatto un confronto con quelli relativi all'argano esistente (argano di ‟Moho‟); nello specifico sono stati comparati gli stress ottenuti nelle varie parti dei due argani (tamburo, spalla, flange, alette di rinforzo, albero) in seguito all‟applicazione di 4 casi di carico differenti (E‟, E‟‟, P‟, P‟‟).

Magma Extraction from the Mantle Wedge at Convergent Margins Through Dikes: A Parametric Sensitivity Analysis

We address under what conditions a magma generated by partial melting at 100 km depth in the mantle wedge above a subduction zone can reach the crust in dikes before stalling. We also address under what conditions primitive basaltic magma (Mg # >60) can be delivered from this depth to the crust. We employ linear elastic fracture mechanics with magma solidification theory and perform a parametric sensitivity analysis. All dikes are initiated at a depth of 100 km in the thermal core of the wedge, and the Moho is fixed at 35 km depth. We consider a range of melt solidus temperatures (800-1100 degrees C), viscosities (10-100 Pa s), and densities (2400-2700 kg m(-3)). We also consider a range of host rock fracture toughness values (50-300 MPa m(1/2)) and dike lengths (2-5 km) and two thermal structures for the mantle wedge (1260 and 1400 degrees C at 100 km depth and 760 and 900 degrees C at 35 km depth). For the given parameter space, many dikes can reach the Moho in less than a few hundred hours, well within the time constraints provided by U series isotope disequilibria studies. Increasing the temperature in the mantle wedge, or increasing the dike length, allows additional dikes to propagate to the Moho. We conclude that some dikes with vertical lengths near their critical lengths and relatively high solidus temperatures will stall in the mantle before reaching the Moho, and these may be returned by corner flow to depths where they can melt under hydrous conditions. Thus, a chemical signature in arc lavas suggesting partial melting of slab basalts may be partly influenced by these recycled dikes. Alternatively, dikes with lengths well above their critical lengths can easily deliver primitive magmas to the crust, particularly if the mantle wedge is relatively hot. Dike transport remains a viable primary mechanism of magma ascent in convergent tectonic settings, but the potential for less rapid mechanisms making an important contribution increases as the mantle temperature at the Moho approaches the solidus temperature of the magma.

Ar-40/Ar-39 Evidence for Middle Proterozoic (1300-1500 Ma) Slow Cooling of the Southern Black Hills, South Dakota, Midcontinent, North America: Implications for Early Proterozoic P-T Evolution and Posttectonic Magmatism

Ar-40/Ar-39 total gas and plateau dates from muscovite and biotite in the southern Black Hills, South Dakota, provide evidence for a period of Middle Proterozoic slow cooling. Early Proterozoic (1600-1650 Ma) mica dates were obtained from metasedimentary rocks located in a synformal structure between the Harney Peak and Bear Mountain domes and also south of Bear Mountain. Metamorphic rocks from the dome areas and undeformed samples of the similar to 1710 Ma Harney Peak Granite (HPG) yield Middle Proterozoic mica dates (similar to 1270-1500 Ma). Two samples collected between the synform and Bear Mountain dome yield intermediate total gas mica dates of similar to 1550 Ma. We suggest two end-member interpretations to explain the map pattern of cooling ages: (1) subhorizontal slow cooling of an area which exhibits variation in mica Ar retention intervals or (2) mild folding of a Middle Proterozoic (similar to 1500 Ma) similar to 300 degrees C isotherm. According to the second interpretation, the preservation of older dates between the domes may reflect reactivation of a preexisting synformal structure (and downwarping of relatively cold rocks) during a period of approximately east-west contraction and slow uplift during the Middle Proterozoic. The mica data, together with hornblende data from the Black Hills published elsewhere, indicate that the ambient country-rock temperature at the 3-4 kbar depth of emplacement of the HPG was between 350 degrees C and 500 degrees C, suggesting that the average upper crustal geothermal gradient was 25 degrees-40 degrees C/km prior to intrusion. The thermochronologic data suggest HPG emplacement was followed by a similar to 200 m.y. period of stability and tectonic quiescence with little uplift. We propose that crust thickened during the Early Proterozoic was uplifted and erosionally(?) thinned prior to similar to 1710 Ma and that the HPG magma was emplaced into isostatically stable crust of relatively normal thickness. We speculate that uplift and crustal thinning prior to HPG intrusion was the result of differential thinning of the subcrustal lithosphere beneath the Black Hills. If so, this process would have also caused an increase in mantle heat flux across the Moho and triggered vapor-absent melting of biotite to produce the HPG magma. This scenario for posttectonic granite generation is supported, in part, by the fact that in the whole of the Black Hills, the HPG is spatially associated with the deepest exposed Early Proterozoic country rock.

A deep crust–mantle boundary in the asteroid 4 Vesta

The asteroid 4 Vesta was recently found to have two large impact craters near its south pole, exposing subsurface material. Modelling suggested that surface material in the northern hemisphere of Vesta came from a depth of about 20 kilometres, whereas the exposed southern material comes from a depth of 60 to 100 kilometres. Large amounts of olivine from the mantle were not seen, suggesting that the outer 100 kilometres or so is mainly igneous crust. Here we analyse the data on Vesta and conclude that the crust–mantle boundary (or Moho) is deeper than 80 kilometres.

Medidas sanitarias para prevenir el ingreso y producción de haloanisoles en bodegas

La contaminación del vino con aromas y/o sabores “fúngicos" o “a moho" es consecuencia de la presencia de ciertos compuestos conocidos como haloanisoles, producidos por hongos filamentosos. El objetivo del presente trabajo es realizar una revisión bibliográfica de la información disponible relacionada con los hongos filamentosos productores de haloanisoles; su permanencia, desarrollo y la bioquímica de las reacciones de síntesis de haloanisoles, a fin de poder establecer las medidas sanitarias preventivas correspondientes.

Major element and mineral composition of gabbroic xenoliths and megacrysts from Jeju Island, South Korea

Gabbroic xenoliths and diverse megacrysts (e.g., clinopyroxenes, amphiboles and plagioclases), which correspond to the lithology ranging from gabbro-norite to gabbro, occur in the Pleisto-Holocene alkali basalts from Jeju Island, South Korea. The gabbroic xenoliths consist primarily of moderate-K2O plagioclase, Ti-Al-rich clinopyroxene and CaO-rich orthopyroxene; additionally, TiO2-rich amphibole (kaersutite) and Ti-Fe oxides might or might not be present. The plagioclase is the most dominant phase (approx. 60-70 vol.%). The xenoliths and megacrysts provide evidence for the modal metasomatism of the lower continental crust by the mafic magmas during the Pleistocene. The coarse grain size (up to 5 mm), moderate Mg# [=100xMg/(Mg+Fe(total)) atomic ratio] of pyroxenes (70-77) and textural features (e.g., poikilitic) indicate that the gabbroic xenoliths are consistent with a cumulus origin. The clinopyroxenes from these xenoliths are enriched in REE with smooth convex-upward MREE patterns, which are expected for cumulus minerals formed from a melt enriched in incompatible trace elements. The strikingly similar major and trace element variations and the patterns of constituent minerals clearly indicate a genetic link between the gabbroic xenoliths (plus megacrysts) and the host basalt, indicating that the xenoliths belong to the Jeju Pleisto-Holocene magma system. On the basis of the textural features, the mineral equilibria and the major and trace element variations, the xenoliths appear to have crystallized from basaltic melts at the reservoir-roof environment within the lower crust (4-7 kbars) above the present Moho estimates beneath Jeju Island, where the xenoliths represent wall rocks. Following the consolidation of the xenolith lithologies, volatile- and incompatible element-enriched melt/fluid, as metasomatic agents, infiltrated through the grain boundaries and/or cracks and reacted with the preexisting anhydrous phases, which produced the metasomatic amphiboles. This volatile-enriched melt/fluid could have evolved from the initially anhydrous compositions to the volatile-saturated compositions by the active fractional crystallization in the Jeju Pleisto-Holocene magma system. This process was significant in that it was a relatively young event and played an important role in the formation of the hydrous minerals and the metasomatization of the lower continental crust, which is a plume-impacted area along the Asian continental margin. The major and trace element analyses of the mineral phases from the xenoliths were performed to define the principal geochemical characteristics of the crustal lithosphere segment represented by the studied xenoliths.

Magnetic and gravity measurements on North-South profiles across the Grosse Meteor Bank, North Atlantic

During the first section of the "Meteor" cruise No. 2 a profile was run from the Azores to the south across the flanks of the Mid-Atlantic Ridge with a chain of seamounts. The profile extended between the Cruiser (living) and the Hyeres seamounts, which, according to our soundings, form a connected massif, and across the centre of the Grosse Meteor Bank (30°N, 28.5 °W). These seamounts rise from a depth of more than 4000 m up till close to the surface of the sea forming there a large almost flat plateau. In the case of the Grosse Meteor Bank, this plateau has a N-S extension of approx. 30 nautical miles and an E-W extension of approx. 20 nautical miles and reaches a height of 275 m in water depth. The gravity measurements yielded a density of the topographic masses of 2.6 g/cm**3 for the Grosse Meteor Bank. Magnitude and shape of the measured free-air anomaly are very well shown in a model computation with this density. The theoretical gravity effects of the seismically detected swell of cristalline rock and of the Moho depression (mountain root) are not indicated by the observational data. It can, therefore, be assumed that the latter two neutralize each other. It seems, accordingly, that there is no local isostatic compensation of the topographic masses. Hence, the density of 2.6 g/cm**3 obtained would be about the true density of rock. In connection with the mean velocity of P waves (Aric et al., 1968) obtained by seismic refraction methods it must be concluded that the material of the 1200-4000 m thick surface layer of the Grosse Meteor Bank consists of consolidated sediments. This finding is supported by the total intensity of the Earth's magnetic field over the Grosse Meteor Bank. On the assumption of a homogeneous magnetization in the direction of the present Earth's field, the computed anomaly of the massif deviates considerably from the measured anomaly while the magnetic field of the seismically detected crystalline body is capable of interpreting the observed data. Deviating from the prevailing interpretation of the seamounts' plateau as a volcanic cone with submarine abrasion, the Grosse Meteor Bank and the seamounts in the vicinity are assumed to be of continental origin. The questions whether these seamounts submerged later on or whether the sealevel has risen subsequently are, therefore, largely nonexistent.

Geochemical composition of olivine, spinel and amphibole in ODP Hole 209-1271B

Dunite and gabbroic materials recovered from Hole 1271B, Ocean Drilling Program (ODP) Leg 209, were examined for mineral chemistry to understand melt flow and melt-mantle reactions in the shallowest upper mantle of the Mid-Atlantic Ridge near the 15°20' Fracture Zone. Hole 1271B was drilled to 103.8 meters below seafloor on the inner corner high along the south wall of the 15°20' Fracture Zone. The total length of core collected was 15.9 m (recovery = ~15%). The dominant rock type in Hole 1271B is dunite, followed by brown amphibole gabbro, olivine gabbro, and troctolite, along with minor amounts of harzburgite and olivine gabbronorite. A large proportion of the dunite is associated with gabbroic rocks in Hole 1271B, similar to those observed in the Mohorovicic (Moho) transition zone of the Oman ophiolite, indicating significant magmatic activity in this region near the 15°20' Fracture Zone. Olivine Fo content varies from 89.2 to 91.2 in impregnated dunite and from 85.6 to 88.6 in troctolite, olivine gabbro, and olivine gabbronorite. Spinel Cr# (= 100 x Cr/[Cr + Al] molar ratio) ranges from 38.9 to 62.7 in dunite and from 46.3 to 57.6 in troctolites, olivine gabbro, and olivine gabbronorite. Compositional trends for spinel from dunite through troctolite toward olivine gabbro/gabbronorite are characterized by increases in TiO2, Cr#, and Fe3+#, very similar to those reported from Hess Deep Site 895. Olivine gabbro, olivine gabbronorite, and troctolite in Hole 1271B are considered to have formed as hybrid rocks between dunite and an evolved melt in the walls of a melt channel in the shallowest upper mantle that is tens of meters wide. The melt trapped in the wall rock crystallized plagioclase and clinopyroxene. On the other hand, dunite in the center of the melt channel became more refractory by melt-mantle reactions, increasing spinel Cr# to 62.5.

El problema isostático inverso de Vening Meinesz. Teoría y desarrollo. Aplicación práctica para la determinación de la profundidad de la discontinuidad de Mohorovičić en la Península Ibérica

La discontinuidad de Mohorovičić, más conocida simplemente como “Moho” constituye la superficie de separación entre los materiales rocosos menos densos de la corteza y los materiales rocosos más densos del manto, suponiendo estas capas de densidad constante del orden de 2.67 y 3.27 g/cm3, y es un contorno básico para cualquier estudio geofísico de la corteza terrestre. Los estudios sísmicos y gravimétricos realizados demuestran que la profundidad del Moho es del orden de 30-40 km por debajo de la Península Ibérica y 5-15 km bajo las zonas marinas. Además las distintas técnicas existentes muestran gran correlación en los resultados. Haciendo la suposición de que el campo de gravedad de la Península Ibérica (como le ocurre al 90% de la Tierra) está isostáticamente compensado por la variable profundidad del Moho, suponiendo un contraste de densidad constante entre la corteza y el manto y siguiendo el modelo isostático de Vening Meinesz (1931), se formula el problema isostático inverso para obtener tal profundidad a partir de la anomalía Bouguer de la gravedad calculada gracias a la gravedad observada en la superficie terrestre. La particularidad de este modelo es la compensación isostática regional de la que parte la teoría, que se asemeja a la realidad en mayor medida que otros modelos existentes, como el de Airy-Heiskanen, que ha sido históricamente el más utilizado en trabajos semejantes. Además, su solución está relacionada con el campo de gravedad global para toda la Tierra, por lo que los actuales modelos gravitacionales, la mayoría derivados de observaciones satelitales, deberían ser importantes fuentes de información para nuestra solución. El objetivo de esta tesis es el estudio con detalle de este método, desarrollado por Helmut Moritz en 1990, que desde entonces ha tenido poca evolución y seguidores y que nunca se ha puesto en práctica en la Península Ibérica. Después de tratar su teoría, desarrollo y aspectos computacionales, se está en posición de obtener un modelo digital del Moho para esta zona a fin de poder utilizarse para el estudio de la distribución de masas bajo la superficie terrestre. A partir de los datos del Moho obtenidos por métodos alternativos se hará una comparación. La precisión de ninguno de estos métodos es extremadamente alta (+5 km aproximadamente). No obstante, en aquellas zonas donde exista una discrepancia de datos significaría un área descompensada, con posibles movimientos tectónicos o alto grado de riesgo sísmico, lo que le da a este estudio un valor añadido. ABSTRACT The Mohorovičić discontinuity, simply known as “Moho” constitutes the division between the rocky and less thick materials of the mantle and the heavier ones in the crust, assuming densities of the orders of 2.67 y 3.27 g/cm3 respectively. It is also a basic contour for every geophysical kind of studies about the terrestrial crust. The seismic and previous gravimetric observations done in the study area show that the Moho depth is of the order of 30-40 km beneath the ground and 5-15 km under the ocean basin. Besides, the different techniques show a good correlation in their results. Assuming that the Iberian Peninsula gravity field (as it happens for the 90% of the Earth) is isostatically compensated according to the variable Moho depth, supposing a constant density contrast between crust and mantle, and following the isostatic Vening Meinesz model (1931), the inverse isostatic problem can be formulated from Bouguer gravity anomaly data obtained thanks to the observed gravity at the surface of the Earth. The main difference between this model and other existing ones, such as Airy- Heiskanen’s (pure local compensation and mostly used in these kinds of works) is the approaching to a regional isostatic compensation, much more in accordance with reality. Besides, its solution is related to the global gravity field, and the current gravitational models -mostly satellite derived- should be important data sources in such solution. The aim of this thesis is to study with detail this method, developed by Helmut Moritz in 1990, which hardly ever has it put into practice. Moreover, it has never been used in Iberia. After studying its theory, development and computational aspects, we are able to get a Digital Moho Model of the Iberian Peninsula, in order to study the masses distribution beneath the Earth’s surface. With the depth Moho information obtained from alternative methods, a comparison will be done. Both methods give results with the same order of accuracy, which is not quite high (+ 5 km approximately). Nevertheless, the areas in which a higher difference is observed would mean a disturbance of the compensation, which could show an unbalanced area with possible tectonic movements or potential seismic risk. It will give us an important additive value, which could be used in, at first, non related fields, such as density discrepancies or natural disasters contingency plans.

Major and rare earth element concentrations of basalts beneath ODP Hole 118-735B (Table 3)

A wide-angle seismic experiment at the Atlantis II Fracture Zone, Southwest Indian Ridge, together with geochemical analyses of dredged basalt glass samples from a site conjugate to Ocean Drilling Program hole 735B has allowed determination of the thickness and the most likely lithological composition of the crust beneath hole 735B. The measured Na, composition of 3.3 +/- 0.1 corresponds to a melt thickness of 3 +/- 1 km, a result consistent with rare earth element inversions which indicate a melt thickness of between 1.5 and 4.5 km. The seismic crustal thickness to the north and south of the Atlantis Platform (on which hole 735B is located) is 4 +/- 1 km, and probably consists largely of magmatic material since the seismic and inferred melt thicknesses agree within experimental uncertainty. Beneath hole 735B itself. the Moho is at a depth of 5 +/- 1 km beneath the seafloor. The seismic model suggests that, on average. about 1 km of upper crust has been unroofed on the Atlantis Platform. However, allowing for the inferred local unroofing of 2 km of upper crust at 735B, the base of the magmatic crust beneath this location is probably about 2 km beneath the seafloor, and is underlain by a 2-3 km thick layer of serpentinised mantle peridotite. The P-wave velocity of 6.9 km/s for the serpentinised peridotite layer corresponds to a 35 +/- 10 vol% serpentine content. The Moho beneath hole 735B probably represents a serpentinisation front.

Hybrid multi-objective optimization and hybridized self-organizing response surface method

Numerical optimization is a technique where a computer is used to explore design parameter combinations to find extremes in performance factors. In multi-objective optimization several performance factors can be optimized simultaneously. The solution to multi-objective optimization problems is not a single design, but a family of optimized designs referred to as the Pareto frontier. The Pareto frontier is a trade-off curve in the objective function space composed of solutions where performance in one objective function is traded for performance in others. A Multi-Objective Hybridized Optimizer (MOHO) was created for the purpose of solving multi-objective optimization problems by utilizing a set of constituent optimization algorithms. MOHO tracks the progress of the Pareto frontier approximation development and automatically switches amongst those constituent evolutionary optimization algorithms to speed the formation of an accurate Pareto frontier approximation. Aerodynamic shape optimization is one of the oldest applications of numerical optimization. MOHO was used to perform shape optimization on a 0.5-inch ballistic penetrator traveling at Mach number 2.5. Two objectives were simultaneously optimized: minimize aerodynamic drag and maximize penetrator volume. This problem was solved twice. The first time the problem was solved by using Modified Newton Impact Theory (MNIT) to determine the pressure drag on the penetrator. In the second solution, a Parabolized Navier-Stokes (PNS) solver that includes viscosity was used to evaluate the drag on the penetrator. The studies show the difference in the optimized penetrator shapes when viscosity is absent and present in the optimization. In modern optimization problems, objective function evaluations may require many hours on a computer cluster to perform these types of analysis. One solution is to create a response surface that models the behavior of the objective function. Once enough data about the behavior of the objective function has been collected, a response surface can be used to represent the actual objective function in the optimization process. The Hybrid Self-Organizing Response Surface Method (HYBSORSM) algorithm was developed and used to make response surfaces of objective functions. HYBSORSM was evaluated using a suite of 295 non-linear functions. These functions involve from 2 to 100 variables demonstrating robustness and accuracy of HYBSORSM.