Tin-bearing magnetite with nanoscale Mg-Si defects: Evidence for the early stages of mineralization in a skarn system

Tin-bearing magnetite is reported from several types of magmatic-hydrothermal ore deposits. The question of whether tin is incorporated within solid solution, as Sn4+, or as nanoinclusions remains open, however. We report a micron- to nanoscale investigation of Sn (Mg, Si)-bearing magnetite from serpentinite in the Dulong Zn-Sn-In skarn, South China, with the dual aims of understanding the mechanisms involved in accommodating Sn and associated elements into the Fe-oxide, and the inferences that this carries for constraining the early stages of skarn formation. Magnetite preserves a range of textures that record the evolution of metasomatism during prograde growth of grain cores and retrograde rim replacement. Observations reveal the presence of chondrodite and sellaite (MgF2) as nanoscale inclusions preserved in magnetite. This implies initiation of the Dulong mineralizing system during a humite-bearing, magnesium skarn stage. Magnesium-Si defects, forming along (110) planes prior to Sn-enrichment, are recognized for the first time. Release of high volatile, F-rich fluids is interpreted to lead to precipitation of cassiterite inclusions along <111*> directions in magnetite.


Introduction
Magnetite is a ubiquitous mineral in a wide range of different types of ore deposit, including Fe-Cu skarns, iron oxide-copper-gold (IOCG) and iron oxide apatite (IOA) systems, banded iron formations (BIF), and magmatic Fe-Ti-V and Ni-Cu-PGE deposits. Trace element geochemistry, which is commonly used to define distinct sub-types or generations of accessory minerals (e.g., Mao et al., 2016;Kontonikas-Charos et al., 2019) is now also extensively applied to magnetite (e.g., Wang et al., 2018). Magnetite is a refractory mineral and represents an untapped record of genetic information that is preserved as incorporated trace elements, mineral inclusions at various scales, as well as structural defects. Thus, nanoscale studies may represent a particularly valuable auxiliary tool to diagnose the earliest stage of evolving hydrothermal systems (Cook et al., 2022a). This is in agreement with the significant advantage to be gained by combining multiple techniques to adequately correlate petrographic and geochemical evidence across different scales of observation (Cook et al., 2017).
The incorporation of trace elements, including Mg, Al, Si, Ca, Ti, V, Cr, Sn and W, into magnetite, has prompted many microbeam studies focusing on fingerprinting of deposit types, thus allowing evaluation of genetic models (e.g., Neumann et al., 2017;Yin et al., 2019;Verdugo-Ihl et al., 2021). Recent nanoscale studies of magnetite from IOCG and IOA systems have used high angle annular dark field scanning transmission electron microscopy (HAADF STEM) imaging to better understand the physical nature of these trace elements. It has been shown that such elements, notably Si, are not always lattice-bound, rather, they can occur as sub-micron scale nanoparticles (Deditius et al., 2018;Ciobanu et al., 2019Ciobanu et al., , 2022Verdugo-Ihl et al., 2020;Huang and Beaudoin, 2021;Cook et al., 2022b).
A comparable approach can be undertaken for Sn-bearing magnetite. Tin, exsolved as cassiterite inclusions in magnetite has been long recognized from micron-scale investigations (Desborough and Sainsbury, 1970;Wang et al., 2012;Wang et al., 2018). Although published data show that average Sn concentrations in magnetite are generally no more than a few tens of ppm (e.g., Nadoll et al., 2014), concentrations commonly reach hundreds and even thousands of ppm in relatively Sn-rich systems (e.g., Huanggang Fe-Sn skarn, China; Mei et al., 2014;Olympic Dam IOCG, Australia,Verdugo-Ihl et al., 2020). As proposed by Desborough and Sainsbury (1970), high concentrations of Sn in magnetite represents a geochemical tracer for, and pathfinder to, Sn mineralization. Synthesis experiments at 1,400°C demonstrate that magnetite can be doped with Sn 4+ giving the nominal composition Sn x Fe 3-x O 4 , where 0<x<0.6 (Hernández- Gómez et al., 2001). Substitution mechanisms for Sn incorporation into magnetite as Sn 4+ may be either require a vacancy (e.g., Barkov et al., 2008), or be part of a coupled substitution with other divalent elements (e.g., Chen et al., 1992;Wang et al., 2018).
Magnetite, containing inclusions of cassiterite, was reported from the Dulong Zn-Sn-In deposit, Southwest China (Xu et al., 2021a). Whereas magnetite and cassiterite are distinct components of the calcic skarn in this deposit, the presence of cassiterite inclusions is limited to serpentine-bearing rocks at the contact to marble. Such Sn-bearing magnetite resembles the cassiterite-in-magnetite reported from small bodies of tactites (rocks of complex mineralogical composition, formed by contact metamorphism and metasomatism of carbonate rocks) in fractured marble near the granite at Lost River tin Mine (Alaska; Desborough and Sainsbury, 1970). Wang et al. (2018) reported stanniferous magnetite from calcic skarn assemblages in the Haobugao Fe-Zn skarn deposit (Great Xing'an Range, China) and contended that Sn-enriched, cassiterite-bearing magnetite forms as temperature decreases. In contrast, the preferential distribution of Sn in magnetite from serpentine rocks at Dulong, but conspicuously not in the calcic skarn, prompts the question of whether such an association can instead be related to an initial skarn stage. This same stage is also characterized by the highest concentrations of indium in sphalerite from pyrrhotite-bearing hornfels (biotite schist) at the contact with marble (Xu et al., 2021b).
Tin-bearing magnetite from Dulong has been investigated down to the nanoscale to understand its origin and association with serpentine rocks. Skarn systems are the result of fluid-rock interaction between a (typically) carbonaceous protolith and hydrothermal fluids, and result in patterns of mineralogical zonation generally defined by calc-silicates such as garnet and pyroxene (e.g., Einaudi et al., 1981;Meinert, 1992;Meinert et al., 2005) but also high-temperature Mg-silicates, such as forsterite, humite, and/or chondrodite, (e.g., Cioflica et al., 1971Cioflica et al., , 1977Gower et al., 1985;Ciobanu and Cook, 2004;Yao et al., 2014). Although not always well preserved, and indeed often not recorded at all, the occurrence of Mg-silicates is convincing evidence for an early magnesian skarn stage and may, in some cases, be associated with Au, Cu, W, Mo or B mineralization (e.g., Zhao et al., 1999;Cook and Ciobanu, 2001;Houzar et al., 2008;Marincea and Dumitraş, 2019;Huang et al., 2022). Intense retrogression (e.g., serpentinization) generally obliterates these minerals, particularly in F-rich systems (e.g., Houzar et al., 2008;Huang et al., 2022). Changes in redox conditions during the switch between prograde (diopside+andradite) and retrograde (hedenbergite+grossular) skarn stages also shows enrichment in elements such as W and Mo (Xu et al., 2016. In this contribution, we bridge micron-and nanoscales of observation using complementary techniques to address Sn-rich magnetite from serpentinite in the Dulong skarn, South China (Xu et al., 2021a;2021b). Our objectives are to 1) characterize Sn, Mg, Si-bearing magnetite down to the atomic scale, and 2) to discuss its petrogenetic significance. We show that atypical Mg-Si defects predate and promote the Sn-enrichment in magnetite that is expressed by inclusions of cassiterite occurring alongside F-and Mg-bearing phases (chondrodite and sellaite) and serpentine. We discuss the findings in the context of an evolving skarn system.

