The Canary Islands swell: a coherence analysis of bathymetry and gravity

The Canary Archipelago is an intraplate volcanic chain, located near the West African continental margin, emplaced on old oceanic lithosphere of Jurassic age, with an extended volcanic activity since Middle Miocene. The adjacent seafloor does not show the broad oceanic swell usually observed in hotspot-generated oceanic islands. However, the observation of a noticeable depth anomaly in the basement west of the Canaries might indicate that the swell is masked by a thick sedimentary cover and the influence of the Canarian volcanism. We use a spectral approach, based on coherence analysis, to determine the swell and its compensation mechanism. The coherence between gravity and topography indicates that the swell is caused by a subsurface load correlated with the surface volcanic load. The residual gravity/geoid anomaly indicates that the subsurface load extends 600 km SSW and 800 km N and NNE of the islands. We used computed depth anomalies from available deep seismic profiles to constrain the extent and amplitude of the basement uplift caused by a relatively low-density anomaly within the lithospheric mantle, and coherence analysis to constrain the elastic thickness of the lithosphere ( T_e ) and the compensation depth of the swell. Depth anomalies and coherence are well simulated with T_e =28–36 km, compensation depth of 40–65 km, and a negative density contrast within the lithosphere of ∼33 kg m⁻³. The density contrast corresponds to a temperature increment of ∼325°C, which we interpret to be partially maintained by a low-viscosity convective layer in the lowermost lithosphere, and which probably involves the shallower parts of the asthenosphere. This interpretation does not require a significant rejuvenation of the mechanical properties of the lithosphere.     

Flexural uplift of the Stara Planina range, central Bulgaria

The Stara Planina is an E–W-trending range within the Balkan belt in central Bulgaria. This topographically high mountain range was the site of Mesozoic through early Cenozoic thrusting and convergence, and its high topography is generally thought to have resulted from crustal shortening associated with those events. However, uplift of this belt appears to be much younger than the age of thrusting and correlates instead with the age of Pliocene–Quaternary normal faulting along the southern side of the range. Flexural modelling indicates the morphology of the range is consistent with flexural uplift of footwall rocks during Pliocene–Quaternary displacement on S-dipping normal faults bounding the south side of the mountains, provided that the effective elastic plate thickness of 12  km under the Moesian platform is reduced to about 3  km under the Stara Planina. This small value of elastic plate thickness under the Stara Planina is similar to values observed in the Basin and Range Province of the western United States, and suggests that weakening of the lithosphere is due to heating of the lithosphere during extension, perhaps to the point that large-scale flow of material is possible within the lower crust. Because weakening is observed to affect the Moesian lithosphere for ≈10  km beyond (north of) the surface expression of extension, this study suggests that processes within the uppermost mantle, such as convection, play an active role in the extension process. The results of this study also suggest that much of the topographic relief in thrust belts where convergence is accompanied by coeval extension in the upper plate (or 'back arc'), such as in the Apennines, may be a flexural response to unloading during normal faulting, rather than a direct response to crustal shortening in the thrust belt.     

Internal structure and eruptive history of a kilometre-scale mud volcano system, South Caspian Sea

We describe the internal structure of a multi‐kilometre scale mud volcano edifice from the South Caspian Sea using three‐dimensional (3D) seismic reflection data leading to a reconstruction of the volcano system's eruptive history. By adapting elements of classic seismic stratigraphy to the study of this volcano, we have found its edifice to consist of a series of stacked mud cones. This internal architecture is most likely to have formed as a result of repeated episodes of expulsion of a fluid‐mud mix. Underlying the stack of cones is an asymmetric fault‐bounded caldera measuring approximately 2 km in diameter. This caldera shows close structural similarity to the trapdoor type of magmatic caldera. Based on the geometrical relationships of individual mud cones to this caldera, we conclude that caldera‐like collapse of the edifice floor initiated following the deposition of the first mud cone (the pioneer cone). Growth of the caldera continued until the later stages of edifice evolution when it eventually abated. This eruptive history shows strong similarities to recent models for magmatic caldera eruption cycles. The study therefore highlights the potential analogue value of mud volcano systems to the study of igneous volcanism. Furthermore, it identifies 3D seismic data as a potentially useful tool in reconstructing the history of mud volcanic eruption and fluid and sediment expulsion from sedimentary basins.     

