# Equation of state and sound wave velocities of fayalite at high pressures and temperatures: implications for the seismic properties of the martian mantle

Frederic Bejina, Misha Bystricky, Nicolas Tercé, Matthew Whitaker, Haiyan Chen

► To cite this version:

Frederic Bejina, Misha Bystricky, Nicolas Tercé, Matthew Whitaker, Haiyan Chen. Equation of state and sound wave velocities of fayalite at high pressures and temperatures: implications for the seismic properties of the martian mantle. European Journal of Mineralogy, 2021, 33 (4), pp.519-535. 10.5194/ejm-33-519-2021. hal-03337044

HAL Id: hal-03337044

https://hal.science/hal-03337044

Submitted on 8 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Equation of state and sound wave velocities of fayalite at high pressures and temperatures: implications for the seismic properties of the martian mantle

Frédéric Béjina<fnref n="1" />, Misha Bystricky<fnref n="1" />, Nicolas Tercé<fnref n="1" />, Matthew L. Whitaker<fnref n="2" />,<fnref n="3" />, and Haiyan Chen<fnref n="2" />,<fnref n="3" />

**Correspondence:** Frédéric Béjina (frederic.bejina@irap.omp.eu)

Received: 9 April 2021 – Revised: 30 July 2021 – Accepted: 2 August 2021 – Published: 30 August 2021

**Abstract.** The elastic properties of a pure, synthetic fayalite aggregate were studied by coupled synchrotron X-ray diffraction and ultrasonic interferometry in a DIA-type multi-anvil press. Measurements at pressures up to about 7 GPa and temperatures up to 873 K yielded an adiabatic bulk modulus, $K_{S0} = 127.2 \pm 0.3$ GPa with $(\partial K_S/\partial P)_{T0} = 6.5 \pm 0.1$, and a shear modulus, $G_0 = 53.3 \pm 0.4$ GPa with $(\partial G/\partial P)_{T0} = 1.25 \pm 0.05$. When fixing $(\partial K_S/\partial P)_{T0} = 5.3$ (after $(\partial K_T/\partial P)_{T0}$ from Nestola et al., 2011), $K_{S0}$ increases to about 130 GPa. These estimates of $(K_{S0}, (\partial K_S/\partial P)_{T0})$ follow a general linear trend, $K = f(dK/dP)$, for fayalite. We define limited ranges for both bulk and shear moduli from previous studies, and we discuss how these variations affect seismic velocities and the determination of a mineralogical model in the context of the Mars InSight SEIS (Seismic Experiment for Interior Structure) experiment.

# 1 Introduction

The successful landing of the Mars InSight mission and more specifically the installation of the seismometer SEIS (Seismic Experiment for Interior Structure) (Lognonné et al., 2019) on the surface of Mars will bring new constraints on the deep interior of the red planet. A depth profile of the seismic velocities is one of the major goals of this mission because (1) it will provide the first direct picture of the inner structure of the planet, in particular the size of its core (Witze, 2021; Stähler et al., 2021), the presence or not of a lower mantle, and possible discontinuities in its upper mantle (a transition zone?); and (2) it will allow, for the first time, mineralogical models to be fitted to real martian data and therefore the mineralogy of Mars’ interior to be constrained.

Building such models requires knowing the thermoelastic properties of candidate minerals, i.e., minerals with compositions close to what is expected for Mars. For a solid solution like (Fe, Mg) olivines, unknown elastic parameters of a given composition are calculated from those of pure endmembers according to their relative amounts (e.g., Eq. 6 in Stixrude

and Lithgow-Bertelloni, 2005). As recently proposed, such a relation may not work for all olivine compositions (Béjina et al., 2019; Pamato et al., 2019). Consequently it is important to measure the elastic properties across the entire solid solution in order to clearly define a compositional trend. With this goal in mind we started performing experimental measurements of equations of state (EoS) and sound velocities on silicate minerals with an iron content higher than on Earth. In the text below we will use a simplified notation for the expressions of bulk and shear moduli. Unless mentioned otherwise the elastic parameters at constant temperature and room pressure, $K_{T0}$, $K_{S0}$, $G_0$ and their derivatives, $(\partial K_T/\partial P)_{T0}$, $(\partial K_S/\partial P)_{T0}$ and $(\partial G/\partial P)_{T0}$, will be denoted as $K_T$, $K_S$, $G$, $dK_T/dP$ (or $K'_T$), $dK_S/dP$ (or $K'_S$) and $dG/dP$ respectively.

The iron content of the martian mantle is estimated to be around 17 wt % (e.g., Dreibus and Wänke, 1985; Taylor, 2013); therefore, olivine, one of its major phases, is thought to have a molar Fe content ($X_{Fe}$) around 25%–30%. Very few values for bulk modulus ($K$) and shear modulus ($G$) exist for olivines with $X_{Fe} > 0.2$. If the trend of $G$ values

<footnote n="1">¹IRAP, Université de Toulouse, CNRS, CNES, UPS, 14 avenue Édouard Belin, 31400 Toulouse, France</footnote>

<footnote n="2">²Dept. of Geosciences, Stony Brook University, Stony Brook, NY 11794, USA</footnote>

<footnote n="3">³Brookhaven National Laboratory, Upton, NY 11973, USA</footnote>
with $X_{\text{Fe}}$ seems well defined and decreases linearly with increasing $X_{\text{Fe}}$ (e.g., Liu et al., 2005), the few results on K appear contradictory. Nestola et al. (2011) performed room-temperature static compression experiments and measured an isothermal bulk modulus, $K_T = 124.7$ GPa and $K'_T = 5.3$ for olivines with $X_{\text{Fe}} = 8\%, 20\%, 29\%, 38\%$, concluding that Fe has little effect on both K and $K'$ in this composition range. By contrast Chung (1971) determined a decreasing linear trend of the adiabatic bulk modulus ($K_S$) between forsterite ($K_S = 128.1$ GPa) and fayalite ($K_S = 122.0$ GPa) from sound velocity measurements. In addition data on pure fayalite are largely discrepant, with $K_T$ (see Fig. S4 and Table S6 in Béjina et al., 2019) ranging from about 104 GPa (Plymate and Stout, 1990, at 400 °C) to 138 GPa (Smyth, 1975) and $K_S$ values of 138 GPa (Sumino, 1979), 128 GPa (Graham et al., 1988) and 134 GPa (Isaak et al., 1993). First-principles calculations by Stackhouse et al. (2010) predicted an even larger bulk modulus for fayalite, above 139 GPa. Therefore, a Fe-compositional trend of the bulk modulus for Fe-rich olivines is difficult to define and may even be nonlinear (Béjina et al., 2019; Pamato et al., 2019).

This lack of constraints on the effect of iron on the elastic properties of the Fe–Mg olivine solid solution, in particular for $X_{\text{Fe}} > 0.2$, makes it difficult to adjust a mineralogical model to SEIS data and decipher the Fe content of the Mars mantle. We recently published high-pressure room-temperature X-ray diffraction (XRD) data (Béjina et al., 2019) and calculated the isothermal bulk modulus ($K_T$) for the following olivine compositions: ($Fe_x$, $Mg_{1-x}$)$_2$SiO$_4$ with $x = 0.45, 0.64, 0.82$ and 1. Because room-temperature measurements at high pressure are subjected to consequent differential stress ($t$), these were estimated by using the method of Singh et al. (1998) that allows for the determination of $d$ spacing values under hydrostatic pressure and the calculation of $t$. Subsequently the NaCl pressures were corrected for $t$. We present here the first results from an experiment at high temperatures and pressures above 1 GPa on a polycrystalline fayalite with coupled XRD and ultrasonic interferometry techniques (with the exception of Liu et al., 2010, but they did not calculate the elastic parameters). XRD measurements show that, as expected, high-temperature data present less differential stress than at room temperature, but $t$ can be either positive or negative. The same correction for stress is therefore applied here.

## 2 Experimental method

The main details of our experiment, performed at the superconducting wiggler beamline X17B2 at the National Synchrotron Light Source (Brookhaven National Laboratory, Uptown, NY), as well as the preparation and characteristics of our fayalite (Fa) aggregate, can be found in our previous paper (Béjina et al., 2019). The acoustic velocity measurements were obtained in conjunction with X-ray diffraction

(XRD) in situ at high pressure and temperature with the setup described by Whitaker et al. (2017). We will therefore only give here a brief description of the experiment.

### 2.1 Sample characteristics

Our $Fe_2SiO_4$ sample was synthesized from nanometric powders of $Fe_2O_3$ and $SiO_2$. The oxide powder mixture was annealed twice for 2 h at 1000 °C and $fO_2$ buffered by a gas mixture Ar/H$_2$O + Ar/H$_2$ close to the Fe/FeO equilibrium. After checking that the run product was pure fayalite by XRD, a dense pellet was prepared by spark plasma sintering (Béjina and Bystricky, 2009; Guignard et al., 2011). Observations of a fragment of this aggregate by scanning electron microscopy (SEM) using energy dispersive X-ray spectroscopy, orientation contrast and electron backscatter diffraction (EBSD) showed a very weak preferred crystallographic orientation and a homogeneous grain size around 2 µm. Finally our 1.2 mm diameter specimen was cored from this pellet, its length was brought to 1.2 mm, and both ends were carefully polished with $Al_2O_3$ powders.

### 2.2 High-pressure and high-temperature experiment

Beamline X17B2 at National Synchrotron Light Source (NSLS) was equipped with a 1000 ton DIA-type multi-anvil apparatus placed in the path of a white X-ray beam and with the DIASCoPE acoustic system for sound velocity measurements (Whitaker et al., 2017). Energy-dispersive XRD patterns (see as an example Fig. 3 in Béjina et al., 2019) were collected at a scattering angle of about 6.5° on the detection system composed of 10 detectors placed in a circle (nine between 0 and 180° azimuth plus one at 270°) allowing for the calculation of differential stress (see description in Weidner et al., 2010). We later used the stress values to correct the pressure obtained from NaCl XRD analysis. Our experiment was run using the cell assembly shown in Fig. 1 containing both the fayalite sample and NaCl as a pressure calibrant. A detailed description of the cell is given by Whitaker et al. (2017), the main characteristics being a graphite heater with a WRe3%-WRe25% thermocouple to measure temperature (not corrected for pressure effects); a mixture $10\text{NaCl} : 1\text{BN}$ at the back of the sample serving as both pressure marker using the NaCl EoS by Brown (1999) and as cushioning to reduce differential stress in the sample; and a BN sleeve along the sample helping to provide a pseudo-hydrostatic environment but also creating low-$fO_2$ conditions needed to keep fayalite stable. In order to maximize the coupling between the different parts of the cell assembly for better sound wave propagation, 1 µm thick gold foils are placed at each interface along the wave path (Fig. 1). The ones around the sample also serve as markers when measuring sample length on X-radiographic images (see Sect. 2.3). The beamline setup allows us to take images of the full sample at once, these being recorded on a Prosilica charge-coupled device (CCD)
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Cell assembly designed for ultrasonic interferometry measurements at high pressure and temperature (for details see Whitaker et al., 2017). The X-ray beam enters the cell via an edge of the boron-epoxy cube at 90° to the thermocouple. The dual-mode 10° Y-cut $LiNbO_3$ piezoelectric transducer, placed at the back of one of the anvils (grey), generates both P and S waves. $Al_2O_3$ serves as a buffer rod. To maximize coupling between the different parts on the wave paths, 1 µ thick gold foils are placed at each interface. They also serve as markers to measure the sample length on X-ray images. (b) Examples of P and S wave echoes of the different interfaces: anvil–$Al_2O_3$, $Al_2O_3$–fayalite and fayalite–NaCl. The travel time across the sample for each wave is calculated by cross-correlation of the echoes reflected from both ends.</figcaption></figure>

