Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Climate changes modulated the history of Arctic iodine during the Last Glacial Cycle

Abstract

Iodine has a significant impact on promoting the formation of new ultrafine aerosol particles and accelerating tropospheric ozone loss, thereby affecting radiative forcing and climate. Therefore, understanding the long-term natural evolution of iodine, and its coupling with climate variability, is key to adequately assess its effect on climate on centennial to millennial timescales. Here, using two Greenland ice cores (NEEM and RECAP), we report the Arctic iodine variability during the last 127,000 years. We find the highest and lowest iodine levels recorded during interglacial and glacial periods, respectively, modulated by ocean bioproductivity and sea ice dynamics. Our sub-decadal resolution measurements reveal that high frequency iodine emission variability occurred in pace with Dansgaard/Oeschger events, highlighting the rapid Arctic ocean-ice-atmosphere iodine exchange response to abrupt climate changes. Finally, we discuss if iodine levels during past warmer-than-present climate phases can serve as analogues of future scenarios under an expected ice-free Arctic Ocean. We argue that the combination of natural biogenic ocean iodine release (boosted by ongoing Arctic warming and sea ice retreat) and anthropogenic ozone-induced iodine emissions may lead to a near future scenario with the highest iodine levels of the last 127,000 years.

Introduction

Atmospheric iodine, primarily emitted from oceans, forms new aerosol particles in the atmosphere1,2,3,4,5,6,7 and efficiently destroys ozone in the troposphere and lower stratosphere8,9,10,11, which reduces ozone radiative forcing12,13,14,15. In the Arctic environment, iodine has recently been identified as a significant source of cloud condensation nucleii (CCN), with the potential to influence clouds formation3. Iodine in the atmosphere is increasing globally (3-fold since the 1950s), as independently evidenced by polar and alpine ice core and tree ring measurements, following anthropogenic ozone pollution and global warming16,17,18. In the Arctic, a coastal ice core from Greenland revealed that ocean primary productivity controlled atmospheric iodine variability during the Holocene (i.e., last 11,700 years)19. Before the climatically stable Holocene, the Arctic underwent abrupt environmental changes featuring vast ice sheet growth and marked climate events during the Last Glacial Period (LGP) (115–11.7 kyr before 2000 CE (b2k)) and experienced temperatures warmer than today during the Last Interglacial Period (Eemian, 130–115 kyr b2k). However, to date, the natural evolution of Arctic iodine before the Holocene, in the absence of anthropogenic forcings, and its coupling with climate variability remain unknown.

In this work, we report the iodine levels in the Arctic during the Last Glacial Cycle (LGC) (last 127 kyr) using records from two Greenland ice cores (NEEM and RECAP) drilled in northwestern and eastern coastal Greenland, respectively (Fig. 1). Both ice cores provide continuous paleoenvironmental records back to the Eemian interglacial period20,21,22. We provide evidence that ocean-ice-atmosphere exchange of biogenic iodine in the Arctic Ocean dominated the atmospheric iodine budget during periods of rapid climatic changes. The results show that abrupt warming and sea ice retreat in the Arctic preceded maximum iodine levels, reaching peak concentrations during interglacial periods.

Fig. 1: Present-day (CE 2014) total iodine emission fluxes from the oceanic regions influencing Greenland.
figure 1

Modelled mean annual total iodine emission fluxes (CH3I + CH2I2 + CH2IBr + CH2ICl + HOI + I2) from the North Atlantic and Arctic Oceans using the 3D chemistry-climate model CAM-Chem. The figure includes the location of the NEEM and RECAP ice cores (red) and other ice cores and marine paleoceanographic archives in the Arctic and the northern North Atlantic discussed in the text and figures.

Results and discussion

Present-day iodine emission sources and depositional fluxes over Greenland

Atmospheric iodine in Greenland is controlled by a complex interplay involving iodine emissions from the oceans, as well as particle-bound iodine compounds related to dust and/or sea spray aerosol (ssa) variability23. The oceans, which are the largest reservoirs of iodine on Earth24, are known to be the dominant source of atmospheric iodine25,26,27. The major oceanic sources of gas-phase iodine are (i) inorganic, ozone-induced, hypoiodous acid and molecular iodine emissions and (ii) organic iodine from the metabolic activity of primary producers (phytoplankton and macro- and micro-algae). Furthermore, atmospheric iodine can also be related to mineral dust by promoting iodate stability in the ice matrix and/or favouring the adsorption of iodine gas-phase molecules to dust particles during atmospheric transport10. On the other hand, ssa expelled from the ocean surface during wave breaking incorporate iodine in their composition and are a substrate where gaseous iodine species undergo heterogeneous recycling reactions28,29,30.

We used a global chemistry climate model17 to understand present-day iodine levels in Greenland (see Methods). Present iodine depositional fluxes at the RECAP drilling site (71°18′ N; 26°43′ W; 2315 m a.s.l.) are ~51% higher than at NEEM (77°45′ N, 51°06′ W; 2484 m a.s.l.) (Supplementary Fig. 1). This is explained by the different ocean iodine emission strengths at the source regions of air masses that influence both sites (the Canadian Arctic, Baffin Bay and Hudson Bay for NEEM31 and the North Atlantic and Arctic Ocean for RECAP32). Note also that RECAP is closer to the coast (30 km from the North Atlantic Ocean) than NEEM (350 kms from the Baffin Bay). Our model results show higher iodine emissions in the RECAP source regions (50N–80N, 45W–12E; mean iodine emissions of 9.1 × 10−14 kg m−2 s−1) than in the NEEM source areas (50N–80N, 120W–45E; mean iodine emissions of 3.5 × 10−14 kg m−2 s−1) (Fig. 1).

The Greenland iodine records

Total iodine and sodium ([I], [Na]) concentrations from NEEM and RECAP ice core samples were determined at 1.10 m and 0.55 m mean intervals of melted ice, respectively (Supplementary Figs. 2-4). Calcium ([Ca]) concentrations were also used as proxies of mineral dust in Greenland ice cores33,34. These geochemical measurements along with other ice cores and marine paleoceanographic archives allow the reconstruction of Arctic iodine variability and environmental sources during the LGC.

