High-harmonic geoid signatures related to glacial isostatic adjustment and their detectability by GOCE

The Earth’s asthenosphere and lower continental crust can regionally have viscosities that are one to several orders of magnitude smaller than typical mantle viscosities. As a consequence, such shallow low-viscosity layers could induce high-harmonic (spherical harmonics 50–200) gravity and geoid anomalies due to remaining isostasy deviations following Late-Pleistocene glacial isostatic adjustment (GIA). Such high-harmonic geoid and gravity signatures would depend also on the detailed ice and meltwater loading distribution and history.ESA’s Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, planned for launch in Summer 2008, is designed to map the quasi-static geoid with centimeter accuracy and gravity anomalies with milligal accuracy at a resolution of 100 km or better. This might offer the possibility of detecting gravity and geoid effects of low-viscosity shallow earth layers and differences of the effects of various Pleistocene ice decay scenarios. For example, our predictions show that for a typical low-viscosity crustal zone GOCE should be able to discern differences between ice-load histories down to length scales of about 150 km.One of the major challenges in interpreting such high-harmonic, regional-scale, geoid signatures in GOCE solutions will be to discriminate GIA-signatures from various other solid-earth contributions. It might be of help here that the high-harmonic geoid and gravity signatures form quite characteristic 2D patterns, depending on both ice load and low-viscosity zone model parameters.     

Regional geoid modeling in the area of subglacial Lake Vostok,Antarctica

We present a geoid model for the area of Lake Vostok, Antarctica, from a combination of local airborne gravity, ice-surface and ice-thickness data and a lake bathymetry model. The topography data are used for residual terrain modeling (RTM) in a remove–restore approach together with GOCE satellite data. The quasigeoid is predicted by least-squares collocation (LSC) and subsequently converted to geoid heights. Special aspects of that method in presence of an ice sheet are discussed.It is well known that a body freely floating in water is in a state of hydrostatic equilibrium (HE). This usually applies, e.g., to ice shelves or sea ice. However, it has been shown that this is valid also for the ice sheet covering the subglacial Lake Vostok. Thus, we demonstrate the use of such a refined regional geoid model for glaciological and geophysical applications by means of the HE surface of that lake. The mean quadratic residual geoid signal (0.56 m) w.r.t. the GOCE background model exceeds the residual variations of the estimated apparent lake level (ALL) (0.26 m) within the central part of the lake. An approach considering the actual geopotential at the ALL has been derived and subsequently applied. In this context, downward continuation of the potential field within the ice sheet as well as the latitudinal tilt of off-geoid equipotential surfaces are discussed. In view of the accuracy of the ice-thickness measurements that dominate the total error budget of the estimated ALL these effects are negligible. Thus, the HE surface of subglacial lakes may safely be described by a constant height bias in small-scale regional applications. However, field continuation is significant with respect to the formal uncertainty of the quasigeoid, which is at the level of 5 cm given that accurate airborne gravity data (±2 mGal) are available.     

A comparison of GOCE gravitational models with EGM2008

Four new gravity field models from GOCE, two of them combined with GRACE, are compared here with EGM2008. The objectives are to look into the differences in consecutive ranges of the spherical harmonic expansion globally as well as in selected geographical regions and in the regions of the various data sources used for EGM2008. In general, GOCE is able to contribute to improved global gravity models in the spherical harmonic range between 120 and 200 (and above). The agreement between EGM2008 and the GOCE models is very good in well-surveyed regions such as North America, Europe and Australia, with geoid RMS-differences on the order of 4–6 cm. In other regions, where the surface gravity data available for the development of EGM2008 were poor, such as South America, Africa, South-East Asia or China the RMS-differences are on a level of 30 cm. Here GOCE leads to a significant improvement. These findings are confirmed by the analysis of the areas of the various EGM2008 data sources. In the regions of the so-called “fill-in” data of EGM2008 RMS-geoid height differences are high. In Antarctica GOCE also gives important improvements in terms of spatial resolution and accuracy. In general, the agreement between EGM2008 and the GOCE-models up to degree and order (d/o) 200 is good, with a global (excluding the polar gaps of GOCE orbits, throughout) geoid difference RMS of 11 cm, in the ocean areas 8 cm and 20 cm in the continental areas. GOCE models are better suited for ocean circulation studies because no prior ocean information enters into the data reduction process, as it is the case when deducing gravity anomalies from an altimetric mean sea surface. On the other hand, the good consistency between GOCE-models and EGM2008 in ocean areas very likely indicates that the influence of ocean circulation information on EGM2008 is rather small. The four tested GOCE models behave similarly except at the highest latitudes where GOCE lacks data due to its orbit inclination of 96.5° and some form of regularization which has to be applied.     