Localization of gravity and topography: constraints on the tectonics and mantle dynamics of Venus

We develop a method for spatio-spectral localization of harmonic data on a sphere and use it to interpret recent high-resolution global estimates of the gravity and topography of Venus in the context of geodynamical models. Our approach applies equally to the simple spatial windowing of harmonic data and to variable-length-scale analyses, which are analogous to a wavelet transform in the Cartesian domain. Using the variable-length-scale approach, we calculate the localized RMS amplitudes of gravity and topography, as well as the spectral admittance between the two fields, as functions of position and wavelength. The observed admittances over 10 per cent of the surface of Venus (highland plateaus and tessera regions) are consistent with isostatic compensation of topography by variations in crustal thickness, while admittances over the remaining 90 per cent of the surface (rises, plains and lowlands) indicate that long-wavelength topography is dominantly the result of vertical convective tractions at the base of the lithosphere. The global average crustal thickness is less than 30 km, but can reach values as large as 40 km beneath tesserae and highland plateaus. We also note that an Earth-like radial viscosity structure cannot be rejected by the gravity and topography data and that, without a mechanical model of the lithosphere, admittance values cannot constrain the thickness of the thermal boundary layer of Venus. Modelling the lithosphere as a thin elastic plate indicates that at the time of formation of relief in highland plateaus and tesserae, the effective elastic plate thickness, T_e , was less than 20 km. Estimates of T_e at highland rises are consistently less than 30 km. Our inability to find regions with T_e > 30 km is inconsistent with predictions made by a class of catastrophic resurfacing models.     

Precise location of unsurveyed seamounts in the Austral archipelago area using SEASAT data

Summary. SEASAT altimetric geoid data are used to detect uncharted sea-mounts in the Austral archipelago area. The various physical parameters which affect the geoid signature of a seamount are inspected to analyse their influence on the precision of the location. These parameters are the shape of the seamount, its density, the crustal model, the effective elastic thickness of the lithosphere and the distance between the seamount and sub-satellite tracks, respectively. The results of the tests performed on synthetic data as well as on charted structures show that when a seamount is detected on at least two neighbouring tracks, it is possible to locate it with good confidence and to give an estimation of its height. If the correct elastic thickness is assumed, the precision on the located is order 15 km. Ten previously unsurveyed seamounts have been located in the Austral archipelago. It appears that they are emplaced along two well-defined azimuths (N 110°E and N 150°E) and that two parallel distinct volcanic chains form the Austral archipelago.     

Small-scale mantle flows and induced lithospheric stress near island arcs

Summary. This paper explores the middle ground between complex thermally-coupled viscous flow models and simple corner flow models of island arc environments. The calculation retains the density-driven nature of convection and relaxes the geometrical constraints of corner flow, yet still provides semianalytical solutions for velocity and stress. A novel aspect of the procedure is its allowance for a coupled elastic lithosphere on top of a Newtonian viscous mantle. Initially, simple box-like density drivers illustrate how vertical and horizontal forces are transmitted through the mantle and how the lithosphere responds by trench formation. The flexural strength of the lithosphere spatially broadens the surface topography and gravity anomalies relative to the functional form of the vertical flow stresses applied to the plate base. I find that drivers in the form of inclined subducting slabs cannot induce self-driven parallel flow; however, the necessary flow can be provided by supplying a basal drag of 1–5 MPa to the mantle from the oceanic lithosphere. These basal drag forces create regional lithospheric stress and they should be quantifiable through seismic observations of the neutral surface. The existence of a shallow elevated phase transition is suggested in two slab models of 300 km length where a maximum excess density of 0.2 g cm⁻³ was needed to generate an acceptable mantle flow. A North New Hebrides subduction model which satisfies flow requirements and reproduces general features of topography and gravity contains a high shear stress zone (75 MPa) around the upper slab surface to a depth of 150 km and a deviatoric tensional stress in the back arc to a depth of 70 km. The lithospheric stress state of this model suggests that slab detachment is possible through whole plate fracture.     

Geomorphological evidence of the influence of pre-volcanic basement structure on emplacement and deformation of volcanic edifices at the Cofre de Perote–Pico de Orizaba chain and implications for avalanche generation