camera. A dual-mode 10° Y-cut $LiNbO_3$ piezoelectric transducer generated both P and S waves at resonance frequencies of 62 and 35 MHz respectively. XRD on NaCl and fayalite, P and S wave signals, and sample images were recorded sequentially at each step of the experiment.

The tonnage-temperature path is presented in Fig. S1 in the Supplement. The Fa aggregate was first cold-pressed to 35 tons and then heated at 1073 K for about 15 min in order to relax the accumulated stress. Because the furnace started to become unstable, the temperature was brought down to 873 K, and this remained the maximum temperature of the experiment. The maximum tonnage-temperature conditions were 50 tons and 873 K. XRD on Fa and NaCl, P and S wave data, and X-ray radiographic images of the sample were collected while decreasing temperature (usually by steps of 100 K). Once at room temperature tonnage was lowered, and temperature increased back to 873 K before collecting a new set of data. Similar steps were repeated until tonnage reached 5 tons. Note that at 5 and 7.5 tons, the maximum temperatures were 573 and 673 K respectively. As discussed in Sect. 3, the maximum pressure estimated from NaCl unit-cell volumes and the EoS of Brown (1999) and corrected for differential stress was about 7 GPa at the highest temperature and about 5 GPa at room temperature because of the relaxation of cell-assembly materials.

## 2.3 Measurements of sample length and ($V_p, V_s$)

The length of the Fa sample was determined from X-ray radiographic images taken during the experiment thanks to gold foils placed at each end of our sample that blocked X-rays and acted as markers. The length variation is therefore known as a function of pressure and temperature (Fig. S2). The uncertainty in length is around 2 pixels (e.g., Li et al., 2004) corresponding to about 5 µm.

The round-trip travel times of P and S waves were obtained by cross-correlation of the echoes reflected from the top and bottom of the Fa specimen. This gives a precise superposition of the two echoes, and their cross-correlation presents a unique maximum corresponding to the round-trip travel time across the sample (Fig. 2). To perform these calculations we built a Python script largely based on the `ObsPY` package (Krischer et al., 2015). Figure 3 shows the P and S velocities obtained with this method as a function of pressure. The main source of uncertainty in the velocities of P and S waves comes from the length measurements described above.

# 3 Results

In this section we will discuss both X-ray and ultrasonic interferometry results. It is well accepted that the latter approach is more reliable, in particular because it does not need an external pressure scale such as NaCl EoS (e.g., Kung et al., 2002). Nevertheless, a comparison of the two datasets provides clues on the consistency of our results.

We removed data acquired during the cold-pressing step of our experiment in both X-rays (already presented in Béjina et al., 2019) and sound velocity sets because they present the highest differential stress levels, and these may affect sound velocity measurements.

The microstructure of our sample was also described in our previous publication, and we will only emphasize that our fayalite remained stable throughout the experiment. XRD measurements and postmortem SEM-EBSD observations showed no signs of transformation in our specimen.

## 3.1 Stress, pressure and unit-cell volume determination

The 10-detector array and other high-pressure equipment we used at the NSLS beamline X17B2 have been described in detail by Weidner et al. (2010). The stress estimation was based on the measurement of lattice strains, a method detailed by Singh et al. (1998). We therefore refer readers to these previous publications (and references therein) for more information on evaluating differential stress in high-pressure experiments. For each fitted *hkl* peak the measured *d* spacings, *d*<sub>m</sub>, were determined from the energy peak position on all available detectors (nine detectors subend 180° and 22.5° apart and one at 270°). The values *d*<sub>p</sub>(*hkl*), representing *d* at hydrostatic pressure and used to estimate all unit-cell volumes, were calculated from the following relation (Singh
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Example of a cross-correlation using our ObsPy-based Python script. (a) Recorded $V_p$ data at 7 GPa and 874 K. The black and red echoes are from both ends of the sample. (b) Superposition of the red and black echoes resulting from the cross-correlation with the indication of the travel time, $\Delta t$. (c) Result of the cross-correlation showing the unique maximum that precisely determines the round-trip travel time of the P wave across the specimen.</figcaption></figure>

et al., 1998, their Eqs. 3 and 4):

$$d_m(hkl) = d_P(hkl) \left[ 1 + (1 - 3\cos^2\psi) Q(hkl) \right], \quad (1)$$

where $d_m(hkl)$ is the measured $d$ spacing of a plane (hkl), $\psi$ is the detector azimuth, and $Q(hkl)$ is a term proportional to differential stress ($t$):

$$Q(hkl) = \frac{t}{3} \left( \alpha [2G_R^X(hkl)]^{-1} + (1 - \alpha) [2G_V]^{-1} \right). \quad (2)$$

In the above equation, $G_R^X$ is the aggregate shear modulus calculated under the Reuss condition (iso-stress), $G_V$ is the shear modulus under the Voigt condition (iso-strain), and their respective expression is to be found in Singh et al. (1998) and Hearmon (1956). The weight factor $\alpha = 1$ (or = 0) corresponds to the pure Reuss (or pure Voigt) bound. Our estimation of $t$ was made assuming pure Reuss condition according to Chen et al. (2006).

Examples of plots, $d_m = f(1 - 3\cos^2\psi)$, are presented in Fig. S3. The value of $d_P(hkl)$ corresponds to $d$ when

$3\cos^2\psi = 1$. To estimate the differential stress both in NaCl and Fa we preferred an individual-peak fitting procedure instead of refining the unit-cell parameters by fitting the full measured diffraction patterns with modeled spectra. NaCl peaks were fitted using the software plot85 (Weidner et al., 2010), while fityk (Wojdyr, 2010) was used for Fa peaks. All unit-cell volumes were calculated using both the UnitCell software (Holland and Redfern, 1997) and our own set of Python scripts that make use of the LMFIT minimization package (Newville et al., 2014), also used for the stress calculations. Both procedures gave identical volumes within uncertainties, and UnitCell diagnostics helped filter poorly fitting peaks. In the end, unit-cell volumes of NaCl at high pressure and high temperature were calculated from the $d$ spacings of peaks (200), (220), (222), (400), (420) and (422), and for Fa we kept peaks (061), (112), (130), (131), (132), (133), (134), (152), (211), (222), (241) and (422).
<figure><img src="image_6.png" /><figcaption>Figure 3. Sound wave velocities as a function of pressure, $P_{Fa}$, for $V_p$ (a) and $V_s$ (b).</figcaption></figure>

Pressure was first estimated from the NaCl volumes using a third-order Birch–Murnaghan (BM) EoS (Birch, 1947) with the isothermal bulk modulus, $K_T$, $dK_T/dP$ and $dK_T/dT$, from Brown (1999), as recommended by Strässle et al. (2014), and thermal expansion from Spetzler et al. (1972). This pressure was then corrected for the differential stresses, $t$, measured in both NaCl and Fa (Fig. 4) following the next relation (Béjina et al., 2019):

$$P_{Fa} = P_{NaCl} - \frac{2}{3}(t_{Fa} - t_{NaCl}), \qquad (3)$$

where $P_{Fa}$ is the pressure felt by the Fa sample, $P_{NaCl}$ is the pressure calculated from NaCl volumes, and $t_{Fa}$ and $t_{NaCl}$ are stresses measured in the Fa sample and in NaCl respectively (detailed data are listed in Table S1). Stress was calculated from $d$ spacings measured at different azimuths of the detector array using Eq. (1) and following the method of Singh et al. (1998). As shown in Fig. 4, differential stress in Fa is mostly compressive at the highest pressures but mostly extensive below about 4 GPa. In comparison NaCl presents less variation in $t$, and consequently $P_{Fa}$ is slightly higher than $P_{NaCl}$ below 4 GPa and lower above (Fig. S4). Table 1 lists unit-cell parameters calculated from the above $d_p$ values as a function of pressure and temperature.

## 3.2 Uncertainties

Uncertainties ($1\sigma$) in individual peak positions, hence on $d_{hkl}$ spacing values ($\Delta d_{hkl}$), were directly obtained from the Levenberg–Marquardt fit of the diffractograms in fityk. Because there is no straightforward way to include error bars on data in LMFIT (one can use weighting factors, but there are different ways to estimate them), we used the $d-\Delta d$ and $d+\Delta d$ of each hkl peak to calculate the minimum and maximum values of cell parameters and volumes. The uncertainty in $V$ ($\Delta V$) for each ($P, T$) condition was defined as the maximum

value between $V-V_{min}$ and $V_{max}-V$, with $V$ being the nominal volume calculated from $d$.

Uncertainties in $P_{NaCl}$, $t_{Fa}$ and $t_{NaCl}$ (Table S1) were obtained from those on unit-cell volumes and $d$ spacings by error propagation. Uncertainties in $P_{Fa}$ were then determined according to Eq. (3). These were finally included in the weighted fitting routine of EosFit7 (Angel et al., 2014) in order to obtain the fayalite EoS parameters.

This approach leads to larger errors in $V$ and subsequently in all parameters derived from either NaCl or fayalite volumes compared to a purely statistical approach (for example a weighted fit) but could very well be more representative of the real $\Delta V$. Indeed, in many high-pressure studies, underestimated errors are often attributed to pressure and temperature mainly because pressure and thermal gradients are poorly known, the position of the thermocouple with respect to the sample could change with pressure and/or temperature, uncertainties in the EoS used to estimate pressure are not taken into account, etc.

## 3.3 Isothermal equation of state of fayalite

Unit-cell volumes (Table 1) as a function of stress-corrected pressure and temperature were fitted simultaneously in EosFit7. We used both the second- and third-order Birch–Murnaghan (BM) EoS with a linear expression of $dK_T/dT$. Because high-temperature data measured at high pressure rarely have the precision of those at 1 atm, the modified Holland–Powell expression was chosen to calculate the thermal expansion (Angel et al., 2014). An acceptable value of a fitted $K'_T$, 5.5 for $K_T = 125.2$ GPa, could be obtained after removing the room-temperature data, but uncertainties in all EoS parameters are large. Therefore, for the third-order BM EoS we fixed $K'_T$ at 5.3 (Nestola et al., 2011) to accommodate the small pressure range. All results are listed in Table 2, and one can see that $K_T$ values are compatible with our previous value, $K_T = 131.9 \pm 1.5$ GPa (Béjina et al., 2019), within the error bars. This $K_T$ was obtained using all room-temperature data, including the pressurization cycle, and therefore with a much higher stress (Fig. 4). The compatibility between room-temperature and high-temperature data, with different stress levels in both, is a hint that our pressure correction leads to reliable results.

