Table 1: Station data for the entire survey, including coordinates and description of BPBASE

During GANOVEX VI new gravity data were collected along an east-west profile in North Victoria Land south of the Drygalski Ice Tongue, extending 150 km across the Transantarctic Mountains, and comprising 21 data points. Thirty five additional data points were collected over a small area near Brimstone Peak, near the western end of the regional profile. The survey south of the Drygalski has been connected to northern gravity data (GANOVEX V) by a survey line of 12 points. All data have been terrain corrected, and are further constrained by satellite elevation (GPS) and radar ice-thickness measurements. A pronounced regional Bouguer gravity gradient decreasing to the west by approximately 3 mgal/km is superimposed over a coast-parallel belt of granitoid basement rock. West of this belt the local gravity fields become mote variable. Over Beta Peak (Ferrar dolerite) a 50 mgal spike is obser- ved. Within this area, the Ferrar sills are exposed at the surface. West of Brimstone Peak (Ferrar/Kirk patrick sequences), a smooth regional gradient appears to reassert itself. We interpret the initial gradient east (oceanward) of the break-in-slope to be representative of the crust/mantle boundary within the study area. We interpret the initial break-in-slope and the apparent flattening of the regional gradient to be an effect of the N-S trending zone of dense Ferrar sills and associated deep crusttil fractionate replacing less dense basement. We attribute the variability of the local field to be the product of sub-glacial density contrasts that cannot be removed. The regional gravity gradient of the profile is steeper than that observed to the north (Mt. Melbourne quadrangle) and shallower than that reported to the south (McMurdo Sound). The absolute values of the coastal points of origin south of the Drygalski and within the Mt. Melbourne quadrangle differ by 60 to 100 mgal. In addition, topographic relief within the regional transect area is subdued relative to the Transantarctic Mountains to the north and south. We speculate that the root structure of the Transantarctic Mountains undergoes a change somewhere between the Mt. Melbourne quadrangle and the region south of the Drygalski Ice Tongue.

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.

Analysis of flexural isostasy of the northern Andes

The mode in which a lithosphere plate supports overlying topography is greatly driven by the strength of the plate. By analyzing the geophysical signature of lithosphere flexure, in the space and spectral domains, the strength of the plates that support the north Andean mountains and adjacent basins, and the topography of Kenya was investigated. In addition, the effect of windowing on elastic thickness estimates obtained via the coherence method was evaluated. ^ The coherence between the topography and Bouguer gravity spectra of northern South America suggests that the average elastic thickness of the lithosphere is 30 km. Although lateral variations were not resolved by the coherence implementation, these became apparent by modeling the foreland stratigraphy of the Llanos, Barinas and Maracaibo sub-Andean basins. Flexural models reveal a zone of lithosphere weakness beneath the eastern flank of the Eastern Cordillera and western flank of the Venezuelan Andes. The gravity anomaly calculated from these models is consistent with the observed Bouguer gravity anomaly. This zone of weakness appears to separate the strong, old Guyana shield lithosphere from the weaker and probably younger Andean lithosphere. The zone of weakness may correspond to a Paleozoic feature at the western margin of cratonic South America, or a Mesozoic rift arm that weakened the proto-Andean lithosphere. ^ Using synthetic data as well as the northern South America topography and gravity, this study demonstrates that lithosphere strength calculated from the coherence of mirrored data may overestimate the true lithosphere strength. As a result, many lithosphere plates may be weaker than currently thought. In light of this observation, gravity and topography data from Kenya were reevaluated using multitaper spectral techniques. The elastic thickness of this plate, currently undergoing rifting, was estimated at 7 to 8 km, a factor of 2 less than previously estimated. These estimates suggest that despite intense fracturing and sustained tensile stresses, continental lithosphere plates undergoing rifting are able to retain some strength. ^

Using Micro-Gravity Techniques to Map Alluvium Thickness and Pleistocene Location of the West Branch of the Susquehanna River Near Muncy, Pennsylvania