The seismically active Andean and Central American margins: Can satellite gravity map lithospheric structures?

The spatial resolution and quality of geopotential models (EGM2008, EIGEN-5C, ITG-GRACE03s, and GOCO-01s) have been assessed as applied to lithospheric structure of the Andean and Central American subduction zones. For the validation, we compared the geopotential models with existing terrestrial gravity data and density models as constrained by seismic and geological data. The quality and resolution of the downward continued geopotential models in the Andes and Central America decrease with increasing topography and depend on the availability of terrestrial gravity data. High resolution of downward continued gravity data has been obtained over the Southern Andes where elevations are lower than 3000 m and sufficient terrestrial gravity data are available. The resolution decreases with an increase in elevation over the north Chilean Andes and Central America. The low resolution in Central America is mainly attributed to limited surface gravity data coverage of the region.To determine the minimum spatial dimension of a causative body that could be resolved using gravity gradient data, a synthetic gravity gradient response of a spherical anomalous mass has been computed at GOCE orbit height (254.9 km). It is shown that the minimum diameter of such a structure with density contrast of 240 kg m⁻³ should be at least ∼45 km to generate signal detectable at orbit height. The batholithic structure in Northern Chile, which is assumed to be associated with plate coupling and asperity generation, is about 60–120 km wide and could be traceable in GOCE data. Short wavelength anomalous structures are more pronounced in the components of the gravity gradient tensor and invariants than in the gravity field.As the ultimate objective of this study is to understand the state of stress along plate interface, the geometry of the density model, as constrained by combined gravity models and seismic data, has been used to develop dynamic model of the Andean margin. The results show that the stress regime in the fore-arc (high and low) tends to follow the trend of the earthquake distributions.     

Constraints on Glacial Isostatic Adjustment from GOCE and Sea Level Data

Current constraints on the glacial isostatic adjustment (GIA) process are mainly provided by relative sea-level data and GPS measurements. Due to a lack of resolving power in the shallow earth (down to about 200 km), these data sets only provide weak constraints on the shallow viscosity structure and the thickness of the lithosphere. Future high-resolution gravity data, as expected from ESA’s Gravity field and steady-state Ocean Circulation Explorer (GOCE) launched on March 17, 2009, are predicted to provide additional information on the shallow earth, more specifically the viscosity structure. Here we present an overview of recent developments in extracting information on rheology and stratification of the shallow earth from high-resolution quasi-steady gravity and geoid data to be obtained from GOCE.     

Improving airborne strapdown vector gravimetry using stabilized horizontal components

Integrating the deflections of the vertical along the flight line can yield geoid profiles which are valuable in the study of geodesy and geophysics, fortunately, the deflections can be measured directly by vector gravimetry. Airborne vector gravimetry using a Strapdown Inertial Navigation System and the Global Navigation Satellite System (SINS/GNSS) has shown promising results in previous studies. However, the quality of the SINS and GNSS is a major limitation; in particular, the attitude errors induced by the gyros will result in large measurement errors to the horizontal components of the gravity disturbance, and these measurement errors represent the behavior of low-frequency trend. An airborne vector gravimetry method used to remove the bias and low-frequency trends in the gravity disturbance estimated for each survey line has been developed. This method uses the horizontal components of the gravity disturbance computed from EGM2008 (Earth Gravitational Model 2008) as a reference. Firstly, the horizontal measurement results obtained from the gravimeter are divided into high- and low-frequency components according to the resolution of the EGM2008, and then, the bias and low-frequency trends of the low-frequency components are corrected using a linear fit to the EGM2008 reference data. Finally, the ultimate results can be acquired after combining the high-frequency components and the corrected low-frequency components. The data used was obtained from the SGA-WZ, which is the first strapdown airborne gravimeter developed in China. The results of this method are promising. The internal accuracy of the gravity disturbance's horizontal components for repeated survey lines exceeds 3.5 mGal, and the corresponding resolution is approximately 4.8 km based on 160-s data smoothing and an airplane averaging speed of approximately 216 km/h. After applying the WCF (Wavenumber Correlation Filter), the internal accuracy of the horizontal components exceeds 2 mGal. This can satisfy the requirement of the application in geodesy and solid earth geophysics.     