As mentioned above high-pressure experiments rarely have precise temperature measurements. Consequently the thermal parameters are not fully reliable and present large uncertainties. Our expression of the thermal expansion tends to overestimate volumes at high temperature compared to Kroll et al. (2012), a point that we discuss in Sect. 4.1. To circumvent this difficulty and test the validity of our values for $K_T$, we fitted data at each temperature independently using a third-order BM EoS with $K'_T = 5.3$ (Fig. 5). Obtained values of $V_0$ and $K_T$ extrapolated to room pressure are listed in Table 3 as a function of temperature. From these values, an overall $K_{T(0,298)} = 129 \pm 5$ GPa was calculated from a lin-
<figure><img src="image_7.png" /><figcaption>**Figure 4.** (a) Differential stress calculated in both NaCl and fayalite using the method of Singh et al. (1998) (compressive stresses are marked as positive). (b) Pressure, $P_{Fa}$ (left axis and grey squares), determined from the NaCl unit-cell volumes and corrected for the differential stress. Right axis and red dots show temperature. Data from the cold-pressing step were not included in the determination of the elastic parameters, as discussed in the text.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Cell volume of our Fa sample as a function of pressure, $P_{Fa}$, at all temperatures. $P_{Fa}$ is corrected from differential stress as described by Eq. (3). Each fit is an isothermal third-order Birch–Murnaghan equation of state (see results in Table 3).</figcaption></figure>

ear fit, $K_T = f(T)$, and $V_{(0,298)} = 307.9 \pm 0.2$ Å³ was determined from the modified Holland-Powell equation (Angel et al., 2014). Despite the dispersion in the data we found a $K_T(0,298)$ value compatible with the EoS Fit 7 results within the uncertainty, therefore validating the simultaneous fit to the whole dataset.

### 3.4 $K_s$, $dK_s/dP$, $G$ and $dG/dT$ for fayalite

Following the procedure described, for example, in Kung et al. (2002) and using finite strain equations from Davies and Dziewoński (1975), we calculated $K_{S0}$, $G_0$ and their pressure derivatives at each temperature (Table 4). This method only uses the density of the sample and sound wave velocities and does not need an external pressure gauge such as NaCl (see also Li and Zhang, 2005). The formulation of Kung et al. (2002), thereafter named K02, is a truncated version of Davies and Dziewoński (1975), DD75, this latter giving access to second derivatives of $K_s$ and $G$ with respect to pressure. Our data do not have the resolution necessary to obtain reliable values for these, but they nevertheless allow for the calculation of pressure. A comparison of the pressure calculated by the finite strain equations with $P_{Fa}$ (Eq. 3) shows a remarkable agreement (Fig. 6). One can also see in Fig. 6 that uncorrected pressure is more discrepant. It is therefore not surprising that fitting $G = \rho V_s^2$ and $K_s = \rho V_p^2 - 4G/3$ to the data at each temperature using the NaCl pressure scale gives results very close to those from the finite strain method (Table 4).

Table 5 lists the values of all elastic parameters at 298 K calculated after a linear fit of data from Table 4 as a function of temperature. While we list the results from both DD75 and K02 formulations, we only retain the values from the K02 equations because of the large uncertainties around the DD75 results. One can nevertheless note that the two for-
<table><caption>**Table 1.** Unit-cell parameters of Fa and velocities $V_p$ and $V_s$ measured as a function of $P_{Fa}$ ($P_{NaCl}$ corrected for differential stress) and temperature. First column is the Fa XRD pattern numbers. Uncertainties are $1\sigma$.</caption><thead><tr><th></th><th>P<sub>Fa</sub><br>(GPa)</th><th>Temperature<br>(K)</th><th>V<br>(&Aring;<sup>3</sup>)</th><th>a<br>(&Aring;)</th><th>b<br>(&Aring;)</th><th>c<br>(&Aring;)</th><th>V<sub>p</sub><br>(km s<sup>-1</sup>)</th><th>V<sub>s</sub><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>20</td><td>4.91(43)</td><td>877(44)</td><td>304.01(18)</td><td>4.8121(37)</td><td>10.3853(54)</td><td>6.0834(36)</td><td>6.90(3)</td><td>3.47(2)</td></tr><tr><td>22</td><td>4.28(34)</td><td>672(34)</td><td>302.77(17)</td><td>4.8044(37)</td><td>10.3745(51)</td><td>6.0744(34)</td><td>6.95(3)</td><td>3.50(2)</td></tr><tr><td>24</td><td>3.69(29)</td><td>471(24)</td><td>301.60(15)</td><td>4.7980(32)</td><td>10.3663(45)</td><td>6.0641(31)</td><td>7.00(3)</td><td>3.53(2)</td></tr><tr><td>26</td><td>3.24(24)</td><td>302(15)</td><td>300.67(16)</td><td>4.7928(36)</td><td>10.3569(51)</td><td>6.0573(31)</td><td>7.06(3)</td><td>3.58(2)</td></tr><tr><td>28</td><td>4.82(41)</td><td>772(39)</td><td>303.37(17)</td><td>4.8074(34)</td><td>10.3767(52)</td><td>6.0814(36)</td><td>6.96(3)</td><td>3.49(2)</td></tr><tr><td>30</td><td>4.12(32)</td><td>571(29)</td><td>302.39(18)</td><td>4.8014(37)</td><td>10.3711(56)</td><td>6.0727(36)</td><td>7.00(3)</td><td>3.52(2)</td></tr><tr><td>32</td><td>3.54(27)</td><td>373(19)</td><td>301.20(19)</td><td>4.7951(32)</td><td>10.3588(56)</td><td>6.0630(40)</td><td>7.05(3)</td><td>3.56(2)</td></tr><tr><td>34</td><td>3.34(27)</td><td>302(15)</td><td>300.81(25)</td><td>4.7919(43)</td><td>10.3599(72)</td><td>6.0597(48)</td><td>7.05(3)</td><td>3.57(2)</td></tr><tr><td>38</td><td>6.76(65)</td><td>874(44)</td><td>300.49(20)</td><td>4.7930(38)</td><td>10.3359(61)</td><td>6.0658(41)</td><td>7.06(4)</td><td>3.49(2)</td></tr><tr><td>39</td><td>6.82(64)</td><td>872(44)</td><td>300.14(18)</td><td>4.7909(35)</td><td>10.3328(57)</td><td>6.0631(39)</td><td>7.06(4)</td><td>3.49(2)</td></tr><tr><td>42</td><td>6.47(56)</td><td>775(39)</td><td>299.65(19)</td><td>4.7886(39)</td><td>10.3284(60)</td><td>6.0588(39)</td><td>7.08(4)</td><td>3.51(2)</td></tr><tr><td>44</td><td>6.08(51)</td><td>672(34)</td><td>299.29(20)</td><td>4.7863(36)</td><td>10.3261(62)</td><td>6.0555(44)</td><td>7.10(4)</td><td>3.52(2)</td></tr><tr><td>46</td><td>5.79(46)</td><td>574(29)</td><td>298.60(19)</td><td>4.7828(42)</td><td>10.3194(61)</td><td>6.0496(37)</td><td>7.12(4)</td><td>3.53(2)</td></tr><tr><td>48</td><td>5.47(41)</td><td>473(24)</td><td>298.11(19)</td><td>4.7805(45)</td><td>10.3125(65)</td><td>6.0459(37)</td><td>7.15(4)</td><td>3.55(2)</td></tr><tr><td>50</td><td>5.18(37)</td><td>373(19)</td><td>297.54(16)</td><td>4.7788(38)</td><td>10.3055(54)</td><td>6.0407(32)</td><td>7.17(4)</td><td>3.58(2)</td></tr><tr><td>52</td><td>4.96(34)</td><td>302(15)</td><td>297.26(17)</td><td>4.7773(40)</td><td>10.3042(58)</td><td>6.0382(33)</td><td>7.20(4)</td><td>3.59(2)</td></tr><tr><td>56</td><td>6.55(64)</td><td>874(44)</td><td>300.70(18)</td><td>4.7965(38)</td><td>10.3363(59)</td><td>6.0651(38)</td><td>7.07(3)</td><td>3.50(2)</td></tr><tr><td>58</td><td>6.17(54)</td><td>774(39)</td><td>300.28(18)</td><td>4.7932(39)</td><td>10.3342(55)</td><td>6.0626(37)</td><td>7.09(3)</td><td>3.52(2)</td></tr><tr><td>60</td><td>5.83(49)</td><td>674(34)</td><td>299.67(17)</td><td>4.7901(42)</td><td>10.3292(57)</td><td>6.0568(32)</td><td>7.09(4)</td><td>3.52(2)</td></tr><tr><td>62</td><td>5.49(45)</td><td>573(29)</td><td>299.16(16)</td><td>4.7863(37)</td><td>10.3261(52)</td><td>6.0535(33)</td><td>7.13(4)</td><td>3.54(2)</td></tr><tr><td>64</td><td>5.12(38)</td><td>472(24)</td><td>298.64(17)</td><td>4.7842(38)</td><td>10.3210(54)</td><td>6.0474(33)</td><td>7.15(4)</td><td>3.56(2)</td></tr><tr><td>66</td><td>4.83(37)</td><td>373(19)</td><td>298.12(17)</td><td>4.7812(44)</td><td>10.3180(59)</td><td>6.0439(31)</td><td>7.19(4)</td><td>3.59(2)</td></tr><tr><td>68</td><td>4.64(33)</td><td>302(15)</td><td>297.79(16)</td><td>4.7801(37)</td><td>10.3132(49)</td><td>6.0403(29)</td><td>7.21(4)</td><td>3.61(2)</td></tr><tr><td>72</td><td>5.96(58)</td><td>873(44)</td><td>301.76(17)</td><td>4.8018(37)</td><td>10.3525(53)</td><td>6.0705(36)</td><td>7.02(3)</td><td>3.49(2)</td></tr><tr><td>74</td><td>5.58(51)</td><td>773(39)</td><td>301.35(20)</td><td>4.7979(41)</td><td>10.3534(62)</td><td>6.0677(42)</td><td>7.03(3)</td><td>3.50(2)</td></tr><tr><td>76</td><td>5.23(46)</td><td>672(34)</td><td>300.63(17)</td><td>4.7961(41)</td><td>10.3432(54)</td><td>6.0611(33)</td><td>7.03(3)</td><td>3.51(2)</td></tr><tr><td>78</td><td>4.89(40)</td><td>573(29)</td><td>300.26(19)</td><td>4.7919(37)</td><td>10.3434(58)</td><td>6.0589(40)</td><td>7.07(3)</td><td>3.53(2)</td></tr><tr><td>80</td><td>4.55(34)</td><td>473(24)</td><td>299.70(17)</td><td>4.7906(41)</td><td>10.3357(56)</td><td>6.0528(30)</td><td>7.10(3)</td><td>3.55(2)</td></tr><tr><td>82</td><td>4.25(33)</td><td>372(19)</td><td>299.17(18)</td><td>4.7871(42)</td><td>10.3321(60)</td><td>6.0490(34)</td><td>7.13(4)</td><td>3.58(2)</td></tr><tr><td>84</td><td>4.03(30)</td><td>301(15)</td><td>299.11(20)</td><td>4.7873(45)</td><td>10.3318(64)</td><td>6.0478(39)</td><td>7.11(4)</td><td>3.58(2)</td></tr><tr><td>88</td><td>5.33(48)</td><td>877(44)</td><td>303.29(19)</td><td>4.8094(41)</td><td>10.3724(57)</td><td>6.0799(38)</td><td>6.98(3)</td><td>3.49(2)</td></tr><tr><td>90</td><td>4.91(44)</td><td>772(39)</td><td>302.81(21)</td><td>4.8065(44)</td><td>10.3690(62)</td><td>6.0764(40)</td><td>7.00(3)</td><td>3.50(2)</td></tr><tr><td>92</td><td>4.57(39)</td><td>672(34)</td><td>302.37(19)</td><td>4.8043(44)</td><td>10.3655(59)</td><td>6.0720(38)</td><td>6.99(3)</td><td>3.50(2)</td></tr><tr><td>94</td><td>4.22(33)</td><td>572(29)</td><td>301.77(21)</td><td>4.8011(44)</td><td>10.3615(65)</td><td>6.0671(43)</td><td>7.02(3)</td><td>3.52(2)</td></tr><tr><td>96</td><td>3.90(32)</td><td>471(24)</td><td>301.25(23)</td><td>4.7980(47)</td><td>10.3600(65)</td><td>6.0620(43)</td><td>7.06(3)</td><td>3.54(2)</td></tr><tr><td>98</td><td>3.60(26)</td><td>372(19)</td><td>300.85(18)</td><td>4.7958(40)</td><td>10.3559(52)</td><td>6.0585(38)</td><td>7.08(3)</td><td>3.56(2)</td></tr><tr><td>100</td><td>3.42(24)</td><td>302(15)</td><td>300.35(19)</td><td>4.7940(44)</td><td>10.3478(59)</td><td>6.0550(36)</td><td>7.10(3)</td><td>3.58(2)</td></tr><tr><td>104</td><td>4.57(48)</td><td>876(44)</td><td>304.90(19)</td><td>4.8158(37)</td><td>10.3983(57)</td><td>6.0900(40)</td><td>6.91(3)</td><td>3.47(2)</td></tr><tr><td>106</td><td>4.20(41)</td><td>773(39)</td><td>304.36(19)</td><td>4.8128(40)</td><td>10.3919(57)</td><td>6.0857(39)</td><td>6.92(3)</td><td>3.48(2)</td></tr><tr><td>108</td><td>3.86(34)</td><td>674(34)</td><td>303.90(17)</td><td>4.8099(36)</td><td>10.3919(50)</td><td>6.0810(36)</td><td>6.95(3)</td><td>3.50(2)</td></tr><tr><td>110</td><td>3.52(31)</td><td>572(29)</td><td>303.32(18)</td><td>4.8073(40)</td><td>10.3840(56)</td><td>6.0764(35)</td><td>6.97(3)</td><td>3.51(2)</td></tr><tr><td>112</td><td>3.23(25)</td><td>472(24)</td><td>302.59(16)</td><td>4.8025(32)</td><td>10.3791(47)</td><td>6.0713(33)</td><td>6.97(3)</td><td>3.52(2)</td></tr><tr><td>116</td><td>2.76(20)</td><td>302(15)</td><td>301.82(18)</td><td>4.7997(37)</td><td>10.3720(53)</td><td>6.0636(36)</td><td>7.03(3)</td><td>3.57(2)</td></tr><tr><td>120</td><td>3.57(33)</td><td>872(44)</td><td>307.07(19)</td><td>4.8264(38)</td><td>10.4250(53)</td><td>6.1035(39)</td><td>6.81(3)</td><td>3.44(2)</td></tr><tr><td>122</td><td>3.22(27)</td><td>773(39)</td><td>306.48(18)</td><td>4.8230(39)</td><td>10.4197(51)</td><td>6.0990(36)</td><td>6.82(3)</td><td>3.45(2)</td></tr><tr><td>124</td><td>2.92(22)</td><td>671(34)</td><td>305.78(18)</td><td>4.8193(39)</td><td>10.4139(53)</td><td>6.0930(35)</td><td>6.82(3)</td><td>3.46(2)</td></tr><tr><td>126</td><td>2.61(20)</td><td>571(29)</td><td>305.14(16)</td><td>4.8153(33)</td><td>10.4095(48)</td><td>6.0880(34)</td><td>6.86(3)</td><td>3.48(2)</td></tr><tr><td>128</td><td>2.31(18)</td><td>472(24)</td><td>304.59(18)</td><td>4.8116(36)</td><td>10.4087(53)</td><td>6.0832(37)</td><td>6.88(3)</td><td>3.50(2)</td></tr><tr><td>130</td><td>2.05(17)</td><td>373(19)</td><td>304.05(17)</td><td>4.8093(34)</td><td>10.4003(49)</td><td>6.0788(36)</td><td>6.88(3)</td><td>3.51(2)</td></tr><tr><td>132</td><td>1.86(14)</td><td>302(15)</td><td>303.58(18)</td><td>4.8075(35)</td><td>10.3961(52)</td><td>6.0744(36)</td><td>6.92(3)</td><td>3.54(2)</td></tr><tr><td>136</td><td>2.55(22)</td><td>675(34)</td><td>306.60(18)</td><td>4.8223(32)</td><td>10.4288(52)</td><td>6.0967(38)</td><td>6.79(3)</td><td>3.45(2)</td></tr><tr><td>138</td><td>2.20(19)</td><td>573(29)</td><td>306.17(18)</td><td>4.8190(37)</td><td>10.4275(53)</td><td>6.0936(37)</td><td>6.80(3)</td><td>3.47(2)</td></tr><tr><td>140</td><td>1.91(16)</td><td>473(24)</td><td>305.48(16)</td><td>4.8157(32)</td><td>10.4187(48)</td><td>6.0889(34)</td><td>6.83(3)</td><td>3.48(2)</td></tr><tr><td>142</td><td>1.60(14)</td><td>373(19)</td><td>304.96(18)</td><td>4.8132(32)</td><td>10.4161(52)</td><td>6.0836(37)</td><td>6.85(3)</td><td>3.51(2)</td></tr><tr><td>144</td><td>1.40(13)</td><td>301(15)</td><td>304.56(18)</td><td>4.8108(35)</td><td>10.4132(54)</td><td>6.0805(37)</td><td>6.86(3)</td><td>3.52(2)</td></tr><tr><td>148</td><td>1.77(14)</td><td>573(29)</td><td>307.31(21)</td><td>4.8239(38)</td><td>10.4434(60)</td><td>6.1006(42)</td><td>6.78(3)</td><td>3.47(2)</td></tr><tr><td>150</td><td>1.44(13)</td><td>473(24)</td><td>306.80(20)</td><td>4.8222(35)</td><td>10.4394(57)</td><td>6.0952(41)</td><td>6.80(3)</td><td>3.48(2)</td></tr><tr><td>152</td><td>1.14(10)</td><td>373(19)</td><td>306.23(20)</td><td>4.8187(37)</td><td>10.4345(59)</td><td>6.0910(41)</td><td>6.82(3)</td><td>3.49(2)</td></tr><tr><td>154</td><td>0.94(12)</td><td>301(15)</td><td>305.81(19)</td><td>4.8159(33)</td><td>10.4319(56)</td><td>6.0875(42)</td><td>6.80(3)</td><td>3.50(2)</td></tr><tr><td>156</td><td>-0.020(52)</td><td>301(15)</td><td>307.84(19)</td><td>4.8248(37)</td><td>10.4604(56)</td><td>6.1007(38)</td><td>not measured</td><td></td></tr></tbody></table>
<table><caption>**Table 2.** Birch–Murnaghan EoS parameters determined from our XRD data at all temperatures using EoSFit7 with the modified Holland–Powell equation for thermal expansion (Angel et al., 2014).</caption><thead><tr><th>V<sub>0</sub><br>(Å<sup>3</sup>)</th><th>K<sub>T</sub><br>(GPa)</th><th>dK<sub>T</sub>/dP</th><th>dK<sub>T</sub>/dT<br>(GPa K<sup>-1</sup>)</th><th>d<sup>2</sup>K<sub>T</sub>/dP<sup>2</sup></th><th>α<sub>0</sub><br>(10<sup>-5</sup> K<sup>-1</sup>)</th><th>α<sub>1</sub><br>(10<sup>-4</sup> K<sup>-1/2</sup>)</th></tr></thead><tbody><tr><td>307.83 ± 0.20</td><td>131.3 ± 5.1</td><td>4<sup>b</sup></td><td>-0.0061 ± 0.0274</td><td>-0.02962<sup>a</sup></td><td>6.541 ± 3.503</td><td>-1.745 ± 3.424</td></tr><tr><td>307.86 ± 0.21</td><td>128.6 ± 5.1</td><td>5.3<sup>b</sup></td><td>-0.0133 ± 0.0273</td><td>-0.05348<sup>a</sup></td><td>6.969 ± 3.639</td><td>-1.555 ± 3.517</td></tr></tbody></table>

