Mg/Ca ratios of Cenozoic benthic foraminifera from the Indian and South Atlantic Ocean

We investigate the evolution of Cenozoic climate and ice volume as evidenced by the oxygen isotopic composition of seawater (delta18Osw) derived from benthic foraminiferal Mg/Ca ratios to constrain the temperature effect contained in foraminiferal delta18O values. We have constructed two benthic foraminiferal Mg/Ca records from intermediate water depth sites (Ocean Drilling Program sites 757 and 689 from the subtropical Indian Ocean and the Weddell Sea, respectively). Together with the previously published composite record of Lear et al. (2002, doi:10.1126/science.287.5451.269) and the Neogene record from the Southern Ocean of Billups and Schrag (2002, doi:10.1029/2000PA000567), we obtain three, almost complete representations of the delta18Osw for the past 52 Myr. We discuss the sensitivity of early Cenozoic Mg/Ca-derived paleotemperatures (and hence the delta18Osw) to assumptions about seawater Mg/Ca ratios. We find that during the middle Eocene (~ 49-40 Ma), modern seawater ratios yield Mg/Ca-derived temperatures that are in good agreement with the oxygen isotope paleothermometer assuming ice-free conditions. Intermediate waters cooled during the middle Eocene reaching minimum temperatures by 40 Ma. The corresponding delta18Osw reconstructions support ice growth on Antarctica beginning by at least 40 Ma. At the Eocene/Oligocene boundary, Mg/Ca ratios (and hence temperatures) from Weddell Sea site 689 display a well-defined maximum. We caution against a paleoclimatic significance of this result and put forth that the partitioning coefficient of Mg in benthic foraminifera may be sensitive to factors other than temperature. Throughout the remainder of the Cenozoic, the temporal variability among delta18Osw records is similar and similar to longer-term trends in the benthic foraminiferal delta18O record. An exception occurs during the Pliocene when delta18Osw minima in two of the three records suggest reductions in global ice volume that are not apparent in foraminiferal delta18O records, which provides a new perspective to the ongoing debate about the stability of the Antarctic ice sheet. Maximum delta18Osw values recorded during the Pleistocene at Southern Ocean site 747 agree well with values derived from the geochemistry of pore waters (Schrag et al., 1996, doi:10.1126/science.272.5270.1930) further highlighting the value of the new Mg/Ca calibrations of Martin et al. (2002, doi:10.1016/S0012-821X(02)00472-7) and Lear et al. (2002, doi:10.1016/S0016-7037(02)00941-9) applied in this study. We conclude that the application of foraminiferal Mg/Ca ratios allows a refined view of Cenozoic ice volume history despite uncertainties related to the geochemical cycling of Mg and Ca on long time scales.

Sr/Ca, U/Ca, and stable isotope data and bimonthly records of Ogasawara coral core OGA-02-1

Instrumental climate observations provide robust records of global land and ocean temperatures during the twentieth century. Unlike for temperature, continuous salinity observations in the surface ocean are scarce prior to 1970, and the magnitude of salinity changes during the twentieth century is largely unknown. Surface ocean salinity is a major component in climate dynamics, as it influences ocean circulation and water mass formation. Here we present an annually resolved reconstruction of salinity variations in the surface waters of the western subtropical North Pacific Ocean since 1873, based on bimonthly records of d18O, Sr/Ca, and U/Ca in a coral from the Ogasawara Islands. The reconstruction indicates that an abrupt regime shift toward fresher surface ocean conditions occurred between 1905 and 1910. Observational atmospheric data suggest that the abrupt freshening was associated with a weakening of the winds that drive the Kuroshio Current system and the associated subtropical gyre circulation. We note that the abrupt early-twentieth-century freshening in the western subtropical North Pacific precedes abrupt climate change in the northern North Atlantic by a few years. The potential for abrupt regime shifts in surface ocean salinity should be considered in climate predictions for the coming decades.

Middle Miocene stable carbon and oxygen isotope record of foraminifera of ODP Hole 120-747A

