_{1}

^{*}

The article contains the results of statistical processing of a large summary of
*δ*
^{18}О-
*δ*
^{13}С isotope values in the primary carbonatites of the world. From literary sources, 1593 paired values
*δ*
^{18}О-
*δ*
^{13}С from 173 carbonatite occurrences of the world were collected. This report exceeds all previously published reports on С-О isotopes in carbonatites by quantity of the used values and carbonatite occurrences. Statistical data analysis is performed on diagrams in the coordinates
* δ*
^{18}О (‰, V-SMOW) -
*δ*
^{13}С (‰, V-PDV). For each carbonatite occurrence, not only the arithmetic mean values are calculated, but also the regression line. Distinct linear trend of
*δ*
^{18}О-
*δ*
^{13}С values is found in half of the carbonatite occurrences. The starting, middle, and ending points of the trend line are determined. The slope of the trend line (angular coefficient) varies over a wide range. The trend is dominated by an average angular coefficient of 0.30 (positive correlation
*δ*
^{18}О-
*δ*
^{13}С). In the literature, it is associated with the Rayleigh high-temperature fractionation of carbonatite melts or with their sedimentary contamination. Half of the carbonatite occurrences do not show a linear trend of
*δ*
^{18}О-
*δ*
^{13}С values, probably due to the combined action of multidirectional trends. The initial ratio
^{87}Sr/
^{86}Sr in the used carbonatite occurrences varies from 0.701 to 0.708. Statistics show no correlation of 87Sr/86Sr with the
*δ*
^{18}О-
*δ*
^{13}С system.