Cenozoic building and deformational processes in the North Patagonian Andes

The Oligocene to present evolution of the North Patagonian Andes is analyzed linking geological and geophysical data in order to decipher the deformational processes that acted through time and relate them to basin formation processes. Seismic reflection profiles reveal the shallow structure of the retroarc area where contractional structures, associated with Oligocene to early Miocene inverted extensional depocenters, are partially onlapped by early to late Miocene synorogenic deposits. From the construction of five structural cross sections along the retroarc area between 40° and 43°30′ S, constrained by surface, gravity and seismic data, a shortening gradient is observed along Andean strike. The highest shortening of 18.7 km (15.34%) is determined near 41°30′ S coincidentally with maximum mean topographic values on the eastern Andean slope, where basement blocks were uplifted in the orogenic front area, and the deepest and broadest synorogenic depocenters were formed towards the foreland. Additionally, eastward shifting of Miocene calc-alkaline rocks occurred at these latitudes, which is interpreted as indicative of a change in the subduction parameters at this time. Deep crustal retroarc structure is evaluated through inversion of gravity models that made possible to infer Moho attenuated zones. These coincide with the occurrence of younger than 5 Ma within-plate volcanics as well as with crustal thermal anomalies suggested by shallowing of the Curie isotherm calculated from magnetic data. Younger volcanism and thermal anomalies are explained by slab steepening since early Pliocene, after a mild-shallow subduction setting in the middle to late Miocene, age of the main compressive event.     

Accurate Establishment of Error Models for the Satellite Gravity Gradiometry Recovery and Requirements Analysis for the Future GOCE Follow-On Mission

Firstly, the new single and combined error models applied to estimate the cumulative geoid height error are efficiently produced by the dominating error sources consisting of the gravity gradient of the satellite-equipped gradiometer and the orbital position of the space-borne GPS/GLONASS receiver using the power spectral principle. At degree 250, the cumulative geoid height error is 1.769 × 10^?1 m based on the new combined error model, which preferably accords with a recovery accuracy of 1.760 ×10^?1 m from the GOCE-only Earth gravity field model GO_CONS_GCF_2_TIM_R2 released in Germany. Therefore, the new combined error model of the cumulative geoid height is correct and reliable in this study. Secondly, the requirements analysis for the future GOCE Follow-On satellite system is carried out in respect of the preferred design of the matching measurement accuracy of key payloads comprising the gravity gradient and orbital position and the optimal selection of the orbital altitude of the satellite. We recommend the gravity gradient with an accuracy of 10^?13?10^?15 /s², the orbital position with a precision of 1-0.1 cm and the orbital altitude of 200-250 km in the future GOCE Follow-On mission.     

GOCE, Satellite Gravimetry and Antarctic Mass Transports

In 2009 the European Space Agency satellite mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) was launched. Its objectives are the precise and detailed determination of the Earth’s gravity field and geoid. Its core instrument, a three axis gravitational gradiometer, measures the gravity gradient components V _xx , V _yy , V _zz and V _xz (second-order derivatives of the gravity potential V) with high precision and V _xy , V _yz with low precision, all in the instrument reference frame. The long wavelength gravity field is recovered from the orbit, measured by GPS (Global Positioning System). Characteristic elements of the mission are precise star tracking, a Sun-synchronous and very low (260 km) orbit, angular control by magnetic torquing and an extremely stiff and thermally stable instrument environment. GOCE is complementary to GRACE (Gravity Recovery and Climate Experiment), another satellite gravity mission, launched in 2002. While GRACE is designed to measure temporal gravity variations, albeit with limited spatial resolution, GOCE is aiming at maximum spatial resolution, at the expense of accuracy at large spatial scales. Thus, GOCE will not provide temporal variations but is tailored to the recovery of the fine scales of the stationary field. GRACE is very successful in delivering time series of large-scale mass changes of the Antarctic ice sheet, among other things. Currently, emphasis of respective GRACE analyses is on regional refinement and on changes of temporal trends. One of the challenges is the separation of ice mass changes from glacial isostatic adjustment. Already from a few months of GOCE data, detailed gravity gradients can be recovered. They are presented here for the area of Antarctica. As one application, GOCE gravity gradients are an important addition to the sparse gravity data of Antarctica. They will help studies of the crustal and lithospheric field. A second area of application is ocean circulation. The geoid surface from the gravity field model GOCO01S allows us now to generate rather detailed maps of the mean dynamic ocean topography and of geostrophic flow velocities in the region of the Antarctic Circumpolar Current.     