<table><caption>**Table 3.** Third-order Birch–Murnaghan EoS parameters obtained from a weighted least-square fit of XRD data at each temperature; $dK_T/dP$ was fixed at 5.3 (Nestola et al., 2011).</caption><tr><td>Temperature (K)</td><td>302</td><td>373</td><td>472</td><td>572</td><td>673</td><td>773</td><td>874</td></tr><tr><td>V<sub>0</sub> (Å<sup>3</sup>)</td><td>307.84 ± 0.05</td><td>308.90 ± 0.13</td><td>310.24 ± 0.09</td><td>311.60 ± 0.16</td><td>312.78 ± 0.17</td><td>314.38 ± 0.32</td><td>315.90 ± 0.22</td></tr><tr><td>K<sub>T</sub> (GPa)</td><td>130.4 ± 1.4</td><td>125.3 ± 2.8</td><td>121.5 ± 1.7</td><td>121.6 ± 2.7</td><td>121.7 ± 2.4</td><td>119.3 ± 3.7</td><td>116.7 ± 2.2</td></tr></table>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Pressure calculated using the finite strain equations ($P_{fs}$) as a function of $P_{Fa}$ obtained from unit-cell volumes of NaCl (using EoS by Brown, 1999) and corrected for differential stress in both Fa and NaCl. The line shows the 1:1 ratio (not a fit). Grey crosses show $P_{NaCl}$ not corrected for stress.</figcaption></figure>

mulations give similar results and that $K_{02}$ values for $K_S$, $dK_S/dP$, G and $dG/dP$ are nearly identical to those obtained when fitting the data using $P_{Fa}$. On Fig. 7 we show the fit of $K_{02}$ data (Table 4) as a function of T. Points at the highest temperatures may seem to diverge from the trends, but removing them does not significantly change the results of the weighted least-square fit. The value of $K_S$ we calculated at 298 K is very close to the one found for $K_T$. The relation $K_S/K_T = 1 + \gamma\alpha T$, with $\gamma$ the Grüneisen parameter, $\alpha$ the thermal expansion and T the temperature, implies a $K_T$ about 1% smaller than $K_S$, a difference too small to resolve here because of our uncertainties.

# 4 Discussion

## 4.1 The hard case of fayalite

### 4.1.1 Oxygen fugacity

Physical properties of Fe-rich minerals are always a challenge to study, in particular at high pressures. For example, variations in the $Fe^{2+}/Fe^{3+}$ ratio may cause volume changes that in turn may perturb measurements. A non-exhaustive search of experimentally measured unit-cell volumes ($V_0$) of Fa at ambient conditions shows values between 306.9 Å$^3$ (natural sample, Richard and Richet, 1990) and 310.1 Å$^3$ (synthetic crystal, Graham et al., 1988), a range more than twice as large as that of forsterite, the pure Mg endmember of olivines. In a 1 atm furnace, performing annealing runs at constant oxygen fugacity ($fO_2$), and therefore at a constant $Fe^{2+}/Fe^{3+}$ ratio, is feasible with proper gas mixtures. At high pressures, solid buffers are used to maintain the $fO_2$ along a chemical equilibrium, often metal/metal oxide, and consequently the $fO_2$ changes with temperature and pressure. It is therefore possible that the range of $V_0$ mentioned above reflects the formation conditions and/or the thermal history of the different Fa specimens.