Oxygen and carbon isotope composition of carbonatites were summarized in a number of previous works. The largest number of δ^{18}О-δ^{13}С values (about 440) was collected and used to construct histograms in [^{18}О-δ^{13}С diagram for the Greenland, Europe, and North and South America regions without a division for individual carbonatite occurrences. The fields of primary igneous carbonatites on the δ^{18}О-δ^{13}С diagram are outlined in the works [

This paper uses 1593 pairs of conjugate values of δ^{18}О-δ^{13}С out of 173 carbonatite occurrences of the world. In addition, a linear regression analysis of the values is performed for most occurrences. This paper exceeds all previously published reports on С-О isotopes in carbonatites by quantity of the used values and carbonatite occurrences. Data on the ratio ^{87}Sr/^{86}Sr in 92 carbonatite occurrences are taken additionally from the sources used. The limited size of the article does not allow to provide a complete database and a list of used references (about 100 titles).

Carbonatite occurrences are represented by bodies of various shapes and sizes (complexes, massifs, dikes, facies zones). Isotope analysis is applied to carbonatite rocks (sovite, alvikite, beforsite, etc.), monofractions of calcite, dolomite, ankerite, siderite. Authors of literature classify the analyzed material as primary carbonatites (PC). This is mainly done on the basis of petrographic studies, in which secondary endo- and exogenous minerals are not detected. Single anomalous values are excluded from the primary category by the author of the article. Numerous δ^{18}О-δ^{13}С values refer to secondary carbonatites in the used literature: carbonate tuffs and lavas, hydrothermal veins, hydrothermally altered and recrystallized carbonatites, weathered and oxidized carbonatites, secondary calcite. Data on secondary carbonatites is not used in this article.

All isotope diagrams have a horizontal x-axis δ^{18}О (‰, V-SMOW) and a vertical y-axis δ^{13}С (‰, V-PDV). The equal scale of both axes, a multiple of 1‰, allows to visually comparing the shape of the point sets (point fields) and the slope of the trend lines. The names of carbonatite occurrences and their identification number (ID) are coordinated with the database [

The diagram in ^{18}О-δ^{13}С from 173 carbonatite occurrences, including various carbonatite facies in one occurrence. The number of points in the individual occurrences varies from 2 to 54. The points fill a very wide field. Analysis of the field is complex and incorrect due to the variable number and large scatter of points that characterize individual carbonatite occurrences. The contours PC-98% and PC-90% presented in the diagram are proved in

Subsequent statistical analysis of isotope data uses summary indicators characterizing carbonatite occurrences. Trend (linear regression) analysis was performed in 140 occurrences that satisfy two conditions: 1) there are three or more points; 2) the arithmetic difference between the maximum and minimum δ^{18}О values is more than 0.5‰. Trend lines under opposite conditions (less than three points and difference δ^{18}О less than 0.5‰) cannot be reliable. The linear regression equation y = kx + b and the trend line are calculated in Microsoft Excel 97-2003. The complete database (it is too large to be here) contains a point diagram for each occurrence, a calculated angular coefficient k, a constant b, a coefficient of determination (approximation) R^{2}.

Examples of point diagrams in order of increasing coefficient R^{2} are shown in ^{2} from 0.00 to 0.29. The number of such occurrences in the database is 70. The linear form of point fields begins to confidently be fixed from R^{2} ≥ 0.30 (occurrence 413. Chetlassky in ^{2} is from 0.30 to 0.99. The trend line is depicted as a vector directed upwards δ^{18}О. Such a direction is taken in the literature on the geochemistry of carbonatites.

Appendix provides a summary ^{18}О-δ^{13}С in occurrences without a linear trend (n) and with a linear trend (n*); 2) middle point arithmetic average δ^{18}О-δ^{13}С from among the values; 3) the starting point of the trend line or nonlinear field of δ^{18}О-δ^{13}С points; 4) minimum initial ratio ^{87}Sr/^{86}Sr. The position of carbonatite occurrences in the diagrams (

Trend line and middle point for 70 occurrences in which a linear trend is revealed are shown in the diagram (^{18}Оmax - δ^{18}Оmin) ranges from 0.5‰ to 11‰, in 90% of cases it does not exceed 7.5‰, on average it is 3.5‰.

The slope of the trend line varies widely, as seen in

The averaged shape of the field of points for three groups of occurrences with a linear trend and for one group without a trend is modeled in

arithmetic average differences δ^{18}Оmax - δ^{18}Оmin and δ^{13}Сmax - δ^{13}Сmin are calculated in each group. Rectangles with sides equal to these averages are shown in Figures 6(a)-(d). The modeled fields of points are inscribed in rectangles along the middle trend line. All fields in accordance with their averages are placed in ^{13}С value. This may be due to the cumulative effect of trends 6a, 6b and 6c.

Each carbonatite occurrence with a linear trend can be characterized by the starting point of the trend. The δ^{18}О of the starting point is equal to the minimum value in the statistical sample. The δ^{13}С value is calculated from the empirical regression equation. The values of δ^{18}О-δ^{13}С starting points of the trends are given in the Table. Occurrences without a linear trend also imply the presence of a starting point. This follows from the previously made assumption that the nonlinear point field 6d in ^{18}О, but in different directions along δ^{13}С. Therefore, the starting point of field 6d must have δ^{18}О equal to the minimum of the statistical sample. The δ^{13}С value in some approximation can be taken equal to the average of the sample (^{18}О-δ^{13}С of the starting point in the occurrences without the identified linear trend is also given in the Table.

The diagram shows two groups of points (^{18}О-δ^{13}С (n = 24) and with a difference δ^{18}Оmax - δ^{18}Оmin < 0.5 (n = 9) that were excluded from the trend analysis. Visual analysis of the diagram allows to delineate the internal contour PC-90% in addition to the PC-98% contour justified in ^{18}О = 7.9‰ divides the contour PC-90% into two parts, each of which is 45% of the total number of starting points.

The contours of primary carbonatites PC-98%, PC-90% and PC-45% (left and right contours) are shown in the diagram (_{2} from melts. Another subordinate trend k = 0.96 is not discussed in the literature. The beginning of the vectors is at the point (5‰ δ^{18}О; −6.5‰ δ^{13}С). The full sector of the angular coefficient (from −0.73 to 1.51) covers almost all occurrences from this point. Perhaps this point is close to the primary mantle source of carbonatites.

The used literature on O-C isotopy also contains data on the isotope composition of strontium. The minimum initial value of ^{87}Sr/^{86}Sr in 92 carbonatite occurrences is given in the Table. The field of minimum values in the coordinates ^{87}Sr/^{86}Sr-δ^{18}О is presented in the diagram (^{87}Sr/^{86}Sr = 0.702; δ^{18}O = 5‰) and the sedimentary contaminant (^{87}Sr/^{86}Sr = 0.710; δ^{18}O = 20‰) at equal concentrations of strontium in the sources. The stable enrichment of carbonatites with strontium in comparison with sedimentary carbonates is known. Under this condition, a band of points above the mixing line may reflect crustal contamination of magmas. However, a wide scatter of points below the mixing line leaves room for other hypotheses, including contamination of the source in the mantle. The PC-98%, PC-90% and PC-45% contours, previously substantiated in the coordinates δ^{18}О-δ^{13}С, are delineated in the diagram. The PC-45% contour is divided by the value ^{87}Sr/^{86}Sr = 0.704 into two fields. The field ^{87}Sr/^{86}Sr < 0.704 and δ^{18}О < 7.75‰ can be considered as the primary mantle field in the O-C-Sr isotope system.

Data on the oxygen and carbon isotope composition of primary carbonatites for 173 carbonatite occurrences of the world were collected (1593 paired values of δ^{18}О-δ^{13}С). Primary carbonatites are rocks without petrographic signs of secondary hydrothermal and exogenous mineral changes. Primary carbonatites demonstrate a wide variation of the δ^{18}О-δ^{13}С values and linear trends, which indicates the isotopic heterogeneity of carbonatite substance.

Linear regression analysis of δ^{18}О-δ^{13}С values reveals linear trends in half of the carbonatite occurrences. The trend with an average angular coefficient of 0.30 (positive correlation δ^{18}О-δ^{13}С) sharply dominates. In the literature, this is explained by the Rayleigh high-temperature fractionation of carbonatite melts or by their sedimentary (crustal) contamination. The trend line span (arithmetic difference δ^{18}Omax - δ^{18}Omin) ranges from 0.5‰ to 11‰, on average it is 3.5‰. Increased trends (over 7.5‰) suggest the action not only of endogenous factors, but also the influence of secondary processes not recorded in petrographic observations.

The second trend with an average angular coefficient of −0.27 (negative correlation δ^{18}О-δ^{13}С) is rarer. This trend is usually associated with the CO_{2} degassing from melts. A rare third trend is not discussed in the literature. It has an average angular coefficient of 0.96 (positive correlation δ^{18}О-δ^{13}С). The linear trend of δ^{18}О-δ^{13}С values is not detected in half of carbonatite occurrences due to increased variation of δ^{13}С. This may be due to the combined action of different factors—contamination, high-temperature fractionation and degassing of melts.

The fields of primary carbonatites (PC) are delineated in the coordinates δ^{18}О-δ^{13}С (‰), including 98%, 90% and 45% of the numbers of occurrences. The PC-90% contour can be considered acceptable for primary carbonatites. In-depth petrographic and other argumentation of the primary nature of carbonates is required for occurrences outside this contour. The PC-45% (δ^{18}О < 7.75‰) contour with a high probability includes only primary carbonatites with a mantle source of a carbonate substance and with minimal effect of isotope fractionation or contamination of melts. A greater influence of these factors is expected for occurrences in the PC-45% (δ^{18}О > 7.75‰) contour.

Strontium in carbonatite occurrences has a wide variation of the initial ^{87}Sr/^{86}Sr ratio from 0.701 to 0.708. This variation and the absence of correlation between ^{87}Sr/^{86}Sr and the δ^{18}О-δ^{13}С allow both mantle and crustal contamination of carbonatite magmas.

The stated statistical data on the O, C and Sr isotope composition in primary carbonatites leave room for additional and alternative judgments.

The author declares no conflicts of interest regarding the publication of this paper.

Bolonin, A.V. (2019) Oxygen and Carbon Isotope Composition in Primary Carbonatites of the World: Data Summary and Linear Trends. Open Journal of Geology, 9, 424-439. https://doi.org/10.4236/ojg.2019.98028

Carbonatite occurrences: δ^{18}О-δ^{13}С values of middle and starting points of trends; minimum initial ratio ^{87}Sr/^{86}Sr (n* – occurrences with a linear trend, n – other occurrences)