Geoid anomalies over Gorringe Ridge,North Atlantic Ocean

Gorringe Ridge is a strong uplifted block of oceanic crust and upper mantle lying at the eastern end of the Azores-Gibraltar plate boundary. The geoid over this structure derived from Seasat altimeter data exhibits a 9-m height anomaly with a north-south lateral extension smaller than 200 km. An attempt is made to interpret this geoid together with the gravity anomalies and with the seismicity, which has been compiled as a function of depth.It is first shown that the flexure of the oceanic lithosphere due to the ridge loading does not provide a good fit of the geoid anomalies and probably should be discarded, as it assumes a continuous unfractured elastic plate.Models involving local heterogeneities are then tested. The comparison of the observed geoid anomalies with the anomalies due to the uncompensated relief indicates that the topographic high has no shallow compensation.Uncompensated models, previously proposed to explain the gravity anomalies, are tested using the geoid. One model (Purdy and Bonnin, in Bonnin [11]), which involves an uplift of upper mantle material at depth, generates too strong geoid anomalies and must be discarded. Another model, which represents a nascent subduction zone (Le Pichon et al. [25]), fits both the gravity and geoid anomalies, but leads to difficulties in explaining the deep seismicity north of Gorringe Ridge.A model in isostatic equilibrium is also able to fit both gravity and geoid anomalies. This model involves a deep root of density 3.0 g cm^?3, as has been previously proposed for many oceanic ridges and plateaus. This model is compatible with the deep seismicity, but the origin of this low-density material at great depth is up to now an unresolved question.More likely, dynamical models taking into account the forces induced by the convection related to the slow plate convergence in this area will have to be considered.     

A data sequence acquired at Mt. Etna during the 2002–2003 eruption highlights the potential of continuous gravity observations as a tool to monitor and study active volcanoes

A 2.5-month long gravity sequence, encompassing the starting period of the 2002–2003 Etna eruption and coming from a summit station only 1 km away from the new fractures, is presented and discussed. The sequence comprises four hours-long anomalies that have a great chance to reflect mass redistributions linked to the ensuing activity. In particular, the start of the eruptive activity on the northeastern flank was marked by a gravity decrease as strong as about 400 μGal, which reverted soon afterwards. This strong decrease/increase anomaly is interpreted as the opening, by tectonic forces, of a fracture system along the Northeastern Rift of Mt. Etna, followed by an intrusion of magma from the central conduit to the new fractures. They were used by the intruding magma as a path to the eruptive vents at lower elevations.Afterwards, on three occasions, in November and December 2002, 6–12 h-lasting gravity decreases, with amplitude ranging between 10 and 30 μGal, were observed simultaneously with increases in the amplitude of the volcanic tremor from four seismic stations. A correlation analysis, between the gravity signal and the overall spectral amplitude of each tremor sequence is performed over the 7 November–9 December period. A marked anti-correlation is found over each contemporaneous gravity decrease/tremor increase, while, over the rest of the investigated period, the correlation is negligible. Accordingly, a joint source is inferred to have acted during the occurrence of the three common anomalies. On the grounds of some volcanological observations spanning the period covered by our analysis, we propose the temporary accumulation of a gas cloud at some level within the plumbing system of the volcano to have acted as a joint source.The present work is a further evidence of the potential of continuous gravity observations as a tool to monitor and study active volcanoes and encourages their employment in spite of the difficulty of running spring gravimeters in a continuous fashion under the adverse conditions normally encountered on the summit zone of an active volcano.     

Electrical imaging of Narmada–Son Lineament Zone,Central India from magnetotellurics

