Skip to main content

Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists
Proceedings of the International Ocean Discovery Program Volume 363
publications.iodp.org

https://doi.org/10.14379/iodp.proc.363.205.2021

Data report: composite depth scale and splice revision for IODP Site U1488 (Expedition 363 Western Pacific Warm Pool) using XRF core scanning data and core images1

Anna Joy Drury,2, 3 Denise K. Kulhanek,4, 5 Thomas Westerhold,2 Ana Christina Ravelo,6 Gregory Mountain,7 Ann Holbourn,5 Yair Rosenthal,7 Anaïs Schmitt,8 Jennifer B. Wurtzel,9, 10 Jian Xu,11 Roy Wilkens,12 and Ursula Röhl2

Keywords: International Ocean Discovery Program, IODP, JOIDES Resolution, Expedition 363, Western Pacific Warm Pool, Site U1488, composite depth scale revision, composite splice revision, composite core images

MS 363-205: Received 5 May 2020 · Accepted 28 June 2021 · Published 29 October 2021

Abstract

The Western Pacific Warm Pool (WPWP) is a major source of heat and moisture to the atmosphere. Small perturbations in WPWP sea-surface temperatures greatly influence local Hadley and Walker cells, thereby affecting global atmospheric circulation patterns. International Ocean Discovery Program (IODP) Expedition 363 sought to document the regional expression and driving mechanisms of WPWP climate variability during the Neogene on millennial, orbital, and geological timescales. Located in the heart of the WPWP, IODP Site U1488 (02°02.59ʹN, 141°45.29ʹE) was drilled in 2604 m water depth on the southern part of the Eauripik Rise in the Caroline Basin. At Site U1488, a continuous shipboard composite stratigraphic section from 0 to ~331 m core composite depth below seafloor (CCSF) was compiled using high-resolution shipboard physical property data from three holes. This section comprises upper Miocene to recent foraminifer-rich nannofossil ooze and foraminifer-nannofossil ooze, making Site U1488 ideally suited to reconstruct the paleoceanographic history of the central WPWP region. However, the high carbonate content (>90% below ~180 m CCSF) of Site U1488 sediments means that the physical property data sets commonly used for splice construction (gamma ray attenuation bulk density, magnetic susceptibility, and natural gamma radiation) were too low amplitude to provide robust constraints on splice tie points below 120 m CCSF. As a result, P-wave data, which are relatively untested as a correlation tool, became critical for correlating between holes. Here, we verify and extend the Site U1488 shipboard composite splice using high-resolution (2 cm) X-ray fluorescence Ba/Sr core scanning data combined with composite linescan images. Overall, using these data at Site U1488 resulted in revised core offsets that differ by up to 0.84 m relative to the shipboard core offsets and a composite depth scale down to 329.33 m revised CCSF. The revised splice will allow optimization of postexpedition research and ensure that high-resolution studies of Site U1488 are conducted on a continuous stratigraphic section.

Introduction

The Western Pacific Warm Pool (WPWP) is the largest reservoir of oceanic heat on Earth and as a result plays a pivotal role in oceanic and atmospheric heat and transport processes (Rosenthal et al., 2018a). Forming part of the larger Indo-Pacific Warm Pool, the WPWP is a major source of heat and moisture to the atmosphere. Small perturbations in regional sea-surface temperatures (SSTs) greatly affect the local Hadley and Walker cells and perturb global atmospheric circulation patterns, thereby significantly affecting regional and global hydrological cycles (Rosenthal et al., 2018a). Prior to Expedition 363, expanded and continuous sedimentary archives documenting WPWP Neogene history were scarce. A key objective of International Ocean Discovery Program (IODP) Expedition 363 (6 October–8 December 2016) was recovert of continuous sedimentary archives to document the regional expression and driving mechanisms of WPWP climate variability throughout the Neogene on millennial, orbital, and geological timescales.