Maximum iodine levels were found during the Eemian in the NEEM ice core, with a total mean iodine depositional flux (Iflux) of 9.6 µg m−2 yr−1 coinciding with the highest Arctic Ocean primary productivity inferred from the phytoplankton biomarker brassicasterol35 (Fig. 2). Our analysis suggests that these high Eemian iodine levels resulted from the interplay of (i) enhanced biogenic iodine emissions driven by primary productivity; (ii) higher sea surface temperatures (SSTs) in the Arctic and the North Atlantic36,37, which facilitates sea-air transfer of volatiles; and (iii) dominant open water conditions in the Arctic, as reflected by the low PBIP25 values. PBIP25 has been used as a proxy of the intensity of sea ice cover in the Arctic38, as well as in Baffin Bay31,39 since it accounts for the algal and phytoplankton (brassicasterol) biomarkers activity, which allows to assess the spatial and temporal extent of the sea ice cover40 (Fig. 2). The Arctic sea ice extent was significantly reduced during the Eemian, and minimum sea ice concentrations towards almost ice-free summers occurred in wide areas as far north as the northern Barents Sea as a consequence of the strong inflow of warm Atlantic waters38, thus facilitating springtime blooms and oceanic iodine production. Indeed, summer North Atlantic SSTs during the Eemian were 2 °C greater than at present36,37. Furthermore, a mean global sea level 6–8 m above the present level41,42,43 (Fig. 2) triggered the flooding of most parts of the shallow Siberian marginal seas, eventually resulting in more productive coastal areas in the Arctic. At NEEM, the highest Iflux was recorded at ~122 kyr b2k concomitant with maximum levels of biological productivity, sea level rise and SST. The coupling of all these drivers during the Eemian would have eventually released large amounts of iodine into the Arctic atmosphere, resulting in the largest iodine emissions on record. The significant decrease in NEEM and RECAP Iflux recorded during the last glacial inception at ~120 kyr b2k is correlated with an increase in sea ice31,38 and a decrease in both SST and brassicasterol (Fig. 2), as also indicated by independent modelling experiments that resemble pre-industrial scenarios41.

Fig. 2: Iodine levels measured in the NEEM and RECAP ice cores together with other paleoenvironmental proxies for the last 130 kyr b2k.
figure 2

From bottom to top: a, b iodine depositional fluxes (Iflux) from NEEM and RECAP ice cores; c, d Arctic primary productivity and sea ice evolution (brassicasterol and PBPIP25 profiles, respectively38,46); e North Atlantic summer sea surface temperature (SST)36; f North Atlantic relative sea level changes88; g Greenland air temperature reconstruction from NGRIP64,65; h solar irradiation reconstruction in July at 77° (NEEM) and 71°N (RECAP)89.

Ice core iodine showed reduced values throughout the LGP, with mean NEEM and RECAP Iflux of 3.8 and 5.1 µg m−2 yr−1, respectively (Fig. 2). The LGP is characterized by abrupt climatic variability marked by Dansgaard/Oeschger (D/O) events. These climatic phases consist of warm Greenland interstadials and cold Greenland stadials that are well recorded in Greenland ice cores and marine sediments from the North Atlantic region and the Nordic seas44,45,46,47,48,49. To investigate in detail the iodine variability across the rapid D/O transitions, additional high-resolution iodine measurements were carried out in the NEEM ice core from 34 to 42 kyr b2k (D/O 7 to 10), with a 3- to 10-year resolution (Fig. 3 and Supplementary Fig. 4). In order to differentiate the main sources of iodine that determine Iflux variability, we have decoupled oceanic iodine emissions (Iocean) from dust-related iodine (Idust) (‘Methods’). The sub-decadal resolution measurements of iodine during this period show the concomitant evolution of Iocean, Arctic Ocean primary productivity and sea ice variability46,48,49 (Fig. 3). Each D/O event showed a similar sequence: (i) iodine emissions increased concurrently with a reduction in sea ice at the onset of the interstadials; (ii) iodine emissions decreased during the late interstadial periods as the climate cooled and sea ice gradually increased; and iii) negligible iodine was released from the oceans with the development of thicker sea ice conditions during the cold stadial periods.

Fig. 3: High-resolution iodine level variability in Greenland NEEM and RECAP ice cores from 34 to 42 kyr b2k (=D/O 7 to 10).
figure 3

From bottom to top: a high-resolution NEEM (dark yellow) and RECAP (red) iodine depositional fluxes (Iflux), b, c dust-related iodine (Idust) and reconstructed ocean iodine emissions (Iocean), with dark yellow lines representing the 11-sample running average in NEEM; d, e Arctic primary productivity and f, g sea ice evolution (brassicasterol and PBPIP25 profiles, respectively) from sediment cores in the Norwegian Sea48,49; h Greenland air temperature reconstruction from NGRIP64. Blue and white bands represent Greenland interstadials (GIs) and stadials (GSs), respectively.

Interestingly, NEEM and RECAP Iflux and Iocean do not follow the same trend during D/O events (Fig. 3). In the context of strongly reduced iodine emissions during the LGP due to the extensive development of perennial sea ice in the Canadian Arctic31, a significant fraction of iodine registered in the NEEM ice core was strongly controlled by dust and, to a lesser extent, by ssa. This is corroborated by the robust statistical correlations between [I] and [Ca] and between [I] and [Na] from D/O 7–10 (r > 0.75 and 0.78, respectively). These correlations are also seen in RECAP during the LGP (r > 0.88 and 0.87, respectively). Sea spray aerosol emissions (indicated by [Na] concentrations in both ice cores) show a 50% increase during the LGP33 (Supplementary Figs. 2 and 3), thereby contributing to higher deposition of ssa-related iodine during this period. On the other hand, [Ca] concentrations in Greenland ice cores show a major increase during stadials33 (Supplementary Figs. 2 and 3) due to enhanced mineral dust associated with increased aridity in East Asia and enhanced atmospheric transport mechanisms33,50,51. The higher Idust values found in NEEM and RECAP during the coldest phases of the LGP are in agreement with the highest iodine fluxes recorded in the Talos Dome ice core (eastern Antarctica) between 16.8 and 33.8 kyr b2k10. This Antarctic ice core reported a correlation between iodate, the most stable iodine species in the atmosphere, and atmospheric dust, indicating uptake of gas-phase iodine molecules into fine dust particles during atmospheric transport10.