Pre-volcanic structure of the basement influences volcanism distribution and avalanche generation in volcanic edifices. Therefore, systematic studies of basement structure below volcanic chains are necessary to understand the deformation effects observed in the surface and vice versa. Based on a compilation of pre-existing data, interpretation of aerial photographs and satellite images, and a collection of structural data we analyzed morphological and structural features of the Cofre de Perote–Pico de Orizaba (CP–PO) volcanic chain and its basement. We have identified three sets of regional lineaments that are related to basement trends. (1) NW 55° SE fractures are parallel to anticline folds observed in Cretaceous rocks that originated during Laramide shortening. These folds present an abrupt morphology observed only in the eastern flank but that is likely to continue below the volcanic chain. (2) NE 55° SW fractures are parallel to normal faults at the basement. We infer that these basement faults confine the CP–PO chain within a stepped graben with a total normal displacement of about 400 m. These faults have been active through time since they have affected volcanic deposits and induced the emplacement of monogenetic vents. Notably, lineaments of monogenetic vents concentrate where the basement is relatively shallow. (3) Another set of faults, oriented N–S, has been observed affecting the scarce basement outcrops at the western flank of the chain covered by lacustrine deposits. Lineaments measured in the volcanic edifice of Pico de Orizaba correlate with the regional trends.In particular, the NE 55° SW alignment of monogenetic vents and fractures at Pico de Orizaba suggest that the same dike trend exists within the volcanic edifice. A normal fault with similar orientation was documented at the NE continuation of an alignment crossing the volcanic edifice along the Jamapa canyon. In the absence of magmatic activity related to collapses, the displacement of NE 55° SW faults represents a potential triggering mechanism for generating avalanches at Pico de Orizaba volcano. Instability is enhanced by the presence of N–S trending fractures crossing the entire volcanic edifice and E–W fractures affecting only the present day cone. We conclude that mechanical instability of the volcanic chain is influenced by the basement structure heterogeneity, but further detailed studies are necessary at individual volcanoes to evaluate their effects on volcano deformation.     

Effect of the viscoelastic lithosphere on polar wander speed caused by the Late Pleistocene glacial cycles

Previous studies of the wander of the rotation pole associated with the Late Pleistocene glacial cycles indicate that the predicted polar wander speed is sensitive to the density jump at the 670 km discontinuity, the thickness of the elastic lithosphere, and the lower mantle viscosity. In particular, the M1 mode related to the density jump at 670 km depth has been shown to contribute a dominant portion of predicted polar wander speed for sufficiently small lower mantle viscosities. In this study, we examine the sensitivity of polar wander to variations in the viscosity of the viscoelastic lithosphere using simplified compressible Maxwell viscoelastic earth models. Model calculations for earth models with a viscoelastic lithosphere of finite viscosity indicate that the contribution of the M1 mode is similar to those associated with the density discontinuity at the core–mantle boundary (C0 mode) and the lithosphere (L0 mode). We speculate that this is due to the interaction between the M1 mode and the transient mode associated with the viscoelastic lithosphere, which reduces the magnitude of polar wander rates. Therefore, the M1 mode does not contribute a dominant portion of the predicted polar wander speed for earth models with a viscoelastic lithosphere of finite viscosity. In this case, predictions of polar wander speed as a function of lower mantle viscosity exhibit the qualitative form of an 'inverted parabola', as predicted for the J ˙₂ curve. We caution, however, that these results are obtained for simplified earth models, and the results for seismological earth models such as PREM may be complicated by the interaction between the M1 mode and the large set of transient modes.     

Modelling of ground deformation and gravity fields using finite element method: an application to Etna volcano

Elastic finite element models are applied to investigate the effects of topography and medium heterogeneities on the surface deformation and the gravity field produced by volcanic pressure sources. Changes in the gravity field cannot be interpreted only in terms of gain of mass disregarding the ground deformation of the rocks surrounding the source. Contributions to gravity changes depend also on surface and subsurface mass redistribution driven by dilation of the volcanic source. Both ground deformation and gravity changes were firstly evaluated by solving a coupled axisymmetric problem to estimate the effects of topography and medium heterogeneities. Numerical results show significant discrepancies in the ground deformation and gravity field compared to those predicted by analytical solutions, which disregard topography, elastic heterogeneities and density subsurface structures. With this in mind, we reviewed the expected gravity changes accompanying the 1993–1997 inflation phase on Mt Etna by setting up a fully 3-D finite element model in which we used the real topography, to include the geometry, and seismic tomography, to infer the crustal heterogeneities. The inflation phase was clearly detected by different geodetic techniques (EDM, GPS, SAR and levelling data) that showed a uniform expansion of the overall volcano edifice. When the gravity data are integrated with ground deformation data and a coupled FEM modelling was solved, a mass intrusion could have occurred at depth to justify both ground deformation and gravity observations.     

