Seismic Properties of a Unique Olivine-Rich Eclogite in the Western Gneiss Region, Norway

Investigating the seismic properties of natural eclogite is crucial for identifying the composition, density, and mechanical structure of the Earth’s deep crust and mantle. For this purpose, numerous studies have addressed the seismic properties of various types of eclogite, except for a rare eclogite type that contains abundant olivine and orthopyroxene. In this contribution, we calculated the ambient-condition seismic velocities and seismic anisotropies of this eclogite type using an olivine-rich eclogite from northwestern Flemsøya in the Nordøyane ultrahigh-pressure (UHP) domain of the Western Gneiss Region in Norway. Detailed analyses of the seismic properties data suggest that patterns of seismic anisotropy of the Flem eclogite were largely controlled by the strength of the crystal-preferred orientation (CPO) and characterized by significant destructive effects of the CPO interactions, which together, resulted in very weak bulk rock seismic anisotropies (AVp = 1.0–2.5%, max. AVs = 0.6–2.0%). The magnitudes of the seismic anisotropies of the Flem eclogite were similar to those of dry eclogite but much lower than those of gabbro, peridotite, hydrous-phase-bearing eclogite, and blueschist. Furthermore, we found that amphibole CPOs were the main contributors to the higher seismic anisotropies in some amphibole-rich samples. The average seismic velocities of Flem eclogite were greatly affected by the relative volume proportions of omphacite and amphibole. The Vp (8.00–8.33 km/s) and Vs (4.55–4.72 km/s) were remarkably larger than the hydrous-phase-bearing eclogite, blueschist, and gabbro, but lower than dry eclogite and peridotite. The Vp/Vs ratio was almost constant (avg. ≈ 1.765) among Flem eclogite, slightly larger than olivine-free dry eclogite, but similar to peridotite, indicating that an abundance of olivine is the source of their high Vp/Vs ratios. The Vp/Vs ratios of Flem eclogite were also higher than other (non-)retrograded eclogite and significantly lower than those of gabbro. The seismic features derived from the Flem eclogite can thus be used to distinguish olivine-rich eclogite from other common rock types (especially gabbro) in the deep continental crust or subduction channel when high-resolution seismic wave data are available.


Introduction
Eclogite is a unique high-pressure to ultrahigh-pressure (HP-UHP) metamorphic rock that commonly originates from the subduction of an oceanic or continental crust or derives from the thickening of an orogenic crust during a continental collision [1][2][3][4]. A typical eclogite consists of a mainly bi-mineralic assemblage (i.e., garnet and omphacite), which is metamorphosed from mafic components, such as gabbro, diabase, and basalt [4]. Hydrous minerals, such as amphibole (e.g., glaucophane and hornblende), epidote, lawsonite, and mica (e.g., phengite and biotite) are also common in the eclogite and indicate different equilibrated pressure-temperature (P-T) conditions and hydrous states of the eclogite facies' metamorphisms [5][6][7]. Magnesium-rich minerals, such as olivine and orthopyroxene, can also be sporadically abundant in the eclogite when its protolith is rich in olivine (i.e., olivine gabbro or troctolite) [8].
The eclogite in the Flem Gabbro from Flemsøya in the Western Gneiss Region (WGR), Norway, is exceptional due to its peculiar mineral assemblage of olivine and orthopyroxene, which were inherited from an olivine gabbro protolith [8,33]. This eclogite body may once have been situated in a continental subduction channel (the interface between the subducted slab and the wedge of an overlying crust and mantle) or a continental root zone (around Moho) owing to continental subduction and collision [2]. In this context, it can provide a valuable opportunity to examine the seismic properties of this rare type of eclogite and its implications for the detectability of eclogite bodies in the deep crust.