Geological background
The Dulong deposit belongs to the southwestern part of the Yangtze Block ( Figure 1A). A U-Pb concordia intercept age of 89.4 ± 1.4 Ma for cassiterite relates mineralization to the first phase of the Laojunshan granite ( Figure 1B), which has a zircon U-Pb age of 87.5-92.9 Ma (Liu et al., 2007;Xu et al., 2015;Zhao et al., 2018). Geochemical data suggest that the Laojunshan granite is a highly fractionated S-type granite derived from partial melting of Paleoproterozoic-Mesoproterozoic pelitic basement rocks (e.g., Xu et al., 2015). The Xinzhaiyan Formation is the main host for ore and consists of interstratified quartz-mica schist, calcareous and dolomitic marble ( Figure 1B). The dolomitic proportion in the marble increases with depth in the sequence (e.g., Ye et al., 2018).
The Dulong orefield contains four ore blocks (Tongjie, Manjiazhai, Jinshupo and Lazizhai) with a total contained resource of 5.0 Mt Zn, 0.4 Mt Sn, and 7 kt In (Zhao et al., 2018). Xu et al. (2021a) provided an overview of deposit zoning, skarn mineralogy and mineralization. The serpentinemagnetite skarn described here is concentrated in the relatively deep and south parts of the Manjiazhai ore block ( Figures 1B,C).
Due to extensive retrogression, no prograde forsterite is observed, although previous publications (e.g., Li et al., 2018) hint at early prograde forsterite at Dulong. Xu et al. (2021a) presented data showing that, apart from cassiterite, skarn silicates (notably garnet) are important carriers of Sn. Subordinate Sn is also contained in iron-oxides and sulfides.
The LA-ICP-MS trace element concentration data for magnetite from serpentinite given by Xu et al. (2021a) show concentrations of Sn, Zn and W that are an order of magnitude Frontiers in Earth Science frontiersin.org 03 greater than in magnetite from marble ( Figure 2). Likewise, cassiterite from serpentinite is richer in W, In and U, relative the same mineral in marble. Such differences imply formation of two distinct stages of magnetite and associated cassiterite during skarn evolution, if they can be linked to the mineralizing event rather than simply being signatures inherited from pre-existing lithologies.