Massif Central (France): new constraints on the geodynamical evolution from teleseismic tomography

The Massif Central, the most significant geomorphological unit of the Hercynian belt in France, is characterized by graben structures which are part of the European Cenozoic Rift System (ECRIS) and also by distinct volcanic episodes, the most recent dated at 20 Ma to 4000 years BP. In order to study the lithosphere-asthenosphere system beneath this volcanic area, we performed a teleseismic field experiment.
During a six-month period, a joint French-German team operated a network of 79 mobile short-period seismic stations in addition to the 14 permanent stations. Inversion of P -wave traveltime residuals of teleseismic events recorded by this dense array yielded a detailed image of the 3-D velocity structure beneath the Massif Central down to 180 km depth. The upper 60 km of the lithosphere displays strong lateral heterogeneities and shows a remarkable correlation between the volcanic provinces and the negative velocity perturbations. The 3-D model reveals two channels of low velocities, interpreted as the remaining thermal signature of magma ascent following large lithospheric fractures inherited from Hercynian time and reactivated during Oligocene times. The teleseismic inversion model yields no indication of a low-velocity zone in the mantle associated with the graben structures proper. The observation of smaller velocity perturbations and a change in the shape of the velocity pattern in the 60–100 km depth range indicates a smooth transition from the lithosphere to the asthenosphere, thus giving an idea of the lithosphere thickness. A broad volume of low velocities having a diameter of about 200 km from 100 km depth to the bottom of the model is present beneath the Massif Central. This body is likely to be the source responsible for the volcanism. It could be interpreted as the top of a plume-type structure which is now in its cooling phase.     

Deglacial land emergence and lateral upper-mantle heterogeneity in the Svalbard Archipelago—II. Extended results for high-resolution load models

In Paper I (Breuer & Wolf 1995), a preliminary interpretation of the postglacial land emergence observed at a restricted set of six locations in the Svalbard Archipelago was given. The study was based on a simple model of the Barents Sea ice sheet and suggested increases in lithosphere thickness and asthenosphere viscosity with increasing distance from the continental margin.
In the present paper, the newly developed high-resolution load model. BARENTS-2, and land-uplift observations from an extended set of 25 locations are used to study further the possibility of resolving lateral heterogeneity in the upper mantle below the northern Barents Sea. A comparison of the calculated and observed uplift values shows that the lithosphere thickness is not well resolved by the observations, although values above 110 km are most common for this parameter. In contrast to this, there are indications of a lateral variation of asthenosphere viscosity. Whereas values in the range 10¹⁸-10²⁰Pas are inferred for locations close to the continental margin, 10²⁰-10²¹ Pa s are suggested further away from the margin.
A study of the sensitivity of the values found for lithosphere thickness and asthenosphere viscosity to modifications of load model BARENTS-2 shows that such modifications can be largely accommodated by appropriate changes in lithosphere thickness, whereas the suggested lateral variation of asthenosphere viscosity is essentially unaffected. An estimate of the influence of the Fennoscandian. ice sheet leads to the conclusion that its neglect results in an underestimation of the thickness of the Barents Sea ice sheet by about 10 per cent.     

Plan-view curvature of foreland basins and its implications for the palaeostrength of the lithosphere underlying the western Alps