Laurentide glaciation during the early Pleistocene (~970 ka) dammed the southeast-flowing West Branch of the Susquehanna River (WBSR), scouring bedrock and creating 100-km-long glacial Lake Lesley near the Great Bend at Muncy, Pennsylvania (Ramage et al., 1998). Local drill logs and well data indicate that subsequent paleo-outwash floods and modern fluvial processes have deposited as much as 30 meters of alluvium in this area, but little is known about the valley fill architecture and the bedrock-alluvium interface. By gaining a greater understanding of the bedrock-alluvium interface the project will not only supplement existing depth to bedrock information, but also provide information pertinent to the evolution of the Muncy Valley landscape. This project determined if variations in the thickness of the valley fill were detectable using micro-gravity techniques to map the bedrock-alluvium interface. The gravity method was deemed appropriate due to scale of the study area (~30 km2), ease of operation by a single person, and the available geophysical equipment. A LaCoste and Romberg Gravitron unit was used to collect gravitational field readings at 49 locations over 5 transects across the Muncy Creek and Susquehanna River valleys (approximately 30 km2), with at least two gravity base stations per transect. Precise latitude, longitude and ground surface elevation at each location were measured using an OPUS corrected Trimble RTK-GPS unit. Base stations were chosen based on ease of access due to the necessity of repeat measurements. Gravity measurement locations were selected and marked to provide easy access and repeat measurements. The gravimeter was returned to a base station within every two hours and a looping procedure was used to determine drift and maximize confidence in the gravity measurements. A two-minute calibration reading at each station was used to minimize any tares in the data. The Gravitron digitally recorded finite impulse response filtered gravity measurements every 20 seconds at each station. A measurement period of 15 minutes was used for each base station occupation and a minimum of 5 minutes at all other locations. Longer or multiple measurements were utilized at some sites if drift or other externalities (i.e. train or truck traffic) were effecting readings. Average, median, standard deviation and 95% confidence interval were calculated for each station. Tidal, drift, latitude, free-air, Bouguer and terrain corrections were then applied. The results show that the gravitational field decreases as alluvium thickness increases across the axes of the Susquehanna River and Muncy Creek valleys. However, the location of the gravity low does not correspond with the present-day location of the West Branch of the Susquehanna River (WBSR), suggesting that the WBSR may have been constrained along Bald Eagle Mountain by a glacial lobe originating from the Muncy Creek Valley to the northeast. Using a 3-D inversion model, the topography of the bedrock-alluvium interface was determined over the extent of the study area using a density contrast of -0.8 g/cm3. Our results are consistent with the bedrock geometry of the area, and provide a low-cost, non-invasive and efficient method for exploring the subsurface and for supplementing existing well data.

Interpretation of gravity data to delineate structural features connected to low-temperature geothermal resources at Northeastern Portugal

A great number of low-temperature geothermal fields occur in Northern-Portugal related to fractured rocks. The most important superficial manifestations of these hydrothermal systems appear in pull-apart tectonic basins and are strongly conditioned by the orientation of the main fault systems in the region. This work presents the interpretation of gravity gradient maps and 3D inversion model produced from a regional gravity survey. The horizontal gradients reveal a complex fault system. The obtained 3D model of density contrast puts into evidence the main fault zone in the region and the depth distribution of the granitic bodies. Their relationship with the hydrothermal systems supports the conceptual models elaborated from hydrochemical and isotopic water analyses. This work emphasizes the importance of the role of the gravity method and analysis to better understand the connection between hydrothermal systems and the fractured rock pattern and surrounding geology. (c) 2013 Elsevier B.V. All rights reserved.

Gravimetric study of a potential mineral deposit in the Itapororoca region, Brazil

A semi-detailed gravity survey was carried out over an area of 650 km(2) localized in the Eo-Neoproterozoic coastal zone of Paraiba State where 548 new gravity stations were added to the existing database. Gravity measurements were made with a LaCoste and Romberg model G meter with a precision of 0.04 mGal. The altitude was determined by barometric levelling with a fixed base achieving a 1.2 m measure of uncertainty, corresponding to an overall accuracy of 0.24 mGal for the Bouguer anomaly. The residual Bouguer map for a 7th degree regional polynomial showed a circumscribed negative anomaly coincident with a localized aero-magnetic anomaly and with hydro-thermally altered outcrops, near the city of Itapororoca. The 3D gravity modelling, constrained by geologic mapping was interpreted as a low density, fractured and/or altered material with a most probable volume of approximately 23 km(3), extending to about 8,500 m depth. This result is in accordance with a volcanic body associated with hydrothermal processes accompanied by surface mineralization evidence, which may be of interest to the mining industry.

