首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 406 毫秒
1.
The elasticity of ferropericlase with a potential mantle composition of (Mg0.83,Fe0.17)O is determined using ultrasonic interferometry in conjunction with in situ X-radiation techniques (X-ray diffraction and X-radiography) in a DIA-type cubic anvil high-pressure apparatus to pressures of 9 GPa (NaCl pressure scale) at room temperature. In this study, we demonstrate that it is possible to directly monitor the specimen length using an X-ray image technique and show that these lengths are consistent with those derived from X-ray diffraction data when no plastic deformation of the specimen occurs during the experiment. By combining the ultrasonic and X-ray diffraction data, the adiabatic elastic bulk (KS) and shear (G) moduli and specimen volume can be measured simultaneously. This enables pressure scale-free measurements of the equation of state of the specimen using a parameterization such as the Birch-Murnaghan equation of state. The elastic moduli determined for (Mg0.83,Fe0.17)O are KS0=165.5(12) GPa, G0=112.4(4) GPa, and their pressure derivatives are KS0′=4.17(20) and G0′=1.89(6). If these results are compared with those for MgO, they demonstrate that KS0 and KS0′ are insensitive to the addition of 17 mol% FeO, but G0 and G0′ are reduced by 14% and 24%, respectively. We calculate that the P and S wave velocities of a perovskite plus ferropericlase phase assemblage with a pyrolite composition at the top of the lower mantle (660 km depth) are lowered by 0.8 and 2.3%, respectively, when compared with those calculated using the elastic properties of end-member MgO. Consequently, the magnitudes of the calculated wave velocity jumps across the 660 km discontinuity are reduced by about 11% for P wave and 20% for S wave, if this discontinuity is considered as a phase transformation boundary only (ringwoodite→perovskite+ferropericlase).  相似文献   

2.
Using acoustic measurement interfaced with a large volume multi-anvil apparatus in conjunction with in situ X-radiation techniques, we are able to measure the density and elastic wave velocities (VP and VS) for both ortho- and high-pressure clino-MgSiO3 polymorphs in the same experimental run. The elastic bulk and shear moduli of the unquenchable high-pressure clinoenstatite phase were measured within its stability field for the first time. The measured density contrast associated with the phase transition OEN → HP-CEN is 2.6-2.9% in the pressure of 7-9 GPa, and the corresponding velocity jumps are 3-4% for P waves and 5-6% for S waves. The elastic moduli of the HP-CEN phase are KS=156.7(8) GPa, G = 98.5(4) GPa and their pressure derivatives are KS′=5.5(3) and G′ = 1.5(1) at a pressure of 6.5 GPa, room temperature. In addition, we observed anomalous elastic behavior in orthoenstatite at pressure above 9 GPa at room temperature. Both elastic wave velocities exhibited softening between 9 and 13-14 GPa, which we suggest is associated with a transition to a metastable phase intermediate between OEN and HP-CEN.  相似文献   

3.
In situ synchrotron X-ray diffraction experiments were conducted using the SPEED-1500 multi-anvil press at SPring-8 on majoritic garnet synthesized from natural mid-ocean ridge basalt (MORB), whose chemical composition is close to the average of oceanic crust, at 19 GPa and 2200 K. Pressure-volume-temperature data were collected using a newly developed high-pressure cell assembly to 21 GPa and 1273 K. Data were fit to the high-temperature Birch-Murnaghan equation of state, with fixed values for the ambient cell volume (V0 = 1574.14(4) Å3) and the pressure derivative of the isothermal bulk modulus (K′T = 4). This yielded an isothermal bulk modulus of KT0 = 173(1) GPa, a temperature derivative of the bulk modulus (∂KT/∂T)P = −0.022(5) GPa K−1, and a volumetric coefficient of thermal expansivity α = a + bT with values of a = 2.0(3) × 10−5 K−1 and b = 1.0(5) × 10−8 K−2. The derived thermoelastic parameters are very similar to those of pyrope. The density of subducted oceanic crust compared to pyrolitic mantle at the conditions in Earth's transition zone (410-660 km depth) was calculated using these results and previously reported thermoelastic parameters for MORB and pyrolite mineral assembledges. These calculations show that oceanic crust is denser than pyrolitic mantle throughout the mantle transition zone along a normal geotherm, and the density difference is insensitive to temperature at the pressures in lower part of the transition zone.  相似文献   