In zones of continental collision, three methods can be used to calculate the strength of the underthrust lithosphere: (1) a forward model approach to the Bouguer gravity field, (2) an inverse model of the gravity and topography using admittance techniques, or (3) a forward model of the stratigraphic infill of the foreland basin to estimate the cross-sectional profile of the downflexed plate. The use of reconstructed stratigraphy has the potential to yield values for the equivalent elastic thickness (T_e) of the cratonic lithosphere at varying slices in geological time, and hence enable an insight into the longer term (10–50 Myr) mechanical behaviour of the continental lithosphere. Calculations of T_e based on isopachs of foreland basin stratigraphy use sea level as a reference line to estimate the basement deflection, and therefore are limited to using stratigraphy which records shallow marine or coastal sedimentation. A new empirical approach is applied to evaluating ancient T_e values using the reconstructed palaeocurvature of the basin in plan view. The radius of curvature of 12 curvilinear foreland basins is plotted against their documented T_e values and shows a linear relationship. The maximum T_e value for a given radius of curvature can also be plotted as a straight line. The palaeocurvature of reconstructed basins can then be compared with the plots, and estimates of likely maximum T_e values may be obtained. During Eocene times, the underfilled foreland basin of the Alps was characterized on its cratonic edge by the deposition of Nummulite-rich limestones. Palaeogeographical reconstructions of the Nummulitic Limestones enable estimates of the palaeocurvature of the cratonic margin of the Alpine foreland basin during the Eocene. By comparing this value with the curvature of documented basins, it is possible to suggest that the European lithosphere underlying the western Alps had an effective elastic thickness of no greater than 17 km during the Eocene. It has been suggested that the transition in the depositional state of the Alpine foreland basin from an underfilled to a filled state during middle Oligocene times was linked to a thickening of the continental lithosphere associated with the effective ramp of the Tethyan passive margin. The T_e value of less than 17 km during the underfilled stage combined with a value of 10±5 km for the later filled stage at 17 Ma does not lend support to this hypothesis.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

Gravity studies of the Rockall and Exmouth Plateaux using SEASAT altimetry

Abstract SEASAT altimetric measurements are used to determine the gravity anomalies across two passive continental margins: the western margin of the Rockall Plateau, UK, and the Exmouth Plateau off north-west Australia. The small gravity anomalies observed over the starved western margin of the Rockall Plateau require the existence of a major density contrast within the crust, as well as the Moho, and show that the elastic thickness is less than 5 km at the time of rifting. The gravity anomaly over the Exmouth Plateau is compared with the gravity anomaly calculated from the sediment loading of a thin elastic plate, taking account of the variation in crustal thickness. This comparison shows that the Exmouth Plateau also has a small effective elastic thickness of 5 km, even for loads emplaced between 60 and 120 Myr after rifting. Elastic thicknesses of about 5 km have also been reported for other sedimentary basins, and are to be expected if the rheological properties of the crust and mantle depend on the ratio of the present temperature to the melting temperature. Flexural effects are therefore likely to be of minor importance in sedimentary basins.     

Isostasy of the Northern Bay of Biscay Continental Margin

Summary. Spectral analysis of eight marine gravity profiles and seven SEASAT profiles, combined with corresponding bathymetric data over the Northern Bay of Biscay origin, yield identical admittance functions for wavelengths greater than 120 km. the resulting admittance function has been interpreted in terms of an Airy model of compensation for wavelengths greater than 250 km and in terms of an elastic plate model of compensation for shorter wavelengths. the Airy model corresponds to a crustal thickness variation across the margin. the plate model with an elastic thickness of 8 km is associated with the regional compensation of a sedimentary load which was probably emplaced during and just after rifting.     

Sediment blanketing and the flexural strength of extended continental lithosphere

The flexural rigidity of the oceanic lithosphere is strongly dependent on its temperature structure at the time of loading. It is commonly assumed that the depth to the 450°C isotherm defines the effective elastic thickness T_e of the lithosphere. However, recent gravity studies across the Baltimore Canyon and Nova Scotian margins suggest that temperature may play a more complicated role in controlling the mechanical strength of extended continental lithosphere. For example, the flexural strength of the Baltimore Canyon margin (with sediment thicknesses of ? 15 km) appears to be controlled by the depth to the 150°C isotherm whereas the strength of the Nova Scotian margin (with sediment thicknesses cf ? 10 km) is controlled by the depth to the 250°C isotherm. The apparent correlation between sediment thickness and controlling isotherm suggests that sediment blanketing may play a role in modifying the flexural strength of extended continental lithosphere. This hypothesis was investigated by simulating the sedimentation history of a margin as a Gaussian function in which sedimentation peak and rate are determined by the mean and standard deviation of the function. The temperature structure of the lithosphere is continually modified as sediments are deposited on, and incorporated into the temperature structure of, the underlying lithosphere. Given a ‘starting’ value of T_e defined by the degree of extension of the lithosphere, the modification of T_e appears to be directly proportional to the sedimentation rate and cumulative sediment thickness, and inversely proportional to the time at which the sedimentation rate is a maximum. The first-order consequence of sediment blanketing is to reduce the cooling rate of the lithosphere relative to cooling in the absence of sediments. At thermal equilibrium, the initial value of T_e is reduced by the cumulative sediment thickness. Local isostatic conditions (i. e. T_e? 0) can only be approached when the sedimentation rate is unrealistically high (> 1000 m/Myr) during the rift or early post-rift phase of basin development. However, while these early loads may be locally compensated, any subsequent loads will be regionally compensated. Thus, it is unlikely that the low present-day flexural strengths interpreted from the Baltimore Canyon and Nova Scotian passive continental margins are a consequence of sediment blanketing.     