The largest iodine variability occurred during the LGP-Holocene transition when iodine emissions strongly fluctuated following the main climatic phases of the Last Glacial Termination. NEEM and RECAP Iocean showed very similar evolution patterns, with minimum Iocean values during Heinrich event 1 (H1, 16.8 kyr b2k) (Fig. 4), when surface waters cooled and a thick perennial sea ice cover reached both the Norwegian Sea and Baffin Bay31,46,52. The Iocean increased in both ice cores at the onset of the Bølling (14.7–14.1 kyr b2k) warming in the Arctic, slightly decreased during the Intra-Allerød Cold Period (14.1–13.9 kyr b2k) and increased again during the Allerød (13.9–12.9 kyr b2k) warm period (Fig. 4). The increase in Iocean values during the Bølling-Allerød period was synchronous with an abrupt increase in ocean productivity as far north as the Fram Strait53. NEEM and RECAP Iocean moderately decreased during the last cold spell of the Last Glacial Termination, the Younger Dryas (12.9–11.7 kyr b2k).

Fig. 4: Iocean variability in Greenland ice cores during the Last Glacial Termination/Holocene transition.
figure 4

From bottom to top: a Iocean values in NEEM (dark yellow) and RECAP (red) ice cores; b, c Arctic primary productivity and sea ice evolution (brassicasterol and PBPIP25 profiles, respectively) in Baffin Bay (core GeoB19927)60 and the Norwegian Sea (core MD99-2284)46; d Greenland air temperature reconstruction from NGRIP64. Bands indicate the main cold (blue) and warm (red) climatic phases of this period (IACP = Intra-Allerød Cold Period).

Oceanic iodine emissions increased in both ice cores at the onset of the Holocene (~11 kyr b2k; Fig. 4) resulting from the interplay of different factors: (i) enhanced marine primary productivity in the Arctic46,54,55,56; (ii) dominant open water conditions in the subpolar North Atlantic46; (iii) abrupt sea level rise (~120 m) during the Last Glacial Termination57,58 resulting in more shallow oceanic domains increasing iodine emission from coastal areas; and iv) maximum levels of solar irradiance in the Arctic that would have also increased oxidative stress in algae, leading to an enhancement of biological iodine production in the ocean and subsequent release into the Arctic atmosphere59. Iflux strongly increased in RECAP from 11 to 9 kyrs b2k, while the NEEM Iflux increase was significantly more gradual at that time (Fig. 2), most likely due to the lower productivity and delayed sea ice retreat in Baffin Bay during the early Holocene60 (Fig. 4). Indeed, this area remained densely sea ice covered until ~7.4 kyr b2k, when warmer North Atlantic waters penetrated Baffin Bay, increasing regional bioproductivity61 (Fig. 4). The mean Holocene iodine flux in NEEM is 57.9% of that in RECAP (mean iodine fluxes of 5.1 and 8.8 µg m−2 yr-1, respectively), in agreement with present-day reactive iodine deposition at both sites obtained from atmospheric chemistry transport model results (Supplementary Fig. 1).

Environmental implications

Ocean-emitted iodine leads to new particle formation (NPF) in the atmosphere5,6. Recent field observations have demonstrated that iodine is a significant driver of NPF in the central Arctic3 and is currently considered a potentially relevant source of CCN3—from which marine clouds originate—that scatter incoming radiation and can contribute a cooling effect to the Earth´s radiation budget62. The significant present-day contribution of iodine to NPF leads us to speculate that, in the absence of other human-induced NPF nucleation drivers during pre-industrial times63, oceanic iodine may have played an even larger relative role in total NPF nucleation in a pristine Arctic during the LGC.

We find that the time lag between temperature changes and iodine maximum emission values recorded in Greenland ice cores throughout the LGC strongly depends on the atmosphere-ice-ocean interactions driving warming SST, sea ice thinning rate and algae metabolic activity that eventually control biogenic iodine emissions. Indeed, the highest iodine emission periods during the LGC were preceded by the highest air temperatures that would have first enhanced SST and ocean biological activity and later favoured the transfer of volatile iodine compounds from the ocean to the atmosphere. Together with the efficiency of iodine to destroy ozone and activate CCN in the polar atmosphere, this points to a possible contribution of iodine to ozone loss and new particle formation, and to their associated radiative impacts, during past abrupt climate change periods, such as D/O events and the onset of interglacials. However, this potential contribution of iodine to past ozone loss and atmospheric NPF is currently unknown and warrants further research.

We now turn to the use of past warm intervals recorded in Greenland ice cores as analogues to evaluate potential feedback and thresholds in the Arctic climate system and their implications for future warming impacts. The rate of ongoing atmospheric warming resembles only the abrupt climate changes and sea ice retreat that occurred during the D/O events64, with climate warming of 5–16.5 °C occurring within a century65. The high resolution NEEM iodine dataset during D/O 7–10 show empirical evidence of a fast response of iodine emissions synchronously with the rapid North Atlantic sea ice retreat preceding abrupt Greenland warming during D/O events48. Therefore, one question arising is if Arctic iodine levels are facing this scenario in the near future?