Materials and methodology
All analytical work was carried out at Adelaide Microscopy, University of Adelaide. A Cameca SX-Five electron probe microanalyser (EPMA) was used to obtain quantitative compositional data and perform qualitative mapping of magnetite. The instrument was operated at an accelerating voltage of 15 kV and a beam current of 20 nA using a slightly defocused beam of 1 µm for quantitative analysis. Analytical conditions for mapping were adjusted to an accelerating voltage of 20 kV, a probe current of 100 nA, and dwell time of 100 ms.
HAADF STEM and energy-dispersive X-ray (EDX) mapping were performed using an ultra-high resolution, probe-corrected, FEI Titan Themis S/TEM (e.g., Ciobanu et al., 2016). FEI software was used for image acquisition (TIA v4.15), with supporting use of the drift-corrected frame integration package (DCFI) included in the Velox (v2.13.0.1138) software. Various filters (Radial Wiener, Highpass, Average and Gaussian blur) were used to eliminate noise and/ or enhance images. All EDS data acquisition and processing was carried out using Velox software.
HR STEM imaging of beam-sensitive minerals was achieved using the integrated differential phase contrast (iDPC) technique, described by Lazić et al., 2016. For thin, non-magnetic samples the iDPC technique images the scalar electrostatic potential field, yielding images of the atomic species with contrast roughly proportional to the atomic number (Z). In contrast to HAADF STEM, the detectability of light elements is significantly improved. The iDPC signal is also superior to the HAADF signal, due to the use of the intense signal in the bright field disc and the noise-suppressing nature of the integration process. As a result, the iDPC imaging technique is ideal for imaging beam-sensitive materials.

Results
Magnetite from serpentinite displays a broad range of intra-grain textures that correlate with measured Sn concentrations ( Figure 3). Magnetite featuring Sn-rich cores (Figures 3A,B) is associated with In-bearing sphalerite (Xu et al., 2021b). These cores can be very small, resorbed, and often overgrown by asymmetric margins ( Figure 3C), or display a patchiness or mottling expressed by submicron-scale inclusions of cassiterite ( Figures 3D,E). Different spatial arrangements of the inclusions define an intertwined lamellar substructure in the patchy cores, as well as inclusion-free subdomains in the mottled cores ( Figures 3F,G).
A representative Sn-bearing magnetite grain from serpentinite sample DL66 featuring patchy core zonation was assessed by EPMA spot analysis and element mapping at the micron scale (Table 1, Figure 4). A profile across the inclusionfree rim ( Figure 4B) shows concentrations of MgO, SiO 2 and SnO 2 (up to several wt% close to the core boundary) that correlate antithetically with FeO*. The patchy texture within the core is reproduced on the Sn map, whereas the rim is marked by a band of relatively higher Mg and Si, and scattered Si-Mgbearing inclusions ( Figure 4C).

Frontiers in Earth Science
frontiersin.org