Paleoceanographic variability at southern high latitude Ocean Drilling Program (ODP) Site 747 was investigated in this study through the interval which spans the Middle Miocene Climate Transition (MMCT). Between 15.0 and 12.2 million years ago (Ma), foraminiferal d18O records derived from both benthic (Cibicidoides spp.) and planktonic taxa (Globorotalia praescitula and Globigerina bulloides) reveal a history of changes in water column thermal and salinity structure and a strong imprint of seasonality. Prior to the MMCT, in the interval between 14.35 and 13.9 Ma, G. bulloides displays relatively high d18O values similar to those of G. praescitula, interpreted to indicate weakening of the thermocline and/or increased seasonality with cooler early-spring and/or late-fall temperatures. Following this interval, G. bulloidesd18O values diverge significantly from benthic and G. praescitula values, with G. bulloides values remaining relatively low for at least 600 kyr following the benthic foraminiferal d18O shift during the MMCT at ~13.9 Ma. This divergence in d18O records occurs in direct association with the Mi3 cooling and glaciation event and may suggest: (1) a strengthening of the vertical temperature gradient, with greater cooling of deep waters than surface waters, (2) changes in the depth habitat of G. bulloides, (3) changes in the dominant season of G. bulloides calcification, (4) modification of surface-water d18O values in association with enhanced sea-ice formation, (5) increased surface-water carbonate ion concentration, and/or (6) a significant decrease in surface-water salinity across the MMCT. The first of these possible scenarios is not likely, particularly in light of recent Mg/Ca evidence for significant surface-water cooling in the Southern Ocean associated with the MMCT. Of the remaining possibilities, we favor a change in surface salinity to explain the observed trends in d18O values and hypothesize that surface salinity may have decreased by up to 2 salinity units at ~13.9 Ma. In this scenario, the development of a lower-salinity Antarctic surface layer coincided with regional cooling of both surface and deep waters of the Southern Ocean during the Mi3 glaciation of East Antarctica, and contributed into the dominance of Neogloboquadrina spp. between 13.8 and 13.2 Ma. Additionally, the distinct patterns observed in planktonic foraminiferal d18O records spanning the MMCT correspond with changes in the vertical d13C gradient between planktonic and benthic foraminiferal records and major changes in planktonic foraminiferal assemblages at Site 747, providing further evidence of the environmental significance of this climatic transition.

Mg/Ca and stable oxygen measurement of the past 27 Ma of ODP Hole 120-747A

We explore the applicability of paired Mg/Ca and 18O/16O measurements on benthic foraminifera from Southern Ocean site 747 to paleoceanographic reconstructions on pre-Pleistocene timescales. We focus on the late Oligocene through Pleistocene (27-0 Ma) history of paleotemperatures and the evolution of the d18O values of seawater (d18Osw) at a temporal resolution of ~100-200 kyr. Absolute paleotemperature estimates depend on assumptions of how Mg/Ca ratios of seawater have changed over the past 27 Myr, but relative changes that occur on geologically brief timescales are robust. Results indicate that at the Oligocene to Miocene boundary (23.8 Ma), temperatures lag the increase in global ice-volume deduced from benthic foraminiferal d18O values, but the smaller-scale Miocene glaciations are accompanied by ocean cooling of -1°C. During the mid-Miocene phase of Antarctic ice sheet growth (~15-13 Ma), water temperatures cool by ~3°C. Unlike the benthic foraminiferal d18O values, which remain relatively constant thereafter, temperatures vary (by 3°C) and reach maxima at ~12 and ~8.5 Ma. The onset of significant Northern Hemisphere glaciation during the late Pliocene is synchronous with an ~4°C cooling at site 747. A comparison of our d18Osw curve to the Haq et al. (1987, doi:10.1126/science.235.4793.1156 ) sea level curve yields excellent agreement between sequence boundaries and times of increasing seawater 18O/16O ratios. At ~12-11 Ma in particular, when benthic foraminiferal d18O values do not support a further increase in ice volume, the d18Osw curve comes to a maximum that corresponds to a major mid-Miocene sea level regression. The agreement between the character of our Mg/Ca-based d18Osw curve and sequence stratigraphy demonstrates that benthic foramaniferal Mg/Ca ratios can be used to trace the d18Osw on pre-Pleistocene timescales despite a number of uncertainties related to poorly constrained temperature calibrations and paleoseawater Mg/Ca ratios. The Mg/Ca record also highlights that deep ocean temperatures can vary independently and unexpectedly from ice volume changes, which can lead to misinterpretations of the d18O record.

Geochemistry of serpentenites from DSDP Hole 84-566C