Our volume $V_0$ fitted from high-temperature data down to room temperature (Table 2) is very close to that of the starting material (Béjina et al., 2019), but, nevertheless, our high-temperature volumes are larger than what is expected from thermal expansion expressions (e.g., Fei, 1995; Kroll et al., 2012). There are very few studies on point defects in minerals in relation to their unit-cell volumes and elastic properties. To our knowledge, among minerals present in planetary mantles, only defects in wüstite, $Fe_xO$, have been studied in detail (McCammon and Liu, 1984; Zhang, 2000; Zhang and Zhao, 2005). These show that the cell parameter of $Fe_xO$ increases as the mineral tends toward perfect stoichiometry, i.e., as $Fe^{2+}$ stoichiometry tends to 1. However, a compar-

<footnote>a implied value; b fixed value.</footnote>
<table><caption>**Table 4.** Adiabatic bulk modulus, shear modulus and their pressure derivatives from fitting finite strain equations at each temperature using the full form of Davies and Dziewoński (1975), DD75, or their truncated version from Kung et al. (2002), K02, or the stress-corrected NaCl pressure scale, $P_{Fa}$.</caption><thead><tr><th>Temperature (K)</th><th>K<sub>S</sub> (GPa)</th><th>dK<sub>S</sub>/dP</th><th>d<sup>2</sup>K<sub>S</sub>/dP<sup>2</sup></th><th>G (GPa)</th><th>dG/dP</th><th>d<sup>2</sup>G/dP<sup>2</sup></th></tr></thead><tbody><tr><td colspan="7">DD75</td></tr><tr><td>302</td><td>125.7 ± 2.8</td><td>7.8 ± 2.3</td><td>-1.2 ± 0.5</td><td>52.4 ± 1.7</td><td>2.2 ± 1.4</td><td>-0.4 ± 0.5</td></tr><tr><td>373</td><td>124.0 ± 4</td><td>7.7 ± 2.8</td><td>-1.2 ± 0.6</td><td>51.5 ± 2.2</td><td>2.1 ± 1.7</td><td>-0.4 ± 0.6</td></tr><tr><td>472</td><td>123 ± 4</td><td>6.6 ± 2.9</td><td>-0.6 ± 0.6</td><td>50.9 ± 2.5</td><td>1.7 ± 1.7</td><td>-0.2 ± 0.5</td></tr><tr><td>572</td><td>117 ± 5</td><td>9.0 ± 3.3</td><td>-1.6 ± 0.7</td><td>50.0 ± 3.1</td><td>1.8 ± 2.0</td><td>-0.3 ± 0.6</td></tr><tr><td>673</td><td>114 ± 10</td><td>9 ± 5</td><td>-1.7 ± 1.1</td><td>47 ± 6</td><td>2.6 ± 3.2</td><td>-0.5 ± 0.9</td></tr><tr><td>773</td><td>105 ± 16</td><td>12 ± 8</td><td>-3.2 ± 2.2</td><td>47 ± 9</td><td>3 ± 5</td><td>-0.6 ± 1.2</td></tr><tr><td>874</td><td>114 ± 17</td><td>6 ± 8</td><td>-0.8 ± 1.2</td><td>45 ± 10</td><td>3 ± 5</td><td>-0.5 ± 1.0</td></tr><tr><td colspan="7">K02</td></tr><tr><td>302</td><td>126.9 ± 1.4</td><td>6.6 ± 0.5</td><td></td><td>53.5 ± 0.8</td><td>1.23 ± 0.27</td><td></td></tr><tr><td>373</td><td>125.6 ± 1.5</td><td>6.4 ± 0.4</td><td></td><td>52.6 ± 0.9</td><td>1.23 ± 0.27</td><td></td></tr><tr><td>472</td><td>123.5 ± 1.6</td><td>6.2 ± 0.4</td><td></td><td>51.7 ± 0.9</td><td>1.17 ± 0.27</td><td></td></tr><tr><td>572</td><td>121.4 ± 1.7</td><td>6.1 ± 0.5</td><td></td><td>51.0 ± 1.0</td><td>1.09 ± 0.27</td><td></td></tr><tr><td>673</td><td>118.5 ± 2.3</td><td>6.2 ± 0.6</td><td></td><td>50.0 ± 1.4</td><td>1.14 ± 0.33</td><td></td></tr><tr><td>773</td><td>116.9 ± 3.2</td><td>5.9 ± 0.7</td><td></td><td>49.6 ± 1.9</td><td>1.0 ± 0.4</td><td></td></tr><tr><td>874</td><td>114.6 ± 3.1</td><td>5.9 ± 0.6</td><td></td><td>49.3 ± 1.9</td><td>0.9 ± 0.4</td><td></td></tr><tr><td colspan="7">P<sub>Fa</sub></td></tr><tr><td>302</td><td>126.9 ± 0.7</td><td>6.4 ± 0.2</td><td></td><td>53.6 ± 0.3</td><td>1.15 ± 0.09</td><td></td></tr><tr><td>373</td><td>126.1 ± 1.1</td><td>6.1 ± 0.3</td><td></td><td>52.7 ± 0.3</td><td>1.18 ± 0.07</td><td></td></tr><tr><td>472</td><td>123.9 ± 0.7</td><td>6.0 ± 0.2</td><td></td><td>51.9 ± 0.2</td><td>1.06 ± 0.06</td><td></td></tr><tr><td>572</td><td>121.2 ± 0.8</td><td>6.0 ± 0.2</td><td></td><td>51.3 ± 0.2</td><td>0.97 ± 0.06</td><td></td></tr><tr><td>673</td><td>119.1 ± 0.9</td><td>5.7 ± 0.2</td><td></td><td>50.3 ± 0.4</td><td>0.99 ± 0.08</td><td></td></tr><tr><td>773</td><td>118.2 ± 1.5</td><td>5.4 ± 0.3</td><td></td><td>49.6 ± 0.4</td><td>0.98 ± 0.09</td><td></td></tr><tr><td>874</td><td>115.7 ± 1.6</td><td>5.3 ± 0.3</td><td></td><td>49.5 ± 0.5</td><td>0.84 ± 0.09</td><td></td></tr></tbody></table>

<table><caption>**Table 5.** Results from a linear least-square fit of data (Table 4) with *T* using the formulations of Davies and Dziewoński (1975), DD75, and Kung et al. (2002), K02, and compared with fitting the data using $P_{Fa}$. Uncertainties are $1\sigma$ calculated from the covariance of the minimization. *K* and *G* are expressed in gigapascals, *T* in kelvin, and $P_{Fa}$ in gigapascals.</caption><thead><tr><th></th><th>DD75</th><th>K02</th><th>P<sub>Fa</sub></th></tr></thead><tbody><tr><td>K<sub>S</sub>(0,298)</td><td>126.2 ± 2.4</td><td>127.16 ± 0.32</td><td>127.1 ± 0.6</td></tr><tr><td>dK<sub>S</sub>/dT (10<sup>-2</sup>)</td><td>-2.9 ± 0.5</td><td>-2.18 ± 0.05</td><td>-2.02 ± 0.10</td></tr><tr><td>dK<sub>S</sub>/dP</td><td>7.5 ± 1.7</td><td>6.52 ± 0.14</td><td>6.36 ± 0.11</td></tr><tr><td>d<sup>2</sup>K<sub>S</sub>/dPdT (10<sup>-3</sup>)</td><td>2.4 ± 3.1</td><td>-1.24 ± 0.22</td><td>-1.80 ± 0.17</td></tr><tr><td>d<sup>2</sup>K<sub>S</sub>/dP<sup>2</sup></td><td>-1.1 ± 0.7</td><td></td><td></td></tr><tr><td>d(d<sup>2</sup>K<sub>S</sub>/dP<sup>2</sup>)/dT (10<sup>-4</sup>)</td><td>-6.0 ± 13.0</td><td></td><td></td></tr><tr><td>G<sub>(0,298)</sub></td><td>52.5 ± 0.6</td><td>53.3 ± 0.4</td><td>53.03 ± 0.33</td></tr><tr><td>dG/dT (10<sup>-2</sup>)</td><td>-1.09 ± 0.11</td><td>-0.80 ± 0.07</td><td>-0.76 ± 0.05</td></tr><tr><td>dG/dP</td><td>2.0 ± 0.5</td><td>1.25 ± 0.05</td><td>1.16 ± 0.07</td></tr><tr><td>d<sup>2</sup>G/dPdT (10<sup>-4</sup>)</td><td>2.0 ± 8.0</td><td>-4.6 ± 0.7</td><td>-5.2 ± 1.1</td></tr><tr><td>d<sup>2</sup>G/dP<sup>2</sup></td><td>-0.35 ± 0.16</td><td></td><td></td></tr><tr><td>d(d<sup>2</sup>G/dP<sup>2</sup>)/dT (10<sup>-5</sup>)</td><td>-8.0 ± 27.0</td><td></td><td></td></tr></tbody></table>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** (a, b) Adiabatic bulk modulus and its pressure derivative as a function of temperature. (c, d) Same for shear modulus. Fits are performed with the formulation by Kung et al. (2002). Solid blue lines are the results of a weighted least-square fit to the data, and dashed lines show the 95% confidence interval. Results (labeled as K02) are listed in Table 5.</figcaption></figure>

ison between wüstite and fayalite should be made with caution because the point-defect populations are very different in the two minerals. FeO is known for its non-stoichiometry allowing for the accommodation of high levels of metal vacancies ($V_{Fe}$) which form clusters of $Fe^{3+}$ surrounded by $V_{Fe}$ (e.g., Minervini and Grimes, 1999), clusters that are not present in fayalite (e.g., Nakamura and Schmalzried, 1983). On the other hand, in a non-stoichiometric $Fe_{2}SiO_{4}$ wadsleyite, Hazen et al. (2000) proposed that the substitution of $2Fe^{3+}$ for $Fe^{2+} + Si$ increases its unit-cell volume. A parallel could be made with the defect $Fe_{Si}^{*}$ ($Fe^{3+}$ on a $Si^{4+}$ site) proposed in fayalite (Nakamura and Schmalzried, 1983; Hillert et al., 1996) that probably increases the size of the tetrahedron. $Fe^{3+}$-related defects could very well lead to larger unit cells in Fe-bearing silicates. Like for most Fe-rich minerals, it is difficult to prepare perfectly stoichiometric fayalite. In our high-pressure experiment, our sample was placed in a BN sleeve, a very reducing material, but also in contact with $Al_{2}O_{3}$ and NaCl, making it difficult to estimate the oxygen fugacity in the cell assembly and therefore the nature of the point-defect population. One must also note that our maximum temperature is, so far, the highest at which Fa unit-cell volumes have been measured at high pressure.

### 4.1.2 Shear modulus

Our value of the shear modulus at room conditions, $G_{(0.298)} = 53.3 \pm 0.4$ GPa, is close to previous data (Fig. 8). This is the first estimate of $G$ and $G'$ obtained from ultra-