Nanoscale characterization
The nanoscale study was conducted on a FIB-prepared TEM foil ( Figure 5A) extracted from across the core-rim boundary within this grain ( Figure 4C). Low-resolution HAADF STEM imaging and STEM EDS mapping of the foil ( Figure 5B) reveal cassiterite particles that increase in density away from the boundary, as well as a smaller population of darker, often acicular Mg-bearing silicate inclusions with contrasting density trend relative to the same boundary. Based on composition, some of these silicate inclusions are attributable to varieties of serpentine, albeit undifferentiated. Frontiers in Earth Science frontiersin.org

Inclusions in magnetite
Cassiterite, the most abundant included phase, varies in size from nanoparticles (<100 nm in diameter; NP) to fine particles (100-2,500 nm) and with variable distribution throughout the magnetite. Many such cassiterite inclusions form binary associations with a Mg-silicate as seen from higher magnification images and element maps ( Figure 6A). Although the silicate is typically far smaller than the cassiterite, their mutual boundaries suggest co-crystallization ( Figure 6B). STEM EDS element maps show measurable amounts of other elements in the cassiterite, including Ti and W (as shown on the element maps), but also U and Pb. EDS STEM integrated spectra yield the empirical formula (Sn 0.997 Ti 0.006 W 0.005 Pb 0.007 U 0.006 ) 1.021 O 2 for cassiterite (Table 2); inclusions with acicular habit ( Figure 6A) are the richest in trace elements. High-resolution images of cassiterite viewed on two main zone axes ( Figures 7A-D) show a very good fit with the crystal models and with STEM simulations (Figures 7E,F). We note that cassiterite displays epitaxial relationships with magnetite, in which <110*> cassiterite is parallel to <111*> magnetite , indicating preferential growth of Sn-bearing inclusions in the host Fe-oxide.
Fluorine was noted in other two Mg-bearing inclusions with a granular habit (Table 2, Figures 8,9). These correspond to the Mg-fluoride sellaite (MgF 2 ), and a coexisting silicate distinct from serpentine. Calculated compositions for the silicate from STEM EDS data yielded the empirical formula (Mg 5.45 Fe 0.26 ) 5.71 (Si 2.14 O 4 ) 2 F 1.16 , indicating a species belonging to the humite family (Gibbs et al., 1970;Ottolini et al., 2000;Friedrich et al., 2001) with a best match to chondrodite, ideally (Mg, Fe) 5 (SiO 4 ) 2 F 2 . The cation excess and F deficit can be attributed, in part, to difficulties in accurate measurement of F concentrations due to matrix effects in humite-group minerals as discussed for several analytical techniques, including EPMA (Ottolini et al., 2000).   The chondrodite grain analyzed by STEM EDS (maps in Figure 8A, data in Table 2) was suitable for imaging only if lowdose HAADF and iDPC STEM techniques were used ( Figures  8B,C). The chondrodite structure displays two different motifs when imaged on the [110] zone axis, with a relatively good correlation between the two types of images. The HAADF image ( Figure 8B) shows better the motif that occurs as bright, waving strips of atoms with~7.5 Å periodicity along the c axis. The iDPC

Frontiers in Earth Science
frontiersin.org 07 image ( Figure 8C), in contrast, shows a lighter, squarish, motif placed between the brighter strips. STEM simulation shows that the lighter motif on the iDPC image is attributable to the SiO 2 tetrahedra, flanked by Mg (Fe) and F cations ( Figure 8D). The waving strips can be related to the overlapping Mg and F cations at the corners of the SiO 2 tetrahedra. Although the match between the images and models is imperfect, the fast Fourier transform (FFT) pattern and electron diffraction simulation (Figures 8B,D) are in excellent agreement with the identification of this mineral from the chemical composition (Table 2).
Likewise, the identification of sellaite, the second F-bearing species, was made using STEM EDS analysis and high-resolution imaging of the stubby grains preserved in the upper part of the foil (Figure 9). Calculated compositions from integrated spectra give the empirical formula (Mg 0.98 Fe 0.02 )F 1.91 (Table 2). Element mapping ( Figure 9A) shows variable content of F and Mg, mostly due to variation in the thickness and exposure of grains rather than the co-existence of different species. As in the case of chondrodite, we note that the sellaite grains are separated from the cassiterite and serpentine inclusions. Identification of sellaite is confirmed by high-resolution iDPC imaging on [100], matching the crystal structure model and STEM simulation on this plane of view ( Figures 9B-E).