ELKO SPINOR FIELDS: LAGRANGIANS FOR GRAVITY DERIVED FROM SUPERGRAVITY

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Investigation of the Goias Alkaline Province, Central Brazil: Application of gravity and magnetic methods

We investigate the strong magnetic and gravity anomalies of the Goias Alkaline Province (GAP), a region of Late Cretaceous alkaline magmatism along the northern border of the Parana Basin, Brazil. The alkaline complexes (eight of which are present in outcrops, two others inferred from magnetic signals) are characterized by a series of small intrusions forming almost circular magnetic and gravimetric anomalies varying from -4000 to +6000 nT and from -10 to +40 mGal, respectively. We used the Aneuler method and Analytical Signal Amplitude to obtain depth and geometry for mapped sources from the magnetic anomaly data. These results were used as the reference models in the 3D gravity inversion. The 3D inversion results show that the alkaline intrusions have depths of 10-12 km. The intrusions in the northern GAP follow two alignments and have different sizes. In the anomaly magnetic map, dominant guidelines correlate strongly with the extensional regimes that correlate with the rise of alkaline magmatism. The emplacement of these intrusions marks mechanical discontinuities and zones of weakness in the upper crust. According to the 3D inversion results, those intrusions are located within the upper crust (from the surface to 18 km depth) and have spheres as the preferable geometry. Such spherical shapes are more consistent with magmatic chambers instead of plug intrusions. The Registro do Araguaia anomaly (similar to 15 by 25 km) has a particular magnetic signature that indicates that the top is deeper than 1500 m. North of this circular anomaly are lineaments with structural indices indicating contacts on their edges and dikes/sills in the interiors. Results of 3D inversion of magnetic and gravity data suggest that the Registro do Araguaia is the largest body in the area, reaching 18 km depth and indicating a circular layered structure. (C) 2011 Elsevier Ltd. All rights reserved.

Development of one-dimensional and two-dimensional computational tools that simulate steady internal condensing flows in terrestrial and zero-gravity environments

This dissertation presents an effective quasi one-dimensional (1-D) computational simulation tool and a full two-dimensional (2-D) computational simulation methodology for steady annular/stratified internal condensing flows of pure vapor. These simulation tools are used to investigate internal condensing flows in both gravity as well as shear driven environments. Through accurate numerical simulations of the full two dimensional governing equations, results for laminar/laminar condensing flows inside mm-scale ducts are presented. The methodology has been developed using MATLAB/COMSOL platform and is currently capable of simulating film-wise condensation for steady (and unsteady flows). Moreover, a novel 1-D solution technique, capable of simulating condensing flows inside rectangular and circular ducts with different thermal boundary conditions is also presented. The results obtained from the 2-D scientific tool and 1-D engineering tool, are validated and synthesized with experimental results for gravity dominated flows inside vertical tube and inclined channel; and, also, for shear/pressure driven flows inside horizontal channels. Furthermore, these simulation tools are employed to demonstrate key differences of physics between gravity dominated and shear/pressure driven flows. A transition map that distinguishes shear driven, gravity driven, and “mixed” driven flow zones within the non-dimensional parameter space that govern these duct flows is presented along with the film thickness and heat transfer correlations that are valid in these zones. It has also been shown that internal condensing flows in a micro-meter scale duct experiences shear driven flow, even in different gravitational environments. The full 2-D steady computational tool has been employed to investigate the length of annularity. The result for a shear driven flow in a horizontal channel shows that in absence of any noise or pressure fluctuation at the inlet, the onset of non-annularity is partly due to insufficient shear at the liquid-vapor interface. This result is being further corroborated/investigated by R. R. Naik with the help of the unsteady simulation tool. The condensing flow results and flow physics understanding developed through these simulation tools will be instrumental in reliable design of modern micro-scale and spacebased thermal systems.

Tab. 2: Gravity measurements of the Inntal underground at Thaur

A cross-section of the Inn-valley has been surveyed by refraction- and refiection-seismic and gravimetrie methods. The thickness of the Inn-va.!ley sediments is 340- 390 m. At the northern edge of the valley an intermediate layer between sediments and basement has been detected, which is up to 300 ITl thick. This zone seems to mark the boundary of the northern calcareous alps.

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.