the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
![](https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-avatar-thumb150.png)
Revising chronological uncertainties in marine archives using global anthropogenic signals: a case study on the oceanic 13C Suess effect
Ulysses S. Ninnemann
Are Olsen
Neil L. Rose
David J. R. Thornalley
Tor L. Mjell
François Counillon
Marine sediments are excellent archives for reconstructing past changes in climate and ocean circulation. Overlapping with instrumental records, they hold the potential to elucidate natural variability and contextualize current changes. Yet, dating uncertainties of traditional approaches (e.g., up to ± 30–50 years for the last 2 centuries) pose major challenges for integrating the shorter instrumental records with these extended marine archives. Hence, robust sediment chronologies are crucial, and most existing age model constraints do not provide sufficient age control, particularly for the 20th century, which is the most critical period for comparing proxy records to historical changes. Here we propose a novel chronostratigraphic approach that uses anthropogenic signals such as the oceanic 13C Suess effect and spheroidal carbonaceous fly-ash particles to reduce age model uncertainties in high-resolution marine archives. As a test, we apply this new approach to a marine sediment core located at the Gardar Drift, in the subpolar North Atlantic, and revise the previously published age model for this site. We further provide a refined estimate of regional reservoir corrections and uncertainties for Gardar Drift.
- Article
(4021 KB) - Full-text XML
-
Supplement
(568 KB) - BibTeX
- EndNote
Among the most prominent features of 20th century climate in the circum-North Atlantic are the observed basin-wide multi-decadal variations in the Atlantic Ocean sea surface temperatures (SSTs) – the Atlantic multi-decadal variability, AMV. This has impacts on the North American and European climate (Sutton and Hodson, 2005), the frequency of Atlantic hurricanes (Goldenberg et al., 2001), the extent of Arctic sea ice (Miles et al., 2014), and rainfall patterns in the African Sahel (Wang et al., 2012). However, instrumental SST records are limited to the last ∼ 150 years (e.g., Kaplan et al., 1998), and only in a few location has widespread coverage existed since the 1950s onwards. Yet longer records of climate and ocean circulation are required to understand and assess the mechanisms behind its variability. For example, it is still debated whether AMV is driven internally, linked to multi-decadal variations in the Atlantic Meridional Overturning Circulation (AMOC) (Zhang et al., 2019), or driven externally, e.g., due to solar and volcanic forcings (Otterå et al., 2010) or the timing of anthropogenic forcings (Booth et al., 2012); it is even debated whether or not such an oscillation exists at all (Mann et al., 2020). Annually laminated mollusk shell archives offer the excellent chronological constraint required to investigate such questions; however, they are limited to shelf locations, and the range of proxies that can be applied in these archives is limited (Reynolds et al., 2016). Also overlapping with and extending the instrumental records further back in time, marine sediments hold the potential to resolve these issues and contextualize current changes. New high-resolution proxy records, particularly from North Atlantic sedimentary drift sites, are now emerging, closing the time gap between modern and paleo-observations (e.g., Boessenkool et al., 2007; Mjell et al., 2016; Thornalley et al., 2018; Spooner et al., 2020). For instance, Mjell et al. (2016) found that AMV and deep-ocean circulation varied on similar timescales over the last 600 years; however, due to age model uncertainties as high as the duration of half an AMV cycle, determining the precise phasing was not possible and required independent age constraints. Hence, integrating nearly continuous but shorter observational records with longer (but with relatively lower resolution) marine archives still poses one of the major challenges for the (paleo)oceanographic community.
Recent marine sediments are dated using an array of approaches, all of which have their own limitations and uncertainties. Radiocarbon (14C) dating is one of the most common methods for dating marine sediment cores. The uncertainties with this method can exceed 50 years and include several caveats and assumptions such as uncertain and variable reservoir effects and confounding influences such as the effect of fossil fuel emissions on atmospheric radiocarbon and the H-bomb 14C spike, which further increases the uncertainties when dating recent sediments (Reimer et al., 2004; Hughen, 2007; Graven, 2015). In the latter case, the 14C bomb spike can serve as an additional high-resolution dating tool in marine settings, yet this requires annually resolved archives (Scourse et al., 2012). Geochemical composition of tephra shards and fingerprinting these to known volcanic eruptions can also provide absolute age markers. The precision of these age markers can be 1–2 years, yet this method is only regionally applicable and the occurrence of multiple, closely spaced eruptions with similar geochemistry can lead to greater uncertainty (Lowe, 2011). Radionuclide dating (210Pb, 137Cs, 241Am) (Appleby, 2008) and more recently the increases in mercury (Hg) concentrations (i.e., as an anthropogenic (pollution) indicator) are used as chronostratigraphic markers on recent marine sediments (Moros et al., 2017; Perner et al., 2019). For instance, 210Pb dating is widely used for dating recent sediments (0–150 years), while chronostratigraphic markers such as fallout from testing nuclear weapons in 1963 and Chernobyl fallout in 1986 can also be determined from the presence of 137Cs (Appleby, 2008). Still, 210Pb-based age models also involve multiple assumptions and are ideally validated using an independent age marker (e.g., 137Cs or 241Am) to assess the influence of post-depositional remobilization or bioturbation. Yet, it remains difficult to confirm to what extent the assumptions for dating are met (Smith, 2001). 137Cs profiles are often used to partially validate 210Pb chronologies, but this can only be undertaken for specific periods (e.g., bomb testing, Chernobyl). In addition, 137Cs is also prone to post-depositional remobilization and is not always above the detection limit – depending on core locations (e.g., Barsanti et al., 2020). Although the application of 210Pb dating in combination with 137Cs in lacustrine environments is well established, delayed input from 137Cs fallouts highlights the need for care in using 137Cs as chronostratigraphic markers even in lake sediments (Appleby et al., 2023). The situation is considerably more difficult in marine environments (e.g., Appleby et al., 2021). Indeed, a recent review highlights the continuing importance of, and need for, independent age control markers to corroborate 210Pb-based age models (Barsanti et al., 2020). Clearly progress is needed to improve age constraints in the 20th century in a way that will allow us to calibrate proxies using observational time series and, ultimately, reliably extend these observational records. Anthropogenic signals, such as the oceanic 13C Suess effect and spheroidal carbonaceous fly-ash particles (SCPs), are evident in high-resolution marine archives and hold the potential to provide a means for improving age control over the 20th century.
Atmospheric CO2 has been increasing due to human activities, such as fossil fuel combustion and deforestation, since the beginning of the industrial period. Due to preferential uptake of the lighter isotope (i.e., 12C), increased anthropogenic CO2 emissions cause the ratio (δ13C) and the 14C/C ratio (Δ14C) to decline. The decreasing trend in the radiocarbon () content of CO2 was first named the “Suess effect” by Suess (1955). In 1979, due to its similarity, Keeling (1979) extended the Suess effect terminology to the shifts in the ratio of the atmospheric CO2. The 13C Suess effect propagates into different reservoirs of the Earth system; for instance, the addition of low-δ13C anthropogenic CO2 from the atmosphere to the surface ocean also affects the natural δ13C gradients (Eide et al., 2017; Olsen and Ninnemann, 2010). Foraminiferal δ13C records (planktonic and benthic) from high-resolution marine archives capture this accelerating decline in δ13C over the last century (e.g., Mellon et al., 2019) and thus hold huge potential for refining age control for recent sediments.
Another new and promising approach for dating recent marine sediments is the use of spheroidal carbonaceous fly-ash particles (SCPs) (Spooner et al., 2020; Thornalley et al., 2018). SCPs are only produced from high-temperature industrial sources, such as coal and oil, and are thus purely anthropogenic in origin. They are emitted to the atmosphere along with combustion flue gases and are therefore transported to and recorded in many natural archives worldwide – including regions that are remote from industrial sources (e.g., Rose et al., 2004, 2012). In lake sediment records, SCPs were first observed during the mid-19th century in the UK, Europe, and North America and show a very distinct concentration profile. The SCP concentration trend starts with a gradual increase from the beginning of the SCP record until the mid-20th century, followed by a rapid increase at ca. 1950 linked to the increased demand for electricity following the Second World War (Rose, 2015). The beginning of the SCP record may vary regionally because it depends on the regional developments in industrial history as well as the sedimentation rates. However, the rapid increase observed in the mid-20th century has been considered to be a global signal (Rose, 2015) – making SCPs a robust and ideal stratigraphic marker for a mid-20th century Anthropocene. The first applications of the SCP method to marine sediment archives (Thornalley et al., 2018; Spooner et al., 2020; Kaiser et al., 2023) have been shown to follow similar trends to those established from lake records (Rose, 2015), providing an independent means to improve marine-based chronologies over the last 150 years.
Here we combine these two novel chronostratigraphic approaches that use anthropogenic signals (i.e., oceanic 13C Suess effect change and spheroidal carbonaceous fly-ash particles – SCPs) to reduce age model uncertainties in high-resolution marine archives. As a test, we apply this new approach to a high-resolution site at the Gardar Drift, off southern Iceland, to revise the previously published age model at this site (i.e., Mjell et al., 2016). We further provide refined regional 14C reservoir corrections and uncertainties for Gardar Drift using a combination of accelerator mass spectrometry (AMS) radiocarbon dates and oceanic 13C Suess effect estimates for our core location.
In this study we use sediment samples from the Gardar Drift multicore, GS06-144-09 MC (60°19′ N, 23°58′ W; 2081 m water depth), recovered during the University of Bergen cruise no. GS06-144 on board the research vessel (R/V) G. O. Sars. Four successful identical cores (GS06-144-09 MC A-D) were recovered at this station. The 44.5 cm long GS06-144-09 MC-D has been sampled at 0.5 cm intervals. Each sample was soaked in distilled water and shaken for 12 h in order to disperse the sediment before they were wet-sieved and separated into size fractions of > 63 and < 63 µm. The fine fractions (< 63 µm) were used for mean sortable silt grain size analysis (Mjell et al., 2016), whereas the > 63 µm fraction was used for selection of foraminifera for stable isotope analysis and 14C AMS dating (Table 1). The 44 cm long GS06-144-09 MC-C was sampled at 0.5 cm intervals. Each sample was dried and weighed. Dry bulk sediment samples from GS06-144-09 MC-C were used for SCP analysis.
Samples from GS06-144-09 MC-D have previously been analyzed for the activity of 210Pb, 226Ra, and 137Cs at the Gamma Dating Centre, Department of Geosciences and Natural Resource Management, University of Copenhagen, Denmark (Mjell et al., 2016). The initial age model of GS06-144-09 MC-D was based on 210Pb excess dates from the top 7.25 cm and two 14C AMS dates (Mjell et al., 2016). The presence of 137Cs in marine sediment cores is often used to validate the 210Pb chronologies and can also provide additional information (e.g., an independent tie point) for the onset of atmospheric weapon testing (e.g., Perner et al., 2018). In core GS06-144-09 MC-D, the content of 137Cs was very low and below the detection limit except in the top 4 cm of the core. This may indicate that the top 4 cm could be younger than ∼ 1950 CE. However, for core GS06-144-09 MC-D, traces (near detection limit) of 137Cs were also episodically present below this depth (Fig. S1 in Mjell et al., 2016). Hence, here we choose not to include the information provided by 137Cs in our age model, and we also do not include the 210Pb dates, as it will not be possible to validate with 137Cs.
In general, an ideal approach to build the best possible chronology is to integrate all available information. However, here we aim to demonstrate the potential utility of two novel approaches, oceanic 13C Suess effect change and SCPs, in building robust marine sediment chronologies. Therefore, we focus on these two novel techniques in a more stand-alone manner to assess their utility independently and their consistency with each other. Our methods include stable carbon isotopes of planktonic foraminifera (δ13C), 14C AMS dates, SCP analysis, and time series of oceanic 13C Suess effect change computed for our core location.
2.1 Stable isotope analysis (δ13C)
Stable isotope analyses (δ13C) were performed on the planktonic foraminifera Globigerina bulloides, Neogloboquadrina incompta, and Globorotalia inflata every 0.5 cm throughout the core. G. bulloides was picked from the 250–300 µm size fraction, while N. incompta was picked from 150–250 µm and G. inflata was picked from the 250–350 µm size fraction. Approximately 5–7 shells of G. bulloides, ∼ 5 shells of G. inflata, and ∼ 10 shells of N. incompta from each sample were used for stable isotope analysis. Foraminifera were ultrasonically rinsed for 20 s in methanol to remove any contaminants prior to analysis. Stable isotope analyses were measured using a Finnigan MAT 251 and a MAT 253 mass spectrometer in the FARLAB (Facility for Advanced Isotopic Research) at the Department of Earth Science, University of Bergen. All samples were run in two replicates whenever foraminifera were sufficiently abundant. The stable isotope results are expressed as the average of the two replicate measurements and reported relative to Vienna Pee Dee Belemnite (VPDB), calibrated using NBS-19. Long-term analytical precision (1σ) of the standards over the analysis period was better than 0.04 ‰ for δ13C.
2.2 13C Suess effect estimates
Recently, Eide et al. (2017) calculated globally gridded surface-to-seabed 13C Suess effect estimates for the industrialized era. These estimates were based on the two-step back-calculation technique of Olsen and Ninnemann (2010) for waters deeper than 200 m, while for waters above they were determined by combining the 200 m level estimate with values of the surface ocean 13C Suess effect as evident in coral and sclerosponge records. The two-step back-calculation approach first takes advantage of the relationships between preformed δ13C and chlorofluorocarbons (CFC-11 or CFC-12) in the ocean to quantify the 13C Suess effect since CFCs first appeared in the atmosphere (the 1940s). In the second step, these estimates are extended to the full industrialized era under the assumption of transient steady state (Gammon et al., 1982; Tanhua et al., 2007), which states that after an initial adjustment period, the response in tracer concentrations at depth will be proportional to the change in boundary concentration in exponentially forced systems. This means that we can expect that the ratio of the 13C Suess effect at any point in the ocean to that in the atmosphere will remain constant in time, i.e.,
where Δt1 and Δt2 represent two time intervals since the preindustrial. In the case of Eide et al. (2017) these are the periods 1940 to 1994 and preindustrial (defined as atmospheric δ13C = −6.5) to 1994.
![https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f01](https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f01-thumb.png)
Figure 1Planktonic δ13C records from Site GS06-144-09 MC-D plotted vs. depth (cm). Yellow highlighting marks the sharp decline in δ13C due to the Suess effect. (a) G. bulloides δ13C record (blue) with five-point mean (bold line), (b) N. incompta δ13C record (green) with five-point mean (bold line), and (c) G. inflata δ13C record (pink) with five-point mean (bold line). The five-point mean is extended into the core top by taking the mean of samples at 0 and 0.5 cm, shown as dashed bold lines to highlight the large abrupt δ13C decrease at the core top.
Here, we use Eq. (1) to derive time series of the Suess effect since the preindustrial at 10 depth layers from the surface to 200 m (e.g, δ13CSE_0, δ13CSE_50) above the Gardar Drift core site. This depth interval covers the depth habitats of the planktonic foraminiferal species we have used for stable isotope analysis. The time series were determined by taking the ratio between the Suess effect determined by Eide et al. (2017) at each of the 10 depth levels we consider in the grid box covering the Gardar Drift (60–61° N, 23–24° W) and the atmospheric δ13C decline until 1994 and multiplying this by the atmospheric δ13C history since the preindustrial provided by Rubino et al. (2013). The thus calculated marine Suess effect time series are presented in Fig. 2. We set the starting point in time to 1800, as an appreciable decline in atmospheric δ13C is only visible after that year.
2.3 SCP analysis
We followed the SCP method outlined by Rose (1994). Approximately 0.2 g of dried bulk sediment was weighed into 15 mL polypropylene tubes. One SCP reference standard (Rose, 2008) and a blank were included for quality control purposes and treated exactly the same as the samples. The SCP extraction method included nitric acid (HNO3), hydrofluoric acid (HF), and hydrochloric acid (HCl) stages to respectively remove organic matter, silicious material, and carbonates. Following the acid digestion stages, a known fraction of the final residue was evaporated onto a cover slip and mounted on a microscope slide using Naphrax mountant. A light microscope with 400× magnification was used to identify and count the total number of SCPs on each slide. SCP identification followed the criteria described in Rose (2008) based on morphology, color, depth, and porosity. SCP concentrations are reported as the number of SCPs per gram of dry sediment (gDM−1). SCP analyses were performed at the Department of Geography, University College London. The concentration of the SCP reference material was 5318 gDM−1 (±1022, 90 % confidence level), close to the reported concentration of 6005 ± 70 gDM−1 (Rose, 2008). No SCPs were observed in the blank.
3.1 Planktonic foraminiferal δ13C vs. the oceanic 13C Suess effect
In the subpolar North Atlantic, G. bulloides calcifies in the upper 50 m of the water column over the late spring and summer, depending on food availability (Jonkers et al., 2013; Schiebel et al., 1997; Spero and Lea, 1996; Chapman, 2010). On the other hand, the habitat depth of N. incompta is highly variable, ranging from the surface to deeper thermocline, most likely calcifying between 50 and 125 m water depth (Chapman, 2010; Field, 2004; Pak and Kennett, 2002; Pak et al., 2004; Von Langen et al., 2005; Nyland et al., 2006; Schiebel et al., 2001). G. inflata is a deep-dwelling foraminiferal species, living at the base of the seasonal thermocline or deeper in the main thermocline if the base of the seasonal thermocline is warmer than 16 °C (Cléroux et al., 2007). In the North Atlantic, G. inflata calcifies between 200 and 400 m south of 57° N and between 100 and 200 m north of 57° N (Ganssen and Kroon, 2000).
To calculate the age estimates based on the 13C Suess effect, we assume a calcification depth of 50 m for G. bulloides and compare our G. bulloides δ13C record with the 13C Suess effect change at 50 m (δ13CSE_50) at our core location. In order to avoid any uncertainties regarding planktonic foraminiferal depth habitats, we also present a stacked planktonic δ13C record (δ13Cstack – i.e., the average of G. bulloides, N. incompta, and G. inflata) and compare it with the average 13C Suess effect change over the top 200 m of the water column (δ13CSE_0–200), which spans the depth habitats of all three planktonic species (i.e., G. bulloides, N. incompta, and G. inflata) used in this study (Figs. S1 and S2 in the Supplement).
3.1.1 G. bulloides δ13C vs. oceanic Suess effect change at 50 m (δ13CSE_50)
The G. bulloides δ13C record shows large natural variability over the 10–44 cm core interval, varying between ∼ 0.08 ‰ and ∼ −0.6 ‰. However, the most prominent feature occurs towards the core top. δ13C values reach a peak of 0.27 ‰ at 7.5 cm, then start to gradually decrease and reach 0.05 ‰ at 1 cm. This is followed by a very sharp decline of ∼ 0.8 ‰ centered at 0.5 cm. The gradual decrease observed in G. bulloides δ13C, with a sharper decline at the core top, indicates the presence of the 13C Suess effect. Compared to the 13C Suess effect change at 50 m, the relative change in G. bulloides δ13C seems to be very similar (Fig. S3 in the Supplement). Does the δ13CSE_50 curve provide a means to narrow down chronological uncertainties over the industrial period? To explore this, we objectively matched our G. bulloides δ13C record with the δ13CSE_50 curve to find the starting point (1800 CE) of the Suess effect curve on the G. bulloides δ13C record.
To objectively place the start of the δ13CSE_50 curve (1800 CE) on the G. bulloides δ13C record, first we computed the curvature of the δ13CSE_50 curve. We use a third-degree polynomial fit using the polyfit function in MATLAB. Secondly, we apply third-degree polynomial curve fits to the G. bulloides δ13C record for different core depth intervals (n=12) starting from 12 cm to cover the whole industrial period. We apply curve fits to 12–0, 11–0, 10–0, 9–0, 8.5–0, 8–0, 7.5–0, 7–0, 6.5–0 6–0, 5.5–0, and 5–0 cm intervals. When applying curve fits, we use G. bulloides δ13C, as well as its three-point running mean and five-point running mean, assuming the overall trends might be better represented in the smoothed data. Goodness-of-fit results for each curve fit are presented in Table S1 in the Supplement. Finally, we compared the curvature of the δ13CSE_50 curve with the various curve fits applied to G. bulloides δ13C records to find which curve fit is the most similar to the curvature of δ13CSE_50. To do this, we calculate the correlation coefficients between our target curve (in this case, the curvature of δ13CSE_50) and each of the third-degree polynomial curves using their individual polynomial coefficients (i.e., p1, p2, p3, p4; Table S1). The curvature of the G. bulloides δ13C record for the 7.5–0 cm interval is the most similar to the curvature of our δ13CSE_50 curve (r = 0.73), suggesting 7.5 cm could be 1800 CE. We further do the same test using the three-point and five-point running mean of the data. Although the correlation is poorer, the same result is also reached (i.e., best fit when the record starts at 7.5 cm) when three-point running mean (r = 0.46) and five-point running mean (r = 0.22) of G. bulloides δ13C are used. Placing the start of the oceanic Suess effect change (∼ 1800 CE) on our G. bulloides δ13C record is one of the main challenges of our approach as the most prominent δ13C decline does not happen until the most recent years or core top. This is also evident from our correlation analysis results (Table S1). For instance, a close second-best fit occurs when we place 1800 CE at 5 cm (r = 0.69) instead of 7.5 cm (r = 0.73). A comparison between the two correlation coefficients using the Fisher's z transformation suggests that the difference between the correlation coefficients is not statistically significant (z = 0.064, p = 0.949). This indicates that 5 cm could also be 1800 CE.
3.1.2 δ13Cstack vs. the 0–200 m average oceanic Suess effect change (δ13CSE_0–200)
N. incompta and G. inflata δ13C also followed a very similar trend as the G. bulloides δ13C record – with the most prominent decline towards the core top, indicating the presence of the 13C Suess effect. To cross-check our approach described in Sect. 3.1.1 and to avoid any uncertainties that may be caused due to habitat depth variability, we use the stacked planktonic δ13C of G. bulloides, N. incompta, and G. inflata (δ13Cstack). Considering the habitat depth range of all three planktonic species, we then compare the δ13Cstack with the 0–200 m average of the 13C Suess effect (δ13CSE_0–200) (Figs. S1 and S2). Similarly, to place the start of the δ13CSE_0–200 curve (1800 CE) on our δ13Cstack record, first we find the curvature of the δ13CSE_0–200 curve. We use a third-degree polynomial fit using the polyfit function in MATLAB. Secondly, we apply third-degree polynomial curve fits to the δ13Cstack record for the same core depth intervals as in Sect. 3.1.2. Finally, we compare the curvature of the Δ13CSE_0–200 curve with the various curve fits applied to our δ13Cstack record and find which curve fit is the most similar to the curvature of Δ13CSE_0–200. To do this, we again calculate the correlation coefficients between our target curve (in this case, the curvature of δ13CSE_0–200) and each of the third-degree polynomial curves using their individual polynomial coefficients (i.e., p1, p2, p3, p4; Table S1). In this case, we get similar results for intervals 5–0, 5.5–0, and 7.5–0 cm (r = −0.60). Although the negative correlation coefficients indicate that the similarity approach used here may not capture the complexity of comparing third-degree polynomials, it gives us a rough estimate of which curve is most similar to our target curve (i.e., Δ13CSE_0–200) and overall agrees with our initial finding based on G. bulloides that 7.5 or 5 cm may in fact be 1800 CE.
![https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f03](https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f03-thumb.png)
Figure 3Overview of core-top age calculation. (a) G. bulloides δ13C record (blue, with five-point mean; bold line) vs. depth (cm). The five-point mean is extended into the core top by taking the mean of samples at 0 and 0.5 cm, shown as dashed bold lines to highlight the large abrupt δ13C decrease at the core top. The dark gray line and gray shading respectively mark the mean and standard deviation of the 1–7.5 cm and 0 and 0.5 cm intervals. (b) δ13CSE_50 curve (pink). The arrow and dashed lines mark when a 0.57 ‰ magnitude decline occurs in the record.
3.2 Core-top age
In paleoceanographic studies it is common to use the year a sediment core was retrieved as the core-top age. However, this is highly dependent on the sedimentation rates of the region and may not always be the case. The core-top (0 cm) 14C AMS date for GS06-144-09 MC-D indicated the presence of bomb carbon, confirming that the top should be younger than ∼ 1957 CE (Mjell et al., 2016). Therefore, based on high sedimentation rates at the site, Mjell et al. (2016) assumed 2006 CE to be the core-top age, i.e., the year core GS06-144-09 MC was retrieved. Here we explore this further considering the new information provided by the relative change in our oceanic 13C Suess effect curve. For this, we use the G. bulloides δ13C record and the δ13CSE_50 curve.
Based on our previous curve fits, we place 1800 CE at 7.5 cm. The most prominent change in the G. bulloides δ13C record occurs at 0.5 cm. Hence, first we find the mean and standard deviation of the 1–7.5 cm interval (−0.05 ± 0.2; n=14), i.e., the mean δ13C over the industrial period and secondly the core top (two data points at 0 and 0.5 cm; −0.62 ± 0.17; n=2), of the G. bulloides δ13C record, i.e., where the sharpest decline due to the Suess effect occurs. We then calculated the difference in means using a t test (0.57 ‰) and found the magnitude of the sharpest decline in G. bulloides δ13C due to the Suess effect. Finally, we used the δ13CSE_50 curve to find when a 0.57 ‰ magnitude decline relative to the preindustrial value occurred. Based on our δ13CSE_50 curve, a decline of 0.57 ‰ occurs in ∼ 1972. This would then place 1972 CE at 0.25 cm (i.e., the mid-point of our two samples at 0 and 0.5 cm), suggesting a much older core-top age than previously assumed for GS06-144-09 MC (Mjell et al., 2016).
We repeated the same approach and evaluated how the core-top age would change if we placed 1800 CE at 5 cm (i.e., our second-best fit). We computed the mean and standard deviation of the 1–5 cm interval (−0.09 ± 0.2; n=9) and calculated the difference in means (0.53 ‰) between the 1–5 cm interval and the core top (0–0.5 cm). Finally, we determined when a magnitude of 0.53 ‰ decline occurred in the δ13CSE_50 curve. Based on our δ13CSE_50 curve, a decline of 0.53 ‰ occurred in ∼ 1969, placing 1969 CE at 0.25 cm. This suggests that placing 1800 CE at 7.5 cm vs. 5 cm changes our core-top age (or our tie point at 0.25 cm) by 3 years. When building our age model, here we choose 7.5 cm as 1800 CE (i.e., based on the best curve fit) and 0.25 cm as 1972 CE, and we introduce 3-year uncertainty to the selection of these tie points.
3.3 Revising regional reservoir corrections (ΔR) at Gardar Drift
To build an age model for the marine sediment cores based on radiocarbon dating it is necessary to convert 14C dates into calendar years. Surface ocean 14C is depleted relative to the atmosphere, which is known as the marine reservoir effect. Global marine radiocarbon calibration curves, e.g., the latest Marine20 curve (Heaton et al., 2020), account for the global average offset between the marine and atmospheric reservoirs; however, there are temporal and spatial deviations from this offset. Marine reservoir ages range from 400 years in the subtropics to more than 1000 years in polar oceans (Key et al., 2004). Therefore, the accurate calibration of 14C ages depends on knowledge of the local radiocarbon reservoir age of the surface ocean, i.e., the regional difference (ΔR) from the global marine radiocarbon calibration curve. The marine reservoir database within CALIB (http://calib.org/marine/, last access: 1 November 2023) is the most extensive and valuable source for ΔR values for the modern ocean (Reimer and Reimer, 2001; Stuiver and Reimer, 1986). This online platform provides the user with an average ΔR value for their core location based on the information provided on coordinates and number of nearest points. The ΔR values within the marine reservoir database are determined based on the known-age approach, i.e., when the death date (in calendar ages) of a pre-bomb marine sample (e.g., a mollusk shell) is known. However, as a consequence of nuclear tests in the 1950s and early 1960s the ΔR calculation with the known-age approach can only be applied to samples collected before 1950 CE; hence, the majority of the samples within the marine reservoir database are not homogenously distributed – making them temporally and spatially limited (Alves et al., 2018). Therefore, deriving a ΔR using the nearest points to a core location is problematic for many regions, where the closest ΔR is either not available or located in a different oceanographic setting (e.g., Hinojosa et al., 2015). When selecting samples for ΔR calculation, it is also important to review the ecological information on the taxa from which the ΔR value is derived, as some studies find species-specific values due to habitat, feeding mechanisms, and food sources. For instance, suspension feeders are thought to be the most suitable for dating, whereas deposit feeders, omnivore species, or carnivorous species are generally excluded due to their greater uncertainty in 14C ages as they incorporate old carbon (Pieńkowski et al., 2021; England et al., 2013; Forman and Polyak, 1997). However, some studies find no difference in 14C ages due to feeding mechanisms when the mollusks are derived from areas with no carboniferous rocks or local freshwater inputs to the surface ocean (Ascough et al., 2005).
Table S2 in the Supplement shows the ΔR values for our core site (GS06-144-09 MC; 60°19′ N, 23°58′ W), located south of Iceland, derived from the nearest points available in the marine reservoir database (Reimer and Reimer, 2001). When the 10 nearest points are used (i.e., based on the distance (km) from core location), the ΔR for our core site is −72 ± 64 14C yr. However, when we exclude carnivore and deposit feeding species, the ΔR value becomes −80 ± 54 14C yr. It is also important to note that even the individual samples have a large range of ΔR values, varying between −23 ± 45 and −220 ± 85 14C yr, suggesting there might be other factors influencing the ΔR. For instance, considering the oceanographic setting, another approach could be to only select samples located around southern Iceland – i.e., those potentially under the influence of the Irminger Current, where our core site lies. Then, the ΔR value would be −92 ± 93 14C yr (or −126 ± 66 14C yr when carnivore and deposit feeding species are excluded). This suggests that the available ΔR values within the CALIB marine reservoir database (Reimer and Reimer, 2001) for the region are highly variable and highly dependent on the selection criteria used by the investigator.
Global Ocean Data Analysis Project (GLODAP) radiocarbon observations (Key et al., 2004) provide an alternative approach to estimate the spatial variations in the reservoir ages (Gebbie and Huybers, 2012; Waelbroeck et al., 2019). For instance, Waelbroeck et al. (2019) extracted the pre-bomb surface mean (upper 250 m) reservoir ages from re-gridded (4° × 4°) GLODAP data. Following the Waelbroeck et al. (2019) approach, we extract the reservoir ages (443 ± 75.8 14C yr) at our core site (60° N, 24° W) from GLODAP data. Waelbroeck et al. (2019) note, however, that the error for their reservoir ages should be at least 100 14C yr if the computed GLODAP standard deviation is less than this value (i.e., in our case 443 ± 100 14C yr). Considering the global average marine reservoir age of ∼ 600 years based on Marine20 (Heaton et al., 2020), this would suggest a ΔR of −157 ± 100 14C yr for our region. The large difference (and/or uncertainties) in regional reservoir corrections extracted using two independent methods (e.g., CALIB vs. GLODAP-based) highlights the need for additional approaches to further constrain regional reservoir ages.
![https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f04](https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f04-thumb.png)
Figure 4G. bulloides δ13C (7.5–0.25 cm, light blue line plotted with the five-point mean; bold dark blue line) vs. the δ13CSE_50 curve spanning the 1800–1972 CE interval (bold pink line). Gray triangles on the depth axis mark the three 14C AMS samples at 2.5, 4, and 5 cm depth intervals, while the dashed gray lines and triangles on the age axis mark their corresponding “known ages” based on the Δ13CSE_50 comparison.
Here we suggest an alternative approach for calculating the ΔR for marine sediment cores that is independent of uncertainties such as the distance between core sites and sample locations, different oceanographic settings (e.g., coastal and fjord regions vs. open ocean), or the feeding ecology of the species used for dating. Based on our comparison of the G. bulloides δ13C record and the δ13CSE_50 curve we obtain two tie points, placing 1972 CE at 0.25 cm and 1800 CE at 7.5 cm. Figure 4 shows the G. bulloides δ13C record on a depth scale spanning the 7.5–0.25 cm core interval, plotted together with the δ13CSE_50 curve spanning the 1800–1972 CE interval. First, we estimate the “known ages” for depths 2.5, 4, and 5.5 cm (i.e., where we have 14C dates) by reading the corresponding ages from the δ13CSE_50 curve. Next, we calculate the ΔR value for each sample using the known-age approach in the online application deltar (Reimer and Reimer, 2017) based on the most recent Marine20 curve (Heaton et al., 2020). Finally, we calculate the weighted mean (Eq. 2) and standard deviation of ΔR following Reimer and Reimer (2001) and provide a revised ΔR estimate for the Gardar Drift. The uncertainty of ΔR is determined as the maximum value of either the weighted uncertainty in the mean of ΔR or the standard deviation of ΔR, as in Eqs. (3) and (5). Our refined ΔR estimate (−69 ± 38 14C yr) is similar to the value obtained from the marine reservoir database when the 10 nearest points are used (−72 ± 64 14C yr) – although with better uncertainty estimates.
3.4 Revised age model for GS06-144-09 MC
We use Bacon (version 2.5.0), the age–depth modeling approach that uses Bayesian statistics (Blaauw and Christen, 2011), operated through R (version 4.0.3) – free software for statistical computing and graphics. A total of 10 14C AMS dates (Table 1) are calibrated through Bacon using the most recent Marine20 curve (Heaton et al., 2020) and a ΔR value of −69 ± 38 (this study) – assuming a constant ΔR value throughout the core. Since our ΔR estimate is based on the comparison with the 13C Suess effect curve, we can only calculate a ΔR value for the last ∼ 200 years with this approach. Although we assume relatively stable conditions over the last millennium (e.g., compared to glacial–interglacial changes), changes in ocean circulation and ventilation before this period will also effect the ΔR in the region (e.g., during the Little Ice Age; Spooner et al., 2020).
![https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f05](https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f05-thumb.png)
Figure 5Age–depth plots of GS06-144-09 MC (a) when additional tie points for 0.25 cm (1972 CE) and 7.5 cm (1800 CE) are used and (b) when “known” calendar ages for samples at 2.5, 4, and 5.5 cm that were derived from the δ13CSE_50 comparison are used as additional tie points.
Additional tie points for 0.25 cm (1972 CE) and 7.5 cm (1800 CE) are used based on the information obtained from the Suess effect curve. Based on the core-top (0 cm) 14C AMS date (> 1950 CE) and the year the core was retrieved (2006 CE) the core-top age should be between ∼ 1950 and 2006 CE. As the core-top age cannot be younger than 2006 CE, we use this information as a prior in Bacon to set a minimum age limit for the core top. According to the revised age model, the date for the core top (0–0.5 cm) is 1977 CE. The average uncertainty for the last ∼ 200 years (i.e., the 0–7.5 cm interval) is ∼ ± 42 years and for the whole core (i.e., the 0–44 cm interval) is ± 90 years. The resulting age–depth plot is provided in Fig. 5a. Although Bacon selects the best age–depth model (i.e., dotted red lines in Fig. 5), considering the sedimentation rate profile based on the prior information, the tie points at the core top and 1800 CE play a crucial role, providing a basis for sedimentation rates. This is also seen from Fig. 5a, illustrated by the large range of 14C AMS dates that exceeds the calibration range of Marine20 due to bomb carbon. This further underscores the need for independent chronological approaches, particularly for the last century.
As a comparison, we also include the “known” calendar ages for samples at 2.5, 4, and 5.5 cm that were derived from the δ13CSE_50 comparison, together with their uncalibrated 14C dates, in the Bacon input file. For all the tie points derived from the δ13CSE_50 comparisons we add a ± 3-year uncertainty. Including the known calendar ages does not change the overall age model, but as expected, it highly decreases the age model uncertainties for the last ∼ 200 years (Fig. 5b). Based on this, the core-top age (0–0.5 cm) is again 1977 CE. The average age model uncertainty for the last ∼ 200 years (i.e., the 0–7.5 cm interval) is ± 17.5 years. Below this point, the uncertainty increases (average of ± 84 years for the 0–44 cm interval) and is highly dependent on the uncertainty of the 14C AMS dates. The average sedimentation rates for the top 0–7.5 cm interval are 43 and 63 cm kyr−1 for the 7.5–44 cm interval. The average sedimentation rate of the core (0–44 cm) is 59 cm kyr−1, giving a sample spacing of ∼ 8.5 years per 0.5 cm sample.
3.5 SCP analysis
To cross-check the validity of our Suess-effect-derived age model, here we use another independent approach: spheroidal carbonaceous fly-ash particles (SCPs).
SCP concentrations at GS06-144-09 MC are generally very low, varying between 152 and 616 gDM−1. Based on our revised age model, SCP concentrations start to gradually increase during the 1930s. A more marked increase in SCP concentrations occurs after 1954 and reaches peak values in 1966, followed by a decline towards the core top. Figure S4 in the Supplement shows a comparison of the GS06-144-09 MC SCP concentrations with previously published SCP profiles from Apavatn Lake, Iceland (Rose, 2015), and Nunatak Lake, Greenland (Bindler et al., 2001; Rose, 2015). Despite similarly low concentrations, both lake records show the same increase after ca. 1950 as the Gardar Drift marine sediment core. This suggests that the SCP concentrations at Gardar Drift follow a similar temporal pattern to the lake sediments in the region. Although the low SCP concentrations at GS06-144-09 MC result in considerable uncertainty for the SCP profile, the rapid increase after the 1950s at GS06-144-09 MC is consistent with the SCP trend in the region and is consistent with our Suess-effect-based revised age model.
![https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f06](https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f06-thumb.png)
Figure 6SCP concentration profile of GS06-144-09 MC plotted vs. the revised age model (as shown in Fig. 5b). The dashed red line marks 1950.
![https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f07](https://gchron.copernicus.org/articles/6/449/2024/gchron-6-449-2024-f07-thumb.png)
Figure 7Age–depth plot for the top 10 cm of GS06-144-09 MC to highlight the differences (e.g., in sedimentation rate) between the original 210Pb-based chronology (Mjell et al., 2016) and the tie points derived based on anthropogenic signals (this study). 14C dates are calibrated with CALIB (version 8.2) (Stuiver and Reimer, 1993) using Marine20 and ΔR = −69 ± 38 14C yr.
Our case study off Gardar Drift demonstrates the utility of two novel chronostratigraphic approaches that use anthropogenic signals (i.e., the oceanic 13C Suess effect change and SCP concentrations) in reducing age model uncertainties of recent high-resolution marine archives. In addition, using a combination of 14C AMS dates and oceanic 13C Suess effect estimates, we further provide refined regional 14C reservoir corrections and uncertainties for Gardar Drift. Despite the similarity of our refined ΔR estimate to those available in the marine reservoir database (Reimer and Reimer, 2001), it is also important to note the shortcomings of our approach. For instance, by reading the corresponding known ages from the 13C Suess effect curve to calculate ΔR, our approach assumes constant sedimentation rates and no bioturbation or reworking at the core top. Although we do not see any visible traces of bioturbation in our core, we acknowledge that this is often not the case, and we rarely have sites with true known ages. One exception to this, with potential to overcome this limitation, would be to use absolute age markers derived by identifying tephra layers and fingerprinting these to known volcanic eruptions. Yet this method is also only applicable in specific geologic settings and can also be affected by bioturbation – a limitation shared by all dating methods.
Bioturbation is one of the main sources of uncertainties of our approach as it will typically influence the age distributions and smooth the record. Generally, the smoothing, or attenuation, is greater when the sedimentation rates are low (∼ 10 cm kyr−1) (Anderson, 2001). For instance, according to Anderson (2001) minimum attenuation (i.e., < 5 %) is observed only when sedimentation rates exceed 50–70 cm kyr−1 – a range often observed at sedimentary drift sites, such as the Gardar Drift. Given the average sedimentation rate of ∼ 43 cm kyr−1 for the top 0–7.5 cm interval of our core (i.e., spanning an interval from ca. 1977 to 1800 CE) and sampling resolution of 0.5 cm, our ultimate chronological precision potentially achievable using these methods would be ∼ 12 years.
We further compare our revised age model based on anthropogenic signals with the previously published age model for Site GS06-144-09 MC-D. Figure 7 shows the 210Pb dates (Mjell et al., 2016), 14C dates, and information provided by the anthropogenic signals (i.e., 13C Suess-effect-derived tie points and the interval where the SCPs are present). The significant mismatch between the 210Pb and 14C dates once again highlights the need for independent approaches, as well as the potential of using anthropogenic signals to improve age model constraints over the last 2 centuries.
One of the main differences between our revised age model and that of Mjell et al. (2016) is the core-top age (1977 vs. 2006, respectively). This once more emphasizes the need to validate 210Pb-based chronologies as well as the common assumption of the year a sediment core was retrieved as the core-top age. Here we suggest and assume that the significant decline in our foraminiferal δ13C records over the last century is mainly caused by the oceanic 13C Suess effect. This is particularly the case for our G. bulloides δ13C, where the actual decline in foraminiferal δ13C is the same as the 13C Suess effect decline at 50 m depth. However, this may also be registered differently in other species. It is important to note that, although difficult to distinguish, our foraminiferal δ13C signals are also subject to natural climate variability. For instance, there are significant changes in the subpolar gyre circulation over the 20th century; more specifically, the observed productivity decline in the region (Spooner et al., 2020) will also be registered by our foraminiferal δ13C. Here, we have focused on the relative difference between average G. bulloides δ13C values over the industrial period vs. the core top (i.e., sharpest 13C decline due to the Suess effect) and demonstrated the potential utility of the 13C Suess effect approach in recent marine sediment chronologies. However, further sensitivity studies are needed to distinguish the effects of natural vs. anthropogenic climate variability in foraminiferal δ13C records.
The scale of the, ongoing, Suess effect is now starting to exceed the entire range of δ13C exhibited through most open-ocean environments (Eide et al., 2017), and, as such, it should be a dominant feature in records able to resolve short timescales. Indeed, the lack of this signal in a core-top record suggests that modern sediments were not recovered and/or that sedimentation rates and bioturbation may confound sub-centennial-scale interpretation of foraminiferal isotope records at a given core site.
Finally, Fig. 8 shows the sortable silt record of Mjell et al. (2016) on its original age model that is based on 210Pb and two 14C dates vs. the revised age model (as shown in Fig. 5b) for GS06-144-09 MC-D, plotted together with the AMV index (Gray et al., 2004), to illustrate how our proxy-based interpretations for the 20th century might change with revised marine sediment chronologies.
Although marine based uncertainties over the last 2 centuries might still be too high (∼ ± 18 years in average) for a significant lead–lag comparison with observational records, our new approach based on anthropogenic signals provides an independent and valuable first step in refining age models and for validation of existing age model approaches and their assumptions.
High-resolution (i.e., decadal to multi-decadal) marine sediment records from North Atlantic sedimentary drift sites are now emerging, with the potential to extend instrumental records further back in time, distinguish natural climate variability from anthropogenic variability, and contextualize current changes. However, age model uncertainties, particularly over the 20th century, pose major challenges, especially for integrating shorter instrumental records with those from extended marine archives. Recent sediments are dated using an array of methodologies, yet all have their own limitations (e.g., bomb carbon, local reservoir corrections for radiocarbon); they are either not applicable to all locations (e.g., tephrochronology) or can be below the detection limits and require another independent approach to confirm (e.g., 210Pb, 37Cs). Here we propose a new chronostratigraphic approach that uses anthropogenic signals to reduce age model uncertainties over the last 2 centuries. As a test application, we use the Gardar Drift sediment core GS06-144-09 MC and revise the age model at this site. Comparing planktonic δ13C records of GS06-144-09 MC with oceanic 13C Suess effect changes above the core location, we assign the beginning of the industrial period (i.e., 1800 CE) in our core and similarly derive the core-top age. We further use a combination of 14C AMS dates and the 13C Suess effect change estimates at our core location to calculate regional reservoir corrections at Gardar Drift. Our refined ΔR estimate for Gardar Drift (−69 ± 38 14C yr) is similar to the value obtained from the marine reservoir database when the 10 nearest points are used (−72 ± 64 14C yr) but with better uncertainty estimates. Furthermore, to validate our 13C Suess-effect-based age model we use another independent approach: spheroidal carbonaceous fly-ash particles (SCPs). The rapid increase in SCP concentrations after the 1950s at GS06-144-09 MC is consistent with the SCP trend in the region and our 13C Suess-effect-based age model. Our new approach, based on anthropogenic signals, provides an independent and valuable first step in refining age models and for validation of existing age model approaches and their assumptions.
Data are available as information files in the Supplement.
The supplement related to this article is available online at: https://doi.org/10.5194/gchron-6-449-2024-supplement.
NI and USN conceptualized the study. NI refined the new age model approach together with USN, FC, and AO. TLM processed the multicore samples and performed stable isotope analysis. NI processed samples for SCP analysis and conducted SCP analysis together with NLR and DJRT. USN led the efforts on stable isotope analysis. AO led the efforts on oceanic 13C Suess effect estimates for Gardar Drift. NI led the writing effort and coordinated input from all co-authors.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
We thank the crew of R/V G. O. Sars, the Institute of Marine Research (IMR), the University of Bergen, and the scientific party of UiB cruise no. GS06-144. We thank Marie Eide for her help with calculating oceanic 13C Suess effect estimates for our core location at Gardar Drift. Stable isotope data were produced at the Facility for Advanced Isotopic Research and Monitoring of Weather, Climate and Biogeochemical Cycling (FARLAB; Research Council of Norway grant 245907).
This study was funded by the Bjerknes Centre for Climate Research (BCCR) – Centre of Climate Dynamics (SKD) Strategic Project PARCIM (Proxy Assimilation for Reconstructing Climate and Improving Model). Nil Irvalı received additional support from the University of Bergen Meltzer Research Fund.
This paper was edited by Richard Staff and reviewed by James David Scourse and two anonymous referees.
Alves, E. Q., Macario, K., Ascough, P., and Bronk Ramsey, C.: The Worldwide Marine Radiocarbon Reservoir Effect: Definitions, Mechanisms, and Prospects, Rev. Geophys., 56, 278–305, https://doi.org/10.1002/2017RG000588, 2018.
Anderson, D. M.: Attenuation of millennial-scale events by bioturbation in marine sediments, Paleoceanography, 16, 352–357, https://doi.org/10.1029/2000PA000530, 2001.
Appleby, P. G.: Three decades of dating recent sediments by fallout radionuclides: a review, Holocene, 18, 83–93, https://doi.org/10.1177/0959683607085598, 2008.
Appleby, P. G., Piliposyan, G., and Hess, S.: Detection of a hot 137Cs particle in marine sediments from Norway: potential implication for 137Cs dating, Geo-Mar. Lett., 42, 2, https://doi.org/10.1007/s00367-021-00727-2, 2021.
Appleby, P. G., Piliposyan, G., Weckström, J., and Piliposian, G.: Delayed inputs of hot 137Cs and 241Am particles from Chernobyl to sediments from three Finnish lakes: implications for sediment dating, J. Paleolimnol., 69, 293–303, https://doi.org/10.1007/s10933-022-00273-6, 2023.
Ascough, P., Cook, G., and Dugmore, A.: Methodological approaches to determining the marine radiocarbon reservoir effect, Prog. Phys. Geogr., 29, 532–547, https://doi.org/10.1191/0309133305pp461ra, 2005.
Barsanti, M., Garcia-Tenorio, R., Schirone, A., Rozmaric, M., Ruiz-Fernández, A. C., Sanchez-Cabeza, J. A., Delbono, I., Conte, F., De Oliveira Godoy, J. M., Heijnis, H., Eriksson, M., Hatje, V., Laissaoui, A., Nguyen, H. Q., Okuku, E., Al-Rousan, S. A., Uddin, S., Yii, M. W., and Osvath, I.: Challenges and limitations of the 210Pb sediment dating method: Results from an IAEA modelling interlaboratory comparison exercise, Quat. Geochronol., 59, 101093, https://doi.org/10.1016/j.quageo.2020.101093, 2020.
Bindler, R., Renberg, I., Appleby, P. G., Anderson, N. J., and Rose, N. L.: Mercury accumulation rates and spatial patterns in lake sediments from west Greenland: a coast to ice margin transect, Environ. Sci. Technol., 35, 1736–1741, 2001.
Blaauw, M. and Christen, J. A.: Flexible paleoclimate age-depth models using an autoregressive gamma process, Bayesian Anal., 6, 457–474, https://doi.org/10.1214/ba/1339616472, 2011.
Boessenkool, K. P., Hall, I. R., Elderfield, H., and Yashayaev, I.: North Atlantic climate and deep-ocean flow speed changes during the last 230 years, Geophys. Res. Lett., 34, L13614, https://doi.org/10.1029/2007gl030285, 2007.
Booth, B. B. B., Dunstone, N. J., Halloran, P. R., Andrews, T., and Bellouin, N.: Aerosols implicated as a prime driver of twentieth-century North Atlantic climate variability, Nature, 484, 228–232, https://doi.org/10.1038/nature10946, 2012.
Chapman, M. R.: Seasonal production patterns of planktonic foraminifera in the NE Atlantic Ocean: Implications for paleotemperature and hydrographic reconstructions, Paleoceanography, 25, PA1101, https://doi.org/10.1029/2008pa001708, 2010.
Cléroux, C., Cortijo, E., Duplessy, J.-C., and Zahn, R.: Deep-dwelling foraminifera as thermocline temperature recorders, Geochem. Geophy. Geosy., 8, Q04N11, https://doi.org/10.1029/2006GC001474, 2007.
Eide, M., Olsen, A., Ninnemann, U. S., and Eldevik, T.: A global estimate of the full oceanic 13C Suess effect since the preindustrial, Global Biogeochem. Cy., 31, 492–514, https://doi.org/10.1002/2016GB005472, 2017.
England, J., Dyke, A. S., Coulthard, R. D., Mcneely, R., and Aitken, A.: The exaggerated radiocarbon age of deposit-feeding molluscs in calcareous environments, Boreas, 42, 362–373, https://doi.org/10.1111/j.1502-3885.2012.00256.x, 2013.
Field, D. B.: Variability in vertical distributions of planktonic foraminifera in the California Current: Relationships to vertical ocean structure, Paleoceanography, 19, PA2014, https://doi.org/10.1029/2003pa000970, 2004.
Forman, S. L. and Polyak, L.: Radiocarbon content of pre-bomb marine mollusks and variations in the 14C Reservoir age for coastal areas of the Barents and Kara Seas, Russia, Geophys. Res. Lett., 24, 885–888, https://doi.org/10.1029/97GL00761, 1997.
Gammon, R. H., Cline, J., and Wisegarver, D.: Chlorofluoromethanes in the northeast Pacific Ocean: Measured vertical distributions and application as transient tracers of upper ocean mixing, J. Geophys. Res.-Oceans, 87, 9441–9454, https://doi.org/10.1029/JC087iC12p09441, 1982.
Ganssen, G. M. and Kroon, D.: The isotopic signature of planktonic foraminifera from NE Atlantic surface sediments: implications for the reconstruction of past oceanic conditions, J. Geol. Soc. London, 157, 693–699, 2000.
Gebbie, G. and Huybers, P.: The Mean Age of Ocean Waters Inferred from Radiocarbon Observations: Sensitivity to Surface Sources and Accounting for Mixing Histories, J. Phys. Oceanogr., 42, 291–305, https://doi.org/10.1175/JPO-D-11-043.1, 2012.
Goldenberg, S. B., Landsea, C. W., Mestas-Nuñez, A. M., and Gray, W. M.: The Recent Increase in Atlantic Hurricane Activity: Causes and Implications, Science, 293, 474–479, https://doi.org/10.1126/science.1060040, 2001.
Graven, H. D.: Impact of fossil fuel emissions on atmospheric radiocarbon and various applications of radiocarbon over this century, P. Natl. Acad. Sci. USA, 112, 9542–9545, https://doi.org/10.1073/pnas.1504467112, 2015.
Gray, S. T., Graumlich, L. J., Betancourt, J. L., and Pederson, G. T.: A tree-ring based reconstruction of the Atlantic Multidecadal Oscillation since 1567 A.D, Geophys. Res. Lett., 31, L12205, https://doi.org/10.1029/2004GL019932, 2004.
Heaton, T. J., Köhler, P., Butzin, M., Bard, E., Reimer, R. W., Austin, W. E. N., Bronk Ramsey, C., Grootes, P. M., Hughen, K. A., Kromer, B., Reimer, P. J., Adkins, J., Burke, A., Cook, M. S., Olsen, J., and Skinner, L. C.: Marine20—The Marine Radiocarbon Age Calibration Curve (0–55,000 cal BP), Radiocarbon, 62, 779–820, https://doi.org/10.1017/RDC.2020.68, 2020.
Hinojosa, J. L., Moy, C. M., Prior, C. A., Eglinton, T. I., McIntyre, C. P., Stirling, C. H., and Wilson, G. S.: Investigating the influence of regional climate and oceanography on marine radiocarbon reservoir ages in southwest New Zealand, Estuar. Coast. Shelf S., 167, 526–539, https://doi.org/10.1016/j.ecss.2015.11.003, 2015.
Hughen, K. A.: Radiocarbon Dating of Deep-Sea Sediments, in: Proxies in Late Cenozoic Paleoceanography, edited by: Hillaire-Marcel, C. and de Vernal, A., Elsevier Science & Technology, 185–210, https://doi.org/10.1016/S1572-5480(07)01010-X, 2007.
Jonkers, L., van Heuven, S., Zahn, R., and Peeters, F. J. C.: Seasonal patterns of shell flux, δ18O and δ13C of small and large N. pachyderma (s) and G. bulloides in the subpolar North Atlantic, Paleoceanography, 28, 164–174, https://doi.org/10.1002/palo.20018, 2013.
Kaiser, J., Abel, S., Arz, H. W., Cundy, A. B., Dellwig, O., Gaca, P., Gerdts, G., Hajdas, I., Labrenz, M., Milton, J. A., Moros, M., Primpke, S., Roberts, S. L., Rose, N. L., Turner, S. D., Voss, M., and Ivar do Sul, J. A.: The East Gotland Basin (Baltic Sea) as a candidate Global boundary Stratotype Section and Point for the Anthropocene series, Anthropocene Review, 10, 25–48, https://doi.org/10.1177/20530196221132709, 2023.
Kaplan, A., Cane, M. A., Kushnir, Y., Clement, A. C., Blumenthal, M. B., and Rajagopalan, B.: Analyses of global sea surface temperature 1856–1991, J. Geophys. Res.-Oceans, 103, 18567–18589, https://doi.org/10.1029/97JC01736, 1998.
Keeling, C. D.: The Suess effect: 13Carbon-14Carbon interrelations, Environ. Int., 2, 229–300, https://doi.org/10.1016/0160-4120(79)90005-9, 1979.
Key, R. M., Kozyr, A., Sabine, C. L., Lee, K., Wanninkhof, R., Bullister, J. L., Feely, R. A., Millero, F. J., Mordy, C., and Peng, T. H.: A global ocean carbon climatology: Results from Global Data Analysis Project (GLODAP), Global Biogeochem. Cy., 18, GB4031, https://doi.org/10.1029/2004GB002247, 2004.
Lowe, D. J.: Tephrochronology and its application: A review, Quat. Geochronol., 6, 107–153, https://doi.org/10.1016/j.quageo.2010.08.003, 2011.
Mann, M. E., Steinman, B. A., and Miller, S. K.: Absence of internal multidecadal and interdecadal oscillations in climate model simulations, Nat. Commun., 11, 49, https://doi.org/10.1038/s41467-019-13823-w, 2020.
Mellon, S., Kienast, M., Algar, C., de Menocal, P., Kienast, S. S., Marchitto, T. M., Moros, M., and Thomas, H.: Foraminifera Trace Anthropogenic CO2 in the NW Atlantic by 1950, Geophys. Res. Lett., 46, 14683–14691, https://doi.org/10.1029/2019GL084965, 2019.
Miles, M. W., Divine, D. V., Furevik, T., Jansen, E., Moros, M., and Ogilvie, A. E. J.: A signal of persistent Atlantic multidecadal variability in Arctic sea ice, Geophys. Res. Lett., 41, 463–469, https://doi.org/10.1002/2013GL058084, 2014.
Mjell, T. L., Ninnemann, U. S., Kleiven, H. F., and Hall, I. R.: Multidecadal changes in Iceland Scotland Overflow Water vigor over the last 600 years and its relationship to climate, Geophys. Res. Lett., 43, 2111–2117, https://doi.org/10.1002/2016GL068227, 2016.
Moros, M., Andersen, T. J., Schulz-Bull, D., Häusler, K., Bunke, D., Snowball, I., Kotilainen, A., Zillén, L., Jensen, J. B., Kabel, K., Hand, I., Leipe, T., Lougheed, B. C., Wagner, B., and Arz, H. W.: Towards an event stratigraphy for Baltic Sea sediments deposited since AD 1900: approaches and challenges, Boreas, 46, 129–142, https://doi.org/10.1111/bor.12193, 2017.
Nyland, B. F., Jansen, E., Elderfield, H., and Andersson, C.: Neogloboquadrina pachyderma (dex. and sin.) Mg/Ca and δ18O records from the Norwegian Sea, Geochem. Geophy. Geosy., 7, Q10P17, https://doi.org/10.1029/2005GC001055, 2006.
Olsen, A. and Ninnemann, U.: Large δ13C Gradients in the Preindustrial North Atlantic Revealed, Science, 330, 658–659, https://doi.org/10.1126/science.1193769, 2010.
Otterå, O. H., Bentsen, M., Drange, H., and Suo, L. L.: External forcing as a metronome for Atlantic multidecadal variability, Nat. Geosci., 3, 688–694, https://doi.org/10.1038/ngeo955, 2010.
Pak, D. K. and Kennett, J. P.: A foraminiferal isotopic proxy for upper water mass stratification, J. Foramin. Res., 32, 319–327, https://doi.org/10.2113/32.3.319, 2002.
Pak, D. K., Lea, D. W., and Kennett, J. P.: Seasonal and interannual variation in Santa Barbara Basin water temperatures observed in sediment trap foraminiferal , Geochem. Geophy. Geosy., 5, Q12008, https://doi.org/10.1029/2004GC000760, 2004.
Perner, K., Moros, M., Jansen, E., Kuijpers, A., Troelstra, S. R., and Prins, M. A.: Subarctic Front migration at the Reykjanes Ridge during the mid- to late Holocene: evidence from planktic foraminifera, Boreas, 47, 175–188, https://doi.org/10.1111/bor.12263, 2018.
Perner, K., Moros, M., Otterå, O. H., Blanz, T., Schneider, R. R., and Jansen, E.: An oceanic perspective on Greenland's recent freshwater discharge since 1850, Sci. Rep.-UK, 9, 17680, https://doi.org/10.1038/s41598-019-53723-z, 2019.
Pieńkowski, A. J., Husum, K., Furze, M. F. A., Missana, A. F. J. M., Irvalı, N., Divine, D. V., and Eilertsen, V. T.: Revised ΔR values for the Barents Sea and its archipelagos as a pre-requisite for accurate and robust marine-based 14C chronologies, Quat. Geochronol., 68, 101244, https://doi.org/10.1016/j.quageo.2021.101244, 2021.
Reimer, P. J. and Reimer, R. W.: A Marine Reservoir Correction Database and On-Line Interface, Radiocarbon, 43, 461–463, https://doi.org/10.1017/S0033822200038339, 2001.
Reimer, P. J., Brown, T. A., and Reimer, R. W.: Discussion: Reporting and Calibration of Post-Bomb 14C Data, Radiocarbon, 46, 1299–1304, https://doi.org/10.1017/S0033822200033154, 2004.
Reimer, R. W. and Reimer, P. J.: An Online Application for ΔR Calculation, Radiocarbon, 59, 1623–1627, https://doi.org/10.1017/RDC.2016.117, 2017.
Reynolds, D. J., Scourse, J. D., Halloran, P. R., Nederbragt, A. J., Wanamaker, A. D., Butler, P. G., Richardson, C. A., Heinemeier, J., Eiríksson, J., Knudsen, K. L., and Hall, I. R.: Annually resolved North Atlantic marine climate over the last millennium, Nat. Commun., 7, 13502, https://doi.org/10.1038/ncomms13502, 2016.
Rose, N. L.: A note on further refinements to a procedure for the extraction of carbonaceous fly-ash particles from sediments, J. Paleolimnol., 11, 201–204, https://doi.org/10.1007/BF00686866, 1994.
Rose, N. L.: Quality control in the analysis of lake sediments for spheroidal carbonaceous particles, Limnol. Oceanogr.-Meth., 6, 172–179, https://doi.org/10.4319/lom.2008.6.172, 2008.
Rose, N. L.: Spheroidal Carbonaceous Fly Ash Particles Provide a Globally Synchronous Stratigraphic Marker for the Anthropocene, Environ. Sci. Technol., 49, 4155–4162, https://doi.org/10.1021/acs.est.5b00543, 2015.
Rose, N. L., Rose, C. L., Boyle, J. F., and Appleby, P. G.: Lake-Sediment Evidence for Local and Remote Sources of Atmospherically Deposited Pollutants on Svalbard, J. Paleolimnol., 31, 499–513, https://doi.org/10.1023/B:JOPL.0000022548.97476.39, 2004.
Rose, N. L., Jones, V. J., Noon, P. E., Hodgson, D. A., Flower, R. J., and Appleby, P. G.: Long-Range Transport of Pollutants to the Falkland Islands and Antarctica: Evidence from Lake Sediment Fly Ash Particle Records, Environ. Sci. Technol., 46, 9881–9889, https://doi.org/10.1021/es3023013, 2012.
Rubino, M., Etheridge, D. M., Trudinger, C. M., Allison, C. E., Battle, M. O., Langenfelds, R. L., Steele, L. P., Curran, M., Bender, M., White, J. W. C., Jenk, T. M., Blunier, T., and Francey, R. J.: A revised 1000 year atmospheric δ13C-CO2 record from Law Dome and South Pole, Antarctica, J. Geophys. Res.-Atmos., 118, 8482–8499, https://doi.org/10.1002/jgrd.50668, 2013.
Schiebel, R., Bijma, J., and Hemleben, C.: Population dynamics of the planktic foraminifer Globigerina bulloides from the eastern North Atlantic, Deep-Sea Res. Pt. I, 44, 1701–1713, https://doi.org/10.1016/S0967-0637(97)00036-8, 1997.
Schiebel, R., Waniek, J., Bork, M., and Hemleben, C.: Planktic foraminiferal production stimulated by chlorophyll redistribution and entrainment of nutrients, Deep-Sea Res. Pt. I, 48, 721–740, https://doi.org/10.1016/S0967-0637(00)00065-0, 2001.
Scourse, J. D., Wanamaker, A. D., Weidman, C., Heinemeier, J., Reimer, P. J., Butler, P. G., Witbaard, R., and Richardson, C. A.: The Marine Radiocarbon Bomb Pulse Across the Temperate North Atlantic: A Compilation of Δ14C Time Histories from Arctica Islandica Growth Increments, Radiocarbon, 54, 165–186, https://doi.org/10.2458/azu_js_rc.v54i2.16026, 2012.
Smith, J. N.: Why should we believe 210Pb sediment geochronologies?, J. Environ. Radioactiv., 55, 121–123, https://doi.org/10.1016/s0265-931x(00)00152-1, 2001.
Spero, H. J. and Lea, D. W.: Experimental determination of stable isotope variability in Globigerina bulloides: implications for paleoceanographic reconstructions, Mar. Micropaleontol., 28, 231–246, https://doi.org/10.1016/0377-8398(96)00003-5, 1996.
Spooner, P. T., Thornalley, D. J. R., Oppo, D. W., Fox, A. D., Radionovskaya, S., Rose, N. L., Mallett, R., Cooper, E., and Roberts, J. M.: Exceptional 20th Century Ocean Circulation in the Northeast Atlantic, Geophys. Res. Lett., 47, e2020GL087577, https://doi.org/10.1029/2020GL087577, 2020.
Stuiver, M. and Reimer, P. J.: A Computer Program for Radiocarbon Age Calibration, Radiocarbon, 28, 1022–1030, https://doi.org/10.1017/S0033822200060276, 1986.
Stuiver, M. and Reimer, P. J.: Extended 14C Data Base and Revised CALIB 3.0 14C Age Calibration Program, Radiocarbon, 35, 215–230, https://doi.org/10.1017/S0033822200013904, 1993.
Suess, H. E.: Radiocarbon Concentration in Modern Wood, Science, 122, 415–417, https://doi.org/10.1126/science.122.3166.415.b, 1955.
Sutton, R. T. and Hodson, D. L. R.: Atlantic Ocean Forcing of North American and European Summer Climate, Science, 309, 115–118, https://doi.org/10.1126/science.1109496, 2005.
Tanhua, T., Körtzinger, A., Friis, K., Waugh, D. W., and Wallace, D. W. R.: An estimate of anthropogenic CO2 inventory from decadal changes in oceanic carbon content, P. Natl. Acad. Sci. USA, 104, 3037–3042, https://doi.org/10.1073/pnas.0606574104, 2007.
Thornalley, D. J. R., Oppo, D. W., Ortega, P., Robson, J. I., Brierley, C. M., Davis, R., Hall, I. R., Moffa-Sanchez, P., Rose, N. L., Spooner, P. T., Yashayaev, I., and Keigwin, L. D.: Anomalously weak Labrador Sea convection and Atlantic overturning during the past 150 years, Nature, 556, 227–230, https://doi.org/10.1038/s41586-018-0007-4, 2018.
von Langen, P. J., Pak, D. K., Spero, H. J., and Lea, D. W.: Effects of temperature on Mg/Ca in neogloboquadrinid shells determined by live culturing, Geochem. Geophy. Geosy., 6, Q10P03, https://doi.org/10.1029/2005GC000989, 2005.
Waelbroeck, C., Lougheed, B. C., Vazquez Riveiros, N., Missiaen, L., Pedro, J., Dokken, T., Hajdas, I., Wacker, L., Abbott, P., Dumoulin, J.-P., Thil, F., Eynaud, F., Rossignol, L., Fersi, W., Albuquerque, A. L., Arz, H., Austin, W. E. N., Came, R., Carlson, A. E., Collins, J. A., Dennielou, B., Desprat, S., Dickson, A., Elliot, M., Farmer, C., Giraudeau, J., Gottschalk, J., Henderiks, J., Hughen, K., Jung, S., Knutz, P., Lebreiro, S., Lund, D. C., Lynch-Stieglitz, J., Malaizé, B., Marchitto, T., Martínez-Méndez, G., Mollenhauer, G., Naughton, F., Nave, S., Nürnberg, D., Oppo, D., Peck, V., Peeters, F. J. C., Penaud, A., Portilho-Ramos, R. d. C., Repschläger, J., Roberts, J., Rühlemann, C., Salgueiro, E., Sanchez Goni, M. F., Schönfeld, J., Scussolini, P., Skinner, L. C., Skonieczny, C., Thornalley, D., Toucanne, S., Rooij, D. V., Vidal, L., Voelker, A. H. L., Wary, M., Weldeab, S., and Ziegler, M.: Consistently dated Atlantic sediment cores over the last 40 thousand years, Scientific Data, 6, 165, https://doi.org/10.1038/s41597-019-0173-8, 2019.
Wang, C., Dong, S., Evan, A. T., Foltz, G. R., and Lee, S.-K.: Multidecadal Covariability of North Atlantic Sea Surface Temperature, African Dust, Sahel Rainfall, and Atlantic Hurricanes, J. Climate, 25, 5404–5415, https://doi.org/10.1175/JCLI-D-11-00413.1, 2012.
Zhang, R., Sutton, R., Danabasoglu, G., Kwon, Y.-O., Marsh, R., Yeager, S. G., Amrhein, D. E., and Little, C. M.: A Review of the Role of the Atlantic Meridional Overturning Circulation in Atlantic Multidecadal Variability and Associated Climate Impacts, Rev. Geophys., 57, 316–375, https://doi.org/10.1029/2019RG000644, 2019.