Geological Background
The Western Gneiss Region (WGR) in Norway is one of the largest, best-exposed, and most studied HP-UHP terranes in the world. It resulted from the deep subduction of the Baltica basement and subsequent collision with the Laurentia continent ≈425-400 Ma [34][35][36][37][38]. The peak metamorphic pressure and temperature (P-T) conditions increase from SE to NW with the highest P-T estimates preserved in three UHP domains (Nordfjord, Sørøyane, and Nordøyane) along the west Norwegian coast (Figure 1a). The Flemsøya is located in the Nordøyane UHP domain, which has the highest peak P-T condition (>800 • C and 3.0-4.0 GPa; see Terry et al. [39], Carswell et al. [40], and Butler et al. [41]) among the three UHP domains. The eclogite samples analyzed in this study were collected from northwestern Flemsøya in the Nordøyane UHP domain (Figure 1b). The Flem Gabbro (or Sandvikhaugane Gabbro) is the largest gabbro block that experienced heterogeneous eclogitization [33]. The studied eclogite (named as Flem eclogite hereafter) is located in the northernmost margin of the Flem gabbro. It is almost completely eclogitized and features variously deformed structures (massive and foliated structures) in the outcrop [42]. Our recent petrological work proposed that the Flem eclogite experienced UHP metamorphism at ≈2.7-3.7 GPa and ≈700-820 • C, in accordance with the previous P-T estimates in the Nordøyane UHP domain [42].  [43]). The dark gray shaded area includes three ultrahigh-pressure (UHP) domains (Nordfjord, Sørøyane, and Nordøyane) along the western coast of the WGR. The dashed lines are isobars of peak metamorphic pressures. The blue arrow indicates the Flemsøya region where our eclogite samples were collected. (b) Geological map of the NW part of Flemsøya showing the distributions of eclogite (modified from Mørk [33] and Terry and Robinson [38]). The blue arrow indicates the locality where our eclogite samples were collected.  [43]). The dark gray shaded area includes three ultrahigh-pressure (UHP) domains (Nordfjord, Sørøyane, and Nordøyane) along the western coast of the WGR. The dashed lines are isobars of peak metamorphic pressures. The blue arrow indicates the Flemsøya region where our eclogite samples were collected. (b) Geological map of the NW part of Flemsøya showing the distributions of eclogite (modified from Mørk [33] and Terry and Robinson [38]). The blue arrow indicates the locality where our eclogite samples were collected.

Sample Descriptions
Based on the field occurrences and deformation microstructures, the Flem eclogite could be divided into massive eclogite (MEC) and foliated eclogite (FEC). In the hand specimen, MEC hardly shows discernable foliation and lineation, whereas FEC developed these structural frameworks conspicuously. The MEC was mainly composed of olivine (Ol), garnet (Grt), omphacite (Omp), orthopyroxene (Opx), phlogopite (Phl), ilmenite (Ilm), amphibole (Amp), minor spinel (Sp), and plagioclase (Pl) ( Table 1). The Ol crystals were mostly (sub)euhedral and often rimmed with small Opx and Omp grains (Figure 2a). The Ol crystals lacked a shape-preferred orientation (SPO) and presented as an irregularly shaped Grt matrix (Figure 2a). In contrast, the FEC often displayed a pronounced foliation characterized by mineral layering, Grt bands, and elongated Ol porphyroclasts (Figure 2b). The mineral assemblage of FEC was similar to but had relatively higher proportions of Omp, Grt, and Amp, but lower Ol content, than the MEC (Table 1). Based on the field occurrences and deformation microstructures, the Flem eclogite could be divided into massive eclogite (MEC) and foliated eclogite (FEC). In the hand specimen, MEC hardly shows discernable foliation and lineation, whereas FEC developed these structural frameworks conspicuously. The MEC was mainly composed of olivine (Ol), garnet (Grt), omphacite (Omp), orthopyroxene (Opx), phlogopite (Phl), ilmenite (Ilm), amphibole (Amp), minor spinel (Sp), and plagioclase (Pl) ( Table 1). The Ol crystals were mostly (sub)euhedral and often rimmed with small Opx and Omp grains (Figure 2a). The Ol crystals lacked a shape-preferred orientation (SPO) and presented as an irregularly shaped Grt matrix (Figure 2a). In contrast, the FEC often displayed a pronounced foliation characterized by mineral layering, Grt bands, and elongated Ol porphyroclasts ( Figure 2b). The mineral assemblage of FEC was similar to but had relatively higher proportions of Omp, Grt, and Amp, but lower Ol content, than the MEC (Table 1).