sonic interferometry and on an aggregate since the study by Chung (1971), all other studies having used Fa single crystals (Graham et al., 1988; Speziale et al., 2004). Except for the $G'$ data point of Chung (1971) that is very low, all $G'$ are within the range [1.2, 1.8]. It is interesting to note that $G$ seems less sensitive to the sample's environment than $K$. For example, it is accepted that S waves are less sensitive to differential stress than P waves, especially if the wave propagation direction is aligned with the maximum principal stress. Keeping in mind that $G = \rho V_{s}^{2}$, it is not surprising that $G$ is less affected by differential stress.

### 4.1.3 Bulk modulus

Comparing bulk moduli between studies is always difficult and more so than for $G$. For Fe–Mg olivines, $G$ decreases linearly with $X_{Fe}$ by about 30 GPa over the entire solid solution, whereas $K$ varies by only a few gigapascals with no clear trend with $X_{Fe}$. As discussed, for example, in Angel (2000), the fitted values of $K$ and $K'$ are often strongly correlated, and therefore, fixing either one affects the estimation of the other. This can be seen in Fig. 8 representing experimental values of $K$ of fayalite plotted against $dK/dP$. The figure shows a clear trend with a negative slope, meaning that as $K'$ increases, $K$ decreases. Many high-pressure experiments, in particular in multi-anvil, only explore a limited pressure range. For such experiments, analyzing data often requires one to fix either $K$ or $K'$ and calculate the other from a fit to the data. In our case, for example, Table 2 shows that for
<figure><img src="image_11.png" /></figure>

**Figure 8.** Trend of two elastic parameters for fayalite with their respective pressure derivatives. (a) Bulk modulus versus d*K*/dP. The blue color is for *K*<sub>S</sub> values, grey for *K*<sub>T</sub> and red for this study. *K*<sub>S</sub> and *K*<sub>T</sub> values are plotted as published, with no conversion from one to the other. The full red square is our data as reported in Table 4 (*P*<sub>Fa</sub>), the open square being *K*<sub>S</sub> fitted to our data with fixed *K*'<sub>S</sub> = 5.3. Arrows at the end of error bars of the data of Yagi et al. (1975) symbolize that they extend outside the plot range. (b) Shear modulus versus d*G*/dP. On both plots, a given symbol refers to the same publication. Dotted lines represent values of *K*<sub>S</sub> and *G* determined at room pressure; therefore, no pressure derivative was determined.

$K'_T$ fixed at 4 (second-order BM EoS) and 5.3 (Nestola et al., 2011), $K_T$ varies by almost 3 GPa.

Another source of difficulty when comparing EoS param-
eters comes from the choice of a pressure scale. Histori-
cally the NaCl Decker scale (Decker, 1971) has been the
most used one, but Strässle et al. (2014) proposed that the
scale of Brown (1999) is more accurate. When measuring
sound wave velocities, the finite strain approach (Davies and
Dziewoński, 1975) allows for the determination of the elastic
parameters without prior knowledge of the pressure, which
can be calculated afterwards (see their Eq. 1c). These differ-
ent methods can lead to small changes in pressure estimates
with a strong impact on the determination of an EoS, in par-
ticular for d*K*/dP.

Among other factors affecting the determination of *K* are,
as mentioned above, the sample's stoichiometry and there-
fore the *f*O<sub>2</sub> conditions of the experiment, as well as the
(often ignored) level of non-hydrostaticity of most high-
pressure experiments.

Our result for *K*<sub>*T*</sub> measured at room temperature,
*K*<sub>*T*</sub>(0.302) = 130±1 GPa for *K*'<sub>*T*</sub> = 5.3 (Table 3), is within the
uncertainties of our previous estimate (Béjina et al., 2019).

Our choice of $K'_T$ comes from the high accuracy of the data
of Nestola et al. (2011) measured for olivines with $X_{Fe}$ be-
tween 8% and 38%, close to estimates, $K'_S = 4.9$ (Speziale
et al., 2004) and 5.2 (Graham et al., 1988), for fayalite. When
recalculating $K_T$ using EoS Fit 7 and high-temperature data
as well, we obtain $K_{T(0.298)} = 129 \pm 5$ GPa for the same
$K'_T = 5.3$ (Table 2). From the relation between $K_S$ and $K_T$
given in Sect. 3.4, one would expect $K_S \approx 131$ GPa with
$K'_S = 5.3$ at 300 K. The values we deduced from our dataset
are $K_S = 127.2 \pm 0.3$ GPa and $K'_S = 6.5 \pm 0.1$ (Table 5). Fix-
ing $K'_S$ at 5.3 would increase $K_S$ to about 130 GPa, a value
very consistent with our estimate of $K_T$ (Table 2). These val-
ues follow the trend $K_S = f(dK_S/dP)$ from previous stud-
ies (Fig. 8) and therefore are in general agreement with these.
From this figure, $K_S$ and $dK_S/dP$ of fayalite seem to fall in
a limited range, [125, 135] GPa and [4, 6.5] respectively.

It is in general difficult to pinpoint a main reason for low
values of K_S even when correlated with high K'_S. Among
a variety of reasons, we possibly have the wrong tempera-
ture reading of the thermocouple, the effect of oxygen fu-
gacity and hence point defects on the elasticity of the olivine
structure, the effect of differential stress (Yagi et al., 1975),
<figure><img src="image_12.png" /><figcaption>**Figure 9.** P wave (top) and S wave (bottom) velocities in fayalite calculated with *burnman* (Cottaar et al., 2014) using the formalism of EoS of Stixrude and Lithgow-Bertelloni (2011) (SLB2011). Full lines are Mars profiles along the hot areotherm A and dashed lines along the cold areotherm C of Bertka and Fei (1997). Dotted lines are along a geotherm (Priestley and McKenzie, 2013) included in *burnman*. Pressure and temperature profiles are shown in the inset. Red lines are velocities calculated using the elastic parameters from this work and grey curves those for parameters from Stixrude and Lithgow-Bertelloni (2011). S wave velocities in ahrensite are also shown to visualize the velocity jump when transiting to the spinel structure.</figcaption></figure>

compositional inhomogeneities or microstructure irregularities, etc. In our case, one could argue that the transition of our Fa specimen to ahrensite (Ahr) may have impacted our estimates of the elastic parameters. However we emphasize that no Ahr was observed by SEM, no distinct Ahr peaks can be seen on the XRD patterns, and sound velocities do not increase when approaching the phase boundary as seen by Liu et al. (2010).

In their study of nanocrystalline forsterite (11 nm grain size prior to the experiment), Couvy et al. (2010) proposed that a small grain size (i.e., many grain boundaries and triple junctions) leads to a lower $K_T$. This effect has been observed in different materials and for grain sizes usually below 100 nm (e.g., see Marquardt et al., 2011, for MgO and references therein for other materials). To our knowledge no systematic studies have been performed, in particular on olivines, to quantify this effect of grain size on the bulk modulus. The grain size of our Fa sample, prepared from nano-size oxide powders, was about 1 µm prior to the experiment and around 3 µm after. Even though much larger than 100 nm, we cannot rule out that the grain size of our sample could be small enough to start affecting its elastic properties. This assumption is difficult to verify because postmortem microstructures are rarely described in most studies.

## 4.2 Seismic velocities in Mars' upper mantle

In order to visualize the effect of fayalite elastic parameters on P and S wave velocities, we calculated $V_p$ and $V_s$ using the Python package *burnman* (Cottaar et al., 2014) with the formalism of Stixrude and Lithgow-Bertelloni (2011) (SLB2011) and data for fayalite ($K_{T0}$ = 135 GPa, $K' = 4.22$, $G = 50.9$ GPa, $G' = 1.46$) and then with our determination of these parameters (Table 5). We kept the temperature-dependent parameters of SLB2011 because of the large uncertainties in ours. Figure 9 represents the P and S velocities in Fa as a function of pressure along two areotherms for Mars (profiles A and C of Bertka and Fei, 1997) and a geotherm from the CAM2016 model (Priestley and McKenzie, 2013). We chose this geotherm in order to have a very different temperature profile compared to the two areotherms. This figure shows that the differences in P and S velocities with the change in parameters are very small. For the geotherm the effect of the hot temperatures prevails over that of pressure, and we observe a decrease in the velocities up to about 5 GPa.

To further test the influence of elastic parameters of Fa on its P and S velocities, we performed the same calculations along the areotherm A from Bertka and Fei (1997), allowing all parameters from the SLB2011 EoS to vary in the inter-
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Velocity profiles in Mars' upper mantle calculated with burnman using the formulation of Stixrude and Lithgow-Bertelloni (2011). The black dots represent the velocities for fayalite using our *Ks*, *dKs/dP*, *G* and *dG/dP*. The red, green and blue domains are the ranges of velocities calculated when varying thermoelastic parameters, as described in the text.</figcaption></figure>

vals: $K_{T0} = [125, 135]$ GPa, $G_0 = [50, 53]$ GPa, $K' = [4, 6]$, $G' = [1, 1.3]$, the Debye temperature $\Theta = [550, 650]$ K, the Grüneisen parameter $\gamma = [1, 1.3]$ and exponent $q = [3, 4]$. Figure 10 shows $V_p$, $V_s$ and $V_p/V_s$ along this areotherm for Fa and two olivines: one with 25 mol% Fe and one with 40 mol% Fe. The EoS for these olivines are calculated by mixing Fe and Mg olivine endmembers according to their molar fraction (Cottaar et al., 2014). In our calculation, only Fa parameters vary, those of forsterite (Fo) being considered as perfectly determined. Consequently mixing Fo and Fa parameters to calculate velocities results in smaller domains of possible velocities compared to pure Fo, but the difference in velocities between these two olivines is rather small. This could be a challenge when fitting data from the SEIS instrument.

Another way to look at this is to test the effect of the bulk modulus on the optimized olivine content resulting from a fit of a simple model of Earth's upper mantle to the Preliminary Reference Earth Model (PREM) (Davies and Dziewoński, 1975) using burnman. We only test here the sensitivity with K because its experimental determination presents a much larger dispersion than G. We built a very simple model of the Earth's upper mantle with three mineral phases: olivine, orthopyroxene (Opx) and garnet (Gt), all three being simple Fe–Mg solid solutions. The molar percentage ratio of Opx to Gt was fixed at 2:1, and only the olivine fraction

<figure><img src="image_14.png" /><figcaption>**Figure 11.** Optimum olivine contents determined by adjusting $V_p$ of a simple upper-mantle model to PREM using burnman for three $K_T$ values for fayalite. Large variations in $K_T$ produce small changes in olivine fractions. See text for calculation details.</figcaption></figure>

was optimized. These calculations were done for three values of the olivine $K_{T0}$, 120, 130 and 140 GPa, which correspond approximately to the range of experimental values determined for this mineral. The properties of this modeled rock were calculated as described by Cottaar et al. (2014)
along a geotherm (Anderson, 1982, provided by burnman). Figure 11 shows that the three optimum olivine contents fall in a tight range despite the large variation in $K_{T0}$, between 61 mol % and 66 mol %. Therefore, considering the uncertainties in the effect of Fe on the bulk modulus of olivine, it is difficult to say at this point if the SEIS data will have the necessary resolution to determine the Fe concentration of the minerals of Mars' upper mantle with an uncertainty less than 15 %.

## 5 Conclusions