Subsurface mass redistributions at Mount Etna (Italy) during the 1995–1996 explosive activity detected by microgravity studies

Gravity changes are presented from a series of field microgravity surveys conducted at Mt Etna between August 1994 and November 1996, a period including the 1995–1996 explosive summit activity. Data were collected along a microgravity network of 69 stations at a monthly to annual sampling rate, depending on each subarray of the network.
  Results show that seasonal changes in water level within the volcano may induce gravity changes of up to 20  μgal on Etna's southern slope, and indicate that significant magma movement occurred within and below Etna's edifice between 1994 and 1996. In particular, between September 1994 and October 1995, a mass increase of 2 × 10¹⁰  kg occurred 2000  m beneath the summit craters. Between October 1995 and July 1996 this mass was lost, while another 2 × 10¹⁰  kg was injected at about 1000  m  a.s.l. into the 1989 fracture system. From the gravity data alone, it is not possible to distinguish whether the first shallow intrusion (1994–1995) was then injected laterally into the 1989 fracture, or summit activity was fed by the first shallow intrusion, while new magma entered the 1989 fracture system.
  While magma was being redistributed within the volcanic edifice, measurements along an E–W-trending profile on the southern slope of the volcano detected some 1.5 × 10¹¹  kg of magma accumulating 2–3  km below sea level between October 1995 and November 1996.     

嫦娥一号月球地形模型的一些应用研究

利用嫦娥一号绕月卫星激光高度计数据获取得到的新的360 阶次的地形模型CLTMs03,开展了初步应用研究,包括对嫦娥一号CCD立体相机数据进行了校对,发现并了解立体相机分析中的约900 m的垂向系统误差;利用撞击坑计数法,基于标校后的照相数据对嫦娥一号地形中新发现的玉兔火山区域的年龄特征进行了分析;结合新的CLTM-s03 地形模型和CEGM-02 重力场模型,新证实了西拉德撞击坑是一个布格质量瘤异常区域。     

Thermal origin of mid-plate hot-spot swells

Summary. Additional evidence supports the idea that the shallow rises surrounding mid-plate, hot-spot volcanoes are caused by a broad-scale reheating of the lithosphere above hot-spots. Firstly, as required by the reheating concept, the rises appear to be supported by a density deficiency within the normal thickness of the lithosphere. The gravity anomalies over the Bermuda, Cape Verde, Hawaii and Cook-Austral swells indicate that the compensation of these swells is only 40 to 100 km deep. The geoid anomaly over the Hawaiian swell is consistent with these depths. Secondly, as also required by the reheating concept; swells and the volcanoes formed on swells subside at the same rate as younger, hotter lithosphere which is at the same ocean depth. Almost all mid-plate swells rise to an ocean depth of 4250 m, the depth of normal 25 Myr-old lithosphere. The Hawaiian Swell, Emperor Guyots, Cook-Austral Swell and Bikini and Enewetok Atolls all subside as 25 Myr-old lithosphere subsides.     

Physical controls and resulting morphological forms of Quaternary ice-contact volcanoes in western Canada

Volcanism in association with large quantities of ice manifests itself in a variety of morphological forms developed under differing physical conditions. These physical conditions include the location, amount, thickness, and type of confining ice, the location and quantity of trapped water, and the surrounding topography. Ice in the form of thick continental ice sheets, thinner alpine glaciers or even perennial accumulations of snow can influence the resulting morphological form of the volcano. The chemical composition of the erupting magma, effusion rate and total erupted volume will also affect the shape of the resulting edifice. Quaternary volcanoes influenced by many of these factors are found in western Canada. The morphological forms range from tuyas to subglacial mounds to stratovolcanoes. This paper provides a review of ice-contact volcanism and volcanic products formed in these environments using Canadian examples for illustration.