We studied a unique chrysotile-antigorite serpentinite, drilled on Deep Sea Drilling Project Leg 84 (Site 566) in the Guatemala forearc. Our in situ major and trace element data provide new constraints on possible reactions and associated trace element mobilisation during shallow serpentinite subduction. Chrysotile of the studied serpentinite, formed by the hydration of an upper mantle peridotite precursor, is partially replaced by antigorite (alone) which also occurs in 0.5 mm wide unoriented veins crosscutting the rock. Based on textural relationships and the P-T-X stability of the rock forming phases, the replacement of chrysotile by antigorite occurred at T < 300 °C, due to interaction between the chrysotile-serpentinite and an aqueous fluid. A comparison of the chemical compositions of reactant and product phases reveals that about 90% of the Cl, more than 80% of the B and about 50% of the Sr hosted originally by chrysotile was lost during fluid-assisted chrysotile-to-antigorite transformation and accompanying partial dehydration, and documents the much lower affinity of antigorite for trace element uptake than that of chrysotile. The fluid-assisted chrysotile-to-antigorite transformation and associated trace element loss documented here can occur in the shallow (< 30 km) region of subduction zones. This transformation decreases notably the Cl and B inventory of subducting serpentinites, which are regarded as one of the most important carriers of these elements into subduction zones. The evolution of serpentinites during initial subduction stages thus appears to be critical in the recycling of specific trace elements such as B or Cl from forearc to subarc depths.

Carbonate system of surface sediments from the South Pacific

B/Ca ratios in Cibicides mundulus and Cibicides wuellerstorfi have been shown to correlate with the degree of calcite saturation in seawater (D[CO32-]). In the South Pacific, a region of high importance in the global carbon cycle, these species are not continuously present in down-core records. Small numbers of epibenthic foraminifera in samples present an additional challenge, which can be overcome by using laser ablation-inductively coupled-mass spectrometry (LA-ICP-MS). We present a laser ablation based core-top calibration for Cibicides cf. wuellerstorfi, a C. wuellerstorfi morphotype that is abundant in the South Pacific and extend the existing global core top calibration for C. mundulus and C. wuellerstorfi to this region. B/Ca in C. cf. wuellerstorfi are linearly correlated with D[CO32-] and possibly display a higher sensitivity to calcite saturation changes than C. wuellerstorfi. Trace element profiles through C. wuellerstorfi and C. mundulus reveal an intra-shell B/Ca variation of ±36% around the mean shell value. Mg/Ca and B/Ca display opposite trends along the shell. Both phenomena likely result from ontogenetic effects. Intra-shell variability equals intra-sample variability, mean sample B/Ca values can thus be reliably calculated from averaged spot results of single specimen. In the global B/Ca-D[CO32-] range, we observe an inverse relationship between water mass age and D[CO32-].

Magnetic susceptibility of deep-sea cores from the Scotia Sea

We studied two deep-sea cores from the Scotia Sea to reconstruct past atmospheric circulation in the southern hemisphere and to resolve a long-standing debate on the interpretation of magnetic susceptibility (MS) records in Southern Ocean (SO) sediment. High-sedimentation sites MD07-3134 (0.2 - 1.2 m/kyr) and MD07-3133 (0.3 - 2 m/kyr) cover the last 92.5 kyr and 36 kyr, respectively. Both exhibit a one-to-one coupling of the MS and Ca2+ signal to the non-sea salt (nss) Ca2+ signal of the EDML ice core, clearly identifying atmospheric circulation as means of distribution. Comparison of additional proxies also excludes major influence by volcanic sources, sea-ice, icebergs, or oceanic current transport. The close resemblance of the dust proxies over the last glacial cycle, in turn, allows for the establishment of an age model of unprecedented resolution and precision for SO deep-sea sediment because atmospheric transport involves no major leads or lags. This is of particular importance because MS is routinely measured on deep-sea cores in the SO but the sediments usually lack biogenic carbonate and therefore had only limited stratigraphic control so far. Southern South America (SSA) is the likely source of eolian material because Site MD07-3133, located closer to the continent, has slightly higher MS values than Site MD07-3134, and also the MS record of Patagonian Site SALSA shows comparable variability. Patagonia was the dust source for both the Scotia Sea and East Antarctica. Dust fluxes were several times higher during glacial times, when atmospheric circulation was either stronger or shifted in latitude, sea level was lowered, shelf surfaces were exposed, and environmental conditions in SSA were dominated by glaciers and extended outwash plains. Hence, MS records of SO deep-sea sediment are reliable tracers of atmospheric circulation, allowing for chronologically-constrained reconstructions of the circum Antarctic paleoclimate history.

