Circulation timescales of Atlantic Water in the Arctic Ocean determined from anthropogenic radionuclides
- 1Laboratory of Ion Beam Physics, Institute for Particle Physics and Astrophysics, ETH Zürich, Zurich, Switzerland
- 2Environmental Physics, Institute of Biogeochemistry and Pollutant Dynamics, ETH Zürich, Zurich, Switzerland
- 3Bedford Institute of Oceanography, Fisheries and Oceans Canada, Dartmouth, NS, Canada
Correspondence: Anne-Marie Wefing (firstname.lastname@example.org) and Núria Casacuberta (email@example.com)
The inflow of Atlantic Water to the Arctic Ocean is a crucial determinant for the future trajectory of this ocean basin with regard to warming, loss of sea ice, and ocean acidification. Yet many details of the fate and circulation of these waters within the Arctic remain unclear. Here, we use the two long-lived anthropogenic radionuclides 129I and 236U together with two age models to constrain the pathways and circulation times of Atlantic Water in the surface (10–35 m depth) and in the mid-depth Atlantic layer (250–800 m depth). We thereby benefit from the unique time-dependent tagging of Atlantic Water by these two isotopes. In the surface layer, a binary mixing model yields tracer ages of Atlantic Water between 9–16 years in the Amundsen Basin, 12–17 years in the Fram Strait (East Greenland Current), and up to 20 years in the Canada Basin, reflecting the pathways of Atlantic Water through the Arctic and their exiting through the Fram Strait. In the mid-depth Atlantic layer (250–800 m), the transit time distribution (TTD) model yields mean ages in the central Arctic ranging between 15 and 55 years, while the mode ages representing the most probable ages of the TTD range between 3 and 30 years. The estimated mean ages are overall in good agreement with previous studies using artificial radionuclides or ventilation tracers. Although we find the overall flow to be dominated by advection, the shift in the mode age towards a younger age compared to the mean age also reflects the presence of a substantial amount of lateral mixing. For applications interested in how fast signals are transported into the Arctic's interior, the mode age appears to be a suitable measure. The short mode ages obtained in this study suggest that changes in the properties of Atlantic Water will quickly spread through the Arctic Ocean and can lead to relatively rapid changes throughout the upper water column in future years.
1.1 The role of Atlantic Water in the Arctic Ocean
The Arctic sea-ice extent in 2019 was one of the lowest in the satellite record but in line with the long-term trend of declining ice cover (e.g., Meredith et al., 2021). In addition to the main contribution coming from the solar heating of the surface mixed layer in summer (Carmack et al., 2015), recent work has also highlighted the contribution of the inflow of warm Atlantic Water (AW) through the Fram Strait to the observed sea-ice decrease (Polyakov et al., 2005, 2017; Årthun et al., 2019). This is especially relevant for the Eurasian Basin of the Arctic Ocean, and the process is now commonly known as “Atlantification”.
AW constitutes about 90 % of the total inflow of waters to the Arctic Ocean with the remainder composed of ∼ 9 % Pacific Water entering through the Bering Strait and ∼ 1 % freshwater that includes runoff and precipitation (e.g., Woodgate, 2013). Both Pacific Water and freshwater largely reside in the upper water column of the Arctic Ocean, including the Polar Mixed Layer and the upper halocline (hereafter referred to as the surface layer). AW enters the Arctic Ocean through either the Fram Strait or the Barents Sea, forming the Fram Strait Branch Water (FSBW) and Barents Sea Branch Water (BSBW), respectively. These two branches encounter one another in the St. Anna Trough and constitute the mid-depth Atlantic layer of the Arctic Ocean, also referred to as Arctic Atlantic Water (AAW), which is found between 250 and 800 m depth. This water mass is characterized by higher temperatures and densities compared to the overlying halocline waters (e.g., Rudels, 2015). Core FSBW is generally found in about 300–400 m depth and can be associated with a potential density of σΘ=27.91. The BSBW core characterized by σΘ=28 is located below, at about 800 m depth (Karcher et al., 2012). Both branches circulate through the Arctic Ocean mainly under a cyclonic regime (Fig. 1a), following the pathways of the Arctic Ocean Boundary Current (Mauldin et al., 2010; Rudels et al., 2012). However, this pattern might be subject to temporal variations (Karcher et al., 2012).
Despite ongoing work, pathways of AW circulation in the Arctic Ocean are still not well understood. Generally, AW can be traced by its T–S (temperature–salinity) properties, atmospherically introduced anthropogenic pollutants such as CFCs, or anthropogenic radionuclides introduced from atmospheric nuclear weapon tests as well as nuclear reprocessing plants (Woodgate, 2013). AW circulation times or water mass ages can only be obtained from transient tracers and are rather poorly constrained. Recent studies suggest circulation times on the order of 15–30 years from the Barents Sea opening through the Arctic Ocean and back to the Fram Strait (Smith et al., 2011; Karcher et al., 2011).
As a consequence of ongoing climate change, increased freshwater input from melting sea ice as well as changes in the Atlantic and Pacific inflows may lead to changes in the circulation pattern in future years (e.g., Woodgate, 2013). Acquiring a better understanding of circulation pathways and their temporal variations will be key to predicting the nature of future changes in the Arctic domain, examples being the propagation of increased AW temperature or the accumulation of anthropogenic CO2 (Terhaar et al., 2020a). Therefore, suitable tracers are required that unambiguously label AW in the Arctic Ocean and at the same time provide information about timescales associated with the circulation of AW to revise present circulation times and trace future changes.
In this study, we will investigate circulation pathways, mixing regimes as well as tracer ages of AW in the surface and Atlantic layers of the Arctic Ocean. This will be done using a novel approach that combines the two long-lived anthropogenic radionuclides 129I and 236U. Obtained ages will be put into the context of available literature data, and different approaches on how to estimate circulation timescales will be compared. Strengths and weaknesses especially of applying the transit time distribution model introduced in Sect. 1.3 to anthropogenic radionuclides will be discussed and implications for the Arctic Ocean will be highlighted.
1.2 129I and 236U as tracers of Atlantic Water in the Arctic Ocean
Several recent studies have employed the combination of two long-lived anthropogenic radionuclides, 129I and 236U, as tracers for the circulation of AW in the Arctic Ocean and the Fram Strait (Casacuberta et al., 2016, 2018; Wefing et al., 2019). In contrast to other established anthropogenic radionuclide tracers such as 137Cs, the half-life (t1∕2) of 236U is around 23 × 106 years and therefore of the same order of magnitude as for 129I with × 106 years. Thus, these two isotopes can be regarded as conservative tracers whose application as age tracers comes from their strongly time-dependent input function (Casacuberta et al., 2018). An important advantage of the two isotopes is their low sampling volume. Measurements of 129I and 236U require only around 250 mL and 3–5 L of seawater, respectively.
129I and 236U have been introduced to the marine environment from two sources: (i) liquid releases from nuclear fuel reprocessing plants (mainly 129I) and (ii) global fallout from atmospheric nuclear weapon tests (mainly 236U) (Casacuberta et al., 2016). In the Arctic Ocean, releases from the two European reprocessing plants, the Sellafield nuclear reprocessing plant, United Kingdom, and the La Hague nuclear fuel reprocessing plant, France, are the dominant source of 129I and 236U followed by the global fallout signal. The natural background of both radionuclides in the Arctic Ocean is several orders of magnitude lower compared to the inputs from anthropogenic activities (Casacuberta et al., 2016).
The liquid releases from the two European reprocessing plants are transported to the Arctic Ocean mainly by the Norwegian Coastal Current (NCC) (Edmonds et al., 2001; Gascard et al., 2004; Christl et al., 2015), which supplies the Arctic shelf seas and the surface layer of the Eurasian Basin (Rudels, 2015; Fig. 1a, dark green). After their circulation through the Arctic Ocean, waters from the surface layer largely exit it via the Fram Strait, forming the upper layer of the East Greenland Current.
The much smaller global fallout signal is carried into the Arctic Ocean by all currents. This includes the input of Pacific Water through the Bering Strait, which remains in the upper water column and dominates the water properties in the Canada and Makarov basins. The input from the Atlantic occurs primarily through the Norwegian Atlantic Current (NwAC). This current also carries a part of the reprocessing plant signal as a result of the partial mixing with the NCC that brings this signal northward from the reprocessing plants. Although only a small part of the NCC is actually mixed with the NwAC resulting in a strongly diluted signal, the radioisotope signal carried by the NwAC is still much larger compared to the pure global fallout signal carried by Pacific Water. This mixture then enters the Arctic Ocean via the FSBW and the BSBW (Casacuberta et al., 2018; Fig. 1a, red), supplying the mid-depth waters. Thus, AW is uniquely tagged with the reprocessing plant signal of 129I and 236U, permitting us to assess the pathways and circulation times of AW in the Arctic Ocean.
1.3 Using 129I and 236U as transient tracers to estimate circulation times
The quantitative use of 129I and 236U as transient tracers for the estimation of AW circulation times or water mass ages requires accurate knowledge of the tracer input functions. These input functions need to be constructed for the location where the tracer signal enters the study area. Consequently, this is the location where the tracer age is defined as zero. Different models can then be applied to determine the water mass age from the tracer measurements, such as a binary mixing or a transit time distribution (TTD) model. The former has been used with anthropogenic radionuclides in earlier studies (e.g., Smith et al., 1998, 2005, 2011; Christl et al., 2015; Wefing et al., 2019), where it has also been referred to as a “tracer age model” or “dual-tracer approach”. The TTD model has been applied widely in the context of ocean interior ventilation studies (e.g., Haine and Hall, 2002; Waugh et al., 2003; Tanhua et al., 2009; Smith et al., 2011; Stöven et al., 2015) and to determine the oceanic uptake of anthropogenic CO2 (e.g., Hall et al., 2002; Waugh et al., 2006; Tanhua et al., 2008; Khatiwala et al., 2009; Olsen et al., 2010; Khatiwala et al., 2013; Stöven and Tanhua, 2014; Stöven et al., 2016; He et al., 2018). In this study the binary mixing model will be applied to samples from the surface layer, whereas the TTD model will be used for the mid-depth Atlantic layer due to the following model characteristics.
In the binary mixing model, only two processes are considered. The first one is lateral advection that carries the tracer signal defined in the input function into the Arctic Ocean and out again without lateral mixing. The second one is mixing with water masses that contain a constant (low) background radionuclide concentration from global fallout (here, e.g., Pacific Water). The relationship between the concentrations of two tracers then permits the definition of a dilution factor and a single tracer age for each water parcel (Smith et al., 1998). The assumptions behind the binary mixing model are reasonable for the surface layer of the Arctic Ocean, as shown in Smith et al. (2011), given its strong lateral confinement. However, these assumptions generally do not hold for the mid-depth Atlantic layer and for most other parts of the world's oceans, since lateral mixing by, e.g., mesoscale processes is a strong feature of the ocean's flow.
In the TTD method this lateral mixing is explicitly considered in addition to the purely advective flow. This leads to a probability density function (PDF) of water mass ages, or a distribution of transit times. The TTD method is an inversion method that aims to solve for the flow characteristics of an advection–diffusion problem given information about the distribution of tracers containing time information. In order for this problem to remain solvable, the flow problem is often limited to 1D and the boundary conditions at the source have to be well known. This is typically given for atmospherically introduced transient tracers such as CFCs, SF6, or 3H with known input functions or natural radionuclides such as 39Ar or 14C that have relatively constant surface concentrations and are only subject to radioactive decay. An important limitation of the TTD method is that it is not straightforward to extend it to the consideration of the mixing of different endmembers, especially when the endmembers are hard to “unmix” (Haine and Hall, 2002). In the Arctic Ocean, only the mid-depth Atlantic layer can be considered to consist of one endmember since no mixing with Pacific Water is assumed to take place, justifying our use of a one endmember TTD model (e.g., Rudels, 2015). This is not the case for the surface layer; hence only the binary mixing model will be applied there.
To our knowledge only one study has applied the TTD model to determine lateral transit times of AW in the Arctic Ocean using anthropogenic radionuclides (Smith et al., 2011). These authors combined radionuclides from nuclear reprocessing (129I and 137Cs) with the atmospherically introduced tracer CFC-11. By introducing those radionuclides with a point-like source function in addition to the classically used CFCs, Smith et al. (2011) substantially expanded the application of the TTD to the lateral flow of AW. The combination of these tracers could be applied to the mid-depth Atlantic layer because the AW is largely isolated from the atmosphere in the eastern Norwegian Sea, proximal to the region where the radionuclides' input function is defined. Hence the input functions for both sets of tracers could be initialized at the same location. This permits the comparison of lateral circulation ages from anthropogenic radionuclides with CFC-derived ventilation ages.
In this study the two models presented above will be applied to revise circulation features and investigate circulation times of AW in the Arctic Ocean using the combination of 129I and 236U. The detailed methodological approach is outlined in the following section.
This study is based on 129I and 236U data from various expeditions in the Arctic Ocean and the Fram Strait between 2011 and 2016 (Fig. 1 and Table S1 in the Supplement). All seawater samples were collected from Niskin bottles during CTD (conductivity–temperature–depth) casts, chemically purified, and measured using accelerator mass spectrometry (AMS). Details on sample collection and chemical treatment can be found in Casacuberta et al. (2016, datasets PS78 2011, PS80 2012), Casacuberta et al. (2018, dataset PS94 2015), and Wefing et al. (2019, dataset PS100 2016). Samples from the Nansen Basin in 2013 (ACCACIA) were also treated according to the protocols in Wefing et al. (2019). All these samples were measured at the Laboratory of Ion Beam Physics at ETH Zurich. Samples collected in the Canada Basin in 2015 (HLY1502) were processed and measured at the A. E. Lalonde AMS Laboratory at the University of Ottawa (129I) and the Centro Nacional de Aceleradores in Sevilla and ETH Zurich (236U). For the scope of this study, we only considered samples from the surface layer and the mid-depth Atlantic layer as defined in Table 1.
2.2 129I and 236U input functions
Two different input functions were defined for 129I and 236U, respectively: one for the surface and one for the mid-depth Atlantic layer (Fig. 1b and c). Both input functions were reconstructed from 1901 to 2019 in the following way: for 1901 to 1949 only the natural background was considered, which was assumed to be zero for both isotopes. From 1950 to 1971, only the global fallout tracer signal was considered, and from 1971 to 2019, the signals from global fallout and from reprocessing plant releases were combined.
For the global fallout signal we assumed a temporally constant 129I concentration of 1 × 107 atoms L−1 but used a time-dependent 236U concentration. This is valid because the atmospheric weapon tests in the 1950s and 60s introduced significantly less 129I compared to 236U (Raisbeck and Yiou, 1999; Sakaguchi et al., 2009). The 236U global fallout signal was calculated from a diffusive global surface ocean model with 1∘ × 1∘ lateral resolution fed by the latitudinally averaged atmospheric depositional flux of 236U from weapon tests (Christl et al., 2015). To determine the relevant concentration for our tracer input function for the Arctic Ocean we used the model output for 0–50 m depth at a location of 74∘ N and 19∘ E. This location corresponds approximately to the region where the NwAC bifurcates into FSBW and BSBW (Fig. 1a). This is also the region where the reprocessing plant signal carried by the NCC is admixed to the NwAC. The modeled global fallout 236U concentration decreases from about 3.5 × 107 atoms L−1 in the 1960s to about 1 × 107 atoms L−1 in the early 1990s. For the input function after 1990 we used a fixed 236U concentration of 1 × 107 atoms L−1, which matches recent measurements in the North Pacific Ocean (Eigl et al., 2017).
For the inputs of 129I and 236U from the two reprocessing plants we used the respective reported or reconstructed radionuclide releases (new data to 2019 were kindly provided by ORANO/IRSN, Christophe Ray, Pascal Bailly du Bois, personal communication, 2018) to model 129I and 236U in the streams emanating from the two sites, following the approach of Christl et al. (2015). In the present study the two streams were treated separately and for each branch, a 3-year moving average was applied to account for mixing and water mass residence time in the North Sea. To combine the two streams and the global fallout signal into input functions for the Arctic Ocean, we adapted the approach of Casacuberta et al. (2018). They used measurements of 129I and 236U concentrations in the Fram Strait and Barents Sea in 2015 to determine the fraction of radionuclides coming from La Hague, Sellafield, and global fallout in each of the three branches entering the Arctic Ocean, i.e., NCC (referred to as Arctic Shelf Break Branch therein) and the FSBW and BSBW, the latter two formed from the NwAC. In this study we used the fractions identified for the NCC for the surface layer input function and the average of the fractions identified for FSBW and BSBW for the Atlantic layer input function. Note that the surface and Atlantic layer input functions of 129I are in reasonable agreement (± 10 %) with those used by Smith et al. (2011) (Fig. S1 in the Supplement). This proves the solidity of the input function as the two studies used a different approach combining the reprocessing plant releases and also defined the surface input function at a different location. The calculated 129I and 236U input functions (Fig. 1b and c, Table S2 in the Supplement) were then used in the binary mixing model and the TTD model to investigate circulation times.
2.3 Models to determine circulation timescales and mixing regimes
2.3.1 Binary mixing model
In our application of the binary mixing model for the surface layer of the Arctic Ocean, we assumed that water labeled with the time-dependent tracer signal entering the Arctic Ocean in the surface AW input function only mixes with water carrying the global fallout signal of 1 × 107 atoms L−1 for 129I and 236U (i.e., Pacific and Atlantic Water not labeled with reprocessing plant releases).
In order to estimate circulation times or tracer ages (both terms have been used synonymously), mixing lines between each year of the input function and the global fallout endmember were constructed in 236U vs. 129I tracer space (Fig. 2a). The tracer age, as well as the corresponding dilution factor (DF) of the input function, were constrained by plotting the measured data in the same tracer space. It is important to note that dilution here corresponds to a dilution with global fallout background. Other dilutions, such as the addition of waters from sea-ice melt and meteoric waters in the surface layer, were not considered. Those contain a negligible radionuclide signal (Casacuberta et al., 2016, 2018) and therefore would represent a third endmember in the 129I and 236U tracer space, albeit one that is close to the global fallout endmember. Thus, even in the presence of a substantial input of freshwater, the effect would be small, justifying our neglect of this input.
The mixing lines converge towards the constant global fallout endmember; hence the determination of a tracer age becomes less precise for a strong dilution of tracer-labeled water with water containing only the global fallout signal. As the concept of the binary mixing model is based on the assumption of purely advective flow, tracer signals from different years are assumed not to mix with one another; i.e., only mixing between the input function signal from a single year with the constant global fallout signal is considered.
2.3.2 Transit time distribution model
In the TTD model that we applied to the study of the mid-depth Atlantic layer, the concentration of a stable tracer at sampling location x and time t is described by
The symbol c0(t) represents the tracer input function and G(x,t) is the Green's function which propagates the tracer signal and therefore describes the properties of the flow. One can consider G(x,t) as the PDF that weighs the tracer signals from different years of the input function for a certain sampling time t at location x. For one-dimensional transport described by the 1D advection–diffusion equation
(tracer concentration c, velocity u, diffusivity k), the PDF has the shape of an inverse Gaussian function:
Substituting and yields
with the two parameters Γ and Δ (Haine and Hall, 2002, both given in years) defining the shape of the TTD. This has been referred to as the inverse Gaussian TTD (IG-TTD), which is essentially a distribution of circulation times or ages. Γ is the first moment of the PDF and represents the mean age. Δ is related to the second moment of the PDF and is a measure of the width of the PDF; i.e., it describes how much a tracer signal disperses during the flow as a result of lateral mixing. The larger this mixing, the more incorporation of waters with different ages occurs within a water parcel sampled downstream from the initialization point of the input function. For Δ=0, the flow is purely advective. In this case, Γ is equal to the tracer age from the binary mixing model at a location x=uΓ. The ratio Δ∕Γ (the Péclet number) is a measure of the advectiveness of the flow (with a smaller Δ∕Γ indicating greater advective flow), and typical values range between 0.4–1.8 (e.g., Stöven et al., 2015, for the Fram Strait). Δ∕Γ is often set to the value of 1, reflecting an equal contribution from advection and lateral diffusion. This value was found to be a good representation in many parts of the world's oceans and reduces the number of free parameters to one (e.g., Waugh et al., 2004).
In the TTD model, we used the AW input function for the mid-depth Atlantic layer including the entire time range from 1901 up to the sample collection year. In order to determine the two moments of the TTD, Δ and Γ isolines were plotted in the 236U vs. 129I tracer space (Fig. 2b). More precisely, for a certain Γ isoline, Γ was fixed and c(t) was calculated for 129I and 236U concentrations according to Eq. (1) using the 129I and 236U AW input functions for c0(t) and varying Δ such that . Plotting the sample data in the same plot allowed for the identification of both Δ and Γ simultaneously; i.e., we did not assume a fixed Δ∕Γ ratio. Using the TTD model with the 129I–236U tracer pair can provide reliable results for Γ≈ 10–70 years and Δ≈ 10–50 years (Fig. 2b). For lower values of both parameters, the isolines intersect due to the shape of the input functions, which are not monotonically increasing or decreasing. For higher values, the isolines converge, which also translates into greater uncertainties.
In addition to the mean age Γ, we also considered the mode age tmode. This age is the circulation time with the highest probability within the PDF and has therefore also been referred to as the “most probable” age (Smith et al., 2011). It is given by
For a fixed Γ, the absolute difference between mode and mean age increases with an increasing Δ∕Γ ratio, i.e., with the amount of lateral mixing involved. Note that in our application of the TTD model, we considered tracer input from a single source only and thus neglected the mixing of tracer signals from different sources (which is the principle point of the binary mixing model).
3.1 Pan-Arctic distribution of 129I and 236U
The compilation of 129I and 236U concentrations presented in this study (data in Table S1) allows for a synoptic view of these two tracers across the Arctic Ocean and Fram Strait. Results will be discussed for the surface layer (10–35 m depth, Fig. 3a and b) and the Atlantic layer (represented by samples from 250–300 m depth, Fig. 3c and d) (see also Table 1). The data points in Fig. 3 represent single measurements, i.e., individual stations and depths.
In the surface layer of the Arctic Ocean, the highest 129I concentrations (> 6000 × 106 atoms L−1) were found in the Amundsen Basin in 2015 and the core Fram Strait outflow along the Greenland shelf break in 2016. These locations reflect the pathway for Arctic surface waters originating from the NCC that carries the highest 129I tracer signal into the Arctic Ocean (Casacuberta et al., 2018). In 2015, high 129I concentrations were clearly restricted to the Eurasian Basin and the sharp decline across the Lomonosov Ridge towards the Makarov Basin implies the recirculation of AW within the Amundsen Basin. In contrast, 129I concentrations decreased gradually towards the Makarov and Canada basins in 2011/12, indicating that AW was probably transported further into the Makarov Basin at that time. Lowest 129I concentrations in the surface layer (around 100 × 106 atoms L−1) were observed in the Canada Basin, suggesting the minimal presence of Atlantic-origin waters. 129I concentrations in the Amundsen Basin in 2015/16 samples (Fig. 3a – diamonds) were generally slightly higher compared to those from 2011/12 (Fig. 3a – circles), which might mirror the shape of the 129I input function which had a sharp increase during the 1990s.
The distribution of 236U in the surface layer was generally similar to that of 129I. Highest 236U concentrations (> 25 × 106 atoms L−1) were found in the Makarov Basin and the Amundsen Basin in 2011/12, followed by the Fram Strait outflow waters in 2016 (about 20 × 106 atoms L−1). Again, this reflects the presence of AW transported from the NCC carrying the high 236U tracer signal of the reprocessing plants. The inflowing FSBW and the Nansen Basin showed lower 236U between 12–15 × 106 atoms L−1, which is consistent with the 236U input function for FSBW in the Atlantic layer (Fig. 1c). As for 129I, samples taken in the Canada Basin exhibited the lowest 236U concentrations within the surface layer of 8 × 106 atoms L−1.
At around 300 m depth in the Atlantic layer, 129I concentrations in the Eurasian Basin and the Fram Strait were significantly lower compared to the surface layer, with highest concentrations of about 3600 × 106 atoms L−1 observed in the Nansen Basin. In the Canada Basin and Makarov Basin, however, 129I concentrations were higher compared to the surface layer, around 1000–3000 × 106 atoms L−1. A similar pattern was observed for 236U, whereby the highest 236U concentration reported in this study was measured in one sample taken at 300 m depth in the Canada Basin in 2015 (around 45 × 106 atoms L−1). In the rest of the Arctic Ocean, all 236U concentrations at the same depth level were below 24 × 106 atoms L−1, with highest concentrations detected in the Makarov Basin and the Fram Strait core outflow.
For the Canada Basin, the extremely high 236U concentrations measured between 250–300 m depth suggest that minimal dilution of the Atlantic layer input occurred during transport. In this regard, it should be noted that the 236U input function used in this study is largely based on a reconstruction of releases from the reprocessing plants, since no discharge data for 236U are available for time periods prior to the 1990s (Christl et al., 2015). A recent study used shell archives to reconstruct the time series for early 236U releases from Sellafield and La Hague (Castrillejo et al., 2020). They found higher 236U releases (by an order of magnitude) from Sellafield during the 1970s compared to those reconstructed by Christl et al. (2015). This finding suggests that the early component of the 236U input function may have been underestimated, which could explain the high 236U concentrations measured in the intermediate water samples from the Canada Basin that are presumed to be older than intermediate water samples from other basins.
The differences observed between tracer distributions in the surface and the Atlantic layer reflect the differing pathways of AW between the two layers: in the surface layer, AW is largely restricted to the Eurasian Basin, whereas in the mid-depth Atlantic layer, it is also transported into the Canada Basin. This circulation pattern confirms hydrographic observations (e.g., Rudels, 2015) and highlights the strength of 129I and 236U as AW tracers.
3.2 Circulation times and dilution factors in the surface layer
The application of the binary mixing model to samples from the surface layer revealed tracer ages on the order of 10–20 years across the Arctic (Fig. 4a). Rather than plotting each individual sample on the binary mixing lines, samples from 2015 and 2016 were here grouped by basins and 129I and 236U concentrations were averaged to provide a general overview.
In the central Arctic, tracer ages ranged from 8 years in the Makarov Basin to 12 years in the Fram Strait. The oldest tracer age of 22 years was found in the Canada Basin (Fig. 4a). This station, however, had also the strongest dilution with the global fallout endmember (DF > 20 or less than 5 % of the input function signal); hence the determination of a tracer age was less precise. The same holds true for the Makarov Basin. Less dilution with global fallout was found in the Amundsen Basin and the Fram Strait (DF of 2–3), which allowed for a more robust determination of tracer ages. A consideration of uncertainties on tracer ages and dilution factors arising from uncertainties in the input function can be found in Appendix A.
To investigate the spatial and temporal pattern of surface water tracer ages in greater detail, ages were constrained for individual stations and for all sampling years (Fig. 4b). Here, samples with a dilution factor larger than 20 were considered to be primarily composed of Pacific-origin water and were not included in the subsequent analyses. Taking the results from all expeditions into account provides a better spatial coverage of the Arctic Ocean, including samples from the eastern Amundsen Basin, close to the Laptev Sea shelf.
An overall increase in tracer ages was observed from close to the Laptev Sea shelf across the Amundsen Basin and towards the Fram Strait outflow. This pattern reflects the transport of surface waters originating from the NCC. Dilution factors constrained for the Arctic Ocean and Fram Strait generally ranged between 2 and 6 (Fig. 4c). The lowest dilution of tracer-labeled Atlantic-origin surface waters with the global fallout background was observed for samples closest to the Laptev Sea, in the Amundsen Basin and in the core Fram Strait outflow (DF between 2 and 3). This pattern is consistent with advective transport of NCC-origin water across the Arctic Ocean. In the Nansen Basin, dilution factors were slightly higher (up to 6.5). This might be due to mixing of the NCC signal with mainly global fallout labeled AW from the NwAC transported into the Nansen Basin through the Fram Strait, as for instance suggested in Karcher et al. (2012).
Regarding differences between sampling years (2011/12 vs. 2015), a slight decrease in tracer ages over time could be observed for the central Arctic Ocean (Fig. 4b). In the Makarov Basin, dilution factors significantly increased from about 2 in 2011/12 to 10–20 in 2015 (Fig. 4c). This implies that the Pacific Water fraction increased over time and can be explained by a shift in the Atlantic–Pacific water front that generally determines the water mass provenance in the Makarov Basin. The alignment of the Atlantic–Pacific front across the Arctic Ocean is influenced by the prevailing atmospheric circulation (indicated by, e.g., the Arctic Oscillation (AO) index; Morison et al., 2012; Alkire et al., 2019, and references therein). However, no clear trend in the AO index was observed between 2011/12 and 2015 and hence the reason behind changes in the Pacific Water fraction in the Makarov Basin remains unclear. Tracer ages should not be affected by the water mass provenance but could change as a result of potentially changing AW pathways in the central Arctic Ocean.
In the Fram Strait, tracer ages of shelf break stations in the core of the EGC (East Greenland Current) were 12–13 years, whereas towards the coast, tracer ages increased up to 17 years (Fig. 4b). A similar pattern was observed for dilution factors, with less dilution in the stations along the shelf break (coinciding with lower tracer ages) and higher dilution factors closer to the coast of Greenland (Fig. 4c). This indicates that a recirculation of waters could have occurred on the shelf, as previously proposed in Wefing et al. (2019). Additionally, meltwaters from Greenland glaciers may have influenced the results from the coastal stations. Closer towards the coast, samples did not seem to dilute with the global fallout endmember but rather seemed to be influenced by a different endmember with a higher 236U concentration (see also Fig. 5 in Wefing et al., 2019), which led to increased tracer ages in the binary mixing model. Due to a lack of 129I and 236U data from Greenland glacier meltwater, however, we could not investigate this hypothesis further.
3.3 Circulation times and mixing regimes in the mid-depth Atlantic layer
The mean ages (Γ) inferred from the TTD method for the mid-depth Atlantic layer were overall larger than the tracer ages inferred from the binary mixing model for the surface layer (Fig. 5a). For an overview about the TTD parameters, stations were first grouped by basins as in Sect. 3.2. Here, all samples from the mid-depth Atlantic layer (250–800 m depth) were averaged. Uncertainties on Γ and Δ are discussed in Appendix A.
Mean ages were lowest (about 20 years) in the Makarov and Amundsen basins and higher (about 35 years) in the Nansen Basin and the Fram Strait outflow. The values of Δ tend to correlate positively with the mean age, with the lowest Δ of 9 years determined for the Makarov Basin, followed by the Amundsen Basin and the Fram Strait (15 years) and the Nansen Basin (32 years). An average Γ and Δ could not be constrained for the Canada Basin, as the average value from 250–800 m depth plots outside the TTD model domain (Fig. 5a). Therefore, the data point from 500 m depth was plotted instead, showing a mean age of about 30 years and a Δ of about 5 years at that depth. The inferred ratio of Δ∕Γ varied substantially across basins but always stayed below 1; i.e., it was indicative of a rather advective flow. The lowest value (0.5) of Δ∕Γ was found in the Makarov Basin and the Fram Strait outflow, followed by the Amundsen Basin (0.75). The highest ratio () was determined for the Nansen Basin. Mode ages (tmode) derived from Γ and Δ were lowest (5–7 years) in the Nansen and Amundsen basins and increased towards the Makarov Basin (11 years) and the Fram Strait (17 years).
From the Δ and Γ values determined for each basin the corresponding IG-TTDs were constructed according to Eq. (4) (Fig. 5b). The shape of these distributions and especially the span between mean and mode ages clearly illustrate the different flow regimes in the different Arctic basins: in the Nansen Basin, the distribution was highly asymmetric and the difference between mean and mode age was largest (almost 30 years). The long tail towards older ages indicated the presence of very old water parcels. In contrast, the Makarov Basin TTD was more symmetric and the two age measures differed by less than 10 years. This suggests that less mixing was involved in the transport of AW to the Makarov Basin.
Overall, the mean age distribution (Fig. 6a; here depicted at 250–300 m depth) of individual stations showed a similar pattern as the tracer ages in the surface layer. Ages increased along the expected AW flow in the mid-depth Atlantic layer (e.g., Rudels, 2015), with younger waters in the Eurasian Basin and older waters in the Canada Basin and the Fram Strait outflow. A similar pattern was also observed in the distribution of mode ages (Fig. 6c). Across the Lomonosov Ridge, mean ages were comparable but significant differences were evident in the mode ages as well as the Δ∕Γ ratios (Fig. 6b), showing greater advective flow in the Makarov Basin compared to the Amundsen Basin.
A high mean age and the highest Δ∕Γ ratio were found for one station close to the Laptev Sea shelf. This result was rather unexpected given the proximity of this station to the initialization point of the input functions and was due to the comparably low 129I and 236U concentrations measured at 300 m depth. A possible explanation could be a higher influence of the FSBW input function at that depth rather than the mixture of FSBW and BSBW that is used for our Atlantic layer input function. Using the FSBW input function instead, the mean age for this station would be on the order of 10 years and the Δ∕Γ ratio would be about 0.5, which is much more reasonable. In the context of this study, the Laptev Sea shelf data point will therefore not be discussed further.
In the Fram Strait, a similar pattern was observed in the mid-depth Atlantic layer as for the surface layer: lower mean ages and Δ∕Γ ratios around 0.5 along the shelf break indicate greater advective flow in the East Greenland Current core, whereas Γ and Δ∕Γ increased towards the coast of Greenland.
Considering individual samples from the entire Atlantic layer (from 250–800 m depth) on a transect across the Arctic Ocean (Fig. S2a and b in the Supplement), we observed that the high mean ages and Δ∕Γ ratios in the Nansen Basin were not limited to 300 m depth but were also obtained for samples from 500 m depth. Whereas mean ages and Δ∕Γ ratios in the Amundsen and Makarov basins did not change with depth, an increase in mean ages of up to more than 50 years at 800 m depth was observed in the Fram Strait. In the Canada Basin, mean ages and Δ∕Γ ratios could be determined for samples from 500 and 800 m depth and lay around 30 years and 0.1, respectively. This implies that across the Atlantic layer of the Arctic Ocean, Δ∕Γ ratios were lowest in the Canada Basin, suggesting very advective AW flow. Mode ages were more or less constant over depth, a slight increase was only observed in the Makarov Basin (Fig. S2c). Largest mode ages of around 30 years were obtained for the Canada Basin at 500 and 800 m depth.
In the first part of the discussion, the use of 236U (in combination with 129I) will be validated as a transient tracer to estimate circulation times in the Arctic Ocean. This will be done in the context of earlier studies that used different models and tracers. The results presented in this paper will be mainly compared to the study by Smith et al. (2011), which is also based on the use of two anthropogenic radionuclides (129I and 137Cs). The TTD model results determined in this work will additionally be compared to other studies conducted in the Arctic Ocean that used CFCs and SF6 as ventilation tracers. A second subsection of the discussion examines the use of the mode age as a measure of circulation times and conceptual differences between mean and mode ages will be discussed. The last part of the discussion highlights the implications of the obtained circulation times for a better understanding of AW circulation in the Arctic Ocean.
4.1 Circulation patterns of Atlantic Water in the Arctic Ocean revealed by 129I and 236U
4.1.1 Circulation times of Atlantic Water in the context of earlier studies
For the central Arctic Ocean, tracer ages and mean ages obtained from the combination of 129I and 236U in the mixing model and the TTD model, respectively, are overall in good agreement with those presented in Smith et al. (2011) (based on samples collected in the mid-1990s), especially in the Amundsen Basin (Table 2). For both, the surface layer and the mid-depth Atlantic layer, AW ages generally increased from the entrance of the Arctic Ocean through the Amundsen and Makarov basins towards the Fram Strait. This pattern is consistent with the general understanding of AW flow in the Arctic Ocean derived from hydrographic measurements (e.g., Rudels, 2015). 129I–236U mean ages for the Fram Strait are in line with mean ages observed slightly upstream (north of Greenland) by Smith et al. (2011). In the Fram Strait, the observed increase in mean ages from 300 to 800 m depth indicates that the mean transport was slower in the BSBW layer compared to the FSBW layer. This result, together with the fact that increased mean ages were not observed in the BSBW layer of the Amundsen Basin, supports the common perception that BSBW tends to follow the path of the Arctic Ocean Boundary Current through the Canada Basin, whereas FSBW instead returns to the Fram Strait within the Eurasian Basin, along the flanks of the Lomonosov Ridge (e.g., Rudels, 2015).
4.1.2 Dilution factors and the Atlantic–Pacific Water interface in the surface layer of the Makarov Basin
Dilution factors obtained from the binary mixing model for the surface layer reflect mixing between tracer-labeled AW and waters carrying the global fallout background. The results determined in this study agree with those reported in Smith et al. (2011) (Table 2). Note that dilution factors calculated by Smith et al. (2011) were referenced to 60∘ N, whereas in the present study, the input function was defined at the entrance to the Arctic Ocean. Hence, dilution factors from Smith et al. (2011) were corrected by the dilution factor of 3, which was determined for one station in the Barents Sea opening.
Low dilution factors in the Amundsen Basin and the core Fram Strait outflow suggest a confined transport of tracer-labeled waters originating from the NCC into the central Arctic Ocean and out through the Fram Strait. The higher dilution factors in the Nansen Basin indicate the presence of waters labeled only by the global fallout signal, either old AW or waters originating from the NwAC (note that dilution with Pacific Water is unlikely to happen in the Nansen Basin).
In the Makarov Basin, the global fallout signal is mainly transported by water of Pacific origin; hence the dilution factors can be interpreted as a measure of the Pacific Water fraction. From 2011/12 to 2015, a significant increase in dilution factors from 2 to 20 was observed in the Makarov Basin, hence indicating an increase in the Pacific Water fraction. Note that further downstream, a change in the spreading of Pacific Water could also be reflected in the composition of outflowing waters in the Fram Strait, but no time series of 129I and 236U data is available for the Fram Strait to date. The finding of a high Pacific Water fraction in the surface of the Makarov Basin in 2015 is supported by estimates from nutrient relationships (N : P ratios) from the same year (Alkire et al., 2019, “West leg” of the GEOTRACES GN01 expedition), which determined Pacific Water fractions greater than 70 %. New tracers that can be used to differentiate between Atlantic- and Pacific derived waters are a valuable tool because N : P ratios, used as the classical Atlantic–Pacific Water tracer, are known to be influenced by processes occurring along the transport, especially denitrification in the broad shelf seas of the Arctic Ocean (Bauch et al., 2011; Alkire et al., 2019). In contrast, 129I and 236U behave more conservatively than nutrients and are not affected by any biogeochemical processes taking place during transport, suggesting the potential of using these two anthropogenic radionuclides as water mass provenance tracers in the surface Arctic Ocean.
4.1.3 Lateral mixing within the Atlantic Water flow and its impact on mean ages
For the mid-depth Atlantic layer, along-flow lateral mixing taking place within the AW flow is described by the Δ∕Γ ratio and influences the mean ages obtained from the TTD model. Δ∕Γ ratios determined from 129I and 236U agree well with the findings from Smith et al. (2011) for all basins (Table 2). Overall, the AW flow was found to be rather advective, especially in the Makarov and Canada basins and the core of the Fram Strait outflow.
An exception is the Nansen Basin, where Δ∕Γ ratios ≥1 were determined for two stations (Fig. 6b), which coincided with high mean ages (Fig. 6a). The high mean ages in the Nansen Basin were an especially unexpected finding, as the stations are located close to the tracer sources in the Fram Strait and the Barents Sea (i.e., one would expect younger waters that have not undergone much lateral mixing). High Δ∕Γ ratios imply a broad TTD and significant mixing of waters with different ages. One possible explanation for the rather high Δ∕Γ ratios found in the Nansen Basin resulting in high mean ages is a mixture of inflowing AW in the FSBW with water that has been transported through the Arctic and re-entered the Nansen Basin as a result of recirculation north of the Fram Strait. In a similar process, the direct recirculation of northward-flowing FSBW within the Fram Strait might have led to the elevated Δ∕Γ ratios observed in the eastern- and southernmost stations of the Fram Strait, which are also accompanied by higher mean ages.
Two recent ventilation studies in the Arctic Ocean using the TTD model with atmospherically introduced transient tracers assumed a fixed (Tanhua et al., 2009; Stöven et al., 2016). Mean ages for the Nansen and Amundsen basins and the Fram Strait determined from these studies are in general agreement with the 129I–236U TTD model results (Table 2). Across the Lomonosov Ridge, however, Tanhua et al. (2009) observed a sharp front for all samples below 200 m depth with an increase in mean ages in the direction of the Makarov Basin, where mean ages were significantly higher compared to the present 129I–236U study. In this regard it should be noted that a recent model study suggested a general overestimation of CFC-derived mean ages, especially in waters ventilated in the Barents Sea, due to an undersaturation of pCFC-12 in the surface waters (Terhaar et al., 2020b).
In the Makarov Basin, the greatest discrepancy between the uniform assumed by Tanhua et al. (2009) and the present data (Δ∕Γ of 0.2–0.4) was observed. Consequently, the assumption of a uniform , which was found to be representative of many parts of the world's oceans, does not appear to be valid for the entire Arctic Ocean, especially for the Makarov Basin. Here the AW flow was found to be rather advective with less lateral mixing occurring compared to the Eurasian Basin. This result is also in line with the recent study by Rajasakaren et al. (2019), who constrained a Δ∕Γ of about 0.6 for the Amerasian Basin, based on the comparison between mean ages derived from CFC-12 and SF6. To provide the best estimate of mean ages, a combination of different tracers is required to constrain Δ∕Γ separately, which supports the findings from Smith et al. (2011).
Generally, some limitations should be considered with respect to the use of the TTD model in the Arctic Ocean and the Fram Strait in particular. In the Fram Strait, highly advective flow (low Δ∕Γ ratios, Fig. 6b) was observed for the East Greenland Current core in addition to slightly higher mean ages compared to those in the central Arctic (Fig. 6a). However, an important consideration for the Fram Strait is that it likely registers a mixture of Atlantic-origin water masses composed of those that circulated through the short Eurasian Basin loop and those that followed a pathway through the Canada Basin, i.e., the long loop. As Smith et al. (2011) noted, the IG-TTD is a unimodal distribution which assumes a single dominant pathway in the water transport. This does not hold for the Fram Strait and the results from the simple TTD approach applied in this study should therefore be regarded with some reservation. In a successive approach, a bimodal TTD which would capture both circulation loops could be considered for the Fram Strait (see, e.g., Stöven and Tanhua, 2014). The functional form of the Green's function G(t), describing the propagation of the tracer signal, introduces an uncertainty into the reported mean and mode ages. Examples of different functional forms have for instance been discussed in Haine and Hall (2002). The inverse Gaussian form of G(t), which solves the 1D advection–diffusion equation, has been shown to perform well for many oceanographic applications and is therefore the most widely used functional form.
4.2 Different age measures determined from the TTD model
4.2.1 Mode ages in the context of Arctic Ocean circulation
To our knowledge, mode ages of the TTD have only been considered in Smith et al. (2011), who suggested that the mode age might be a better comparison to the tracer age for pulse-like input functions. They only constrained the mode age for a single station, however, where they found it was similar to the tracer ages obtained from the binary mixing model.
A somewhat similar approach is the study by Mauldin et al. (2010), who used the combination of 3H and 3He to estimate the boundary current transport of BSBW. They applied a leaky pipe model with an advective core and considered mixing with the exterior region (based on Waugh and Hall, 2005) instead of the classical TTD model. The assumed distribution does not have the shape of an inverse Gaussian but peaks at the “advective time” tadv, which is calculated from the modeled current velocity of the core region. The overall shape of the distribution with a peak at low transit times and a long tail towards higher transit times is however similar to the IG-TTD, and as a first approximation, this advective time can be compared to the mode ages derived in the present study. The advective times reported by Mauldin et al. (2010) for the BSBW pathway along the marginal slopes of the Amundsen, Makarov, and Canada basins generally agree with the mode ages obtained in the present study (Table 3). For the Canada Basin, the 129I–236U mode ages are about 25 years higher (i.e., implying longer circulation times) compared to the Mauldin et al. (2010) study.
Here it should be noted that the latter employed data from the 1990s when the Arctic Oscillation was in a positive phase and exceptionally strong cyclonic boundary current conditions prevailed in the Arctic Ocean, associated with accelerated boundary current flow (Karcher et al., 2012). This implies low advective ages. The Canada Basin mode age obtained from 129I and 236U for 2015 cannot unambiguously be attributed to the boundary current. Additionally, AW pathways to the Canada Basin probably changed as a consequence of a weakening of the boundary current due to the transition of the AO to a more negative phase (Karcher et al., 2012).
Mode age results for the Arctic Ocean can also be compared to estimates of circulation times based on the propagation of hydrographic properties. Several studies have for example used temperature anomalies as tracers for AW flow (Swift et al., 1997; Rudels et al., 2000; Dmitrenko et al., 2008; Polyakov et al., 2011). Warm pulses were observed in inflowing AW on hydrographic sections or at moorings in the Fram Strait and subsequently tracked through the Arctic Ocean and re-analyzed at downstream sampling sites. Circulation times deduced from several studies are summarized in Table 3. In addition, model experiments (e.g., Karcher et al., 2011) suggest that warm water pulses entering the Arctic Ocean through the Fram Strait take about 15 years to recirculate back to the Fram Strait via the Eurasian Basin loop and about 20–30 years via the Canada Basin loop. These circulation times, both from observations and models, are overall within the same range as mode ages obtained from the combination of 129I and 236U.
4.2.2 Comparison of mean and mode ages
Mean and mode ages provide different circulation timescales and also impart a different understanding of the flow of AW in the mid-depth Atlantic layer. The mean age of the TTD is mainly influenced by the long tail of the distribution that extends towards old ages. The asymmetry depends on the Δ∕Γ ratio as a higher influence of mixing leads to a stronger dispersion of the TTD peak. Mean ages are classically used as the age measure derived from the TTD, especially in ocean ventilation studies that often aim at an estimation of the uptake of anthropogenic carbon (Cant) from the atmosphere. In this regard, the focus lies on the Cant inventory accumulating over time rather than the determination of tracer signal transport times. In contrast, mode ages reflecting the transit time with the highest probability may better reflect the propagation of biogeochemical pulses (e.g., temperature-salinity anomalies, plankton blooms).
In general, there is no unique definition of a circulation time once the effect of mixing is taken into account. With respect to water mass ages, tracer ages derived from the simple mixing model are a valid approximation if the mixing of ages within a water parcel can be assumed to be limited to a short time range, that is if the Δ∕Γ ratio is low. In this case, the different age measures deduced from the TTD are also similar and can be compared to tracer ages. If significant mixing is involved, only the TTD model provides meaningful age estimates. We found that overall, mode ages derived from the TTD model for the Atlantic layer are in the same range as tracer ages from the surface layer.
4.3 Future perspectives on the use of 129I and 236U in the Arctic Ocean
Climate change will continue triggering large-scale circulation changes in the Arctic Ocean, fostering the use of transient tracers which can unambiguously label AW and can provide information about circulation pathways and timescales. The present study highlights the combination of 129I and 236U in the Arctic Ocean as a powerful tool to understand such processes.
In addition to circulation times and mixing regimes obtained from tracer measurements, simulations using regional ocean models can predict changes in AW circulation patterns in the coming years. A numerical model used to simulate 129I transport through the Arctic Ocean, for example, revealed significant circulation changes during the mid-1990s and 2000s resulting from changing AO phases (Karcher et al., 2012). The inclusion of different tracers (e.g., 236U) into these models, together with a better data coverage, will help to improve our understanding of the Arctic Ocean circulation and ventilation as well as potential responses to climate change.
The properties of AW flow also play a key role in understanding and predicting the transport of heat or Cant to the Arctic Ocean. A recent study suggests that the amount of Cant that entered the Arctic Ocean through lateral transport, i.e., inflowing AW, is about 3 times larger compared to what entered through air–sea gas exchange (Terhaar et al., 2019). Latest estimates predict a Cant inventory of about 9 Pg carbon at the end of the 21st century, which exceeds projections from earlier studies (Terhaar et al., 2020a). The consequent decrease in calcium carbonate saturation states is suggested to occur in the mid-depth Atlantic layer of the Arctic Ocean rather than the surface waters, in depths between 400–800 m, according to Terhaar et al. (2020a). A study by Ulfsbo et al. (2018) found a strong linear correlation between measurements of Cant and 129I in the Atlantic layer of the Eurasian Basin in 2011/12 and 2015, also pointing to AW as the main source of both substances. Circulation times obtained from the combination of 129I and 236U can therefore help understanding the fate of Cant and predict its propagation into the different Arctic basins.
129I and 236U measurements on samples collected in the Arctic Ocean and Fram Strait between 2011 and 2016 were used to investigate circulation timescales and flow characteristics of AW in the Arctic Ocean. Two different models were utilized: a binary mixing model for the surface layer and the TTD model for Atlantic Water in the mid-depth Atlantic layer.
The binary mixing model, which takes into account the mixing between AW labeled with a radionuclide signal from European nuclear reprocessing plants and waters carrying background concentrations from global fallout, provided AW tracer ages which support findings from previous studies (summarized in Fig. 7, green). This is a useful finding as it is problematic to use atmospherically introduced ventilation tracers to study circulation timescales of surface waters which are constantly exposed to atmospheric gas exchanges. The 129I and 236U tracer pair has been shown to be conservative with a well-defined, geographically specific input function that incidentally provides a method for the accurate resolution of Atlantic-origin from Pacific-origin water in the surface layer of the Makarov Basin. For this location, dilution factors obtained from the mixing model suggested an increase in the Pacific Water fraction between 2011/12 and 2015.
In the mid-depth Atlantic layer (250–800 m depth), the TTD model was used to evaluate the flow of AW through the major Arctic basins. Mean ages (Γ, Fig. 7, red) inferred from the TTD model were in good agreement with those obtained from ventilation tracer studies for Arctic Ocean basins where Δ∕Γ ratios from 129I and 236U matched the assumed from other studies. This assumption did not appear to be valid for the entire Arctic Ocean, however, especially in the Makarov and Canada basins, where Δ∕Γ ratios were significantly smaller than 1.
Mode ages, which correspond to the most probable age of the TTD, were in the same range as circulation times deduced from the propagation of temperature anomalies (tmode, Fig. 7, red). For the lateral propagation of tracer signals, especially those of a pulsed nature, the mode age may be a better age measure compared to the mean age, as it is not influenced by the asymmetric shape of the TTD and the long tail towards old ages. For a comparison between tracer ages, mean, and mode ages, the different underlying models for the surface and the Atlantic layers need to be considered. In the case of lateral mixing, tracer ages are not meaningful and mean ages will always exceed mode ages. Here, the choice of the age estimate generally depends on the context.
Overall, AW pathways and circulation times constrained from the combination of 129I and 236U add to a better understanding of the present state of Arctic Ocean circulation. The Arctic is considered a region to experience major changes as a consequence of global warming in future years, including higher freshwater input due to sea-ice melt, increased heat transported into the Arctic Ocean by AW, or ocean acidification due to the uptake of anthropogenic carbon. Hence transient tracers labeling AW are valuable tools to track the future evolution of AW circulation in the Arctic Ocean.
Uncertainties for basin-averaged tracer ages and dilution factors as well as the TTD parameters result from both the standard deviation of the averaged samples (error bars in Figs. 4a and 5) and from uncertainties in the 129I and 236U input functions (Fig. 1b and c, about 10 %–20 % for the last 30 years). Due to the irregular shape of the input functions the resulting uncertainties are highly asymmetric.
To estimate the influence of input function uncertainties, a sensitivity analysis was performed as an example for the Amundsen Basin average from 2015. In the binary mixing model (surface layer, tracer age of 10 years and dilution factor of 2.2), the uncertainty in the 236U input function was the main driver for tracer age uncertainties (Fig. A1a), whereas dilution factors were mainly affected by uncertainties in the 129I input function (Fig. A1b). In 2005 (corresponding to tracer age of 10 years for samples from 2015), the 129I and 236U surface layer input functions had associated uncertainties of about 20 % and 10 %, respectively. This translated to approx. years uncertainty in the tracer age and in the dilution factor.
In the TTD model (mid-depth Atlantic layer, Γ = 20 years, Δ = 15 years), the uncertainty in the 129I input function was the main driver for uncertainties in Γ (Fig. A1c) and the uncertainty in the 236U input function was the main driver for uncertainties in Δ (Fig. A1d). The actual Atlantic layer input function uncertainties of around 10 % for 129I and 5 % for 236U resulted in uncertainties of the order of ± 5 years for both Γ and Δ.
Data used in this study can be found in Table S1 in the Supplement.
The supplement related to this article is available online at: https://doi.org/10.5194/os-17-111-2021-supplement.
AMW performed the conceptualization, investigation, formal analysis, and wrote the original draft. NC supported the conceptualization, supervised the study, acquired funding, and performed writing (review and editing). MC supported the investigation and formal analysis, and performed writing (review and editing). NG supervised and performed writing (review and editing). JNS supported the conceptualization and performed writing (review and editing).
The authors declare that they have no conflict of interest.
The authors acknowledge all people involved in the sampling activities across the Arctic Ocean and the subsequent sample processing and measurements, providing the available 129I and 236U data across the Arctic Ocean. Special thanks go to Mercedes López-Lora and Elena Chamizo Calvo from CNA Sevilla and Tim Kenna from the Lamont-Doherty Earth Observatory for sampling and processing 236U samples from the USCGC Healy cruise in 2015. We also thank Gideon Henderson from the University of Oxford, who provided seawater samples from the ACCACIA cruise in 2013. Anne-Marie Wefing received funding from the ETH Doctoral Grant ETH-06 16-1 “Combining U-236 with other multi-sourced anthropogenic tracers as a new tool to study ocean circulation”. Núria Casacuberta's research was supported by the Swiss National Science Foundation (AMBIZIONE PZ00P2_154805) and the Laboratory of Ion Beam Physics at ETH Zurich, which is partially funded by its consortium partners EAWAG, EMPA, and PSI. Nicolas Gruber acknowledges support from the European Union's Horizon 2020 research and innovation program under grant agreement no. 820989 (COMFORT).
We thank two anonymous reviewers for their comments that contributed to a substantial improvement of this paper.
This research has been supported by the Eidgenössische Technische Hochschule Zürich (grant no. ETH-06 16-1), the Swiss National Science Foundation (AMBIZIONE PZ00P2_154805), and the European Union's Horizon 2020 research and innovation program (grant agreement no. 820989, COMFORT).
This paper was edited by Arvind Singh and reviewed by two anonymous referees.
Alkire, M. B., Rember, R., and Polyakov, I.: Discrepancy in the Identification of the Atlantic/Pacific Front in the Central Arctic Ocean: NO Versus Nutrient Relationships, Geophys. Res. Lett., 46, 3843–3852, https://doi.org/10.1029/2018GL081837, 2019. a, b, c
Bauch, D., van der Loeff, M. R., Andersen, N., Torres-Valdes, S., Bakker, K., and Abrahamsen, E. P.: Origin of Freshwater and Polynya Water in the Arctic Ocean Halocline in Summer 2007, Prog. Oceanogr., 91, 482–495, https://doi.org/10.1016/j.pocean.2011.07.017, 2011. a
Carmack, E., Polyakov, I., Padman, L., Fer, I., Hunke, E., Hutchings, J., Jackson, J., Kelley, D., Kwok, R., Layton, C., Melling, H., Perovich, D., Persson, O., Ruddick, B., Timmermans, M.-L., Toole, J., Ross, T., Vavrus, S., and Winsor, P.: Toward Quantifying the Increasing Role of Oceanic Heat in Sea Ice Loss in the New Arctic, B. Am. Meteorol. Soc., 96, 2079–2105, https://doi.org/10.1175/BAMS-D-13-00177.1, 2015. a
Casacuberta, N., Masqué, P., Henderson, G., Rutgers van-der-Loeff, M., Bauch, D., Vockenhuber, C., Daraoui, A., Walther, C., Synal, H.-A., and Christl, M.: First 236U Data from the Arctic Ocean and Use of 236U/238U and 129I/236U as a New Dual Tracer, Earth Planet Sc. Lett., 440, 127–134, https://doi.org/10.1016/j.epsl.2016.02.020, 2016. a, b, c, d, e
Casacuberta, N., Christl, M., Vockenhuber, C., Wefing, A.-M., Wacker, L., Masqué, P., Synal, H.-A., and Rutgers van der Loeff, M.: Tracing the Three Atlantic Branches Entering the Arctic Ocean With 129I and 236U, J. Geophys. Res.-Oceans, 123, 6909–6921, https://doi.org/10.1029/2018JC014168, 2018. a, b, c, d, e, f, g, h
Castrillejo, M., Witbaard, R., Casacuberta, N., Richardson, C. A., Dekker, R., Synal, H.-A., and Christl, M.: Unravelling 5 Decades of Anthropogenic 236U Discharge from Nuclear Reprocessing Plants, Sci. Total Environ., 717, 137094, https://doi.org/10.1016/j.scitotenv.2020.137094, 2020. a
Christl, M., Casacuberta, N., Vockenhuber, C., Elsässer, C., Bailly du Bois, P., Herrmann, J., and Synal, H.-A.: Reconstruction of the 236U Input Function for the Northeast Atlantic Ocean: Implications for 129I/236U and 236U/238U-Based Tracer Ages, J. Geophys. Res.-Oceans, 120, 7282–7299, https://doi.org/10.1002/2015JC011116, 2015. a, b, c, d, e, f, g
Dmitrenko, I. A., Polyakov, I. V., Kirillov, S. A., Timokhov, L. A., Frolov, I. E., Sokolov, V. T., Simmons, H. L., Ivanov, V. V., and Walsh, D.: Toward a Warmer Arctic Ocean: Spreading of the Early 21st Century Atlantic Water Warm Anomaly along the Eurasian Basin Margins, J. Geophys. Res., 113, C05023, https://doi.org/10.1029/2007JC004158, 2008. a
Edmonds, H. N., Zhou, Z. Q., Raisbeck, G. M., Yiou, F., Kilius, L., and Edmond, J. M.: Distribution and Behavior of Anthropogenic 129I in Water Masses Ventilating the North Atlantic Ocean, J. Geophys. Res.-Oceans, 106, 6881–6894, https://doi.org/10.1029/1999JC000282, 2001. a
Eigl, R., Steier, P., Sakata, K., and Sakaguchi, A.: Vertical Distribution of 236U in the North Pacific Ocean, J. Environ. Radioactiv., 169–170, 70–78, https://doi.org/10.1016/j.jenvrad.2016.12.010, 2017. a
Gascard, J.-C., Raisbeck, G., Sequeira, S., Yiou, F., and Mork, K. A.: The Norwegian Atlantic Current in the Lofoten Basin Inferred from Hydrological and Tracer Data (129I) and Its Interaction with the Norwegian Coastal Current, Geophys. Res. Lett., 31, L01318, https://doi.org/10.1029/2003GL018303, 2004. a
Haine, T. W. N. and Hall, T. M.: A Generalized Transport Theory: Water-Mass Composition and Age, J. Phys. Oceanogr., 32, 1932–1946, https://doi.org/10.1175/1520-0485(2002)032<1932:AGTTWM>2.0.CO;2, 2002. a, b, c, d
Hall, T. M., Haine, T. W. N., and Waugh, D. W.: Inferring the Concentration of Anthropogenic Carbon in the Ocean from Tracers, Global Biogeochem. Cy., 16, 78-1–78-15, https://doi.org/10.1029/2001GB001835, 2002. a
He, Y.-C., Tjiputra, J., Langehaug, H. R., Jeansson, E., Gao, Y., Schwinger, J., and Olsen, A.: A Model-Based Evaluation of the Inverse Gaussian Transit-Time Distribution Method for Inferring Anthropogenic Carbon Storage in the Ocean, J. Geophys. Res.-Oceans, 123, 1777–1800, https://doi.org/10.1002/2017JC013504, 2018. a
Karcher, M., Beszczynska-Möller, A., Kauker, F., Gerdes, R., Heyen, S., Rudels, B., and Schauer, U.: Arctic Ocean Warming and Its Consequences for the Denmark Strait Overflow, J. Geophys. Res., 116, C02037, https://doi.org/10.1029/2010JC006265, 2011. a, b
Karcher, M., Smith, J. N., Kauker, F., Gerdes, R., and Smethie, W. M.: Recent Changes in Arctic Ocean Circulation Revealed by Iodine-129 Observations and Modeling, J. Geophys. Res.-Oceans, 117, C08007, https://doi.org/10.1029/2011JC007513, 2012. a, b, c, d, e, f
Khatiwala, S., Tanhua, T., Mikaloff Fletcher, S., Gerber, M., Doney, S. C., Graven, H. D., Gruber, N., McKinley, G. A., Murata, A., Ríos, A. F., and Sabine, C. L.: Global ocean storage of anthropogenic carbon, Biogeosciences, 10, 2169–2191, https://doi.org/10.5194/bg-10-2169-2013, 2013. a
Mauldin, A., Schlosser, P., Newton, R., Smethie, W. M., Bayer, R., Rhein, M., and Jones, E. P.: The Velocity and Mixing Time Scale of the Arctic Ocean Boundary Current Estimated with Transient Tracers, J. Geophys. Res., 115, C08002, https://doi.org/10.1029/2009JC005965, 2010. a, b, c, d
Meredith, M., Sommerkorn, M., Cassotta, S., Derksen, C., Ekaykin, A., Hollowed, A., Kofinas, G., Mackintosh, A., Melbourne-Thomas, J., Muelbert, M. M. C., Ottersen, G., Pritchard, H., and Schuur, E. A. G.: Polar Regions, in: 2019 IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, edited by: Pörtner, H.-O., Roberts, D. C., Masson-Delmotte, V., Zhai, P., Tignor, M., Poloczanska, E., Mintenbeck, K., Alegría, A., Nicolai, M., Okem, A., Petzold, J., Rama, B., and Weyer, N. M., in press, 2021. a
Morison, J., Kwok, R., Peralta-Ferriz, C., Alkire, M., Rigor, I., Andersen, R., and Steele, M.: Changing Arctic Ocean Freshwater Pathways, Nature, 481, 66–70, https://doi.org/10.1038/nature10705, 2012. a
Olsen, A., Omar, A. M., Jeansson, E., Anderson, L. G., and Bellerby, R. G. J.: Nordic Seas Transit Time Distributions and Anthropogenic CO2, J. Geophys. Res., 115, C05005, https://doi.org/10.1029/2009JC005488, 2010. a
Polyakov, I. V., Beszczynska, A., Carmack, E. C., Dmitrenko, I. A., Fahrbach, E., Frolov, I. E., Gerdes, R., Hansen, E., Holfort, J., Ivanov, V. V., Johnson, M. A., Karcher, M., Kauker, F., Morison, J., Orvik, K. A., Schauer, U., Simmons, H. L., Skagseth, O. y., Sokolov, V. T., Steele, M., Timokhov, L. A., Walsh, D., and Walsh, J. E.: One More Step toward a Warmer Arctic, Geophys. Res. Lett., 32, L17605, https://doi.org/10.1029/2005GL023740, 2005. a
Polyakov, I. V., Alexeev, V. A., Ashik, I. M., Bacon, S., Beszczynska-Möller, A., Carmack, E. C., Dmitrenko, I. A., Fortier, L., Gascard, J.-C., Hansen, E., Hölemann, J., Ivanov, V. V., Kikuchi, T., Kirillov, S., Lenn, Y.-D., McLaughlin, F. A., Piechura, J., Repina, I., Timokhov, L. A., Walczowski, W., and Woodgate, R.: Fate of Early 2000s Arctic Warm Water Pulse, B. Am. Meteorol. Soc., 92, 561–566, https://doi.org/10.1175/2010BAMS2921.1, 2011. a
Polyakov, I. V., Pnyushkov, A. V., Alkire, M. B., Ashik, I. M., Baumann, T. M., Carmack, E. C., Goszczko, I., Guthrie, J., Ivanov, V. V., Kanzow, T., Krishfield, R., Kwok, R., Sundfjord, A., Morison, J., Rember, R., and Yulin, A.: Greater Role for Atlantic Inflows on Sea-Ice Loss in the Eurasian Basin of the Arctic Ocean, Science, 356, 285–291, https://doi.org/10.1126/science.aai8204, 2017. a
Rajasakaren, B., Jeansson, E., Olsen, A., Tanhua, T., Johannessen, T., and Smethie, W. M.: Trends in Anthropogenic Carbon in the Arctic Ocean, Prog. Oceanogr., 178, 102177, https://doi.org/10.1016/j.pocean.2019.102177, 2019. a
Rudels, B.: Arctic Ocean Circulation, Processes and Water Masses: A Description of Observations and Ideas with Focus on the Period Prior to the International Polar Year 2007–2009, Prog. Oceanogr., 132, 22–67, https://doi.org/10.1016/j.pocean.2013.11.006, 2015. a, b, c, d, e, f, g
Rudels, B., Meyer, R., Fahrbach, E., Ivanov, V. V., Østerhus, S., Quadfasel, D., Schauer, U., Tverberg, V., and Woodgate, R. A.: Water mass distribution in Fram Strait and over the Yermak Plateau in summer 1997, Ann. Geophys., 18, 687–705, https://doi.org/10.1007/s00585-000-0687-5, 2000. a
Rudels, B., Björk, G., Nilsson, J., Winsor, P., Lake, I., and Nohr, C.: The Interaction between Waters from the Arctic Ocean and the Nordic Seas North of Fram Strait and along the East Greenland Current: Results from the Arctic Ocean-02 Oden Expedition, J. Marine Syst., 55, 1–30, https://doi.org/10.1016/j.jmarsys.2004.06.008, 2005. a
Rudels, B., Anderson, L., Eriksson, P., Fahrbach, E., Jakobsson, M., Jones, E. P., Melling, H., Prinsenberg, S., Schauer, U., and Yao, T.: Observations in the Ocean, in: Arctic Climate Change, vol. 43, edited by: Lemke, P. and Jacobi, H.-W., Springer Netherlands, Dordrecht, https://doi.org/10.1007/978-94-007-2027-5_4, 117–198, 2012. a
Sakaguchi, A., Kawai, K., Steier, P., Quinto, F., Mino, K., Tomita, J., Hoshi, M., Whitehead, N., and Yamamoto, M.: First Results on 236U Levels in Global Fallout, Sci. Total Environ., 407, 4238–4242, https://doi.org/10.1016/j.scitotenv.2009.01.058, 2009. a
Smith, J. N., Jones, E. P., Moran, S., Smethie, W. M., and Kieser, W.: Iodine 129/CFC 11 Transit Times for Denmark Strait Overflow Water in the Labrador and Irminger Seas, J. Geophys. Res., 110, C05006, https://doi.org/10.1029/2004JC002516, 2005. a
Smith, J. N., McLaughlin, F. A., Smethie, W. M., Moran, S. B., and Lepore, K.: Iodine-129, 137Cs, and CFC-11 Tracer Transit Time Distributions in the Arctic Ocean, J. Geophys. Res., 116, C04024, https://doi.org/10.1029/2010JC006471, 2011. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t
Stöven, T., Tanhua, T., Hoppema, M., and Bullister, J. L.: Perspectives of transient tracer applications and limiting cases, Ocean Sci., 11, 699–718, https://doi.org/10.5194/os-11-699-2015, 2015. a, b
Stöven, T., Tanhua, T., Hoppema, M., and von Appen, W.-J.: Transient tracer distributions in the Fram Strait in 2012 and inferred anthropogenic carbon content and transport, Ocean Sci., 12, 319–333, https://doi.org/10.5194/os-12-319-2016, 2016. a, b, c
Swift, J., Jones, E., Aagaard, K., Carmack, E., Hingston, M., MacDonald, R., McLaughlin, F., and Perkin, R.: Waters of the Makarov and Canada Basins, Deep-Sea Res. Pt. II, 44, 1503–1529, https://doi.org/10.1016/S0967-0645(97)00055-6, 1997. a
Tanhua, T., Jones, E. P., Jeansson, E., Jutterström, S., Smethie, W. M., Wallace, D. W. R., and Anderson, L. G.: Ventilation of the Arctic Ocean: Mean Ages and Inventories of Anthropogenic CO2 and CFC-11, J. Geophys. Res., 114, C01002, https://doi.org/10.1029/2008JC004868, 2009. a, b, c, d, e
Terhaar, J., Orr, J. C., Gehlen, M., Ethè, C., and Bopp, L.: Model constraints on the anthropogenic carbon budget of the Arctic Ocean, Biogeosciences, 16, 2343–2367, https://doi.org/10.5194/bg-16-2343-2019, 2019. a
Terhaar, J., Kwiatkowski, L., and Bopp, L.: Emergent Constraint on Arctic Ocean Acidification in the Twenty-First Century, Nature, 582, 379–383, https://doi.org/10.1038/s41586-020-2360-3, 2020a. a, b, c
Terhaar, J., Tanhua, T., Stöven, T., Orr, J. C., and Bopp, L.: Evaluation of Data-Based Estimates of Anthropogenic Carbon in the Arctic Ocean, J. Geophys. Res.-Oceans, 125, e2020JC016124, https://doi.org/10.1029/2020JC016124, 2020b. a
Ulfsbo, A., Jones, E. M., Casacuberta, N., Korhonen, M., Rabe, B., Karcher, M., and van Heuven, S. M.: Rapid Changes in Anthropogenic Carbon Storage and Ocean Acidification in the Intermediate Layers of the Eurasian Arctic Ocean: 1996–2015, Global Biogeochem. Cy., 32, 1254–1275, https://doi.org/10.1029/2017GB005738, 2018. a
Waugh, D. W., Haine, T. W., and Hall, T. M.: Transport Times and Anthropogenic Carbon in the Subpolar North Atlantic Ocean, Deep-Sea Res. Pt. I, 51, 1475–1491, https://doi.org/10.1016/j.dsr.2004.06.011, 2004. a
Waugh, D. W., Hall, T. M., Mcneil, B. I., Key, R., and Matear, R. J.: Anthropogenic CO2 in the Oceans Estimated Using Transit Time Distributions, Tellus B, 58, 376–389, https://doi.org/10.1111/j.1600-0889.2006.00222.x, 2006. a
Wefing, A.-M., Christl, M., Vockenhuber, C., Rutgers van der Loeff, M., and Casacuberta, N.: Tracing Atlantic Waters Using 129I and 236U in the Fram Strait in 2016, J. Geophys. Res.-Oceans, 124, 882–896, https://doi.org/10.1029/2018JC014399, 2019. a, b, c, d, e, f
Woodgate, R.: Arctic Ocean Circulation: Going Around At the Top Of the World, Nature Education Knowledge, 4, 8, available at: https://www.nature.com/scitable/knowledge/library/arctic-ocean-circulation-going-around-at-the-102811553/ (last access: 13 June 2018), 2013. a, b, c