4.
First-principles calculations have been used to determine the equation of state of Fe3C in both its low-pressure magnetically ordered and high-pressure non-magnetically ordered states; at 0 K the ferromagnetic transition was found to occur at about 60 GPa. In the high pressure, non-magnetically ordered regime at 0 K the material may be described by a Birch-Murnaghan third-order equation of state with V0=8.968(7) Å3 per atom, K0=316.62(2) GPa and K′=4.30(2). At atmospheric pressure the ferromagnetic phase transition in Fe3C occurs at ∼483 K; preliminary measurements of the thermal expansion by powder neutron diffraction show that this transition produces a large effect on thermoelastic properties. The volumetric thermal expansion coefficient in the paramagnetic phase was found to be 4.34×10−5 K−1 at T∼550 K. By applying a thermal expansion correction to the calculated equation of state at 0 K, predicted values for the density and adiabatic incompressibility of this material at core pressures and temperatures were obtained. These results appear to be sufficiently different from seismological data so as to preclude Fe3C as the major inner core-forming phase.  相似文献   

5.
We present new elasticity measurements on single-crystal fayalite and combine our results with other data from resonance, pulse superposition interferometry, and Brillouin scattering to provide a set of recommended values for the adiabatic elastic moduliC ij and their temperature variations. We use a resonance method (RPR) with specimens that were previously investigated by pulse superposition experiments. The nineC ij of fayalite are determined from three new sets of measurements. One set of our newC ij data is over the range 300–500 K. We believe that the relatively large discrep ncies found in someC ij are due in large part to specimen inhomogeneities (chemical and microstructural) coupled with differences in the way various techniques sample, rather than only systematic errors associated with experimental procedures or in the preparations of the specimens.Our recommendeaC ij's (GPa) and (C ij/T) p (GPa/K) are: The resulting values for the isotropic bulk and shear moduli,K s and , and their temperature derivatives are:K s=134(4) GPa; =50.7(0.3) GPa; (K s/T) p =–0.024(0.005) GPa/K; and (/T) p =–0.013(0.001) GPa/K. An important conclusion is thatK s increases as the Fe/(Fe+Mg) ratio in olivine is increased.  相似文献   

6.
Ringwoodite Mg2SiO4 with spinel structure is an important phase in the earth's mantle transition zone. Controlled deformation experiments showed that ringwoodite underwent ductile deformation when compressed axially at 6-10 GPa and at room temperature in a multianvil D-DIA deformation apparatus. Texture evolution during cyclic compression has been recorded in situ using X-ray transparent anvils with monochromatic synchrotron X-ray diffraction and a two-dimensional detector. Quantitative analysis of the images with the Rietveld method revealed a 1 1 0 fiber texture. By comparing this texture pattern with polycrystal plasticity simulations, it is inferred that {1 1 1}〈 1 0〉 slip is the dominant deformation mechanism in ringwoodite, consistent with high temperature mechanisms observed in other spinel-structured materials. Although strong ringwoodite textures may develop in the transition zone, the contribution to bulk anisotropy is minimal due to the weak single-crystal anisotropy.  相似文献   

7.
Iron sulfide (FeS) was investigated using first-principles calculations up to a pressure of 400 GPa. A number of new phase transitions were found. An antiferromagnetic MnP-type structure, FeS II, was confirmed to be stable at low pressures, whereas at high pressures (40–135 GPa) we find a new stable phase, with a non-magnetic MnP-type structure, FeS VI. The observed first-order change in the cell shape between the two phases can be explained by the difference in magnetic configurations. The calculated cell parameters, atomic coordinates, and bulk modulus of non-magnetic MnP-type phase are consistent with those determined from experiment. The upper pressure limit of the stability of the non-magnetic MnP-type phase was calculated to be 135 GPa. A hitherto unsuspected phase transition from the non-magnetic MnP-type to a phase with Pmmn symmetry, FeS VII, was identified using the evolutionary crystal structure prediction (USPEX) method. The structure of the Pmmn phase has no known analogues, but can be described as a distortion of the NaCl-type structure. The Pmmn phase with the distorted NaCl-type structure is stable from 135 GPa at least up to 400 GPa. According to previous experiments and the present study, the transition sequence of FeS at low temperatures is as follows: troilite ➔ antiferromagnetic MnP-type phase ➔ monoclinic phase ➔ non-magnetic MnP-type phase ➔ Pmmn phase. The calculated volume reduction from the monoclinic to the non-magnetic MnP-type phase is 1.0% at 40 GPa, which is in good agreement with experimental observations. The calculated volume reduction from the non-magnetic MnP-type to the Pmmn phase is 3.7% at 135 GPa.  相似文献   