Mg-Si-bearing defects in magnetite
One of the most outstanding nanoscale features in the studied magnetite is the presence of dense fields of parallel, linear defects throughout the Sn-bearing, inclusion-rich domains, and to a lesser extent, within the margin of the foil (Figure 10). Such defects are clearly crosscut by needles of

Frontiers in Earth Science
frontiersin.org serpentine attached to grains of cassiterite ( Figure 10A). In closer detail, the linear defects observed on HAADF images consist of either: 1) arrays of hexagonal-shaped nanoparticles (~5 nm in width) that appear darker than magnetite; or 2) thin, still darker lines ( Figure 10B). Defect-free nanodomains are surrounded by cassiterite and serpentine inclusions ( Figure 10C), and magnetite from the cassiterite-free margin of the foil is densely populated by serpentine inclusions ( Figure 10D). The two types of defects are best displayed on that [110] zone axis in magnetite and may, or may not, overlap with one another ( Figure 10E). Fast Fourier  Figure 10F). In contrast, FFT patterns obtained from areas with the hexagonal-shaped darker NPs show no modifications from ideal magnetite ( Figures 10G,H). Loopshaped defects and hexagonal to rounded,~25 nm-diameter inclusions of serpentine located along the core-rim boundary in magnetite ( Figures 10I,J) resemble the defects within the Sn-rich part of the same grain, but the smallest serpentine inclusions are five times larger than those forming parallel arrays along <110*> directions in magnetite ( Figures 10E,J). Although the serpentine is not sufficiently stable under the electron beam to allow highresolution imaging, FFT patterns indicate that it is a monoclinic variety of antigorite ( Figure 10K; Dodony et al., 2002;Capitani and Mellini, 2004).
Higher magnification STEM EDS maps of areas with <110*> defects within magnetite show the presence of Si and Mg (Figure 11). The same two elements are measurable along both the lines and the hexagonal-shaped blebs ( Figures 11A,C). Profiles across a set of such linear defects indicate a positive correlation between Si and Mg ( Figure 11D). An antithetic distribution between Si and Mg with Fe is shown on element maps for a single line defect ( Figure 12A).
High-resolution imaging shows changes in intensity and typical arrangement of atoms on [110] magnetite ( Figures  12B-D). Interpretation of the HAADF STEM images illustrating this type of defect is made considering the atomfill crystal structure model for Fe atoms in [110] magnetite  Frontiers in Earth Science frontiersin.org defect of~6 Å width, representing 2d 110 in magnetite, can be attributed to partial or total replacement of Fe by Mg and Si ( Figure 12C). Using the same atomic model, the atomic arrangement across the defect is represented by the sequence: 3/2123/23′2/3212/3, whereby 3′ is a 3-type array with ½ (110) gaps in 2 M atoms (vacancies?) ( Figure 12D). The same sequence can also be considered as produced by twinning along the 3-type arrays with ½(110) shifts. Analogous to this are the (110) twins along the 2-type arrays, leading to the sequence 2,123/2/3212/, without chemical changes, which have been reported from silician magnetite at Olympic Dam, South Australia (Cook et al., 2022b).