Iodine emissions and depositional fluxes during the Eemian and/or Holocene Thermal Maximum shed light on the possible environmental conditions in an expected ice-free Arctic Ocean and thus remain useful as an observational constraint on projections of future impacts. Thus, while atmospheric iodine depositional fluxes recorded in RECAP during the Holocene Thermal Maximum19 mirror present-day conditions17, iodine levels recorded in NEEM during the mid-Eemian tripled present-day values at this location. Similar trends might be expected in near future environmental scenarios including (i) predicted ice-free summertime conditions in the Arctic Ocean by 2050 CE66 and (ii) 3–4 °C warmer global temperatures with a three- to four-fold amplification in the Arctic by 2100 CE67. However, one important point to consider is that the predicted iodine emissions in the Arctic in the near future based on past analogues represent only the biogenic fraction of the total foreseen atmospheric iodine concentrations in the Arctic. Ozone-induced inorganic iodine emissions during these past warm periods could be considered negligible in the context of the very low ozone concentrations before industrialization19,68. In contrast, at present, ocean release of inorganic iodine is estimated to account for 75% of the total source of atmospheric iodine, with the remaining 25% coming from biogenic emissions69. Furthermore, ozone-related inorganic iodine emissions are predicted to increase during the 21st century15. Therefore, we argue that the future combined (inorganic + biogenic) ice-free Arctic Ocean iodine emissions might not have analogues in the past, which could potentially lead to a near future scenario with the highest iodine levels of the last 127,000 years.

In summary, this study uncovers how past climate changes modulated the variability of the ocean-atmosphere exchange of iodine in the Arctic during the LGC. We show that ocean productivity was the dominant source of Arctic iodine during interglacial periods, leading to highest levels during the Holocene and the Eemian, while dust-related iodine constituted a significant source of iodine in the Arctic atmosphere during the coldest phases of the Last Glacial Period. The mirrored trends of iodine emissions and ocean primary productivity proxies during different past environmental scenarios highlight the key role of ocean biology and sea ice cover in driving iodine variability in the Arctic region. Finally, we conclude that understanding the past variability of iodine as a key environmental element is important to further comprehending the evolution of Arctic atmospheric chemistry on centennial to millennial timescales and its coupling to climate changes.

Methods

Ice core chronologies and sampling resolution

NEEM age-depth model was achieved by transferring the annual layer counting from the NGRIP to the NEEM ice cores using 787 tie points (mainly volcanic layers)21. NEEM samples were collected at 1.10 m mean intervals of melted ice where 10 mL of meltwater was collected in an acid-cleaned polyethylene bottle and immediately frozen. Chronological resolution between samples (n = 84) ranges from 4 to 267 years during the Holocene. There is a sampling gap between 3-9 ky b2k where the ice was not analysed due to the poor quality of the core in the “brittle ice” section31. The mean sampling frequency during the Glacial Period (n = 247) is ~400 years, where the highest resolution achieved was between 25 and 50 kyrs b2k with sampling resolution ranging from 60 to 150 years during warmer interstadial periods and colder glacial stadial periods, respectively. Additional sampling was carried out during the intervals 34.5–35.8, 37.6–38.8 and 39.4–41.6 (n = 1304) providing a time resolution ranging from 3 to 10 years during those intervals to better resolve the rapid stadial-interstadial transition.

RECAP age-depth model is based on (i) annual layer counting from the last four millennia using the StratiCounter algorithm70 constrained by volcanic eruption markers and synchronized to the GICC05 Greenland Ice Core Chronology framework71; (ii) a modified Dansgaard-Johnsen ice flow model72 constrained to well-dated age markers from 4048 to 11,703 ka b2k and (ii) linear interpolation between 73 GICC05-modelext age markers between 11,703 and 120,215 years b2k32. Each RECAP sample integrated a 55 cm depth interval, where the resolution in the Holocene (n = 1035) ranged from sub-annual in the upper metres to decadal or centennial resolution during the mid to early Holocene (average resolution of 12 years per sample). Due to the thinning of annual layers according to depth, each sample integrates approximately 470 years during the glacial period (n = 218)32. For details on the age model we refer the readers to Simonsen et al.22.

Geochemical analyses

Samples from NEEM and RECAP ice cores were collected using continuous flow analysis systems. NEEM samples were sent to the Environmental Analytical Chemistry laboratory of the CNR-ISP and Ca’ Foscari University of Venice for iodine and sodium analyses. Total iodine and sodium concentrations were determined by inductively coupled plasma mass spectrometry (ICP-MS). Sodium and iodine detection limits were 1 ppb and 0.005 ppb, respectively. RECAP samples were sent to ISP-CNR and to Curtin University of Technology (Perth, Australia). Total iodine and sodium concentrations in Perth were determined by inductively coupled plasma-sector field mass spectrometry (ICP-SFMS) with sodium and iodine detection limits of 1.1 ppb and 0.002 ppb, respectively. The analytical procedure for both records was performed following Vallelonga et al.73 and Maffezzoli et al.32. Instrumental errors for iodine and sodium concentrations were 5%.

Iflux, Idust and Iocean calculation

Total iodine depositional fluxes (Iflux) in the ice cores were calculated following Corella et al.19. To estimate Iflux in the RECAP ice core, we re-scaled NGRIP accumulation rates65 by a factor of 436/174 since it is the modern accumulation ratio at the NGRIP and RECAP drilling sites74,75,76. Dust-related iodine (Idust) and oceanic gas-phase iodine emissions (Iocean) were calculated according to Eqs. (1) and (2)

$${{{{{{\rm{I}}}}}}_{{{{{\rm{dust}}}}}}}={{{{{\rm{nssCa}}}}}}\,{{{{{\rm{X}}}}}}([{{{{{\rm{I}}}}}}]{/}[{{{{{\rm{Ca}}}}}}])_{\left.{{{{{\rm{crust}}}}}}\right)}{/}[{{{{{\rm{I}}}}}}]$$
(1)
$${{{{{{\rm{I}}}}}}_{{{{{\rm{ocean}}}}}}}=({{{{{\rm{nssI}}}}}}-({{{{{\rm{nssCa}}}}}}\,{{{{{\rm{X}}}}}}{([{{{{{\rm{I}}}}}}]{/}[{{{{{\rm{Ca}}}}}}])}_{{{{{{\rm{crust}}}}}}})){/}[{{{{{\rm{I}}}}}}]$$
(2)

