Influence of water on the physical properties of olivine, wadsleyite, and ringwoodite

The incorporation of water in nominally anhydrous minerals plays a crucial role in many geodynamic processes and evolution of the Earth and affects the physical and chemical properties of the main constituents of the Earth’s mantle. Technological advances now allow the transport properties of minerals to be precisely measured under extreme conditions of pressure and temperature (P and T ) that closely mimic the P –T conditions throughout much of the Earth’s interior. This contribution provides an overview of the recent progress in the experimental studies on the influence of water on physical properties (i.e., diffusivity, electrical conductivity, thermal conductivity, sound velocity, and rheology) of olivine, wadsleyite, and ringwoodite together with their applications. In particular, consistency among various experimental data is investigated, discrepancies are evaluated, and confusions are clarified. With such progress in the experimental determination of transport properties of major mantle minerals, we can expect new insights into a broad range of geoscience problems. Many unresolved issues around water inside Earth require an integrated approach and concerted efforts from multiple disciplines.

Olivine [α-(Mg,Fe) 2 SiO 4 ] is one of the most abundant minerals of the Earth's upper mantle (UM), and it transforms to wadsleyite [β-(Mg,Fe) 2 SiO 4 ] and ringwoodite [γ -(Mg,Fe) 2 SiO 4 ] in the mantle transition zone (MTZ). Previous studies showed that wadsleyite and ringwoodite can incorporate up to approximately 3 wt. % H 2 O within their crystal structures (e.g., Inoue et al., 1995;Kohlstedt et al., 1996;Ohtani et al., 2004), whereas the solubility of water in olivine is limited (< 0.9 wt. % H 2 O) over all the pressure and temperature (P -T ) conditions of the UM (e.g., Kohlstedt et 40 B.-H. Zhang and Q.-K. Xia: Physical properties of olivine, wadsleyite, and ringwoodite al., 1996;Mosenfelder et al., 2006;Bali et al., 2008;Padrón-Navarta and Hermann, 2017;Fei and Katsura, 2020). A comprehensive dataset of the water content in diverse types of natural samples, including peridotite and granulite xenoliths hosted by alkali basalts, ultrahighpressure metamorphic eclogites, and Mesozoic-Cenozoic alkali basalts, has been built not only for eastern China (e.g., Xia et al., 2006Xia et al., , 2013Xia et al., , 2019Gu et al., 2019), but also for worldwide localities . This dataset provides an integrated picture of the temporal and spatial distribution of water in the lower crust and UM at a continental scale. In addition, the recently discovered ringwoodite inclusion with ∼ 1.4 wt. % H 2 O and ice-VII inclusions in diamonds (Pearson et al., 2014;Tschauner et al., 2018) clearly demonstrated a water-rich MTZ, at least locally. Considering water solubility in the crystal structure of wadsleyite and ringwoodite (Inoue et al., 1995;Kohlstedt et al., 1996), the MTZ is expected to contribute greatly to global water circulation (Karato, 2011). In this sense, the water content of the Earth's mantle is a key parameter of the Earth's water budget. Thus, the precise determination of water content in the mantle is essential for understanding various geodynamic processes in the deep Earth, such as mantle convection and water circulation (Karato, 2011).
In the past several decades, although numerous experimental data on water solubility and incorporation mechanism of water in NAMs have been accumulated (e.g., Inoue et al., 1995;Kohlstedt et al., 1996;Demouchy and Mackwell, 2006;Withers et al., 2012;Jollands et al., 2016), determining the real water content and its distribution from such information alone is difficult. Thus, the actual abundance of water in the UM and MTZ is still unknown. Measuring the physical properties (i.e., diffusivity, electrical conductivity, thermal conductivity, sound velocity, and rheology) of mantle minerals as a function of water content and combining multidisciplinary constraints from seismology, magnetotelluric survey, petrology, geochemistry, and theoretical simulation can be used to directly determine the water content and distribution in the Earth's present mantle. Although many studies have investigated the transport properties of NAMs (e.g., Jacobsen and Smyth, 2006;Yoshino, 2010;Yoshino and Katsura, 2013;Gardés et al., , 2015Karato, 2011Karato, , 2019, there are still some controversies and confusion. For example, the effect of water on the electrical conductivity of olivine is often contradictory between different groups (D. Yoshino et al., 2006;Karato and Dai, 2009;, 2013Gardés et al., 2015;Dai and Karato, 2015). With respect to the causes of these experimental discrepancies or controversies, one has to fingerprint the type of defect that is involved in each process (from fast electrical conductivity properties to extremely slow rheological processes) (e.g., Berry et al., 2005;Walker et al., 2007;Purevjav et al., 2014Purevjav et al., , 2016Caracas and Panero, 2017;Padrón-Navarta and Hermann, 2017;Tollan et al., 2017Tollan et al., , 2018Le Losq et al., 2019). All these processes related to physical properties of NAMs cannot be the consequence of a single defect. Especially, the microscopic substitution mechanisms of water in lattice and type and concentration of point defects needs to be considered carefully.
In the present study, we mainly focus on the influence of water on the physical properties of olivine, wadsleyite, and ringwoodite. Firstly, we critically review the recent advances in diffusivity, electrical conductivity, thermal properties, sound velocity, and rheology; evaluate the reliability of experimental results; and clarify the controversies and discrepancies. Secondly, we present possible links between various transport properties and discuss the geophysical applications. Finally, we provide a summary of the limitations and possibilities of the methodologies and available data and suggest some future directions.

Diffusivity
Diffusion is due to thermally activated atomic-scale random motion of particles (atoms, ions, and molecules) in minerals, glasses, melts, and gases (Zhang, 2010). The random motion leads to a net flux when the concentration of a component is not uniform. Diffusion studies in silicate minerals provide important constraints on numerous geodynamics, petrological, and geochemical processes (e.g., Watson and Baxter, 2007;Brady and Cherniak, 2010). Several researchers have made a very extensive and detailed review of diffusion coefficients of different elements (such as H, O, Fe-Mg, and other cations) in mantle minerals (e.g., Béjina et al., 2003;Ingrin and Blanchard, 2006;Brady and Cherniak, 2010;Chakraborty, 2010;Zhang, 2017). Herein, we present an overview of the latest diffusion data in olivine, wadsleyite, and ringwoodite after 2010 and mainly emphasize the effect of water on the diffusion coefficient rather than diffusion theory.
If a single mechanism of diffusion is involved, the variation of diffusivity (D) with temperature and pressure obeys the Arrhenius equation for a given species: where D 0 , R, H , E, and V are the pre-exponential factor, ideal gas constant, activation enthalpy, activation energy, and activation volume, respectively. If more influencing factors (such as water content, oxygen fugacity, and composition) are simultaneously taken into account, the diffusion coefficient, D, can be rewritten as follows (Dohmen et al., 2007;Zhang, 2017): where f O 2 is oxygen fugacity; C H 2 O is oxygen fugacity; X is composition (e.g., iron content); n, r, and p are constants; and α is a parameter related to the activation energy. It is worth noting that the diffusion in crystalline solids is strongly dependent upon point defect (vacancy, interstitial, and substitutional defects) concentrations and temperature. Various atomic mechanisms of diffusion in crystals have been identified and catalogued: vacancy mechanism, interstitial mechanism, divacancy mechanism, interstitialcy mechanism, and interstitial-substitutional exchange mechanism (e.g., Mehrer, 2007;Chakraborty, 2010;de Koker and Stixrude, 2010;Zhang, 2017). In this case, we mainly focus on the experimental data on diffusion in olivine and its highpressure polymorphs.

Olivine
Over the long history of diffusion studies in olivine, many sets of data have been accumulated. However, not all of these are consistent with each other. A summary of the new advances in diffusion in olivine is provided below.