The Narmada–Son Lineament (NSL) Zone is the second most important tectonic feature after Himalayas, in the Indian geology. Magnetotelluric (MT) studies were carried out in the NSL zone along a 130 km long NNE-SSW trending profile. The area of investigation extends from Edlabad (20°46′16″; 75°59′05″) in the South to Khandwa (21°53′51″; 76°18′05″) in the North. The data shows in general the validity of a two-dimensional (2D) approach. Besides providing details on the shallow crustal section, the 2D modeling results resolved four high conductive zones extending from the middle to deep crust, spatially coinciding with the major structural features in the area namely the Gavligarh, Tapti, Barwani-Sukta and Narmada South faults. The model for the shallow section has brought out a moderately resistive layer (30–150 Ω m) representing the exposed Deccan trap layer, overlying a conductive layer (10–30 Ω m) inferred to be the subtrappean Gondwana sediments, the latter resting on a high resistive basement/upper crust. The Deccan trap thickness varies from around a few hundred meters to as much as 1.5 km along the traverse. A subtrappean sedimentary basin like feature is delineated in the northern half of the traverse where a sudden thickening of subtrappean sediments amounting to as much as 2 km is noticed. The high resistive upper crust is relatively thick towards the southern end and tends to become thinner towards the middle and northern part of the traverse. The lower crustal segment is conductive over a major part of the profile. Considering the generally enhanced heat flow values in the NSL region, coupled with characteristic gravity highs and enhanced seismic velocities coinciding with the mid to lower crustal conductors delineated from MT, presence of zones of high density mafic bodies/intrusives with fluids, presumably associated with magmatic underplating of the crust in the zone of major tectonic faults in NSL region are inferred.     

Imaging Pn velocities beneath the Pannonian basin

The Pannonian depression is an extensional back-arc basin in central Europe and is an integral part of the Alpine–Carpathian orogenic mountain belts. It can be characterized by thinned lower crust, shallow Moho discontinuity, high surface heat flow and Moho temperature, implying recent active tectonic processes. Imaging the velocity structure of the upper mantle may help us to better understand the structure and formation of the Pannonian region.In this paper, P_n traveltimes from regional earthquakes are used to tomographically image the lateral velocity variations in the uppermost mantle beneath the Pannonian basin. The set of linear tomographic equations, built up of the time term equation for each source–receiver pair, is solved by a truncated singular value decomposition algorithm. The explicit computation of the generalized inverse of the tomographic equations makes it possible to deduce both the resolution matrix and the model covariance matrix, allowing us to estimate the resolution and reliability of the solution.The mean compressional wave velocity in the uppermost mantle beneath the Pannonian basin is 7.9 km/s, substantially lower than the average continental P_n velocity of 8.1 km/s. It is mostly due to the high Moho temperature having values on average 400–500 °C more than those in the surrounding areas. The velocity anomalies range from −0.3 to 0.3 km/s relative to the mean velocity of 7.9 km/s. Due to high Moho temperature, below the North Hungarian range low (7.6–7.7 km/s) velocities can be found. High-velocity anomalies of around 8.1 km/s can be detected along the W-SW boundaries of Hungary and at the junction of the Pannonian basin and the Southern Carpathians. The Great Hungarian Plain shows average (7.9 km/s) P_n velocities.     

Elastic thickness structure of South America estimated using wavelets and satellite-derived gravity data

We used a wavelet formulation of the classical spectral isostatic analysis to invert satellite-derived gravity and topography/bathymetry for elastic thickness (T_e) over South America and its surrounding plates. To provide a homogeneous representation of the gravity field for this vast region, we corrected free-air anomalies derived from a combination of terrestrial/marine gravity data with data from the GRACE and CHAMP satellite missions (model EIGEN-CG03C) by a simple Bouguer slab using a smoothed representation of surface relief (wavelengths > 125 km). The resulting Bouguer anomaly compares well with terrestrial data acquired in the Central Andes and allows T_e to be confidently estimated for values greater than 10 km. The T_e map resolves regional-scale features that are well-correlated with known surface structures and shows maximum values of 100 ± 15 km over the Archean–Neoproterozoic core of the continent, decreasing to less than 30 km around continental margins. Several regions of the oceanic plates and continental margins have an elastic thickness less than 10 km. We performed a quantitative analysis by comparing the elastic thickness with the thermal structure predicted from the age of oceanic crust and igneous–metamorphic rocks. This demonstrates that oceanic plates have been weakened by thermal interaction with hotspots and locally by fracturing and hydration near the trench. We observe that only the nucleus of the continent has resisted the thermomechanical weakening induced by the rifting of Africa and South America along the passive margin and the Andean orogeny along the active margin. This latter region shows along-strike variations in T_e that correlate with the geotectonic segmentation of the margin and with the pattern of crustal seismicity. Our results reveal that the rigidity structure follows the segmentation of the seismogenic zone along the subduction fault, suggesting a causal relationship that should be investigated in order to improve the understanding and predictability of great earthquakes and tsunamis.     

