Abstract
Neodymium (Nd) isotopes are considered a valuable tracer of modern and past ocean circulation. However, the promise of Nd isotope as a water mass tracer is hindered because there is not an entirely self-consistent model of the marine geochemical cycle of rare earth elements (REEs, of which Nd is one). That is, the prevailing mechanisms to describe the distributions of elemental and isotopic Nd are not completely reconciled. Here, we use published [Nd] and Nd isotope data to examine the prevailing model assumptions, and further compare these data to emergent alternative models that emphasize benthic processes in controlling the cycle of marine REEs and Nd isotopes. Our conclusion is that changing from a “top-down” driven model for REE cycling to one of a “bottom-up” benthic source model can provide consistent interpretations of these data for both elemental and isotopic Nd distributions. We discuss the implications such a benthic flux model carries for interpretation of Nd isotope data as a tracer for understanding modern and past changes in ocean circulation.
Introduction
Neodymium (Nd) isotopes (denoted εNd, which reflects the 143Nd/144Nd ratio normalized to a Chondritic Uniform Reservoir; Jacobsen and Wasserburg, ) are widely recognized as a valuable tracer for ocean circulation (e.g., Piepgras and Wasserburg, , ; Elderfield, ; von Blanckenburg, ; Frank, ; Goldstein and Hemming, ; Piotrowski et al., ; van de Flierdt and Frank, ). The potential for Nd isotopes to serve as a circulation tracer lies with the observation that εNd distributions appear to mirror the pattern of global deep-water thermohaline circulation. These observations support the assumption that εNd is conservative or “quasi-conservative” in the oceans and will trace water masses in a manner similar to temperature and salinity (Frank, ; Goldstein and Hemming, ). Furthermore, elemental Nd is assumed to have negligible bioactivity, although recent work suggests that there are important caveats to this assumption (Shiller et al., ). However, because of the concerted efforts invested into measuring neodymium in the oceans as well as within sedimentary archives, inconsistencies between the observed distributions and the assumptions inherent in the use of the Nd as a tracer are becoming increasingly apparent. The best known of these inconsistencies is the “Nd paradox” (Jeandel et al., , ; Tachikawa et al., ; Lacan and Jeandel, ). Described in many ways, this paradox fundamentally describes the contradictory observations that while εNd appears to behave conservatively, the water column profile of dissolved elemental Nd ([Nd]) appears to reflect the behavior of a reactive element, showing a nutrient-like distribution (see discussion in Goldstein and Hemming, ). In short, εNd distributions imply a residence time of ≤103 years, while [Nd] distributions imply a residence time of ≥104 years (Bertram and Elderfield, ; Jones et al., ). This discrepancy appears in paired geochemical mixing models (Goldstein and Hemming, ) and models that attempt to reproduce the global distribution of εNd (Tachikawa et al., ; Jones et al., ; Arsouze et al., ; Rempfer et al., ). In addition to the “Nd paradox,” inconsistencies have arisen regarding past records of authigenic εNd, which sometimes prove to be difficult to explain when assuming water mass mixing of fixed end-members (Osborne et al., ; Stewart et al., ), or show exceptional disparities between adjacent sites (Stumpf et al., ; Roberts and Piotrowski, ; Howe et al., ; Hu et al., ).
The problem with εNd as a circulation proxy is fundamentally that while εNd appears to trace water masses effectively, it is not particularly potent at discerning geochemical processes that might impact its distribution. For example, it is difficult to use εNd to determine how a water mass acquires its isotope signature. Although the pattern of global εNd mirrors that of thermohaline circulation the cause of that apparent coincidence is not straightforward. For example, that the deep North Pacific is a Nd isotope end-member for mixing is difficult to justify given that no deep water forms within this basin. If isotopic end members can develop unassociated with preformed properties of the water mass, then there is clearly a need to explain and constrain the effects of these processes on εNd, especially for interpretation of εNd as ancient records of circulation.
Relatively recently, several possible mechanisms for non-conservative behavior of εNd have been forwarded: boundary exchange, Submarine Groundwater Discharge (SGD) and a benthic Nd flux (Lacan and Jeandel, ; Johannesson and Burdige, ; Abbott et al., ). These three hypotheses are all similar in that bottom water can be altered along the sediment-water boundary. Although these hypotheses may have been articulated as distinctly different, the differences among these hypotheses may be simply semantic or reflect our incomplete understanding of the processes involved. Boundary exchange, as described by Lacan and Jeandel () is typically thought to reflect an isotope exchange of Nd with net zero concentration change that occurs between the sediments and bottom water predominantly at ocean margins (Jeandel, ). While this model is highly cited, the specific mechanism or processes are not well defined. In contrast, both the SGD and benthic flux models propose more distinct mechanisms but lack the robust global data sets needed for verification. All three of these hypotheses invoke a potential source of dissolved [Nd] to the oceans from sedimentary fluids; SGD waters are considered here to be originally fresh waters, while the benthic flux model considers seawater as the source fluid origin. These latter hypotheses will differ in the quantity and location of the Nd fluxes they predict and the εNd that they carry. Clearly, the distinctions among all three models is nuanced, and again, probably simply reflects the need for further investigation or even simple definition. However, all three models challenge the underlying assumptions around the conservative nature of εNd in the water column.
A more comprehensive understanding of the marine geochemical cycle of REEs in conjunction with εNd is needed to build more accurate geochemical models and interpret down core records with respect to ocean circulation. The problems arising from the notion of εNd as a paleoproxy perhaps reflect the superficial, and perhaps data limited synthesis between analyses of εNd and the rare earth elements (REEs) throughout the ocean basins. Fortunately, the comprehensive nature of the GEOTRACES water column sampling will go far to improve our understanding of water column processes that affect the REEs and εNd. Here, we discuss published εNd and Nd data from the Atlantic Ocean, with the intent to clarify our understanding of these tracers where they appear to be most robust. While we do not present any new data, we attempt to illustrate that a holistic view of published REE and εNd distributions in the oceans is consistent with predominant control from a benthic sedimentary source. Our conclusion is that processes occurring within the water column (e.g., reversible scavenging) are likely a secondary control on εNd or [REE] in the oceans, save where particle fluxes are exceptionally high, and that benthic fluxes exert primary control over the distribution of these elements and their isotopes.
Modern ocean REEs and Nd isotopes
It was originally held that the REEs are input to the oceans predominantly via rivers, with potential additions via dust where such inputs are relatively high (Byrne and Sholkovitz, ; Greaves et al., ; Tachikawa et al., ; Goswami et al., ; Dunlea et al., ; Stichel et al., 2015; also see Bayon et al., for the alternative view that dust is a sink term). More recently, the emphasis of REE input to the oceans has shifted to greater consideration of a benthic continental margin source for REEs, as we will discuss further later (Spivack and Wasserburg, ; Jeandel et al., ; Goldstein and Hemming, ; Tachikawa et al., ; Arsouze et al., ; Rickli et al., , , ; Carter et al., ; Grasse et al., ; Singh et al., ; Grenier et al., ; Wilson et al., ; Garcia-Solsona et al., ).
In the open ocean, REEs are thought to be transferred from shallow to deep ocean via reversible scavenging on particles sinking through the water column (Elderfield, ; Byrne and Kim, ; Sholkovitz et al., ; Byrne and Sholkovitz, ). However, the heavier atomic mass REEs (HREEs, such as Yb, Lu) tend to have vertical profiles that are more “Si-like,” which is thought to reflect the tendency for greater complexation across the REE series (Cantrell and Byrne, ; Elderfield, ; Byrne and Sholkovitz, ). Finally, removal of REEs from the ocean is thought to be predominantly via metal-oxides, although the role of organic matter, phosphates, or even carbonates, is likely to be important (Byrne and Kim, Byrne and Sholkovitz, ; Schijf et al., ). This model is considered the standard of marine REEs. Nd isotopes should also behave in a manner that is consistent with this model, although as pointed out previously this does not appear to be the case (Goldstein and Hemming, ; Jones et al., ; Arsouze et al., ). Successful efforts to reconcile the inconsistent behavior between elemental Nd and its isotopes, i.e., the Nd paradox, will likely begin with an examination of these first principles of marine REE cycling.
The shift from considering riverine point sources to broad marginal inputs improves model results (Tachikawa et al., ; Arsouze et al., ) and offers an explanation for some of the irregularities observed when comparing Nd and its isotopes such as the constancy of surface ocean water εNd (−7.8 upstream vs. −7.9/−8.0 downstream; Osborne et al., ) around the mouth of the Mississippi river that carries a far less radiogenic signal (at εNd < −11; Goldstein and Jacobsen, , ; Bayon et al., . Budgetary estimates from diffuse benthic REE fluxes, such as SGD or from sedimentary pore water (Sholkovitz et al., ; Greaves et al., ; Haley and Klinkhammer, ; Tachikawa et al., ; Johannesson and Burdige, ; Schacht et al., ; Johannesson et al., , ; Abbott et al., ; Fröllje et al., ), suggest that they may dominate the flux of REEs to the oceans. Importantly, such diffuse, or at least hard to quantify, sources are consistent with deep-water isopycnal mixing of εNd signals off margins that cannot be explained through surface water sources (Grasse et al., ; Grenier et al., ). However, while a diffuse benthic source of Nd appears to better describe the dominant flux of REEs to the ocean, outstanding questions remain. Arguably the most important of these is: what drives this flux? Any mechanism must be reconciled with the remarkably consistent REE pattern of seawater in light of diverse sediment compositions and sedimentary/diagenetic environments. As such, it seems unlikely that simple, unidirectional particulate dissolution can be the mechanism, because the dissolved REE pattern and εNd would to some extent reflect such an input, following arguments posed by Sholkovitz () and reiterated by Bau et al. ().
Another supposition in the current understanding of marine REE cycling is that reversible scavenging controls the vertical distribution of these elements in the water column. In his deeply insightful paper, Elderfield () suggested cause for concern over this conjecture: not for shallow REE adsorption, but in light of a mechanism for desoprtion at depth. For many trivalent cations the pH-edge/Langmuir-front for adsorption onto particle surfaces lies well under a pH of 7 (Morel and Hering, ; Ngwenya et al., ); thus, over a typical oceanic pH range of 7.6 and 8.2 desorption is unlikely. The thermodynamics of desorption become even less favorable as the increasing dissolved concentration at depth would tend to push the equilibrium toward further adsorption (reiterating arguments posed by Elderfield, ). Moreover, given free Nd is <16% of the total dissolved Nd (Millero, ; Schijf et al., ), and surface adsorption coefficients are typically not exceptionally high (in the range of 3 to 5 onto carboxyl and moncarboxylic groups; Smith and Martell, ; Byrne and Kim, ; Ngwenya et al., ), the potential for REE transfer from the shallow to deep ocean via particle surface scavenging should be limited (Elderfield and Greaves, ; Stichel et al., 2015).
In spite of such arguments, lab and field observations attest that the REEs are indeed highly susceptible to adsorption (Sholkovitz, ), and that adsorption onto particles in the ocean certainly does happen (e.g., Sholkovitz et al., ; Tachikawa et al., ). One positive test for REE adsorption in the oceans may be seen in conservative Arctic Ocean REE profiles where particles are near absent (Yang and Haley, ). We do note, however, that a counter-argument can be made in that there does not appear to be an appreciable difference in the [Nd] profiles relating to marginal sites vs. distal sites, where particle loading may be very different (German and Elderfield, ; Grasse et al., ; Goswami et al., ; Haley et al., ; Abbott et al., ; Stichel et al., 2015).
Considering all the arguments and data, a preferred model is that REEs are most likely transferred from the ocean surface to depth through incorporation onto complex organics, such as humic acids, polysaccharides, or other complex organic molecules that have binding coefficients higher than simple carboxyl groups (logK > 9; Byrne and Kim, ; Stanley and Byrne, ; Davranche et al., ; Pourret et al., ; Schijf et al., ). In this way, scavenging is actually achieved via an organic coating carrier phase that is remineralized, thus offering a mechanism for transfer of REEs from the particulate to dissolved phase without subsequent re-adsorption (i.e., not a desorptive process, but a remineralization; Sholkovitz et al., ). This model implies that [Nd] should behave much like δ13C in the water column, wherein their distributions reflect a preformed signal modified by top-down uni-directional input additions related to POC remineralization. Available GEOTRACES data allow us to make an initial comparison of water column [Nd] and δ13C (Mawji et al., ; Stichel et al., 2015; Figure 1). We consider the comparison inconclusive; these data appear to share many similarities (e.g., at Stations 12, 24, 10, especially in the upper 1,000 m) but not everywhere (e.g., at Stations 1, 3, and 7). Comparison of these [Nd] and δ13C data yields no correlation (r2 < 0.05; Figure S1), although the fit of all the data from the upper 50 to 500 m, ostensibly under highest remineralization, is better (r2 ~0.4; Figure S2), which would support scavenging of Nd on isotopically light-C POM. The form of the [Nd] profiles in Figure 1, especially in the upper 1,000 m, does appear to share many of the characteristics of POM remineralization in the oceans (e.g., Suess, ; Martin et al., ) and of AOU (discussed by Stichel et al., 2015). Unfortunately, global compilations, such as presented by Tachikawa et al. (), produce poor correlations between δ13C and [Nd] (r2 < 0.001; not shown), and the similarity between deep (>1,000 m) [Nd] and P* (van de Flierdt et al., ) becomes difficult to justify under the hypothesis of POM scavenging of Nd.
Figure 1
A similar uni-directional release of REEs at depth can also potentially come from deterioration of the lithogenic minerals themselves (Rousseau et al., ; Abadie et al., ). This mechanism would differ markedly in its effect on water column [REE] and εNd in comparison to a model based on an organic surface coating. Unfortunately, the evaluation of [Nd] and εNd in light of a lithogenic source of Nd at depth has been ambiguous in the eastern North Atlantic (Stichel et al., 2015; van de Flierdt et al., ).
Evaluation of modes of scavenging may be complicated through multiple substrates for scavenging, each of which can deliver REEs to depth (e.g., complex organic molecules as well as lithogenic and biogenic silicates; Akagi et al., ; Akagi, ). Adding additional ambiguity to the assessment of reverse scavenging is the issue of preformed REEs at depth that are likely conservative in nature (Haley et al., ; Zheng et al., ). Such preformed REEs must be an important component of deep water REEs (at least on basinal scales), such that REE profiles cannot be taken at face value: de-convolution of the conservative and non-conservative REEs is needed for accurate interpretation (Bertram and Elderfield, ; Haley et al., ; Zheng et al., ). Such efforts, either using a water mass approach (e.g., Bertram and Elderfield, ; Grasse et al., ; Stichel et al., , 2015 or a more statistical approach (Singh et al., ; Haley et al., ; Zheng et al., ), are critical in evaluating the importance of scavenging on water column REEs. It goes beyond the scope of this paper to explore solutions for resolving preformed vs. “reactive” REEs. Our goal here is simply to suggest that much of the difficulty in resolving the “Nd paradox” is caused by undue emphasis placed on the function of reverse scavenging, which likely plays only a minor role in controlling the distribution of REEs and εNd in the water column.
Our exemplar here is the outstanding work of Bertram and Elderfield (). These authors carefully evaluated the possible influence of scavenging and reverse scavenging, largely under an assumption of “top down” forcing, because of the similarities to Si and a lack of sedimentary data (Bertram and Elderfield, ). Their conclusion is that particle-water exchange is necessary to explain deep ocean εNd, and give two possible mechanisms: a “top down” possibility wherein scavenged Nd exchanges over time in the deep ocean and a “…second option [that] simply places the site of exchange close to or within the sediment-porewater system with exchange within the nepheloid zone or transport via diffusion into overlying seawater” (Bertram and Elderfield, ). [Interestingly, the authors also note that “the simple (water density) & Si-REE correlations seen in bottom and near-bottom waters may also imply involvement of a benthic source. Unfortunately, no satisfactory data are available on REEs in oxic porewaters.”].
Finally, regardless of the impact on the water column, we emphasize that all of the potential “top down” mechanisms for remobilizing surface REEs at depth would be greatly enhanced within seafloor sediments (Freslon et al., ); i.e., manifest as a pore water flux. In summary, while scavenging certainly happens on particles, we argue that the REE mass transfer is likely relatively small and subsequent desorption at depth in the water column is unlikely. On the other hand, any such desorption or remineralization mechanism is likely to be significant in the sediments, where the chemical environment can change significantly.
Benthic control
Regardless of the discussions presented above, mechanisms driven via a “top-down” flux of REEs return to the disparity between εNd and [Nd] as discussed by Goldstein and Hemming (). The question is fundamentally how can such processes occur that account for [Nd] distributions, yet does not obviate the conservative qualities of εNd? That a water mass such as North Atlantic Deep Water (NADW) can maintain its signature on ocean-basin scales is arguably less justifiable in a “top-down” explanation of marine REEs than it is in the “bottom-up” explanation that will be discussed below.
If we consider our prior argument that reversible scavenging does not significantly impact Nd's water column distribution, what then can explain the distribution of [Nd]/εNd at depth? A “bottom up” hypothesis argues that Nd and εNd are truly conservative below the permanent thermocline, except when exposed to a benthic flux. This benthic flux will certainly impact overlying bottom water, but its influence may also be carried laterally along isopycnals, for example off the continental slope (as indicated by Grasse et al., ; Singh et al., ; Stichel et al., 2015). The interaction with the sediments is then what makes εNd “quasi-conservative,” vs. being truly conservative, as has been suggested previously (Lacan and Jeandel, ; Abbott et al., ). The two important aspects of this suggestion are that (1) in the absence of a positive benthic flux, Nd and εNd will behave conservatively and (2) the ocean biogeochemistry of Nd and εNd is thus defined by benthic diagenetic processes. That is, this benthic control model suggests that the REEs and εNd are controlled by the nature of the benthic flux, which can be a positive, negative, or neutral term.
The idea of a benthic source for REEs, as described previously (Abbott et al., ), would certainly help resolve model and observational inconsistencies of REE in the Pacific Ocean (Jones et al., ; Arsouze et al., ), but is such a source necessary anywhere else-in the Atlantic, Indian, Arctic or Southern Oceans? The North Atlantic, for example, has been well described in “top-down” models (Jones et al., ), or even through simple conservative mixing (van de Flierdt et al., ). Although we note here that the Denmark Straits is the ‘birthplace’ of the very non-conservative boundary exchange hypothesis (Lacan and Jeandel, ,, ) and that other models have required benthic source fluxes in the North Atlantic (Tachikawa et al., ).
Why then invoke a benthic model? First, if correct, the benthic model would provide a mechanism to reconcile the apparent conflicts in our interpretation of Nd and εNd; and second, the benthic model makes different predictions about how changes in εNd and REEs records can be interpreted, when compared to a reversible scavenging model, as will be discussed later.
The concept of benthic control of [REE]/εNd on Pacific waters is still a nascent idea (Abbott et al., , ; Du et al., ). Here, we expand on this idea in a general sense. First, benthic control of [REE]/εNd will directly influence only the bottom water, such that a water mass detached from the bottom will maintain its εNd so long as it maintains its temperature and salinity (i.e., all its conservative properties). Second, benthic control represents a diffuse boundary flux, which implies that deep water [REE]/εNd distributions may be more disposed to being controlled by eddy diffusivity. Because eddy diffusion coefficients in the vertical are 10−5-10−4 m2/s (or greater; Ku and Luo, ), over hundreds to thousands of meters per year (horizontal coefficients are even greater), the impact of a benthic flux may potentially integrate over significant depths within the deep ocean. Furthermore, eddy diffusivity is also an order of magnitude higher on the seafloor relative to the water column because of the turbulence created by topography (Waterhouse et al., ); therefore, a benthic flux will be dissipated with great efficiency into the water column. This dissipation is the rationale for the observation that there is an apparent jump in REE concentrations between pore water and overlying water (Abbott et al., , ).
Implications of a benthic flux in the atlantic
There is little direct observational data regarding the benthic REE flux to the oceans (Sholkovitz et al., ; Haley et al., ; Abbott et al., ), and even less that tracks the εNd of this flux. However, comparison with radium offers some support for the hypothesis. Again, the holistic nature of GEOTRACES sampling has provided rare, directly comparable Ra-REE-εNd data (Charette et al., ; Stichel et al., 2015; USGT10/GA03 section from the Gulf of Cadiz to the Mauritanian coast). All isotopes of Ra positively correlate with [Nd] in the bottom water, although we plot only 228Ra, as this isotope is most directly associated with a benthic source influence (Figure 2A). The correlation of 228Ra with [Nd] in bottom water is strong (r2 = 0.7 including all the data, and by omitting Site 7, r2 = 0.9). Such a positive correlation is consistent with both 228Ra and [Nd] sharing a common benthic source. In contrast, benthic 228Th does not covary with [Nd], which might otherwise be expected if the source of deep Nd was delivered via top-down particle-based mechanisms.
Figure 2
Bottom water εNd does not covary with Ra isotopes (Figure 2B). While such a correlation is not obligatory under the benthic control model, work with Pacific sediments suggests that certain components of the sediments may be more reactive during diagenesis, and such reactivity is also likely element specific (Abbott et al., ; Du et al., ). Thus, we might expect that reaction of these sedimentary elements imposes a co-variation between the magnitude and isotopic composition of the flux.
We reiterate that there is undeniable evidence in support of the idea that seawater εNd traces basin scale modern water mass circulation, and the benthic control hypothesis suggested here does not detract or contradict these observations. On the contrary, we suggest that records of εNd may potentially be used to actually quantify current velocities, given that the “exposure time” of a water mass to the bottom (Abbott et al., ) implies a temporal sensitivity necessary to establish velocity. As an example, the bottom water data of Lambelet et al. () for the North Atlantic, can be plotted vs. latitude (Figure 3). A benthic control based interpretation of these data may be that there are three diagenetic regimes in the North Atlantic: north of 55°N, between 45 and 55°N and south of 40°N. We will first focus on the increasing [Nd] and decreasing εNd of the middle regime (~43°-55°N). In this region, there is a total change in εNd of ~2 units, of which conservative mixing with LSW can account for ~1 unit (Figure 4; using water properties defined in van de Flierdt et al. (). Applying the simple exposure time model of Abbott et al., ; assuming flux of 20 pmol/cm2/yr with εNd of −15, which is reasonable but arbitrarily chosen, and a deep water [Nd] of 15 pM, which is measured), we might estimate that the residual 1-unit corresponds to ~35 years of benthic flux exposure on a 1 km thick bottom water layer, occurring over ~12° latitude. This combination results in a bottom water current velocity of 0.2 cm/s. This estimate of effective transport rate, despite its obvious simplifications and assumptions, is not incomparable to the 1.8–2 cm/s estimates derived from CFC and tritum/He3 (Smethie et al., ; Steinfeldt and Rhein, ). This simple model scenario also predicts [Nd] of the water mass increases to 23 pM, similar to that observed (Figure 3). Other scenarios are obviously possible; for example, a 5 pmol/cm2/yr flux with εNd of −20 offers a similar 40 year exposure estimate, but predicts a final [Nd] of only 17 pM. The only existing pore waters from the Atlantic (Buzzard's Bay; Sholkovitz et al., ) show pore water [Nd] of >700 pM at ~13 cm depth in the sediment, which compares to fluxes as large as any seen in the Pacific (up to 30 pmol/cm2/yr; Haley and Klinkhammer, ; Abbott et al., ). Such a flux of 30 pmol/cm2/yr with an εNd of −25, as implied might be the case in this region (Lambelet et al., ), predicts only 4 years of exposure for ~1 εNd unit change, which equates to 1.1 cm/s velocity, but an increase of [Nd] to just 16 pM. Again, these calculations are obviously crude, and poorly constrained; we use them simply to illustrate sensitivity of deep water to benthic fluxes, and the feasibility of using the concept for quantitative reconstruction of circulation.
Figure 3
Figure 4
The Lambelet et al. () data show constancy of εNd in the Northwest Atlantic Bottom Water (NWABW) after it detaches (or is displaced) from contact with the bottom south of ~40°N. At these lower latitudes, Antarctic Bottom Water (AABW) is in contact with the sediments, and it is this southern-sourced water mass that appears to be modified during its transit northwards. While such a pattern can be explained with a benthic control model for [REE]/εNd, a “top-down” process that clearly alters bottom water but not the deep water overlying bottom water seems difficult to reconcile. Moreover, the potential density of AABW is remarkably stable along its flow path: σ = 27.9 kg/m3 at ~60°S (Stichel et al., ) vs. σ = 27.903 kg/m3 at ~37°N (from Lambelet et al., data; Figure 4). This constancy indicates that little water mass mixing occurs during this bottom water transit, yet εNd changes significantly: at least one unit in just the northern hemisphere data of Lambelet et al. (). Again, if we argue that conservation of NADW εNd precludes a top down driver for AABW εNd change, the logical conclusion is that benthic alteration impacts bottom water (i.e., AABW) εNd. Through the same approximations done above (assuming an arbitrarily chosen 5 pmol/cm2/yr flux with εNd of −15 and measured 30 pM [Nd] of bottom water), we might estimate that the 1-unit ΔεNd of AABW observed corresponds to ~150 years of exposure. Over 25° of latitude, this corresponds to a current velocity of 0.06 cm/s, implying that AABW has about one third the velocity as NWABW. (Note that because this calculation is based on a concentration gradient, the exact positions of where the slope is calculated can differ, making such estimations amenable to paleo-reconstructions where end-members may not be well constrained.) This effective transport rate estimate is, as was the case for our NADW estimate, lower than the CFC estimate of ~1.2 cm/s for AABW (Haine et al., ). We can predict that the benthic fluxes affecting AABW must be more complex, because linear extrapolation of the rate of change in εNd over latitude as shown in Figures 3, 4 would imply that AABW has an εNd of −6 at 60°S, which it does not (Stichel et al., ). The likely rationale for this difference is that the diagenetic fluxes along the greater AABW flow path are more nuanced in both their magnitude and composition than can be expressed by our simplified approximations.
There are three further points relevant to this discussion. First, the benthic control hypothesis does not invoke water mass mixing to explain patterns in εNd such as those shown in Figures 3, 4) That is, there is no obligate change in conservative properties, such as temperature or salinity, as the bottom water εNd is modified. Again, mid-depth water masses detached from the sediment will likely maintain their εNd signatures, or mix them as they would other conservative tracer properties along isopycnals, for example.
Second, the North Atlantic appears to have a propensity to mask the impacts of a benthic flux because ventilation rates are high, and thus exposure times are low. Instances where benthic influences may be directly observed in the North Atlantic (e.g., Lacan and Jeandel, ) would thus reflect rather extreme cases of high flux or divergent isotopic composition. In contrast, AABW in the South Atlantic (and PDW in the Pacific) circulate rather more slowly and thus would be prone to longer benthic exposure, resulting in more obvious impacts of a benthic source flux.
Third, the benthic flux mechanism seems to be intimately connected to authigenic metal oxides in the sediments (Abbott et al., ), which appear to modulate the REE flux and its εNd signature. Authigenic metal oxides appear to act as a “capacitor” of REE: i.e., both a source and sink of pore water and bottom water REE. Analogous to electrical capacitance, more abundant sedimentary metal oxides will produce greater Nd fluxes with less temporally variable εNd signatures (Abbott et al., ). While this mechanism is clearly far from resolved, it offers an explanation for the apparent lack of influence from the non-radiogenic lithogenic sediments of the central North Pacific (Jones et al., ). Specifically, we predict that there is no positive benthic flux from these sediments, or that reactive diagenetic materials are not present in these sediments (Abbott et al., ; Du et al., ), or both. We also note that the non-radiogenic Nd data presented by Jones et al. () reflect the post-leach lithogenic fraction of the sediments, which may not play a role at all in defining the benthic flux as defined by our proposed mechanism.
Paleoproxy Nd
Arguably the most important implication of benthic control on marine [REE]/εNd is how we interpret paleoceanographic records of water mass change. For these studies, the difference between a “top-down” or “bottom-up” control on bottom water and on the marine budget of REEs becomes paramount for two reasons:
Firstly, most interpretations of paleo-data treat εNd as fundamentally conservative, discounting the possibility of changes in REE supply that is not the result of water mass mixing or circulation. In a reversible scavenging model, this may be a valid assumption if deep-water sensitivity to surface ocean change is small. However, this supposition then begs the question of why reversible scavenging controls the modern deep [REE]/εNd? On the other hand, in a benthic control scenario the εNd of bottom water may evolve through non-mixing processes: such as a change in the benthic flux or a change in the “exposure time” of the bottom water to this flux (Abbott et al., ). As such, an εNd record interpreted through a benthic control model may reflect a change in benthic flux or exposure time, the latter of which may well be the result of changing circulation patterns or velocities. Despite the ostensible complications, we argue that interpreting εNd records purely as a conservative tracer discounts the “pseudo-conservative” nature of the tracer. Our benthic control model offers an internal cycling constraint on [Nd] and εNd compared to a reverse scavenging model. Both models require the ability to reconstruct fluxes through time to accurately use εNd as a tracer: the reversible scavenging model requires estimation of vertical particle fluxes and particulate type over time, whereas the benthic control model requires estimation of the benthic flux over time. Determining either of these fluxes will be non-trivial tasks, but may be resolvable.
Secondly, similarity between surface sediment archival phases (dispersed oxides, coatings, fish teeth, forams) and bottom water is often cited as evidence for record fidelity under the top down model for εNd. In the benthic control model, such correspondence will result from bottom water that is either dominated by a strong positive benthic flux, wherein the preformed bottom water εNd signature is overwhelmed by a pore water signature supported by the authigenic phases, or from a negative or net-zero pore water flux, such that pore water and authigenic phases passively record bottom water signatures. There is also a range of intermediate conditions where εNd of the archival phase may differ from bottom water, and where we would predict there is a significant, but not overwhelming, positive benthic flux into bottom water (e.g., Du et al., in the Pacific, or Gutjahr et al., ; Huck et al., in the Atlantic). The relationship between authigenic phases (used as archival records) and contemporaneous bottom water is far more complex in the benthic control model, because the authigenic phases associated with diagenesis are intimately associated with the nature of the benthic flux that impact the bottom water signal (Abbott et al., ). Moreover, the concentrations of sedimentary reactive particulate phases, pore water, and bottom water differ by orders of magnitude, so we stress the importance of looking at both εNd and the [REEs].
Our benthic model may give an initially grim outlook for paleoproxy work, but the situation is probably not at all intractable: in fact, these challenges are likely to be fully resolvable with further observations and models. Most pressing is the need to better determine the mechanisms and sources that define the Nd sink terms within the sediment column: these processes are ultimately what will define the proxy records. For example, Abbott et al. () did not discuss the deep (~10 cm) sink term seen in pore water [Nd] profiles; if this sink records contemporaneous pore water εNd, an apparent temporal offset could be generated in records of εNd from diagenetic Fe-Mn oxides and other isotope and trace element proxies archived in the calcite lattice of foraminifera (Piotrowski et al., ).
Theory and observations all indicate that some aspect of bottom water εNd is maintained in authigenic phases, and is likely identifiable in all sediment (Gutjahr et al., ; Wilson et al., ; Blaser et al., ; Du et al., ). Furthermore, the additional Nd generated through diagenesis should be far less variable and readily constrained as “baseline” variations in sedimentary records, vs. higher-order changes that are derived from seawater; analogous to Sr isotope records. A compilation of global surface sediment authigenic-phase εNd vs. bottom water εNd (Figure 5) illustrates that authigenic (archival) phases do trend with bottom water globally, but there is a remarkably uniform 1- εNd unit offset toward authigenic phases being more radiogenic: an offset equivalent in the Atlantic and the Pacific. Furthermore, we note that the Pacific data actually have a more robust correlation compared to the Atlantic data, consistent with the preceding discussion that benthic fluxes in the Atlantic are more difficult to identify from bottom water data. Bearing in mind that the concentration of Nd in the authigenic phases (at ppm Nd) is nominally six orders of magnitude greater than the bottom water (at pM Nd), it is difficult to explain these data without invoking a Nd source from within the sediments that is globally persistent. Of course, there are many outstanding questions to answer before better understanding the benthic flux of REEs and εNd: for example, can a flux be more non-radiogenic than the sediment host? How do these fluxes change over time?
Figure 5
In summary, while complex, diagenesis does follow predictable rules that make it amenable to modeling, which, in turn, can provide a robust estimation of past changes in the benthic control on bottom water Nd. In turn, these modeled boundary conditions can be used to calculate current direction and velocities as suggested in the modern examples above. Moreover, such studies may provide further evidence for changes in pore waters that can relate to changes in ocean redox or acid-base chemistry. The benthic model adds complexity, but also potential for the application of [REE]/εNd.
Beyond the REEs and εNd
The REEs and εNd are often touted as being excellent tracers of processes that impact marine geochemical cycles in a much broader sense. As such, the hypothesis of benthic control on deep ocean REE cycling may support suggestions for the importance of these diagenetic processes in marine geochemical cycles in general. For example, the same authigenic Fe/Mn oxides apparently involved in the marine cycling of [REE]/εNd are also rich in other trace metals (e.g., Cu, Zn, Ni, Co). These authigenic metal oxides may be acting in an expansive way as a capacitor for geochemical fluxes; i.e., as both a source and sink of these trace metals.
Statements
Author contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Acknowledgments
This work was supported by NSF grants OCE-1147407 to JM and BH, OCE-1357529 to BH, and OCE-1715106 to JM. We thank Drs. V. Goswami and T. van de Flierdt for their comments, and Catherine Jeandel for editorial handling of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2017.00426/full#supplementary-material
References
1
AbadieC.LacanF.RadicA.PradouxC.PoitrassonF. (2017). Iron isotopes reveal distinct dissolved iron sources and pathways in the intermediate versus deep Southern Ocean. Proc. Natl. Acad. Sci. U.S.A.114, 858–863. 10.1073/pnas.1603107114
2
AbbottA. N.HaleyB. A.McManusJ. (2015a). Bottoms up: sedimentary control of the deep North Pacific Ocean's εNd signature. Geology43, 1035–1035. 10.1130/G37114.1
3
AbbottA. N.HaleyB. A.McManusJ. (2016). The impact of sedimentary coatings on the diagenetic Nd flux. Earth Planet. Sci. Lett. 449, 217–227. 10.1016/j.epsl.2016.06.001
4
AbbottA. N.HaleyB. A.McManusJ.ReimersC. E. (2015b). The sedimentary flux of dissolved rare earth elements to the ocean. Geochim. Cosmochim. Acta154, 186–200. 10.1016/j.gca.2015.01.010
5
AkagiT. (2014). Diatoms spread a high εNd-signature in the North Pacific Ocean. Geochem. J.48, 121–131. 10.2343/geochemj.2.0292
6
AkagiT.FuF.HongoY.TakahashiK. (2011). Composition of rare earth elements in settling particles collected in the highly productive North Pacific Ocean and Bering Sea: implications for siliceous-matter dissolution kinetics and formation of two REE-enriched phases. Geochim. Cosmochim. Acta75, 4857–4876. 10.1016/j.gca.2011.06.001
7
ArsouzeT.DutayJ.-C.LacanF.JeandelC. (2009). Reconstructing the Nd oceanic cycle using a coupled dynamical – biogeochemical model. Biogeosciences6, 2829–2846. 10.5194/bg-6-2829-2009
8
BauM.TepeN.MohwinkelD. (2013). Siderophore-promoted transfer of rare earth elements and iron from volcanic ash into glacial meltwater, river and ocean water. Earth Planet. Sci. Lett. 364, 30–36. 10.1016/j.epsl.2013.01.002
9
BayonG.GermanC. R.BurtonK. W.NesbittR. W.RogersN. (2004). Sedimentary Fe–Mn oxyhydroxides as paleoceanographic archives and the role of aeolian flux in regulating oceanic dissolved REE. Earth Planet. Sci. Lett. 224, 477–492. 10.1016/j.epsl.2004.05.033
10
BayonG.ToucanneS.SkoniecznyC.AndréL.BermellS.CheronS.et al. (2015). Rare earth elements and neodymium isotopes in world river sediments revisited. Geochim. Cosmochim. Acta170, 17–38. 10.1016/j.gca.2015.08.001
11
BertramC.ElderfieldH. (1993). The geochemical balance of the rare earth elements and neodymium isotopes in the oceans. Geochim. Cosmochim. Acta57, 1957–1986.
12
von BlanckenburgF. (1999). Tracing past ocean circulation?Science286, 1862–1863. 10.1126/science.286.5446.1862b
13
BlaserP.LippoldJ.GutjahrM.FrankN.LinkJ. M.FrankM. (2016). Extracting foraminiferal seawater Nd isotope signatures from bulk deep sea sediment by chemical leaching. Chem. Geol. 439, 189–204. 10.1016/j.chemgeo.2016.06.024
14
ByrneR. H.KimK.-H. (1990). Rare earth element scavenging in seawater. Geochim. Cosmochim. Acta54, 2645–2656. 10.1016/0016-7037(90)90002-3
15
ByrneR. H.SholkovitzE. R. (1996). Chapter 158 Marine chemistry and geochemistry of the lanthanides. Handb. Phys. Chem. Rare Earths22, 497–593. 10.1016/S0168-1273(96)23009-0
16
CantrellK. J.ByrneR. H. (1987). Rare earth element complexation by carbonate and oxalate ions. Geochim. Cosmochim. Acta51, 597–605. 10.1016/0016-7037(87)90072-X
17
CarterP.VanceD.HillenbrandC. D.SmithJ. A.ShoosmithD. R. (2012). The neodymium isotopic composition of waters masses in the eastern Pacific sector of the Southern Ocean. Geochim. Cosmochim. Acta79, 41–59. 10.1016/j.gca.2011.11.034
18
CharetteM. A.MorrisP. J.HendersonP. B.MooreW. S. (2015). Radium isotope distributions during the US GEOTRACES North Atlantic cruises. Mar. Chem. 177, 184–195. 10.1016/j.marchem.2015.01.001
19
DavrancheM.PourretO.GruauG.DiaA.Le Coz-BouhnikM. (2005). Adsorption of REE(III)-humate complexes onto MnO2: experimental evidence for cerium anomaly and lanthanide tetrad effect suppression. Geochim. Cosmochim. Acta69, 4825–4835. 10.1016/j.gca.2005.06.005
20
DuJ.HaleyB.MixA. (2016). Neodymium Isotopes in authigenic phases, bottom waters, and detrital sediments in the Gulf of Alaska and their implications for paleo-circulation reconstruction. Geochim. Cosmochim. Acta193, 14–35. 10.1016/j.gca.2016.08.005
21
DunleaA. G.MurrayR. W.SauvageJ.SpivackA. J.HarrisR. N.D'HondtS. (2015). Dust, volcanic ash, and the evolution of the South Pacific Gyre through the Cenozoic. Paleoceanography30, 1078–1099. 10.1002/2015PA002829
22
ElderfieldH. (1988). The oceanic chemistry of the rare-earth elements [and discussion]. Philos. Trans. R. Soc. Lond. Math. Phys. Eng. Sci. 325, 105–126. 10.1098/rsta.1988.0046
23
ElderfieldH.GreavesM. J. (1982). The rare earth elements in seawater. Nature296, 214–219. 10.1038/296214a0
24
FrankM. (2002). Radiogenic isotopes: tracers of past ocean circulation and erosional input. Rev. Geophys. 40, 1-1–1-38. 10.1029/2000RG000094
25
FreslonN.BayonG.ToucanneS.BermellS.BollingerC.ChéronS.et al. (2014). Rare earth elements and neodymium isotopes in sedimentary organic matter. Geochim. Cosmochim. Acta140, 177–198. 10.1016/j.gca.2014.05.016
26
FrölljeH.PahnkeK.SchnetgerB.BrumsackH.-J.DulaiH.FitzsimmonsJ. N. (2016). Haawiian imprint on dissolved Nd and Ra isotopes and rare earth elements in the central North Pacific: local survey and seasonal variability. Geochim. Cosmochim. Acta189, 110–131. 10.1016/j.gca.2016.06.001.
27
Garcia-SolsonaE.JeandelC.LabatutM.LacanF.VanceD.ChavagnacV.et al. (2014). Rare earth elements and Nd isotopes tracing water mass mixing and particle-seawater interactions in the SE Atlantic. Geochim. Cosmochim. Acta125, 351–372. 10.1016/j.gca.2013.10.009
28
GermanC. R.ElderfieldH. (1990). Rare earth elements in the NW Indian Ocean. Geochim.Cosmochim. Acta54, 1929–1940.
29
GoldsteinS. J.JacobsenS. B. (1987). The Nd and Sr isotopic systematics of river-water dissolved material: Implications for the sources of Nd and Sr in seawater. Chem. Geol. Isot. Geosci. Sect. 66, 245–272. 10.1016/0168-9622(87)90045-5
30
GoldsteinS. J.JacobsenS. B. (1988). Nd and Sr isotopic systematics of river water suspended material: implications for crustal evolution. Earth Planet. Sci. Lett. 87, 249–265. 10.1016/0012-821X(88)90013-1
31
GoldsteinS. L.HemmingS. R. (2003). Long-lived isotopic tracers in oceanography, paleoceanography, and ice-sheet dynamics, in Treatise on Geochemistry, eds HollandH. D.TurekianK. K. (Pergamon; Oxford), 453–489.
32
GoswamiV.SinghS. K.BhushanR. (2014). Impact of water mass mixing and dust deposition on Nd concentration and eNd of the Arabian Sea water column. Geochim. Cosmochim. Acta145, 30–49. 10.1016/j.gca.2014.09.006
33
GrasseP.StichelT.StumpfR.StrammaL.FrankM. (2012). The distribution of neodymium isotopes and concentrations in the Eastern Equatorial Pacific: water mass advection versus particle exchange. Earth Planet. Sci. Lett. 353–354, 198–207. 10.1016/j.epsl.2012.07.044
34
GreavesM. J.ElderfieldH.SholkovitzE. R. (1999). Aeolian sources of rare earth elements to the Western Pacific Ocean. Mar. Chem.68, 31–38.
35
GrenierM.JeandelC.LacanF.VanceD.VenchiaruttiC.CrosA.et al. (2013). From the subtropics to the central equatorial Pacific Ocean: neodymium isotopic composition and rare earth element concentration variations. J. Geophys. Res. Oceans118, 592–618. 10.1029/2012JC008239
36
GutjahrM.FrankM.StirlingC. H.KeigwinL. D.HallidayA. N. (2008). Tracing the Nd isotope evolution of North Atlantic Deep and Intermediate Waters in the western North Atlantic since the Last Glacial Maximum from Blake Ridge sediments. Earth Planet. Sci. Lett. 266, 61–77. 10.1016/j.epsl.2007.10.037
37
GutjahrM.FrankM.StirlingC. H.KlemmV.van de FlierdtT.HallidayA. N. (2007). Reliable extraction of a deepwater trace metal isotope signal from Fe–Mn oxyhydroxide coatings of marine sediments. Chem. Geol. 242, 351–370. 10.1016/j.chemgeo.2007.03.021
38
HaineT. W. N.WatsonA. J.LiddicoatM. I.DicksonR. R. (1998). The flow of Antarctic bottom water to the southwest Indian Ocean estimated using CFCs. J. Geophys. Res. Oceans103, 27637–27653. 10.1029/98JC02476
39
HaleyB. A.FrankM.HathorneE.PisiasN. (2014). Biogeochemical implications from dissolved rare earth element and Nd isotope distributions in the Gulf of Alaska. Geochim. Cosmochim. Acta126, 455–474. 10.1016/j.gca.2013.11.012
40
HaleyB. A.KlinkhammerG. P. (2003). Complete separation of rare earth elements from small volume seawater samples by automated ion chromatography: method development and application to benthic flux. Mar. Chem. 82, 197–220. 10.1016/S0304-4203(03)00070-7
41
HaleyB. A.KlinkhammerG. P.McManusJ. (2004). Rare earth elements in pore waters of marine sediments. Geochim. Cosmochim. Acta68, 1265–1279. 10.1016/j.gca.2003.09.012
42
HoweJ. N. W.PiotrowskiA. M.OppoD. W.HuangK.-F.MulitzaS.ChiessiC. M.et al. (2016). Antarctic intermediate water circulation in the South Altantic over the past 25,000 years. Paleoc31, 1302–1314. 10.1002/2016PA002975
43
HuR.PiotrowskiA. M.BostockH. C.CrowhurstS.RennieV. (2016). Variability of neodymium isotopes associated with planktonic foraminifera in the Pacific Ocean during the Holocene and last glacial maximum. Earth Planet. Sci. Lett. 447, 130–138. 10.1016/j.epsl.2016.05.011
44
HuckC. E.van de FlierdtT.Jiménez-EspejoF. J.BohatyS. M.RöhlU.HammondS. J. (2016). Robustness of fossil fish teeth for seawater neodymium isotope reconstructions under variable redox conditions in an ancient shallow marine setting: robustness of fossil fish tooth nd. Geochem. Geophys. Geosyst.17, 679–698. 10.1002/2015GC006218
45
JacobsenS. B.WasserburgG. J. (1980). Sm-Nd isotopic evolution of chondrites. Earth Planet. Sci. Lett. 50, 139–155. 10.1016/0012-821X(80)90125-9
46
JeandelC. (2016). Overview of the mechanisms that could explain the ‘Boundary Exchange’ at the land–ocean contact. Philos. Trans. R. Soc. Math. Phys. Eng. Sci. 374:20150287. 10.1098/rsta.2015.0287
47
JeandelC.BishopJ. K.ZindlerA. (1995). Exchange of neodymium and its isotopes between seawater and small and large particles in the Sargasso Sea. Geochim. Cosmochim. Acta59, 535–547. 10.1016/0016-7037(94)00367-U
48
JeandelC.ThouronD.FieuxM. (1998). Concentrations and isotopic compositions of neodymium in the eastern Indian Ocean and Indonesian straits. Geochim. Cosmochim. Acta62, 2597–2607. 10.1016/S0016-7037(98)00169-0
49
JohannessonK. H.BurdigeD. J. (2007). Balancing the global oceanic neodymium budget: Evaluating the role of groundwater. Earth Planet. Sci. Lett. 253, 129–142. 10.1016/j.epsl.2006.10.021
50
JohannessonK. H.ChevisD. A.BurdigeD. J.CableJ. E.MartinJ. B.RoyM. (2011). Submarine groundwater discharge is an important net source of light and middle REEs to coastal waters of the Indian River Lagoon, Florida, USA. Geochim. Cosmochim. Acta75, 825–843. 10.1016/j.gca.2010.11.005
51
JohannessonK. H.PalmoreC. D.FackrellJ.ProutyN. G.SwarzenskiP. W.ChevisD. A.et al. (2017). Rare earth element behavior during groundwater–seawater mixing along the Kona Coast of Hawaii. Geochim. Cosmochim. Acta198, 229–258. 10.1016/j.gca.2016.11.009
52
JonesC. E.HallidayA. N.ReaD. K.OwenR. M. (1994). Neodymium isotopic variations in North Pacific modern silicate sediment and the insignificance of detrital REE contributions to seawater. Earth Planet. Sci. Lett. 127, 55–66. 10.1016/0012-821X(94)90197-X
53
JonesK. M.KhatiwalaS. P.GoldsteinS. L.HemmingS. R.van de FlierdtT. (2008). Modeling the distribution of Nd isotopes in the oceans using an ocean general circulation model. Earth Planet. Sci. Lett. 272, 610–619. 10.1016/j.epsl.2008.05.027
54
KuT.-L.LuoS. (1994). New appraisal of radium 226 as a large-scale oceanic mixing tracer. J. Geophys. Res. Oceans99, 10255–10273. 10.1029/94JC00089
55
LacanF.JeandelC. (2001). Tracing Papua New Guinea imprint on the central equatorial Pacific Ocean using neodymium isotopic compositions and rare earth element patterns. Earth Planet. Sci. Lett. 186, 497–512. 10.1016/S0012-821X(01)00263-1
56
LacanF.JeandelC. (2004a). Subpolar Mode Water formation traced by neodymium isotopic composition. Geophys. Res. Lett. 31:L14306. 10.1029/2004GL019747
57
LacanF.JeandelC. (2004b). Neodymium isotopic composition and rare earth element concentrations in the deep and intermediate Nordic Seas: constraints on the Iceland Scotland overflow water signature. Geochem. Geophys. Geosystems5:Q11006. 10.1029/2004GC000742
58
LacanF.JeandelC. (2004c). Denmark Strait water circulation traced by heterogeneity in neodymium isotopic compositions. Deep Sea Res. I Oceanogr. Res. Pap. 51, 71–82. 10.1016/j.dsr.2003.09.006
59
LacanF.JeandelC. (2005a). Acquisition of the neodymium isotopic composition of the North Atlantic Deep Water. Geochem. Geophys. Geosyst.6:Q12008. 10.1029/2005GC000956
60
LacanF.JeandelC. (2005b). Neodymium isotopes as a new tool for quantifying exchange fluxes at the continent–ocean interface. Earth Planet. Sci. Lett. 232, 245–257. 10.1016/j.epsl.2005.01.004
61
LambeletM.van de FlierdtT.CrocketK.RehkämperM.KreissigK.ColesB.et al. (2016). Neodymium isotopic composition and concentration in the western North Atlantic Ocean: results from the GEOTRACES GA02 section. Geochim. Cosmochim. Acta177, 1–29. 10.1016/j.gca.2015.12.019
62
MartinJ. H.KnauerG. A.KarlD. M.BroenkowW. W. (1987). VERTEX: carbon cycling in the northeast Pacific. Deep Sea Res. Oceanogr. Res. Pap.34, 267–285.
63
MawjiE.SchlitzerR.DodasE. M.AbadieC.AbouchamiW.AndersonR. F.et al. (2015). The GEOTRACES Intermediate Data Product 2014. Mar. Chem. 177, 1–8. 10.1016/j.marchem.2015.04.005
64
MilleroF. J. (1992). Stability constants for the formation of rare earth-inorganic complexes as a function of ionic strength. Geochim. Cosmochim. Acta56, 3123–3132. 10.1016/0016-7037(92)90293-R
65
MorelF. M. M.HeringJ. G. (1993). Principles and Applications of Aquatic Chemistry, 1st Edn. New York, NY: Wiley-Interscience.
66
NgwenyaB. T.MagennisM.OliveV.MosselmansJ. F. W.EllamR. M. (2010). Discrete Site surface complexation constants for lanthanide adsorption to bacteria as determined by experiments and linear free energy relationships. Environ. Sci. Technol. 44, 650–656. 10.1021/es9014234
67
OsborneA. H.HaleyB. A.HathorneE. C.FlögelS.FrankM. (2014). Neodymium isotopes and concentrations in Caribbean seawater: Tracing water mass mixing and continental input in a semi-enclosed ocean basin. Earth Planet. Sci. Lett. 406, 174–186. 10.1016/j.epsl.2014.09.011
68
PiepgrasD. J.WasserburgG. J. (1980). Neodymium isotopic variations in seawater. Earth Planet. Sci. Lett. 50, 128–138. 10.1016/0012-821X(80)90124-7
69
PiepgrasD. J.WasserburgG. J. (1987). Rare earth element transport in the western North Atlantic inferred from Nd isotopic observations. Geochim. Cosmochim. Acta51, 1257–1271. 10.1016/0016-7037(87)90217-1
70
PiotrowskiA. M.GoldsteinS. L.HemmingS. R.FairbanksR. G. (2004). Intensification and variability of ocean thermohaline circulation through the last deglaciation. Earth Planet. Sci. Lett. 225, 205–220. 10.1016/j.epsl.2004.06.002
71
PiotrowskiA. M.GoldsteinS. L.HemmingS. R.FairbanksR. G. (2005). Temporal relationships of carbon cycling and ocean circulation at glacial boundaries. Science307, 1933–1938. 10.1126/science.1104883
72
PourretO.DavrancheM.GruauG.DiaA. (2007). Organic complexation of rare earth elements in natural waters: Evaluating model calculations from ultrafiltration data. Geochim. Cosmochim. Acta71, 2718–2735. 10.1016/j.gca.2007.04.001
73
RempferJ.StockerT. F.JoosF.DutayJ.-C.SiddallM. (2011). Modelling Nd-isotopes with a coarse resolution ocean circulation model: Sensitivities to model parameters and source/sink distributions. Geochim. Cosmochim. Acta75, 5927–5950. 10.1016/j.gca.2011.07.044
74
RickliJ.FrankM.BakerA. R.AciegoS.de SouzaG.GeorgR. B.et al. (2010). Hafnium and neodymium isotopes in surface waters of the eastern Atlantic Ocean: Implications for sources and inputs of trace metals to the ocean. Geochim. Cosmochim. Acta74, 540–557. 10.1016/j.gca.2009.10.006
75
RickliJ.FrankM.HallidayA. N. (2009). The hafnium–neodymium isotopic composition of Atlantic seawater. Earth Planet. Sci. Lett. 280, 118–127. 10.1016/j.epsl.2009.01.026
76
RickliJ.GutjahrM.VanceD.Fischer-GöddeM.HillenbrandC.-D.KuhnG. (2014). Neodymium and hafnium boundary contributions to seawater along the West Antarctic continental margin. Earth Planet. Sci. Lett. 394, 99–110. 10.1016/j.epsl.2014.03.008
77
RobertsN. L.PiotrowskiA. M. (2015). Radiogenic Nd isotope labeling of the northern NE Atlantic during MIS 2. Earth Planet. Sci. Lett. 423, 125–133. 10.1016/j.epsl.2015.05.011
78
RousseauT. C. C.SonkeJ. E.ChmeleffJ.van BeekP.SouhautM.BoaventuraG.et al. (2015). Rapid neodymium release to marine waters from lithogenic sediments in the Amazon estuary. Nat. Commun. 6:7592. 10.1038/ncomms8592
79
SchachtU.WallmannK.KutterolfS. (2010). The influence of volcanic ash alteration on the REE composition of marine pore waters. J. Geochem. Explor.106, 176–187. 10.1016/j.gexplo.2010.02.006
80
SchijfJ.ChristensonE. A.ByrneR. H. (2015). YREE scavenging in seawater: a new look at an old model. Mar. Chem. 177, Part 3, 460–471. 10.1016/j.marchem.2015.06.010
81
ShillerA. M.ChanE. W.JoungD. J.RedmondM. C. Kessler, J. D. (2017). Light rare earth element depletion during Deepwater Horizon blowout methanotrophy. Sci. Rep.7:10389. 10.1038/s41598-017-11060-z
82
SholkovitzE. R. (1989). Artifacts associated with the chemical leaching of sediments for rare-earth elements. Chem. Geol.77, 47–51.
83
SholkovitzE. R. (1993). The geochemistry of rare earth elements in the Amazon River estuary. Geochim. Cosmochim. Acta57, 2181–2190. 10.1016/0016-7037(93)90559-F
84
SholkovitzE. R.LandingW. M.LewisB. L. (1994). Ocean particle chemistry: the fractionation of rare earth elements between suspended particles and seawater. Geochim. Cosmochim. Acta58, 1567–1579. 10.1016/0016-7037(94)90559-2
85
SholkovitzE. R.PiepgrasD. J.JacobsenS. B. (1989). The pore water chemistry of rare earth elements in Buzzards Bay sediments. Geochim. Cosmochim. Acta53, 2847–2856. 10.1016/0016-7037(89)90162-2
86
SinghS. P.SinghS. K.GoswamiV.BhushanR.RaiV. K. (2012). Spatial distribution of dissolved neodymium and eNd in the Bay of Bengal: role of particulate matter and mixing of water masses. Geochim. Cosmochim. Acta, 94, 38–96, 10.1016/j.gca.2012.07.017.
87
SmethieW. M.FineR. A.PutzkaA.JonesE. P. (2000). Tracing the flow of North Atlantic Deep Water using chlorofluorocarbons. J. Geophys. Res. Oceans105, 14297–14323. 10.1029/1999JC900274
88
SmithR. M.MartellA. E. (1989). Critical Stability Constants. Boston, MA: Springer.
89
SpivackA. J.WasserburgG. J. (1988). Neodymium isotopic composition of the Mediterranean outflow and the eastern North Atlantic. Geochim. Cosmochim. Acta52, 2767–2773.
90
StanleyJ. K.ByrneR. H. (1990). The influence of solution chemistry on REE uptake by Ulva lactuca L. in seawater. Geochim. Cosmochim. Acta54, 1587–1595. 10.1016/0016-7037(90)90393-Y
91
SteinfeldtR.RheinM. (2004). Spreading velocities and dilution of North Atlantic Deep Water in the tropical Atlantic based on CFC time series. J. Geophys. Res. Oceans109:C03046. 10.1029/2003JC002050
92
StewartJ. A.GutjahrM.JamesR. H.AnandP.WilsonP. A. (2016). Influence of the Amazon River on the Nd isotope composition of deep water in the western equatorial Atlantic during the Oligocene–Miocene transition. Earth Planet. Sci. Lett. 454, 132–141. 10.1016/j.epsl.2016.08.037
93
StichelT.FrankM.RickliJ.HaleyB. A. (2012). The hafnium and neodymium isotope composition of seawater in the Atlantic sector of the Southern Ocean. Earth Planet. Sci. Lett.317–318, 282–294. 10.1016/j.epsl.2011.11.025
94
StichelT.HartmanA. E.DugganB.GoldsteinS. L.ScherH.PahnkeK. (2015). Separating biogeochemical cycling of neodymium from water mass mixing in the Eastern North Atlantic. Earth Planet. Sci. Lett. 412, 245–260. 10.1016/j.epsl.2014.12.008
95
StumpfR.FrankM.SchoenfeldJ.HaleyB. A. (2010). Late Quaternary variability of Mediterranean Outflow Water from radiogenic Nd and Pb isotopes. Quat. Sci. Rev. 29, 2462–2472, 10.1016/j.quatscirev.2010.06.021
96
SuessE. (1980). Particulate organic carbon flux in the oceans - surface productivity and oxygen utilization. Nature288, 260–263.
97
TachikawaK.ArsouzeT.BayonG.BoryA.ColinC.DutayJ.-C.et al. (2017). The large-scale evolution of neodymium isotopic composition in the global modern and Holocene ocean revealed from seawater and archive data. Chem. Geol. 457, 131–148. 10.1016/j.chemgeo.2017.03.018
98
TachikawaK.AthiasV.JeandelC. (2003). Neodymium budget in the modern ocean and paleo-oceanographic implications. J. Geophys. Res. Oceans108:3254. 10.1029/1999JC000285
99
TachikawaK.JeandelC.Roy-BarmanM. (1999b). A new approach to the Nd residence time in the ocean: the role of atmospheric inputs. Earth Planet. Sci. Lett. 170, 433–446.
100
TachikawaK.JeandelC.VangriesheimA.DupréB. (1999a). Distribution of rare earth elements and neodymium isotopes in suspended particles of the tropical Atlantic Ocean (EUMELI site). Deep Sea Res. Oceanogr. Res. Pap. 46, 733–755.
101
van de FlierdtT.FrankM. (2010). Neodymium isotopes in paleoceanography. Quat. Sci. Rev. 29, 2439–2441. 10.1016/j.quascirev.2010.07.001
102
van de FlierdtT.GriffithsA. M.LambeletM.LittleS. H.StichelT.WilsonD. J. (2016). Neodymium in the oceans: a global database, a regional comparison and implications for palaeoceanographic research. Phil. Trans. R. Soc. A374:20150293. 10.1098/rsta.2015.0293
103
WaterhouseA. F.MacKinnonJ. A.NashJ. D.AlfordM. H.KunzeE.SimmonsH. L.et al. (2014). Global patterns of diapycnal mixing from measurements of the turbulent dissipation rate. J. Phys. Oceanogr. 44, 1854–1872. 10.1175/JPO-D-13-0104.1
104
WilsonD. J.PiotrowskiA. M.GalyA.CleggJ. A. (2013). Reactivity of neodymium carriers in deep sea sediments: Implications for boundary exchange and paleoceanography. Geochim. Cosmochim. Acta109, 197–221. 10.1016/j.gca.2013.01.042
105
YangJ.HaleyB. (2016). The profile of the rare earth elements in the Canada Basin, Arctic Ocean. Geochem. Geophys. Geosystems17, 3241–3253. 10.1002/2016GC006412
106
ZhengX.-Y.PlancherelY.SaitoM. A.ScottP. M.HendersonG. M. (2016). Rare earth elements (REEs) in the tropical South Atlantic and quantitative deconvolution of their non-conservative behavior. Geochim. Cosmochim. Acta177, 217–237. 10.1016/j.gca.2016.01.018
Summary
Keywords
rare earth elements, neodymium isotopes, marine geochemistry, benthic flux
Citation
Haley BA, Du J, Abbott AN and McManus J (2017) The Impact of Benthic Processes on Rare Earth Element and Neodymium Isotope Distributions in the Oceans. Front. Mar. Sci. 4:426. doi: 10.3389/fmars.2017.00426
Received
02 October 2017
Accepted
12 December 2017
Published
21 December 2017
Volume
4 - 2017
Edited by
Catherine Jeandel, UMR5566 Laboratoire d'études en Géophysique et Océanographie Spatiales (LEGOS), France
Reviewed by
Tina Van De Flierdt, Imperial College London, United Kingdom; Vineet Goswami, Colorado State University, United States
Updates
Copyright
© 2017 Haley, Du, Abbott and McManus.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Brian A. Haley bhaley@coas.oregonstate.edu
This article was submitted to Marine Biogeochemistry, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.