with non-sea salt iodine (nssI) representing the excess iodine production beyond the iodine production related to ssa. It is calculated as nssI = [I] − [Na] × ([I]/[Na])seawater with a [I]/[Na] seawater concentration ratio of 5.6 × 10−6 77. Non-sea salt calcium (nssCa) represents the non-sea spray aerosol fraction of calcium and is representative of mineral dust input. It is calculated as nssCa = [Ca] − ssCa, where ssCa = [Na] × ([Ca]/[Na])seawater with a [Ca]/[Na] seawater concentration ratio of 0.0391. ([I]/[Ca])crust represents the ratio of both elements in the continental crust (4.75 × 10−5 78). Both, Iocean and Idust correspond to unitless fractional proxies.

Atmospheric chemistry modelling

We also used the global 3-D chemistry-climate Community Atmospheric Model with chemistry (CAM-Chem) version 4 to estimate (i) oceanic iodine emissions from different source areas in the Arctic and North Atlantic Oceans (Fig. 1) and (ii) the deposition fluxes of reactive iodine (HI + HOI + I2O2 + I2O3 + I2O4 + INO2 + IONO2) in Greenland (Supplementary Fig. 1). CAM-Chem includes a comprehensive benchmark chemistry scheme to simulate the evolution of trace gases and aerosols in the troposphere and stratosphere79,80. The model implements a halogen chemistry scheme for chlorine, bromine and iodine12,81,82. This includes the photochemical breakdown of five very short-lived bromocarbons (VSLBr = CHBr3, CH2Br2, CH2BrCl, CHBrCl2, CHBr2Cl) and four iodocarbons (VSLI = CH3I, CH2ICl, CH2IBr, CH2I2), which are naturally emitted from the ocean into the atmosphere83. Additionally, abiotic oceanic sources of HOI and I2 were included in the lowest layer of the model26 based on laboratory studies of the oxidation of aqueous iodide by atmospheric ozone deposited on the ocean surface25,84. In this simulation, the model was configured with a horizontal resolution of 1.9° latitude by 2.5° longitude and 26 levels, from the surface to 40 km (with eight levels above 100 hPa)17,29,82. At the lower boundary, the time-varying zonally averaged distributions of CO2, CH4, H2, N2O and long-lived halocarbons (CFC-11, CFC-12, CFC-113, HCFC-22, H-1211, H-1301, CCl4, CH3CCl3, CH3Cl and CH3Br) are specified following their observed distribution for 2000 CE85. Monthly mean time variations observed for sea surface temperature and sea ice distribution are also prescribed. To obtain a reasonable representation of the overall stratospheric circulation, the integrated momentum that would have been deposited above the model top is specified by an upper boundary condition. We used the output of a previous work17 for CE 2010. This simulation was run in free-running mode79 with prescribed sea surface temperatures, sea ice and meteorological fields from 1950 to 201086. The simulated dynamics and transport therefore represent the daily synoptic conditions of the observations in CE 2010, allowing the online coupling between the ocean, ice and atmospheric modules of the CESM model79. Abiotic oceanic sources of iodine (HOI and I2) are modelled according to the meteorology (temperature, winds, sea surface temperature and surface pressure) in CE 2010 based on MacDonald et al.84 parameterization. The authors would like to remark that a quantitative assessment of the role of iodine in aerosol radiative forcing (through its influence in NPF in the past) is yet not feasible mainly because there are still remaining chemical mechanistic gaps in our knowledge of the iodine particle formation and growth4,87.

Data availability