The bulk and shear moduli of fayalite presented here are the first to be determined by coupled synchrotron X-ray radiation and ultrasonic interferometry. A fit using the third-order Birch–Murnaghan EoS with $K'_T = 5.3$ (fixed) to the data gives $K_T = 128.6 \pm 5.1$ GPa at room conditions, whereas $K_S = 127.2 \pm 0.3$ GPa and $K'_S = 6.5$ were obtained from sound wave velocities. Fixing $K'_S$ to 5.3 increases $K_S$ to 130 GPa, therefore showing that our $K_T$ and $K_S$ are compatible within uncertainties. These values of $K$ and $K'$ follow the linear trend $K = f(K')$ defined by previously published data. The shear modulus values of fayalite show less dispersion than the bulk modulus; therefore, the latter will be critical when fitting a mineralogical model to Mars InSight SEIS data. The range of $K$ values for the olivine Fe–Mg solid solution is defined between about 125 and 135 GPa. With such variation, a precise compositional and mineralogical model of the Mars mantle will be challenging to obtain. Because experiments on Fe-rich minerals are always difficult, more experiments are needed to better define the trend of $K$ with $X_{Fe}$ so that more robust values of $K$ can be established for minerals with an Fe content consistent with Mars' mantle composition.

*Data availability.* Data are provided in the main text and in the Supplement.

*Supplement.* The supplement related to this article is available online at: https://doi.org/10.5194/ejm-33-519-2021-supplement.

*Author contributions.* FB and MB designed the experiment. FB analyzed the data and wrote the manuscript. MB analyzed ultrasonic data and corrected the manuscript. All authors participated in running the experiment.

*Competing interests.* The authors declare that they have no conflict of interest.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

*Special issue statement.* This article is part of the special issue "Probing the Earth: experiments and mineral physics at mantle depths". It is a result of the 17th International Symposium on Experimental Mineralogy, Petrology and Geochemistry, Potsdam, Germany, 21–24 June 2020.

*Acknowledgements.* The high-pressure experiments were performed at hutch X17B2 of NSLS (Brookhaven National Lab.) thanks to the support of COMPRES. We would also like to thank the two anonymous reviewers for their constructive remarks.

*Financial support.* This research has been supported by the TS8 CNES, focusing on Mars InSight SEIS, and by ANR MAGIS (grant no. ANR-19-CE31-0008-08).

*Review statement.* This paper was edited by Monika Koch-Müller and reviewed by two anonymous referees.

## References

Anderson, O. L.: The Earth's core and the phase diagram of iron, Philos. T. R. Soc. A, 306, 21–35, https://doi.org/10.1098/rsta.1982.0063, 1982.

Andrault, D., Bouhifd, M., Itié, J., and Richet, P.: TS3 Compression and amorphization of (Mg,Fe)2SiO4 olivines: An X-ray diffraction study up to 70 GPa, Phys. Chem. Miner., 22, 99–107, https://doi.org/10.1007/BF00202469, 1995.

Angel, R. J.: Equations of State, Rev. Mineral. Geochem., 41, 35–59, https://doi.org/10.2138/rmg.2000.41.2, 2000.

Angel, R. J., Alvaro, M., and González-Platas, J.: EosFit7c and a Fortran module (library) for equation of state calculations, Z. Kristallogr., 229, 405–419, https://doi.org/10.1515/zkri-2013-1711, 2014.

Béjina, F. and Bystricky, M.: Use of the spark plasma sintering technique for the synthesis of dense mineral aggregates suitable for high-pressure experiments, High Pressure Res., 29, 630–634, https://doi.org/10.1080/08957950903425249, 2009.

Béjina, F., Bystricky, M., Tercé, N., Whitaker, M. L., and Chen, H.: Bulk modulus of Fe-rich olivines corrected for non-hydrostaticity, C. R. Geosci., 351, 86–94, https://doi.org/10.1016/j.crte.2018.06.002, 2019.

Bertka, C. M. and Fei, Y.: Mineralogy of the Martian interior up to core-mantle boundary pressures, J. Geophys. Res., 102, 5251–5264, https://doi.org/10.1029/96JB03270, 1997.

Birch, F.: Finite elastic Strain of cubic Crystals, Phys. Rev., 71, 809–824, https://doi.org/10.1103/PhysRev.71.809, 1947.

Brown, J. M.: The NaCl pressure standard, J. Appl. Phys., 86, 5801–5808, https://doi.org/10.1063/1.371596, 1999.
Chen, J., Li, L., Yu, T., Long, H., Weidner, D., Wang, L., and Vaughan, M.: Do Reuss and Voigt bounds really bound in high-pressure rheology experiments?, J. Phys., 18, S1049–S1059, https://doi.org/10.1088/0953-8984/18/25/S11, 2006.

Chung, D. H.: Elasticity and equations of state of olivines in the Mg₂SiO₄-Fe₂SiO₄ system, Geophys. J. Int., 25, 511–538, https://doi.org/10.1111/j.1365-246x.1971.tb02201.x, 1971.

Cottaar, S., Heister, T., Rose, I., and Unterborn, C.: BurnMan: A lower mantle mineral physics toolkit, Geochem. Geophys. Geosy., 15, 1164–1179, https://doi.org/10.1002/2013GC005122, 2014.

Couvy, H., Chen, J., and Drozd, V.: Compressibility of nanocrystalline forsterite, Phys. Chem. Miner., 37, 343–351, https://doi.org/10.1007/s00269-009-0337-8, 2010.

Davies, G. F. and Dziewoński, A. M.: Homogeneity and constitution of the Earth's lower mantle and outer core, Phys. Earth Planet. Int., 10, 336–343, https://doi.org/10.1016/0031-9201(75)90060-6, 1975.

Decker, D. L.: High-Pressure Equation of State for NaCl, KCl, and CsCl, J. Appl. Phys., 42, 3239–3244, https://doi.org/10.1063/1.1660714, 1971.

Dreibus, G. and Wänke, H.: Mars, a volatile-rich planet, Meteoritics, 20, 367–381, 1985.

Fei, Y.: Thermal Expansion, in: Mineral Physics & Crystallography: A Handbook of Physical Constants, edited by: Ahrens, T. J., AGU, Washington, DC, 29–44, https://doi.org/10.1029/RF002p0029, 1995.

Graham, E. K., Schwab, J. A., Sopkin, S. M., and Takei, H.: The pressure and temperature dependence of the elastic properties of single-crystal Fayalite Fe₂SiO₄, Phys. Chem. Miner., 16, 186–198, https://doi.org/10.1007/BF00203203, 1988.

Guignard, J., Bystricky, M., and Béjina, F.: Dense fine-grained aggregates prepared by spark plasma sintering (SPS), an original technique in experimental petrology, Eur. J. Mineral., 23, 323–331, https://doi.org/10.1127/0935-1221/2011/0023-2099, 2011.

Hazen, R., Yang, H., and Prewitt, C.: High-pressure crystal chemistry of Fe³⁺-wadsleyite, β-Fe₂.₃₃Si₀.₆₇O₄, Am. Mineral., 85, 778–783, https://doi.org/10.2138/am-2000-5-618, 2000.

Hearmon, R. F. S.: The elastic constants of anisotropic materials – II, Adv. Phys., 5, 323–382, https://doi.org/10.1080/00018732.1956.tADP0323, 1956.

Hillert, M., Selleby, M., and Sundman, B.: A Reassessment of the Non-Stoichiometry of Fayalite, Phys. Chem. Miner., 23, 387–390, https://doi.org/10.1007/BF00199504, 1996.

Holland, T. J. B. and Redfern, S. A. T.: Unit cell refinement from powder diffraction data: the use of regression diagnostics, Mineral. Mag., 61, 65–77, https://doi.org/10.1107/S0021889896011673, 1997.

Isaak, D. G., Graham, E. K., Bass, J. D., and Wang, H.: The elastic properties of single-crystal fayalite as determined by dynamical measurement techniques, Pure Appl. Geophys., 141, 393–414, https://doi.org/10.1007/BF00998337, 1993.

Krischer, L., Megies, T., Barsch, R., Beyreuther, M., Lecocq, T., Caudron, C., and Wassermann, J.: ObsPy: a bridge for seismology into the scientific Python ecosystem, Comput. Sci. Discov., 8, 014003, https://doi.org/10.1088/1749-4699/8/1/014003, 2015.

Kroll, H., Kirfel, A., Heinemann, R., and Barbier, B.: Volume thermal expansion and related thermophysical parameters in the Mg,Fe olivine solid-solution series, Eur. J. Mineral., 24, 935–956, https://doi.org/10.1127/0935-1221/2012/0024-2235, 2012.

Kudoh, Y. and Takeda, H.: Single crystal X-ray diffraction study on the bond compressibility of Fayalite, Fe₂SiO₄ and Rutile, TiO₂ under high pressure, Physica B+C, 139–140, 333–336, https://doi.org/10.1016/0378-4363(86)90591-7, 1986.

Kung, J., Li, B., Weidner, D., Zhang, J., and Liebermann, R.: Elasticity of (Mg₀.₈₃,Fe₀.₁₇)O ferropericlase at high pressure: ultrasonic measurements in conjunction with X-radiation techniques, Earth Planet. Sc. Lett., 203, 557–566, https://doi.org/10.1016/S0012-821X(02)00838-5, 2002.

Li, B. and Zhang, J.: Pressure and temperature dependence of elastic wave velocity of MgSiO₃ perovskite and the composition of the lower mantle, Phys. Earth Planet. Int., 151, 143–154, https://doi.org/10.1016/j.pepi.2005.02.004, 2005.

Li, B., Kung, J., and Liebermann, R. C.: Modern techniques in measuring elasticity of Earth materials at high pressure and high temperature using ultrasonic interferometry in conjunction with synchrotron X-radiation in multi-anvil apparatus, Phys. Earth Planet. Int., 143-144, 559–574, https://doi.org/10.1016/j.pepi.2003.09.020, 2004.

Liu, Q., Liu, W., Whitaker, M. L., Wang, L., and Li, B.: In situ ultrasonic velocity measurements across the olivine-spinel transformation in Fe₂SiO₄, Am. Mineral., 95, 1000–1005, https://doi.org/10.2138/am.2010.3369, 2010.

Liu, W., Kung, J., and Li, B.: Elasticity of San Carlos olivine to 8 GPa and 1073 K, Geophys. Res. Lett., 32, L16301, https://doi.org/10.1029/2005GL023453, 2005.

Lognonne, P., Banerdt, W. B., Giardini, D., Pike, W. T., TS10 Christensen, U., Laudet, P., de Raucourt, S., Zweifel, P., Calcutt, S., Bierwirth, M., Hurst, K. J., Ijpelaan, F., Umland, J. W., Llorca-Cejudo, R., Larson, S. A., Garcia, R. F., Kedar, S., Knapmeyer-Endrun, B., Mimoun, D., Mocquet, A., Panning, M. P., Weber, R. C., Sylvestre-Baron, A., Pont, G., Verdier, N., Kerjean, L., Facto, L. J., Gharakanian, V., Feldman, J. E., Hoffman, T. L., Klein, D. B., Klein, K., Onufer, N. P., Paredes-Garcia, J., Petkov, M. P., Willis, J. R., Smrekar, S. E., Drilleau, M., Gabsi, T., Nebut, T., Robert, O., Tillier, S., Moreau, C., Parise, M., Aveni, G., Ben Charef, S., Bennour, Y., Camus, T., Dandonneau, P. A., Desfoux, C., Lecomte, B., Pot, O., Revuz, P., Mance, D., tenPierick, J., Bowles, N. E., Charalambous, C., Delahunty, A. K., Hurley, J., Irshad, R., Liu, Huafeng, Mukherjee, A. G., Standley, I. M., Stott, A. E., Temple, J., Warren, T., Eberhardt, M., Kramer, A., Kühne, W., Miettinen, E. P., Monecke, M., Aicardi, C., André, M., Baroukh, J., Borrien, A., Bouisset, A., Boutte, P., Brethomé, K., Brysbaert, C., Carlier, T., Deleuze, M., Desmarres, J. M., Dilhan, D., Doucet, C., Faye, D., Faye-Refalo, N., Gonzalez, R., Imbert, C., Larigauderie, C., Locatelli, E., Luno, L., Meyer, J. R., Mialhe, F., Mouret, J. M., Nonon, M., Pahn, Y., Paillet, A., Pasquier, P., Perez, G., Pérez, R., Perrin, L., Pouilloux, B., Rosak, A., de Larclause, I. Savin, Sicre, J., Sodki, M., Toulemont, N., Vella, B., Yana, C., Alibay, F., Avalos, O. M., Balzer, M. A., Bhandari, P., Blanco, E., Bone, B. D., Bousman, J. C., Bruneau, P., Calef, F. J., Calvet, R. J., D'Agostino, S. A., de los Santos, G., Deen, R. G., Denise, R. W., Ervin, J., Ferraro, N. W., Gengl, H. E., Grinblat, F., Hernandez, D., Hetzel, M., Johnson, M. E., Khachikyan, L., Lin, J. Y., Madzunkov, S. M., Marshall, S. L., Mikellides, I. G., Miller, E. A., Raff, W., Singer, J. E., Sunday, C. M., Villalvazo, J. F., Wallace, M. C., Banfield, D., Rodriguez-Manfredi, J. A., Russell, C. T., Trebi-Ollennu, A., Maki, J. N., Beucler, E., Böse, M., Bonjour, C., Berenguer, J.
L., Ceylan, S., Clinton, J., Conejero, V., Daubar, I., Dehant, V., Delage, P., Euchner, F., Estève, I., Fayon, L., Ferraioli, L., Johnson, C. L., Gagnepain-Beyneix, J., Golombek, M., Khan, A., Kawamura, T., Kenda, B., Labrot, P., Murdoch, N., Pardo, C., Perrin, C., Pou, L., Sauron, A., Savoie, D., Stähler, S., Stutzmann, E., Teanby, N. A., Tromp, J., van Driel, M., Wieczorek, M., Widmer-Schnidrig, R., and Wookey, J.: SEIS: Insight's Seismic Experiment for Internal Structure of Mars, Space Sci. Rev., 215, 12, https://doi.org/10.1007/s11214-018-0574-6, 2019.

Marquardt, H., Gleason, A., Marquardt, K., Speziale, S., Miyagi, L., Neusser, G., Wenk, H.-R., and Jeanloz, R.: Elastic properties of MgO nanocrystals and grain boundaries at high pressures by Brillouin scattering, Phys. Rev. B, 84, 064131, https://doi.org/10.1103/PhysRevB.84.064131, 2011.

McCammon, C. A. and Liu, L.-G.: The effects of pressure and temperature on nonstoichiometric wüstite, FeₓO: The iron-rich phase boundary, Phys. Chem. Miner., 10, 106–113, https://doi.org/10.1007/BF00309644, 1984.

Minervini, L. and Grimes, R. W.: Defect clustering in wüstite, J. Phys. Chem. Solids, 60, 235–245, https://doi.org/10.1016/S0022-3697(98)00273-X, 1999.

Nakamura, A. and Schmalzried, H.: On the Nonstoichiometry and Point Defects of Olivine, Phys. Chem. Miner., 10, 27–37, https://doi.org/10.1007/BF01204323, 1983.

Nestola, F., Pasqual, D., Smyth, J. R., Novella, D., Secco, L., Manghnani, M. H., and Dal Negro, A.: New accurate elastic parameters for the forsterite-fayalite solid solution, Am. Mineral., 96, 1742–1747, https://doi.org/10.2138/am.2011.3829, 2011.

Newville, M., Stensitzki, T., Allen, D. B., and Ingargiola, A.: LMFIT: Non-Linear Least-Square Minimization and Curve-Fitting for Python, Zenodo [data set], https://doi.org/10.5281/zenodo.4516651, 2014.

Pamato, M. G., Nestola, F., Novella, D., Smyth, J. R., Pasqual, D., Gatta, G. D., Alvaro, M., and Secco, L.: The High-Pressure Structural Evolution of Olivine along the Forsterite-Fayalite Join, Minerals, 9, 790, https://doi.org/10.3390/min9120790, 2019.

Plymate, T. and Stout, J.: Pressure-volume-temperature behavior of fayalite based on static compression measurements at 400°C, Phys. Chem. Miner., 17, 413–420 https://doi.org/10.1007/BF00201452, 1990.

Priestley, K. and McKenzie, D.: The relationship between shear wave velocity, temperature, attenuation and viscosity in the shallow part of the mantle, Earth Planet. Sc. Lett., 381, 78–91, https://doi.org/10.1016/j.epsl.2013.08.022, 2013.

Richard, G. and Richet, P.: Room-temperature amorphization of fayalite and high-pressure properties of Fe₂SiO₄ liquid, Geophys. Res. Lett., 17, 2093–2096, https://doi.org/10.1029/GL017i012p02093, 1990.

Singh, A., Balasingh, C., Mao, H.-K., Hemley, R., and Shu, J.: Analysis of lattice strains measured under nonhydrostatic pressure, J. Appl. Phys., 83, 7567–7575, https://doi.org/10.1063/1.3678782, 1998.

Smyth, J. R.: High temperature crystal chemistry of Fayalite, Am. Mineral., 60, 1092–1097, 1975.

Spetzler, H., Sammis, C. G., and O'Connell, R. J.: Equation of state of NaCl: Ultrasonic measurements to 8 kbar and 800 °C and static lattice theory, J. Phys. Chem. Solids, 33, 1727–1750, https://doi.org/10.1016/S0022-3697(72)80468-2, 1972.

Speziale, S., Duffy, T. S., and Angel, R. J.: Single-crystal elasticity of fayalite to 12 GPa, J. Geophys. Res., 109, B12202, https://doi.org/10.1029/2004JB003162, 2004.

Stackhouse, S., Stixrude, L., and Karki, B.: Determination of the high-pressure properties of fayalite from first-principles calculations, Earth Planet. Sc. Lett., 289, 449–456, https://doi.org/10.1016/j.epsl.2009.11.033, 2010.

Stähler, S. C., Khan, A., Banerdt, W. B., Lognonné, P., Giardini, D., and Ceylan, S., Drilleau, M., Duran, A. C., Garcia, R. F., Huang, Q., Kim, D., Lekic, V., Samuel, H., Schimmel, M., Schmerr, N., Sollberger, D., Stutzmann, É., Xu, Z., Antonangeli, D., Charalambous, C., Davis, P. M., Irving, J. C. E., Kawamura, T., Knapmeyer, M., Maguire, R., Marusiak, A. G., Panning, M. P., Perrin, C., Plesa, A.-C., Rivoldini, A., Schmelzbach, C., Zenhäusern, G., Beucler, É., Clinton, J., Dahmen, N., van Driel, M., Gudkova, T., Horleston, A., Pike, W. T., Plasman, M., and Smrekar, S. E.: Seismic Detection of the Martian Core, Science, 373, 443–448, https://doi.org/10.1126/science.abi7730, 2021.

Stixrude, L. and Lithgow-Bertelloni, C.: Thermodynamics of mantle minerals - I. Physical properties, Geophys. J. Int., 162, 610–632, https://doi.org/10.1111/j.1365-246X.2005.02642.x, 2005.

Stixrude, L. and Lithgow-Bertelloni, C.: Thermodynamics of mantle minerals – II. Phase equilibria, Geophys. J. Int., 184, 1180–1213, https://doi.org/10.1111/j.1365-246X.2010.04890.x, 2011.

Strässle, T., Klotz, S., Kunc, K., Pomjakushin, V., and White, J. S.: Equation of state of lead from high-pressure neutron diffraction up to 8.9 GPa and its implication for the NaCl pressure scale, Phys. Rev. B, 90, 014101, https://doi.org/10.1103/PhysRevB.90.014101, 2014.

Sumino, Y.: The elastic constants of Mn₂SiO₄, Fe₂SiO₄ and Co₂SiO₄, and the elastic properties of olivine group minerals at high temperature, J. Phys. Earth, 27, 209–238, https://doi.org/10.4294/jpe1952.27.209, 1979.

Taylor, G. J.: The bulk composition of Mars, Chem. Erde, 73, 401–420, https://doi.org/10.1016/j.chemer.2013.09.006, 2013.

Weidner, D. J., Vaughan, M. T., Wang, L., Long, H., Li, L., Dixon, N. A., and Durham, W. B.: Precise stress measurements with white synchrotron x rays, Rev. Sci. Instrum., 81, 013903, https://doi.org/10.1063/1.3263760, 2010.

Whitaker, M. L., Baldwin, K. J., and Huebsch, W. R.: DIAS-CoPE: Directly integrated acoustic system combined with pressure experiments – A new method for fast acoustic velocity measurements at high pressure, Rev. Sci. Instrum., 88, 034901, https://doi.org/10.1063/1.4977596, 2017.

Williams, Q., Knittle, E., Reichlin, R., Martin, S., and Jeanloz, R.: Structural and electronic properties of Fe₂SiO₄-Fayalite at ultrahigh pressures: Amorphization and gap closure, J. Geophys. Res., 95, 21549–21563, https://doi.org/10.1029/JB095iB13p21549, 1990.

Witze, A.: Seismic ripples reveal size of Mars's core, Nature, 591, 514–515, https://doi.org/10.1038/d41586-021-00696-7, 2021.

Wojdyr, M.: Fityk: a general-purpose peak fitting program, J. Appl. Crystallogr., 43, 1126–1128, https://doi.org/10.1107/S0021889810030499, 2010.

Yagi, T., Ida, Y., Sato, Y., and Akimoto, S.-I.: Effect of hydrostatic pressure on the lattice parameters of Fe₂SiO₄ olivine up to 70 kbar, Phys. Earth Planet. Int., 10, 348–354, https://doi.org/10.1016/0031-9201(75)90062-X, 1975.
Zhang, J.: Effect of defects on the elastic proper-
ties of wüstite, Phys. Rev. Lett., 84, 507–510,
https://doi.org/10.1103/PhysRevLett.84.507, 2000.

Zhang, L.: Single crystal hydrostatic compression of (Mg, Mn, Fe, Co)2SiO4 olivines, Phys. Chem. Miner., 25, 308–312, https://doi.org/10.1007/s002690050119, 1998.

Zhang, J. and Zhao, Y.: Effects of defect and pressure on the ther-
mal expansivity of Fe<sub>x</sub>O, Phys. Chem. Miner., 32, 241–247,
https://doi.org/10.1007/s00269-005-0447-x, 2005.

Zhang, J. S., Hu, Y., Shelton, H., Kung, J., and Dera, P.: Single-crystal X-ray diffraction study of Fe2SiO4 Fayalite up to 31 GPa, Phys. Chem. Miner., 44, 171–179, https://doi.org/10.1007/s00269-016-0846-1, 2017.