Three holes at IODP Site U1488 (02°02.59ʹN, 141°45.29ʹE; Figure F1) were drilled in 2604 m water depth on the southern part of the Eauripik Rise in the Caroline Basin. Because of its location in the heart of the WPWP, Site U1488 is ideally suited to reconstruct the Neogene paleoceanographic history of the central WPWP region (Rosenthal et al., 2018c). Three holes were cored to a maximum depth of ~315 m core depth below seafloor, method A (CSF-A; hereafter referred to as meters below seafloor [mbsf]), comprising upper Miocene to recent foraminifer-rich nannofossil ooze and foraminifer-nannofossil ooze. Shipboard studies show very good to good preservation of benthic and planktonic foraminifers throughout most of the sequence, and because the site remained above the calcite compensation depth throughout its history, dissolution effects on planktonic foraminifers are minimal (Rosenthal et al., 2018c). Based on shipboard nannofossil and planktonic foraminiferal biostratigraphic and paleomagnetic datums, sedimentation rates were ~2 cm/ky in the Pleistocene and ~4 cm/ky throughout the late Miocene and Pliocene. During the expedition, high-resolution shipboard physical property data were used to construct a continuous shipboard composite depth scale and stratigraphic splice from 0 to ~331 m core composite depth below seafloor (CCSF), which spanned from recent to ~9.4 Ma based on shipboard biostratigraphy (Rosenthal et al., 2018c).

Figure F1. Site map.

At Site U1488, P-wave velocity data, which are relatively untested as a stratigraphic correlation tool, were critical for determining the hole-specific cumulative core offsets (Figure F2) (Rosenthal et al., 2018b). This was necessary because the high carbonate content at Site U1488 (>90% below ~180 m CCSF) meant that the physical property data sets commonly used for splice construction (gamma ray attenuation [GRA] bulk density, magnetic susceptibility [MS], and natural gamma radiation [NGR]) had amplitude variability that was too low to provide a robust means for stratigraphic correlation between holes (Figure F2) (Rosenthal et al., 2018b; Rosenthal et al., 2018c). The shipboard composite depth scale and splice were predominantly generated using MS in the upper ~65 m CCSF (Figure F2A). Below 65 m CCSF, the amplitude of the MS data greatly decreased, and between 65 and ~120 m CCSF the composite depth scale and splice were constructed using a combination of MS, P-wave velocity, and GRA bulk density data. Below 120 m CCSF, the amplitude of the variability in the MS data was too low to use for stratigraphic correlation, so only P-wave velocity and GRA bulk density data were used (Figure F2B) (Rosenthal et al., 2018c).

Figure F2. MS, P-wave, and XRF Ba/Sr data.

The stratigraphic splice for Site U1488 must be robust to facilitate postexpedition research, which includes the generation of high-resolution paleoclimatological records to contribute to the overarching objectives of Expedition 363. However, research on other high-carbonate content settings (e.g., Ocean Drilling Program [ODP] Site 982 in the North Atlantic) has shown that physical property data can be highly unreliable for generating robust, continuous composite depth scales and stratigraphic splices (Drury et al., 2018). As such, it is crucial to verify shipboard composite depth scales and splices from high-carbonate sites such as Site U1488.

Here, we present linescan core images and high-resolution (2 cm) Ba/Sr X-ray fluorescence (XRF) core scanning data for Site U1488. We combine these two data sets with the available shipboard physical property data to verify and revise the shipboard composite depth scale and splice for Site U1488.

Methods and materials

XRF core scanning data

XRF core scanner Ba and Sr data were collected every 2 cm downcore over a 1.2 cm2 area with a downcore slit size of 10 mm and a cross-core slit size of 12 mm (50 kV; 0.5 mA; 7 s count time) from the split core surface of Site U1488 archive halves with the XRF Core Scanner III (AVAATECH Serial No. 12) at MARUM (the Center for Marine Environmental Sciences), University of Bremen (Germany). Where needed, the split core surface of the archive half was lightly scraped with a glass slide to achieve a flat surface suitable for XRF core scanning. The split core surface was then covered with a 4 μm thick SPEXCerti Prep Ultralene1 foil to avoid contamination of the XRF measurement unit and desiccation of the sediment. The data reported here utilize the SGX Sensortech Silicon Drift Detector (Model SiriusSD D65133Be-INF with 133eV X-ray resolution), the Topaz-X High-Resolution Digital MCA, and an Oxford Instruments 100W Neptune X-Ray tube with rhodium (Rh) target material. The raw X-ray data spectra were processed using the Iterative Least Square software (WIN AXIL) package from Canberra Eurisys.

Linescan core and spliced images

