Revision 1 Magnetic Anisotropy in Natural Amphibole Crystals Andrea
1 Revision 1 2 Magnetic anisotropy in natural amphibole crystals 3
4 Andrea R. Biedermann1,a, Christian Bender Koch2, Thomas Pettke3 and Ann M. Hirt1,*
5 1 Institute of Geophysics, ETH Zurich, Sonneggstrasse 5, 8092 Zurich, Switzerland; 6 [email protected], [email protected]
7 a Now at Department of Geology and Mineral Resources Engineering, Norwegian University 8 of Science and Technology, Sem Sælands vei 1, 7491 Trondheim, Norway
9 2 Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 10 Copenhagen Ø, Denmark; [email protected]
11 3 Institute of Geological Sciences, University of Bern, Baltzerstrasse 1-3, 3012 Bern, 12 Switzerland; [email protected]
13 Address for Correspondence: Prof. Dr. Ann M. Hirt, Institute of Geophysics, ETH Zurich, 14 Sonneggstrasse 5, 8092 Zurich, Switzerland, phone: +41446332705, fax: +41446331065
15 Abstract 16 Anisotropy of magnetic susceptibility (AMS) is often used as a proxy for mineral fabric in 17 deformed rocks. In order to do so quantitatively, it is necessary to quantify the intrinsic 18 magnetic anisotropy of single crystals of rock-forming minerals. Amphiboles are common in 19 mafic igneous and metamorphic rocks and often define rock texture due to their general 20 prismatic crystal habits. Amphiboles may dominate the magnetic anisotropy in intermediate to 21 felsic igneous rocks and in some metamorphic rock types, because they have a high Fe 22 concentration and they can develop a strong crystallographic preferred orientation. In this 23 study, the AMS is characterized in 28 single crystals and one crystal aggregate of 24 compositionally diverse clino- and ortho-amphiboles. High-field methods were used to isolate 25 the paramagnetic component of the anisotropy, which is unaffected by ferromagnetic 26 inclusions that often occur in amphibole crystals. Laue imaging, laser ablation inductively 27 coupled plasma mass spectrometry and Mössbauer spectroscopy were performed to relate the 28 magnetic anisotropy to crystal structure and Fe concentration. The minimum susceptibility is 29 parallel to the crystallographic a*-axis and the maximum susceptibility is generally parallel to 30 the crystallographic b-axis in tremolite, actinolite, and hornblende. Gedrite has its minimum 31 susceptibility along the a-axis, and maximum susceptibility aligned with c. In richterite, 32 however, the intermediate susceptibility is parallel to the b-axis and the minimum and 33 maximum susceptibility directions are distributed in the a-c-plane. The degree of anisotropy, 34 k’, increases generally with Fe concentration, following a linear trend described by: k’ = 35 1.61×10-9 Fe – 1.17×10-9 m3/kg. Additionally, it may depend on the Fe2+/Fe3+ ratio. For most 36 samples, the degree of anisotropy increases by a factor of approximately 8 upon cooling from 37 room temperature to 77 K. Ferroactinolite, one pargasite crystal and riebeckite show a larger 38 increase, which is related to the onset of local ferromagnetic (s.l.) interactions below about 39 100 K. This comprehensive data set increases our understanding of the magnetic structure of 40 amphiboles, and it is central to interpreting magnetic fabrics of rocks whose AMS is 41 controlled by amphibole minerals.
42 1
43 Keywords: AMS (Anisotropy of magnetic susceptibility), magnetic properties, single crystal, 44 amphibole, hornblende, actinolite, richterite, tremolite
45 1. Introduction 46 Members of the amphibole group are common rock forming minerals occurring in a
47 wide range of igneous and metamorphic rocks. Amphiboles crystallize generally in
48 idiomorphic, prismatic to needle-like habits; hence they often display preferential orientation
49 in a deformed rock. This, combined with their intrinsic magnetic anisotropy, often makes
50 amphiboles, together with phyllosilicates, the main carriers of magnetic anisotropy in igneous
51 and metamorphic rocks. Amphiboles can be responsible for the magnetic fabric of a rock in
52 two ways. Firstly, the amphibole minerals themselves can dominate the paramagnetic
53 anisotropy (e.g. Borradaile et al. 1993; Schulmann and Ježek 2011; Zak et al. 2008).
54 Secondly, the shape of magnetite inclusions can be controlled by the crystallographic
55 preferred orientation (CPO) of amphiboles, which in turn is responsible for the magnetic
56 anisotropy (Archanjo et al. 1994).
57 Because amphiboles possessing a CPO can be an important carrier of the magnetic
58 anisotropy in a rock, it is essential to quantify their intrinsic anisotropy of magnetic
59 susceptibility (AMS). Until now, only a few studies have been conducted on the AMS of
60 amphibole single crystals, returning inconsistent results. Finke (1909) measured one
61 hornblende crystal and found the maximum susceptibility at an angle of -21°55’ to the
62 crystallographic c-axis. Wagner et al. (1981) measured the magnetic anisotropy in six crystals
63 from the hornblende group. In addition, they cited an unpublished study by Parry (1971), who
64 examined high-field AMS in 18 hornblende crystals. Both studies concluded that the
65 maximum principal susceptibility is sub-parallel to the crystallographic c-axis. However, in
66 Wagner’s hornblende, the minimum susceptibility aligns with the crystallographic a-axis,
67 whereas it is parallel to b in Parry’s study. This difference was attributed to the presence of
68 ferromagnetic inclusions that influenced Wagner’s low-field measurements (Wagner et al.
2
69 1981). Borradaile et al. (1987) measured five aggregates of amphibole crystals, including two
70 actinolites and one sample each of hornblende, crocidolite (fibrous riebeckite) and
71 glaucophane. Due to imperfect alignment of the individual grains within the aggregates, this
72 study gives an estimate of the lower limit of the AMS. The authors provide no directional
73 dependency. Lagroix and Borradaile (2000) measured two pargasite crystals and suggested
74 that the maximum susceptibility is sub-parallel to the b-axis.
75 Differences in orientations of the principal axes of the AMS ellipsoid reported in these
76 previous studies illustrate the importance of systematically investigating the magnetic
77 anisotropy of amphiboles. In the present study, the intrinsic magnetic anisotropy of a series of
78 amphiboles having a wide range of chemical compositions is characterized. The magnetic
79 anisotropy was measured in low and high magnetic fields and at different temperatures in
80 order to isolate the paramagnetic AMS. The paramagnetic AMS is then interpreted based on
81 the general crystal structure of the amphiboles and their chemical composition. A main focus
82 is put on the dependence of AMS on the Fe concentration.
83 2. Material and methods
84 Amphibole is an inosilicate that has the general formula A0-1B2C5T8O22(OH,F,Cl)2,
85 where A = Na, K; B = Ca, Na, Fe2+, Mn, Li, Mg; C = Mg, Fe2+, Fe3+, Al, Ti, Mn, (Ni, Cr, V,
6- 86 Li, Zn); and T = Si, Al. The main structural element defining amphiboles are [(T4O11) ]n
87 chains, i.e. double chains of tetrahedrally coordinated silica or aluminum. All tetrahedra of the
88 same double chain point in the same direction, whereas two double chains with oppositely
89 pointing tetrahedra are bonded by a band of octahedrally coordinated cations (C) (Figure 1).
90 These cations occupy one of three sites (labeled M1, M2 and M3), which possess variable
91 distortions of the octahedra, depending on the local environment. The M3 octahedra share 6
92 edges with adjacent octahedra, the M2 sites share 3 edges with octahedra and one with each of
93 the two M4 polyhedra. M1 shares 5 edges with octahedra and 1 with the M4 polyhedron. The
3
94 sizes of the M1, M2 and M3 sites depend on the radius of the cation a given site hosts.
95 Additionally, the sizes of M1 and M3 octahedra are also influenced by the amount of (OH)-
96 substitution by F- and Cl-. The two tetrahedral double chains and the octahedral strip form so-
97 called I-beams. Neighboring I-beams are bonded by cations in M4 (B) and A sites. The
98 cations in the M4 sites have usually a higher than 6-fold coordination, and the coordination
99 number is determined by the size of the cation. The A-site is coordinated by 12 surrounding
100 oxygen atoms and can be vacant, partially occupied or filled.
101 Amphiboles are subdivided into 4 main groups, according to the main M4- or B-cation.
102 Amphiboles that are classified as calcic (1), calcic-sodic (2) or sodic (3) have large cations,
103 Ca and Na, in varying proportions, in the M4 sites. The M4 site in the Fe-Mg-Mn amphiboles
104 (4) is occupied by smaller cations. Depending on the size of the cations in the M4 sites, the
105 stacking sequence changes and thus amphiboles can possess a monoclinic or orthorhombic
106 unit cell. Clinoamphiboles (space group C2/m) are more common than orthoamphiboles
107 (Pnma). The symmetry elements as defined by the space groups can dictate the orientation of
108 principal susceptibility directions according to Neumann’s principle (Neumann 1885), which
109 states that any physical property of a crystal has to include all symmetry elements of its space
110 group. Therefore, each principal susceptibility has to be parallel to one of the crystallographic
111 axes for the orthoamphiboles, while for clinoamphiboles one principal susceptibility has to be
112 parallel to the crystallographic b-axis (Nye 1957).
113 With respect to the magnetic properties, it is the location and arrangement of Fe atoms
114 that will be of greatest importance. This is due to the large magnetic moment of Fe in
115 combination with a relatively high abundance in the crystal lattice. Iron can be present as Fe2+
116 or Fe3+, whereby the Fe3+/Fe2+ ratio rarely exceeds 1/2. Exceptions to this general rule can be
117 found in Fe-rich hornblende (often referred to as oxy-hornblende) and hastingsite. In
118 actinolite, Fe2+ prefers M1 and M3 over M2, and some Fe2+ can also be located at M4 (Deer et
4
119 al. 1997). Hornblende generally shows a similar Mg/Fe ratio in each of the M1, M2 and M3
120 sites, and Fe2+ can be located in M4. Fe3+, like other small cations (e.g. Ti, Al), is
121 preferentially located at the M2 sites (Deer et al. 1997). Metamorphic and skarn hornblende
122 are different; in these minerals, Mg shows a preference for M2 and Fe has a tendency to be
123 located in M1 or M3. In richterite, Fe2+ is enriched in the M2 sites compared to the M1 or M3
124 sites, where the Fe/Mg ratios are similar. If Mn is present, it occupies the M4 or M2 sites. In
125 contrast, Fe2+ prefers M1 and M3 in riebeckite, where some Fe2+ can also enter M4. Fe3+ in
126 riebeckite is located in M2 sites. In the anthophyllite-gedrite series, Fe2+ shows a preference
127 for M4 over M1, M2 and M3, which possess similar Fe2+ concentrations. If there is a
128 difference between M1, M2 and M3, the Fe2+ concentration is lowest in M2. When OH is
129 substituted by F or Cl, this can also influence the site occupancy; Fe-F avoidance forces a
130 larger proportion of Fe2+ into the M2 sites, whereas Cl prefers bonds to Fe2+ over bonds to
131 Mg, which results in a preferential ordering of Fe2+ into the M1 and M3 sites (Deer et al.
132 1997).
133 2.1 Sample description
134 A collection of 28 natural single crystals and one aggregate was analyzed, chosen to
135 cover a broad range of chemical compositions (Table 1). These include crystals from the
136 calcic amphiboles (tremolite, actinolite and hornblende groups), sodic-calcic (richterite) and
137 sodic (riebeckite) (clino-)amphiboles and one Fe-Mn-Mn orthoamphibole (gedrite). The
138 individual crystals and the aggregate were cleaned with ethanol in an ultrasonic cleaner and
139 oriented prior to measurements. Crystals were oriented based on Laue X-ray diffraction,
140 performed at the Laboratory of Crystallography, ETH Zurich. Laue images were analyzed
141 with the OrientExpress 3.4 crystal orientation software (Laugier and Filhol 1983). The
142 oriented crystals were glued into diamagnetic plastic cylinders with no magnetic anisotropy.
143 The crystallographic a* (or a for the orthoamphibole), b and c directions corresponded to the
5
144 sample x, y and z-axes, respectively for magnetic measurements. The accuracy of crystal
145 orientation is ± 5°. For the riebeckite aggregate, only the c-axis, which corresponds to the
146 long axis of individual fibrous crystals, was oriented.
147 2.2 Chemical analysis
148 Bulk chemical composition was determined using laser ablation inductively coupled
149 plasma mass spectrometry (LA-ICP-MS) at the Institute of Geological Sciences, University of
150 Bern. Because most crystals could not be modified (e.g. for producing polished sections), the
151 analyses were performed on crystal surfaces or cleavage planes, and for some crystals also on
152 polished cross-sections (cf. Table A). LA-ICP-MS was preferred over electron probe
153 microanalysis (EPMA) because LA-ICP-MS samples a bigger volume (here a cylinder of 90
154 µm diameter and about 60 µm depth). Hence minute impurities that may greatly affect the
155 magnetic properties, such as exsolutions, melt or mineral inclusions, are also comprised in the
156 mineral analysis. Notably metamorphic amphiboles often contain mineral inclusions, referred
157 to as poikilitic growth. Per sample, four to six analyses were performed to obtain information
158 on element homogeneity, expressed as one standard deviation (SD) in Table A (online
159 supplementary).
160 The LA-ICP-MS used consists of a GeoLas Pro system coupled with an Elan DRC-e
161 ICP-MS, optimized following procedures detailed in Pettke et al. (2012). SRM610 from NIST
162 was used as the external standard material, with preferred values reported by Spandler et al.
163 (2011). Data reduction employed SILLS (Guillong et al. 2008), and internal standardization
164 was done by normalizing to 98 wt% total major oxides, thus allowing for 2 wt% of OH, F or
165 Cl that cannot be measured with LA-ICP-MS. Data are considered to be accurate to better
166 than 2 % 1 SD. The analytical accuracy is thus insufficient to reliably quantify Fe3+ and Fe2+
167 by LA-ICP-MS.
6
168 In order to determine the relative proportions of Fe2+ and Fe3+ and help in assigning the
169 Fe to the various crystallographic sites Mössbauer spectra were measured on selected samples
170 at the Department of Chemistry, University of Copenhagen. Absorbers were prepared by
171 mixing powdered mineral samples and boron nitride into PerspexR sample holders. Spectra
172 were obtained at room temperature using a conventional constant acceleration spectrometer
173 with the absorber perpendicular to the γ-ray direction and samples FAkt1 and FAkt4 were
174 also measured at the magic angle (tilted at 54.7⁰ relative to the γ-rays), to remove effects of
175 mineral alignment in the powder. The spectrometer was calibrated using the spectrum of a
176 thin foil of natural Fe at room temperature, and isomer shifts are given relative to the center of
177 this absorber. All spectra exhibit three absorption lines which are interpreted as being
178 composed of three overlapping doublets: one due to Fe3+ and two due to Fe2+. The maximum
179 intensities of the absorption lines were between 3% and 7% and baseline counts in the folded
180 spectra were between 2 and 10 million counts. Lines were deconvoluted assuming Lorentzian
181 lineshape and constraining the lines of each component to be identical. This constraint could
182 not be applied to samples FAkt1 and FAkt4 due to preferred orientation of the grains in the
183 powdered sample. New absorbers were prepared and measured using the magic angle and
184 these spectra were suitable for the constraints (equal width and intensity of the two lines in
185 the doublet). Assuming identical f-factors for each of the Fe components, the relative spectral
186 areas were converted into abundances. The samples were scanned using amplitudes of 5 and
187 12 mm/s to achieve good spectral resolution of the amphibole components and to allow for
188 checking for the presence of inclusions of magnetically ordered Fe oxides.
189 Site occupancies for individual cations were then determined based on the general
190 formula A0-1B2C5T8O22(OH,F)2. Because the A site can be vacant or partially filled, the
191 recalculation is not straightforward, unless the Fe2+/Fe3+ ratio is known. For each sample, one
192 of the following three models was used:
7
193 (1) When Mössbauer data were available, i.e., the Fe3+ and Fe2+ concentrations are
194 known, the site distribution was defined based on a charge balance, setting the total
195 number of cation charges to 46.
196 (2) For orthoamphiboles, it was assumed that no Na is located in M4. In this case, the
197 number of cations minus Na and K equals 15.
198 (3) For clinoamphiboles, it was assumed that no Mg, Mn or Fe2+ occupy the M4 site and
199 the total number of cations minus Ca, Na and K is equal to 13.
200 2.3 Magnetic measurements
201 Characterization of ferromagnetic inclusions. All magnetic measurements were made
202 at the Laboratory for Natural Magnetism, ETH Zurich. Acquisition of isothermal remanent
203 magnetization (IRM) was measured on selected samples in order to check for the presence of
204 and identify ferromagnetic inclusions within the crystals. Samples were magnetized along the
205 –z direction with an ASC Scientific IM-10-30 Pulse Magnetizer in a field of 2 T.
206 Subsequently, the sample was remagnetized along the +z direction in increasing fields
207 between 20 mT and 2 T. After each magnetization step, the IRM was measured on a 2G
208 Enterprises, three-axis cryogenic magnetometer, Model 755.
209 Low-field AMS and mean mass susceptibility. Magnetic susceptibility can be
210 described by a symmetric second-order tensor with eigenvalues , or | |
211 | | | | in the case of diamagnetic samples, i.e. corresponds to the most and to the
212 least negative susceptibility for diamagnetic samples (Hrouda 2004). The corresponding
213 eigenvectors define the directions of the principal susceptibilities. This can be represented by
214 a magnitude ellipsoid of susceptibility, often referred to as anisotropy ellipsoid. The ellipticity
215 can be described by the AMS degree P = k1/k3, or k’