Convergent flow of ice within the Astrolabe subglacial basin, Trre Adélie, East Antarctica: an hypothesis derived from nummerical modelling experiments

The existence of a large subglacial lake beneath the antarctic Ice Sheet at Terre Adélie indicates the presence of basal ice at its pressure-melting temperature. A numerical model of the ice-sheet thermal regime is employed using the balance velocity of the ice sheet as an initial model input in order to calculate ice-sheet basal temperatures. However, the results from this model show the Terre Adélie area to be characterised by basal freezing. Heat in addition to that accounted for in the model is thus required at the ice-sheet base in order for pressure melting temperatures to be attained. The sources for such heat are (1) an enhanced geothermal heat flux and (2) an increase in frictional heating caused by the flow of ice. In this paper the latter possibility is expanded by hypothesising that subglacial topography induces convergent ice flow around Terre Adélie, causing enhanced basal ice velocities. Model experiments indicate that an increase in ice velocity (from 7 to at least 42 m yr⁻¹) is required to raise the temperature of the basal ice to the pressure melting value. Increased ice velocity, and consequent frictional heat production due to convergent ice flow, may therefore be important in explaining the location of the subglacial lake in this region. These results allow the process of convergent ice flow within a contemporary ice sheet to be quantified. A verification (or otherwise) of the model results may be possible if ice surface velocity measurements from modem GPS methods are made.     

The Zihuatanejo, Mexico, earthquake of 1994 December 10 (M= 6.6): source characteristics and tectonic implications

An analysis of the Zihuatanejo, Mexico, earthquake of 1994 December 10 ( M = 6.6), based on teleseismic and near-source data, shows that it was a normal-faulting, intermediate-depth ( H = 50 ± 5 km) event. It was located about 30 km inland, within the subducted Cocos plate. The preferred fault plane has an azimuth of 130°, a dip of 79° and a rake of −86°. The rupture consisted of two subevents which were separated in time by about 2 s, with the second subevent occurring downdip of the first. The measured stress drop was relatively high, requiring a Δσ of about a kilobar to explain the high-frequency level of the near-source spectra. A rough estimate of the thickness of the seismogenic part of the oceanic lithosphere below Zihuatanejo, based on the depth and the rupture extent of this event, is 40 km.
This event and the Oaxaca earthquake of 1931 January 15 ( M = 7.8) are the two significant normal-faulting, intermediate-depth shocks whose epicentres are closest to the coast. Both of these earthquakes were preceded by several large to great shallow, low-angle thrust earthquakes, occurring updip. The observations in other subduction zones show just the opposite: normal-faulting events precede, not succeed, updip, thrust shocks. Indeed, the thrust events, soon after their occurrence, are expected to cause compression in the slab, thus inhibiting the occurrence of normal-faulting events. To explain the occurrence of the Zihuatanejo earthquake, we note that the Cocos plate, after an initial shallow-angle subduction, unbends and becomes subhorizontal. In the region of the unbending, the bottom of the slab is in horizontal extension. We speculate that the large updip seismic slip during shallow, low-angle thrust events increases the buckling of the slab, resulting in an incremental tensional stress at the bottom of the slab and causing normal-faulting earthquakes. This explanation may also hold for the 1931 Oaxaca event.     

Palaeosecular variation over the past 5 Myr based on a new generalized database

A generalized database of global palaeomagnetic data from 3719 lava flows and thin dykes of age 0–5 Ma has been constructed for use with a relational database management system. The database includes all data whose virtual geomagnetic poles (VGP) lie within 45 of the spin axis and can be used for studies of palaeosecular variation and for geomagnetic field modelling. Because many of these data were collected and processed more than 15–20 years ago, each result has been characterized according to the demagnetization procedures carried out. Analysis of these data in terms of the latitude variation of the angular dispersion of VGPs (palaeosecular variation from lavas) strongly suggests that careful data selection is required and that many of the older studies may need to be redone using more modern methods. Differences between the angular dispersions for separate normal- and reverse-polarity data sets confirm that many older studies have not been adequately cleaned magnetically. Therefore, the use of the database for geomagnetic field modelling should be carried out with some caution. Using a VGP cut-off angle that varies with latitude, the best data set consists of 2636 results that show a smooth increase of VGP angular dispersion with latitude. Model G for palaeosecular variation, which is based on modelling of the antisymmetric (dipole) and symmetric (quadrupole) dynamo families, provides a good fit to these results.     

