The origin and composition of carbonatite-derived carbonate-bearing fluorapatite deposits

Carbonate-bearing fluorapatite rocks occur at over 30 globally distributed carbonatite complexes and represent a substantial potential supply of phosphorus for the fertiliser industry. However, the process(es) involved in forming carbonate-bearing fluorapatite at some carbonatites remain equivocal, with both hydrothermal and weathering mechanisms inferred. In this contribution, we compare the paragenesis and trace element contents of carbonate-bearing fluorapatite rocks from the Kovdor, Sokli, Bukusu, Catalão I and Glenover carbonatites in order to further understand their origin, as well as to comment upon the concentration of elements that may be deleterious to fertiliser production. The paragenesis of apatite from each deposit is broadly equivalent, comprising residual magmatic grains overgrown by several different stages of carbonate-bearing fluorapatite. The first forms epitactic overgrowths on residual magmatic grains, followed by the formation of massive apatite which, in turn, is cross-cut by late euhedral and colloform apatite generations. Compositionally, the paragenetic sequence corresponds to a substantial decrease in the concentration of rare earth elements (REE), Sr, Na and Th, with an increase in U and Cd. The carbonate-bearing fluorapatite exhibits a negative Ce anomaly, attributed to oxic conditions in a surficial environment and, in combination with the textural and compositional commonality, supports a weathering origin for these rocks. Carbonate-bearing fluorapatite has Th contents which are several orders of magnitude lower than magmatic apatite grains, potentially making such apatite a more environmentally attractive feedstock for the fertiliser industry. Uranium and cadmium contents are higher in carbonate-bearing fluorapatite than magmatic carbonatite apatite, but are much lower than most marine phosphorites.