42
B.-H. Zhang and Q.-K. Xia: Physical properties of olivine, wadsleyite, and ringwoodite significantly higher water content (150-1648 ppm) than the ones from Novella et al. (2017) (74 ppm). The smaller activation enthalpies determined by Sun et al. (2019) might be due to a lowering of activation energy of proton migration as water content increases. The anisotropy between the [010] and [001] directions reported by Sun et al. (2019) is comparable to that of Zhang (2012) but in contrast to the result from Novella et al. (2017). However, the reason for this difference is unknown, which cannot be explained by difference in the measured temperature/or pressure range for each study.
Dehydration experiment is another effective way to obtain a hydrogen diffusion coefficient in single-crystal or polycrystalline olivine. Padrón-Navarta et al. (2014) and Jollands et al. (2016) studied hydrogen diffusion in Ti-doped forsterite as a function of the type of point-defect mechanism by which the hydroxyl is incorporated. In the case of Ti-bearing forsterite, the defect chemistry is dominated by [Ti] (Ti 4+ is in an octahedral site charge-balanced by a substitution of Si 4+ by two protons on a neighboring tetrahedral site) and [Si] (four H atoms charge-balance a Si vacancy in a tetrahedral site). Demouchy and Mackwell (2006) and  investigated the incorporation and diffusion of hydrogen in San Carlos olivine single crystals by performing hydration/dehydration experiments under hydrothermal conditions at high pressure. Their results showed that hydrogen diffusion in olivine is the fastest along the [001] axis, and hydrogen incorporation is controlled by proton-metal vacancy (PV mechanism) associated defects. Ferriss et al. (2018) carried out sequential dehydration experiments on natural olivine single crystal to determine hydrogen diffusivities at 1 atm. Their experiments revealed that the apparent diffusivities of total H + vary both over time and for olivine with different defect concentrations. Consequently, hydrogen diffusion rates are slower than proton-polaron redox (in experiments at low temperatures for short times, hydrogen diffuses into olivine by exchange of protons with polarons and electron holes localized on iron atoms occupying octahedrally coordinated metal cation sites) rates but faster than proton-metal vacancy diffusion and with a fast [100] direction (Fig. 1b). Jollands et al. (2019) studied coupled inter-site reaction and hydrogen diffusion in natural olivine. They demonstrated that hydrogen loss from the Si-vacancy defect can be extremely rapid and approaches the fastest known mechanism of proton-polaron diffusion in olivine. Additionally, hydrogen diffusion through grain boundaries in olivine aggregates has been investigated by Demouchy (2010) at 0.3 GPa. She found that the effective hydrogen diffusivity, which includes the grain-boundary effect in olivine aggregates, is the fastest and at least 1 order of magnitude higher than hydrogen in an olivine single crystal. Figure 1b shows a comparison of new data on hydrogen diffusion in olivine from dehydration experiments. It is found that (1) the fastest mechanism for hydrogen diffusion is redox-associated (proton-polaron exchange) through grain boundaries in polycrystalline olivine (Demouchy, 2010), and the slowest dif-fusion is the tetrahedral site mechanism via silicon vacancy defects [Si] in forsterite single crystal (Padrón-Navarta et al., 2014). (2) Hydrogen diffusion in Ti-doped forsterite is slightly faster than that in pure forsterite. Similarly, hydrogen diffusion is faster in olivine (Ferriss et al., 2018;Jollands et al., 2019) than in forsterite (Padrón-Navarta et al., 2014;Jollands et al., 2016). (3) Hydrogen diffusion in olivine dominated by a proton-polaron mechanism (Demouchy, 2010;Jollands et al., 2019) is faster than the proton-vacancy mechanism (Demouchy and Mackwell, 2006;Padrón-Navarta et al., 2014;Jollands et al., 2016;Ferriss et al., 2018). (4) Despite the underlying cause being unclear, the anisotropy of hydrogen diffusion reported by different groups was significantly inconsistent. For example, the fastest diffusion is along the [001] orientation in Demouchy and Mackwell (2006) and Jollands et al. (2016), but the fastest diffusion is along the [100] orientation in Ferriss et al. (2018) and Jollands et al. (2019). Understanding these issues is necessary for interpreting hydrogen incorporation in olivine and designing future experiments.
Mg and Fe-Mg. Fei et al. (2018) measured Mg lattice diffusion coefficients in iron-free olivine aggregates as a function of pressure (1-13 GPa), temperature (1100-1300 K), and bulk water content (1-350 ppm) using secondary ion mass spectrometry (SIMS) in depth-profiling mode. They demonstrated that the Mg lattice diffusion coefficient was significantly enhanced by dissolved water in olivine (r = 1.2 in Eq. 2). For example, at the same pressure, the diffusion coefficient of Mg in hydrous olivine (350 ppm) is at least 2 orders of magnitude faster than that in dry olivine (Fig. 1c). The enhancement of Fe-Mg interdiffusion in olivine by the presence of water has been investigated by Hier-Majumder et al. (2005). They found a large water exponent (r = 0.9) but a slightly lower activation energy (220 kJ/mol) compared with that (280-343 kJ/mol) reported in anhydrous olivine (Tachibana et al., 2013;Zhang and Shan, 2015a).
O and Si. Costa and Chakraborty (2008) investigated the effect of water on O diffusion rate in Fe-bearing olivine. They found that, under hydrous conditions, the anisotropy of O diffusion was weak to non-existent. Fei et al. (2014) studied the influence of water content (< 1-800 ppm) on O selfdiffusion coefficients in forsterite along the b axis at 8 GPa and 1600-1800 K. Surprisingly, their experimental results indicated that water has no significant effect on O self-diffusion rate in forsterite (r ≈ 0.05). Very similar to O diffusion in olivine, Costa and Chakraborty (2008) found that water has little effect on Si diffusion rate in olivine. The diffusion of Si in forsterite single crystal along the b axis reported by Fei et al. (2012Fei et al. ( , 2013 under anhydrous and hydrous conditions indicated a small C H 2 O exponent (r = 1/3). They concluded that the effect of water on the upper-mantle rheology is very small. As shown in Fig. 1c, the O diffusivity in olivine measured by Costa and Chakraborty (2008) is consistent with that in forsterite reported by Fei et al. (2014) and Si diffusion in olivine (Costa and Chakraborty, 2008;Zhang and Shan, 2015b) but at least 1 order of magnitude higher than O diffusion calculated by Zhang et al. (2011) and Si diffusion in forsterite by Fei et al. (2012Fei et al. ( , 2013. Analysis of the point defect structure indicated that O diffused by an interstitial mechanism under hydrous conditions, whereas Si diffused by a mechanism involving vacancy complexes (Costa and Chakraborty, 2008;Fei et al., 2012Fei et al., , 2013. It is suggested that, when C H 2 O is high, the defect chemistry could be changed by incorporation of protons in Si vacancies and the hydrated Si vacancies (i.e., H Si , (2H) Si , (3H) Si , and (4H) x Si ) could dominate Si diffusion, possibly leading to a hydrolytic weakening of olivine (Costa and Chakraborty, 2008;Fei et al., 2013).

Wadsleyite
Diffusion experiments in wadsleyite are much more difficult to perform because crystals large enough to measure diffusion coefficients are difficult to obtain. As a result, diffusion coefficients that are known for these materials come from measurements in coarse polycrystals, while diffusion experiments in wadsleyite single crystals are rare. An additional difficulty associated with diffusion experiments is the in situ control of oxygen fugacity and water fugacity under high P -T conditions.
H. Hae et al. (2006) were the first to present hydrogen diffusion in polycrystalline wadsleyite at 15-16 GPa and 1173-1473 K using infrared (IR) spectroscopy. They found that the hydrogen diffusion rate in wadsleyite is essentially the same as the average diffusivity of hydrogen in olivine. The refined diffusion coefficient of hydrogen and activation energy in (Fe,Mg) 2 SiO 4 wadsleyite derived from the cB model (the cB model is a thermodynamic model , which suggests the defect Gibbs free energy g act is proportional to the isothermal bulk modulus (B) and the mean volume per atom ( ) through the relation g act = c act B , where c act is a dimensionless factor) agree well with the experimental data reported by Hae et al. (2006) (Fig. 2a). A much more extensive study of hydrogen self-diffusion in wadsleyite was performed by Sun et al. (2018) from H-D interdiffusion experiments in wadsleyite single-crystal couples. Their results demonstrated that the anisotropy of hydrogen self-diffusion in wadsleyite gradually decreases with increasing temperature.
Fe-Mg. Kubo et al. (2004) first investigated the influence of water on Fe-Mg interdiffusion in wadsleyite at 16-17 GPa and 1503-1803 K under the Ni-NiO (NNO) buffered condition. Their results showed that the Fe-Mg interdiffusion rate in wadsleyite with ∼ 330 wt. ppm H 2 O is about 1 order of magnitude faster than that containing 50-90 wt. ppm H 2 O at 1503 K, suggesting a large water effect on Fe-Mg interdiffusivity. Unfortunately, dehydrated samples in Kubo et al. (2004) during their experiments rendered the relationship between water and Fe-Mg interdiffusion in wadsleyite unreliable. By combining all existing results, Holzapfel et al. (2009) proposed a new expression to constrain the effects of pressure and composition on D Fe−Mg and found that the interdiffusivity of Fe-Mg in wadsleyite decreases with increasing pressure. Recently, Zhang et al. (2021) measured Fe-Mg interdiffusion rates in polycrystalline wadsleyite aggregates as a function of water content (up to ∼ 0.44 wt. % H 2 O) at 16 GPa and 1373-1773 K by using the diffusion couple method. Their results indicated that water significantly enhanced the rates of Fe-Mg interdiffusion in wadsleyite (r = 0.54) (Fig. 2a). In addition, the activation energy (∼ 113 kJ/mol) obtained by Zhang et al. (2021) was much lower than that (229-289 kJ/mol) determined by Holzapfel et al. (2009) and Zhang and Shan (2015a).
O and Si. Shimojuku et al. (2009) showed that O diffusion is faster than Si diffusion in wadsleyite, and the activation energy (291 kJ/mol) for O diffusion is smaller than that (513 kJ/mol) derived from a thermodynamic cB model by Zhang et al. (2011). After considering the degradation of the depth resolution by the roughness in SIMS measurement and correcting the convolution effect, the Si diffusion rates in Mg 2 SiO 4 wadsleyite under nominally dry conditions (20-60 ppm H 2 O) obtained by Shimojuku et al. (2010) were comparable to those determined in (Mg 0.9 Fe 0.1 ) 2 SiO 4 wadsleyite with similar water contents (10-82 ppm H 2 O) by Shimojuku et al. (2009) but approximately half an order of magnitude smaller than those reported by Shimojuku et al. (2004) Fig. 2a). (In microbeam analysis, the measured concentration profile is always continuous (even the element is discontinuous at the interface of two samples) owing to the spreading of the ion beam and spatial averaging of compositions on both sides of the interface. It is a convolution result of both the finite abrasion of the sample volume and the interface of two samples.) This finding indicates that the effect of water is significantly larger than the convolution effect. Enhancement of the Si diffusivity by incorporating hydrogen possibly leads to water weakening in wadsleyite.

Ringwoodite
H. The diffusion of hydrogen in synthetic polycrystalline ringwoodite with a grain size of approximately 10 µm was measured by Kudo et al. (2006) in the temperature range from 1173 to 1373 K at 19 GPa. The diffusion rate of hydrogen in ringwoodite is expressed by using the following equation: D = 1.6 × 10 −4 exp(−140 (kJ/mol)/RT ). Given that the diffusion coefficient in polycrystalline samples is contributed by grain-boundary diffusion and the grain size in real mantle may be greater (∼ millimeter or more), the experimental results of Kudo et al. (2006) are likely to be the maximum value of the diffusion coefficient of hydrogen in the real MTZ. As shown in Fig. 2b for hydrogen diffusion in (Mg,Fe) 2 SiO 4 ringwoodite, the experimental data reported by Kudo et al. (2006) are comparable with that calculated by the cB model (Zhang, 2012). Sun et al. (2015) reported a hydrogen self-diffusion coefficient in Fe-bearing ringwood-  ; Z12 (Zhang, 2012); S18 (Sun et al., 2018); K06 (Kudo et al., 2006); K04 ; H09 (Holzapfel et al., 2009); S04 ; S09 (Shimojuku et al., 2009); S10 (Shimojuku et al., 2010); Z11 (Zhang et al., 2011); Z15a (Zhang and Shan, 2015a); Z15b (Zhang and Shan, 2015b); Z19b (B. H. Zhang et al., 2019b); Z21  ite through the interdiffusion of hydrogen and deuterium in a pair of synthesized single crystals with 5500-7700 ppm water at 21 GPa and 1000-1300 K. Note that the results obtained by Sun et al. (2015) are 2 orders of magnitude lower than those reported by Kudo et al. (2006) and Zhang (2012), despite the activation energies being comparable to each other (Fig. 2b).
Fe-Mg. To date, only B. H. Zhang et al. (2019b) have determined the kinetics of Fe-Mg interdiffusion in ringwoodite aggregates as a function of water content (40-6000 ppm H 2 O) at 20 GPa and 1373-1673 K by using the diffusion couple method. A relatively small water exponent (r = 0.25) suggests the non-negligible role of water in enhancing Fe-Mg interdiffusion in ringwoodite. The activation energy (140 kJ/mol) determined by B. H. Zhang et al. (2019b) was much lower than that (345 kJ/mol) reproduced in terms of bulk elastic and expansivity data through a thermodynamic model that interconnects point defect parameters with bulk properties (Zhang and Shan, 2015a).
O and Si. To the best of our knowledge, only Shimojuku et al. (2009) have reported O diffusion rates in ringwoodite to date. The experiments were conducted for (Mg,Fe) 2 SiO 4 ringwoodite at 22 GPa and 1673-1873 K using SIMS measurement. With similar water contents (130-220 ppm H 2 O), the diffusion rates of O in ringwoodite were slightly higher than those of Si. Meanwhile, the diffusivity and activation energy for Si diffusion in ringwoodite are consistent with those derived from the cB model on the basis of the observed compensation effect for Si diffusion in silicates under anhydrous conditions (Zhang and Shan, 2015b). However, the activation energy for O diffusion in ringwoodite reported by Shimojuku et al. (2009) is inconsistent with that determined by Zhang et al. (2011).

Effect of water on diffusivity
On the basis of the foregoing comparisons (Figs. 1 and 2) and discussions, the diffusion rates of different elements in olivine and its high-pressure polymorphs, wadsleyite and ringwoodite, yield the order of D H > D Fe−Mg > D O (≈ D Si ). Their water content exponents are in the order of r Fo Mg > r Ol Fe−Mg > r Wad Fe−Mg > r Ring Fe−Mg (Fig. 3). This finding implies that water can significantly enhance the D Fe−Mg in olivine compared with wadsleyite and ringwoodite, which needs to be verified by further experiments. The magnitude of the enhancement of interdiffusivity in wadsleyite by water is larger than that in ringwoodite (D Wad Fe−Mg > D Ring Fe−Mg ) when C H 2 O > 0.2 wt. %. These different characteristics of diffusion response to water may be related to the hydrogen incorporation mechanism in host minerals. Experimental observations (Purevjav et al., 2016;Zhang et al., 2021) and theoretical calculations (Wright and Catlow, 1996;Caracas and Panero, 2017) demonstrated that, for Fe-free hydrous wadsleyite, all hydrogen cations in the crystal should have been fully concentrated into a vacant M3 octahedral sites bonded with oxygen O1. In the case of Fe-bearing hydrous wadsleyite, reduction from ferric to ferrous iron probably drives the formation of Mg vacancies in wadsleyite, which in turn contributes to Fe-Mg interdiffusion. In the case of ringwoodite, the mobility of hydroxyl defects related to an octahedral vacancy is greater than that related to a tetrahedral vacancy (Mrosko et al., 2013). Furthermore, although both tetrahedral and octahedral sites can be simultaneously hydrated in hydrous ringwoodite (Purevjav et al., 2014), Caracas and Panero (2017) The smaller water content exponent obtained here (r ≈ 0.25 < 1) implies that not all hydrogen is involved in the diffusion process.

Electrical conductivity
Electrical conductivity is very sensitive to temperature and water; therefore, electrical conductivity can be used to constrain and quantify the chemical composition, temperature, and water content and distribution in the Earth's interior (Xu et al., 2000;Yoshino, 2010). In most cases, the T and P dependences of electrical conductivity can be described by the Arrhenius equation: where σ 0 is the pre-exponential factor. The electrical conductivity is generally dominated by one or more conduction mechanisms. For example, the electrical conductivity of hydrous iron-bearing silicate minerals can be expressed by the summation of "dry" (sum of ionic σ i and hopping conduction σ h ) and "wet" (proton conduction σ p ) mechanisms Yoshino, 2010;Yoshino and Katsura, 2013): Since the first suggestion by Karato (1990) of the enhancement of electrical conductivity by water in olivine, numer-ous experimental studies have been carried out on the role of water in electrical conductivity in olivine (e.g., Yoshino et al., 2006Poe et al., 2010;Yang, 2012; and its high-pressure polymorphs (e.g., Huang et al., 2005;Yoshino et al., 2008;Romano et al., 2009;Manthilake et al., 2009;Dai and Karato, 2009;Yoshino and Katsura, 2012). Because several review papers have provided a very detailed and comprehensive overview of this topic (Yoshino, 2010;Yoshino and Katsura, 2013;Pommier, 2013;Karato, 2015Karato, , 2019, here we will only give a brief introduction and comparison to the influence of water on the electrical conductivity of olivine, wadsleyite, and ringwoodite.

Olivine
Single crystal. The electrical conductivity anisotropy of hydrous olivine is an important issue for understanding the anisotropic σ beneath the East Pacific Rise (Evans et al., 2005). The effect of water on the electrical conductivity of single-crystal olivine has been reported by several groups (Fig. 4). Yoshino et al. (2006) first measured the electrical conductivity of hydrous single-crystal olivine with 0.01-0.02 wt. % H 2 O at 3 GPa under low-temperature conditions (< 1000 K) to minimize possible dehydration of olivine. Nevertheless, they did not determine how conductivity varies with varying water along any of the three crystallographic directions, but rather assumed that conductivity increased linearly with water content without change in activation energy. The extrapolation of their experimental data to higher temperatures suggests that the electrical conductivity of hydrous olivine at the top of the asthenosphere should be nearly isotropic. Therefore, Yoshino et al. (2006) concluded that the hydration of olivine cannot account for the high-conductivity anisotropy in the UM (Evans et al., 2005). Poe et al. (2010) investigated the electrical conductivity of single-crystal San Carlos olivine as a function of dissolved water content (up to 0.22 wt. % H 2 O) at 8 GPa. They found that σ p is largely anisotropic, the conductivity along [010] is approximately 1 order of magnitude greater than that along other directions when the water content is larger than 0.1 wt. % at high temperature, whereas the electrical anisotropy is negligible when the water content is lower than 0.1 wt. %. Similarly, Yang (2012) reported negligible anisotropy for the proton conduction of a single-crystal olivine with low water content (∼ 40 ppm) at 1 GPa. By contrast,  presented new experimental results on the electrical conductivity of hydrous single-crystal olivine at temperatures up to 1373 K and pressure of 4 GPa. They argued that their new results agree well with previous studies at low temperature (< 900 K) and significantly deviate at high temperature where the temperature dependence of conductivity increases strongly (Fig. 4). An important point of the investigations of  is that the electrical conductivity anisotropy ( temperature, and a change in conduction mechanism is observed because the activation energy substantially increases at high temperature. Recently, Fei et al. (2020) measured the ionic conductivity of olivine single crystal as a function of pressure from 2 to 10 GPa, temperature from 1450 to 2180 K, and water content from 20 to 580 wt. ppm. They revealed that olivine ionic conductivity has negative pressure and positive temperature dependences and is significantly enhanced by H 2 O incorporation. It is clear from Fig. 4 that the anisotropy of proton conduction of olivine reported by different groups disagree with each other. The experimental results from Yang (2012) are consistent with those measured by Yoshino et al. (2006) and Poe et al. (2010). Note that Poe et al. (2010) demonstrated that the conductivity along [010] is the highest in olivine single crystal, whereas the highest conductivity is along the [100] direction in other studies (Yoshino et al., 2006;.  found that the magnitude of electrical conductivity anisotropy increases with an increase in temperature, which is the opposite of early experimental observations (Yoshino et al., 2006;Yang, 2012). In addition, the anisotropic ionic conductivity and activation enthalpies (337-396 kJ/mol) reported by Fei et al. (2020) are significantly higher than those of proton conductivity (Yoshino et al., 2006;Poe et al., 2010;Yang, 2012;. Polycrystalline. D.  measured the electrical conductivity of hydrous olivine aggregates at temperature up to 1273 K, which showed high conductivity and strong temperature dependence, and proposed the following relation to fit their experimental data: where C H 2 O is the water content (in wt. %) and r is the water exponent. However, their results may have been biased because of the significant dehydration of olivine at high temperature, and the contribution of hopping conduction (σ h ) to bulk conductivity was not considered.  identified three conduction mechanisms by measuring the electrical conductivity of olivine aggregates with various amounts of water at a pressure of 10 GPa. They demonstrated the significant effect of water on the electrical conductivity of hydrous olivine and suggested that the behavior of proton conductivity with water content is similar to that of N-type semiconductors through the equation where α is a constant. As shown in Fig. 5a, the electrical conductivity of anhydrous olivine aggregates is consistent with others (Xu et al., 2000;Constable, 2006;. However, at a similar water content, the electrical conductivity determined by D. Wang et al. (2006) is at least 2 orders of magnitude higher than that of . Figure 5b shows a compilation of the influence of water content on the bulk conductivity of olivine aggregates at 1573 K. It is found that D.  reported the highest conductivities, while Du Frane and Tyburczy (2012) had the lowest ones that were calculated from H-D interdiffusion in a single crystal. When the water content is low (< 0.01 wt. %), the difference in conductivity data among Yoshino's group (Yoshino et al., 2006), Poe's group (Poe et al., 2010), and Frane and Tyburczy (2012) is small (within 0.2 log unit), while the results vary widely at higher water content. The possible reasons for these discrepancies will be discussed later.

Wadsleyite
Xu et al. (1998) reported that the electrical conductivity of wadsleyite under MTZ conditions is 2 orders of magnitude higher than that of olivine. Although the results of Xu et al. (1998) were widely accepted as the standard data for hopping conduction in wadsleyite, Huang et al. (2005) attributed the results of Xu et al. (1998) to proton conduction because they found a significant amount of water in Xu et al.'s samples. On the other hand, Huang et al. (2005) and Romano et al. (2009) only reported the proton conduction, whereas other studies (Yoshino et al., 2008;Dai and Karato, 2009;Manthilake et al., 2009;Yoshino and Katsura, 2012) distinguished between small polaron and proton conduction mechanisms of wadsleyite in a wider temperature range. Figure 6a shows a compilation of electrical conductivity data for anhydrous and hydrous wadsleyite at various water contents. The published conductivity data regarding the effect of water on σ p in wadsleyite disagree among different groups (Fig. 6), which can be roughly divided into three categories: Karato's group (Huang et al., 2005;Dai and Karato, 2009), Yoshino's group (Yoshino et al., 2008;Manthilake et al., 2009;Yoshino and Katsura, 2012;Sun et al., 2018), and Romano's dataset (Romano et al., 2009). Each set of experimental data is internally consistent in terms of the effect of water on conductivity (σ and r) and activation enthalpy (H ). For anhydrous wadsleyite, Yoshino's group (Yoshino et al., 2008;Yoshino and Katsura, 2012) reported higher conductivity values (0.5 log unit) than Karato's group ), although these two groups found similar values for H (∼ 145 kJ/mol) of hopping conduction. As for the hydrous case, the conductivity values measured by Karato's group (Huang et al., 2005;Dai and Karato, 2009) are much higher than those of Yoshino's group (Yoshino et al., 2008;Yoshino and Katsura, 2012) at the same water content ( Fig. 3a and b). Romano et al. (2009) obtained the highest absolute conductivity values for hydrous wadsleyite among the published data ( Fig. 6a), but they reported the lowest H value (∼ 64 kJ/mol). Note that the trend of bulk conductivity (σ i + σ h + σ p ) with water content at 1773 K is largely different among Karato's group (Huang et al., 2005;Dai and Karato, 2009), Yoshino's group (Yoshino et al., 2008;Yoshino and Katsura, 2012;Sun et al., 2018), and Romano et al. (2009) (Fig. 6b) because each group obtained different values for H and r and used different fitting equations. Consequently, although all of these studies suggest that water can significantly enhance the electrical conductivity, the measured absolute conductivity values and activation enthalpies vary widely between studies.

Ringwoodite
Huang et al. (2005) first measured the electrical conductivity of hydrous ringwoodite, but they did not determine the hopping conduction at high temperatures. Similar to wadsleyite, Yoshino et al. (2008) successfully separated the contributions of hopping and proton conductions in ringwoodite (Fig. 7a). Several important features were revealed by their experiments: (1)   approximately 1 order of magnitude higher than that of wadsleyite; (2) the activation energy of proton conduction largely decreases from 95 to 43 kJ/mol with increasing water content from 0.01 wt. % to 1 wt. %; (3) the electrical conductivity increases with increasing water content. For example, the contribution of proton conduction to the bulk conductivity (σ i + σ h + σ p ) and calculated conductivity from H-D interdiffusion in single-crystal ringwoodite is negligible at 1873 K and water content below 0.1 wt. % (Sun et al., 2015), but it is much larger than that of Huang et al. (2005) at water content above 0.5 wt. % (Fig. 7b). (4) Small polaron conduction becomes dominant at high temperatures. Overall, the dependence of conductivity on water content measured by Huang et al. (2005) is much weaker than that reported by Yoshino et al. (2008). It is worth noting that, when the water content increases from 0.1 wt. % to 0.5 wt. %, the electrical conductivity increases by only 0.5 log units in Huang et al. (2005); in contrast, it goes up by 0.6-2.5 log units in Yoshino et al. (2008). Indeed, this discrepancy may reflect the different water dependence of conductivity. The water exponent is estimated to be about 0.69 by Huang et al. (2005), while it is equal to 1 in the study of Yoshino et al. (2008).

Discrepancies and interpretations
On the basis of the above discussions, available experimental data on the magnitude of water's effect on olivine, wadsleyite, and ringwoodite are largely inconsistent as illustrated in Figs. 4-7. These discrepancies between different studies have stimulated several controversies and ongoing debates (e.g., Karato and Dai, 2009;, 2013Jones, 2014Jones, , 2016Gardés et al., 2015;Dai and Karato, 2015). Importantly, the inconsistencies in the experimental data reported by different researchers have led to large differences in the estimated water content of the UM and MTZ and also created great confusion for the geoscience community. Meanwhile, these controversies and disagreements have the unfortunate and unsatisfactory consequence that those who wish to make use of the experimental results of the effects of water on NAMs must make a decision as to which laboratory's values to choose over the others because they do not know whose experimental data are more reliable. The possible causes of these existing discrepancies can be simply attributed to three aspects. (1) Experimental technique (low-frequency versus impedance spectroscopy measurements with or without a metal shield in the assembly). Karato and Dai (2009) criticized the use of the low-frequency method by Yoshino's group (Yoshino et al., 2006(Yoshino et al., , 2008 as it tends to produce conductivity tens of a percent smaller than that obtained by impedance spectroscopy. In fact, later studies on hydrous ringwoodite (Yoshino et al., 2008), olivine , and wadsleyite  have proven that the difference in conductivity measured by low-frequency and impedance spectroscopy is very small and that the large data inconsistency between different groups cannot be explained by this technical reason. Another technical difference is that use of the metal guard ring by Karato and colleagues (e.g., Huang et al., 2005;Poe et al., 2010) in their assembly would prevent water purge and lead to the formation of fluid or supercritical liquid by the released free water due to the dehydration of water-rich samples at high temperatures (> 900 K) as suggested by Yoshino (2010), Yoshino andKatsura (2012, 2013), and Yang (2012).
(2) Water content determination (unpolarized Fouriertransformation infrared (FTIR) spectroscopy using Paterson calibration versus polarized FTIR by Bell calibration for determining water loss/gain). The determination of sample water content before and after conductivity measurement is a key procedure to judge whether the sample is changed (water loss or gain). Most previous studies between Yoshino's group (e.g., Yoshino et al., 2006Yoshino et al., , 2008Yoshino and Katsura, 2012) and Karato's group (e.g., Huang et al., 2005;Dai and Karato, 2009) used the unpolarized FTIR spectroscopy and the calibration from Paterson (1982) to determine the water content. However, Bell et al. (2003) claimed that the unpolarized FTIR method with the calibration of Paterson (1982) may underestimate the water content in NAMs (especially olivine) by a factor of 3 compared with the application of the Bell et al. (2003) calibration. If this is the case, the conductivity values of NAMs obtained using the unpolarized calibration of Paterson (1982) may also be underestimated by the same factor, except for the study of Poe et al. (2010).
(3) Sample dehydration during conductivity measurement. As for the dehydration of hydrous samples, another very confusing point is that Yoshino's group (Yoshino et al., 2006(Yoshino et al., , 2008Yoshino and Katsura, 2012;B. H. Zhang et al., , 2019cZhao and Yoshino, 2016) and others (Poe et al., 2010;Yang, 2012) found that water-bearing NAM samples (e.g., olivine, wadsleyite, ringwoodite, orthopyroxene, omphacite, and clinopyroxene) were easily dehydrated during annealing at temperatures higher than 900 K, whereas no sample dehydration occurred in the study of  even at 1373 K. In fact, the dehydration temperature of a hydrous sample should not vary from laboratory to laboratory.
(4) Formalism to fitting experimental data (Eq. 5 versus Eq. 6). Considering the different dependences of activation enthalpy on water content, Karato's group used Eq. (5) to fit their conductivity data, whereas Yoshino's group used Eq. (6). Obviously, different fitting equations will lead to systematic differences in predicted conductivity data between different laboratories. Consequently, the laboratorybased laws of NAMs' electrical conductivity predict contrasting effects of water, precluding the interpretation of geophysical data in terms of mantle hydration. We found that several researchers Jones, 2014;Karato, 2019) have attempted to develop a unified law to reconcile the large discrepancies in available data on olivine conductivity and test the internal consistency of the database. They argued that the experimental data from different laboratory studies are mostly consistent when the uncertainties and biases in the water contents of the olivine samples are considered. In addition, Karato and Dai (2009) criticized the result from small polaron conduction in wadsleyite and ringwoodite obtained by Yoshino et al. (2008) as not being small polaron conductivity because the samples are not truly dry ones. In fact, preparing truly dry NAM samples under high pressure is very difficult because samples prepared without adding water often contain a certain amount of water (e.g., Yoshino et al., 2008;B. H. Zhang et al., 2019c). From a purely mathematical point of view, the contributions of hopping and proton conductions can be completely separated and predicted by the use of Eq. (4) associated with Eq. (5) or (6). Therefore, the conductivity data from a "truly dry" sample are unnecessary as long as the data of a hydrous sample are reliable .
To sum up, one solution to reconcile or eliminate the abovementioned discrepancies between different laboratories is to use the same water content determination methods and the same formula to fit the experimental data, which will surely confirm which laboratory's experimental data are more reliable. In addition, it is more important to understand the fundamental physical reasons explaining why different experimental methods lead to such differences. We look forward to the early realization of this goal.
where ρ is density and C P is heat capacity at constant pressure. Heat is transmitted within the Earth's crust and mantle mainly by diffusion (in the steady state: conduction) and radiation. As a result, the total thermal conductivity consists of contributions from scattering of quantized lattice vibrations (phonons, termed κ p ) and radiative attenuation of photons (κ r ). Generally, in most polycrystalline materials, heat diffusion dominates heat radiation at low temperatures, whereas heat radiation starts to dominate at higher temperatures (> 1500 K). At present, direct measurement of radiative heat transfer is still very difficult because of the challenge of experimental technology. Thus, only phonon conduction is presented and discussed here.
In the past decades, substantial collections of mineral thermal properties were compiled by many researchers (e.g., Clauser and Huenges, 1995;Hofmeister, 2007;Clauser, 2011;Hofmeister et al., 2007;Xiong and Zhang, 2019). In this review, we do not present a complete compilation for all available thermal transport property data published to date. Instead, the latest experimental data on the effect of water on the thermophysical properties of olivine, wadsleyite, and ringwoodite are presented.

Olivine
The most prominent methods for measuring thermal diffusivity and thermal conductivity of geomaterials are steadystate (e.g., Ångström method) and transient techniques (e.g., pulse heating method and laser-flash analysis) (Clauser and Huenges, 1995;Hofmeister et al., 2007). The λ and κ of anhydrous olivine have been reported by several researchers (Katsura, 1995;Xu et al., 2004;Osako et al., 2004). Recently, Y. Zhang et al. (2019) applied a pulse heating method to explore the influence of iron on the λ and κ of olivine aggregates. Their results showed that the κ of olivine first decreases and then slightly increases with increasing iron content in olivine. Ge et al. (2021) measured the λ and κ of dunite up to 3 GPa and 823 K. They showed that λ and κ have a very weak temperature dependence ( Fig. 8a and b), which is likely due to the presence of some amounts of serpentine in the dunite sample because the λ and κ of serpentine are low and almost independent of temperature (Osako et al., 2010;Ge et al., 2021).
Compared with the available data on anhydrous olivine, only two studies have reported on the effect of water on the thermal properties of olivine so far. Chang et al. (2017) first determined the κ of hydrous single-crystal Fo 90 at room temperature and pressures up to 15 GPa by using an ultrafast optical pump-probe technique. They found that the κ of hydrous single-crystal Fo 90 has a weaker pressure depen-dence below 5 GPa. Surprisingly, the values of κ for hydrous olivine are even larger than those of anhydrous olivine under similar P -T conditions (Fig. 8b), which may be due to the large uncertainty of the experimental method itself. By contrast, B. H. Zhang et al. (2019a) systematically investigated the λ and κ of polycrystalline Fo 90 as a function of pressure, temperature, and water content (up to 0.2 wt. %). Their experimental results demonstrated that water can significantly reduce the λ and κ of olivine aggregate ( Fig. 8a  and b). With the increase in water content from 0.08 wt. % to 0.2 wt. %, the absolute values of the λ and κ for olivine samples decrease by 5 %-13 % and 17 %-33 % and by 3 %-8 % and 14 %-21 % (Fig. 8c and d), respectively.
For a given pressure, the overall trends of decreasing the λ and κ with temperature for olivine measured using different experimental methods are roughly consistent as shown in Fig. 8a and b. However, we must emphasize two issues. First, the effects of temperature and pressure on λ and κ were considered separately in most previous studies (e.g., Katsura, 1995;Osako et al., 2004). For example, λ and κ decrease with temperature and are fitted by an empirical form, a + b/T + c/T 2 ; at constant temperature, λ and κ increase with increasing pressure and can be linearly fitted by a 0 + d × P . Such separated fitting equations further hinder the extrapolation of the available experimental results to the real mantle P -T conditions. Fortunately, Xu et al. (2004) proposed a unified power-law formula to describe the combined influence of temperature and pressure on the λ and κ of olivine: λ = λ 300 K,1 bar 300 T n (1 + a × P ) , where λ 300 K, 1 bar and κ 300 K, 1 bar are the fitting results of λ and κ at room pressure and 300 K, a and b are pressure factors, and n is a fitting parameter. Eqs. (8)-(9) can better fit the λ and κ results obtained from experimental measurements (Xu et al., 2004;Y. Zhang et al., 2019;B. H. Zhang et al., 2019a) and can also conveniently scale laboratory results to real mantle conditions. Second, besides the effect of pressure and temperature, other factors (Tommasi et al., 2001;Y. Zhang et al., 2019;Ge et al., 2021) such as composition, water, crystal structure, and grain size can also affect the λ and κ of minerals including olivine. For example, even at the same pressure (3 GPa) and using a similar impulse heating technique, the thermal conductivity of anhydrous olivine reported by B. H. Zhang et al. (2019a) is higher than that determined by Y.   (Katsura, 1995;Xu et al., 2004;Osako et al., 2004;Y. Zhang et al., 2019;Ge et al., 2021) and hydrous olivine (Chang et al., 2017;B. H. Zhang et al., 2019a) (Fig. 8a and b). This difference may originate from the presence of a small amount of water in Xu et al.'s samples because protonation increases structure disorder of NAMs and adds new vibration modes. As a result, the number of scattering events is more significant in the wet phase, which leads to reduced mean free paths of phonons and decreased λ and κ (B. H. Zhang et al., 2019a).

Wadsleyite and ringwoodite
Considering the substantial difficulties in performing thermal property measurements under P -T conditions of the MTZ, Xu et al. (2004) first reported the lattice λ and κ of anhydrous (Mg 0.9 Fe 0.1 ) 2 SiO 4 wadsleyite and ringwoodite at 14 and 20 GPa by using the Ångström method on cylindri-cal polycrystalline specimens. With increasing pressure, the olivine phase transforms to its high-pressure polymorphs, wadsleyite and ringwoodite, and the corresponding λ and κ increase in the order of (λ or κ) olivine < (λ or κ) wadsleyite < (λ or κ) ringwoodite (Fig. 9), which may suggest a discontinuity of κ occurring at 410 and 520 km. Recently, Marzotto et al. (2020) investigated the effect of water on the lattice thermal conductivity of ringwoodite (0.11 wt. %-1.73 wt. % H 2 O) at pressures up to 25 GPa and room temperature by combing DAC experiments with ultrafast optics. Their results indicated that the incorporation of 1.73 wt. % water substantially reduces the ringwoodite thermal conductivity by more than 40 % at the MTZ pressures. However, Marzotto et al. (2020) did not measure the temperature dependence on the thermal conductivity of hydrous ringwoodite because their experiments were carried out under high pressure and 300 K. Given the importance of wadsleyite and ringwoodite in the MTZ, the effect of water on the λ and κ of wadsleyite and ringwoodite needs to be systematically investigated in the future.
Recently, simultaneous in situ Brillouin spectroscopy and synchrotron single-crystal X-ray diffraction measurements of the sound velocities and density of single-crystal minerals in an externally or internally heated diamond anvil cell (DAC) were extended to P -T conditions of the MTZ (e.g., Mao et al., 2012Mao et al., , 2015Zhang and Bass, 2016). The successful application of these techniques to the study of the behavior and physical properties of Earth materials, including unquenchable mantle phases, liquids, and multi-phase aggregates, will greatly improve our interpretation of seismic observations and our understanding of the structure, composition, and dynamics of the Earth's interior.

Olivine
So far, only Mao et al. (2010) have reported measurements of the single-crystal elasticity and velocity of forsterite with 0.9 wt. % H 2 O to 14 GPa at 300 K by Brillouin scattering. Their results demonstrated that hydration could decrease the magnitude of the velocity contrast at 410 km depth and the V P /V S ratio of forsterite. Figure 10a and b show the V P and V S of forsterite with 0.9 wt. % H 2 O compared with previous data obtained under dry conditions at 300 K (Zha et al., 1996;Li et al., 2004;Darling et al., 2004;Liu et al., 2005;Zhang and Bass, 2016), respectively. It is found that the V P of hydrous forsterite (Mao et al., 2010) is almost identical to that of anhydrous forsterite reported by Li et al. (2004) but slightly higher (∼ 0.8 %) than that measured by Zha et al. (1996). The same trend is also observed in the V S (∼ 1 %) of hydrous forsterite (Fig. 10b), which suggests that water can even enhance the seismic velocity. Compared to that of the H 2 O effect, Fe can significantly reduce the V P (∼ 1.6 %) and V S (∼ 2.6 %) of anhydrous Fe-bearing olivine (Darling et al., 2004;Liu et al., 2005;Zhang and Bass, 2016) at 1 bar pressure. Note that the V P and V S of anhydrous Febearing olivine measured by Zhang and Bass (2016) showed a marked decrease at pressure more than 7 GPa ( Fig. 10a and  b). Nevertheless, the low gradients of V P and V S versus pressure on anhydrous olivine reported by Mao et al. (2015) and Zhang and Bass (2016) suggest that linear extrapolations to high pressure/temperature based on experiments performed over a limited pressure/temperature range may overestimate the velocities of olivine at 410 km depth (Zha et al., 1996;Li et al., 1996Li et al., , 2004Darling et al., 2004;Liu et al., 2005;Mao et al., 2010).
The V P /V S ratio is commonly used to identify compositional variations in the Earth's mantle. Figure 10c shows the effect of hydration on the V P /V S ratio of Fe-free and Fe-bearing olivine. At ambient pressure, the V P /V S ratio (1.7111) of anhydrous forsterite is almost the same (∼ 1.7058) as that of the hydrous one. With increasing pressure, the V P /V S ratio gradually increases. At pressure corresponding to 410 km depth, the V P /V S of hydrous forsterite with 0.9 wt. % H 2 O increases to 1.7725 compared with 1.7889 for the anhydrous phase, a difference of 0.92 %. At the same time, the V P /V S ratios of anhydrous Fe-bearing olivine (Darling et al., 2004;Liu et al., 2005;Zhang and Bass, 2016) are at least 0.79 % higher than those of anhydrous forsterite. This finding suggests that increasing iron content will increase the V P /V S of forsterite, whereas hydration will decrease the V P /V S . Similarly, the effects of hydration and iron on Poisson's ratio in Fig. 10d are almost the same as those of the V P /V S of hydrous and anhydrous olivine.

Wadsleyite
In addition to Fe, hydration is known to decrease the velocities and elastic properties of olivine polymorphs as a result of bond-weakening defects that accompany hydrogen incorporation as a structural hydroxyl (e.g., Inoue et al., 1995;Kohlstedt et al., 1996;Jacobsen et al., 2004;Jacobsen and Smyth, 2006  With the available experimental data on V P and V S , the effect of iron and hydration on the sound velocities of wadsleyite can be evaluated quantitatively (Fig. 11). Under ambient conditions, compared with the V P (Fig. 11a) and V S (Fig. 11b) of anhydrous pure-Mg wadsleyite reported by Li et al. (1996), the presence of 8 mol %-13 mol % Fe in wadsleyite lowers the V P and V S by 1.4 % and 2.3 %, respectively (Li and Liebermann, 2000;Liu et al., 2009;Wang et al., 2014). When adding 0.24 wt. % and 1.93 wt. % H 2 O in Fe-bearing wadsleyite (Mao et al., 2011;Buchen et al., 2018), the V P and V S were further lowered by water, leading to total reductions of 3.3 % and 4.7 % in V P and 3.7 % and 5.5 % in V S , respectively. Similarly, adding 1.63 wt. % H 2 O in Fe-free wadsleyite lowers the V P and V S by 3.7 % and 4.2 %, respectively , which are almost equivalent to those of Fe-bearing wadsleyite with 0.24 wt. % H 2 O (Buchen et al., 2018). In combining the results on anhydrous Fe-free and Fe-bearing wadsleyite (Li et al., 1996;Li and Liebermann, 2000;Liu et al., 2009;Wang et al., 2014) with those of slightly hydrous and largely hydrous wadsleyite (Mao et al., 2011;Buchen et al., 2018;Wang et al., 2019), iron and water incorporation lowers the sound velocities of wadsleyite. However, the two effects are difficult to distinguish using sound velocities alone. Figure 11c and d show the V P /V S ratio and Poisson's ratio of anhydrous Fe-bearing wadsleyite (Li and Liebermann, 2000;Liu et al., 2009;Wang et al., 2014), which are indistinguishable from those of hydrous Fe-bearing wadsleyite (Mao et al., 2011;Buchen et al., 2018) and hydrous pure-Mg wadsleyite . Among them, the calculated values for the V P /V S ratio and Poisson's ratio in anhydrous pure-Mg wadsleyite (Li et al., 1996) are the lowest. This observation indicates that an increase in the iron and water contents in wadsleyite will increase the V P /V S ratio and Poisson's ratio, which is contrary to hydration decreasing the V P /V S ratio and Poisson's ratio in Fe-free olivine ( Fig. 10c and d).

Ringwoodite
J.  measured the sound velocities and single-crystal elastic moduli of hydrous Mg 2 SiO 4 ringwoodite containing 2.3 wt. % H 2 O using Brillouin spectroscopy at pressure up to 23.4 GPa and at room temperature. They speculated that the MTZ seismic velocity gradients are not due to "wet" ringwoodite. Jacobsen et al. (2004) and Jacobsen and Smyth (2006) reported the effects of 1 wt. % H 2 O on the sound velocities and elastic moduli of singlecrystal ringwoodite with Fo 90 composition up to 9 GPa. Both experiments were performed using a gigahertz ultrasonic interferometer to monitor P-and S-wave travel times under hydrostatic compression in a DAC. Their results indicated that water can significantly reduce the sound velocities and shear moduli of hydrous Fe-bearing ringwoodite; thus, the V P /V S ratios were elevated by 2.3 % at pressure less than 18 GPa. Mao et al. (2012) investigated the hydration effect on the single-crystal elastic moduli of Febearing ringwoodite with 1.1 wt. % H 2 O at pressure up to 12 GPa and temperature up to 673 K. They found that the presence of 1.1 wt. % H 2 O lowers the elastic moduli of Febearing ringwoodite by 5 %-9 %, whereas the sound velocities of hydrous Fe-bearing ringwoodite were modeled from a thermal Birch-Murnaghan equation of state but not directly measured. Their results suggest that the observed seismic velocity anomalies and related depth depression of the 660 km discontinuity could be attributed to thermal variations together with the presence of ∼ 0.1 wt. % H 2 O. Recently, Schulze et al. (2018) measured the V P and V S of hydrous (Mg 0.89 Fe 0.11 ) 2 SiO 4 (with 0.21 wt. %, 1.04 wt. % and 1.71 wt. % H 2 O) and Mg 2 SiO 4 (with 0.42 wt. % H 2 O) ringwoodite at pressure up to 22.1 GPa and at room temperature. They showed that the effect of water on the V P and V S decreases with increasing pressure, and the hydration-induced reduction of seismic velocities almost vanishes under conditions of the MTZ.
To evaluate the effect of water on the sound velocities of ringwoodite at high pressure, Fig. 12a and b compare the V P and V S of anhydrous and hydrous ringwoodite (Li, 2003;Sinogeikin et al., 2003;Jacobsen and Smyth, 2006;J. Wang et al., 2006;Higo et al., 2008;Schulze et al., 2018) at MTZ pressures and 300 K. If the V P and V S of pure-Mg anhydrous ringwoodite reported by Li (2003) are considered a benchmark, the presence of 9 mol % Fe in ringwoodite (Sinogeikin et al., 2003;Higo et al., 2008) could reduce the V P and V S by 2.2 % and 2.54 % under ambient conditions, respectively. As shown in Fig. 12a, the magnitude of the effect of Fe content on the V P of anhydrous ringwoodite (Sinogeikin et al., 2003;Higo et al., 2008) is comparable to that of the slightly hydrous ringwoodite (Schulze et al., 2018; 0.21 wt. % H 2 O) but lower than those of the large hydration effect (J.  and a combination of the Fe and H 2 O effect (Jacobsen and Smyth, 2006; Schulze et al., 2018). Unlike V P , the influence of Fe and H 2 O on the V S of ringwoodite is great in Fig. 12b. Under ambient conditions, the V S of anhydrous ringwoodite with 9 mol % Fe (Sinogeikin et al., 2003;Higo et al., 2008) or Fe-free ringwoodite with 2.3 wt. % H 2 O (J.  is about 2.54 % or 4.09 % lower, respectively, than that of the pure-Mg anhydrous ringwood-ite reported by Li (2003). However, the V S of hydrous Febearing ringwoodite with 1.0 wt. % H 2 O reported by Jacobsen and Smyth (2006) is 9.02 % lower than Li's results. Notably, the V P and V S from J.  and Jacobsen and Smyth (2006) indicate obvious crossover with other studies because of the relatively high gradients of V P and V S versus pressure. Schulze et al. (2018) demonstrated that the effect of water on the V P and V S is large at low pressure (< 12 GPa), while this effect decreases with increasing pressure and is even negligible under the conditions of the MTZ. They concluded that water might not be detectable by seismic velocity variations in the MTZ. Figure 12c shows the calculated V P /V S ratio of anhydrous and hydrous ringwoodite as a function of pressure at room temperature. For Fe-free ringwoodite, water can decrease the V P /V S ratio of pure-Mg ringwoodite (J. . In contrast, both Fe and H 2 O content greatly increase the V P /V S ratio of Fe-bearing ringwoodite (Sinogeikin et al., 2003;Jacobsen and Smyth, 2006;Higo et al., 2008;Schulze et al., 2018). As shown in Fig. 12d, the effects of Fe and H 2 O on Poisson's ratio of ringwoodite are almost the same as those in the V P /V S ratio. Water reduces Poisson's ratio of Fe-free ringwoodite, whereas Poisson's ratio of Fe-bearing ringwoodite was enhanced by Fe and H 2 O. These observations clearly show that the influence of water on the V P /V S ratio and Poisson's ratio is distinctly different between pure-Mg and Fe-bearing ringwoodite.

Rheological properties
The results from theoretical calculations and experimental observations demonstrate that the rheological property of rocks and minerals is frequently described by a power-law creep equation of the form (Mei and Kohlstedt, 2000a, b;Karato and Jung, 2003;Hirth and Kohlstedt, 2003) whereε is the axial strain rate; σ is the differential stress; d is the grain size; A is a constant; C H 2 O is the water content (in wt. % H 2 O); E and V are the activation energy and activation volume for creep, respectively; n is the stress exponent; p is the grain-size exponent; and r is the water exponent. When a deviatoric stress is applied to a polycrystalline sample at low stresses and high temperatures, plastic deformation occurs at a rate that is sensitive to grain size and only weakly dependent on the applied stress. Thus, the deformation mechanism described by Eq. (10) in materials with small grain sizes and/or at low stresses is dominated by diffusion creep (ε diff ) (n = 1, p = 2-3). By contrast, dislocation creep (ε disl , another important mechanism of deformation) is operated for materials with coarse grain size at high stress, in which stress sensitivity is higher and strain rate does not depend on grain size in Eq. (10) (n = 3-5, p = 0). Especially for olivine-rich  Sinogeikin et al., 2003;polycrystal: Li, 2003;Higo et al., 2008); hydrous ringwoodite single crystal (Jacobsen and Smyth, 2006;J. Wang et al., 2006;Schulze et al., 2018). Closed and open symbols denote the experimental data, and solid and dashed lines stand for the calculated results from third-order finite strain theory.
rocks with and without melt under anhydrous and hydrous conditions, p ≈ 3 in the diffusion creep region and n ≈ 3.5 in the dislocation creep regime (e.g., Mei and Kohlstedt, 2000a, b;Hirth and Kohlstedt, 2003). In addition, grain-boundary sliding (GBS) accommodated by dislocation motion has been identified as an important deformation mechanism in olivine (Hirth and Kohlstedt, 2003;Hansen et al., 2011). In most cases in the mantle, diffusion creep, dislocation creep, and GBS occur in parallel. Thus, the total strain rate of a sample is a simple sum of strain rates due to each mechanism: In general, deformation experiments are conducted under conditions close to a transition between these two mechanisms. For each mechanism, the flow law will depend on the water content. The nature of the transition between the "wet" flow law and the "dry" flow law is not very clear because it depends on the details of the microscopic mechanism of deformation, which includes diffusion creep, dislocation creep, and GBS. Therefore, separating the contributions from dif-ferent mechanisms by simultaneous inversion is important (Hirth and Kohlstedt, 2003).
sion creep and n = 3.5 for dislocation creep) and grain-size exponent (p = 3 for diffusion creep and p = 0 for dislocation creep) reported by Mei and Kohlstedt (2000a, b) are in good agreement with other deformation experiments (Karato et al., 1986;Kohlstedt, 1995, 2003;Zimmerman and Kohlstedt, 2004;Katayama and Karato, 2008a, b;Tasaka et al., 2015Tasaka et al., , 2016Faul et al., 2016;Ohuchi et al., 2017) and theoretical predictions (Coble, 1963). Compared with deformation of anhydrous olivine, the strain rate of hydrous olivine aggregates greatly increases with increasing water content (Mei and Kohlstedt, 2000a, b;Hirth and Kohlstedt, 2003;Katayama and Karato, 2008a, b;Faul et al., 2016) (Fig. 13b). Nonetheless, the values for the stress exponent and grain-size exponent obtained in hydrous olivine are comparable to those determined under dry conditions. In addition, the water exponent in the diffusion creep regime is 1, but it is 1.2 in the dislocation creep regime. This difference may reflect the importance of various waterderived point defects for the creep behavior of olivine (Mei and Kohlstedt, 2000a, b;Faul et al., 2016). It is worth noting that the creep rate during deformation is limited by the diffusivity of the slowest ionic species diffusing along its fastest path. Intrinsic point defects in olivine include tetrahedral (Si), octahedral (M), and oxygen vacancies and interstitials. Previous studies demonstrated that, under anhydrous conditions, metal vacancies and ferric iron on metal sites are the majority point defects with the dominant charge neutrality condition [Fe and Kohlstedt, 2000a, b;Hirth and Kohlstedt, 2003). Deformation of olivine is ratecontrolled by Si as the slowest diffusing ionic species (Hirth and Kohlstedt, 2003;Kohlstedt, 2006 In the case of the latter charge neutrality condition, dislocation climb is rate-limited by diffusion of Si occurring by an interstitial mechanism (Mei and Kohlstedt, 2000a, b). Remarkably, Fei et al. (2012Fei et al. ( , 2013Fei et al. ( , 2014 proposed that the effect of water on O and Si self-diffusion in olivine is very limited or even negligible and that water does not significantly affect the creep rate in olivine assuming the diffusioncontrolled deformation mechanism in the UM. However, as to the water effect on olivine rheology, Fei et al.'s observations based on diffusion experiments are completely opposite to the results of the deformation experiments (e.g., Mei and Kohlstedt, 2000a, b;Hirth and Kohlstedt, 2003;Katayama and Karato, 2008a, b;Faul et al., 2016).
Except for the creep flow law experiments of olivine discussed above, the fabric transition and slip system of olivine under anhydrous and hydrous conditions have been investigated extensively (e.g., Jung and Karato, 2001;Katayama et al., 2004;Jung et al., 2006;Ohuchi et al., 2012;Ohuchi and Irifune, 2013;Wallis et al., 2019). These experiments revealed that olivine has various types of lattice-preferred orientation and slip system under the different levels of water content and stress. In addition, Girard et al. (2013) and Tielke et al. (2017Tielke et al. ( , 2019 carried out shear deformation experiments to quantify the hydrolytic weakening of olivine single crystals. The brittle behavior of hydrous olivine was investigated by Ohuchi et al. (2018) through pure shear deformation at 1-1.8 GPa and 670-1250 K. Note that all of these experimental findings have important applications for the interpretation of seismic anisotropy, mechanism of shear localization, strength of the mantle lithosphere, and intermediate-depth earthquakes in subduction zones and in the deep UM (e.g., Jung and Karato, 2001;Katayama et al., 2004;Jung et al., 2006;Ohuchi and Irifune, 2013;Ohuchi et al., 2012Ohuchi et al., , 2018Girard et al., 2013;Tielke et al., 2017Tielke et al., , 2019Wallis et al., 2019).

Wadsleyite and ringwoodite
Because wadsleyite and ringwoodite are only stable at 410-520 km depth (13-18 GPa) and 520-660 km depth (16-23 GPa) of the MTZ, performing deformation experiments on wadsleyite and ringwoodite is very difficult under the P -T conditions of the MTZ. Generally, a conventional deformation apparatus such as the Griggs or Paterson apparatus can be operated only at low pressures (P < 3 GPa or P < 0.5 GPa, respectively). Yamazaki and Karato (2001) designed the rotational Drickamer apparatus, in which deformation experiments can be performed to large strain in simple shear geometry. Because anvils are better supported in this design, deformation experiments can be performed even under lower mantle conditions in conjunction with synchrotron X-radiation (Girard et al., 2016). Subsequently, the deformation DIA was developed and has been used up to ∼ 20 GPa and ∼ 1700 K (e.g., Wang et al., 2003;Kawazoe et al., 2010).
Several deformation experiments (e.g., Xu et al., 2003;Nishihara et al., 2008;Kawazoe et al., 2010) have been conducted to identify the operating deformation mechanisms and determine the flow law for each mechanism for water-poor polycrystalline wadsleyite and ringwoodite and to constrain the stress exponent n (Fig. 14a). However, due to the very high level of stress and limited experimental data points, the flow law for parameters in Eq. (10) was not well constrained quantitatively by these studies. With the improvement of experimental techniques, more data from a wide range of strain rates, temperature, pressure, grain size, and water content were reported to investigate the rheological behavior of the Earth's deep mantle. For example, Kawazoe et al. (2013) and Farla et al. (2015) demonstrated that water can significantly reduce the strength of wadsleyite under MTZ conditions. The slope of the stress-strain rate data for wadsleyite with 2100 wt. ppm H 2 O, n = 6 ± 3 reported by Hustoft et al. (2013), suggests that dislocation creep is the dominant deformation mechanism under wet conditions (Fig. 14a). Similar to wadsleyite, the dislocation creep strain rates in hydrous ringwoodite (Kawazoe et al., 2016) were markedly enhanced by water incorporation compared with that of anhydrous ringwoodite (Xu et al., 2003). In addition, the inferred strain rates for dislocation creep in wadsleyite and ringwoodite from experimental measurements of Si and O self-diffusion rates in these minerals (Shimojuku et al., 2009) under P -T conditions similar to those in deformation experiments are also plotted in Fig. 14a and b for comparison. The predicted dislocation creep strain rates for ringwoodite by Shimojuku et al. (2009) agree reasonably well with the experimental observations of Kawazoe et al. (2016). However, the results of dislocation creep strain rates between Shimojuku et al. (2009) and others (Hustoft et al., 2013;Farla et al., 2015) are inconsistent in wadsleyite (Fig. 14a), possibly due to the large uncertainties in the selection of the appropriate model parameters as suggested by Hustoft et al. (2013). Recently, Fei et al. (2017) estimated the viscosity of ringwoodite and bridgmanite under the conditions of the MTZ and lower mantle by measuring their dislocation mobilities within their crystal structures.

Links between transport properties
In this section, we will try to explore the quantitative relations between the different transport processes in silicate minerals. As discussed above, the propagation of elastic wave (compressional and shear wave velocity) and heat transfer (thermal conductivity/thermal diffusivity) are primarily caused by atomic displacements in the vicinity of equilibrium position (e.g., Hofmeister et al., 2007;Karato, 2008). However, three other transport properties, namely electrical conductivity, diffusivity, and rheology, are basically controlled by atomic or ionic transport, in which atomic displacements on spatial and temporal scales are larger than the typical displacements due to atomic vibration about a mean position (e.g., Chakraborty, 1995;Mehrer, 2007;Zhang, 2010;Ni et al., 2015). Here, we do not dig into the underlying physical principles but summarize the quantitative links between them.

Diffusivity and electrical conductivity
The electrical conductivity in minerals is used in studies of mineral structure and electronic and ionic transport processes. Ionic motion dominates electrical transport in many kinds of silicate glasses. Therefore, studying electrical transport is an important alternative approach to understanding diffusion behavior. If the electrical conduction of a material is dominated by the motion of a single ionic species, the relationship between electrical conductivity (σ ) and diffusion coefficient (D) is given by the Nernst-Einstein relation (Karato, 1990;Chakraborty, 1995): where z is the valence of the conducting species, e is the charge of the electron, c is the concentration of the migrating species, k B is Boltzmann's constant, and H R is the Haven ratio determined by the diffusion mechanism and usually has a value between 0.1 and 1. Hydrogen diffusion in NAMs may provide a different perspective on how proton conduction affects bulk electrical conductivity with information exclusively on hydrogen mobility. Combined with the Nernst-Einstein relation (Eq. 12) and experimentally determined hydrogen diffusivity in olivine Novella et al., 2017;Sun et al., 2019), wadsleyite Sun Figure 14. Plots of strain rate versus differential stress for (a) wadsleyite and (b) ringwoodite deformed under dry and wet conditions. Data sources: K13 (Kawazoe et al., 2013); H13 (Hustoft et al., 2013); F15 (Farla et al., 2015); S09 (Shimojuku et al., 2009); X03 (Xu et al., 2003); K16 (Kawazoe et al., 2016, and ringwoodite (Sun et al., 2015), these hydrogen diffusion results can be used to constrain the contribution of the total electrical conductivity of corresponding NAMs. As a result, a comparison between the indirect calculated conductivity and direct measured electrical conductivity as a function of water content is shown in Fig. 15a for olivine, Fig. 15b for wadsleyite, and Fig. 15c for ringwoodite, respectively. In this calculation, the isotropic electrical conductivity of polycrystalline olivine and its high-pressure polymorphs was approximately estimated by a geometric mean of the hydrogen diffusivities along each principle axis Novella et al., 2017). Owing to the relatively larger activation energy determined by hydrogen diffusion, the calculated proton conductivities are usually higher than those from laboratory conductivity measurements, except for the studies of Du Frane and Tyburczy (2012) for olivine and Sun et al. (2015) for ringwoodite. This is because electrical conductivity in a material is the sum of conduction of each charge carrier (or defect) type acting in parallel. In general, all defects are present in some amount in every crystal, and each contributes to the total conductivity. In most cases, only one or two types of defects dominate under a given set of thermodynamic conditions. If Eq. (12) is used to calculate the diffusion coefficient from conductivity measurements, one has to separate partial conductivities due to individual ions from total conductivity. By contrast, the electrical conductivity can be estimated from the diffusion coefficient in a similar way if the predominant ionic species (or defect) are the same between diffusion process and electrical conduction. For example, the electrical conductivities of olivine, wadsleyite, and ringwoodite calculated from hydrogen selfdiffusion in olivine and its high-pressure polymorphs were applied to constrain the possible water contents in the UM and MTZ (e.g., Karato, 1990;Sun et al., 2015Sun et al., , 2018Sun et al., , 2019.

Diffusivity and rheology
When a differential stress is applied to a polycrystalline material, a heterogeneous state is created due to heterogeneous mechanical properties. Grain boundaries are weaker than the grains themselves, and hence the mutual sliding of grains along grain boundaries occurs upon the application of stress. Because grain boundaries are a good sink/source for defects (Karato, 2008), the concentration of point defects near grain boundaries will be dependent upon the stress state at grain boundaries. When diffusion along grain boundaries is important, the effective diffusion coefficient should involve the contributions from volume (D V ) and grain-boundary (D GB ) diffusion: where δ and d are the thickness of the grain boundary and grain size, respectively. The diffusion creep (including Coble creep and Nabarro-Herring creep) and dislocation creep are believed to be driven by the diffusion of the slowest element, that is, Si in silicate minerals, on grain boundaries (for Coble creep) and within grain interiors (for Nabarro-Herring creep and dislocation creep) (Frost and Ashby, 1982;Weertman, 1999). Thus, plastic deformation rates at high temperature in minerals can be calculated from the slowest diffusing species in the major constituent elements. For example, creep due to diffusional mass transport through the bulk of grains (Nabarro-Herring creep) can be expressed as follows (Frost and Ashby, 1982): where α is a constant (≈ 13.3), σ is the stress, D is the diffusion coefficient, and is the molar volume. If the diffusion Figure 15. Comparison of laboratory conductivity data with the calculated conductivity from hydrogen diffusion results through Eq. (12) as a function of 1/T for different water contents in (a) olivine, (b) wadsleyite, and (c) ringwoodite. Figure 16. Comparison of laboratory strain rate versus C H 2 O with those calculated from Si diffusion data. All data are normalized to 8 GPa, 1600 K, and a stress of 300 MPa for olivine (a) and 16 GPa, 1873 K, and a stress of 10 MPa for wadsleyite (b), respectively. Si self-diffusivities in olivine (Costa and Chakraborty, 2008;Fei et al., 2013) and wadsleyite (Shimojuku et al., 2009(Shimojuku et al., , 2010Druzhbin, 2019) were converted to strain rate using Eqs. (14)- (17) with l g /l c = 1 and parameters reported therein. Note that the water exponent (≈ 0.5) in the calculation of strain rate (b) by Shimojuku et al. (2009Shimojuku et al. ( , 2010 is based on the geometric mean of Si self-diffusion in wadsleyite single crystal (Druzhbin, 2019). (14) is replaced with an effective diffusion coefficient (D eff ) given by Eq. (13), one obtains the following: If (π δ/d)D GB D V at small grain size and/or at low T , theṅ This type of creep mechanism in Eq. (16) is referred to as the Coble creep (Coble, 1963). In the case of dislocation creep, dislocations move continuously, responding to applied stress until they cannot glide due to the mutual interaction of the edge portions. Dislocation climb overcomes the obstacle by which the dislocation glide can operate again. Consequently, the plastic strain is accomplished by dislocation glide. In the Weertman model, the strain rate is written as follows (Weertman, 1999): where G is the shear modulus, b is the Burgers vector, l g is the glide distance, and l c is the climb distance. The ratio of l g /l c is one of the important parameters for describing the dislocation creep in the Weertman model. In general, the D V and D GB are the slowest diffusion rates of Si, which are known from diffusion experiments measured under P -T conditions similar to that in deformation experiments. Therefore, plastic deformation rates for diffusion creep and dislocation creep at high temperature and high pressure can be calculated from Eqs. (14) to (17) using the Si self-diffusion coefficient in olivine, wadsleyite, and ringwoodite (e.g., Karato, 2008;Shimojuku et al., 2009). As an example, Fig. 16 shows a comparison of strain rate from deformation experiments with the calculated strain rate by the Weertman model (Eq. 17). When all data are normalized to a pressure of 8 GPa, a temperature of 1600 K, and a stress of 300 MPa for olivine (Fig. 16a), it is apparent that the dislocation creep rates reported by deformation experiments (Mei and Kohlstedt, 2000a, b;Hirth and Kohlstedt, 2003) are about 1 to 3 orders of magnitude higher than those calculated from Si diffusion in olivine (Costa and Chakraborty, 2008;Fei et al., 2013) using Eqs. (14)- (17) with l g /l c = 1. In addition, owing to the difficulties in performing diffusion experiments under the MTZ conditions, few experiments have been reported on the diffusion of silicon in wadsleyite (Shimojuku et al., , 2009(Shimojuku et al., , 2010Druzhbin, 2019) and ringwoodite (Shimojuku et al., 2009). Among them, only Druzhbin (2019) measured the Si diffusion in wadsleyite single crystal as a function of water content. In this calculation, a water exponent (≈ 0.5) for Si diffusion was obtained by the geometric mean of Druzhbin's results from three principal crystal orientations, which is also comparable to that for Si diffusion in olivine (Costa and Chakraborty, 2008). Consequently, the calculated dislocation creep rates in wadsleyite (Fig. 16b) using a water exponent of 0.5 and parameters of Si diffusion from Shimojuku et al. (2009Shimojuku et al. ( , 2010 are lower than those determined by plastic deformation experiment (Kawazoe et al., 2013). This discrepancy could be due to the different stress conditions in deformation experiment (Mei and Kohlstedt, 2000a, b;Kawazoe et al., 2013;stress ∼ 100-540 MPa) and diffusion experiment (e.g., Costa and Chakraborty, 2008;Shimojuku et al., 2009Shimojuku et al., , 2010Fei et al., 2013;Druzhbin, 2019; stress is negligibly small because of hydrostatic pressure). In fact, the stress range in deformation experiment may exceed estimated stress conditions of the MTZ (∼ 10 MPa or less) by more than 1 to 2 orders of magnitude. Such high-stress conditions should have produced high dislocation densities, thus resulting in an apparently high water-content dependence of creep rate in olivine (ε ∝ (C H 2 O ) 1.2 , Hirth and Kohlstedt, 2003) and wadsleyite (ε ∝ (C H 2 O ) 2.7 , Kawazoe et al., 2013). If this is the case, the high water-content dependence of wadsleyite creep obtained by deformation experiments cannot be applied to the MTZ due to the low stresses and low dislocation density.

High-conductivity and low-velocity anomalies in the UM
In the UM, some studies have shown that the electrical conductivity in the asthenosphere just below the lithosphereasthenosphere boundary (LAB) is approximately 10 −2 to 10 −1 S/m, known as the high conductive layer (HCL, e.g., Evans et al., 2005;Naif et al., 2013). The conductivity largely increases with increasing depth from 50 to 100 km and also becomes strongly anisotropic in this depth range. Initially, the origin of the HCL has been attributed to the presence of hydrous olivine (e.g., Karato, 1990;. However, conductivity anisotropy in olivine single crystals appears to have an insufficient effect on electrical conductivity to account for the observed HCL (e.g., Yoshino et al., 2006;Poe et al., 2010;Yang, 2012), which instead may be explained by the presence of partial melt elongated in the direction of plate motion (e.g., Caricchi et al., 2011;Zhang et al., 2014;Pommier et al., 2015;Stern et al., 2015;Zhang and Yoshino, 2020). In addition, water dissolved in olivine is unlikely to produce the conductivity anomalies observed in the UM  because the concentration (more than 1000 ppm H 2 O) of water in minerals required to reach UM conductivities would lead to partial melting of the mantle rocks (Hirschmann, 2010). Another striking geophysical anomaly identified in the UM is the low-velocity zone (LVZ, e.g., Hirschmann, 2010), which is characterized by low seismic velocities and high attenuation and located in the asthenosphere near the LAB. Under the oceanic plate, the depth range of LVZ appears to coincide with that of the HCL in some regions, with a 20 to 50 km thick layer (e.g., Evans et al., 2005;Naif et al., 2013). As discussed earlier, the LVZ of the UM cannot be explained by the sound velocities of hydrous olivine. Alternatively, the reduced seismic velocity (LVZ) and elevated electrical con-62 B.-H. Zhang and Q.-K. Xia: Physical properties of olivine, wadsleyite, and ringwoodite ductivity (HCL) have been attributed to the presence of partial melting (e.g., Yoshino et al., 2010;Zhang et al., 2014;Chantel et al., 2016;Freitas et al., 2019), although large discrepancies still remain in the laboratory estimations of the amount of melt volume fraction present in the asthenosphere (e.g., Chantel et al., 2016;Freitas et al., 2019;Zhang and Yoshino, 2020).

Water in the MTZ
Considering the high water storage capacities of wadsleyite and ringwoodite (e.g., Inoue et al., 1995;Kohlstedt et al., 1996), the MTZ has been thought to be a potentially large water reservoir in the Earth's mantle. Contrast of water partition coefficient between wadsleyite and olivine at the UM boundary and between ringwoodite and bridgmanite at the lower mantle boundary would lead to dehydration melting near 410 and 660 discontinuities since excess water from the MTZ cannot be incorporated substantially by minerals with a smaller partition coefficient and lower water solubilities. These partial melts could produce low-velocity and highconductivity anomalies that can be imaged using seismic or electromagnetic observations (e.g., Bercovici and Karato, 2003;Song et al., 2004;Toffelmier and Tyburczy, 2007;Kawakatsu et al., 2009;Tauzin et al., 2010;Schmerr et al., 2012). Indeed, quantifying the amount of water in the Earth's mantle has been a long-standing challenge in deep Earth research integrating seismology, mineral physics, and geodynamics. On the basis of the above discussion in Sects. 2-6, water can dramatically affect the physical properties of mantle minerals including diffusion, electrical conduction, thermal conduction, sound velocity, and rheology. Therefore, the water content in the MTZ could be estimated according to the various transport properties of olivine, wadsleyite, and ringwoodite and their sensitivity to water.

From electrical conductivity
Electromagnetic studies demonstrated that high-conductivity anomaly is an important feature of the MTZ between depths of 400 and 660 km at many localities (Utada et al., 2003;Tarits et al., 2004;Ichiki et al., 2006;Shimizu et al., 2010), even globally (Kuvshinov and Olsen, 2006;Kelbert et al., 2009). For example, in the region beneath the North Pacific (Utada et al., 2003;Shimizu et al., 2010) and beneath northeastern China (Ichiki et al., 2001), the high-conductivity anomaly observed in the MTZ might be associated with water released from the stagnant slab. Recent electromagnetic studies have shown that the lateral variation in electrical conductivity is largely comparable to the depth variation (e.g., Ichiki et al., 2006;Shimizu et al., 2010), suggesting the strong sensitivity of electrical conductivity to some parameters that change laterally such as temperature and water content.
To better understand the nature of the high-conductivity zone and constrain the water content in the MTZ from the perspective of electrical conductivity, we construct a conductivity-depth profile in the Earth's mantle by extrapolating the experimental conductivity (olivine: Xu et al., 2000;wadsleyite: Huang et al., 2005;Dai and Karato, 2009;Yoshino et al., 2008;Yoshino and Katsura, 2012;ringwoodite: Huang et al., 2005;Yoshino et al., 2008) and hydrogen diffusion data (olivine: Sun et al., 2019;wadsleyite: Sun et al., 2018;ringwoodite: Sun et al., 2015) to the deep Earth in terms of temperature, pressure, and water content. In the present calculation, two sets of conductivity data (Yoshino and coworkers, termed the Yoshino model; Karato and coworkers -the Karato model) and hydrogen diffusion data (diffusion model) were used to calculate the electrical conductivity of the mantle by using Eqs. (1)- (4) and (12). The temperature in the MTZ is based on the average adiabatic mantle geotherm from . The electrical conductivity of dry olivine aggregates was referenced from Xu et al. (2000). Figure 17 shows a comparison of laboratory-based conductivity-depth profiles of the Earth's mantle for olivine, wadsleyite, and ringwoodite as a function of water content with geophysically observed electrical conductivity in the mantle. The conductivity-depth profiles obtained from the geophysical observations show significant differences from each other. For example, the semi-global reference model beneath the North Pacific Ocean (Utada et al., 2003;Shimizu et al., 2010) indicates that the electrical conductivity sharply jumps at the top of the MTZ (410 km discontinuity) and gradually increases to the top of the lower mantle. Remarkably, the recent global average model by Kelbert et al. (2009) is comparable to that beneath central Europe revealed by Tarits et al. (2004) in the UM and MTZ. However, both models are lower than that beneath northeastern China (Ichiki et al., 2006). By contrast, the laboratory-based conductivitydepth profiles are strongly influenced by the quantitative relationships between mineral conductivity and water content reported by different laboratories. For example, at 410 km discontinuity, the difference between olivine and wadsleyite is negligibly small in the Yoshino model (Fig. 17a), whereas a large drop in electrical conductivity at constant water content can be seen between olivine and wadsleyite according to the Karato model (Fig. 17b). Conversely, the conductivity-depth profiles derived from the diffusion model (Fig. 17c) show a positive conductivity jump at a depth of 410 km. However, when we compare the three laboratory-based models with the global average model by Kelbert et al. (2009) and the model beneath northeastern China (Ichiki et al., 2006), the estimated water content in the UM and MTZ is approximately 0.1 wt. %-0.5 wt. % H 2 O from the Yoshino and diffusion models, whereas the water content inferred from the Karato model is much less than 0.01 wt.% H 2 O in the UM and 0.1 wt.% H 2 O in the lower part of the MTZ. This discrepancy is mainly due to the difference in the dependence of water on the electrical conductivity of olivine, wadsleyite, Figure 17. Comparison of laboratory-based conductivity-depth profiles of the mantle with olivine composition as a function of water content (wt. %) for the adiabatic mantle geotherm  with the geophysically observed electrical conductivity in the mantle. (a) Profiles constructed from data by Yoshino and coworkers (olivine: wadsleyite: Yoshino and Katsura, 2012;ringwoodite: Yoshino et al., 2008). (b) Profiles constructed from data by Karato and coworkers (olivine: D. Wang et al., 2006;wadsleyite: Dai and Karato, 2009;ringwoodite: Huang et al., 2005). (c) Profiles constructed from H self-diffusion data (olivine: Sun et al., 20019;wadsleyite: Sun et al., 2018;ringwoodite: Sun et al., 2015). For dry olivine, data from Xu et al. (2000) are used. Global average (Kelbert et al., 2009); central Europe (Tarits et al., 2004); North Pacific (Utada et al., 2003;Shimizu et al., 2010); northeastern China (Ichiki et al., 2006). and ringwoodite investigated experimentally (Yoshino, 2010;Yoshino and Katsura, 2013).

From seismological observations
Seismic wave propagation is sensitive to water incorporation and has thus been proposed as promising geophysical observables to quantify hydration (e.g., Jacobsen and Smyth, 2006;Mao et al., 2011Mao et al., , 2012Houser, 2016). Seismological methods provide high three-dimensional resolution and offer a variety of observables that have been proposed to be sensitive to hydration, including lateral wave speed variations as depicted by seismic tomography and characteristics of seismic discontinuities (e.g., Houser and Williams, 2010;Schmandt et al., 2014;Houser, 2016;Buchen et al., 2018).
In this study, we attempt to constrain the water content and interpret seismic signatures in the potential water-rich region of the UM and MTZ by using the measured sound velocities of olivine, wadsleyite, and ringwoodite at simultaneous high pressure and high temperature by following the method used by Jacobsen and Smyth (2006) and Mao et al. (2012). The water content for hydrous (Mg 0.9 Fe 0.1 ) 2 SiO 4 olivine, wadsleyite, and ringwoodite is assumed to be 0.4 wt. %, 2.0 wt. % and 1.0 wt. % H 2 O, respectively, based on previous data (Li et al., 1996;Li, 2003;Sinogeikin et al., 2003;Jacobsen and Smyth, 2006;Higo et al., 2008;Wang et al., 2014;Mao et al., 2010Mao et al., , 2011Mao et al., , 2012Mao et al., , 2015Buchen et al., 2018). Figure 18 shows the modeled mantle velocity structures for olivine and its high-pressure polymorphs in the Earth's mantle as a function of depth and water content along an adiabatic mantle geotherm . As seen in the global seismic models (PREM and AK135 in Fig. 18a, b), the compressional and shear wave velocities exhibit discontinuous increases at certain depths (Dziewonski and Anderson, 1981;Kennett et al., 1995); some distinct features include the lowvelocity zone around 80-150 km, jumps at the top of 410 and 660 km depths, and a high-velocity gradient in the MTZ. Importantly, seismic low-velocity zones have been detected above the 410 km discontinuity in various regions (e.g., Song et al., 2004;Tauzin et al., 2010). Clearly, up to ∼ 3 % and ∼ 8 % decrease in V P and V S of olivine in the UM (Fig. 18c) are unlikely to be explained by olivine hydration because the sound velocities of olivine are insensitive to water as illustrated in Figs. 10 and 18a, b. In fact, these features were predicted by the "water filter" model (Bercovici and Karato, 2003), in which partial melting atop the 410 km discontinuity would occur with vertical transport of water-rich minerals from the MTZ to the UM. This hypothesis has been verified by recent experimental observations on seismic wave velocity of hydrous peridotites (β → α transition) (Freitas et al., 2017).
In contrast to the 410 and 660 km discontinuities, which are sharp features with reasonably well-defined velocity increases of approximately 5 % in V P and V S arising from the olivine-wadsleyite and ringwoodite- Figure 18. Influence of water on the mantle velocity structure. (a) Compressional (V P ) and (b) shear wave (V S ) velocities of olivine polymorphs in the Earth's mantle along an adiabatic mantle geotherm  compared with seismic models (PREM: Dziewonski and Anderson, 1981;AK135: Kennett et al., 1995). Red lines denote the velocity of anhydrous (Mg 0.9 Fe 0.1 ) 2 SiO 4 olivine polymorphs. Light blue lines represent the velocity of hydrous (Mg 0.9 Fe 0.1 ) 2 SiO 4 olivine, wadsleyite, and ringwoodite with 0.4 wt. %, 2.0 wt. % and 1.0 wt. % H 2 O, respectively, based on the data sources (Li et al., 1996(Li et al., , 2003Sinogeikin et al., 2003;Jacobsen and Smyth, 2006;Higo et al., 2008;Wang et al., 2014;Mao et al., 2010Mao et al., , 2011Mao et al., , 2012Mao et al., , 2015Buchen et al., 2018). (c) Comparison of water-induced reduction in velocity of (Mg 0.9 Fe 0.1 ) 2 SiO 4 olivine polymorphs in the Earth's mantle with the low-velocity anomalies in the mantle revealed by seismic observations in the UM (e.g., Song et al., 2004;Kawakatsu et al., 2009;Tauzin et al., 2010;Schmerr, 2012) and MTZ (Sieminski et al., 2003;Kustowski et al., 2008;Houser and Williams, 2010;Schmandt et al., 2014). bridgmanite + ferropericlase phase transitions (Dziewonski and Anderson, 1981;Kennett et al., 1995), the 520 km discontinuity is comparatively small, and there is less agreement on the magnitude of the velocity jump. The 520 km discontinuity may not be sharp (e.g., less than 10 km; Benz and Vidale, 1993), but it is perhaps a rather broad (25-50 km) velocity gradient (Revenaugh and Jordan, 1991). As shown in Figs. 11 and 18a, b, the influence of water on the sound velocities of wadsleyite is very limited. Even if wadsleyite contains 2.0 wt. % H 2 O, the V P and V S under the P -T conditions of the MTZ were reduced by only 0.9 %, and still higher those of the PREM and AK135 models (Fig. 18a, b). This finding suggests that the low-velocity anomalies with 1 %-3 % decrease in V S in the upper part (410-520 km) of the MTZ (e.g., Sieminski et al., 2003;Kustowski et al., 2008) cannot be explained by water alone, which may be caused by partial melting or lateral heterogeneities of composition or temperature in some locations. However, the magnitude of the reduction in seismic wave velocities of ringwoodite by a given amount of water is much greater than previously thought when the combined effect of hydration, temperature, and pressure is taken into consideration. In contrast, Schulze et al. (2018) demonstrated that the influence of hydration on the V P and V S of (Mg 0.89 Fe 0.11 ) 2 SiO 4 ringwoodite decreases with increasing pressure and vanishes under conditions of the MTZ. This observation suggests that hydrous ringwoodite (Fig. 18c) would be unable to account for the low-velocity anomaly in the lower part (520-660 km) of the MTZ derived from seismological studies (e.g., Sieminski et al., 2003;Kustowski et al., 2008;Houser and Williams, 2010;Schmandt et al., 2014). In other words, it is hard to constrain water content in the UM and MTZ by employing experimental measurements of V P and V S of olivine, wadsleyite, and ringwoodite (Fig. 18).

Rheology of the UM and MTZ
Diffusion and dislocation creep are considered to be the two main mechanisms of plastic deformation in the Earth's interior. An important question in geodynamics at present is which creep mechanism dominates in the UM and MTZ. If it is diffusion creep, shape-preferred orientation will be produced in olivine, and strain rate proportional to stress will result in Newtonian viscosity. By contrast, if dislocation creep dominates, lattice-preferred orientation will appear, and the strain rate proportional to stress to the power of 3.0-3.5 will lead to non-Newtonian viscosity (Fei et al., 2016). The mechanisms of weakening effects by water (hydrogen) were reviewed by Karato (2008). Given that high-temperature creep is controlled by diffusion-related processes, the enhancement of creep by water is partly due to enhanced diffusion (Costa and Chakraborty, 2008;Kawazoe et al., 2013).
Although deformation experiments can provide direct constraints on the details of creep mechanisms, strain rates attained in laboratories are usually on the order of 10 −5 s −1 , Figure 19. Grain-size/stress-deformation mechanism maps for (a) (Mg 0.9 Fe 0.1 ) 2 SiO 4 olivine (P = 7 GPa, T = 1500 K, and water-saturated conditions), (b) (Mg 0.9 Fe 0.1 ) 2 SiO 4 wadsleyite (P = 16 GPa, T = 1600 K, and C H 2 O = 1.0 wt. % H 2 O) and (c) (Mg 0.9 Fe 0.1 ) 2 SiO 4 ringwoodite (P = 20 GPa, T = 1800 K, and C H 2 O = 0.02 wt. % H 2 O). The power-law creep constitutive relations in olivine were established by Mei and Kohlstedt (2000a, b) and Hirth and Kohlstedt (2003) and in wadsleyite and ringwoodite calculated from the Si diffusion data (Shimojuku et al., 2009(Shimojuku et al., , 2010Druzhbin, 2019) using Eqs. (14)-(17) with l g /l c = 1. Note that dis GBS represents dislocation grain-boundary sliding. The thick light blue lines represent boundaries between different creep mechanisms. The orange, green, and dark purple boxes show the possible variations of grain size and stress in the laboratory, UM, MTZ and the subducting slabs, respectively. which is 10 orders of magnitude faster than the geological strain rates (∼ 10 −15 s −1 ) (see the colored boxes in Fig. 19). This discrepancy in strain rates is unavoidable because deformation experiments must be conducted on a timescale of hours. To enhance the strain rates, grain sizes are reduced to microns, and deviatoric stresses are increased to ∼ 100 MPa or more; by contrast, the grain size of mantle olivine is typically on the order of millimeters, and deviatoric stresses are expected to be on the order of 1 MPa in most of the mantle (e.g., Karato, 2008). Therefore, considerable extrapolation is involved when using experimentally derived mantle rheology in numerical modeling. However, the extrapolation of experimental data is valid only when extrapolation is within the same deformation mechanism where the same constitutive relationship determined in the laboratory would apply. At the same time, the published results on the plastic deformation of the UM and MTZ minerals should be critically evaluated, and various key parameters (temperature, water content, grain size) in the deep Earth should be estimated in order to assess the mechanism.
When several mechanisms of deformation operate, the use of a deformation mechanism map is a convenient way to understand which mechanism is dominant in a certain parameter space (Karato, 2008). On the basis of Eq. (10) and previous experimental data on olivine rheology (Mei and Kohlstedt, 2000a, b;Karato and Jung, 2003;Hirth and Kohlstedt, 2003;Ohuchi et al., 2017;Gouriet et al., 2019), Fig. 19a shows a deformation mechanism map of olivine on a stress versus grain-size space at a pressure of 7 GPa, a temperature of 1500 K, and water-saturated conditions. This deformation mechanism map for hydrous olivine demonstrates that diffusion creep dominates at lower stresses and smaller grain sizes, and dislocation creep becomes dominant at higher stresses and larger grain sizes. However, Fig. 19a indicates that dislocation GBS is one of the most effective deformation mechanisms at higher pressures and/or lower temperatures (e.g., subducting slabs), as suggested by Ohuchi et al. (2017).
To understand the dominant deformation mechanism in the MTZ, we have constructed the deformation mechanism maps for wadsleyite (Fig. 19b) at a pressure of 16 GPa, a temperature of 1600 K, and a fixed water content of 1.0 wt. % and for ringwoodite (Fig. 19c) at a pressure of 20 GPa, a temperature of 1800 K, and a fixed water content of 0.02 wt. %, respectively. The strain rates for diffusion and dislocation creep were calculated as a function of grain size and stress in terms of Eqs. (14)-(17) and Si self-diffusion data in wadsleyite (Shimojuku et al., , 2009(Shimojuku et al., , 2010Druzhbin, 2019) and ringwoodite (Shimojuku et al., 2009). As discussed in Sect. 7.2, the water exponent in the calculation of strain rate is assumed to be 0.5 for wadsleyite based on the geometric mean of Si self-diffusion in wadsleyite single crystal (Druzhbin, 2019). In this case, it is impossible to calculate the strain rate of ringwoodite as a function of water content, because only Shimojuku et al. (2009) measured Si diffusivities in ringwoodite at a single water content (∼ 0.02 wt. % H 2 O). Therefore, the influence of water on the deformation mechanism map for ringwoodite was not considered in this calculation. Our results imply that the dislocation creep mechanism may be dominant in the MTZ (Carrez et al., 2006;Ritterbex et al., 2015). Transitions between diffusion creep and dislocation creep mechanisms may easily occur through slight changes or heterogeneities in grain size and stress at the MTZ (Shimojuku et al., , 2009(Shimojuku et al., , 2010Druzhbin, 2019). In addition, the possible deformation mechanisms for subducting slabs are also illustrated in Fig. 19b, c. The rheological structure of the subducting slab depends strongly on the subduction parameters through temperature that controls the grain size of ringwoodite and the magnitude of forces acting on a slab . Furthermore, numerical modeling by Karato et al. (2001) shows that, regardless of warm or cold slabs, the subducting slabs lose their strength if the grain-size-sensitive diffusion creep operates because of the grain-size reduction after the olivine-ringwoodite transformation. However, the reliability of the conclusions from the deformation mechanism map depends on the reliability of the data used to construct such a map and also on the estimated values of parameters such as the stress and water-content exponent. Especially, the rheology of the UM, MTZ, and subducting slabs may strongly be affected by water content and chemical environment in the deep Earth. These issues need to be further constrained by deformation experiments with better controlled water fugacity.

Concluding remarks and perspectives
Despite its small abundance, water plays an important role in the properties of the Earth's interior, such as mantle convection, chemical differentiation, and seismic signature. In this article, we review the recent advances in the study of the transport properties of olivine, wadsleyite, and ringwoodite under high temperature and high pressure, specifically focusing on the influence of water. We have seen how water is involved in transport processes such as diffusion and deformation, the effect of water being to enhance the kinetics of transport properties (diffusion, electrical conductivity, and rheology) in olivine, wadsleyite, and ringwoodite through the increase in point defect concentrations, especially vacancies. In addition, the dissolution of water in the crystal lattice will weaken chemical bonding and increase structural disorder. Hence, one may expect that seismic wave velocities and thermal conductivity/diffusivity will be reduced by the incorporation of water.
Direct laboratory experimental measurements on the transport properties of NAMs can provide meaningful constraints water content in the mantle. However, these data have shown large discrepancies between research groups, especially for proton conduction in olivine and its high-pressure polymorphs (e.g., D. Yoshino et al., 2006;Huang et al., 2005;Yoshino et al., 2008;Gardés et al., 2015). As we discussed in Sect. 3, the primary reason for this inconsistency may originate from the different experimental techniques. Another typical example is to estimate the water content in the UM and MTZ by combining geophysical observations and the results of high-pressure experiments (see . However, different conclusions are obtained in Sect. 8. In most geophysical studies on the composition (especially for water) of the Earth's interior, the focus is seismological observations such as seismic wave velocities, which are rather insensitive to water content (e.g., Mao et al., 2010Mao et al., , 2011Buchen et al., 2018) but sensitive to other parameters such as the major element chemistry (e.g., pyrolite versus harzburgite) as shown by Karato (2011). In contrast to seismic wave velocity, electrical conductivity is more sensitive to water content compared with the composition (e.g., Karato, 1990;Yoshino et al., 2010;Yoshino and Katsura, 2013). Meanwhile, the temporal and spatial resolution of magnetotelluric surveys is much worse than that of seismology. Consequently, one cannot obtain robust constraints on water content from geophysical investigations (e.g., seismological observations or magnetotelluric surveys) or laboratory measurements (e.g., electrical conductivity or sound velocity) alone (e.g., Houser, 2016;Wang et al., 2019). This kind of difference between the different methods is not due to the experimental and/or geophysical methods themselves are wrong, but because each method has its own shortcomings or limitations.
The ultimate goal of laboratory measurements on the transport properties of mantle minerals is to explore universal law and construct a quantitative relationship for each transport property as a function of various thermodynamic parameters (T , P , f O 2 , H 2 O, etc.). Considering the substantial difficulties of performing experiments under extreme conditions, most experiments were conducted at very limited temperature and pressure conditions. Thus, different equations for the same property (e.g., Eq. 5 versus Eq. 6) were established by different groups, even using the same equation but obtaining different parameter values (see . As a result, not surprisingly, the fitting quality and model applicability are frequently unsatisfactory. In many cases, great caution must be exercised with the extrapolation/or interpolation of laboratory-based models to real Earth interior conditions. There is one point that we must emphasize: the connection between geophysical observations and physical properties of minerals is provided by the thermodynamics of mantle minerals and links between transport properties. As discussed above, enormous progress has been made over the last 2 decades in the ability to measure the physical properties of NAMs under high P -T conditions. Unlike other properties (i.e., diffusivity, electrical conductivity, and rheology), a unified law involving various thermodynamic parameters (T , P , f O 2 , etc.) has not been established for thermal conductivity and sound velocity yet. The parameterization of these two physical properties in terms of water content and pressure (and temperature) is urgently needed. Further exper-imental or theoretical studies are required to uncover these unresolved issues to better understand the nature of various geodynamic processes in the mantle.
Importantly, a comprehensive understanding of the composition, structure, and dynamics in the Earth, even a thorough interpretation of magnetotelluric or seismological profiles, needs a systematic multidisciplinary approach, including not only geophysical observations and geochemical analyses, but also laboratory measurements of the physical properties of mantle minerals under high temperature and pressure. Nonetheless, with each method being sensitive to specific parameters and probing different processes, multidisciplinary approaches are becoming more popular within the geoscience community and will cover a wide range of applications in the foreseeable future.
Data availability. All data derived from this research are presented in the figures and additional requests for raw data should be addressed to Bao-Hua Zhang.
Author contributions. BHZ and QKX designed the research, BHZ compiled and analyzed data, and BHZ and QKX wrote the paper.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Probing the Earth: reviews of OH groups in anhydrous and hydrous minerals". It is not associated with a conference.