To assist with composite depth scale and splice verification, core images for Site U1488 were compiled using Code for Ocean Drilling Data (CODD; https://www.codd-home.net) (Wilkens et al., 2017) from the digital linescan images taken during the expedition. For each hole, cropped linescan section images were downloaded from the IODP laboratory information management system (LIMS) database (https://web.iodp.tamu.edu/LORE) while on board the R/V JOIDES Resolution. Images from Sections 363-U1488C-3H-2 through 7H-CC and 14H-5 through 17H-CC were slightly lightened because the original linescan images from these sections were obviously artificially darker than the actual sediments, as well as the preceding and succeeding section images. The corrected section-half images were then compiled into entire core images and scaled to the CCSF depth scale using the Includes_Core_Image_Assembly functions within CODD. A complete image of the entire stratigraphic splice was also compiled. All images compiled here are available on the PANGAEA database (https://www.pangaea.de).

Composite depth scale and stratigraphic splice revision

The core images and Ba/Sr data were combined with the shipboard physical property data to verify and, where needed, revise the shipboard composite depth scale and splice for Site U1488. The construction of a revised composite depth scale and splice generally followed the original shipboard methodology (Rosenthal et al., 2018b; Rosenthal et al., 2018c), which builds on the basic strategy first applied during ODP Leg 138 (Ruddiman et al., 1987; Hagelberg et al., 1992, 1995). During the expedition, Hole U1488C was used as the "backbone" of the upper section of the splice and Hole U1488B was used to cover core gaps. For the lower section, Hole U1488B was the backbone and Hole U1488A was used to cover core gaps (Rosenthal et al., 2018c). This strategy enables the use of off-splice core intervals to make a secondary near-continuous off-splice succession suitable for low-resolution sampling that avoids the more heavily sampled splice (see the OFFSPLICE folder in Supplementary material; Rosenthal et al., 2018d). Where appropriate, we maintained the shipboard interhole relationships and tie points when revising the CCSF depth scale and splice. Where the new data showed the shipboard offsets and splice required revision, we chose interhole ties and splice intervals that avoided the top and bottom ~50 cm and disturbed intervals of a core. We chose splice intervals from cores that best represented the stratigraphic characteristics of the section in question.

Here, we report depths on the core depth below seafloor, method A (CCSF-A) depth scale as meters below seafloor (mbsf) in line with the terminology used in the Expedition 363 Proceedings (Rosenthal et al., 2018b). The shipboard composite depth scale "core composite depth below seafloor" is referred to as CCSF. The revised composite depth scale is referred to as "revised CCSF." The CODD Includes_SplicingRevise_Menu was used to assist with the composite verification and revision process (Wilkens et al., 2017). As with the shipboard CCSF depth scale and splice, features used for correlating between holes and creating the revised CCSF depth scale and splice are only precisely correlated between holes at the exact location of the splice and/or offset tie point. Away from these ties, the correlation of individual stratigraphic features within a core may be less accurate between holes. Studies from other sites have provided an adjusted CCSF depth scale where the off-splice intervals are adjusted or "mapped" onto the splice to ensure as accurate a correlation between holes as possible (e.g., adjusted CCSF [Wilkens et al., 2013]; CCSF-M [van Peer et al., 2017)]. This could not be done for Site U1488 because the physical property data and core images were not distinct enough to use accurately and the XRF data were not collected for all holes.

Results

Revised composite depth scale

The XRF core scanning ratio of Ba/Sr provided a robust proxy for stratigraphic correlation between all three holes (Figure F3). The shipboard cumulative core offsets were initially rounded to the nearest centimeter. In total, the offsets for 26 cores were changed from the shipboard offsets (Table T1; Figure F3):

  • Cores 363-U1488A-14H, 21H, 23H, 26H, 31H, 34F, and 35F;
  • Cores 363-U1488B-2H, 8H, 10H–11H, 18H, 21H, 23H, 27H– 28H, 30H–31H, and 33H; and
  • Cores 363-U1488C-3H, 9H–10H, 12H–14H, and 16H.

Figure F3. Linescan and spliced core images.

Most changes are small, and the revised offsets generally differ 0.02–0.17 m relative to the shipboard offsets. A further three offsets were changed by 0.30–0.32 m (Cores 363-U1488A-26H and 34F– 35F; Figure F3K, F3N). Two of these changes were to the successive core gaps between Cores 363-U1488B-33H, 363-U1488A-34F, and 363-U1488A-35F, where the original growth factor–derived offset was increased by 30 cm to make it clearer that these cores are not part of a continuous succession (Table T1).

The largest changes to the shipboard offsets occurred in Cores 363-U1488C-12H (+0.74 m; Figure F3E) and 363-U1488B-27H (−0.75 m; Figure F3K) and 363-U1488B-31H (−0.84 m; Figure F3L, F3M). Changes to offsets were applied to every core in all holes below the depth of the core that was shifted, resulting in cumulative changes of 0.03–1.59 m in the revised CCSF scale relative to the shipboard CCSF scale (Table T1).

All preconditioned physical property data ("cleaned" GRA, MS, P-wave, and NGR) for Holes U1488A, U1488B, and U1488C were revised to the updated CCSF. These data are available as supplementary tables ST1–ST12 in the DATA_HOLE folder in the Supplementary material for this data report. During the expedition, preconditioning treatments were applied to the raw records to "clean up" and aid in the interpretation of noisy or spiky data (for details, see Physical properties in the Methods chapter [Rosenthal et al., 2018b]). The preconditioned physical property data ("cleaned" GRA, MS, P-wave, and NGR) provided along the revised stratigraphic splice are available as supplementary tables ST1–ST4 in the DATA_SPLICE folder in the Supplementary material for this data report.

Revised stratigraphic splice

The revised composite depth scale results in an updated stratigraphic splice down to 329.33 m revised CCSF (Table T2). Appended Cores 363-U1488A-34F and 35F extend the total sedimentary succession to 339.69 m revised CCSF overall (Table T2). The revised composite depth scale required adjustment of 25 splice tie points. Mostly, these changes constituted small revisions to optimize the location of the splice tie point between two holes within a sedimentary cycle. In addition, the exact interval and mbsf depth of some of the shipboard ties were adjusted by up to 1 cm to accommodate the rounding of the shipboard offsets to the nearest centimeter prior to revision of the composite depth scale (Tables T1, T2; both also available in the TABLES folder in the Supplementary material for this data report).

The largest revision occurred around ~134 m revised CCSF between Cores 363-U1488B-14H and 363-U1488C-14H, where a small section of Core 363-U1488A-14H was originally included in the shipboard splice. Following revision of the composite depth scale, Cores 363-U1488B-14H and 363-U1488C-14H were correlated at 134.13 m revised CCSF and Core 363-U1488A-14H was removed from the splice (Figure F3F). Generally, most splice ties are robust because there is excellent agreement between the XRF data, composite core photos, and physical property data among holes.

The splice tie point from Core 363-U1488A-32H to Core 363-U1488B-33H at 323.20 m revised CCSF remains ambiguous and should be treated cautiously because the XRF data did not provide a good correlation between the cores from different holes (Figure F3M). This interhole variability could be due to increased coring disturbance at the base of both holes. This tie point cannot be better resolved with the current data. However, because this tie point is located at the base of the revised splice, the consequences for high-resolution studies should be minimal.

Summary

We present a revised composite depth scale and stratigraphic splice for Site U1488, which was generated using XRF Ba/Sr core scanning data, linescan core images, and physical property data. The revised cumulative offsets differ by up to 0.84 m relative to the shipboard offsets. Splice tie points were adjusted where needed. Site U1488 now has a continuous and expanded revised composite depth scale down to 329.13 m revised CCSF. This revised splice will especially benefit postexpedition research that utilizes high-resolution studies on Site U1488 sediments and will contribute toward understanding orbital-scale climate dynamics in the central WPWP over the last ~9.4 My.

Acknowledgments

We thank the Captain, crew, and shipboard science party for their contributions to the recovery and preliminary analysis of the material recovered during IODP Expedition 363. We gratefully acknowledge the kind assistance of Holger Kuhlmann, Alex Wülbers, and Patrizia Geprägs for their help with core handling and Vera Lukies for her help with XRF core scanning. We thank Peter Blum for his thorough review, which helped us to improve the data report.

This research used samples and data provided by IODP, which is sponsored by the US National Science Foundation and participating countries. A.J. Drury and the XRF core scanning were funded by a MARUM Incentive Fund Award to A.J. Drury. A.J. Drury was additionally funded by DFG (Project Number 408101468) and by European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 796220. D. Kulhanek was funded by IODP JRSO NSF grant Number 1326927. Data are archived as supplementary tables in this volume and are available with the entire core and splice images in the PANGAEA database: https://www.pangaea.de.

References

Drury, A.J., Westerhold, T., Hodell, D., and Röhl, U., 2018. Reinforcing the North Atlantic backbone: revision and extension of the composite splice at ODP Site 982. Climate of the Past, 14(3):321–338. https://doi.org/10.5194/cp-14-321-2018

Hagelberg, T., Shackleton, N., Pisias, N., and the Shipboard Scientific Party, 1992. Development of composite depth sections for Sites 844 through 854. In Mayer, L., Pisias, N., Janecek, T., et al., Proceedings of the Ocean Drilling Program, Initial Reports. 138: College Station, TX (Ocean Drilling Program), 79–85. https://doi.org/10.2973/odp.proc.ir.138.105.1992

Hagelberg, T.K., Pisias, N.G., Shackleton, N.J., Mix, A.C., and Harris, S., 1995. Refinement of a high-resolution, continuous sedimentary section for studying equatorial Pacific Ocean paleoceanography, Leg 138. In Pisias, N.G., Mayer, L.A., Janecek, T.R., Palmer-Julson, A., and van Andel, T.H. (Eds.), Proceedings of the Ocean Drilling Program, Scientific Results. 138: College Station, TX (Ocean Drilling Program). https://doi.org/10.2973/odp.proc.sr.138.103.1995

Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists, 2017. Expedition 363 Preliminary Report: Western Pacific Warm Pool: International Ocean Discovery Program https://doi.org/10.14379/iodp.pr.363.2017

Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., Aiello, I.W., Babila, T.L., Bayon, G., Beaufort, L., Bova, S.C., Chun, J.-H., Haowen, D., Drury, A.J., Dunkley Jones, T., Eichler, P.P.B., Fernando, A.G.S., Gibson, K.A., Hatfield, R.G., Johnson, D.L., Kumagai, Y., Tiegang, L., Linsley, B.K., Meinicke, N., Mountain, G.S., Opdyke, B.N., Pearson, P.N., Poole, C.R., Ravelo, A.C., Sagawa, T., Schmitt, A., Wurtzel, J.B., Jian, X., Yamamoto, M., and Zhang, Y.G., 2018a. Expedition 363 methods. In Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists, Western Pacific Warm Pool. Proceedings of the International Ocean Discovery Program, 363: College Station, TX (International Ocean Discovery Program). https://doi.org/10.14379/iodp.proc.363.102.2018

Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., Aiello, I.W., Babila, T.L., Bayon, G., Beaufort, L., Bova, S.C., Chun, J.-H., Haowen, D., Drury, A.J., Dunkley Jones, T., Eichler, P.P.B., Fernando, A.G.S., Gibson, K.A., Hatfield, R.G., Johnson, D.L., Kumagai, Y., Tiegang, L., Linsley, B.K., Meinicke, N., Mountain, G.S., Opdyke, B.N., Pearson, P.N., Poole, C.R., Ravelo, A.C., Sagawa, T., Schmitt, A., Wurtzel, J.B., Jian, X., Yamamoto, M., and Zhang, Y.G., 2018b. Expedition 363 summary. In Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists, Western Pacific Warm Pool. Proceedings of the International Ocean Discovery Program, 363: College Station, TX (International Ocean Discovery Program). https://doi.org/10.14379/iodp.proc.363.101.2018

Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., Aiello, I.W., Babila, T.L., Bayon, G., Beaufort, L., Bova, S.C., Chun, J.-H., Haowen, D., Drury, A.J., Dunkley Jones, T., Eichler, P.P.B., Fernando, A.G.S., Gibson, K.A., Hatfield, R.G., Johnson, D.L., Kumagai, Y., Tiegang, L., Linsley, B.K., Meinicke, N., Mountain, G.S., Opdyke, B.N., Pearson, P.N., Poole, C.R., Ravelo, A.C., Sagawa, T., Schmitt, A., Wurtzel, J.B., Jian, X., Yamamoto, M., and Zhang, Y.G., 2018c. Site U1488. In Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists, Western Pacific Warm Pool. Proceedings of the International Ocean Discovery Program, 363: College Station, TX (International Ocean Discovery Program). https://doi.org/10.14379/iodp.proc.363.109.2018

Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists, 2018d. Supplementary material, https://doi.org/10.14379/iodp.proc.363supp.2018. Supplement to Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists, Western Pacific Warm Pool. Proceedings of the International Ocean Discovery Program, 363: College Station, TX (International Ocean Discovery Program). https://doi.org/10.14379/iodp.proc.363.2018

Ruddiman, W.F., Cameron, D., and Clement, B.M., 1987. Sediment disturbance and correlation of offset holes drilled with the hydraulic piston corer - Leg 94. In Ruddiman, W.F., Kidd, R. B., Thomas, E., et al., Initial Reports of the Deep Sea Drilling Project, 94: Washington, DC (US Government Printing Office), 615–634. https://doi.org/10.2973/dsdp.proc.94.111.1987

van Peer, T.E., Liebrand, D., Xuan, C., Lippert, P.C., Agnini, C., Blum, N., Blum, P., Bohaty, S.M., Bown, P.R., Greenop, R., Kordesch, W.E.C., Leonhardt, D., Friedrich, O., and Wilson, P.A., 2017. Data report: revised composite depth scale and splice for IODP Site U1406. In Norris, R.D., Wilson, P.A., Blum, P., and the Expedition 342 Scientists, Proceedings of the Integrated Ocean Drilling Program, 342: College Station, TX (Integrated Ocean Drilling Program). https://doi.org/10.2204/iodp.proc.342.202.2017

Wilkens, R.H., Dickens, G.R., Tian, J., Backman, J., and the Expedition 320/321 Scientists, 2013. Data report: revised composite depth scales for Sites U1336, U1337, and U1338. Proceedings of the Integrated Ocean Drilling Program, 320/321. https://doi.org/10.2204/iodp.proc.320321.209.2013

Wilkens, R.H., Westerhold, T., Drury, A.J., Lyle, M., Gorgas, T., and Tian, J., 2017. Revisiting the Ceara Rise, equatorial Atlantic Ocean: isotope stratigraphy of ODP Leg 154 from 0 to 5 Ma. Climate of the Past, 13(7):779–793. https://doi.org/10.5194/cp-13-779-2017


1 Drury, A.J., Kulhanek, D.K., Westerhold, T., Ravelo, A.C., Mountain, G., Holbourn, A., Rosenthal, Y., Schmitt, A., Wurtzel, J.B., Xu, J., Wilkens, R., and Röhl, U., 2021. Data report: composite depth scale and splice revision for IODP Site U1488 (Expedition 363 Western Pacific Warm Pool) using XRF core scanning data and core images. In Rosenthal, Y., Holbourn, A.E., Kulhanek, D.K., and the Expedition 363 Scientists, Western Pacific Warm Pool. Proceedings of the International Ocean Discovery Program, 363: College Station, TX (International Ocean Discovery Program). https://doi.org/10.14379/iodp.proc.363.205.2021

2 MARUM—Center for Marine Environmental Sciences, University of Bremen, Leobener Strasse, DE-28359 Bremen, Germany.

3 Department of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United Kingdom. Correspondence author: a.j.drury@ucl.ac.uk

4 International Ocean Discovery Program, Texas A&M University, 1000 Discovery Drive, College Station TX 77845, USA.

5 Institute of Geosciences, Christian-Albrechts-Universität zu Kiel, Ludewig-Meysnstresse 14, 24118 Kiel, Germany.

6 Ocean Sciences Department, University of California, Santa Cruz, 1156 High Street, Santa Cruz CA 95064, USA.

7 Department of Earth and Planetary Sciences, Rutgers University, 610 Taylor Road, Piscataway NJ 08854, USA.

8 Faculty of Science and Technology, Université de Nantes, 2 rue de la Houssiniére, 44322 Nantes, France.

9 Research School of Earth Sciences, The Australian National University, Building 142 Mills Road, Acton ACT 2601, Australia.

10 Now at New South Wales Department of Primary Industries, Orange Agricultural Institute, 1447 Forest Road, Orange NSW 2800, Australia.

11 Department of Geology, Northwest University, 229 North Taibai Road, Xi'an Shaanxi 710069, China.

12 School of Ocean and Earth Science and Technology (SOEST), University of Hawaii at Manoa, 2525 Correa Road, HIG 153, Honolulu HI 96822, USA.

This work is distributed under the Creative Commons Attribution 4.0 International (CC BY 4.0) license.