8.
The new hydrothermal diamond anvil cell (HDAC) has been designed for optical microscopy and X-ray diffraction at pressures up to 10 GPa and temperatures between –190°C and 1200°C. Laser light reffected from the top and bottom anvil faces and the top and bottom solid sample faces produce interference fringes that provide a very sensitive means of monitoring the volume of sample chamber and for observing volume and refractive index changes in solid samples due to transitions and reactions. Synchrotron radiation has been used to make X-ray diffraction patterns of samples under hydrothermal conditions. Individual heaters and individual thermocouples provide temperature control with an accuracy of ±0.5°C. Liquid nitrogen directly introduced into the HDAC has been used to reduce the sample temperature to –190°C. The - phase boundary of quartz has been used to calculate the transition pressures from measured transition temperatures. With this method we have redetermined 5 isochores of H2O up to 850°C and 1.2 GPa at which the solution rate of the quartz became so rapid that the quartz dissolved completely before the - transition could be observed. When silica solutions were cooled, opal spherules and rods formed.  相似文献   

9.
The evolution with pressure of the unit-cell parameters brownmillerite (Ca2Fe2O5), a stoichiometric defect perovskite structure, has been determined to a maximum pressure of 9.46 GPa, by single-crystal X-ray diffraction measurements at room temperature. Brownmillerite does not exhibit any phase transitions in this pressure range. A fit of a third-order Birch–Murnaghan equation-of-state to the P–V data yields values of KT0=127.0(5) GPa and K0′=5.99(13). Analysis of the unit-cell parameter data shows that the structure compresses anisotropically. Compressional moduli for the axes are Ka0=141(1) GPa, Kb0=118(3) GPa and Kc0=122.2(2) GPa, with Ka0′=8.9(3), Kb0′=6.2(6) and Kc0′=4. The stiffest direction (i.e. along a) coincides with the direction of the FeO4 tetrahedral chains. Comparison of these data with the elasticity systematics of Ca-perovskites shows that the presence of oxygen vacancies in the brownmillerite structure softens the structure by ∼25% and that the ordering of vacancies in the perovskite structure increases the anisotropy of compression.  相似文献   

10.
In situ X-ray diffraction measurements on a calcium aluminosilicate (CAS) phase have been carried out using a laser-heated diamond anvil cell up to a pressure of 44 GPa, employing a synchrotron radiation source. CAS is the major mineral formed from sediments subducted into the Earth's mantle. The sample was heated using a YAG laser after each pressure increment to relax the deviatoric stress in the sample. X-ray diffraction measurements were carried out at T = 300 K using an angle-dispersive technique. The pressure was calculated using an internal platinum metal pressure calibrant. The Birch–Murnaghan equation of state for the CAS phase obtained from the experimental unit cell parameters showed a density of ρ0 = 3.888 g/cm3 and a bulk modulus of K0 = 229 ± 9 GPa for K′0 = 4.7 ± 0.7. When the first pressure derivative of the bulk modulus was fixed at K′0 = 4, then the value of K0 = 239 ± 2 GPa. From the experimental compressibility, the density of the CAS phase was observed to be lower than the density of co-existing Al-bearing stishovite, calcium perovskite, calcium ferrite-type phases, and (Fe,Al)-bearing Mg-perovskite in subducted sediments in the lower mantle. Therefore, the density of subducted sediments in the lower mantle decreases with increasing mineral proportion of the CAS phase.  相似文献   