Geochemistry of last glacial sediments of Isla de los Estados, southeastern Tierra del Fuego

The island of Isla de los Estados is situated at 54.5°S, 64°W, east of Argentinian Tierra del Fuego, and is located in a sensitive geographic position in relation to the zonal circulation between Antarctica and South America. Its terrestrial records of the last deglaciation, recording atmospheric conditions but within an oceanic setting, can help to clarify changes of regional circulation patterns, both atmospheric and marine. Here, we present geochemical analyses from 16-10 ka cal BP of a peat core from Lago Galvarne Bog at the northern coast of the island, and a lake sediment core from Laguna Cascada 3 km further south. The data comprise TC, TN, loss on ignition analyses and continuous XRF scanning on both cores as well as age-depth modeling based on AMS-14C dating. Deglaciation and onset of peat formation in the coastal areas began before 16 ka cal BP followed by a rapid glacial retreat and the start of lacustrine sedimentation further inland. Data suggest initially windy conditions with permafrost succeeded by gradually warmer and wetter conditions until ca 14.5 ka cal BP. The warming trend slows down until ca 13.5 ka cal BP, followed by arid conditions culminating around 12.8 ka cal BP. Our data suggest fairly warm conditions and the establishment of denser peat and forest vegetation ca 10.6 ka cal BP, contemporaneous with the onset of the Antarctic thermal optimum. This indicates large-scale shifts in the placement of zonal flow and the Westerlies at the beginning of the Holocene.

Variability of the Brazil Current during the late Holocene analysed on two marine sediment cores

Our understanding of the centennial-scale variability of the Brazil Current (BC) during the late Holocene is elusive because of the lack of appropriate records. Here we used the Mg/Ca and oxygen isotopic composition of planktonic foraminifera from two marine sediment cores collected at 27° S and 33° S off southeastern South America to assess the late Holocene variability in the upper water column of the BC. Our results show in phase fluctuations of up to 3 °C in sea surface temperatures (SST), and 0.8 per mil in oxygen isotopic composition of surface sea water, a proxy for relative sea surface salinity (SSS). Time-series analyses of our records indicate a cyclicity with a period of ca. 730 yr. We suggest that the observed cyclicity reflects variability in the strength of the BC associated to changes in the Atlantic meridional overturning circulation (AMOC). Positive (negative) SST and SSS anomalies are related to a strong (weak) BC and a weak (strong) AMOC. Moreover, periods of peak strength in the BC occur synchronously to a weak North Brazil Current, negative SST anomalies in the high latitudes of the North Atlantic, and positive (negative) precipitation anomalies over southeastern South America (equatorial Africa), further corroborating our hypothesis. This study shows a tight coupling between the variability of the BC and the high latitudes of the North Atlantic mediated by the AMOC even under late Holocene boundary conditions.

Mg/Ca ratios of Neogloboquadrina pachyderma (s) and Globigerina bulloides in Irminger Sea sediment trap time series

iven the importance of high-latitude areas in the ocean-climate system, there is need for a paleothermometer that is reliable at low temperatures. Here we assess the applicability of the Mg/Ca-temperature proxy in colder waters (5-10?°C) by comparing for the first time the seasonal Mg/Ca and d18O cycles of N. pachyderma (s) and G. bulloides using a sediment trap time-series from the northern North Atlantic. While both species show indistinguishable seasonal d18O patterns that clearly track the near surface temperature cycle, their Mg/Ca are very different. G. bulloides Mg/Ca is high (2.0-3.1 mmol/mol), but varies in concert with the seasonal temperature cycle. The Mg/Ca of N. pachyderma (s), on the other hand, is low (1.1-1.5 mmol/mol) and shows only a very weak seasonal cycle. The d18O patterns indicate that both species calcify in the same depth zone. Consequently, depth habitat differences cannot explain the contrasting Mg/Ca patterns. The elevated Mg/Ca in pristine G. bulloides might be due to the presence of high Mg phases that are not preserved in fossil shells. The contrasting absence of a seasonal trend in the Mg/Ca of N. pachyderma (s) confirms other studies where calcification temperatures were less well constrained. The reason for this absence is not fully known, but may include species-specific vital effects. The very different seasonal patterns of both species' Mg/Ca underscore the importance of parameters other than temperature in controlling planktonic foraminiferal Mg/Ca. Our results therefore lend further caution in the interpretation of Mg/Ca-temperature reconstructions from high northern latitudes.