Magnetite records early-stage mineralization-cassiterite associated with formation of humite skarn
The core-rim zonation patterns in magnetite (Figures 3, 4) could have formed during a single growth episode if we consider the micron-scale zonation patterns showing correlations between Sn, Mg, and Si with trends showing a gradual decrease from grain core to rim (Figure 4).
The magnetite grain core studied here shows up to 1.4 wt% SnO 2 (Table 1) which is within the same range (1-2 wt% Sn) that we have measured in other grains from serpentinite at Dulong (Figure 2). Although a little lower, these values are still within the same order of magnitude of Sn concentration (up to 2.82 wt% for 'magnetite 2') as those reported for low-temperature, Si-bearing magnetite with cores similarly rich in cassiterite inclusions from the Haobugao Fe-Zn polymetallic skarn deposit (Great Xing'an Range, China; Wang et al., 2018). At Haobugao, cassiterite was considered to have formed via exsolution from the late magnetite, whereas an earlier, high temperature, Sn-poor magnetite without µm-scale cassiterite inclusions was interpreted to contain this element in solid solution based on various geochemical trends. Such Sn-enrichment in magnetite is however interpreted by Wang et al. (2018) in the context of a calcic skarn, with different generations of magnetite, and thus distinct to those known from the Dulong skarn (Xu et al., 2021a). Formation of 2-5 µm-sized cassiterite inclusions in skarn magnetite from the Lost River tin mine, Seward Peninsula (Alaska) was also Frontiers in Earth Science frontiersin.org attributed to exsolution from an initial solid solution based on observed epitactic/crystallographic orientations between inclusions and host oxide (Desborough and Sainsbury, 1970). Our study shows the boundary between the core and rim in magnetite preserves evidence of multi-stage formation, rather than a single growth event. One important discovery, only revealed through a nanoscale approach, is the preservation of a Mg-and F-bearing stage (formation of sellaite and chondrodite; Figures 8, 9) since humite group minerals are key indicators of magnesium skarn predating formation of (pyroxene-garnet) calcic skarn at Dulong (Xu et al., 2021a).
This humite skarn is initiated by Mg-Si metasomatism forming fields of defects within magnetite and culminates with the discharge of volatile-rich, Sn-bearing fluids leading to formation of cassiterite inclusions within magnetite (Figures  4-6). The co-crystallization of cassiterite with serpentine from such fluids is supported by two-phase inclusions in apparent equilibrium ( Figure 6B). Collectively, the defects and inclusion associations within the core are attributable to a hightemperature, prograde skarn stage. The inclusion-free rim is, however, interpreted as the result of a superimposed, retrograde stage, as evidenced from orientation of serpentine inclusions but without cassiterite underlining the scalloped contact between the magnetite core and rim ( Figure 5). Retrograde stages are inherent to skarn evolution and are manifested by fluids being drawn backwards at the end of the prograde stage (Meinert, 1992;Dipple and Gerdes, 1998). Such an interpretation can explain the lack of cassiterite inclusions in the grain rims ( Figures 3A,B) as Sn is dissolved by these waning fluids and reprecipitated as coarser cassiterite outside the magnetite. Likewise, a small-scale overprint is recognized from the nanodomains that are free of Mg-Si defects ( Figure 10C). At this stage advanced  Dodony et al. (2002).

Frontiers in Earth Science
frontiersin.org serpentinization of the humite skarn is likely to occur as volatile extraction always enhances the fluid flux (e.g., Dipple and Gerdes, 1998).

Exsolutions from solid solution: Mg-Si defects predating Sn-enrichment in magnetite
The change in crystallographic directions, from defects along <110*> to epitaxial inclusion nucleation along <111*> (Figures 7A,C; 10e, f) indicates the structural control during magnetite growth and interaction with inbound fluids.
Rare spinel phases such as ringwoodite (SiMg 2 O 4 ) and ahrensite (SiFe 2 O 4 ), which have proven useful in petrological studies of the Earth's mantle, are considered high pressure polymorphs of minerals from the olivine group formed in the transitional zone at depths of 400-650 km (Biagioni and Pasero, 2014, and references therein). Such scenarios are clearly unrealistic to explain the formation of the Mg-Si defects shown here (Figures 11, 12). These could be attributed to diffusion-assisted exsolution processes, thus accounting for the spotted appearance (i.e., arrays of hexagonal-shaped Mg-Si-NPs) followed by fluid infiltration along (110) planes acting as twin dislocations in magnetite. They are likely to be more common than presently known in skarn magnetite.   -Gómez et al., 2006). In the titanomagnetite system, Sn incorporation is promoted at higher temperatures via vacancy-mediated long-range substitution of Sn ions within boundaries between magnetic domains, albeit with Sn 4+ -Fe 2+ pairs showing lower rates of diffusion than Ti 4+ -Fe 2+ pairs in magnetite (Hernández-Gómez et al., 2001, 2006). An example of cassiterite particles attached to ilmenite lamellae has been shown for magmatic titanomagnetite (Wang et al., 2012).
The present study shows no unequivocal evidence for cassiterite formation from exsolution processes, despite the preferential orientation between cassiterite and serpentine inclusions along <111*> directions in host magnetite. Such orientations are rather controlled by crystal structural parameters matching each other along these directions in magnetite ( Figures 7A,C, 10J,K).

Inclusions in magnetite: A record of evolving fluids and other petrogenetic implications
Studies of melt inclusions in quartz and fluid inclusions in quartz and cassiterite from topaz-albite granite from Zinnwald (Germany) have shown that F-rich hydrous melt fractions and hydrothermal fluids can form during differentiation of highly evolved Sn-granite magmas (Thomas et al., 2005). This demonstrates that F is extremely efficient in extraction and transport of Sn via hydrothermal fluids formed from a F-H 2 O rich granitic magma, even if the Laojunshan granite is recognized as a highly fractionated, metaluminous S-type granite (Xu et al., 2015).
Such an interpretation is consistent with the presence of grossular and vesuvianite in skarns at Dulong (Xu et al., 2021a), typical of reduced assemblages in Sn-W skarn deposits (Bowman, 1998) affiliated to S-type granites (Meinert, 1997). This implies that the initial fluids, although F-rich, are not the main contributor to Sn mineralization. Magnetite within marble or calcic skarn is not a major host for Sn ( Figure 2). Cassiterite from the calcic skarn and marble (Xu et al., 2021a) is thus considered to represent a later, geochemically distinct generation from the cassiterite occurring as inclusions in magnetite.
The structure and composition of magnetite is adept at preserving evidence of inherited protoliths, stages of alteration, and superimposed events, information that is critical for genetic interpretation in terranes with intensively altered lithologies. Expanding on the information obtained from traditional micronscale techniques, an extraordinary level of detail can be accessed using HAADF STEM imaging (Xu et al., 2014;Ciobanu et al., 2019Ciobanu et al., , 2022Huang and Beaudoin, 2021;Verdugo-Ihl et al., 2021), a technique that provides a more direct visualization of atoms in the crystal structures than traditional TEM methods.

Conclusion and implications
In summary, the magnetite studied here preserves a range of textures that reveal initiation of the skarn system as a humite, Overlays in (D) show interpretation in terms of partial to total replacement of Fe by Mg and Si. Atomic arrangement across the defect in terms of the sequence: 3/2123/23′2/3212/3, whereby 3′ is a 3-type array with ½ (110) gaps in 2 M atoms (v represents vacancies?).

Frontiers in Earth Science
frontiersin.org magnesium skarn stage (chondrodite). Magnesium-Si defects form along (110) planes prior to Sn-enrichment. The release of high-volatile, F-rich fluids can lead to precipitation of cassiterite inclusions along <111*> directions in magnetite. This is an alternative mechanism to "exsolution from solid solution" showing epitaxial growth between cassiterite and host magnetite. In its entirety, the magnetite records the history of metasomatism during prograde core growth and retrograde rim replacement. The Dulong magnetite is the first case in which (110) Mg-Si defects are documented. This is likely to be a common feature of many skarn systems since these include Mg-bearing (calc)silicates. Such characteristics can be useful to fingerprint initiation of metasomatic processes at the atomicscale, particularly important during changes in redox conditions at the switch between prograde and retrograde skarn regimes which also show enrichment in elements such as W and Mo, commonly associated with Sn. Magnetite that coexists with such regime changes, e.g., from more oxidizing diopside + andradite to reducing hedenbergite + grossular assemblages in calcic skarn, may record comparable complex defects and inclusion associations.

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions
JX and CC conceptualized this work with advice and guidance from NC and KE. FIB-SEM work was undertaken by JX and LC-D, EPMA analysis by JX and BW. STEM analysis and interpretation was performed by CC, AS and JX. The manuscript was written and revised by JX, CC and NC with contributions to figure preparation and editing of the manuscript by all other authors.

Funding
This research was funded by Scientific Funds of Fuzhou University, the National Natural Science Foundation of China (42072094), and Key Laboratory of Mineral Resources (KLMR 2017-13). Additional support was made available from the project "Trace elements in Fe-oxides: deportment, distribution and application in ore genesis, geochronology, exploration and mineral processing," supported by BHP Olympic Dam and the South Australian Government Mining and Petroleum Services Centre of Excellence.