11.
The compressional behavior of the MgCr_2O_4 spinel has been investigated with the CASTEP code using density functional theory and planewave pseudopotential technique. We treated the exchange-correlation interaction by both the local density approximation(LDA) and generalized gradient approximation(GGA) with the Perdew-Burker-Ernzerhof functional. Our simulation was conducted for the pressure range of 0–19 GPa. We obtained the isothermal bulk modulus(K_T) of the MgCr_2O_4 spinel as 181.46(48) GPa(GGA; low boundary) or 216.1(11) GPa(LDA; high boundary), with its first derivative(K'_T) as 4.41(6) or 4.5(1), respectively. The oxygen parameter u is not constant but negatively correlated with P, and decreases by about 0.5–0.6% for the investigated P range. The component polyhedra have different compressibilities, increasing in the order of(O_4)_1CrO_6(O_4)_2O_6MgO_4. The Mg-O bond in the MgO_4 tetrahedron is much more compressible than the Cr-O bond in the CrO_6 octahedron.  相似文献   

12.
  相似文献   

13.
The melting curve of forsterite has been studied by static experiment up to a pressure of 15 GPa. Forsterite melts congruently at least up to 12.7 GPa. The congruent melting temperature is expressed by the Kraut-Kennedy equation in the following form: Tm(K)=2163 (1+3.0(V0 ? V)/V0), where the volume change with pressure was calculated by the Birch-Managhan equation of state with the isothermal bulk modulus K0 = 125.4 GPa and its pressure derivative K′ = 5.33. The triple point of forsterite-β-Mg2SiO4-liquid will be located at about 2600°C and 20 GPa, assuming that congruent melting persists up to the limit of the stability field of forsterite. The extrapolation of the previous melting data on enstatite and periclase indicates that the eutectic composition of the forsterite-enstatite system should shift toward the forsterite component with increasing pressure, and there is a possibility of incongruent melting of forsterite into periclase and liquid at higher pressure, although no evidence on incongruent melting has been obtained in the present experiment.  相似文献   

14.
Superposed epoch studies have been carried out in order to determine the ionospheric response at mid-latitudes to southward turnings of the interplanetary magnetic field (IMF). This is compared with the geomagnetic response, as seen in the indices Kp, AE and Dst. The solar wind, IMF and geomagnetic data used were hourly averages from the years 1967–1989 and thus cover a full 22-year cycle in the solar magnetic field. These data were divided into subsets, determined by the magnitudes of the southward turnings and the concomitant increase in solar wind pressure. The superposed epoch studies were carried out using the time of the southward turning as time zero. The response of the mid-latitude ionosphere is studied by looking at the F-layer critical frequencies, foF2, from hourly soundings by the Slough ionosonde and their deviation from the monthly median values, foF2. For the southward turnings with a change in Bz of Bz > 11.5 nT accompanied by a solar wind dynamic pressure P exceeding 5 nPa, the F region critical frequency, foF2, shows a marked decrease, reaching a minimum value about 20 h after the southward turning. This recovers to pre-event values over the subsequent 24 h, on average. The Dst index shows the classic storm-time decrease to about –60 nT. Four days later, the index has still to fully recover and is at about –25 nT. Both the Kp and AE indices show rises before the southward turnings, when the IMF is strongly northward but the solar wind dynamic pressure is enhanced. The average AE index does register a clear isolated pulse (averaging 650 nT for 2 h, compared with a background peak level of near 450 nT at these times) showing enhanced energy deposition at high latitudes in substorms but, like Kp, remains somewhat enhanced for several days, even after the average IMF has returned to zero after 1 day. This AE background decays away over several days as the Dst index recovers, indicating that there is some contamination of the currents observed at the AE stations by the continuing enhanced equatorial ring current. For data averaged over all seasons, the critical frequencies are depressed at Slough by 1.3 MHz, which is close to the lower decile of the overall distribution of foFl values. Taking 30-day periods around summer and winter solstice, the largest depression is 1.6 and 1.2 MHz, respectively. This seasonal dependence is confirmed by a similar study for a Southern Hemisphere station, Argentine Island, giving peak depressions of 1.8 MHz and 0.5 MHz for summer and winter. For the subset of turnings where Bz > 11.5 nT and P 5 nPa, the response of the geomagnetic indices is similar but smaller, while the change in foF2 has all but disappeared. This confirms that the energy deposited at high latitudes, which leads to the geomagnetic and ionospheric disturbances following a southward turning of the IMF, increases with the energy density (dynamic pressure) of the solar wind flow. The magnitude of all responses are shown to depend on Bz. At Slough, the peak depression always occurs when Slough rotates into the noon sector. The largest ionospheric response is for southward turnings seen between 15–21 UT.  相似文献   

15.
Sub-surface characterization in fractured aquifers is challenging due to the co-existence of contrasting materials namely matrix and fractures. Transient hydraulic tomography (THT) is proved to be an efficient and robust technique to estimate hydraulic (Km, Kf) and storage (Sm, Sf) properties in such complex hydrogeologic settings. However, performance of THT is governed by data quality and optimization technique used in inversion. We assessed the performance of gradient and gradient-free optimizers with THT inversion. Laboratory experiments were performed on a two-dimensional, granite rock (80 cm × 45 cm × 5 cm) with known fracture pattern. Cross-hole pumping experiments were conducted at 10 ports (located on fractures), and time-drawdown responses were monitored at 25 ports (located on matrix and fractures). Pumping ports were ranked based on weighted signal-to-noise ratio (SNR) computed at each observation port. Noise-free, good quality (SNR > 100) datasets were inverted using Levenberg–Marquardt: LM (gradient) and Nelder–Mead: NM (gradient-free) methods. All simulations were performed using a coupled simulation-optimization model. Performance of the two optimizers is evaluated by comparing model predictions with observations made at two validation ports that were not used in simulation. Both LM and NM algorithms have broadly captured the preferential flow paths (fracture network) via K and S tomograms, however LM has outperformed NM during validation ( ). Our results conclude that, while method of optimization has a trivial effect on model predictions, exclusion of low quality (SNR ≤ 100) datasets can significantly improve the model performance.  相似文献   

16.
  相似文献   

17.
In this paper we evaluate the present state of the seismic regime in Southern California using the concentration parameter of seismogenic faults (K sf ,Sobolev andZavyalov, 1981). The purpose of this work is to identify potential sites for large earthquakes during the next five or ten years. The data for this study derived from the California Institute of Technology's catalog of southern California earthquakes, and spanned the period between 1932 to June 1982. We examined events as small asM L 1.8 but used a magnitude cutoff atM L =3.3 for a detailed analysis. The size of the target earthquakes (M M ) was chosen as 5.3 and 5.8.The algorithm for calculatingK sf used here was improved over the algorithm described bySobolev andZavyalov (1981) in that it considered the seismic history of each elementary seismoactive volume. The dimensions of the elementary seismoactive volumes were 50 km×50 km and 20 km deep. We found that the mean value ofK sf within 6 months prior to the target events was 6.1±2.0 for target events withM L 5.3 and 5.41.8 for targets withM L 5.8. Seventy-three percent of the targets withM L 5.8 occurred in areas whereK sf was less than 6.1. The variance of the time between the appearance of areas with lowK sf values and the following main shocks was quite large (from a few months to ten years) so this parameter cannot be used here for accurate predictions of occurrence time.Regions where the value ofK sf was below 6.1 at the end of our data set (June, 1982) are proposed as the sites of target earthquakes during the next five to ten years. The most dangerous area is the area east of San Bernardino whereK sf values are presently between 2.9 and 3.7 and where there has been no earthquake withM L 5.3 since 1948.  相似文献   

18.
This study concludes that the elongation axis (K 1) of the ellipsoid of anisotropic magnetic susceptibility (AMS) is a suitable proxy for flow axis in ashflow tuffs. 153 oriented samples (176 specimens) were studied from 18 sites in the 1.1 Ma Tshirege member of the Bandelier Tuff. These sites are distributed around the Valles caldera at distances of 5–25 km outside of the rim.K 1 axes correlate well with postulated radial flow axes at 13 sites.K 1 also agrees with measured geological flow indicators, mainly imbricated larger clasts, at 7 sites. At 2 of the 5 sites where significant disagreement is seen between theoretical radial flow directions and measuredK 1 axes, theK 1 axes correspond well with geological flow indicators, indicating that the divergence of flow from the predicted radial flow pattern is real. Two major topographic buttresses are suggested as the cause of flow divergence for the Tshirege ash flows: the San Pedro buttress northwest of the caldera, and the San Miguel buttress in the southeast. In situK 1 axes plunge about 7° toward the source at two-thirds of the sites; therefore the plunge ofK 1 is a plausible in situ indicator for thedirection of flow. Multiple flow zones in sections of several meters thickness indicate changes of flow direction that are both rapid and large during ash-flow emplacement. These observations raisre the question of how best to represent mean flow directions in ash-flow sheets: by eigenvector methods, by vector-sum methods, or by modes. A method for measuring imbrication of larger clasts using apparent dips in vertical joints is outlined. Imbrication, determined in this way at one-third of the sites, dips toward the source, i.e., up-flow. The minimum (K 3) axis of the AMS ellipsoid correlates with the flow foliation rather than with the larger clast imbrication. The flow axes of ash flows correspond with theK 1 axes, not with the declination ofK 3 axes as suggested by some authors. Initial dip of the sampled ash flows is not large and does not affect the paleomagnetic remanence direction, which is reversed with a mean ofD=173.5°,I=-38.4°, 95=3.4°N=18. This mean is not different at the 95% confidence level from that of earlier workers. The mean pole, at 098.0°E, 74.8°N,A 95=3.3°,N=18, is about 15° far-sided relative to the expected time-averaged geomagnetic pole, suggesting a history of emplacement too short to adequately average secular variation.  相似文献   

19.
We review earthquake distributions associated with hydrocarbon fields in the context of pore pressure diffusion models, poroelastic stress transfer and isostasy theory. These three mechanisms trigger or induce seismic instabilities at both local scale (D5 km) and at regional scale (D20 km). The modeled changes in stress are small (1 MPa), whatever the tectonic setting. Each mechanism corresponds to different production processes. (1) Local hydraulic fracturing due to fluid injection induces seismic-slip on cracks (M L3) within the injected reservoir through decreasing the effective stress. (2) Pure fluid withdrawal causes pore pressure to decrease within the reservoir. It triggers adjustments of the geological structure to perturbations related to the reservoir response to depletion. Poroelastic mechanisms transfer this stress change from the reservoir to the surrounding levels whereM L5 seismic instabilities occur either above or below the reservoir. (3) Massive hydrocarbon recovery induces crustal readjustments due to the removal of load from the upper crust. It can induce larger earthquakes (M L6) at greater distance from the hydrocarbon fields than the two other mechanisms.Due to the mechanical properties of the shallow rock matrices involved, seismic slip triggered either by mechanism (1) or (2), is a second-order process of the main elastoplastic deformation. for a minimum of 80% of commercially productive basins, most of the local deformation is reported as aseismic, i.e., there is no evidence forM L3 earthquakes. Nevertheless, the induced stresses vary as a function of time in a manner that depends on the hydraulic diffusivity (i.e., permeability) of the reservoir and surrounding rocks. Because small earthquakes (M L3) indicate changes in stress and pore pressure, monitoring of seismicity is a means of assessingin situ reservoir behavior.The less constrained seismic response to hydrocarbon recovery is the possible connection between local fluid manipulations, triggered earthquakes and major regional earthquakes. Positive feedback mechanisms suggest that the region of seismic hazard changes is much larger than the area where hydrocarbons are extracted. These observations and models testify that fluid movement and pore pressure changes (increase or decrease) play important roles in the mechanics of earthquakes and in the triggering of natural earthquakes.  相似文献   

20.
Experiments on MgSiO3 enstatite were conducted in the pressure range from 13 to 18 GPa under hydrous conditions in order to clarify the effect of water on the melting phase relations of enstatite at pressures corresponding to the Earth’s mantle transition zone. In some previous experiments [Geol. Soc. Am. Bull. 79 (1968) 1685; Phys. Earth Planet. Inter. 85 (1994) 237], incongruent melting behavior to form Mg2SiO4 forsterite and SiO2 enriched liquid up to 5 GPa was observed, and congruent melting behavior at pressures up to 12 GPa was observed. Under hydrous conditions, we found that the melting reaction changes from congruent to incongruent at around 13.5 GPa. Liquid formed above 13.5 GPa is enriched in MgO component relative to MgSiO3 because it coexists with stishovite (SiO2). Moreover, the solidus temperature decreases drastically at around 13.5 GPa, in unison with the change in the melting reaction. The solidus temperature is about 1400 °C at 13 GPa, but approximately 900 °C at 15 GPa. Our results show that the liquidus phase changes from clinoenstatite to stishovite with increasing pressure and water content above 13.5 GPa. MgSiO3 enstatite is one of the major constituent minerals in the Earth’s mantle, and it is expected that MgO-enriched liquid will be generated in the transition zone if water is present.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号