Mg/Ca and d18O measurements of planktonic foraminifera from Sargasso Sea

Paired Mg/Ca and d18O measurements on planktonic foraminiferal species (G. ruber white, G. ruber pink, G. sacculifer, G. conglobatus, G. aequilateralis, O. universa, N. dutertrei, P. obliquiloculata, G. inflata, G. truncatulinoides, G. hirsuta, and G. crassaformis) from a 6-year sediment trap time series in the Sargasso Sea were used to define the sensitivity of foraminiferal Mg/Ca to calcification temperature. Habitat depths and calcification temperatures were estimated from comparison of d18O of foraminifera with equilibrium calcite, based on historical temperature and salinity data. When considered together, Mg/Ca (mmol/mol) of all species, except two, show a significant (r = 0.93) relationship with temperature (T °C) of the form Mg/Ca = 0.38 (±0.02) exp 0.090 (±0.003)T, equivalent to a 9.0 ± 0.3% change in Mg/Ca for a 1°C change in temperature. Small differences exist in calibrations between species and between different size fractions of the same species. O. universa and G. aequilateralis have higher Mg/Ca than other species, and in general, data can be best described with the same temperature sensitivity for all species and pre-exponential constants in the sequence O. universa > G. aequilateralis = G. bulloides > G. ruber = G. sacculifer = other species. This approach gives an accuracy of ±1.2°C in the estimation of calcification temperature. The 9% sensitivity to temperature is similar to published studies from culture and core top calibrations, but differences exist from some literature values of pre-exponential constants. Different cleaning methodologies and artefacts of core top dissolution are probably implicated, and perhaps environmental factors yet understood. Planktonic foraminiferal Mg/Ca temperature estimates can be used for reconstructing surface temperatures and mixed and thermocline temperatures (using G. ruber pink, G. ruber white, G. sacculifer, N. dutertrei, P. obliquiloculata, etc.). The existence of a single Mg thermometry equation is valuable for extinct species, although use of species-specific equations will, where statistically significant, provide more accurate evaluation of Mg/Ca paleotemperature.

Geochemistry of ocean floor and forearc serpentinites

We provide new insights into the geochemistry of serpentinites from mid-ocean ridges (Mid-Atlantic Ridge and Hess Deep), passive margins (Iberia Abyssal Plain and Newfoundland) and fore-arcs (Mariana and Guatemala) based on bulk-rock and in situ mineral major and trace element compositional data collected on drill cores from the Deep Sea Drilling Project and Ocean Drilling Program. These data are important for constraining the serpentinite-hosted trace element inventory of subduction zones. Bulk serpentinites show up to several orders of magnitude enrichments in Cl, B, Sr, U, Sb, Pb, Rb, Cs and Li relative to elements of similar compatibility during mantle melting, which correspond to the highest primitive mantle-normalized B/Nb, B/Th, U/Th, Sb/Ce, Sr/Nd and Li/Y among subducted lithologies of the oceanic lithosphere (serpentinites, sediments and altered igneous oceanic crust). Among the elements showing relative enrichment, Cl and B are by far the most abundant with bulk concentrations mostly above 1000 µg/g and 30 µg/g, respectively. All other trace elements showing relative enrichments are generally present in low concentrations (µg/g level), except Sr in carbonate-bearing serpentinites (thousands of µg/g). In situ data indicate that concentrations of Cl, B, Sr, U, Sb, Rb and Cs are, and that of Li can be, increased by serpentinization. These elements are largely hosted in serpentine (lizardite and chrysotile, but not antigorite). Aragonite precipitation leads to significant enrichments in Sr, U and B, whereas calcite is important only as an Sr host. Commonly observed brucite is trace element-poor. The overall enrichment patterns are comparable among serpentinites from mid-ocean ridges, passive margins and fore-arcs, whereas the extents of enrichments are often specific to the geodynamic setting. Variability in relative trace element enrichments within a specific setting (and locality) can be several orders of magnitude. Mid-ocean ridge serpentinites often show pronounced bulk-rock U enrichment in addition to ubiquitous Cl, B and Sr enrichment. They also exhibit positive Eu anomalies on chondrite-normalized rare earth element plots. Passive margin serpentinites tend to have higher overall incompatible trace element contents than mid-ocean ridge and fore-arc serpentinites and show the highest B enrichment among all the studied serpentinites. Fore-arc serpentinites are characterized by low overall trace element contents and show the lowest Cl, but the highest Rb, Cs and Sr enrichments. Based on our data, subducted dehydrating serpentinites are likely to release fluids with high B/Nb, B/Th, U/Th, Sb/Ce and Sr/Nd, rendering them one of the potential sources of some of the characteristic trace element fingerprints of arc magmas (e.g. high B/Nb, high Sr/Nd, high Sb/Ce). However, although serpentinites are a substantial part of global subduction zone chemical cycling, owing to their low overall trace element contents (except for B and Cl) their geochemical imprint on arc magma sources (apart from addition of H2O, B and Cl) can be masked considerably by the trace element signal from subducted crustal components.