Assimilation of geodetic dynamic ocean topography using ensemble based Kalman filter

Estimation of ocean circulation is investigated via assimilation of satellite measurements of the dynamic ocean topography (DOT) into the global finite-element ocean model (FEOM). The DOT was obtained by means of a geodetic approach from carefully cross-calibrated multi-mission altimeter data and GRACE gravity fields. The spectral consistency was achieved by consistently filtering both, the sea surface and the geoid. The filter length is determined by the spatial resolution of the gravity field and corresponds to approximately 241 km half width for the GRACE-based gravity field model ITG-Grace03s.The assimilation of the geodetic DOT was performed by employing a local singular evolutive interpolated Kalman (SEIK) filter in combination with the method of weighting of observations. It is shown that this approach leads to a successful assimilation technique that reduced the RMS difference between the model and the data from 16 cm to 5 cm during one year of assimilation. The ocean model returns an optimized mean dynamic ocean topography. The effects of assimilation on transport estimates across several hydrographic World Ocean Circulation Experiment (WOCE) sections show improvements compared to the FEOM run without data assimilation. As a result of the assimilation, DOT estimates are available in the polar or coastal regions where the geodetic estimates from satellite data alone are not adequate. Furthermore, more realistic features of the ocean can be seen in these areas compared to those obtained using the filtered data fields.     

Density modeling of the Escollos Alijos Seamount from inversion of its geoid undulation anomaly

Escollos Alijos is a large seamount located in the NE Pacific Ocean about 300 km off the Baja California Peninsula. Geochronology and geochemical analysis of volcanic rocks capping the seamount indicate recent magmatism that resulted from extensive differentiation of a mildly alkalic basalt parent magma.Escollos Alijos is located towards the eastern edge of a long-wavelength geoid undulation minimum, of up to -47 m with respect of the WGS84 ellipsoid, which extends over the northeastern Pacific Ocean. Subtracting from the geoid undulation its long-wavelength component and the undulation due to the seamount topography itself, a negative undulation anomaly persists that indicates a mass deficit at depth. Linear inversion of the undulation anomaly yields a region characterized by a negative density contrast, localized under the seamount at a depth between 9 and 13 km.The age and chemical composition of Escollos Alijos, and the inferred mass deficit suggest magma trapped between the oceanic crust and the uppermost mantle, which explains the magmatic activity in recent times.     

Using ground gravity to improve ice mass change estimation from GOCE gravity gradients in mid-west Greenland

Vertical gravity gradient anomalies from the Gravity and steady-state Ocean Circulation Explorer (GOCE) DIR-3 model have been used to determine gravity anomalies in mid-west Greenland by using Least-Squares Collocation (LSC) and the Reduced Point Mass (RPM) method. The two methods give nearly identical results. However, compared to LSC, the RPM method needs less computational time as the number of equations to be solved in LSC equals the number of observations. The advantage of the LSC, however, is the acquired error estimates. The observation periods are winter 2009 and summer 2012. In order to enhance the accuracy of the calculated gravity anomalies, ground gravity data from West Greenland is used over locations where the gravity change resulting from ice mass changes is negligible, i.e. over solid rock. In the period considered, the gravity anomaly change due to changes in ice mass varies from ?5 mGal to 4 mGal. It is negative over the outlet glacier Jacobshavn Isbræ, where the mass loss corresponds to a gravity change of approximately ?4 mGal. When using only GOCE vertical gravity gradients, the error estimates range from 5 mGal at the coast to 17 mGal over the ice sheet. Introducing the ground gravity data from West Greenland in the prediction reduces the errors to range from 2 to 10 mGal.     

Volcanic systems and calderas in the Vatnajökull region,central Iceland: Constraints on crustal structure from gravity data

