the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 27 Jun 2024
Mg/Ca and δ18O in multiple species of planktonic foraminifera from 15 Ma to Recent
Abstract. The ratio of the trace element Mg over Ca (Mg/Ca) and the oxygen isotopic composition (δ18O) of foraminiferal calcite are widely employed for reconstructing past ocean temperatures, although geochemical signals are also influenced by several other factors that vary temporally and spatially. Here, we analyze a global dataset of Mg/Ca and δ18O data of 59 middle Miocene to Holocene species of planktonic foraminifera from a wide range of depth habitats, many of which have never been analyzed before for Mg/Ca. We investigate the extent to which Mg/Ca and δ18O covary through time and space, and identify several sources of mismatch between the two proxies. Once the data are adjusted for long term non-thermal factors, Mg/Ca and δ18O are overall positively correlated in a way consistent with temperature being the dominant controller of both through space and time and across many different species, including deep-dwellers. However, we identify several species with systematic offsets in Mg/Ca values, to which multispecies calibrations should be applied with caution. We can track the appearance of such offsets through ancestor-descendent species over the last 15 million years and propose that the emergence of these offsets may be the geochemical expression of evolutionary innovations. We find virtually all of the Mg/Ca and δ18O-derived temperatures from the commonly used genera Globigerinoides and Trilobatus are within uncertainty of each other, highlighting the utility of these species for paleoceanographic reconstructions. Our results highlight the potential of leveraging information from species lineages to improve sea surface temperature reconstruction from planktonic foraminifera over the Cenozoic.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(2306 KB)
-
Supplement
(115 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(2306 KB) - Metadata XML
-
Supplement
(115 KB) - BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-1608 by LJ de Nooijer', Lennart de Nooijer, 04 Sep 2024
Review Boscolo-Galazzo et al., EGUSphere
Dear editor,
I have carefully read the manuscript by Boscolo-Galazzo and co-workers you asked me to review. Although I am very enthusiastic about their approach and think that the manuscript is clearly written, there are two main issues that need to be addressed before this manuscript can be accepted for publication in BG. The second of these major comments, may lead to substantial changes in the manuscript, including the conclusions of this paper. In short, the authors need to substantiate their current claim that phylogeny and calcite chemistry are related by a more rigorous approach. In addition, I also listed some things that are less urgent, but do need attention.
Looking forward to see the authors’ response to my comments,
Lennart de Nooijer
Major comments:
- The statistical approaches should be explicitly described. This includes the exact way the regression analyses were done and the way the MC simulation was applied. Please provide the necessary details. Line 254: do the ‘calibration coefficients’ refer to the uncertainty in the estimated constants from the regression analysis? It would be interesting to know which of the sources of the uncertainty contribute most to the overall CI.
- The evolutionary perspective on differences in Mg/Ca-temperature relationships is very interesting (and slightly annoying for myself: I was just working on something similar…). However, the analysis of offsets in Mg/Ca and del18O versus evolutionary kinships is not making things clearer. Figures 4 and 5 actually suggest that offsets occur here and there, with no single ‘branch’ containing all the offset results. It may help to show which lineages fall towards the left-upper side of the average line and which towards the right-lower side (Figure 6). This would mean that three groups can be compared. But with the current figures, it remains very difficult to see if the species with a larger offset from the average (Figure 6) cluster together. I think this needs to be shown more rigorously if the conclusion of this exercise (i.e. that evolution is somehow responsible for the variability in Mg/Ca versus del18O) can hold. Also: the black color of some lines in figure 4 and 5 suggest that they are not offset from the average Mg/Ca-del18O, but with the red question mark beside them, we actually don’t know, right? This gives the false impression that more species than is maybe the case, fall around the average Mg/Ca-del18O relationship. Can figures 4 and 5 not be combined? If done, it will become clear that the spinose/ non-spinose divide doesn’t match the offset/ non-offset divide. This actually makes me most suspicious of the conclusion that there is a link between phylogeny and calcite chemistry (at least on this level).
Minor comments:
Line 30: remove ‘of’
Line 32: The authors argue in the discussion that deviations from the average Mg/Ca-del18O correlation is mainly due to inter-species variability in Mg/Ca (rather than variability in del18O). This may be the case, but interspecies differences in oxygen isotopes may still contribute to the variability shown in figures 6 and 7. Therefore, it may be more accurate here to say ‘systematic offsets in the Mg/Ca-del18O relationship’ or something similar.
Line 175: I guess you combined the ‘Mg/Ca-temperature calibrations’
Lines 205-208: it may be easier to simply state that you corrected the Mg/Cacc for past seawater Mg/Ca and pH (see 2.4.2). I guess it doesn’t matter for this procedure what Mg/Ca-temperature calibration you use (since you use it twice in ‘opposite’ directions).
Line 239: I don’t understand this: it seems that the effect of pH, salinity and temperature were taken into account twice. In lines 218-234 it was also described that the long-term changes in salinity and pH were corrected for. Conversion to temperature was already mentioned in line 238.
Line 240: what exponential coefficient? For the G. ruber Mg/Ca-pH sensitivity?
Line 250: what does this re-fit mean precisely? I understood that the Mg/Cacc were already corrected for salinity and pH…
Figure 2: so, since the different corrections are plotted ‘on top’ of each other, the order in which they are stacked matters for the suggestion that only part of the data could be accounted for with these corrections. The order could also be such that the upper left part of the data cloud falls within the calibration lines. It is not important for the paper, but perhaps the authors could try to change this? The applied corrections are actually vectors pushing the original calibration in different directions…
Figure 2: I don’t understand the Mg/Casw =2.5. That is (almost) half that of today, but in the methods there lacks a justification for this. Instead, the methods describe how the reconstructed change in Mg/Ca is used to correct the Mg/Cacc. It makes me wonder in figure 2, then what this line actually signifies: the caption suggests that the change in Mg/Casw from 15 Ma till now was used to correct the Mg/Ca of the foraminifera, but this would not (necessarily) result in a constant offset (as the orange line here suggests).
Figure 2B: what does the ’m = -2.08’ mean?
Line 284-285: ‘out of’ should be ‘from’.
Line 311: I tried to look it up, but couldn’t find how many datapoints were generated in total. The 62% sounds a higher percentage than it looks to me, but there may simply be many datapoints obscured by the density of the center of the cloud.
Line 352: should be ‘addition’
Line 434 and on: what do the ratios between parentheses mean?
Line 442: reference to Fig 5 is incorrect here. Also, references to figure 5 in the text of the manuscript appear after references to figure 6.
Figure 4: what are ‘offset’ spinose species? If it refers to species that have a del18O-MgCa relationship that deviates too much from the average relationship (figure 6), then what exactly makes it offset? What is the threshold value for such a distinction and what is this based on?
Line 539: should be ‘as in Figure 6’.
Citation: https://doi.org/10.5194/egusphere-2024-1608-RC1 -
AC1: 'Reply on RC1', Flavia Boscolo Galazzo, 25 Oct 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-1608/egusphere-2024-1608-AC1-supplement.pdf
-
RC2: 'Comment on egusphere-2024-1608', Anonymous Referee #2, 10 Sep 2024
General comments
This paper makes a positive contribution to our understanding of d18O and Mg/Ca as temperature indicators. The authors make use of a suitable dataset to address a number of knowledge gaps from both proxies. However, their final conclusions rely on a number of underlying assumptions some of which are still subject to uncertainties that are not clearly disclosed. My main concern is explained under "specific comments" below but once this is covered or acknowledged in some more detail I recommend that the manuscript be accepted for publication.
Specific Comments
d18O-temperatures: I’m aware that the online tool of Gaskell and Hull (2023) allows for a limited number of d18Osw corrections and that Rohling et al. (2021) (R_21) is a perfectly valid option. However, not only is d18Osw beyond the Pleistocene (or even just the LGM) quite uncertain, but also a recent benthic d18O deconvolution from Rohling et al. (2022) (R_22 https://doi.org/10.1029/2022RG000775) suggests heavier d18Osw during most of the last 15 million years compared to R_21, with the independent model output from De Boer et al. (2010) sitting somewhere in between R_21 and R_22. Meanwhile, clumped isotope derived d18Osw suggests even heavier values for parts of the Miocene (e.g. https://doi.org/10.1029/2020PA003927).
For these reasons, I wonder if it’d be possible to at least account for some of this uncertainty in your d18O-temperatures (larger vertical error bars?) and whether this might bring some of your “discordant" samples back into the uncertainty envelope of Figure 3. Alternatively, you could replot a few more versions of Figure 3 with different d18Osw corrections as offered in Gaskell and Hull (2023). Either way, I think it would be important to at least acknowledge these uncertainties in d18Osw before giving the reader the impression that a high degree of confidence exists behind the non-thermal corrections mentioned in the text. Moreover, this has implications for e.g. your hypothesis in Section 4.1.
Technical corrections
Title: The title is too general you could consider making it a bit more specific and refer to some of your main results.(Line) 220: The sea level record of Spratt and Liesecki only reaches 800 kyrs and doesn’t go as far as 8 Ma as stated on line 221.Also, I understand that in order to use the tool of Gray and Evans (2019) in samples older than 800 kyrs you probably had to modify the original code to include longer datasets for the salinity and and pH corrections (I only have experience with the R version so I’m not sure if the matlab version is a bit more flexible in this sense). If that were the case you could be more explicit with this in the text.221: Given that your choice of sea level/ benthic d18O deconvolution is Rohling et al., 2021 for the conversion of d18O into temperature (Line 259) I suggest using this same sea level record instead of Miller et al., 2005 for this step as well (from 800 kyrs up to 15 Ma, see comment for Line 220) to be as consistent as possible with the corrections applied to both proxies.229: This line suggests a very similar (or equivalent?) Mg/Casw than that one derived in Rosenthal, Bova and Zhou (2022) where they also combined Zhou et al., 2021 with the fluid inclusion data from Brennan et al., 2013. If it is then it might be more clear for the reader to cite the Rosenthal and Bova (2022) record directly (https://www.ncei.noaa.gov/access/paleo-search/study/36413).
Figures 3B-8B: These plots colored by calcification depth are so interesting it would be nice if you could increase the color contrast a bit more so that the 200-0 m range is more visible.Line 626: It goes from ~4 mol/mol at 5 Ma to ~5.25 mol/mol at present (I’m looking at Figure 10 in Rosenthal, Bova and Zhou, 2022).. did you perhaps mean 15 Ma?Citation: https://doi.org/10.5194/egusphere-2024-1608-RC2 -
AC2: 'Reply on RC2', Flavia Boscolo Galazzo, 25 Oct 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-1608/egusphere-2024-1608-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Flavia Boscolo Galazzo, 25 Oct 2024
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-1608 by LJ de Nooijer', Lennart de Nooijer, 04 Sep 2024
Review Boscolo-Galazzo et al., EGUSphere
Dear editor,
I have carefully read the manuscript by Boscolo-Galazzo and co-workers you asked me to review. Although I am very enthusiastic about their approach and think that the manuscript is clearly written, there are two main issues that need to be addressed before this manuscript can be accepted for publication in BG. The second of these major comments, may lead to substantial changes in the manuscript, including the conclusions of this paper. In short, the authors need to substantiate their current claim that phylogeny and calcite chemistry are related by a more rigorous approach. In addition, I also listed some things that are less urgent, but do need attention.
Looking forward to see the authors’ response to my comments,
Lennart de Nooijer
Major comments:
- The statistical approaches should be explicitly described. This includes the exact way the regression analyses were done and the way the MC simulation was applied. Please provide the necessary details. Line 254: do the ‘calibration coefficients’ refer to the uncertainty in the estimated constants from the regression analysis? It would be interesting to know which of the sources of the uncertainty contribute most to the overall CI.
- The evolutionary perspective on differences in Mg/Ca-temperature relationships is very interesting (and slightly annoying for myself: I was just working on something similar…). However, the analysis of offsets in Mg/Ca and del18O versus evolutionary kinships is not making things clearer. Figures 4 and 5 actually suggest that offsets occur here and there, with no single ‘branch’ containing all the offset results. It may help to show which lineages fall towards the left-upper side of the average line and which towards the right-lower side (Figure 6). This would mean that three groups can be compared. But with the current figures, it remains very difficult to see if the species with a larger offset from the average (Figure 6) cluster together. I think this needs to be shown more rigorously if the conclusion of this exercise (i.e. that evolution is somehow responsible for the variability in Mg/Ca versus del18O) can hold. Also: the black color of some lines in figure 4 and 5 suggest that they are not offset from the average Mg/Ca-del18O, but with the red question mark beside them, we actually don’t know, right? This gives the false impression that more species than is maybe the case, fall around the average Mg/Ca-del18O relationship. Can figures 4 and 5 not be combined? If done, it will become clear that the spinose/ non-spinose divide doesn’t match the offset/ non-offset divide. This actually makes me most suspicious of the conclusion that there is a link between phylogeny and calcite chemistry (at least on this level).
Minor comments:
Line 30: remove ‘of’
Line 32: The authors argue in the discussion that deviations from the average Mg/Ca-del18O correlation is mainly due to inter-species variability in Mg/Ca (rather than variability in del18O). This may be the case, but interspecies differences in oxygen isotopes may still contribute to the variability shown in figures 6 and 7. Therefore, it may be more accurate here to say ‘systematic offsets in the Mg/Ca-del18O relationship’ or something similar.
Line 175: I guess you combined the ‘Mg/Ca-temperature calibrations’
Lines 205-208: it may be easier to simply state that you corrected the Mg/Cacc for past seawater Mg/Ca and pH (see 2.4.2). I guess it doesn’t matter for this procedure what Mg/Ca-temperature calibration you use (since you use it twice in ‘opposite’ directions).
Line 239: I don’t understand this: it seems that the effect of pH, salinity and temperature were taken into account twice. In lines 218-234 it was also described that the long-term changes in salinity and pH were corrected for. Conversion to temperature was already mentioned in line 238.
Line 240: what exponential coefficient? For the G. ruber Mg/Ca-pH sensitivity?
Line 250: what does this re-fit mean precisely? I understood that the Mg/Cacc were already corrected for salinity and pH…
Figure 2: so, since the different corrections are plotted ‘on top’ of each other, the order in which they are stacked matters for the suggestion that only part of the data could be accounted for with these corrections. The order could also be such that the upper left part of the data cloud falls within the calibration lines. It is not important for the paper, but perhaps the authors could try to change this? The applied corrections are actually vectors pushing the original calibration in different directions…
Figure 2: I don’t understand the Mg/Casw =2.5. That is (almost) half that of today, but in the methods there lacks a justification for this. Instead, the methods describe how the reconstructed change in Mg/Ca is used to correct the Mg/Cacc. It makes me wonder in figure 2, then what this line actually signifies: the caption suggests that the change in Mg/Casw from 15 Ma till now was used to correct the Mg/Ca of the foraminifera, but this would not (necessarily) result in a constant offset (as the orange line here suggests).
Figure 2B: what does the ’m = -2.08’ mean?
Line 284-285: ‘out of’ should be ‘from’.
Line 311: I tried to look it up, but couldn’t find how many datapoints were generated in total. The 62% sounds a higher percentage than it looks to me, but there may simply be many datapoints obscured by the density of the center of the cloud.
Line 352: should be ‘addition’
Line 434 and on: what do the ratios between parentheses mean?
Line 442: reference to Fig 5 is incorrect here. Also, references to figure 5 in the text of the manuscript appear after references to figure 6.
Figure 4: what are ‘offset’ spinose species? If it refers to species that have a del18O-MgCa relationship that deviates too much from the average relationship (figure 6), then what exactly makes it offset? What is the threshold value for such a distinction and what is this based on?
Line 539: should be ‘as in Figure 6’.
Citation: https://doi.org/10.5194/egusphere-2024-1608-RC1 -
AC1: 'Reply on RC1', Flavia Boscolo Galazzo, 25 Oct 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-1608/egusphere-2024-1608-AC1-supplement.pdf
-
RC2: 'Comment on egusphere-2024-1608', Anonymous Referee #2, 10 Sep 2024
General comments
This paper makes a positive contribution to our understanding of d18O and Mg/Ca as temperature indicators. The authors make use of a suitable dataset to address a number of knowledge gaps from both proxies. However, their final conclusions rely on a number of underlying assumptions some of which are still subject to uncertainties that are not clearly disclosed. My main concern is explained under "specific comments" below but once this is covered or acknowledged in some more detail I recommend that the manuscript be accepted for publication.
Specific Comments
d18O-temperatures: I’m aware that the online tool of Gaskell and Hull (2023) allows for a limited number of d18Osw corrections and that Rohling et al. (2021) (R_21) is a perfectly valid option. However, not only is d18Osw beyond the Pleistocene (or even just the LGM) quite uncertain, but also a recent benthic d18O deconvolution from Rohling et al. (2022) (R_22 https://doi.org/10.1029/2022RG000775) suggests heavier d18Osw during most of the last 15 million years compared to R_21, with the independent model output from De Boer et al. (2010) sitting somewhere in between R_21 and R_22. Meanwhile, clumped isotope derived d18Osw suggests even heavier values for parts of the Miocene (e.g. https://doi.org/10.1029/2020PA003927).
For these reasons, I wonder if it’d be possible to at least account for some of this uncertainty in your d18O-temperatures (larger vertical error bars?) and whether this might bring some of your “discordant" samples back into the uncertainty envelope of Figure 3. Alternatively, you could replot a few more versions of Figure 3 with different d18Osw corrections as offered in Gaskell and Hull (2023). Either way, I think it would be important to at least acknowledge these uncertainties in d18Osw before giving the reader the impression that a high degree of confidence exists behind the non-thermal corrections mentioned in the text. Moreover, this has implications for e.g. your hypothesis in Section 4.1.
Technical corrections
Title: The title is too general you could consider making it a bit more specific and refer to some of your main results.(Line) 220: The sea level record of Spratt and Liesecki only reaches 800 kyrs and doesn’t go as far as 8 Ma as stated on line 221.Also, I understand that in order to use the tool of Gray and Evans (2019) in samples older than 800 kyrs you probably had to modify the original code to include longer datasets for the salinity and and pH corrections (I only have experience with the R version so I’m not sure if the matlab version is a bit more flexible in this sense). If that were the case you could be more explicit with this in the text.221: Given that your choice of sea level/ benthic d18O deconvolution is Rohling et al., 2021 for the conversion of d18O into temperature (Line 259) I suggest using this same sea level record instead of Miller et al., 2005 for this step as well (from 800 kyrs up to 15 Ma, see comment for Line 220) to be as consistent as possible with the corrections applied to both proxies.229: This line suggests a very similar (or equivalent?) Mg/Casw than that one derived in Rosenthal, Bova and Zhou (2022) where they also combined Zhou et al., 2021 with the fluid inclusion data from Brennan et al., 2013. If it is then it might be more clear for the reader to cite the Rosenthal and Bova (2022) record directly (https://www.ncei.noaa.gov/access/paleo-search/study/36413).
Figures 3B-8B: These plots colored by calcification depth are so interesting it would be nice if you could increase the color contrast a bit more so that the 200-0 m range is more visible.Line 626: It goes from ~4 mol/mol at 5 Ma to ~5.25 mol/mol at present (I’m looking at Figure 10 in Rosenthal, Bova and Zhou, 2022).. did you perhaps mean 15 Ma?Citation: https://doi.org/10.5194/egusphere-2024-1608-RC2 -
AC2: 'Reply on RC2', Flavia Boscolo Galazzo, 25 Oct 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-1608/egusphere-2024-1608-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Flavia Boscolo Galazzo, 25 Oct 2024
Peer review completion
Journal article(s) based on this preprint
Data sets
Study Dataset published as Boscolo-Galazzo et al. (2021) F. Boscolo-Galazzo et al. http://science.sciencemag.org/content/suppl/2021/03/10/371.6534.1148.DC1
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 366 | 180 | 174 | 720 | 42 | 23 | 22 |
- HTML: 366
- PDF: 180
- XML: 174
- Total: 720
- Supplement: 42
- BibTeX: 23
- EndNote: 22
Viewed (geographical distribution)
| Country | # | Views | % |
|---|---|---|---|
| United States of America | 1 | 253 | 34 |
| Germany | 2 | 58 | 7 |
| United Kingdom | 3 | 56 | 7 |
| China | 4 | 44 | 5 |
| Netherlands | 5 | 32 | 4 |
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
- 253
Flavia Boscolo-Galazzo
David Evans
Elaine Mawbey
William Gray
Paul Pearson
Bridget Wade
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(2306 KB) - Metadata XML
-
Supplement
(115 KB) - BibTeX
- EndNote
- Final revised paper