Crystal-Preferred Orientation (CPO) Measurement and Analysis
The CPOs of the Flem eclogite were measured by investigating finely polished thin sections using an electron backscatter diffraction (EBSD) system. The thin sections, parallel to the lineation and perpendicular to the foliation (i.e., XZ plane), were used for the FEC, whereas the thin sections, parallel to the foliation of the neighboring FEC in the outcrop (named the apparent XY plane or X'Y' plane), were used for the MEC owing to its vague foliation and lineation. The EBSD system we used was installed in a Quanta 450 Field Emission Scanning Electron Microscope (FESEM) at the State Key Laboratory of Geological Processes and Mineral Resources in the China University of Geosciences (Wuhan), China. The operational settings of the EBSD system were: an acceleration voltage of 20 kV, a working distance of 25 mm, and a spot size of 60 in a low-vacuum mode. The CPOs were presented in the pole figures with marked foliation and/or lineation, and their strengths were quantified using the J-index (ranging from 1 to infinity; see [44]) and the M-index (ranging from 0 to 1; see [45]). The CPO analyses were completed using the MTEX toolbox (ver. 5.2.beta2) in MATLAB (ver. 2017b) (http://mtex-toolbox.github.io/) [46,47]. More details about the EBSD measurement and data treatment of studied samples are available in Cao et al. [42].

Seismic Property Calculation
To calculate the seismic property (velocity and anisotropy) of the MEC and FEC, the elastic stiffness, CPOs, densities, and volume proportions of their constituent minerals are required before inserting them into the Christoffel equation to solve for the seismic velocities (i.e., Vp, Vs1, and Vs2) along every direction in 3D space [48]. The volume proportions of the constituent minerals were obtained directly from the EBSD phase maps ( Table 1). The densities of the bulk rock were averaged from the densities of each mineral weighted by their volume proportions. The elastic stiffness of the single-crystal (Cij) Amp (kataphorite), Opx (enstatite80), Phl, Pl (albite), and Sp were adopted from Brown and Abramson [49], Webb and Jackson [50], Alexandrov and Ryzhova [51], Brown et al. [52], and Duan et al. [53], respectively. Because Omp in MEC and FEC has the composition of Jd20-40Di50-75 [42], the (Cij) of the Omp that has a similar composition (Jd30Di70) recently reported by Hao et al. [54] is used. The elastic stiffness of single-crystal Ol (avg. Fo# ≈ 68) was estimated using the (Cij) of their forsterite [55] and fayalite [56] endmembers, which were linearly averaged by their molar proportions. Likewise, the elastic stiffness of single-crystal Grt (Grs10-15, Prp40-45, and Alm+Sps45) was calculated using the (Cij) of grossular [57], pyrope [58], and alamdine-spessartite [59] endmembers. Since the single-crystal elastic stiffness of ilmenite (FeTiO 3 ) was not reported, we constructed a fictive isotropic tensor using the bulk and shear modulus of ilmenite (MgTiO 3 ) [60]. A list of single-crystal elastic stiffness tensors for common minerals was recently provided by Almqvist and Mainprice [61]. The elastic stiffness tensors <Cij> of mono-mineralic polycrystal and bulk rock were calculated by integrating the elastic stiffness tensor (Cij) of all grains (and phases) over their individual orientations using the Vogit-Reuss-Hill (VRH) averaging scheme. In this study, the calculation was accomplished using the MTEX toolbox (ver. 5.2.beta2) in MATLAB (ver. 2017b) [47,62]. Since this calculation method yields seismic velocities along discrete orientations in 3D space, we could derive the average Vp and Vs by simply averaging their maximum and minimum values ( Table 2).
To understand the contributions of different mineral phases and/or their CPOs to the bulk rock seismic properties, we also calculated the seismic properties of major mono-mineralic aggregates (e.g., olivine, orthopyroxene, garnet, omphacite, and amphibole). These results were compared with the seismic properties of bulk rock, providing a straightforward approach to characterizing the constructive and destructive effects among different mineral aggregates.

Crystal-Preferred Orientation
Both MEC and FEC exhibited obvious but weak CPOs for most of their constituent minerals (Figures 3 and 4). Notably, in most FEC samples, Opx, Omp, and Amp tended to align the maximum of their [001] axes subparallel to the lineation and the maximum of the [010] or [100] axes sub-perpendicular to the foliation (Figures 3 and 4). These CPO patterns are common and are also observed in many other deformed Opx, Omp, and Amp bearing rocks. The weak CPO strengths agreed with their fairly low J-(<2.5) and M-indices (<0.06). For more details of the EBSD analytical method, descriptions of the microstructures and CPOs in the studied MEC and FEC samples, readers are referred to Cao et al. [42].

Olivine Polycrystals
The Ol polycrystals in the MEC showed the maxima of Vp at low angles to the apparent foliation, subparallel either to the X' or Y' direction, while the Vp was low when a seismic ray propagated at moderate to high angles to the apparent foliation ( Figure S1a,b). The degree of shear wave splitting or S-wave polarization anisotropy and velocities of the fast shear wave (Vs1) displayed no regular patterns among MEC samples. In the FEC, the Vp, AVs, and Vs1 all lacked consistent distributions between the samples (Figure S1c-g). Regardless of the MEC and FEC samples, both the fast P-wave and the polarization direction of the fast shear wave (S1) tended to align parallel to the maxima of the [100] axes of the Ol CPOs ( Figure 3). Notably, the P-wave anisotropy (AVp) and maximum polarization anisotropy of the S-wave (max. AVs) were higher in the MEC (AVp = 2.04-3.42%, max. AVs = 1.69-2.53%) than the FEC (AVp = 0.86-1.65%, max. AVs = 0.84-1.28%), which is in agreement with the stronger Ol CPO strength in the MEC than the FEC ( Figure 3).

Olivine Polycrystals
The Ol polycrystals in the MEC showed the maxima of Vp at low angles to the apparent foliation, subparallel either to the X' or Y' direction, while the Vp was low when a seismic ray propagated at moderate to high angles to the apparent foliation ( Figure S1a,b). The degree of shear wave splitting or S-wave polarization anisotropy and velocities of the fast shear wave (Vs1) displayed no regular patterns among MEC samples. In the FEC, the Vp, AVs, and Vs1 all lacked consistent distributions between the samples (Figure S1c-g). Regardless of the MEC and FEC samples, both the fast P-wave and the polarization direction of the fast shear wave (S1) tended to align parallel to the maxima of the [100] axes of the Ol CPOs ( Figure 3). Notably, the P-wave anisotropy (AVp) and maximum polarization anisotropy of the S-wave (max. AVs) were higher in the MEC (AVp = 2.04-3.42%, max. AVs = 1.69-2.53%) than the FEC (AVp = 0.86-1.65%, max. AVs = 0.84-1.28%), which is in agreement with the stronger Ol CPO strength in the MEC than the FEC (Figure 3).

Orthopyroxene Polycrystals
The Vp, AVs, and Vs1 patterns of Opx polycrystals were irregular in the MEC, but they were more consistent in the FEC ( Figure S2). Except for one FEC sample ( Figure S2f), the max. Vp was aligned subparallel to the lineation and the low Vp tended to distribute in a girdle perpendicular to the lineation with their minima mostly subnormal to the foliation ( Figure S2c-e,g). Despite the lack

Orthopyroxene Polycrystals
The Vp, AVs, and Vs1 patterns of Opx polycrystals were irregular in the MEC, but they were more consistent in the FEC ( Figure S2). Except for one FEC sample ( Figure S2f), the max. Vp was aligned subparallel to the lineation and the low Vp tended to distribute in a girdle perpendicular to the lineation with their minima mostly subnormal to the foliation ( Figure S2c-e,g). Despite the lack of consistency in the distributions of AVs and Vs1 in the FEC, the fast S-wave always tended to polarize parallel to the lineation just like the fast P-wave ( Figure S2c-g). The magnitudes of the AVp and max. AVs were low and varied in similar ranges for both the MEC (AVp = 1.30%, max. AVs = 1.20-1.72%) and FEC (AVp = 0.84-1.52%, max. AVs = 0.82-1.70%).

Garnet Polycrystals
In both the MEC and FEC, the Grt polycrystals displayed orthogonal patterns for the max. and min. Vp ( Figure S3). The max. Vp mainly corresponded to the <100> axes in the Grt CPO (Figure 3), which were the fastest P-wave directions for the Grt single crystal. The AVs and Vs1 also displayed some regular clustering and girdling without consistent spatial correlation with the foliation and lineation. The AVp and max. AVs were very low in both the MEC (AVp = 0.03-0.04%, max. AVs = 0.06-0.09%) and FEC (AVp = 0.01-0.03%, max. AVs = 0.03-0.05%).

Omphacite Polycrystals
The Omp polycrystals exhibited disparate and irregular Vp, AVs, and Vs1 patterns in the MEC, whereas they were more consistent and regular in the FEC ( Figure S4). Specifically, the max. Vp was always subparallel to the lineation, complying with the maxima of the Omp [001] axes in the FEC (Figure 4). In contrast, the slow P-wave was perpendicular to the lineation with minima mainly at high angles to the foliation, matching the distribution of the Omp [010] pole ( Figure 4). The fast S-wave tended to polarize subparallel to the lineation when the seismic ray traveled at high angles to the lineation. The intensities of the AVp and max. AVs were also larger in the MEC (AVp = 3.09-4.98%, max. AVs = 2.09-5.29%) than in the FEC (AVp = 1.87-2.69%, max. AVs = 0.83-2.12%).

Amphibole Polycrystals
No consistent Vp, AVs, and Vs1 patterns were observed in the Amp polycrystals in the MEC ( Figure S5a,b). In contrast, these patterns were more regular in the FEC ( Figure S5). Excluding two FEC samples that showed XZ-girdled ( Figure S5f) and off-lineation clustering ( Figure S5g) of the fast Vp, the max. Vp was aligned subparallel to the lineation in the other FEC samples ( Figure S5c-e). Furthermore, the min. Vp was always found in the plane sub-perpendicular to the lineation. Notably, the max. and min. Vps were well correlated with the distributions of the [001] and [100] axes of Amp, respectively ( Figure 4). The AVs was large when the seismic ray propagated along the moderate Vp directions and with S1 polarizing parallel to the max. Vp directions. Except for one FEC sample that had the lowest AVp (1.58%) and max. AVs (1.45%) ( Figure S5g), the seismic anisotropies of the FEC samples (AVp = 3.29-16.95%, max. AVs = 2.56-12.00%) were much higher than those of the MEC (AVp = 2.29-2.80%, max. AVs = 2.00-2.44%).

Bulk Rocks
Overall, the bulk rock Vp distribution patterns of the MEC were similar to those of the Omp polycrystals (cf. Figures S4a,b and 5a,b). In contrast, three FEC samples (subgroup 1) presented Vp, AVs and Vs1 patterns resembling those of the Amp polycrystals (cf. Figures S5c-e and 5c-e); the other two FEC samples (subgroup 2) showed similar Vp, AVs, and Vs1 patterns to those of the Omp or Ol polycrystals (Figures S1f,g, S4f,g and 5f,g). Apart from one MEC sample (Figure 5a), the max. Vp and polarization direction of the fast S-wave were aligned subparallel (occasionally oblique) to the (apparent) lineation in other eclogite samples (Figure 5b (Figures 5f,g, 6 and 7), which exhibit the lowest AVp (0.99-1.09%) and max. AVs (0.63-0.65%) ( Table 2). Taking all samples together, the intensities of the AVp and max. AVs increased markedly with the volume proportion of Amp ( Figure 8) and CPO strength (Figure 9).

Factors Controlling the Seismic Anisotropy of Flem Eclogite
Several factors can affect the seismic anisotropy of a natural rock, including microcracks, mineral assemblage, and the strength of CPOs and their interactions [13,19,[63][64][65][66]. The role of microcracks is nontrivial in natural rocks; however, its contribution to seismic anisotropy decreases dramatically with depth, where the effect is negligible at pressures above ≈300 MPa (≈10 km in depth) owing to the fast closure of cracks under an increasing confining pressure [20,21,23], though a few studies suggest that some spherical micropores may still exist up to 1 GPa [67,68]. Nevertheless, these pressures are much lower than the case of our high-pressure MEC and FEC that were formed at ≈1.0-3.7 GPa or ≈30-110 km depths [42], where the effect of microcrack should be minimized. Because the single-crystal seismic anisotropies can be highly different between minerals, their volume proportions can tremendously influence the seismic anisotropies of the bulk rock that they constitute. This factor is generally the main contributor to the disparate seismic anisotropies between different rock types, such as blueschist and eclogite [12][13][14]. However, the bulk rock seismic anisotropies (AVp and max. AVs) of Flem eclogite show no obvious correlations with the volume proportions of major phases, except for Amp (Figure 8), indicating that bulk rock seismic anisotropies are not mainly

Factors Controlling the Seismic Anisotropy of Flem Eclogite
Several factors can affect the seismic anisotropy of a natural rock, including microcracks, mineral assemblage, and the strength of CPOs and their interactions [13,19,[63][64][65][66]. The role of microcracks is nontrivial in natural rocks; however, its contribution to seismic anisotropy decreases dramatically with depth, where the effect is negligible at pressures above ≈300 MPa (≈10 km in depth) owing to the fast closure of cracks under an increasing confining pressure [20,21,23], though a few studies suggest that some spherical micropores may still exist up to 1 GPa [67,68]. Nevertheless, these pressures are much lower than the case of our high-pressure MEC and FEC that were formed at ≈1.0-3.7 GPa or ≈30-110 km depths [42], where the effect of microcrack should be minimized. Because the single-crystal seismic anisotropies can be highly different between minerals, their volume proportions can tremendously influence the seismic anisotropies of the bulk rock that they constitute. This factor is generally the main contributor to the disparate seismic anisotropies between different rock types, such as blueschist and eclogite [12][13][14]. However, the bulk rock seismic anisotropies (AVp and max. AVs) of Flem eclogite show no obvious correlations with the volume proportions of major phases, except for Amp (Figure 8), indicating that bulk rock seismic anisotropies are not mainly

Factors Controlling the Seismic Anisotropy of Flem Eclogite
Several factors can affect the seismic anisotropy of a natural rock, including microcracks, mineral assemblage, and the strength of CPOs and their interactions [13,19,[63][64][65][66]. The role of microcracks is nontrivial in natural rocks; however, its contribution to seismic anisotropy decreases dramatically with depth, where the effect is negligible at pressures above ≈300 MPa (≈10 km in depth) owing to the fast closure of cracks under an increasing confining pressure [20,21,23], though a few studies suggest that some spherical micropores may still exist up to 1 GPa [67,68]. Nevertheless, these pressures are much lower than the case of our high-pressure MEC and FEC that were formed at ≈1.0-3.7 GPa or ≈30-110 km depths [42], where the effect of microcrack should be minimized. Because the single-crystal seismic anisotropies can be highly different between minerals, their volume proportions can tremendously influence the seismic anisotropies of the bulk rock that they constitute. This factor is generally the main contributor to the disparate seismic anisotropies between different rock types, such as blueschist and eclogite [12][13][14]. However, the bulk rock seismic anisotropies (AVp and max. AVs) of Flem eclogite show no obvious correlations with the volume proportions of major phases, except for Amp (Figure 8), indicating that bulk rock seismic anisotropies are not mainly controlled by the mineral assemblage.
The weak positive relations between the seismic anisotropies and the volume proportions of Amp, which is relatively stronger for AVp (R 2 = 0.583) than max. AVs (R 2 = 0.157), reveal that the high Amp content can partly account for the higher seismic anisotropies in the Amp-rich FEC subgroup 1 (Figure 8).
The stronger CPO strength (denoted by larger J-and M-indices) means that the crystal lattices were aligned more toward a similar orientation, resulting in greater seismic anisotropies. Despite a weak CPO strength, positive correlations between the bulk rock seismic anisotropies and the bulk rock CPO strength were still recognizable in the MEC and FEC (Figure 9). This result suggests that the variations of bulk rock seismic anisotropies in the Flem eclogite were largely affected by their CPO strength, which is also corroborated by previous studies, especially the variations of seismic anisotropies in the same rock type [13,14].
The interactions of CPO can either reinforce or weaken the bulk rock seismic anisotropies, depending on the distribution patterns of the constituent mono-mineralic polycrystals. If the fast velocities (Vp or Vs) of different mono-mineralic polycrystals are aligned in similar orientations (same for the slow velocities), the CPOs are thus thought to coherently contribute to the bulk rock seismic anisotropies (i.e., the constructive effect of CPO interactions). Otherwise, the seismic anisotropies of bulk rock are much weaker than those of mono-mineralic polycrystals, when the fast and slow velocities are subparallel and thus counteract with each other (i.e., destructive effect of CPO interactions) [13,69,70]. In the FEC subgroup 1, the Omp, Amp, and Opx polycrystals presented overall similar distribution patterns for Vp (constructive effect; Figures S2c-e, S4c-e and S5c-e), thus partly contributing to the highest bulk rock seismic anisotropies for the FEC subgroup 1. In contrast, the Vp distribution patterns were less consistent in the FEC subgroup 2 (destructive effect; Figures S2f,g, S4f,g and S5f,g), hence leading to the lowest bulk rock seismic anisotropies in collaboration with the lowest bulk rock CPO strength (Figure 9). The Vp patterns of Ol polycrystals mostly differed from those of Omp, Amp, and Opx ( Figure S1), which is indicative of the destructive effect of Ol CPO. The Grt polycrystals showed almost negligible seismic anisotropies ( Figure S3) and thus weakened the bulk rock seismic anisotropies. A comparison of the seismic anisotropies between bulk rock and constituent mono-mineralic polycrystals illustrated that the bulk rock seismic anisotropies were indeed the combined results of each of the mono-mineralic polycrystals ( Figure 6). Especially for the FEC subgroup 1, their highest seismic anisotropies were mainly due to the Amp polycrystals. The calculated seismic anisotropies were all smaller than the arithmetically averaged seismic anisotropies that approximated the purely constructive contributions of all constituent mono-mineralic polycrystals (Figure 7). This result clearly reflects the existence of destructive effects of CPOs on weakening the bulk rock seismic anisotropies [14]. The degree of deviations (calculated using 100 × AVp avg − AVp cal /AVp cal and 100 × max. AVs avg − max. AVs cal /max. AVs cal ) from the dividing line (k = 1) were much larger for the MEC (90-126% for AVp and 45-74% for max. AVs) and FEC subgroup 1 (34-45% for AVp and 25-54% for max. AVs) than FEC subgroup 2 (15-26% for AVp and 6-13% for max. AVs), implying that the destructive effects of the CPOs were not constant in the Flem eclogite, i.e., they were more remarkable in the MEC and FEC subgroup 1 than in FEC subgroup 2 ( Figure 7).
It is noteworthy that the three factors above are probably not mutually exclusive. The highest seismic anisotropies in FEC subgroup 1 were contributed jointly by their higher Amp contents, larger bulk rock CPO strength, and constructive CPO interactions. The enrichment of Amp in the FEC subgroup 1 may favor the development of a stronger Amp CPO (especially in the Amp-rich layers where strain is likely more localized) and result in a larger bulk rock CPO strength (combined consequences of stronger Amp CPO and higher Amp volume proportion). The stronger CPOs of Amp, Omp, and Opx in FEC subgroup 1 may also imply that more coherent CPO developments and resultant constructive CPO interactions reinforce the seismic anisotropies.

Seismic Properties and the Implications of Flem Eclogite
The Flem eclogite is a unique type of eclogite because of its uncommon mineral assemblage of abundant Ol and Opx. FEC subgroup 1 showed an AVp (1.78-2.46%) and max. AVs (1.30-1.99%) similar to previously reported dry eclogite that consisted of dominantly bi-mineralic Omp and Grt (AVp = 1.2-2.9%, max. AVs = 0.74-2.02%; Bascou et al. [22]), while the seismic anisotropies of the MEC (AVp = 1.06-1.08%, max. AVs = 1.00-1.38%) and FEC subgroup 2 (AVp = 0.99-1.09%, max. AVs = 0.63-0.65%) were even lower ( Figure 10). The seismic anisotropies of the Flem eclogite were also strikingly lower than the seismic anisotropies of other types of eclogite that bear hydrous phases, such as epidote (Ep), glaucophane (Gln), and lawsonite (Lws) (AVp = 1.4-10.2%, max. AVs = 1.35-8.09%; Figure 10). Therefore, compared to the eclogite free of Ol and Opx, and despite the large strain and high Amp content (as reflected by the Amp-rich FEC), the presence of Ol and Opx actually weakened the bulk rock seismic anisotropies, mainly by reducing the mineral CPO strength (Figure 9) through non-dislocation creep mechanisms, such as diffusion creep, grain or phase boundary sliding, and rigid-body-like rotations [42]. The Epand Lws-bearing blueschist, which is the hydrous precursor of eclogite, displayed the highest seismic anisotropies, implying a gradual decrease of the elastic anisotropies in the subducted oceanic crust with increasing depth [12][13][14]. Peridotite is the rock that is typically enriched with Ol and Opx. However, peridotite had a higher AVp and max. AVs (Figure 10), indicating that it can be distinguished from Flem eclogite and dry eclogite, but not Ep-, Gln-, and Lws-bearing eclogite, in terms of the magnitudes of seismic anisotropies [71]. The Flem eclogite is surrounded by non-or weakly metamorphosed gabbro (Figure 1b; see Mørk [33]). Although the seismic anisotropies of the surrounding gabbro are unknown, an examination of the seismic anisotropies of the gabbro from other localities revealed similar intensities to the peridotite [72] (Figure 10). This fact thus raises the feasibility of using the intensities of seismic anisotropies to distinguish Flem eclogite and dry eclogite from their gabbroic country rocks in the deep continental crust.  [22]), while the seismic anisotropies of the MEC (AVp = 1.06-1.08%, max. AVs = 1.00-1.38%) and FEC subgroup 2 (AVp = 0.99-1.09%, max. AVs = 0.63-0.65%) were even lower ( Figure 10). The seismic anisotropies of the Flem eclogite were also strikingly lower than the seismic anisotropies of other types of eclogite that bear hydrous phases, such as epidote (Ep), glaucophane (Gln), and lawsonite (Lws) (AVp = 1.4-10.2%, max. AVs = 1.35-8.09%; Figure 10). Therefore, compared to the eclogite free of Ol and Opx, and despite the large strain and high Amp content (as reflected by the Amp-rich FEC), the presence of Ol and Opx actually weakened the bulk rock seismic anisotropies, mainly by reducing the mineral CPO strength ( Figure  9) through non-dislocation creep mechanisms, such as diffusion creep, grain or phase boundary sliding, and rigid-body-like rotations [42]. The Ep-and Lws-bearing blueschist, which is the hydrous precursor of eclogite, displayed the highest seismic anisotropies, implying a gradual decrease of the elastic anisotropies in the subducted oceanic crust with increasing depth [12][13][14]. Peridotite is the rock that is typically enriched with Ol and Opx. However, peridotite had a higher AVp and max. AVs ( Figure 10), indicating that it can be distinguished from Flem eclogite and dry eclogite, but not Ep-, Gln-, and Lws-bearing eclogite, in terms of the magnitudes of seismic anisotropies [71]. The Flem eclogite is surrounded by non-or weakly metamorphosed gabbro (Figure 1b; see Mørk [33]).
Unlike the seismic anisotropies, the average seismic velocities (approximating the isotropic seismic velocities) and Vp/Vs ratio of Flem eclogite were controlled only by the mineral assemblage. Specifically, both Vp and Vs correlated positively and negatively with the volume proportions of Omp and Amp, respectively (Figure 11a,b). These salient relations were the manifestation of the faster seismic velocities of Omp than the Amp single crystal and the formation of Amp by replacing mainly Omp and Ol. As suggested by Cao et al. [42], the Amp in MEC and FEM has both low-and high-pressure origins, where the high-pressure Amp is the most abundant and is formed by deformation-enhanced hydration reactions under peak-to-early-retrograde UHP metamorphic conditions. Furthermore, despite narrow ranges of seismic velocities and densities, weak positive correlations between Vp, Vs, and bulk rock densities were still displayed (Figure 12), which obeyed Birch's law [9,29,80]. The seismic velocities of Flem eclogite were lower than those of peridotites and dry eclogite; larger than those of Ep-, Gln-, and Lws-eclogite; much greater than those of blueschist and gabbro ( Figure 13). The Vp and Vs of Flem eclogite were plotted in a straight line, yielding an average Vp/Vs ratio of 1.765, which was lower and higher than the Ep-/Gln-and Lws-bearing eclogites, respectively ( Figure 13). Despite the narrow range of the Vp/Vs ratio (≈1.758-1.770; see Table 2), a strong positive correlation (R 2 = 0.922) between the Vp/Vs ratio and Ol volume proportion suggested that Ol could significantly increase the Vp/Vs ratio of eclogite, whereas Amp and Grt decreased the bulk rock Vp/Vs ratio, as indicated by the negative correlations between the Vp/Vs ratio and their volume proportions (Figure 11c). The Vp/Vs ratio of the Flem eclogite was much larger than the non-retrograded eclogite (1.69-1.71) and slightly larger than the retrograded eclogites (1.73-1.77) that contain appreciable amounts of retrograde Amp and Pl [12,28]. This Vp/Vs ratio was also slightly higher than the dry eclogites (1.73-1.75) that are free of Ol and Opx and similar to those of peridotite [81], further indicating that Ol was the source of high Vp/Vs ratio in the Flem eclogite ( Figure 13). Owing to the abundance of Gln, which was characterized by the low Vp/Vs ratio, Ep-blueschist mostly had lower Vp/Vs ratios (1.70-1.73) than the Flem eclogite. However, the occurrence of a high-Vp/Vs-ratio Lws can increased the Vp/Vs ratio of Lws-blueschist to a wide range (1.71-1.81). The enrichment of high Vp/Vs ratio Pl caused the Vp/Vs ratios of gabbro (1.75-1.85) to be higher than the Flem eclogite. Therefore, we propose that the much larger Vp and Vs and relatively smaller Vp/Vs ratios of the olivine-rich eclogite (e.g., Flem eclogite) can be used to differentiate them from the country rocks (e.g., metastable and non-eclogized gabbro) in the deep continental crust or subduction channel when high-resolution seismic waves data are available since the detectable body size depends on the resolution (i.e., a wavelength that is determined by velocity and frequency) of seismic waves at different depths.     [13], lawsonite blueschsit [13,69,70], epidote/glaucophane eclogite [14,73,74], epidote blueschsit [14,70,73,74], dry eclogite [22], peridotite [81], (olivine) gabbro [82].

Conclusions
The Flem eclogite from NW Flemsøya in the Nordøyane UHP domain of the Western Gneiss Region in Norway is a unique type of eclogite that was enriched with olivine and orthopyroxene. Based on the detailed examinations on their seismic properties, several main findings were as follows.

•
The seismic anisotropies of the Flem eclogite were largely controlled by the CPO strength and significant destructive effects of CPO interactions on weakening the seismic anisotropies were manifested in both weakly and strongly foliated eclogite. Amphibole CPOs were the main contributors to the higher seismic anisotropies in some amphibole-rich and strongly foliated eclogite.

•
The average seismic velocities correlated positively with the rock density and omphacite content, but negatively with the amphibole volume proportion, which is indicative of the dominant role of a mineral assemblage in governing the bulk rock seismic velocities.

•
The Vp/Vs ratios were almost constant (≈1.765) between the Flem eclogite samples. These values were similar to the Vp/Vs ratios of peridotite and slightly higher than the dry eclogite that was free of olivine, suggesting that the abundance of olivine was the source of their high Vp/Vs ratios.

•
Compared to other common rock types, including hydrous-phase-bearing eclogite, blueschist, peridotite, and gabbro, the AVp and max. AVs of Flem eclogite were very low and similar to those of dry eclogite. In contrast, the average Vp and Vs of the Flem eclogite were remarkably larger than the hydrous-phase-bearing eclogite, blueschist, and gabbro, but lower than the dry eclogite and peridotite. The Vp/Vs ratios of the Flem eclogite were relatively higher than those of dry, non-retrograded, and retrograded eclogite, but significantly lower than gabbro. These results thus highlight the feasibilities of utilizing seismic anisotropies, velocities, and Vp/Vs ratios for differentiating olivine-rich eclogite (e.g., Flem eclogite) bodies from their country rocks (e.g., gabbro) in the deep continental crust or subduction channel when high-resolution seismic wave data are available.
Supplementary Materials: The following are available online at http://www.mdpi.com/2075-163X/10/9/774/s1, Figure S1: Seismic properties of the olivine aggregates in the MEC and FEC and olivine single crystal at ambient condition; Figure S2: Seismic properties of the orthopyroxene aggregates in the MEC and FEC and orthopyroxene single crystal at ambient condition; Figure S3: Seismic properties of the garnet aggregates in the MEC and FEC and garnet single crystals in MEC and FEC at ambient condition; Figure S4: Seismic properties of the omphacite aggregates in the MEC and FEC and omphacite single crystal at ambient condition; Figure