Time-domain electromagnetic migration in the solution of inverse problems

Time-domain electromagnetic (TDEM) migration is based on downward extrapolation of the observed field in reverse time. In fact, the migrated EM field is the solution of the boundary-value problem for the adjoint Maxwell's equations. The important question is how this imaging technique can be related to the solution of the geoelectrical inverse problem. In this paper we introduce a new formulation of the inverse problem, based on the minimization of the residual-field energy flow through the surface or profile of observations. We demonstrate that TDEM migration can be interpreted as the first step in the solution of this specially formulated TDEM inverse problem. However, in many practical situations this first step produces a very efficient approximation to the geoelectrical model, which makes electromagnetic migration so attractive for practical applications. We demonstrate the effectiveness of this approach in inverting synthetic and practical TDEM data.     

Source investigation of a small event using empirical Green's functions and simulated annealing

We propose a two-step inversion of three-component seismograms that (1) recovers the far-field source time function at each station and (2) estimates the distribution of co-seismic slip on the fault plane for small earthquakes (magnitude 3 to 4). The empirical Green's function (EGF) method consists of finding a small earthquake located near the one we wish to study and then performing a deconvolution to remove the path, site, and instrumental effects from the main-event signal.
The deconvolution between the two earthquakes is an unstable procedure: we have therefore developed a simulated annealing technique to recover a stable and positive source time function (STF) in the time domain at each station with an estimation of uncertainties. Given a good azimuthal coverage, we can obtain information on the directivity effect as well as on the rupture process. We propose an inversion method by simulated annealing using the STF to recover the distribution of slip on the fault plane with a constant rupture-velocity model. This method permits estimation of physical quantities on the fault plane, as well as possible identification of the real fault plane.
We apply this two-step procedure for an event of magnitude 3 recorded in the Gulf of Corinth in August 1991. A nearby event of magnitude 2 provides us with empirical Green's functions for each station. We estimate an active fault area of 0.02 to 0.15 km² and deduce a stress-drop value of 1 to 30 bar and an average slip of 0.1 to 1.6 cm. The selected fault of the main event is in good agreement with the existence of a detachment surface inferred from the tectonics of this half-graben.     

Ccte d'Ivoire-Ghana margin: seismic imaging of passive rifted crust adjacent to a transform continental margin

During May 1990 and January-February 1991, an extensive geophysical data set was collected over the Côte d'Ivoire-Ghana continental margin, located along the equatorial coast of West Africa. The Ghana margin is a transform continental margin running subparallel to the Romanche Fracture Zone and its associated marginal ridge—the Côte d'Ivoire-Ghana Ridge. From this data set, an explosive refraction line running ∼ 150 km, ENE-WSW between 3°55'N, 3°21'W and 4°23'N, 2°4'W, has been modelled together with wide-angle airgun profiles, and seismic reflection and gravity data. This study is centred on the Côte d'Ivoire Basin located just to the north of the Côte d'Ivoire-Ghana Ridge, where bathymetric data suggest that a component of normal rifting occurred, rather than the transform motion observed along the majority of the equatorial West African margin.
Traveltime and amplitude modelling of the ocean-bottom seismometer data shows that the continental Moho beneath the margin rises in an oceanward direction, from ∼ 24 km below sea level to ∼ 17 km. In the centre of the line where the crust thins most rapidly, there exists a region of anomalously high velocity at the base of the crust, reaching some 8 km in thickness. This higher-velocity region is thought to represent an area of localized underplating related to rifting. Modelling of marine gravity data, collected coincident with the seismic line, has been used to test the best-fitting seismic model. This modelling has shown that the observed free-air anomaly is dominated by the effects of crustal thickness, and that a region of higher density is required at the base of the crust to fit the observed data. This higher-density region is consistent in size and location with the high velocities required to fit the seismic data.