A Bouguer anomaly map is presented of southern central Iceland, including the western part of Vatnajökull and adjacent areas. A complete Bouguer reduction for both ice surface and bedrock topography is carried out for the glaciated regions. Parts of the volcanic systems of Vonarskarð-Hágöngur, Bárðarbunga-Veiðivötn, Grímsvötn-Laki, and to a lesser extent Kverkfjöll, show up as distinct features on the gravity map. The large central volcanoes with calderas: Vonarskarð, Bárðarbunga, Kverkfjöll and Grímsvötn, are associated with 15–20 mGal gravity highs caused by high density bodies in the uppermost 5 km of the crust. Each of these bodies is thought to be composed of several hundred km³ of gabbros that have probably accumulated over the lifetime of the volcano. The Skaftárkatlar subglacial geothermal areas are not associated with major anomalous bodies in the upper crust. The central volcanoes of Vonarskarð and Hágöngur belong to the same volcanic system; this also applies to Bárðarbunga and Hamarinn, and Grímsvötn and Þórðarhyrna. None of the smaller of the two volcanoes sharing a system (Hágöngur, Hamarinn and Þórðarhyrna) is associated with distinct gravity anomalies and clear caldera structures have not been identified. However, ridges in the gravity field extend between each pair of central volcanoes, indicating that they are connected by dense dyke swarms. This suggests that when two central volcanoes share the same system, one becomes the main pathway for magma, forming a long-lived crustal magma chamber, a caldera and large volume basic intrusive bodies in the upper crust. Short residence times of magma in the crust beneath these centres favour essentially basaltic volcanism. In the case of the second, auxillary central volcano, magma supply is limited and occurs only sporadically. This setting may lead to longer residence times of magma in the smaller central volcanoes, favouring evolution of the magma and occasional eruption of rhyolites. The eastern margin of the Eastern Volcanic Zone is marked by a NE–SW lineation in the gravity field, probably caused by accumulation of low density, subglacially erupted volcanics within the volcanic zone. This lineation lies 5–10 km to the east of Grímsvötn.     

Structure and development of the Andean system between 36° and 39°S

One of the main morphological changes along the Southern Central Andes occurs from 36° to 39°S. The northern portion is characterized by prominent basement structures and a thick-skinned orogenic front with relief of over 2000 m with a deep level of exhumation where more than 4 km of section has been eroded. Contrastingly, the southern part is formed by mildly inverted basement structures restricted mainly to the hinterland zone, which reaches only 1500–1700 m relief. We quantify the variable contributions of two main contractional stages through the construction of three regionally balanced sections across the Andes, constrained by field and geophysical data. Extensional re-activation described for this segment in late Oligocene-early Miocene and Pliocene to Quaternary times, after the two main contractional episodes, suggests only 3 km of stretching that represents 30–10% of the original longitude. We, therefore, conclude that while initial Late Cretaceous to Eocene compression was similar along strike (∼10–7 km), it is the contrasting degrees of Neogene shortening (∼16–6 km) that have played the largest role in the along strike differences in structure and morphology along this portion of the southern Andes. Variable Neogene arc expansion could be responsible for the contrasting contractional deformation: In the north, late Miocene arc-related rocks cover most of the retroarc zone (>200 km with respect to the late Miocene arc front in the south), presumably driven by a shallow subduction episode in the area, whereas to the south they remain restricted to the continental drainage divide. Other factors involving architecture of previous rift structures, are proposed as additional mechanisms that accommodated variable shortening magnitudes through inversion.     

Gravity wave observations in the summertime polar mesosphere from the Cloud Imaging and Particle Size (CIPS) experiment on the AIM spacecraft

We present the first results of gravity wave signatures on polar mesospheric clouds (PMCs) during the summer of 2007, in the northern hemisphere polar region. The Cloud Imaging and Particle Size (CIPS) experiment has one of the three instruments on board the NASA Aeronomy of Ice in the Mesosphere (AIM) spacecraft, which was launched into a sun-synchronous orbit on April 25, 2007. CIPS is a four-camera, wide-field (120°×80°) imager designed to measure PMC morphology and particle properties. One of the objectives of AIM is to investigate gravity wave effects on PMC formation and evolution. CIPS images show distinct wave patterns and structures in PMCs that are similar to ground-based photographs of noctilucent clouds (NLCs). The observed horizontal wavelengths of the waves were found to vary between 15 and 320 km, with smaller-wavelength structures of less than 50 km being the most common. In this paper we present examples of individual quasi-monochromatic wave events observed by CIPS and statistics on the wave patterns observed in the northern hemisphere during the summer months of 2007, together with a map showing the geographic locations of gravity wave events observed from CIPS.