The ice core iodine data generated in this study have been deposited in the Zenodo database (https://doi.org/10.5281/zenodo.5721369).

Code availability

The software code for the CESM model is available from http://www.cesm.ucar.edu/models/.

References

  1. Allan, J. et al. Iodine observed in new particle formation events in the Arctic atmosphere during ACCACIA. Atmos. Chem. Phys. 15, 5599–5609 (2015).

    ADS  CAS  Google Scholar 

  2. Roscoe, H. K. et al. Particles and iodine compounds in coastal Antarctica. J. Geophys. Res.: Atmos. 120, 7144–7156 (2015).

    ADS  CAS  Google Scholar 

  3. Baccarini, A. et al. Frequent new particle formation over the high Arctic pack ice by enhanced iodine emissions. Nat. Commun. 11, 1–11 (2020).

    Google Scholar 

  4. He, X.-C. et al. Role of iodine oxoacids in atmospheric aerosol nucleation. Science 371, 589–595 (2021).

    ADS  CAS  PubMed  Google Scholar 

  5. O’Dowd, C. D. et al. Marine aerosol formation from biogenic iodine emissions. Nature 417, 632–636 (2002).

    ADS  PubMed  Google Scholar 

  6. Sipilä, M. et al. Molecular-scale evidence of aerosol particle formation via sequential addition of HIO3. Nature 537, 532–534 (2016).

    ADS  PubMed  PubMed Central  Google Scholar 

  7. Gómez Martín, J. C. et al. A gas-to-particle conversion mechanism helps to explain atmospheric particle formation through clustering of iodine oxides. Nat. Commun. 11, 4521 (2020).

    ADS  PubMed  PubMed Central  Google Scholar 

  8. Saiz-Lopez, A. et al. Boundary layer halogens in coastal Antarctica. Science 317, 348–351 (2007).

    ADS  CAS  PubMed  Google Scholar 

  9. Sherwen, T. et al. Iodine’s impact on tropospheric oxidants: a global model study in GEOS-Chem. Atmos. Chem. Phys. 16, 1161–1186 (2016).

    ADS  CAS  Google Scholar 

  10. Spolaor, A. et al. Halogen species record Antarctic sea ice extent over glacial–interglacial periods. Atmos. Chem. Phys. 13, 6623–6635 (2013).

    ADS  Google Scholar 

  11. Koenig, T. K. et al. Quantitative detection of iodine in the stratosphere. Proc. Natl Acad. Sci. USA 117, 1860–1866 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  12. Saiz-Lopez, A. et al. Estimating the climate significance of halogen-driven ozone loss in the tropical marine troposphere. Atmos. Chem. Phys. 12, 3939–3949 (2012).

    ADS  CAS  Google Scholar 

  13. Hossaini, R. et al. Efficiency of short-lived halogens at influencing climate through depletion of stratospheric ozone. Nat. Geosci. 8, 186 (2015).

    ADS  CAS  Google Scholar 

  14. Sherwen, T., Evans, M. J., Carpenter, L. J., Schmidt, J. A. & Mickley, L. J. Halogen chemistry reduces tropospheric O 3 radiative forcing. Atmos. Chem. Phys. 17, 1557–1569 (2017).

    ADS  CAS  Google Scholar 

  15. Iglesias-Suarez, F. et al. Natural halogens buffer tropospheric ozone in a changing climate. Nat. Clim. Change 10, 147–154 (2020).

    ADS  CAS  Google Scholar 

  16. Legrand, M. et al. Alpine ice evidence of a three-fold increase in atmospheric iodine deposition since 1950 in Europe due to increasing oceanic emissions. Proc. Natl Acad. Sci. USA 115, 12136–12141 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  17. Cuevas, C. A. et al. Rapid increase in atmospheric iodine levels in the North Atlantic since the mid-20th century. Nat. Commun. 9, 1452 (2018).

    ADS  PubMed  PubMed Central  Google Scholar 

  18. Zhao, X., Hou, X. & Zhou, W. Atmospheric iodine (127I and 129I) record in spruce tree rings in the Northeast Qinghai-Tibet Plateau. Environ. Sci. Technol. 53, 8706–8714 (2019).

    ADS  CAS  PubMed  Google Scholar 

  19. Corella, J. P. et al. Holocene atmospheric iodine evolution over the North Atlantic. Climate 15, 2019–2030 (2019).

    Google Scholar 

  20. Dahl-Jensen, D. et al. Eemian interglacial reconstructed from a Greenland folded ice core. Nature 493, 489 (2013).

    ADS  CAS  Google Scholar 

  21. Rasmussen, S. O. et al. A first chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core. Climate 9, 2713–2730 (2013).

    Google Scholar 

  22. Simonsen, M. F. et al. East Greenland ice core dust record reveals timing of Greenland ice sheet advance and retreat. Nat. Commun. 10, 1–8 (2019).

    CAS  Google Scholar 

  23. Saiz-Lopez, A. et al. Atmospheric chemistry of iodine. Chem. Rev. 112, 1773–1804 (2012).

    CAS  PubMed  Google Scholar 

  24. Fuge, R. & Johnson, C. C. Iodine and human health, the role of environmental geochemistry and diet, a review. Appl. Geochem. 63, 282–302 (2015).

    CAS  Google Scholar 

  25. Carpenter, L. J. et al. Atmospheric iodine levels influenced by sea surface emissions of inorganic iodine. Nat. Geosci. 6, 108 (2013).

    ADS  CAS  Google Scholar 

  26. Prados-Roman, C. et al. A negative feedback between anthropogenic ozone pollution and enhanced ocean emissions of iodine. Atmos. Chem. Phys. 15, 2215–2224 (2015).

    ADS  CAS  Google Scholar 

  27. Carpenter, L. J. et al. Marine iodine emissions in a changing world. Proc. R. Soc. A: Math., Phys. Eng. Sci. 477, 20200824 (2021).

    ADS  Google Scholar 

  28. Tham, Y. J. et al. Direct field evidence of autocatalytic iodine release from atmospheric aerosol. Proc. Natl Acad. Sci. USA 118, e2009951118 (2021).

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Saiz-Lopez, A. et al. Iodine chemistry in the troposphere and its effect on ozone. Atmos. Chem. Phys. 14, 13119–13143 (2014).

    ADS  Google Scholar 

  30. Gómez Martín, J. C. et al. Spatial and temporal variability of iodine in aerosol. J. Geophys. Res.: Atmos. 126, e2020JD034410 (2021).

    ADS  Google Scholar 

  31. Spolaor, A. et al. Canadian Arctic sea ice reconstructed from bromine in the Greenland NEEM ice core. Sci. Rep. 6, 33925 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  32. Maffezzoli, N. et al. A 120 000-year record of sea ice in the North Atlantic? Climate 15, 2031–2051 (2018).

    Google Scholar 

  33. Schüpbach, S. et al. Greenland records of aerosol source and atmospheric lifetime changes from the Eemian to the Holocene. Nat. Commun. 9, 1476 (2018).

    ADS  PubMed  PubMed Central  Google Scholar 

  34. Simonsen, M. F. et al. Particle shape accounts for instrumental discrepancy in ice core dust size distributions. Climate 14, 601–608 (2018).

    Google Scholar 

  35. Hu, A. et al. Influence of Bering Strait flow and North Atlantic circulation on glacial sea-level changes. Nat. Geosci. 3, 118 (2010).

    ADS  CAS  Google Scholar 

  36. Kandiano, E. S., Bauch, H. A. & Müller, A. Sea surface temperature variability in the North Atlantic during the last two glacial–interglacial cycles: comparison of faunal, oxygen isotopic, and Mg/Ca-derived records. Palaeogeogr., Palaeoclimatol., Palaeoecol. 204, 145–164 (2004).

    Google Scholar 

  37. Bauch, H. A. & Kandiano, E. S. Evidence for early warming and cooling in North Atlantic surface waters during the last interglacial. Paleoceanography 22, PA1201 (2007).

  38. Stein, R., Fahl, K., Gierz, P., Niessen, F. & Lohmann, G. Arctic Ocean sea ice cover during the penultimate glacial and the last interglacial. Nat. Commun. 8, 373 (2017).

    ADS  PubMed  PubMed Central  Google Scholar 

  39. Andrews, J. et al. Land/ocean correlations during the last interglacial/glacial transition, Baffin Bay, northwestern North Atlantic: A review. Quat. Sci. Rev. 4, 333–355 (1985).

    ADS  Google Scholar 

  40. Müller, J. et al. Towards quantitative sea ice reconstructions in the northern North Atlantic: a combined biomarker and numerical modelling approach. Earth Planet. Sci. Lett. 306, 137–148 (2011).

    ADS  Google Scholar 

  41. Otto-Bliesner, B. L. et al. How warm was the last interglacial? New model–data comparisons. Philosophical Transactions of the Royal Society A: Mathematical. Phys. Eng. Sci. 371, 20130097 (2013).

    Google Scholar 

  42. McFarlin, J. M. et al. Pronounced summer warming in northwest Greenland during the Holocene and Last Interglacial. Proc. Natl Acad. Sci. 115, 6357–6362 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  43. Siddall, M. et al. Sea-level fluctuations during the last glacial cycle. Nature 423, 853–858 (2003).

    ADS  CAS  PubMed  Google Scholar 

  44. Dansgaard, W. et al. Evidence for general instability of past climate from a 250-kyr ice-core record. Nature 364, 218–220 (1993).

    ADS  Google Scholar 

  45. Bond, G. et al. Correlations between climate records from North Atlantic sediments and Greenland ice. Nature 365, 143–147 (1993).

    ADS  Google Scholar 

  46. Hoff, U., Rasmussen, T. L., Stein, R., Ezat, M. M. & Fahl, K. Sea ice and millennial-scale climate variability in the Nordic seas 90 kyr ago to present. Nat. Commun. 7, 1–10 (2016).

    Google Scholar 

  47. Rasmussen, S. O. et al. A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy. Quat. Sci. Rev. 106, 14–28 (2014).

    ADS  Google Scholar 

  48. Sadatzki, H. et al. Rapid reductions and millennial-scale variability in Nordic Seas sea ice cover during abrupt glacial climate changes. Proc. Natl Acad. Sci. USA 117, 29478–29486 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  49. Sadatzki, H. et al. Sea ice variability in the southern Norwegian Sea during glacial Dansgaard-Oeschger climate cycles. Sci. Adv. 5, eaau6174 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  50. Bory, A. M., Biscaye, P. E., Piotrowski, A. & Steffensen, J. P. Regional variability of ice core dust composition and provenance in Greenland. Geochem. Geophys., Geosyst. 4, 1107 (2003).

  51. Svensson, A., Biscaye, P. E. & Grousset, F. E. Characterization of late glacial continental dust in the Greenland Ice Core Project ice core. J. Geophys. Res.: Atmos. 105, 4637–4656 (2000).

    ADS  CAS  Google Scholar 

  52. Jennings, A. E. et al. Baffin Bay paleoenvironments in the LGM and HS1: resolving the ice-shelf question. Mar. Geol. 402, 5–16 (2018).

    ADS  CAS  Google Scholar 

  53. Müller, J., Massé, G., Stein, R. & Belt, S. T. Variability of sea-ice conditions in the Fram Strait over the past 30,000 years. Nat. Geosci. 2, 772–776 (2009).

    ADS  Google Scholar 

  54. Kolling, H. M., Stein, R., Fahl, K., Perner, K. & Moros, M. Short-term variability in late Holocene sea ice cover on the East Greenland Shelf and its driving mechanisms. Palaeogeogr. Palaeoclimatol. Palaeoecol. 485, 336–350 (2017).

    Google Scholar 

  55. Müller, J. et al. Holocene cooling culminates in sea ice oscillations in Fram Strait. Quat. Sci. Rev. 47, 1–14 (2012).

    ADS  Google Scholar 

  56. Werner, K. et al. Holocene sea subsurface and surface water masses in the Fram Strait–Comparisons of temperature and sea-ice reconstructions. Quat. Sci. Rev. 147, 194–209 (2016).

    ADS  Google Scholar 

  57. Denton, G. H. et al. The last glacial termination. Science 328, 1652–1656 (2010).

    ADS  CAS  PubMed  Google Scholar 

  58. Thompson, W. G. & Goldstein, S. L. A radiometric calibration of the SPECMAP timescale. Quat. Sci. Rev. 25, 3207–3215 (2006).

    ADS  Google Scholar 

  59. Saiz-Lopez, A., Blaszczak-Boxe, C. S. & Carpenter, L. A mechanism for biologically induced iodine emissions from sea ice. Atmos. Chem. Phys. 15, 9731–9746 (2015).

    ADS  CAS  Google Scholar 

  60. Saini, J. et al. Holocene variability in sea ice and primary productivity in the northeastern Baffin Bay. arktos 6, 55–73 (2020).

    Google Scholar 

  61. Gibb, O. T., Steinhauer, S., Fréchette, B., de Vernal, A. & Hillaire-Marcel, C. Diachronous evolution of sea surface conditions in the Labrador Sea and Baffin Bay since the last deglaciation. Holocene 25, 1882–1897 (2015).

    ADS  Google Scholar 

  62. Slingo, A. Sensitivity of the Earth’s radiation budget to changes in low clouds. Nature 343, 49–51 (1990).

    ADS  Google Scholar 

  63. Gordon, H. et al. Causes and importance of new particle formation in the present‐day and preindustrial atmospheres. J. Geophys. Res.: Atmospheres 122, 8739–8760 (2017).

    ADS  Google Scholar 

  64. Jansen, E. et al. Past perspectives on the present era of abrupt Arctic climate change. Nat. Clim. Change 10, 714–721 (2020).

    ADS  Google Scholar 

  65. Kindler, P. et al. Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core. Climate 10, 887–902 (2014).

    Google Scholar 

  66. Overland, J. E. & Wang, M. When will the summer Arctic be nearly sea ice free? Geophys. Res. Lett. 40, 2097–2101 (2013).

    ADS  Google Scholar 

  67. Miller, G. H. et al. Arctic amplification: can the past constrain the future? Quat. Sci. Rev. 29, 1779–1790 (2010).

    ADS  Google Scholar 

  68. Volz, A. & Kley, D. Evaluation of the Montsouris series of ozone measurements made in the nineteenth century. Nature 332, 240 (1988).

    ADS  CAS  Google Scholar 

  69. Prados-Roman, C. et al. Iodine oxide in the global marine boundary layer. Atmos. Chem. Phys. 15, 583–593 (2015).

    ADS  Google Scholar 

  70. Winstrup, M. et al. An automated approach for annual layer counting in ice cores. Clim 8, 1881–1895 (2012).

    ADS  Google Scholar 

  71. Vinther, B. M. et al. A synchronized dating of three Greenland ice cores throughout the Holocene. J. Geophys. Res.: Atmos. 111, https://doi.org/10.1029/2005JD006921 (2006).

  72. Dansgaard, W. & Johnsen, S. J. A Flow model and a time scale for the ice core from camp century, Greenland. J. Glaciol. 8, 215–223 (1969).

    ADS  Google Scholar 

  73. Vallelonga, P. et al. Sea-ice-related halogen enrichment at Law Dome, coastal East Antarctica. Clim. Past. 13, 171–184 (2017).

    Google Scholar 

  74. Andersen, K. K. et al. High-resolution record of Northern Hemisphere climate extending into the last interglacial period. Nature 431, 147–151 (2004).

    ADS  CAS  PubMed  Google Scholar 

  75. Hansson, M. E. The Renland ice core. A Northern Hemisphere record of aerosol composition over 120,000years. Tellus B 46, 390–418 (1994).

    ADS  Google Scholar 

  76. Hughes, A. G. et al. High-frequency climate variability in the Holocene from a coastal-dome ice core in east-central Greenland. Clim 16, 1369–1386 (2020).

    ADS  Google Scholar 

  77. Turekian, K. K. Oceans. 120 (Prentice-Hall, 1968).

  78. Hans Wedepohl, K. The composition of the continental crust. Geochimica et. Cosmochimica Acta 59, 1217–1232 (1995).

    Google Scholar 

  79. Lamarque, J. F. et al. CAM-chem: description and evaluation of interactive atmospheric chemistry in the Community Earth System Model. Geosci. Model Dev. 5, 369–411 (2012).

    ADS  Google Scholar 

  80. Tilmes, S. et al. Representation of the Community Earth System Model (CESM1) CAM4-chem within the Chemistry-Climate Model Initiative (CCMI). Geosci. Model Dev. 9, 1853–1890 (2016).

    ADS  CAS  Google Scholar 

  81. Saiz-Lopez, A. et al. Injection of iodine to the stratosphere. Geophys. Res. Lett. 42, 6852–6859 (2015).

    ADS  CAS  Google Scholar 

  82. Fernandez, R. P., Salawitch, R. J., Kinnison, D. E., Lamarque, J. F. & Saiz-Lopez, A. Bromine partitioning in the tropical tropopause layer: implications for stratospheric injection. Atmos. Chem. Phys. 14, 13391–13410 (2014).

    ADS  Google Scholar 

  83. Ordóñez, C. et al. Bromine and iodine chemistry in a global chemistry-climate model: description and evaluation of very short-lived oceanic sources. Atmos. Chem. Phys. 12, 1423–1447 (2012).

    ADS  Google Scholar 

  84. MacDonald, S. et al. A laboratory characterisation of inorganic iodine emissions from the sea surface: dependence on oceanic variables and parameterisation for global modelling. Atmos. Chem. Phys. 14, 5841–5852 (2014).

    ADS  Google Scholar 

  85. Meinshausen, M. et al. The RCP greenhouse gas concentrations and their extensions from 1765 to 2300. Clim. Change 109, 213–241 (2011).

    ADS  CAS  Google Scholar 

  86. Rienecker, M. M. et al. MERRA: NASA’s modern-era retrospective analysis for research and applications. J. Clim. 24, 3624–3648 (2011).

    ADS  Google Scholar 

  87. Martín, J. C. G. et al. A gas-to-particle conversion mechanism helps to explain atmospheric particle formation through clustering of iodine oxides. Nat. Commun. 11, 1–14 (2020).

    ADS  Google Scholar 

  88. Bates, S. L., Siddall, M. & Waelbroeck, C. Hydrographic variations in deep ocean temperature over the mid-Pleistocene transition. Quat. Sci. Rev. 88, 147–158 (2014).

    ADS  Google Scholar 

  89. Laskar, J. et al. A long-term numerical solution for the insolation quantities of the Earth. Astron. Astrophys. 428, 261–285 (2004).

    ADS  Google Scholar 

Download references

Acknowledgements

This study received funding from the European Research Council (ERC) Executive Agency under the European Union´s Horizon 2020 Research and Innovation programme (Project ‘ERC-2016-COG 726349 CLIMAHAL’). The RECAP ice coring effort was financed by the Danish Research Council through a Sapere Aude grant, the NSF through the Division of Polar Programmes, the Alfred Wegener Institute, and the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement 610055 through the Ice2Ice project and the Early Human Impact project (267696).

Author information

Authors and Affiliations

Authors

Contributions

A.S.-L., A.S. and J.P.C. designed the research. N.M., A.S. and F.S. prepared the ice core samples and ran the chemical analysis. C.A.C. conducted the CAM-Chem model simulations. A.S., N.M., F.S., G.C., P.V., C.B., B.V., H.A.K., J.M. and R.E. helped with data interpretation. J.P.C. and A.S.-L. wrote the manuscript with contributions from all authors.

Corresponding authors

Correspondence to Juan Pablo Corella or Alfonso Saiz-Lopez.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review information

Nature Communications thanks Melinda Nicewonger and the other anonymous reviewers for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Corella, J.P., Maffezzoli, N., Spolaor, A. et al. Climate changes modulated the history of Arctic iodine during the Last Glacial Cycle. Nat Commun 13, 88 (2022). https://doi.org/10.1038/s41467-021-27642-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1038/s41467-021-27642-5

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing