In this work, many global tomographic inversions and resolution tests are carried out to investigate the influence of various mantle and core phase data from the International Seismological Center (ISC) data set on the determination of 3D velocity structure of the Earth's interior. Our results show that, when only the direct P data are used, the resolution is good for most of the mantle except for the oceanic regions down to about 1000 km depth and for most of the D″ layer, and PP rays can provide a better constraint on the structure down to the middle mantle, in particular for the upper mantle under the oceans. PcP can enhance the ray sampling of the middle and lower mantle around the Pacific rim and Europe, while Pdiff can help improve the spatial resolution in the lowermost mantle. The outer core phases (PKP, PKiKP and PKKP) can improve the resolution in the lowermost mantle of the southern hemisphere and under oceanic regions. When finer blocks or grid nodes are adopted to determine a high-resolution model, pP data are very useful for improving the upper mantle structure. The resulting model inferred from all phases not only displays the general features contained in the previous global tomographic models, but also reveals some new features. For example, the image of the Hawaiian mantle plume is improved notably over the previous studies. It is imaged as a continuous low velocity anomaly beneath the Hawaiian hotspot from the core-mantle boundary (CMB) to the surface, implying that the Hawaiian mantle plume indeed originates from the CMB. Low-velocity anomalies along some mid-oceanic ridges extend down to about 600 km depth. Our results suggested that later seismic phases are of great importance in better understanding the structure and dynamics of the Earth's interior. 相似文献
Although many geochemical, geophysical and seismological studies have suggested that the Hawaiian mantle plume originates from the core–mantle boundary (CMB), so far no tomographic model shows a continuous image of the Hawaiian plume in the entire mantle because of the few seismic stations on the narrow Hawaiian island chain. Here we present a new tomographic image beneath Hawaii determined by using simultaneously 10 kinds of seismic phases, P, pP, PP, PcP, Pdiff, PKPab, PKPbc, PKiKP, PKKPab and PKKPbc, extracted from the data set compiled by the International Seismological Center. Of these phases, PKiKP, PKKPab and PKKPbc are, for the first time, attempted to use in the global seismic tomography. Our results show a slow anomaly beneath Hawaii ascending continuously from the CMB to the surface, implying that the Hawaiian plume indeed originates from the CMB. This image is improved notably over the previous results in the whole mantle, particularly in and below the middle mantle, suggesting that later phases, PP, Pdiff, PKP and particularly PKiKP, are of great importance for better imaging the Hawaiian plume. This slow anomaly is considered to be a plume conduit being tilted, which is likely caused by the mantle flow. This indicates that the position of the Hawaiian hotspot on the surface is not stationary, as evidenced by the recent paleomagnetic and numerical modeling studies. 相似文献
Precursor and coda portions of short-period PcP waves (reflected P wave from the core-mantle boundary, CMB) recorded at J-array stations in Japan were analyzed in order to extract weak scattered signals originating from small-scale heterogeneities in the lowermost mantle beneath northeastern China. Two nuclear explosions at Lop Nor in China detonated on 21 May 1992 (Mb=6.5) and 8 June 1996 (Mb=5.9) were used for our analysis.Three-dimensional grids above the CMB were defined in the area around the PcP bounce points beneath northeastern China to calculate theoretical travel times of scattered waves which propagate from the sources to each grid point and arrive at each station based on the IASP91 model. Subsequently the waveforms were aligned with respect to the theoretical travel times and the semblance (an amplitude dependent measure of coherency) was calculated for each grid point. In order to obtain a more accurate travel time correction, we applied a cross correlation method to PcP waveforms in order to reduce picking error of the PcP onset time. A cross convolution method was also applied so that the two events could be analyzed simultaneously without using unstable deconvolutions.We could identify regions with relative high semblance values in semblance contour maps at about 200 and 375 km above the CMB. Stacking waveforms with respect to the theoretical travel times for the grid points with relative high semblance values indicate coherent wavelets originating at those grid points, that is, they correspond to scattered waves originating from small-scale heterogeneities in the lowermost mantle. Our results indicate the existence of small-scale scattering objects in the D″ layer, especially in the depth range of 200 and 375 km above the CMB beneath northeastern China. Considering recent tomographic images of high velocity anomalies in this area, these scattering objects could be fragments of old oceanic crusts which have subducted through the lower mantle and have accumulated in the D″ layer beneath northeastern China. 相似文献
A three-dimensional (3D) ray tracing technique is used to investigate ray path variations of P, PcP, pP and PP phases in a global tomographic model with P wave velocity changing in three dimensions and with lateral depth variations of the Moho, 410 and 660 km discontinuities. The results show that ray paths in the 3D velocity model deviate considerably from those in the average 1D model. For a PcP wave in Western Pacific to East Asia where the high-velocity (1-2%) Pacific slab is subducting beneath the Eurasian continent, the ray path change amounts to 27 km. For a PcP ray in South Pacific where very slow (−2%) velocity anomalies (the Pacific superplume) exist in the whole mantle, the maximum ray path deviation amounts to 77 km. Ray paths of other phases (P, pP, PP) are also displaced by tens of kilometers. Changes in travel time are as large as 3.9 s. These results suggest that although the maximal velocity anomalies of the global tomographic model are only 1-2%, rays passing through regions with strong lateral heterogeneity (in velocity and/or discontinuity topography) can have significant deviations from those in a 1D model because rays have very long trajectories in the global case. If the blocks or grid nodes adopted for inversion are relatively large (3-5°) and only a low-resolution 3D model is estimated, 1D ray tracing may be feasible. But if fine blocks or grid nodes are used to determine a high-resolution model, 3D ray tracing becomes necessary and important for the global tomography. 相似文献
We investigate the utility of PKP coda waves for studying weak scattering from small-scale heterogeneity in the mid-mantle. Coda waves are potentially a useful probe of heterogeneity in the mid-mantle because they are not preferentially scattered near the CMB, as PKP precursors are, but are sensitive to scattering at all depths. PKP coda waves have not been used for this purpose historically because of interference with other late-arriving energy due to near-surface resonance and scattering. Any study of deep mantle scattering using coda waves requires the removal of near-surface effects from the data. We have analyzed 3624 recordings of PKP precursors and coda made by stations in the Incorporated Research Institutions for Seismology (IRIS) Global Seismographic Network (GSN). To study the range and time dependence of the scattered waves, we binned and stacked envelopes of the recordings. We have considered precursors that arrive within a 20 s window before PKP and coda waves in a 60 s window after PKP. The PKP scattered waves increase in amplitude rapidly with range as predicted by scattering theory. At ranges below 125°, we predict and observe essentially no scattered energy preceding PKP. Coda amplitudes at these ranges are independent of range and provide an estimate of energy due to near-surface effects that we can expect at all ranges. We use the average coda amplitude at ranges from 120 to 125° to correct coda amplitudes at other ranges. PKP coda waves show a strong dependence on time and range and are clearly influenced by scattering in the lower mantle. PKP coda waves, however, do not provide a tighter constraint on the vertical distribution of mantle heterogeneity than is provided by precursors. This is due, in part, to relatively large scatter in coda amplitudes as revealed by a resampling analysis. Modeling using Rayleigh–Born scattering theory and an exponential autocorrelation function shows that PKP coda amplitudes are not highly sensitive to the vertical distribution of heterogeneity in the mantle. To illustrate this we consider single-scattering in two extreme models of mantle heterogeneity. One allows heterogeneity just at the CMB; the other includes heterogeneity throughout the mantle. The amplitudes of precursors are tightly constrained by our stack and support our earlier conclusion that small-scale heterogeneity is uniformly distributed throughout the lower mantle. The best-fit model includes 8 km scale length heterogeneity with an rms velocity contrast throughout the mantle of 1%. 相似文献
I have examined precisely the differential travel times and waveforms of SmKS seismic phases propagating under the southern Indian Ocean obtained from African broadband seismic arrays. The SmKS phases analyzed in this study travel in the mantle with weak heterogeneity confirmed by a global tomographic study for the distance range of 115-135°. The SmKS differential times were obtained from a vespagram (a stack intensity on a time-slowness diagram), and comparison with the vespagram created from synthetic waveforms with PREM gives the travel-time residual for each event-array pair. Although the residuals of S3KS-S2KS times exhibit apparently a systematic dependence on epicentral distance, this is likely due to small-scale heterogeneity beneath the Oceania where is covered by the SmKS ray entering points at the CMB. Waveform modeling was applied to a record section with a small travel-time residual that suggests a small effect from the mantle heterogeneity on the data set, I found that a low-velocity zone in the outermost 50 km in the core rather than PREM can explain an additional arrival detected just after the S3KS phase. This result is still inconclusive because of the small number of data and non-uniqueness of the model and ambiguity due to mantle structure. However, accumulation of the precise measurement described in this study may help the reduction of uncertainty and trade-offs. 相似文献
The differential axial and equatorial rotations of both cores associated with the Quaternary glacial cycles were evaluated based on a realistic earth model in density and elastic structures. The rheological model is composed of compressible Maxwell viscoelastic mantle, inviscid outer core and incompressible Maxwell viscoelastic inner core. The present study is, however, preliminary because I assume a rigid rotation for the fluid outer core. In models with no frictional torques at the boundaries of the outer core, the maximum magnitude of the predicted axial rotations of the outer and inner cores amounts to ∼2° year−1 and ∼1° year−1, respectively, but that for the secular equatorial rotations of both cores is ∼0.0001° at most. However, oscillating parts with a period of ∼225 years are predicted in the equatorial rotations for both cores. Then, I evaluated the differential rotations by adopting a time-dependent electromagnetic (EM) torque as a possible coupling mechanism at the core-mantle boundary (CMB) and inner core boundary (ICB). In a realistic radial magnetic field at the CMB estimated from surface magnetic field, the axial and equatorial rotations couple through frictional torques at the CMB, although these rotations decouple for dipole magnetic field model. The differential rotations were evaluated for conductivity models with a conductance of 108 S of the lowermost mantle inferred from studies of nutation and precession of the Earth and decadal variations of length of day (LOD). The secular parts of equatorial rotations are less sensitive to these parameters, but the magnitude for the axial rotations is much smaller than for frictionless model. These models, however, produce oscillating parts in the equatorial rotations of both cores and also in the axial rotations of the whole Earth and outer and inner cores. These oscillations are sensitive to both the magnitude of radial magnetic field at the CMB and the conductivity structure. No sharp isolated spectral peaks are predicted for models with a thin conductive layer (∼200 m) at the bottom of the mantle. In models with a conductive layer of ∼100 km thickness, however, sharp spectral peaks are predicted at periods of ∼225 and ∼25 years for equatorial and axial rotations, respectively, although these depend on the strength of radial magnetic field at the CMB. While the present study is preliminary in modelling the fluid outer core and coupling mechanism at the CMB, the predicted axial rotations of the whole Earth may be important in explaining the observed LOD through interaction between the equatorial and axial rotations. 相似文献
The attenuation factor QP at the top of the inner core is evaluated by using the amplitude spectral ratio of PKPdf and PKPbc phases observed at African
stations (BGCA mostly), from strong deep earthquakes in the Pacific Ocean area. The maximum depth of penetration of the PKPdf
phase into the inner core (IC) is roughly 377 km, and the sampled region of IC is centered beneath the Southern Indian Ocean.
The derived mean value of QP is 249 ± 31 (95% confidence level) in the frequency range 0.2–2 Hz, where no frequency dependence of attenuation has been
reliably observed. By using Student’s t-test, we show that the value is statistically significantly different (with a probability
greater than 95%) from other mean values of Q derived by using the same method, for both the western (180 °W to 40 °E) and
eastern (40 °E to 180 °E) hemispheres of the IC. The decrease of Q with the radius of the turning point (denoted by rTP), according to QP = 840 − 0.62 rTP, has a moderate statistical support (the R-squared value is 38%). A slightly increase of Q as a function of the angle of
the PKPdf path within the inner core with respect to the Earth’s spin axis is observed, in agreement with various investigations
performed in the time domain. However, the value of the anisotropy, if any, is suggested to be around 3%. 相似文献
Scattering by a slightly-rough core-mantle boundary (CMB) with small-scale radial variations of up to a few hundred metres, has been an attractive (though non-unique) interpretation of at least part of the precursors to PKIKP. Here it is shown that a slightly-rough CMB has an observable effect on PKKP as well, if the signal-to-noise ratio is sufficiently high. The effect may be observed as precursive arrivals and is due to back-scattering at CMB. This work was prompted by observations by Chang and Cleary at LASA of “PKKP” and precursors from the Novaya Zemlya explosions. NORSAR data from several source regions are presented here; small-scale radial variations of 100–200 metres are inferred from these data, although in some regions the CMB appears to be much smoother. On the other hand, the LASA data are anomalous and suggest much larger topography in the sampled region of the CMB. Both large- and small-scale topography must be dynamically produced, if current estimates of the viscosity of the lower mantle (～1022 Poise) are correct. 相似文献
Indirect observations and theoretical predictions for the period of the free core nutation (FCN) differ by anywhere from 15 to 30 days, and various effects have been invoked in attempts to explain this difference. The favored explanation remains as much as 5% departure in the flattening of the core-mantle boundary (CMB) from that of its hydrostatic reference figure. This 5% ‘extra-flattening’ of the CMB is not seen at the Earth's surface, where the difference is only about 0.5%. In contrast to the a posteriori model adjustments used to determine this up to 5% value, and the kinematic results available from viscous flow modeling using the seismically determined lateral heterogeneity in density data, we consider this problem from the perspective of a forward-modeling dynamical study. More specifically, we investigate the related problem of flow-induced surface and CMB topography, arising from convection in the mantle. As such, we have completed a comparative and systematic study of relative surface and CMB topography resulting from numerical models of mantle convection. When effects resulting from boundary curvature are isolated, it appears that the magnitude of CMB topography produced is insufficient in producing a significant extra-flattening of the CMB. However, results concerning effects solely resulting from a depth-dependent mantle viscosity profile, indicate that this factor may indeed lead to enhanced topography at the CMB of the magnitude required to produce the extra-flattening there. 相似文献
Computing synthetic seismograms for media with localized heterogeneous regions can be performed using hybrid methods. Here, a combination of a finite-difference (FD) technique and a frequency-wavenumber (ω − k) filtering is applied to model wave reflection at different kinds of core-mantle boundary (CMB) topography. The FD method is only applied in the neighbourhood of the CMB, while the ω − k filter is used to continue the reflected wavefield to the Earth's surface. Synthetic SH-seismograms for ScS with a dominant frequency of 0.5 Hz are computed at epicentral distances from 44° to 69°. The topography varies in amplitude (maximum amplitude of 1.0–2.7 km) and in its wavenumber spectrum; it is either monochromatic (wavelengths from 55 to 270 km) or statistical (coloured noise). The seismograms for a CMB with topography are compared with those for a plane CMB. We observe that monochromatic topography with short wavelengths (less than 100 km) results in amplitude reduction and shorter travel times than in the case of a plane CMB, but no variations with epicentral distance appear, whereas greater wavelengths exhibit amplitude variations with distance as well as travel time residuals, which both correlate with the CMB topography. Statistical models show amplitude variations with epicentral distance, while the travel time residuals are very small (less than 0.1 s). All synthetics illustrate that wavefront healing occurs along the ray path from the CMB to the Earth's surface. While the seismograms at the CMB exhibit strong fluctuations, the fluctuations at the surface are smoothed and reduced. This demonstrates that it is necessary to use wave theoretical methods for computing synthetic seismograms for complicated structures at greater depth. It also follows that travel times are less sensitive to the structure than the amplitudes. 相似文献
We compare lateral variations at the base of the mantle as inferred from a global dataset of PcP-P travel time residuals, measured on broadband records, and existing P and S tomographic velocity models, as well as ScS-S travel time data in some selected regions. In many regions, the PcP-P dataset implies short scale lateral variations that are not resolved by global tomographic models, except under eastern Eurasia, where data and models describe a broad region of fast velocity anomalies across which variations appear to be of thermal origin. In other regions, such as central America and southeastern Africa, correlated short scale lateral variations (several hundred kilometers) are observed in PcP and ScS, implying large but not excessive values for the ratio R=∂ ln Vs/∂ ln Vp (∼2.5). On the other hand, in at least two instances, in the heart of the African Plume and on the edge of the Pacific Plume, variations in P and S velocities appear to be incompatible, implying strong lateral gradients across compositionally different domains, possibly also involving topography on the core-mantle boundary. One should be cautious in estimating R at the base of the mantle from global datasets, as different smoothing and sampling of P and S datasets may result in strong biases and meaningless results. 相似文献
Deep earthquakes located in the Tonga-Kermadec region produce exceptionally clear and sharp short-period P, S, PcP, ScP, and ScS phases which are recorded at many stations at distances of less than 60°. The data used in this study are produced by short-period stations located in oceanic-type regions (Fiji and New Caledonia), a mobile continental region (eastern Australia) and a shield region (central Australia). Differential travel-time residuals of the above phases at these stations are investigated to determine the contribution to the differential residuals from: (1) the upper part of the mantle (S-P residuals); (2) the core-to-station portion of the mantle (ScS-ScP residuals); and (3) the hypocenter-to core portion of the mantle (ScP-PcP residuals). The use of differential travel-time residuals considerably reduces near-station effects and effects due to inaccurate determination of the source parameters, and hence the results can be interpreted as due to variations along the propagation paths. The results show that (S-P) residuals from phases traveling along event-to-station paths are about 7 s smaller at the shield station than at the oceanic stations. This correlation with surface tectonic environments is equally strong for the (ScS-ScP) residuals, with the shield/oceanic station difference being about 4 s. Moreover, the data suggest that this correlation between differential residuals and surface tectonic environments is caused by variations in shear velocity within the upper part of the mantle. However, the data cannot uniquely resolve the required depth of these variations within the mantle. For example, if the shear velocity variations extend to a depth of 400 km beneath the recording stations, then the average shear velocity difference between shield- and oceanic-type environments is about 4%. However, if the variations extend only to a depth of 200 km, this difference is more than 8%.(ScP-PcP) and (ScS-PcS) residuals vary from about +1 to about +4 s at the different stations, apparently because of compressional velocity variations in the mantle along the Pc path. If the variation in compressional velocity within the mantle below a depth of about 600 km is about 10% and occurs near the source region, these results suggest that, in the vicinity of deep earthquake zones, variations in compressional velocity extend to a depth of about 1000 km. However, these results can equally be explained by a 1% variation in compressional velocity, evenly distributed along the entire Pc path. An estimate of Q determined from the observed predominant frequency of ScS waves, as recorded at the shield station, suggests that the average 〈Qs〉 of the mantle beneath about 600 km is about 1050 at frequencies of about 1 Hz. 相似文献
An ScP phase reflected and converted at the core–mantle boundary (CMB) beneath the region east of the Philippine Islands shows clear pre- and postcursors, recorded on short-period seismic networks in Japan. These waveform variations can be explained by interaction of the ScP wavefield with thin layers at the CMB. The results of forward modeling of double-array stacks reveal two different structural heterogeneities in the lowermost mantle beneath the region east of the Philippine Islands. One of the structures represents a decreased velocity, and increased density across the reflector at the lowermost ～10 km of the mantle, with P- and S-wave velocity reductions of 5–10% and ～30%, respectively, and an increase in density of 5–10%. Another structure consists of a pair of reflectors at ～10 km and ～5 km above the CMB, both of which are characterized by reduced P- and S-wave velocities. The upper reflector is the interface of a low-velocity zone in which P- and S-wave velocities decrease of 10% and 30%, respectively, accompanied by an extremely large increase in density (20–25%). The lower reflector is characterized by a 25% reduction in S-wave velocity relative to the above low-velocity layer, as well as a 5% decrease in P-wave velocity and no change in density. The nature of the low-velocity zone detected locally at the CMB is comparable with that of ultra-low-velocity zones (ULVZs) observed by various seismic probes in the South Pacific and Central America. Extensive observations of the ULVZ beneath the region east of the Philippine Islands indicate massive partial melting at the bottom of the mantle. Low-S-velocity basal layer partly detected within the ULVZ may be resulting from core–mantle chemical interactions, driven by massive partial melting. 相似文献
The origin of large low shear-wave velocity provinces (LLSVPs) in the lowermost mantle beneath the central Pacific and Africa is not well constrained. We explore numerical convection calculations for two proposed hypotheses for these anomalies, namely, thermal upwellings (e.g., plume clusters) and large intrinsically dense piles of mantle material (e.g., thermochemical piles), each of which uniquely affects the topography on Earth's core–mantle boundary (CMB). The thermochemical pile models predict a relatively flat but elevated CMB beneath piles (presumed LLSVPs), with strong upwarping along LLSVP margins. The plume cluster models predict CMB upwarping beneath upwellings that are less geographically organized. Both models display CMB depressions beneath subduction related downwelling. While each of the two models produces a unique, characteristic style of CMB topography, we find that seismic models will require shorter length scales than are currently being employed in order to distinguish between the end-member dynamic models presented here. 相似文献
In 1983, Lay and Helmberger [Geophys. J. R. Astron. Soc. 75 (1983) 799–837] reported the detection of a precursor to the seismic phase ScS. They attributed this precursor to a sharp seismic discontinuity located several hundred kilometers above the core–mantle boundary. Such a lowermost mantle discontinuity implies the existence of a sharp phase change or a chemical boundary. Precursors to ScS and, less frequently, PcP have since been observed in numerous locations, but are not a global phenomenon. Frequently, PcP precursors are weak or absent when ScS precursors are observed in the same location, and vice versa. There can be significant variations in the amplitude and arrival time of the precursor relative to the main phase. The presence or absence of these precursors has led to speculations about the nature of the lowermost mantle. Here we demonstrate that ScS or PcP precursors may be produced by gradients in seismic wave speed associated with large-scale lowermost mantle heterogeneity. Rather than a phase or chemical boundary with substantial topography, such gradients require lateral variations in temperature and, close to the core–mantle boundary, composition. 相似文献
The derivation of P and S velocities at the core-mantle boundary (CMB) from long-period diffracted waves by the use of the simple ray-theoretical formulavCMB=rc/p (vCMB=velocity at the CMB;rc=core radius;p=ray parameter) yields apparent velocity values which differ from the true velocities. Using a dominant period of about 20 sec for calculating theoretical seismograms, we found a linear relation between the apparent velocity and the average velocity in a transition zone at the base of the mantle with fixed velocity on top.The ray parameters determined from long-period earthquake data are found to be 4.540±0.035 and 8.427±0.072 sec/deg for Pdiff and Sdiff, respectively. These values yield apparent velocities of 13.378±0.103 for P and 7.207±0.062 km/sec for S waves. By means of the theoretical relation between apparent and average velocity and under the assumption of linear variation of velocity with depth, one can invert the apparent velocities into true CMB velocities of 13.736±0.170 and 7.320±0.124 km/sec. These results imply positive velocity gradients at the base of the mantle and hence no significant departures from adiabaticity and homogeneity.Contribution No. 211 of the Geophysical Institute, University of Karlsruhe. 相似文献
The amplitudes of the core reflection PcP are sensitive to the wave velocities and densities in the neighborhood of the core-mantle boundary (CMB). We study the amplitude ratio of the long-period phases PcP and P from two South American deep-focus earthquakes with favorable fault-plane solution, depth and magnitude, as recorded by WWNSS and CSN stations in North America.Comparison is made with long-period PcP/P amplitude ratios, derived from theoretical seismograms for a variety of CMB models. Models from previous studies, which were mainly derived from short-period PcP observations and which are characterized by discrete layers above the CMB, are almost all inconsistent with the long-period data. The data also discriminate against low nonzero S velocities below the CMB. Simple first-order-discontinuity models of the CMB, for instance according to the Jeffreys-Bullen earth model or according to recent models based mainly on free oscillations, explain the data reasonably well.Model improvements are attempted by varying the P-velocity gradient above the CMB. The best amplitude fit is obtained for a rather strong decrease in P velocity with depth in this zone which, however, gives no acceptable traveltime fit for PcP. The scatter in body-wave amplitudes is considerable even for long-period waves and may prevent the correct assessment of that part of the amplitude variation of a phase with distance that is due to the variation of velocities and densities with depth alone. 相似文献
The ∼0.2 mm/yr uplift of Hawaiian islands Lanai and Molokai and Hawaiian swell topography pose important constraints on the structure and dynamics of mantle plumes. We have formulated 3-D models of mantle convection to investigate the effects of plume-plate interactions on surface vertical motions and swell topography. In our models, the controlling parameters are plume radius, excess plume temperature, and upper mantle viscosity. We have found that swell height and swell width constraints limit the radius of the Hawaiian plume to be smaller than 70 km. The additional constraint from the uplift at Lanai requires excess plume temperature to be greater than 400 K. If excess plume temperature is 400 K, models with plume radius between 50 and 70 km and upper mantle viscosity between 1020 and 3×1020 Pa s satisfy all the constraints. Our results indicate that mantle plume in the upper mantle may be significantly hotter than previously suggested. This has important implications for mantle convection and mantle melting. In addition to constraining plume dynamics, our models also provide a mechanism to produce the observed uplift at Lanai and Molokai that has never been satisfactorily explained before. 相似文献