Geochemistry of peridotites from ODP Leg 125 diapiric serpentinites

Refractory spinel peridotites were drilled during Leg 125 from two diapiric serpentinite seamounts: Conical Seamount in the Mariana forearc (Sites 778-780) and Torishima Forearc Seamount (Sites 783-784) in the Izu-Ogasawara forearc. Harzburgite is the predominant rock type in the recovered samples, with subordinate dunite; no lherzolite was found. The harzburgite is diopside-free to sparsely diopside-bearing, with modal percentages of diopside that range from 0% to 2%. Spinels in the harzburgites are chrome-rich (Cr/[Cr + Al] = 0.38-0.83; Fe3+/[Fe3+ + Cr + Al] = 0.01-0.07). Olivine and orthopyroxene are magnesian (Mg# = 0.92). Discrete diopsides reveal extreme depletion of light rare earth elements. Primary hornblende is rare. The bulk major-element chemistry shows low average values of TiO2 (trace), Al2O3 (0.55%) and CaO (0.60%), but high Mg# (0.90). These rocks are more depleted than the abyssal peridotites from the mid-oceanic ridge. They are interpreted as residues of extensive partial melting (= 30%), of which the last episode was in the mantle wedge, probably associated with the generation of incipient island-arc magma, including boninite and/or arc-tholeiite. These depleted peridotites probably represent the residues of melting within mantle diapirs that developed within the mantle wedge.

Skeletal density of Porites

Paleotemperature estimates based on coral Sr/Ca have not been widely accepted because the reconstructed glacial-Holocene shift in tropical sea-surface temperature (~4-6°C) is larger than that indicated by foraminiferal Mg/Ca (~2-4°C). We show that corals over-estimate changes in sea-surface temperature (SST) because their records are attenuated during skeletogenesis within the living tissue layer. To quantify this process, we microprofiled skeletal mass accumulation within the tissue layer of Porites from Australasian coral reefs and laboratory culturing experiments. The results show that the sensitivity of the Sr/Ca and d18O thermometers in Porites will be suppressed, variable, and dependent on the relationship between skeletal growth rate and mass accumulation within the tissue layer. Our findings help explain why d18O-SST sensitivities for Porites range from -0.08 per mil/°C to -0.22 per mil/°C and are always less than the value of -0.23 per mil/°C established for biogenic aragonite. Based on this observation, we recalibrated the coral Sr/Ca thermometer to determine a revised sensitivity of -0.084 mmol/mol/°C. After rescaling, most of the published Sr/Ca-SST estimates for the Indo-Pacific region for the last ~14,000 years (-7°C to +2°C relative to modern) fall within the 95% confidence envelope of the foraminiferal Mg/Ca-SST records. We conclude that two types of calibration scales are required for coral paleothermometry; an attenuated Porites-specific thermometer sensitivity for studies of seasonal to interannual change in SST and, importantly, the rescaled -0.084 mmol/mol/°C Sr/Ca sensitivity for studies of 20th-century trends and millennial-scale changes in mean SST. The calibration-scaling concept will apply to the development of transfer functions for all geochemical tracers in corals.