Introduction
Phosphorus is essential to all life and is a limiting nutrient for plant growth (Cordell and White 2014). It is one of three macro-nutrients in fertiliser, but is a non-renewable resource principally mined from phosphorites, carbonatites and alkaline igneous rocks (Northolt et al. 1989). As of 2015, carbonatites and alkaline rocks account for~10% of global phosphate production and make up~5% of the world's phosphorus reserves (Pufahl and Groat 2017). Extraction from igneous sources is predominantly from large, low-grade (5-15% P 2 O 5 ) but easily beneficiated intrusions (e.g. Siilinjärvi, Finland; Khibiny, Russia), where apatite (Ca 5 [PO 4 ] 3 [OH,Cl,F]) crystallised as a magmatic phase. However, non-magmatically derived carbonate-bearing fluorapatite rocks can also occur in association with carbonatite and these rocks are actively mined for their phosphate content at, for example, Catalão I (Brazil), Editorial handling: R. Linnen Electronic supplementary material The online version of this article (https://doi.org/10.1007/s00126-020-01010-7) contains supplementary material, which is available to authorized users. and have been mined out at the Cargill (Canada) and Glenover (South Africa) carbonatites.
Carbonate-bearing fluorapatite rocks (also termed 'staffelite rocks' or 'francolite rocks') are known from dozens of carbonatite complexes (Fig. 1). They are principally composed of massive carbonate-bearing fluorapatite, which accommodates the carbonate ion through substitution into the apatite Cl/F/OH site (type A substitution) or phosphate site (type B;Fleet 2015). In comparison with igneous apatite deposits, these rocks can reach very high grades. Martison, Canada, for example, has an indicated resource of 54.3 Mt at 23.4% P 2 O 5 (Horner et al. 2016), compared with igneous apatite mined at Siilinjärvi with a grade of 3-6% P 2 O 5 (Puustinen and Kauppinen 1989). In addition, in areas where phosphorites are scarce and commercial fertiliser is expensive, such as sub-Saharan Africa, direct application of local carbonate-bearing fluorapatite can deliver increased crop yields (Appleton 1994(Appleton , 2002Van Straaten 2002;Sanchez 2002). In these cases, substitution of carbonate into the apatite structure is particularly beneficial as the reactivity and agronomic potential of apatite dramatically increases with higher carbonate content (Lehr and McClellan 1972).
The competing hypotheses outlined above invoke contrasting mineralisation environments, but the geological evidence supporting each model is equivocal. For instance, brecciation, which is commonly associated with carbonate-bearing fluorapatite deposits, can be interpreted as a hightemperature hydraulic fracturing/explosive degassing process (e.g. at Kovdor : Kapustin 1983;Krasnova 2003), or lowtemperature karstic collapse (e.g. Glenover; Verwoerd 1966). Commonly, the genesis of a deposit is determined by comparison with better-understood localities. For example, deposits considered to have formed by palaeo-tropical weathering processes, which are currently located in arctic climates (e.g. Sokli, Finland, and Martison, Canada;Vartiainen 1989;Sage 1991), bear similarity to deposits formed through modern tropical weathering in Brazil (e.g. at Catalão I, Toledo et al. 2004). However, extensive petrographic and geochemical analyses of many of these deposits are not commonly reported in detail, and compositional data for carbonate-bearing fluorapatite is typically limited to major elements from electron probe microanalysis (EPMA) and lacks CO 3 contents. This scarcity of geochemical data limits the robustness of inter-deposit comparison as means of interpreting rock genesis.

Carbonate-bearing fluorapatite occurrences in carbonatites
Occurrences of carbonate-bearing fluorapatite in carbonatites are often unpublished or are presented in publications which may be of limited circulation or difficult to obtain. Moreover, a myriad of terms are used for carbonate-bearing fluorapatite rocks from carbonatites, including 'carbonate-fluorapatite', 'staffelite', 'dahlite', 'collophane' or 'francolite' rocks, as well as 'weathered carbonatite' and (erroneously) 'phosphorite', making direct comparisons challenging. Several of these terms are widely used but are ambiguous, poorly defined or are prefixed by minerals which are discredited by the International Mineralogical Association. For example, the most commonly used term, francolite, is defined as a carbonate-bearing fluorapatite with greater than 1% fluorine and 'an appreciable amount of CO 2 ' (Sandell et al. 1939). However, McClellan and Van Kauwenbergh (1990) highlight that while most carbonate-bearing fluorapatite contains high fluorine contents (>3.77%), fluorine contents of <1% F can also occur in such apatite. Moreover, the term is commonly applied even where compositional data are not available, based on historical interpretations or textural data. An alternative term, 'carbonate-fluorapatite', has been discredited by the International Mineralogical Association as carbonate substitution is typically a subordinate component and 'carbonate-fluorapatite' lacks a distinct XRD pattern (Pasero et al. 2010). The term 'carbonate-bearing fluorapatite' is recommended instead and is used herein. However, we note that magmatic apatite from carbonatites, which is not the focus of this study, also contains a small amount of structurally bound CO 3 (Prins 1973). Consequently, the term 'carbonate-bearing fluorapatite' is preferentially used in this manuscript to describe fluorapatite with a measurable amount of carbonate substitution formed after the cessation of magmatic crystallisation. Where an alternative term has previously been used to name a local rock type (e.g. 'apatite francolite regolith'), this is also included for consistency. The challenges outlined above mean that detailed accounts of carbonate-bearing fluorapatite rocks are uncommon and can be difficult to obtain. While it is not the intention of this contribution to provide a thorough review of all examples, Fig. 1 indicates where carbonate-bearing fluorapatite has been described at a carbonatite and includes key references for each locality. It is evident that carbonate-bearing fluorapatite rocks are globally distributed and are not solely localised to areas of modern tropical weathering.
In carbonatites, carbonate-bearing fluorapatite typically forms the matrix of a breccia cementing clasts of altered carbonatite, fenite or fragments of other locally derived alkaline rocks and minerals. It occurs in two geological environments: (1) as a regolith covering the carbonatite and (2) as localised breccia units within a carbonatite complex. Regolith-hosted examples are commonly interpreted as a weathering product of the underlying carbonatite host rock, and these are highlighted in Fig.  1. Such an interpretation is relatively straightforward, especially when the weathering is ongoing, as the weathering zone typically has a fixed depth, covers a large part of the complex, lacks associated carbonate minerals which have been removed through dissolution and contains residually enriched magmatic apatite grains and other minerals which are resistant to chemical weathering, such as pyrochlore. Many of these localities are also located in equatorial areas, where deep weathering profiles are common (Fig. 1 ). Where carbonate-bearing fluorapatite is localised to a breccia unit, its genesis is less clear and two interpretations are common: (a) karstic collapse-in a (palaeo)-weathering environment or (b) a hydrothermal vent. In both instances, mixed interpretations are common (Fig. 1).

Study sites and sample selection
Five, globally distributed, study sites were selected for this study in order to assess a representative suite of carbonatebearing fluorapatite deposits. Bukusu and Catalão I are representative of sites where modern weathering processes are active, and samples were acquired from the weathering profiles at each of these locations. At Sokli, Kovdor and Glenover, the role of modern weathering is less apparent and, as outlined above, their genesis is less clear.
All of the study sites have some form of associated economic activity. Sokli and Bukusu are mineral resources at different stages of development, the carbonate-bearing fluorapatite rock at Kovdor is presently stockpiled while igneous apatite is mined as a co-product from the Zheleznyi pit, Glenover had a mineral reserve but was mined out between 1950 and 1970, and phosphate has been mined by at least two different companies at Catalão I for the past 30 years.

Sokli, Finland
The Sokli carbonatite complex is located approximately 145 km NE of Kemijärvi and is part of the Devonian Kola Alkaline Province. It is~18 km 2 in size, is ellipsoidal and composed of concentric rings (towards the centre) of fenite, metacarbonatite and a core comprising phoscorite, calcite and dolomite carbonatites, cross-cut by late Sr-REE and Babearing veins (Vartiainen and Paarma 1979;Vartiainen 1980; O'Brien and Hyvönen 2015; Fig. 1). The complex is covered by 5-10 m of glacial till (Vartiainen and Paarma 1979), with the underlying bedrock heavily weathered to a typical depth of 5-30 m, but to a maximum of 107 m (Vartiainen 1989; Seppo Gehör, pers comm 2019).
Phosphate mineralisation is hosted in the weathered rock and principally occurs in two localities: Loitso, in the northwestern part of the complex, and Pierkuli, along the eastern contact between carbonatite and fenite (Vartiainen 1989;Vartiainen et al. 1990; Fig. 1). The phosphate mineralisation is divided into three main types termed: (1) apatite francolite regolith, (2) apatite silicate residue and (3) apatite residue above ring dykes (Vartiainen 1989). The apatite francolite regolith is the most extensive ore type, and the only one which contains substantial quantities of carbonate-bearing fluorapatite. The latter two types are smaller in extent and are composed of eluvial apatite grains from the dissolution of carbonate with little/no development of secondary minerals.
The apatite francolite regolith is unconsolidated, loose and heterogeneous (Fig. 2a) and comprises several subtypes, as defined by Vartiainen (1989), including hard, granular, soft, mica-rich, Fe-rich and Mn-rich varieties. These are principally defined on the abundance of residual magmatic minerals and neo-formed carbonate-bearing fluorapatite, goethite and Mn oxides. Hard and granular sub-types are cemented by carbonate-bearing fluorapatite, to different degrees. The samples analysed for this study fit the description of the hard apatite francolite regolith subtype (e.g. Fig. 2b) and were acquired from near to the test pit in the Loitso area (Fig. 1). These rocks are creamy-white in colour, hard, porous and predominantly composed of massive apatite, with no clear fragments of carbonatite visible on the hand-specimen scale.

Kovdor, Russia
Kovdor is located~50 km ESE of Sokli and is also part of the Kola Alkaline Province. It is~42 km 2 in size and principally composed of alkali-ultramafic rocks (olivinite, pyroxenite), surrounded by a fenite aureole. Carbonatite and phoscorite make up a much smaller component of the complex, limited to a concentrically zoned plug in the southwestern part of the intrusion, surrounded by calcite carbonatite dykes (Ivanyuk et al. 2002;Krasnova et al. 2004).
The francolite deposit at Kovdor is located to the south of the Zheleznyi iron ore mine, which currently exploits the phoscorite for magnetite, apatite and baddeleyite (Fig. 1). The carbonate-bearing fluorapatite mineralisation forms structurally complex vein-like bodies of approximately 15-300 m wide, occurring in a semi-circular zone of about 3.5 km in diameter. The structure of these carbonate-bearing fluorapatite units roughly matches the surrounding carbonatite bodies, suggesting an association between the two rock types (Lapin and Lyagushkin 2014). The carbonate-bearing fluorapatite units form in funnel-shaped depressions, typically pinching out at a depth of 70-100 m, but locally extending tõ 200 m. The contact between unaltered and carbonatebearing fluorapatite rocks is commonly sharp but disintegrated (Lapin and Lyagushkin 2014). Notably, some of the deposits are not connected with the surface (Krasnova 2003). Krasnova (2003) identifies three cross-cutting types of carbonate-bearing fluorapatite rocks at Kovdor: (1) loose rocks of carbonate-bearing fluorapatite and vermiculite, cross-cut by (2) massive brecciated rocks which are locally cross-cut by stringers of psilomelane, which in-turn is associated with (3) psilomelane-vermiculite 'brecciaform' rocks. Breccia clasts in all these instances comprise residual fragments of carbonatite and other country rocks (Krasnova 2003;Lapin and Lyagushkin 2014). Three samples of material from Kovdor were acquired for this study, all of which predominantly fit the description of the loose carbonate-bearing fluorapatite and vermiculite rocks of Krasnova (2003). They were sourced from approximately 3-4 m depth, from a site of excavation to the south of the Zheleznyi mine (Fig. 1). The samples are soft and friable and composed of residual mica and magnetite, cemented by a groundmass of massive apatite and cross-cut by colloform apatite veins (Fig. 2c).

Bukusu, Uganda
Bukusu is the largest Cenozoic carbonatite complex in south east Uganda, spanning 50 km 2 . It consists of rings of syenite, ijolite and pyroxenite, with veins of magnetite-apatitephlogopite rocks and carbonatite, centred on a central hill of vent agglomerate (Davies 1956;Baldock 1971; Fig. 1). Outcrop, however, is poor due to the formation of up to 45 m of laterite (Davies 1947(Davies , 1956. Phosphate rock occurs on Busumbu ridge, to the south of the central vent (Davies 1947(Davies , 1956 Fig. 1). These rocks form in the regolith above unweathered carbonatite and magnetite-apatite-phlogopite rocks. The phosphate rock is heterogeneous and consists of a mixture of hard phosphatic rock in a mass of soft limonitic material. A single sample of the hard phosphatic rock was investigated for this study, sourced from the British Geological Survey collections (Fig. 2d). Little information is available about the environment from which this sample was collected, except that it is a surface sample. The hard phosphatic rock is grey-brown, friable and concretionary. It is composed of spongey, silvery-white carbonate-bearing fluorapatite locally forming stellate or fibrous aggregates around cores of colourless apatite crystals and, to a lesser extent, mica or magnetite.

Glenover, South Africa
The Glenover carbonatite complex is located approximately 170 km NNW of Johannesburg. It is approximately 15 km 2 in size, ellipsoidal and principally composed of pyroxenite, with variable amounts of phlogopite, magnetite and apatite (Verwoerd 1967). Dolomite carbonatite is the predominant carbonatite type and occurs as dykes, irregular veins and sills which dip towards the centre of the complex. Calcite carbonatite is subordinate and occurs We ath eri ng ho riz on s Approximate location of BRR1 Fig. 2 Examples of carbonate-bearing fluorapatite rocks from the studied localities. a Excavated unconsolidated, loose, apatite francolite regolith, and b sampled 'hard apatite francolite regolith subtype' (terminology after Vartiainen 1989) from Sokli. c Carbonate-bearing fluorapatite at Kovdor, with cream-white colloform apatite forming on top of a groundmass of altered mica and iron oxides. d Example of a clast in breccia from Bukusu, predominantly composed of apatite mixed with Fe oxides. e Contact between ferruginous apatite-rich breccia and apatite-bearing pyroxenite at the Glenover open pit, interpreted as a karstic collapse structure (Verwoerd 1967;van der Walt et al. 2012). Profile height is approximately 20 m. f Sample taken from waste pit at Glenover, representative of the apatite-rich breccia. Despite the brown colouration, the rock is predominantly made up of apatite, with darker areas formed of massive apatite with Fe impurities and lighter areas comprising late euhedral apatite in a pseudo-colloform habit. g Profile of weathered carbonatite at Catalão I, showing (A) barren material / laterite, (B) saprolite and (C) unweathered rock. Bench height is 5 m. The location of sample BRR1 is indicated with a close-up in (h), demonstrating the apatite-rich (cream/white in colour) composition towards the edge of the complex. Minor, late-stage, REEbearing carbonatite stringer veins occur throughout the complex. Exposure is poor, and weathering is extensive, with karst-like dissolution cavities reported down to depths of approximately 30 m (Verwoerd 1967) and oxidation of pyrite in carbonatite occurring to depths of 150 m (van der Walt et al. 2012).
The Glenover phosphate deposit is located near the centre of the complex (Fig. 1), in an area of elevated radioactivity, and is composed of a ferruginous apatiterich breccia body, approximately 200 by 250 m in size. The breccia body is likely to be bowl-shaped, and shallow in depth, with a sharply defined contact with the carbonatite (Fig. 2e) occurring at a maximum depth of approximately 150-170 m (Verwoerd 1967(Verwoerd , 1986van der Walt et al. 2012). The apatite-rich breccia is a dark red heterogeneous rock consisting of iron oxides, haematite after magnetite and small apatite grains, termed 'consolidated ferruginous earth' by Verwoerd (1967;Fig. 2f). Clasts in the breccia include fine-grained pale-red/purple fragments, predominantly composed of apatite mixed with Fe oxides, cemented by very fine-grained colloform white apatite (Verwoerd 1967). The apatite-rich material is now mined out, but samples matching the above description were acquired from adjacent dumps. For comparison, a sample of weathered carbonatite was also collected from a drill core (HTAC_008).

Catalão I, Brazil
Catalão I is located about 10 km NE of the town of Catalão in Goiás state and is part of the Late Cretaceous Alto Paranaíba Alkaline Province (Gomes and Comin-Chiaramonti 2005). The intrusion is approximately 6 km in diameter and is composed of a core of dolomite carbonatite, nelsonite and phoscorite dykes, surrounded by pyroxenite, glimmerite and serpentinised peridotite (Cordeiro et al. 2010(Cordeiro et al. , 2011Fig. 1).
Phosphate mineralisation is hosted in both laterite and saprolite, concentrated above the carbonatite, and may range between 15 and 250 m in thickness (Morteani and Preinfalk 1996;Fig. 2g). Economic mineralisation is located in the north and central parts of the complex, with grades reaching 30% P 2 O 5 (Carvalho and Bressan 1989). The laterite is divided into different zones by different authors (cf. Morteani and Preinfalk 1996;Oliveira and Imbernon 1998;Toledo et al. 2004), but a simplified sequence from the surface down can be summarised as follows: (a) barren material / laterite, (b) saprolite and (c) fresh rock (Fig. 2g). The thickness and extent of these zones is very irregular. The top-most zone is considered overburden, which includes alluvium and soil, and apatite grades are low. The saprolite is somewhat heterogeneous, with lenses of mica-or haematite-rich material, 'silicicrete' and fragments of unweathered rock. The sample analysed in this study is from the saprolite (Fig. 2g) and is brilliant white in colour, being predominantly composed of massive apatite, with minor residual grains of mica (Fig. 2h). For comparison, a sample of weathered carbonatite was also collected from a drill core (BRR3).

Analytical techniques
Textural characterisation of the samples was undertaken using a standard petrographic microscope, CITL Mk3 and Mk5 cold-cathodoluminescence (CL) equipment, and an FEI Quanta 600 scanning electron microscope (SEM). The CL was typically operated with a beam current of~300-500 μA and a corresponding voltage of 8-12 kV.
Electron probe micro analyses (EPMA) were undertaken at Camborne School of Mines using a JEOL JXA 8200 instrument, operated with a 10 nA, 15 kV beam defocussed to a 5-μm spot. Peak counting times were 20 s for CaO, Cl, SO 3 , La 2 O 3 , Ce 2 O 3 , MnO and FeO; 30 s for P 2 O 5 , F and SrO; 40 s for SiO 2 and Y 2 O 3 ; and 60 s for Na 2 O, with background times half those of the peaks. The primary standards used for La, Ce, Nd and Y were synthetic REE glasses from the University of Edinburgh; Ca, P and F were standardised against MKII Durango fluorapatite, while all other elements were standardised using Astimex mineral standards. X-ray counts were converted to oxide concentrations using the inbuilt Jeol φρz correction programme. Astimex fluorapatite and inhouse Durango and Odegaarden apatite grains were used as secondary standards, with results within 2% of published data from these localities.
Trace elements were analysed, using LA-ICP-MS, at the University of Leeds and Camborne School of Mines. Ablation at Leeds was carried out using a Lambda Physik 193 nm ArF excimer laser with a 100 μm spot and a repetition rate of 5 Hz, maintaining a fluence of approximately 4 J cm −2 . This was coupled to an Agilent 7500c ICP-MS operating in reaction mode, with a flow rate of 2.5 ml min −1 H 2, to suppress 40 Ar + interference on 40 Ca. Ablation at Camborne School of Mines utilised a New Wave Research 213 nm Nd-YAG laser coupled to an Agilent 7700 ICP-MS. The laser was operated with a 100 μm spot and a repetition rate of 10 Hz, resulting in a fluence of 3-4 J cm −2 . In both cases, median Ca concentrations for the apatite generation analysed, obtained by EPMA, were used as the internal standard composition and NIST SRM 610 was used for calibration. SRM 612, and at Camborne School of Mines, in-house Durango and Odegaarden apatites were used as a secondary standards with the concentrations of all analysed elements falling within 10% of published values.
Infrared spectra were collected for polished surfaces using a ThermoNicolet i10 FTIR spectrometer fitted with a Ge-tipped attenuated total reflectance head, housed at the University of Bristol. In total, 128 scans were collected from 650 to 4000 cm −1 at a resolution of 4 cm −1 using a mercury cadmium telluride detector and a KBr beamsplitter in absorbance mode. The effective collection area ranged in size from approximately 30 to 10 μm squares.

Apatite textures
Textural details for apatite from each locality are presented in Figs. 3, 4, 5, 6 and 7, grouped by locality. In total, five different apatite habits are recognised. These habits can be summarised as follows: (1) ovoid grains, (2) epitactic overgrowths, (3) massive, (4) colloform and (5) late euhedral. Not all habits occur in apatite from each example, with apatite from Bukusu and Glenover notably not exhibiting epitactic overgrowths or colloform types. The different habits are interpreted as 'primary' (ovoid grains residually derived from a cognate carbonatite melt) and 'secondary' (overgrowths, massive, colloform and late euhedral), which could be derived from either hydrothermal or weathering processes. The term 'secondary' is not intended to convey information about the apatite genesis, but is used here to describe apatite forming after ovoid grains.
Ovoid grains are typically within the centre of 'clasts' and are locally associated with magnetite and biotite. Oscillatory zoning within the ovoid grains is evident in CL images, ranging in colour from yellowish to purple. Zoned, ovoid grains are ubiquitous in plutonic carbonatites and are considered to crystallise early in a carbonatite melt with rounding likely forming due to subsequent mechanical or chemical abrasion during emplacement (Chakhmouradian et al. 2017). Based on this common association, ovoid apatite is almost-certainly derived from a carbonatite melt. Zoning within the grains is not used here to differentiate the ovoid apatite into further sub-types. Grains from Glenover and Catalão I are partially fragmented and cemented by calcite and Fe oxide minerals, un-associated with carbonate-bearing apatite (Figs. 5a and 7b).   Fig. 3 CL (a, b) and BSE images (c-f) of different primary and secondary apatite generations from Sokli. a, b Primary apatite grains hosted in a massive apatite groundmass, cut by fractures cemented with colloform apatite and haematite. Note the muted blue luminescence in the massive apatite, and the lack of luminescence in the colloform apatite. c A primary apatite grain overgrown by euhedral zoned secondary apatite, hosted in a massive apatite groundmass. d, e Primary apatite grains with spongey secondary apatite overgrowths, possibly after euhedral zoned secondary apatite (d) is fractured and infilled with a darker generation of spongey secondary apatite. f Primary apatite with a small epitactic secondary apatite rim, overgrown by colloform apatite formed of 20-30 μm wide parallel apatite needles Epitactic overgrowths are distinct from zoning within the ovoid grains as indicated by a clear change from an ovoid habit to euhedral hexagonal overgrowths. This change is coupled with a decrease in BSE intensity, substantially reduced CL intensity, and a drop in RI compared with the original ovoid grain (Figs. 3c, 4b and 6b). Epitactic overgrowths typically exhibit fine euhedral oscillatory zoning in BSE images, with no corresponding CL variation apparent. Epitactic overgrowths are commonly overprinted by massive apatite, with the original overgrowth only partly preserved (Figs. 3c, f and 4b), or with just the crystal outline retained (Fig.  3d, e).
Massive apatite forms the groundmass of all the carbonate-bearing fluorapatite rocks studied. It grows around and cross-cuts earlier ovoid and euhedral grains (Figs. 3d and 6c). Locally, the massive apatite retains  (Fig. 4a, c). The massive structure is turbid in thin sections and composed of many fine-grained interlocking needles between 1 and 10 μm wide, but up to several hundred micrometres in length (Fig. 3f). Locally, the massive apatite is spongey, with spongey areas commonly corresponding to lower BSE intensity (Figs. 3c and 7c). The spongey texture is likely to be indicative of partial apatite dissolution (Putnis 2002(Putnis , 2009). In most samples, a locally precipitated reaction product is not evident but massive apatite from Glenover is also associated with Ca-rich monazite-(Ce) and cerianite-(Ce) (Fig. 7c). Samples from Bukusu and Sokli do not exhibit any luminescence in CL images, but massive apatite from Kovdor luminesces a muted blue colour while massive apatite from Glenover luminesces yellow.
Colloform apatite is the most distinctive megascopic apatite habit (Fig. 2). The samples in this study are cream-white in colour, but turquoise coloured samples have been reported (Chakhmouradian et al. 2017). In thin section, colloform apatite is first formed of late euhedral apatite needles radiating perpendicular to a growth surface (Fig. 3a, b). The needles vary in thickness, between 1 and 50 μm, although they most commonly are towards the lower end of this range. Euhedral apatite from Catalão I is considered here as part of the late euhedral stage, although it is subtly different from the other examples (Fig. 5b). For instance, original ovoid grains are not commonly found in the cores of the 'overgrowths' although their presence was demonstrated by Toledo et al. (2004) in their type 5 apatite. Additionally, while fine oscillatory zoning is not present, similar changes in the crystallisation environment are manifested by lamellar, or skeletal, growth. Euhedral apatite overgrowths at Catalão I also exhibit subtle colour variation in CL images.
Up to three generations of colloform apatite occur at the different sites and these vary in turbidity, needle thickness and colour. Early generations at Sokli are composed of thinner needles than later stages and are commonly darker in colour (Fig. 3f). In the final growth stages at Glenover, the needlelike form is lost, and the apatite is more equant, forming late euhedral apatite (Fig. 7d-f). The needle-like structure is also absent in later colloform generations at Kovdor, and apatite forms a blister-like texture (Fig. 4d).
Colloform apatite is commonly followed by the growth of REE-bearing crandallite group minerals as well as iron oxides (Figs. 3b and 4d). Additionally, at Glenover, cerianite and Ca-rich monazite also occur in pore spaces, locally associated with quartz, after late euhedral apatite. Luminescence is absent in colloform apatite with the exceptions of Kovdor, where it is muted blue to green in colour (Fig. 4a), and Glenover where the yellow luminescence present in the massive apatite is retained. However, yellow luminescence at Glenover is only found in spongey, altered apatite, and where

Compositional variation
FTIR spectra indicate similarities between carbonatebearing fluorapatite types from the studied localities. Spectra from all samples indicate the presence, to varying extents, of CO 3 2− in the apatite structure, as indicated by V 3 CO 3 peaks between 1400 and 1550 cm −1 (Fig. 8). Nonetheless, there is considerable variety between the magnitude and shape of the CO 3 peaks. Ovoid apatite grains exhibit the lowest magnitude V 3 CO 3 peaks, compared with other apatite types analysed. They are additionally characterised by a distinctive peak at approximately 720-730 cm −1 , and a small, but consistent, OH peak at~3530 cm −1 (Fig. 8). The 720-730 cm −1 peak is probably related to vibrational OH bands that are weakly hydrogen bonded to F ions (Fowler 1973). In contrast to the grains, colloform apatite does not exhibit peaks at 720-730 and~3530 cm −1 , but does show a broad stretch between 2800 and 3600 cm −1 , related to molecular (i.e. not structurally bound) H 2 O, as well as a peak at around 2900 cm −1 (Fig. 8). The V 3 CO 3 peaks are consistently higher than those in the ovoid grains. Euhedral overgrowths from Bukusu and massive apatite from Kovdor exhibit the same spectral features as colloform apatite (Fig. 8). In grains from Catalão I and Kovdor, as well as in massive apatite at Kovdor, a small peak is apparent at 2350 cm −1 , corresponding to molecular CO 2 , which also correlates with a broad H 2 O peak (2800-3600 cm −1 ). These peaks are probably due to fluid inclusions in the area analysed.
Apatite major and trace element compositions are presented in ESM Tables S1-4. The composition of the ovoid grains is typical of apatite derived from a carbonatite, with relatively high F, Sr, Na and REE contents, and low Cl and Mn, compared to apatite from other rock types (e.g. granitoids ;Hogarth 1989;Belousova et al. 2002;Broom-Fendley et al. 2017). All grains exhibit almost-straight REE patterns, with a negative Y anomaly, although grains from Glenover exhibit a small hump Fig. 7 CL (a, d, e) and BSE (b, c, e) image of primary and secondary apatite from Glenover. a-c Primary rounded grains, hosted in secondary apatite (a, c) and in a rock composed of secondary Fe oxides and calcite (b). Note the breakdown and replacement of primary apatite grains in a, c to secondary, yellow-/non-luminescent apatite and the presence of bright grains of Ca-rich monazite associated with the massive apatite. d Colloform apatite, transitioning into late euhedral apatite. e, f Late euhedral apatite, associated with small euhedral grains of Ca-rich monazite. Note the yellow luminescence corresponds with spongey apatite while nonspongey areas do not luminesce around Ce-Pr (Fig. 9). Thorium and uranium form a positive correlation, with grains from Glenover reaching the highest Th contents, up to~300 μg/g (Fig. 10a). The concentration of U in some Sokli grains is decoupled from the Th contents, possibly as a result of the ablation of inclusions or contamination from later surrounding apatite generations.
Secondary apatite has notably lower EPMA totals than the primary grains (ESM Table S2), with the lowest from epitactic overgrowths from Kovdor and Sokli, reaching~90%. The decreasing total corresponds to an under-full PO 4 site (Fig. 10b), and almost-certainly reflects the substitution of CO 3 into the structure, as supported by the presence of CO 3 peaks from the FTIR data (Fig. 8).
Minor and trace element concentrations are notably lower in secondary apatite than primary apatite grains, typically below the EPMA detection limit and, in some cases, the LA-ICP-MS detection limit too. Although only a limited number of analyses are available, epitactic overgrowths are the closest in trace element contents to primary grains, with similar REE, Na and Sr concentrations (Fig. 10c, ESM Table S4). REE, Na and Sr contents in massive and colloform apatite drop steadily towards concentrations below the LA-ICP-MS detection limit (Fig. 10c, d). However, some trace element contents are elevated relative to earlier ovoid apatite. A key example is U, which can reach contents of over 1000 μg/g and, in samples from Kovdor, exhibits a weak negative correlation with REE contents (R 2 = 0.7; Fig. 10e). Conversely, U and REE positively correlate in the euhedral overgrowths from Catalão I (R 2 = 0.75; Fig. 10e). Transition metals, such as V, Cu, Zn and Cd, also exhibit elevated contents in the secondary apatite generations, relative to ovoid apatite. Moreover, with the exception of Cu, these elements also exhibit a weak positive correlation with U (Figs. 10f-h), although data for Zn are limited to only a few analyses. Similarly, data for Mo and Cr are restricted to a small number of data points and are close to the detection limit, but also exhibit a tentative positive correlation with U.
The REE distribution of non-magmatic apatite from Sokli, Kovdor and Glenover emulates that of the magmatic apatite grains in each example, but with lower overall REE contents (Fig. 11). The REE patterns for Catalão I and Bukusu, however, are notably flat, contrasting with the LREE-enriched pattern of the primary apatite grains from these localities. Secondary apatite from Catalão I and Bukusu both exhibit a negative Ce anomaly, and a positive Y anomaly, which is most pronounced in samples from Catalão I. Massive and colloform apatite from Sokli and Kovdor generally also have a strong negative Ce anomaly, but some analyses exhibit a small positive anomaly, or none at all. In addition, at Sokli and Fig. 8 IR spectra of apatite from the study sites. Note the presence of V 3 CO 3 peaks in all samples, but the elevated size of those from secondary apatite. Primary apatite has a distinct OH peak at 3530 cm −1 and an OH-F peak around 720 cm −1 . Secondary apatite has a broad stretch between 2800 and 3600 cm −1 . Spectra are normalised to the maximum PO 4 intensity, and intensity units are arbitrary. PO 4 peaks are trimmed in order to emphasise detail in smaller peaks. Dashed spectra indicate possible mixed signal with fluid inclusion(s). Data in ESM Table S5 L a C e P r N d Kovdor a Y anomaly is typically absent, or slightly negative (Fig. 11). Compositional differences between different secondary apatite generations are limited to analyses from Kovdor, Sokli and Glenover. Massive and colloform secondary apatite from Sokli is compositionally identical, within the range of data for each apatite type (ESM Table S4). While there is also substantial overlap between different generations at Kovdor, some differences between the different apatite generations are apparent (Fig. 10). Massive apatite from Kovdor has the highest REE, Na and Sr contents, most reminiscent of ovoid apatite, compared to colloform apatite. Conversely, a later generation of colloform apatite has noticeably the lowest REE, Na and Sr contents, and the highest U and Cd concentration (Fig. 10g, h). There is no noticeable compositional change from the core to rim of colloform apatite, nor any substantial variation in composition with luminescence colour. At Glenover, massive apatite has a marked negative Ce anomaly which is absent in the late euhedral type (Fig. 11e).

Interpretation and discussion
Generalised paragenetic evolution Figure 12 is a generalised paragenetic sequence of the different apatite generations. Ovoid magmatic fluorapatite grains (i) are fragmented and overgrown by epitactic carbonate-bearing fluorapatite (ii) which is enveloped in massive, fine-grained carbonate-bearing fluorapatite (iii). These rocks, when subsequently fragmented, show additional late euhedral (iv) and colloform (v) generations of apatite occurring along fractures and in vugs. Multiple colloform generations may also occur. REE minerals (vi) typically crystallise last. Locally, magmatic grains exhibit resorption / dissolution (i-a), relating to changing magmatic-hydrothermal conditions during carbonatite emplacement. Dissolution and resorption of epitactic carbonate-bearing fluorapatite (ii-a) and massive carbonate-bearing fluorapatite (iii-a) is common in most samples.    Fig. 12 A generalised paragenetic sequence for carbonate-bearing fluorapatite deposits associated with carbonatites. Roman numerals denote the successive stages of apatite growth. The columns to the right of the paragenetic stages indicate deposits where the different stages have been identified. Black lines denote evidence from this study, while grey lines are inferred from other studies (Oliveira and Imbernon 1998;Tazaki et al. 1986Tazaki et al. , 1987Subasinghe 1988, 1989a, b;Potapoff 1989;Sage 1991;Araújo 2015). Numbers in square boxes refer to figures which show evidence for the interpreted stage. Numbers in circles indicate the equivalent stage as interpreted by previous researchers (Sokli : Vartiainen 1989;Kovdor: Krasnova 2003;Catalão I: Toledo et al. 2004;Mabounié: Boulingui 1997;Matongo: Decrée et al. 2016;Juquiá: Walter et al. 1995a) Previous workers (Vartiainen 1989;Krasnova 2003;Toledo et al. 2004) have recognised the different successive generations of secondary apatite at the Sokli, Kovdor and Catalão I deposits, and these are included in Fig. 12 for comparison. At Sokli, Vartiainen (1989) recognised three generations of carbonate-bearing fluorapatite. Stage 1 and stage 2 carbonate-bearing fluorapatites form a cement, with stage 1 comprising a microcrystalline aggregate with goethite and crandallite and stage 2 forming small (<0.02 mm) individual fibrous crystals of carbonate-bearing fluorapatite, which typically occur as radial growths. Stage 3 is milky-white and occurs in fissure cracks and cavities with abundant dust-like inclusions. At Kovdor, Krasnova (2003) divides carbonatebearing fluorapatite into four sub-types. These consist of (1) overgrowths on relict residual magmatic apatite grains, (2) fine-grained aggregates intermixed with Fe and Mn oxides, (3) crusts/kidney-like aggregates and (4) porcelaneous aggregates. At Catalão I, apatite in the laterite is present both as residual magmatic grains and as neoformed apatite and is subdivided into 7 different subtypes by Toledo et al. (2004). Residual grains (1) are commonly strongly corroded and fractured, with corrosion typically located along fracture planes. Locally, the residual grains are rimmed (2), replaced or cross-cut (3) by secondary apatite. In places, overgrowths form small euhedral tabular prismatic crystals (4) or larger, lamellar grains (5). Sub-types 6 and 7 grow in fissures or aggregates in weathered material. They consist of fibrous crystals in aggregates growing perpendicular to the growth surface (6) or globules, locally linked together in botryoidal growth (7). Other phosphate minerals in the laterite include gorceixite, crandallite, florencite and, in 'silicicrete', monazite (Oliveira and Imbernon 1998).
In general, the paragenetic sequence presented in Fig. 12 is coincident with that of Toledo et al. (2004). Nonetheless, some subtle differences include the exclusion of their stage 3 (intracrystalline bands of acicular secondary apatite across grains of primary apatite), which we consider to be coincident with the formation of massive apatite as part of our stage iii (Figs. 3d and 6c). We also only recognise their stage 5 (lamellar secondary apatite surrounding primary apatite) in our sample from Catalão I and tentatively associate this with our late euhedral stage iv.
Localised dissolution of carbonate-bearing fluorapatite is prominent in some samples, especially from Sokli. This has the effect of masking some important paragenetic stages, such as the formation of epitactic overgrowths (stage ii; Fig. 3c-e), which may explain the absence of this stage in previous work (Vartiainen 1989). With the recognition of the presence of epitactic overgrowths at Sokli, it is clear that the different apatite stages from Kovdor, Sokli and Catalão I are texturally near-identical.
While epitactic overgrowths at Bukusu and Glenover are not readily apparent, it is quite possible that these are either obscured by subsequent alteration or are missing as a consequence of the limited availability of sample material.
Some of the aspects of the paragenetic sequence presented in Fig. 12 have been previously noted by workers at several different deposits. In cases where work has focussed on different aspects of a deposit, information on the order of crystallisation is limited to noting the presence of residual apatite hosted in a secondary apatite groundmass (Cargill, Canada, [Sage 1991]; Ngualla, Tanzania, [Witt et al. 2019]; Araxá, Brazil [Mariano 1989]), the formation of secondary apatite on the rim of a primary apatite grain (Dorowa, Zimbabwe [Fernandes 1978]; Mt Weld, Australia [Lottermoser 1990]), or the formation of fibrous apatite in fractures (Anitapolis, Brazil, [Girard et al. 1993]). However, paragenetic descriptions of secondary apatite are available for the Mabounié, Matongo and Juquiá complexes (Walter et al. 1995a, b;Boulingui 1997;Decrée et al. 2016), and sufficiently detailed descriptions of apatite textures are reported for the Eppawala, Salitre I and Martison deposits such that a crystallisation sequence can be inferred (Tazaki et al. 1986(Tazaki et al. , 1987Subasinghe 1988, 1989a, b;Potapoff 1989;Sage 1991;Araújo 2015). These localities have also been summarised in Fig. 12, with reference to ascribed paragenetic stages.
Notably, the textural evolution of other carbonatebearing fluorapatite deposits bears a strong similarity to those studied in this contribution. For instance, epitactic apatite, equivalent to our stage ii (Fig. 12), occurs at Mabounié (described by Boulingui 1997 as type 2 'overgrowing apatite'), Matongo (as 'prismatic apatite'; Decrée et al. 2016) and Juquiá (as 'type 2 apatite'; Walter et al. 1995a, b). 'Massive apatite', equivalent to stage iii, occurs in most comparison examples with the possible exception of Juquiá where it is not explicitly described. At Mabounié, massive apatite exhibits similar remnant cleavage after carbonate, as is present at Kovdor (Boulingui 1997). Commonly, this stage of apatite growth is also associated with the formation of iron oxyhydroxide minerals (Salitre I, Mabounié, Martison), or siderite (Matongo, Tomtor;Kravchenko and Pokrovsky 1995). The development of late euhedral apatite (type iv) is ubiquitous, although at some localities it is unclear if this stage post-dates massive mineralisation, or if massive mineralisation is a further manifestation of this mineralisation type, with no pore space remaining. The development of subsequent colloform apatite generations (type v) is quite variable. At Juquiá, for instance, Walter et al. (1995b) report several generations of colloform and late euhedral apatite (their apatite types 3-6), with later bands cross-cutting earlier colloform generations.
Similarly to the studied localities, the occurrence of different REE minerals late in the paragenesis, where reported, is quite variable. Neo-formed monazite (or rhabdophane) is commonly associated with carbonatebearing fluorapatite (Matongo, Salitre I), but does not always occur at the end of the paragenetic sequence (Decrée et al. 2016). The presence of crandallite group minerals, such as florencite, notably varies depending on weathering depth (Kravchenko and Pokrovsky 1995;Oliveira and Imbernon 1998;Araújo 2015).
While, as highlighted above, there are subtle local textural differences between the carbonate-bearing fluorapatite deposits in this study and those previously investigated, we consider the paragenetic sequence in Fig. 12 to likely be applicable to carbonatite-hosted apatite deposits in general. Considering the textural similarities alone, a common process for the formation of these rocks, as suggested by other authors (Lapin and Lyagushkin 2014), seems highly probable. As several of the examples listed here (Catalão I, Mabounié, Salitre I) unequivocally formed by weathering, by extension, we propose that such a mechanism can account for carbonatebearing fluorapatite deposits related to most carbonatites, even where modern weathering is not active (e.g. Glenover, Kovdor).

Compositional similarities and evolution
Compositionally, primary grains from each locality are identical to apatite in unaltered carbonatites (e.g. Kovdor, Zaitsev et al. 2015;Sokli, Hornig-Kjarsgaard 1998). While some grains texturally appear partially altered, their composition does not diverge from that of the primary grains (Fig. 10), suggesting that this small-scale alteration is within the compositional bounds of magmatic carbonatite-apatite.
Textural similarities between carbonate-bearing fluorapatite from the studied deposits are strongly borne out in the compositional data, especially those where samples are available for the full paragenetic sequence. The FTIR spectra and major element compositions for non-magmatic apatite from each locality are near-identical, with the exception of epitactic overgrowths ( Fig. 8; ESM Table S2), demonstrating the presence of CO 3 in the apatite structure, as inferred from low EPMA totals and noted in previous studies (Davies 1947). While the trace element contents vary somewhat between different deposits, many of the compositional trends are the same. For example, all secondary apatite varieties trend towards higher U, and lower REE, Na and Sr contents, than primary apatite from the same deposit (Fig. 10c-e). Importantly, this trend is broadly in harmony with the textural evolution of the apatite, with later paragenetic stages containing the highest U, and lowest REE, Na and Sr contents (Fig. 10c-e).

REE distribution and Ce anomalies
The shape of the REE distribution for colloform and massive carbonate-bearing fluorapatite from Sokli, Kovdor and Glenover broadly matches that of primary ovoid grains from these localities (Fig. 11a-c). In early stages from Kovdor, the distribution matches that of primary apatite from the same deposit, but with a negative Ce anomaly (Fig. 11a). In later generations, the REE concentration decreases linearly for each REE dropping off to below the LA-ICP-MS detection limit.
The lack of REE fractionation in later generations of carbonate-bearing fluorapatite at Sokli, Kovdor and Glenover is intriguing. For instance, it contrasts with the extreme REE fractionation observed in hydrothermal apatite from Tundulu, Kangankunde and Songwe (Malawi), which are enriched in the HREE relative to earlier ovoid igneous apatite from the same complexes (Broom-Fendley et al. 2016. Such a contrast circumstantially supports the textural evidence that the carbonate-bearing fluorapatite from the studied localities is not derived from a hydrothermal environment. However, in the Malawian examples, the change in REE distribution is considered to be caused by the differing stability of some REE anion complexes between the LREE and HREE (Migdisov and Williams-Jones 2014). Such a process is not unique to the hydrothermal environment and, depending on the REE complexes formed, can be absent in hydrothermally formed minerals. While it is beyond the scope of this study to comment on the composition of the REEtransporting agent in the studied examples, the lack of fractionation rules out a ligand which fractionates the REE, such as Cl (at elevated temperature; Migdisov et al. 2009), F, and, at ambient temperature, CO 3 Byrne 2000, 2004).
With few exceptions, all secondary apatite exhibits a negative Ce anomaly (Fig. 11), a feature also noted at the Eppawala, Juquiá and Kovdor complexes (Tazaki et al. 1987;Walter et al. 1995a, b;Chakhmouradian et al. 2017). The presence and magnitude of the Ce anomaly does not change with different paragenetic stages.
Anomalous Ce contents, relative to the other REE, are caused by the transition of Ce 3+ to Ce 4+ in highly oxidizing and alkaline pH conditions and, consequently, the different geochemical behaviour of Ce 4+ compared with the other REE (Henderson 1984;Akagi and Masuda 1998). Such conditions are typically met during rock weathering and low-temperature surface processes (Braun et al. 1990(Braun et al. , 1993Leybourne and Johannesson 2008;Janots et al. 2015;Sanematsu and Watanabe 2016). For example, during the weathering of granitic rocks, acidic fluids, sourced from the breakdown of organic matter in the soil horizon, transport the REE as aqueous complexes down the weathering horizon. This process splits the weathering profile of a granite into an upper zone of REE depletion and a deeper zone of REE accumulation, where the REE precipitate due to an increase in pH caused by a change in mineralogy or the presence of groundwater (Sanematsu and Watanabe 2016). Under most ambient conditions, however, Ce 4+ is less soluble than other REE 3+ ions, resulting in Ce 4+ preferentially accumulating in the upper REE depletion zone as a solid phase, such as cerianite-(Ce), or adsorbing onto Mn oxides, Fe oxides or clay minerals (Bau 1999;Takahashi et al. 2000;Bau and Koschinsky 2009). These solid phases, therefore, exhibit a positive Ce anomaly, while redeposited minerals in the REE accumulation zone, lower in the profile, exhibit a negative Ce anomaly. The negative Ce anomaly in the secondary carbonatebearing fluorapatite is, therefore, consistent with formation from a fluid where Ce 4+ has already partitioned into Ce-rich minerals higher in the weathering horizon. Elevated Ce 4+ contents in the upper part of the weathering horizon have been indicated by the presence of CeO 2 -rich pyrochlore from the upper weathering levels at Mt Weld (Lottermoser and England 1988) and Lueshe (Wall et al. 1996). However, the processes of weathering and REE redistribution above a relatively P 2 O 5 -poor granite are most likely not representative of those above a P 2 O 5 -rich carbonatite. For instance, previous studies of the REE profiles of weathered carbonatites do not show Ce anomalies in the REE minerals or bulk rock from the upper levels of a carbonatite weathering zone (e.g . Lottermoser 1990;Walter et al. 1995b;Morteani and Preinfalk 1996;Oliveira and Imbernon 1998;Broom-Fendley et al. 2020). Furthermore, the remarkably low REE concentration of carbonate-bearing fluorapatite (ESM Table S4) contradicts the notion that it is forming in a REE accumulation zone, as apatite typically readily incorporates the REE. Moreover, at Kovdor, Bukusu, Glenover and Catalão I, Fe oxide and Ce 4+ -bearing minerals commonly occur after the formation of the final carbonate-bearing fluorapatite stage. There is, thus, no immediately obvious phase which is likely to retain and fractionate Ce 4+ early in the genesis of these rocks, and, therefore, another mechanism to remove Ce 4+ may be in operation.
While uncommon, negative Ce anomalies can occur in residual solid phases, and preferential transport of Ce 4+ over the trivalent REEs has been described in certain conditions. A key example of a rock with a negative Ce anomaly is a marine phosphorite (McArthur and Walsh 1985). However, this is inherited from the REE pattern of seawater which, in turn, is derived from the retention of Ce 4+ in weathered solids, akin to the mechanism discounted above. As an alternative example, organic ligands, especially siderophores, have been shown experimentally and in natural studies to preferentially complex with Ce 4+ over the REE 3+ , and retain Ce 4+ in solution (Pourret et al. 2007;Tanaka et al. 2010;Loges et al. 2012;Bau et al. 2013;Kraemer et al. 2015Kraemer et al. , 2017. Thus, preferential transport of Ce 4+ could result in a negative Ce anomaly in the residual solid phase and may explain the negative anomaly in the carbonate-bearing fluorapatite deposits. Alternatively, anomalously positive Ce anomalies have also been documented in highly alkaline waters where Ce 4+ forms stable complexes with carbonate (Möller and Bau 1993;Johannesson et al. 1994). Given the extreme dissolution and replacement of carbonate in the secondary apatite-rich rocks, high carbonate concentrations in the REE-transporting fluids are feasible and thus carbonate complexation may also explain the negative Ce anomalies in the residual secondary apatite. The presence of a Ce anomaly in these rocks is, therefore, not unequivocal evidence for formation in a weathering environment, as other processes can lead to Ce fractionation. However, the ubiquity of the Ce anomaly in the secondary carbonate-bearing fluorapatite does point towards formation in an oxic environment, and in most circumstances, this requires the presence of organic ligands which can only occur in the surface environment.

Implications of elevated U contents in secondary carbonate-bearing fluorapatite
In addition to the presence of a Ce anomaly, the increasing concentration of U in the secondary carbonate-bearing fluorapatite also points towards a strong redox control on the formation of this mineral. Oxidised U 6+ is relatively soluble, compared to U 4+ , and can be transported as phosphate or carbonate complexes in neutral and alkaline solutions, respectively (Langmuir 1978). The transport and reprecipitation of uranium in weathering profiles is well documented, with uranium commonly adsorbing onto Fe oxide phases or reprecipitating as uranium-bearing phosphate minerals in the weathering horizon (Murakami et al. 1997;Jerden and Sinha 2006;Göb et al. 2013;Boekhout et al. 2015). Extensive experimental and natural studies have shown that hydroxylapatite and fluorapatite can be highly effective at immobilising aqueous U from groundwater (Fuller et al. 2002;Jerden and Sinha 2003;Ohnuki et al. 2004;Simon et al. 2008;Lammers et al. 2017). In these cases with high U contents in the groundwater, a dissolution reprecipitation process leads to the precipitation of secondary uranium phosphate minerals. However, for solutions with lower U (< 7000 ppm) and elevated Ca contents, which are more appropriate for a weathering carbonatite, the precipitation of uranium phosphates is suppressed, and instead U 6+ adsorbs onto the surface of apatite (Fuller et al. 2002;Simon et al. 2008). Importantly, adsorption is partly reversible, such that successive batches of weathering fluids would be capable of removing U 6+ from earlier apatite and re-fixing it onto later generations of apatite. With the continued input of newly weathered U from the breakdown of primary minerals from the host carbonatite, this process would progressively increase the U content of later apatite generations.

Deleterious and beneficial element contents
The content of U, Th and heavy metals, such as Cd, is an important criterion in the economic viability of potential phosphate sources for fertiliser as these elements are typically passed into the soil and, consequently, into the food chain (Roberts 2014). Regulation of the deleterious element content of fertiliser is becoming increasingly restrictive with the EU, for example, progressively limiting Cd contents from <60 to <20 mg/kg P 2 O 5 by 2032 (Ulrich 2019). Phosphate from igneous sources, such as carbonatites, is low in U and Cd, but relatively high in Th (van Kauwenbergh 1997), and the composition of the primary grains analysed in this study matches this finding (Fig. 10).
As a reference point, the concentration of Moroccan Phosphate Rock (MPR) standard BCR-032 (Serrini and Haemers 1980) is included in Fig. 10. Since this is a bulkrock standard, it is not a perfect comparison with apatite trace element contents but it provides an indication of the differing trace element content of phosphorites compared with igneous apatite sources. Notably, the Th concentration in secondary apatite from carbonatite is several orders of magnitude lower than that of the primary igneous grains from the same samples and is lower than the MPR phosphorite example (Fig. 10a). While this drop in Th content corresponds to an increase in U, the increase in U concentration is only approximately 1 order of magnitude, and it is only in the latest and most volumetrically minor stage that U contents exceed that of MPR (Fig. 10e). Increasing U contents in the carbonate-bearing fluorapatite also correspond to an increase in the concentration of other potentially deleterious elements such as Zn and Cd (Fig. 10). In particular, Cd in secondary apatite approaches and exceeds the mean MPR concentration although, as for U, this is only in the latest and most volumetrically minor stage (Fig. 10h).
Apatite is commonly viewed as a potential economic source of the REE, typically as a by-product or co-product of other commodities (Emsbo et al. 2015). In particular, hydrothermally formed carbonate-bearing fluorapatite can contain elevated contents of the HREE (Broom-Fendley et al. 2016). Morteani and Preinfalk (1996) suggested secondary carbonate-bearing apatite could be a major REE carrier in carbonatites. However, the samples analysed in this study indicate that such apatite has REE concentrations several orders of magnitude below the median content of phosphorites (Fig. 10c) and therefore appears not to represent an attractive source of the REE at this time.

Conclusions and implications
1 The textural, paragenetic and compositional similarity of carbonate-bearing fluorapatite from carbonatites indicates a common origin for the five, globally distributed examples investigated in this study. By textural comparison, an additional six carbonate-bearing fluorapatite deposits are also likely to have formed via the same mechanism. 2 The presence of a negative Ce anomaly and elevated contents of U in carbonate-bearing fluorapatite indicates formation in an oxidised environment which, based on comparison with other weathered rocks, is highly likely to indicate formation within the surface environment. 3 Considering the above two points, we infer that carbonatebearing fluorapatite deposits which are not evidently caused by modern weathering (such as Glenover and Kovdor) formed instead during palaeo-weathering. 4 Progressive formation of later generations of carbonatebearing fluorapatite increases the concentration of potentially deleterious elements for fertiliser production, such as Cd and U. However, only in the latest generation, comprising the smallest mass fraction of apatite, do these approach and exceed the corresponding mean contents observed in Moroccan phosphorites. Thus, carbonatebearing fluorapatite could be an attractive feedstock for fertiliser production, especially in light of increasing regulatory limits on the deleterious element content of phosphate rock. 5 The concentration of potentially economic REE is low overall and decreases in successive generations of carbonate-bearing fluorapatite, suggesting such rocks have limited potential as a source of REE.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.