text:
Deep-water circulation changes lead North Atlantic climate during deglaciation
Latest Pliocene Northern Hemisphere glaciation amplified by intensified Atlantic meridional overturning circulation
Tatsuya Hayashi, Toshiro Yamanaka, … Masao Ohno
Enhanced climate instability in the North Atlantic and southern Europe during the Last Interglacial
P. C. Tzedakis, R. N. Drysdale, … G. Zanchetta
Abrupt ice-age shifts in southern westerly winds and Antarctic climate forced from the north
Christo Buizert, Michael Sigl, … Eric J. Steig
Opposite response modes of NADW dynamics to obliquity forcing during the late Paleogene
Hojun Lee, Kyoung-nam Jo & Sangmin Hyun
Natural drivers of multidecadal Arctic sea ice variability over the last millennium
Paul R. Halloran, Ian R. Hall, … Freya Garry
Abrupt high-latitude climate events and decoupled seasonal trends during the Eemian
J. Sakari Salonen, Karin F. Helmens, … Miska Luoto
Split westerlies over Europe in the early Little Ice Age
Hsun-Ming Hu, Chuan-Chou Shen, … Robert Korty
Atlantic deep water circulation during the last interglacial
Yiming Luo, Jerry Tjiputra, … Jörg Lippold
Observed interannual changes beneath Filchner-Ronne Ice Shelf linked to large-scale atmospheric circulation
Tore Hattermann, Keith W. Nicholls, … Torsten Kanzow
Francesco Muschitiello1,2,3,
William J. D'Andrea2,
Andreas Schmittner4,
Timothy J. Heaton5,
Nicholas L. Balascio6,
Nicole deRoberts2,
Marc W. Caffee ORCID: orcid.org/0000-0002-6846-89677,8,
Thomas E. Woodruff7,
Kees C. Welten9,
Luke C. Skinner10,
Margit H. Simon3 &
Trond M. Dokken3
Palaeoceanography
Palaeoclimate
Constraining the response time of the climate system to changes in North Atlantic Deep Water (NADW) formation is fundamental to improving climate and Atlantic Meridional Overturning Circulation predictability. Here we report a new synchronization of terrestrial, marine, and ice-core records, which allows the first quantitative determination of the response time of North Atlantic climate to changes in high-latitude NADW formation rate during the last deglaciation. Using a continuous record of deep water ventilation from the Nordic Seas, we identify a ∼400-year lead of changes in high-latitude NADW formation ahead of abrupt climate changes recorded in Greenland ice cores at the onset and end of the Younger Dryas stadial, which likely occurred in response to gradual changes in temperature- and wind-driven freshwater transport. We suggest that variations in Nordic Seas deep-water circulation are precursors to abrupt climate changes and that future model studies should address this phasing.
Precise reconstructions that resolve the relative timing of changes in North Atlantic Ocean circulation, climate, and carbon cycling are necessary to anticipate the mechanisms initiating and propagating abrupt global climate changes. During the last deglaciation (∼18,000-11,000-years ago), the climate system underwent numerous abrupt changes that have been attributed to variations in the strength of the Atlantic Meridional Overturning Circulation (AMOC)1,2. Through changes in high-latitude North Atlantic Deep Water (NADW) formation and export, AMOC exerts an important control on the global climate system by redistributing heat near the surface and regulating carbon storage at depth. In particular, the partitioning of carbon between the surface and deep ocean is thought to play a critical role in centennial-to-millennial-scale variations of atmospheric CO2 (refs 3,4,5). However, reconciling the deglacial history of changes in overturning circulation as recorded in marine records with North Atlantic climate and pCO2 as inferred by Greenlandic and Antarctic ice cores, respectively, remains challenging. First, highly resolved records from deep convection sites sensitive to NADW that monitor the descending branch of AMOC are still lacking. Secondly, large uncertainties in high-latitude marine reservoir ages6 limit the precision of marine 14C-based chronologies. Thirdly, direct alignment of marine records to far afield Greenland ice-core stratigraphies hinders testing hypotheses of synchronicity. Lastly, precise comparisons between marine and ice-core climate records are hampered by inconsistencies between the radiocarbon and ice-core timescales7,8.
Here we present a new synchronization of high-latitude NADW, climate, and pCO2 records for the last deglaciation based on new marine and ice core data that allows us to conclude for the first time that changes in deep-water circulation in the Nordic Seas led rapid shifts in North Atlantic climate and changes in carbon cycling.
Site location, 14C ventilation and chronology
We generated a continuous record of deep/intermediate- and surface-water 14C ventilation age from 14C measurements on planktic and benthic foraminifera (Methods) in sediment core MD99-2284 (62° 22.48 N, 0° 58.81 W, 1500 m water depth) from the Norwegian Sea (Fig. 1).
Site location. a Location of core MD99-2284 (star), NGRIP ice cores (black circle)101, and other sediment cores (yellow circles) and deep-sea corals (yellow triangles) discussed in this study. The white areas indicate the extent of Northern Hemisphere ice sheets at 12500 years BP102. Red and white arrows show warm Atlantic inflow to the Nordic Seas and main bottom current pathways in the northern North Atlantic. Simulated winter (March, purple) and summer (September, orange) 50% sea-ice cover fraction during GS-1 (12,500 years BP)103 is also shown. b Meridional section of radiocarbon concentration within the Atlantic Ocean (averaged over 0–40°W) from a pre-industrial control simulation using a coupled climate-biogeochemical model104
Site MD99-2284, which is characterised by exceptionally high sedimentation rates (>400 cm kyr-1), is located at the gateway of the Faroe-Shetland Channel (FSC), where warm surface Atlantic water flows into the Nordic Seas and cold dense water overflows into the North Atlantic. Critically, this overflow water is one of two main NADW pathways flowing into the deep North Atlantic and a key constituent of the AMOC9. During the last glacial period and deglaciation, overflow through the FSC remained a continuous source of NADW10,11. Hence, because deep-water 14C activity reflects the circulation-driven exchange of carbon between the atmosphere and deep-ocean reservoir, bottom-surface water 14C age differences of the FSC directly inform past changes in Nordic Seas deep convection, NADW production, and its southward export12,13.
The age model for the core was established using a combination of tephrochronology and alignment between sea-surface temperature records from core MD99-2284 and a high-resolution hydroclimate reconstruction from a relatively closely located terrestrial sequence in southern Scandinavia (Methods; Supplementary Figs. 1–5 and Note 1–2). The approach enables us to precisely place our marine proxies on the IntCal13 timescale14 and to use the foraminiferal radiocarbon data (Methods) to calculate the marine 14C ventilation age. Our estimate was determined using a random walk model (RWM) (Methods and Supplementary Methods) fitted via Markov chain Monte Carlo (MCMC) that took into account uncertainty structures in both calendar age modelling and 14C measurements. To allow detailed comparison with Greenlandic and Antarctic ice-core records, we synchronized the ice-core GICC05 (ref. 15) and WD2014 (ref. 16), and 14C timescales using previously published and new 10Be records from GRIP17 and WAIS Divide ice cores, respectively (Methods; Supplementary Fig. 6). Ages are hereafter reported as IntCal13 years before 1950 AD ± 1σ (BP).
Deglacial ventilation history of the deep Nordic Seas
Surface and bottom water mass reconstructions at our site are consistent with existing paleoceanographic records of water properties, transport and exchange between the Norwegian Sea and the northern North Atlantic (Supplementary Figs. 7–8), indicating that our reconstructions are representative of regional oceanographic conditions. Benthic-planktic (B-P) ventilation ages in MD99-2284 decreased by ∼700 years seemingly shortly preceding the abrupt warming transition from Greenland Stadial (GS) 2 (equivalent to Heinrich Stadial 1, HS-1, and the Last Glacial Maximum) into Greenland Interstadial (GI) 1 (equivalent to the Bølling-Allerød interstadial, BA; 14581 ± 16 years BP) (Fig. 2a–d; Supplementary Fig. 9). Although relatively older deep/intermediate water during GS-2 (HS-1) at our site is inferred from only one B-P estimate, this is confirmed by other regional B-P records (Supplementary Fig. 9). B-P ages gradually decreased during GI-1 (BA) by an additional ∼700 years, reaching or surpassing present-day values (0–50 years) as late as 14,150 years BP (Fig. 2d). This trend suggests that strengthening of deep convection across the onset of GI-1 (BA) was slower than previously thought, which can be attributed to a gradual northward re-initiation of NADW production in the North Atlantic. This interpretation is supported by transient simulations of the last deglaciation, which reveal a northward time-transgressive recovery of NADW formation during GI-1 (BA), with a ~400-year long re-initiation of deep convection in the Nordic Seas relative to convection sites south of Greenland18. The century-scale northward order of NADW re-initiation is consistent with other model results19. It is also in line with paleoceanographic reconstructions suggesting a delayed resumption of deep convection in the high-latitude Nordic Seas starting midway through GI-1 (BA)20.
Synchronized Antarctic and North Atlantic climate records on the IntCal13 timescale. a CO2 concentration from WAIS Divide ice cores (WDC)3. b δ18O values from NGRIP ice cores101. c Bottom and surface 14C ventilation (R) histories from core MD99-2284 based on benthic and planktic foraminifera, respectively. Dots indicate individual measurements together with their 2σ error bars. Shading reflects the 95% posterior credible interval (2σ) of 14C ventilation as a function of age found by MCMC using random walk model (Methods) and taking into account both the analytical and chronological uncertainty in our observed data. Orange dot reflects individual early GS-1 measurement from nearby core JM11-FI-19PC (1179 m)12 (not incorporated in the random walk model). Red and blue arrows indicate the modern surface (0–100 m) and bottom (1200–1500 m) R values at the study site104, respectively. Grey arrow indicates modern mean marine R. d Benthic-planktic (B-P) offset based on ventilation estimates in (c) and reflecting the strength of deep convection and NADW formation in the Nordic Seas. Red line and shading denote the posterior median value and pointwise 95% credible intervals. Purple square indicates a B-P estimate from the Vøring Plateau (1048 m) based on one solitary U/Th dated deep-sea coral105 (average of three B-Atm measurements—Supplementary Fig. 9) and presented on its independent time scale. e Second principal component of all foraminifera counts in core MD99-2284 dominated by N. labradoricum, which is an indicator for the proximity of the Polar Front and sea-ice edge106 (Supplementary Fig. 10). f Coarse-grained (>150 μm) ice-rafted debris count and PIP25 index in core MD99-2284 reflecting iceberg rafting and occurrence of spring sea-ice cover, respectively. g Abundance of long-chain n-alkanes in core MD99-2284 derived from terrestrial higher plants107 and indicating meltwater discharge. All records are presented on the IntCal13 timescale14. Greenland stratigraphic events relative to the IntCal13 timescale are displayed at the top and cold events are highlighted with grey bars (GS: Greenland Stadial; GI: Greenland Interstadial; YD: Younger Dryas Stadial; BA: Bølling-Allerød Interstadial; HS1: Heinrich Stadial 1). Dashed vertical lines show tephra horizons identified in core MD99-2284
Prior to the onset of GS-1 (12870 ± 26 years BP), which corresponds to the Younger Dryas stadial (YD), the B-P age offset rapidly increased by ∼900 years starting at 13,250 years BP. Ahead of the termination of GS-1 (11577 ± 16 years BP), the B-P offset returned to modern values beginning at 11,890 years BP (Figs. 2d–3). GS-1 (YD) onset and termination are also preceded by surface oceanographic changes inferred from sedimentological and lipid biomarker signatures in MD99-2284. Specifically, increasing (decreasing) B-P ages are virtually synchronous with a southward (northward) migration of the polar front across the coring site inferred from faunal assemblages (Fig. 2e; Supplementary Fig. 10), relatively more (less) frequent sea-ice edge conditions inferred from PIP25 (Methods) and iceberg rafting inferred from IRD (Fig. 2f), and relatively higher (lower) inputs of terrestrial meltwater inferred from sedimentary n-alkane concentrations (Methods; Fig. 2g), ultimately pointing at a tight link between freshwater dynamics and bottom-water ventilation. Although the benthic-atmosphere (B-Atm) offset dominates the B-P signal, the P-Atm offset also decreased prior to the onset of GS-1 (YD) and in phase with the rise in B-P values (Fig. 2c, d). Specifically, P-Atm values slowly declined by ∼500 years starting ∼13,250 years BP and approached the age of the contemporaneous atmosphere by ∼12900 years BP, suggesting that the climatic transition into GS-1 (YD) was preceded by a major slowdown of the Atlantic Inflow that resulted in enhanced isotopic equilibration between surface waters and the atmosphere21.
North Atlantic deep/intermediate water 14C ventilation during deglaciation. a High-pass filtered (cut-off frequency 1/200 years) atmospheric radiocarbon content (Δ14C) corrected for changes in natural 14C production27, reflecting changes in global ocean ventilation and NADW formation28,29. b Comparison of bottom 14C ventilation histories (B-Atm) from core MD99-2284 (grey dots and dashed lines) and B-Atm records from the Faroe-Shetland Channel (1179 m)12 (red) and south of the Iceland-Scotland Ridge (1237 m)12 (orange) using benthic foraminifera (excluding Pyrgo spp. and Miliolids) (note the broken y-axis and scale change at ∼2500 14C years). c Same as in (b) using intermediate coral 14C records from Orphan Knoll108 and the New England Seamounts22,23,24,25 and deep/intermediate coral 14C records from the New England Seamounts22,23,24,25 in the northwest Atlantic Ocean (Fig. 1). d Same as in (b) using deep/intermediate coral 14C records from the Equatorial Atlantic26. Arrows show modern depth-averaged ventilation ages at each location. Histograms display depth represented by individual data points. Note that the 14C record from ref. 24 used here includes one measurement from the northeast Atlantic. Bars reflect 2σ errors of individual measurements. All records are presented on their independent time scale. Greenland stratigraphic events relative to the IntCal13 timescale are displayed at the top (GS: Greenland Stadial; GI: Greenland Interstadial; YD: Younger Dryas Stadial; BA: Bølling-Allerød Interstadial; HS1: Heinrich Stadial 1)
The timing and magnitude of changes in our B-Atm record are supported by other independent reconstructions from the FSC and the Iceland Basin, i.e., regions that monitor ventilation of Nordic Seas overflow waters12 (Fig. 3). The data are also in agreement with deep ventilation records from U–Th dated corals in the deep/intermediate northwest and equatorial Atlantic, which monitor the strength of downstream NADW transport22,23,24,25,26. Specifically, the slow decrease in B-Atm ages at the end of GS-2 (HS-1) and during the early GI-1 (BA) is observed in data from both the northeast and equatorial Atlantic, whereas early changes prior to the start and end of GS-1 (YD) are evident in each of the records. We hypothesise that these parallel changes in radiocarbon across sites likely reflect downstream propagation of North Atlantic deep/intermediate water 14C signatures associated with high-latitude NADW ventilation rates. In addition, the timing and pattern of the reconstructed NADW changes are consistent with shifts in production-corrected atmospheric Δ14C (ref. 27) (Figs. 3, 4), which primarily reflects changes in ocean ventilation and high-latitude NADW formation28,29, and thus support an oceanic origin for the variations in atmospheric radiocarbon concentrations during the last deglaciation. Overall, the agreement across the range of independent deep ventilation records supports the view that our B-Atm data capture large-scale fluctuations in NADW circulation. A strong coupling between AMOC and deep ventilation in the Nordic Seas is also corroborated by transient climate model experiments (Supplementary Fig. 11 and Note 3). The simulations, which reproduce a shutdown and subsequent recovery of the AMOC, identify and validate sign and magnitude of the shifts observed in our B-P age reconstruction.
Detailed view of synchronized CO2, climate and ocean circulation records during the Younger Dryas stadial. a CO2 concentration from WDC ice cores. b Greenland temperature reconstruction109 (average of NEEM, GISP2 and NGRIP ice cores). c Nordic Sea deep water formation rate based on B-P ages from core MD99-2284 (this study). Triangles indicate the location of individual planktic (blue) and (benthic) 14C measurements. d Atmospheric radiocarbon content (Δ14C) corrected for changes in natural 14C production27. The record was high-pass filtered (cut-off frequency 1/200 years) to facilitate comparison with long-term overturning strength variations (c). All the records were generated by MCMC using the same random walk model approach used to estimate 14C ventilation, but the model was fitted directly to the observations. Grey envelopes reflect the 68% credible intervals (1σ) associated with both analytical and chronological uncertainty in the raw data. Coloured segments represent the 68% confidence intervals of the fitted piecewise linear regression functions accounting for both analytical and chronological uncertainties. The regression models were estimated using a modified version of an algorithm for breakpoint analysis100. Note that only the sharp decline in Δ14C during the second half of GS-1 (YD) can be attributed to a resumption of NADW formation27, whereas the preceding minor decline is likely a dynamical response in the Southern Ocean associated with venting of 14C-depleted carbon110. e 95% confidence intervals of the timing and duration of the transitions (Supplementary Table 1). Inset black lines show the 2σ uncertainty of the estimated start and end of each transition. Greenland stratigraphic events relative to the IntCal13 timescale are displayed at the top (GS: Greenland Stadial; GI: Greenland Interstadial; YD: Younger Dryas Stadial; BA: Bølling-Allerød Interstadial)
A detailed comparison between our B-P record and ice-core-based temperature records from Greenland using breakpoint analysis (Methods; Supplementary Table 1; Fig. 4), reveals a complex relationship between Nordic Seas ventilation and abrupt climate change. The analysis reveals that changes in Nordic Seas NADW formation occurred before the climate shifts into and out of GS-1 (YD) by 385 ± 32 (1σ bounds) and 447 ± 27 years, respectively, and that weakening of Nordic Seas NADW occurred 437 ± 79 years prior to the first signs of pCO2 rise near the start of GS-1 (YD). Importantly, the latter finding substantiates the hypothesised30,31 lag between AMOC reduction and pCO2 rise during early deglaciation and is observed in climate simulations as a transient response of the global efficiency of the biological pump to AMOC slowdown (Supplementary Figs. 12–14 and Note 3).
The lag between our high-latitude NADW records and Greenland temperatures across the transitions into and out of GS-1 (YD) requires further consideration. Cooling of sea-surface temperatures and sea-ice expansion in the Norwegian Sea preceding GS-1 (YD) have been documented in terrestrial temperature and isotope records32,33, and are consistent with increasing B-P offsets observed in MD99-2284 beginning at 13,250 years BP. Beyond the last deglaciation, records from the last glacial and Late Pleistocene34,35 have also implied increased input of ice-sheet meltwater, gradual surface cooling, and southward migration of the polar front in the Nordic Seas preceding transitions from warm interstadials to cold stadial conditions, analogous to the oceanographic changes inferred from MD99-2284 that preceded GS-1 (YD).
On the other hand, ocean warming in the Nordic Seas, a northward diversion of the polar front, and NADW resumption during the second half of the GS-1 (YD) have been widely reported36,37,38, and further support decreasing B-P offsets in MD99-2284 starting at 11,890 years BP. A lead of a few centuries of changes in high-latitude NADW formation ahead of Greenland temperatures is compatible with the timing of AMOC slowdown and recovery during deglaciation inferred from several Pa/Th circulation39 and other AMOC-sensitive proxies40,41. Despite uncertainties with the 14C chronologies and coarse temporal resolution, it has been demonstrated that a century-scale response of Pa/Th to AMOC changes should be additionally accounted for when interpreting this proxy as a function of changes in ocean circulation42,43. The existence of a significant lag time for Greenland temperature and Pa/Th records behind changes in NADW would represent an important observational constraint to account for, for example when attempting to infer the physical mechanisms of rapid climate and AMOC change from climate model simulations44,45. Equally abrupt and synchronous changes in the AMOC, Pa/Th and Greenland temperature (i.e., to within a few centuries) should not, therefore, be expected in numerical model simulations.
On the basis of timing inferred from our 14C ventilation records, we propose that the lagged Greenland temperature response behind changes in Nordic Seas NADW formation can be understood as a threshold response to gradual changes in surface temperature-driven freshwater transport, which ultimately modulates the strength of the northward oceanic heat transport. Both empirical34,35,46 and climate modelling studies47 have proposed a nonlinear salt oscillator in the North Atlantic system, whereby meltwater production rates are greater (smaller) toward the end of warm (cold) episodes. This is consistent with evidence for southward (northward) migration of the polar front, more (less) extensive sea-ice cover, and increased (decreased) iceberg discharge at our coring site shortly before the termination of GI-1 (GS-1). We further hypothesise that shifts to higher (lower) rates of meltwater fluxes, although initially prompted by gradual climate warming (cooling), are amplified by concomitant changes in surface winds, which play a critical role in the rearrangement of water masses in the Nordic Atlantic. For instance, model studies suggest that sustained cold events like GS-1 (YD)48 and the Little Ice Age49,50 can be triggered by a positive feedback loop, whereby increasing export of meltwater to the subpolar gyre in response to warming leads to surface ocean cooling and consequent strengthening of atmospheric blocking over the North Atlantic–i.e. large-scale quasi-stationary anticyclonic circulation. This feedback can impede northward heat transport and suppress NADW formation for centuries. Moreover, reduced NADW can lead to an additional freshening through basin-wide subsurface warming51,52, which causes ice-shelf thinning and iceberg discharge52, further bolstering NADW formation weakening. However, we cannot rule out the potential amplifying role of a catastrophic drainage of freshwater from the Baltic Ice Lake into the Nordic Seas53, which occurred precisely at the start of GS-1 (YD), and could have further impacted regional NADW formation.
Atmospheric blocking could be equally important as a mechanism for NADW resumption in the Nordic Seas at the transition out of GS-1 (YD). New studies suggest that under full GS-1 (YD) (ref. 54) and glacial conditions55, the cold Fennoscandian Ice Sheet induces a strong southwesterly wind flow over the Norwegian Sea associated with blocking circulation over the ice sheet. This configuration is more conducive to reduced sea-ice cover through export of sea ice out of the Norwegian Sea and increased surface warming and salinity in response to northward oceanic transport of heat and salt, which in turn gradually promotes deep-water formation in the Norwegian Sea55. Our data and other existing records all lend support to the occurrence of this proposed mechanism midway through GS-1 (YD) (i.e., stronger southwesterly winds36,37, warmer surface ocean conditions56,57, and gradually less frequent sea-ice occurrence58,59) (Supplementary Fig. 7).
In conclusion, our results suggest that gradual changes in high-latitude NADW formation are precursors of rapid climate shifts in the North Atlantic and emphasize the central role of ocean circulation in abrupt climate change, as well as its sensitivity to atmosphere and cryosphere dynamics. Given recent evidence that the current AMOC has been slowing down for several decades60,61, our findings broach the question as to whether the current decline in deep-water circulation may herald a new phase of abrupt change.
The chronology of core MD99-2284 was established by aligning variations in downcore sea-surface temperatures (SST)36,62 with synchronous changes in hydroclimate as recorded in sedimentary leaf-wax hydrogen isotope (δD) records from the ancient lake of Atteköp, southern Sweden63 (Supplementary Figs. 1–5 and Note 1). The terrestrial site is located closely downwind of our marine core location in an area where most of the precipitation is sourced from the North and Norwegian Seas64, and where hydroclimate shifts have a demonstrated modern and past linear relationship with upwind, near-field upper ocean temperatures via changes in sea-to-air-moisture fluxes (Supplementary Figs. 1–2). Leaf-wax δD records in northern Europe have been successfully used as sensitive indicators of annual changes in the isotopic composition of the marine precipitation source and moisture availability33,65. Therefore, precise synchronization of marine SST and terrestrial δD records in this region affords an accurate chronology for core MD99-2284 that circumvents alignment to the far-afield Greenland ice-core isotope stratigraphy, which is not necessarily representative of near-surface paleoceanographic conditions in the North and Norwegian Seas33.
Synchronization between the marine and terrestrial proxy time series was obtained using an automated stratigraphic alignment algorithm driven by a Markov chain Monte Carlo method33,66 run for 106 iterations. The approach involves nonlinear deformation of the SST time series onto the reference δD time series to deliver an optimal alignment accounting for uneven compaction and/or expansion of sediments over time, as well as for analytical errors associated with the proxy data (Supplementary Fig. 3). An account of the mathematical formulation associated with the algorithm is presented in ref. 67.
Atteköp's chronology was constructed using a Bayesian age model based on the Hässeldalen Tephra and 37 AMS 14C dates from selected terrestrial plant macrofossils63 (Supplementary Fig. 3) calibrated with the IntCal13 curve14, which in turn allows placing MD99-2284 records on the atmospheric IntCal13 timescale. However, the alignment between the marine and terrestrial records was done only up to ∼12,600 years BP due to the scarcity of datable plant macrofossils in Atteköp in this portion of the core (Supplementary Fig. 3). To secure the chronology upcore, we thus employed three well-dated tephra layers identified in MD99-2284 sediments–the Vedde Ash, the Abernethy Tephra, and the Saksunarvatn Ash (Supplementary Note 2)—and fit the tephra-based age constraints using a Gaussian Regression model68.
Tephra analysis
Two important regional isochrones, i.e. the Vedde Ash and the Saksunarvatn Ash (both present in Greenland ice cores69), have been previously reported in core MD99-2284 (ref. 36). These horizons have been precisely and accurately radiocarbon dated in Lake Kråkenes, Norway70. The markers are here supplemented by the first reported marine occurrence of the Abernethy Tephra71,72 (Supplementary Fig. 4). The finding constitutes the outcome of a targeted search for late YD tephra based on the counting of glass shards within the interval 200–360 cm in the marine core. A peak in cryptotephra shards was identified in samples 247.5 cm and 249.5 cm. The shards were isolated from the sediment following the methods of ref. 73. Samples were treated with 10% H2O2 and heated to remove organic material, then washed with deionized water over a 63-µm sieve. A series of heavy liquid density separations using sodium polytungstate were then performed to isolate material between 2.2 and 2.5 g cm−3, which was mounted on 27 × 46 mm glass slides in epoxy resin. Slides were polished to expose grain interiors and analyzed at the Concord University Microanalytical Laboratory (WV, USA) using an ARL SEMQ electron microprobe equipped with six wavelength-dispersive spectrometers and a Bruker 5030 SDD energy-dispersive spectrometer. Instrument conditions included a 14 kV accelerating voltage, 10 nA beam current, and beam size of 4 to 6 µm. Results were reported as non-normalized major oxide concentrations (Supplementary Data 1).
Results from the analysis of tephra grains in samples 247.5 cm and 249.5 cm show a single geochemical population of rhyolitic composition similar to the widespread Vedde Ash erupted from the Katla volcano in Iceland during the middle of the YD. Recently, it has been shown that there was a second tephra-producing eruption of Katla within the YD, but closer to the YD-Holocene boundary named the Abernethy Tephra72. Tephra layers from Loch Etteridge (LET-5)72 and Abernethy Forest (AF555)71, in Scotland, were attributed to this eruption, which are both similar in stratigraphy and geochemistry to tephra that we have identified (Supplementary Fig. 4). Based on precise varve counting72 and other lines of evidence74,75, it has been argued75 that the Abernethy Tephra is likely 320 ± 20 years younger than the Vedde Ash (12,064 ± 48 years BP)70. Therefore we here assigned to the layer an age of 11,744 ± 50 years BP. Radiocarbon dating of this layer based on one terrestrial plant macrofossil71 from a sediment sequence in Abernethy Forest yields a calibrated age range of 12,380–11770 years BP, which is consistent with, or slightly underestimates our age assignment.
14C dating
Well-preserved monospecific shells of the planktic foraminifer Neogloboquadrina pachyderma (s) (calcification depth ∼30 to 200 m (ref. 76)), two samples of the benthic foraminifer Cibicidoides wuellerstorfi, one shell of Pyrgo murrhina, and a number of mixed benthic foraminifera were hand-picked from slices of core MD99-2284 for AMS 14C dating for a total of 88 measurements (41 planktic and 47 benthic), including 8 previously reported dates36 (Supplementary Data 2). Mixed benthic samples are mainly composed by C. neoteretis excluding Pyrgo spp and miliolid spp. Samples were mainly selected from the >150 μm fraction, except for 31 samples, which were picked from the 63-150 μm fraction. All shells were cleaned and dried prior to preparation for accelerator mass spectrometry (AMS) measurements, whereas large samples were additionally rinsed with dilute, organic-free HCl.
Large samples were converted to graphite and analysed using standard AMS 14C measurement procedures within the laboratories of Beta Analytic Miami (FL), USA, and ETH Zurich, Switzerland. All the small fraction samples were processed and analysed using a newly developed method77. The approach involves direct analysis of CO2 from ultra-small amounts of carbonate (∼0.5 mg) using a compact AMS facility equipped with a gas ion source at the Laboratory of Ion Beam Physics, ETH Zurich. To test for contamination by secondary carbonates on the small fraction samples, we performed leaching experiments on the sample material surface using HCl 0.02 M and following the procedure outlined in ref. 78. All radiocarbon ages are here reported according to the standard protocol of ref. 79. A subset of 7 samples were prepared at the University of Cambridge using methods detailed in ref. 80 and subsequently dated by AMS at the 14 Chrono Centre, University of Belfast.
The resulting 14C ages served as a basis to generate subsurface (∼100 m) and bottom water (1500 m) ventilation records employing the random walk model described below, and ultimately to reconstruct benthic-planktic (B-P) ventilation ages81. We identified seven outliers that were dismissed from our bottom water ventilation record (one in GS-1 at 400.5 cm; three in GI-1 at 415, 445.5, 534.5 cm; three in GS-2 at 576.5, 581.5, 601 cm). These were all ages based on benthic foraminifera, including one Pyrgo spp. age, which has been shown to yield old 14C ages12. The other benthic dates, which are all based on small shell fragments <1 mg in weight and all associated with the 63-150 μm fraction, were either significantly younger than the corresponding planktic ages (where available), or significantly younger than the contemporaneous atmosphere. Acid leaching tests indicate that four out of the seven outliers are likely affected by modern CO2 carbon contamination, as evidenced by the high 14C content in the respective leach fractions (Supplementary Data 2). Given that sedimentation rates at our site are remarkably high (>400 cm kyr−1), we argue that complications associated with bioturbation and differential mixing effects are likely negligible82,83. Rather, we argue that the systematically younger ages measured in these samples may be the result of preferential modern-carbon contamination among the finer carbonate fraction through secondary growth of calcite within the sediment pore waters and exchange between calcite and the carbonate ions in pore waters83.
Random walk model
Regional reservoir (ΔR) and reservoir (B-Atm and P-Atm) estimates were inferred using a random walk model that incorporates both the uncertainty in the calendar age modelling and our 14C measurements (Supplementary Data 3). B-P estimates were obtained by subtracting the modelled B-Atm and P-Atm curves. This approach was preferred to the projection method, which is more influenced by atmospheric Δ14C variations in the North Atlantic84. In the following, we provide a brief overview of the random walk model and the construction of the reservoir age curves.
For each of our N ocean objects, we observe an estimate ti of its true calendar age θi; and a radiocarbon determination zi that is also subject to noise. This radiocarbon determination can be decomposed into the site-general marine radiocarbon age m(θi); and \({\mathrm{\Delta }}R\left( {\theta _{\mathrm{i}}} \right)\) the local reservoir variation i.e., difference between the reservoir age in the Norwegian Sea. We aim to estimate the Norwegian Sea \({\mathrm{\Delta }}R\left( {\theta _{\mathrm{i}}} \right)\) given these paired observations \((t_{\mathrm{i}},z_{\mathrm{i}})_{{\mathrm{i}} = 1, \ldots ,N}\) and our prior knowledge about the value of \(m\left( {\theta _{\mathrm{i}}} \right)\) provided by the Marine13 radiocarbon calibration curve (ref. 14). We then infer the value of the absolute Norwegian Sea reservoir correction \(R\left( \theta \right)\) by comparison with the atmospheric radiocarbon curve at the same \(\theta\). Specifically we model
$$t_{\mathrm{i}} = \theta _{\mathrm{i}} + \epsilon _{\mathrm{i}}$$
$$z_{\mathrm{i}} = m\left( {\theta _{\mathrm{i}}} \right) + {\mathrm{\Delta }}R\left( {\theta _{\mathrm{i}}} \right) + \eta _{\mathrm{i}},$$
where \(\epsilon _{\mathrm{i}}\) and \(\eta _{\mathrm{i}}\) are the uncertainties in our measurement of calendar age and radiocarbon age assumed here to be independent and identically distributed (IID) with mean 0 and variances \(\tau _{\mathrm{i}}^2\) and \(\sigma _{\mathrm{i}}^2\), respectively. Note in particular that the true calendar ages \(\theta _{\mathrm{i}}\) are only observed subject to noise. This should be incorporated into any estimation procedure. If we ignore this additional uncertainty we are likely to introduce a bias in our estimate for \({\mathrm{\Delta }}R\left( {\theta _{\mathrm{i}}} \right)\) and oversmooth.
We take a Bayesian approach for our estimation of \({\mathrm{\Delta }}R\left( {\theta _{\mathrm{i}}} \right)\). This allows us to incorporate our prior beliefs about the value of \(m\left( {\theta _{\mathrm{i}}} \right)\); model the evolution of \({\mathrm{\Delta }}R\left( {\theta _{\mathrm{i}}} \right)\) over time in a physically interpretable way; and update these prior beliefs in a consistent framework. To make our problem identifiable we require strong prior information on the site-general \(m(\theta )\). This is provided by Marine13, the marine radiocarbon calibration curve14. At any calendar age \(\theta\) this provides pointwise estimates of the value
$$m\left( \theta \right) \sim N\left( {\mu \left( \theta \right),\sigma _{\mathrm{c}}^2\left( \theta \right)} \right).$$
Additionally, we require a prior model for the evolution of the Norwegian Sea's reservoir age \({\mathrm{\Delta }}R\left( \theta \right)\). While we do not expect this to be constant over time, we would expect that knowing \({\mathrm{\Delta }}R\left( \theta \right)\) at a particular \(\theta\) provides some information about its likely value at nearby times. We, therefore, place a Wiener process (random walk) prior on \({\mathrm{\Delta }}R\left( \theta \right)\) similar to that used in the modelling of the IntCal atmospheric calibration curve (see refs. 85,86,87) whereby
$${\mathrm{\Delta }}R\left( \theta \right)|{\mathrm{\Delta }}R\left( {\theta ^ \ast } \right) \sim N\left( {{\mathrm{\Delta }}R\left( {\theta ^ \ast } \right),\rho \left| {\theta - \theta ^ \ast } \right|} \right),$$
The above random walk model is then combined with our prior on \(m(\theta )\) and our observed \((t_{\mathrm{i}},z_{\mathrm{i}})_{{\mathrm{i}} = 1, \ldots ,N}\) through a Metropolis-within-Gibbs sampler. This provides a posterior estimate for \({\mathrm{\Delta }}R\left( \theta \right)\) that allows us to borrow strength from the surrounding observations (i.e., neighbouring \(\theta\)'s). It also updates our site-general marine \(m\left( \theta \right)\). Combining both the estimate of \({\mathrm{\Delta }}R\left( \theta \right)\) and the updated \(m\left( \theta \right)\) and comparing with the atmospheric radiocarbon levels provided by IntCal13 (ref. 14) then allows us to estimate the absolute Norwegian Sea reservoir correction \(R\left( \theta \right)\), i.e., B-Atm and P-Atm. Details on the implementation of this sampler can be found in Supplementary Methods.
10Be measurements
A total of 170 samples from 1793 to 2279 m depth in the WAIS Divide 06A ice core (WDC-06A) were analysed for 10Be concentrations at UC Berkeley and at the Purdue University PRIME Lab (Supplementary Data 4). Samples typically represent continuous ice core sections of ~3 m length (although they varied from 1.9 to 3.9 m) corresponding to ~30 years of snow accumulation. Ice samples of 300–700 g were weighed, melted, and acidified with a solution containing ~0.18 mg Be carrier. The samples were passed through 30 micron Millipore filter, and loaded on cation exchange columns from which the Be fraction was eluted following established procedures88,89. The 10Be/9Be ratios of samples and blanks were measured by AMS at PRIME Lab. Absolute ratios were obtained by normalizing the measurements to well-documented 10Be/9Be standards90. Results were corrected for an average blank 10Be/9Be ratio of 15 ± 3, which corresponds to typical blank corrections of 1–3% of the measured 10Be/9Be ratios. The 10Be results up to 12,000 years BP were used by ref. 16 for the synchronization of the WDC core with the Intcal13 14C record.
Timescale synchronization
Solar activity cycles modulate the amount of Earth bound cosmic rays that controls the production rates of cosmogenic radionuclides 10Be and 14C in the upper atmosphere. The former are integrated in ice cores; the latter are incorporated in a number of absolutely dated records (e.g., tree rings) that take up carbon directly from the atmosphere at the time of formation. Hence, matching of the globally synchronous common short-term 10Be and 14C variations as reconstructed in ice cores and IntCal13 records14, respectively, provides a unique tool for estimating timescale differences in a nearly continuous fashion7,8.
Here we aligned relative changes in GRIP 10Be flux17 (on the GICC05 timescale15) and new WAIS Divide 06A ice core (WDC) 10Be concentration data (on the WD2014 timescale16) to 14C from IntCal13 using a Bayesian wiggle match approach to estimate the likely time offset between the ice and radiocarbon timescales (Supplementary Data 5). The 10Be data are compared to the atmospheric 14C production rate (p14C) simulated from IntCal13 Δ14C data using a Box-Diffusion carbon-cycle model91 initialized with pre-industrial boundary conditions assuming a constant carbon cycle. Even though this assumption does not necessarily hold for climatically unstable period such as the YD, sensitivity experiments8 have shown that centennial-scale variations in atmospheric 14C (such as those analysed in this study) are largely unaffected by the dynamics of the carbon cycle and predominantly reflect solar cycle modulations. To estimate the uncertainty associated with carbon cycle effects on the simulated p14C, we performed four additional sensitivity experiments where the carbon cycle was perturbed by increasing or decreasing the air-sea gas exchange and the ocean diffusivity parameter by 50%, individually. These tests resulted in an uncertainty estimate of ± 2%, which was added to the nominal IntCal13 Δ14C error.
Prior to wiggle matching, the cosmogenic nuclide records were bandpass filtered between cut-off frequencies of 1/100 and 1/500 years to focus on the common (likely solar-induced) secular changes in production. Low-pass filtering eliminates short-term weather noise and sampling inconsistencies between the data sets. On the other hand, high-pass filtering minimizes systematic errors associated with millennial 10Be production rate variability and uncertainties accompanying longer-term climate and carbon cycle changes92. We also apply a 1-year lag for 10Be to account for the average delay between production and deposition, and ultimately scaled the 10Be time series by a factor of 0.7 to optimize the fit with the p14C data.
Bayesian wiggle matching was performed on the filtered and scaled data following the mathematical formulations of ref. 93 and methodology of ref. 8 (Supplementary Fig. 6). We sequentially offset-minimized the 10Be data against p14C with moving time windows of length of 1000 years at 20-year time steps allowing for leads/lags ranging ±150 years applied to the ice core timescales relative to IntCal13. The probabilities of each individual lead-lag measurement were combined to obtain an overall timescale shift likelihood at each time step between GICC05 and IntCal13, and between WD2014 and IntCal13, respectively. The results allowed us to extend the GICC05-IntCal13 and WD2014-IntCal13 timescale transfer functions of ref. 8 and ref. 16, respectively, beyond the Holocene until ∼14,500 years BP. The estimated timescale offsets are well within the stated uncertainties of the ice core timescales (Supplementary Fig. 6), while the associated uncertainties estimates are consistent with the previous results8. We also note that our timescale transfer functions are in excellent agreement with previously published estimates over the interval of overlap8,16. Finally, Greenlandic and Antarctic ice core records presented in this study were consistently placed on the IntCal13 timescale, whereas beyond the limits of the synchronization described above, the cumulative age uncertainties associated with each ice-core timescale were reset. Note that when applying the timescale correction to WDC CO2 data3, ice-age–gas-age differences16 were fully accounted for.
Lipid biomarkers
One hundred and sixty-nine sediment samples were collected from core MD99-2284 at 1-cm intervals. The samples (∼3 cm3) were freeze-dried, homogenized, and lipids were extracted from the sediments via sonication with a 9:1 mixture (v/v) of dichloromethane:methanol for 20 min. The procedure was repeated three times and supernatants were combined. Total lipid extracts (TLE) were evaporated in a Biotage TurboVap under a stream of nitrogen gas. TLE were separated into constituent lipids by flash silica-gel chromatography (pre-cleaned 100% active silica gel) in hexane-saturated columns using as eluting solvents hexane (F1), dichloromethane (F2), 3:1 hexane:ethyl acetate (F3), and methanol (F4), separately.
To monitor changes in terrestrial meltwater discharge and paleo sea-ice conditions at the coring site, we quantified the abundance of long-chain n-alkanes and mono-unsaturated C25 highly branched isoprenoid alkene (IP25)94, respectively, contained in the F1 fraction of 169 samples (Fig. 2 and Supplementary Fig. 7). The latter biomarker is biosynthesised by Arctic sea ice diatoms during spring and it has been widely used for paleoceanographic reconstructions in this sector of the North Atlantic58,95. The F1 fraction was evaporated and transferred to GC vials. Both n-alkanes (C13-C39) and IP25 were identified and quantified on an Agilent 7890A gas chromatograph (GC) equipped with a mass selective detector (MSD). Samples were measured with an Agilent HP-5ms column (30 m, 0.25 mm ID, 0.25 μm film thickness). The GC oven was held for 1.5 min at 60 °C and ramped to 300 °C at 10 °C min-1. After 1 min it was ramped to 320 °C at 2 °C min−1 and held for 10 min. We used a PTV injector in splitless mode, with samples injected at 60 °C and the injector temperature immediately ramped to 320 °C at 900 °C min−1. All the compounds were identified in SIM mode based on mass spectra and retention times from the literature (monitoring ions m/z 57.1, 85.1, 99.1, 245.3, 346.3, 348.3, 350.3). A standard mixture of n-alkanes (C10-C40) was run after every sixth sample injection, and was used to evaluate instrument drift and correct for mass-dependent differences in the response factor of the MSD. Concentrations of individual n-alkanes were calculated based on the comparison of the adjusted peak areas relative to both C36 alkane and 5α-androstane used as an internal standard, and each added to samples before GC analysis. Similarly, IP25 was quantified based on comparison to the internal standard of 9-octyl-heptadec-8-ene (9-OHD)96 added before GC analysis. The measurement accuracy of both n-alkanes and IP25 was determined by monitoring the relative response factor for n-C36 and 5α-androstane, and for IP25 and 9-OHD (over a range of concentrations), respectively, every seven samples (1σ = ± 3.2%).
To quantify changes in marine surface water phytoplankton productivity, we measured the concentration of C37 unsaturated alkenones97 contained in the F2 fraction of 154 samples (Supplementary Fig. 7). Alkenones were recovered by liquid-phase separation using toluene, and diluted for analysis with a Thermo TRACE Ultra Gas Chromatography Flame Ionization Detector (GC-FID) equipped with PTV injector operated in splitless mode. Samples were measured with a DB-1 column (60 m, 0.25 μm ID, 0.1 μm film thickness). Lastly, an in-house alkenone standard with known relative concentrations of C37 unsaturated alkenones was injected every six samples to monitor the instrument performance and analytical precision of the alkenone-unsaturated index UK37 (ref. 98) (1σ = ± 0.0072).
Changes in sea-ice cover were finally reconstructed by calculating a phytoplankton-IP25 index (PIP25)99 (Fig. 2) following the equation
$$PIP_{25} = [IP_{25}]\frac{{[IP_{25}]}}{{\left[ {IP_{25}} \right] + ([P] \times c)}}$$
whereby [IP25] and [P] represent the respective concentrations of IP25 and the phytoplankton biomarker P, while c is a balance factor calculated from the ratio of the mean IP25 to P concentrations.
Breakpoint analysis
To estimate the timing of the transitions recorded in the proxy data, we used a Monte Carlo version of the function Segmented in the R package 'Segmented'100, i.e., a piecewise linear fitting regression method that determines the breakpoints of two lines. To combine the analytical and age uncertainties associated with each proxy record, while at the same time enforcing monotonicity in the age estimates, the data sets were first modelled using the Random Walk Model described above, whereby the model was fitted directly to the observations. We subsequently randomly resampled (5000 times) the modelled data within its uncertainty, assuming a Gaussian error distribution. We then estimated the breakpoint structure of each Monte Carlo realization using the Segmented function run with 100 bootstrap iterations. The inferred likelihood distribution of each breakpoint allowed us quantifying the timing uncertainty for the start and end of the transitions under investigation (Fig. 4; Supplementary Table 1).
All the climate model output and R codes used for the numerical procedures are available from the main author upon reasonable request.
The source data underlying Fig. 2c, d, Supplementary Figures 4 and 6d, f and h (tephra geochemistry, 14C ages, WDC 10Be data, ocean-atmosphere 14C disequilibrium estimates, IntCal13-GICC05 and IntCal13-WD2014 age offsets) are provided as a Source Data file along the online version of this article on the publisher's website and in the PANGAEA paleoclimate data archive
Broecker, W. S., Peteet, D. M. & Rind, D. Does the ocean–atmosphere system have more than one stable mode of operation? Nature 315, 21–26 (1985).
McManus, J. F., Francois, R., Gherardi, J.-M., Keigwin, L. D. & Brown-Leger, S. Collapse and rapid resumption of Atlantic meridional circulation linked to deglacial climate changes. Nature 428, 834–837 (2004).
Marcott, S. A. et al. Centennial-scale changes in the global carbon cycle during the last deglaciation. Nature 514, 616–619 (2014).
Ito, T. & Follows, M. J. Preformed phosphate, soft tissue pump and atmospheric CO2. J. Mar. Res. 63, 813–839 (2005).
Schmittner, A. & Galbraith, E. D. Glacial greenhouse-gas fluctuations controlled by ocean circulation changes. Nature 456, 373–376 (2008).
Waelbroeck, C. et al. The timing of the last deglaciation in North Atlantic climate records. Nature 412, 724–727 (2001).
Muscheler, R., Adolphi, F. & Knudsen, M. F. Assessing the differences between the IntCal and Greenland ice-core time scales for the last 14,000 years via the common cosmogenic radionuclide variations. Quat. Sci. Rev. 106, 81–87 (2014).
Article ADS Google Scholar
Adolphi, F. & Muscheler, R. Synchronizing the Greenland ice core and radiocarbon timescales over the Holocene-Bayesian wiggle-matching of cosmogenic radionuclide records. Clim. Past. 12, 15–30 (2016).
Hansen, B. & Østerhus, S. North Atlantic-Nordic Seas exchanges. Progress. Oceanogr. 45, 109–208 (2000).
Meland, M. Y., Dokken, T. M., Jansen, E. & Hevrøy, K. Water mass properties and exchange between the Nordic seas and the northern North Atlantic during the period 23-6 ka: Benthic oxygen isotopic evidence. Paleoceanography 23, https://doi.org/10.1029/2007PA001416(2008).
Crocket, K. C., Vance, D., Gutjahr, M., Foster, G. L. & Richards, D. A. Persistent Nordic deep-water overflow to the glacial North. Atl. Geol. 39, 515–518 (2011).
Ezat, M. M. et al. Ventilation history of Nordic Seas overflows during the last deglacial period revealed by species-specific benthic foraminiferal 14C dates. Paleoceanography 32, 172–181 (2017).
Thornalley, D. J. R. et al. A warm and poorly ventilated deep Arctic Mediterranean during the last glacial period. Science 349, 706–710 (2015).
Reimer, P. J. et al. IntCal13 and Marine13 Radiocarbon Age Calibration Curves 0–50,000 Years cal BP. Radiocarbon 55, 1869–1887 (2013).
Rasmussen, S. O. et al. A new Greenland ice core chronology for the last glacial termination. J. Geophys. Res. Atmos. 111, https://doi.org/10.1029/2005JD006079 (2006).
Sigl, M. et al. The WAIS Divide deep ice core WD2014 chronology - Part 2: Annual-layer counting (0-31 ka BP). Clim. Past. 12, 769–786 (2016).
Adolphi, F. et al. Persistent link between solar activity and Greenland climate during the Last Glacial Maximum. Nat. Geosci. 7, 662–666 (2014).
Cheng, J. et al. Simulated two-stage recovery of Atlantic meridional overturning circulation during the last deglaciation. Abrupt Clim. Chang. Mech. Patterns, Impacts, Am. Geophys. Union 75–92 (2011).
Renold, M., Raible, C. C., Yoshimori, M. & Stocker, T. F. Simulated resumption of the North Atlantic meridional overturning circulation—Slow basin-wide advection and abrupt local convection. Quat. Sci. Rev. 29, 101–112 (2009).
Bauch, H. A. et al. A multiproxy reconstruction of the evolution of deep and surface waters in the subarctic Nordic seas over the last 30,000 yr. Quat. Sci. Rev. 20, 659–678 (2001).
Ritz, S. P., Stocker, T. F. & Müller, S. A. Modeling the effect of abrupt ocean circulation change on marine reservoir age. Earth. Planet. Sci. Lett. 268, 202–211 (2008).
Wilson, D. J., Crocket, K. C., Van De Flierdt, T., Robinson, L. F. & Adkins, J. F. Dynamic intermediate ocean circulation in the North Atlantic during Heinrich Stadial 1: A radiocarbon and neodymium isotope perspective. Paleoceanography 29, 1072–1093 (2015).
Adkins, R. L., Cheng, C. D., Boyle, H., Druffel & Edwards Deep-Sea coral evidence for rapid change in ventilation of the deep north atlantic 15,400 years Ago. Science 280, 725–728 (1998).
Eltgroth, S. F., Adkins, J. F., Robinson, L. F., Southon, J. & Kashgarian, M. A deep-sea coral record of North Atlantic radiocarbon through the Younger Dryas: Evidence for intermediate water/deepwater reorganization. Paleoceanography 21, https://doi.org/10.1029/2005PA001192(2006).
Thiagarajan, N., Subhas, A. V., Southon, J. R., Eiler, J. M. & Adkins, J. F. Abrupt pre-Bølling-Allerød warming and circulation changes in the deep ocean. Nature 511, 75–78 (2014).
Chen, T. et al. Synchronous centennial abrupt events in the ocean and atmosphere during the last deglaciation. Science 349, 1537–1541 (2015).
Hain, M. P., Sigman, D. M. & Haug, G. H. Distinct roles of the Southern Ocean and North Atlantic in the deglacial atmospheric radiocarbon decline. Earth. Planet. Sci. Lett. 394, 198–208 (2014).
Singarayer, J. S. et al. An oceanic origin for the increase of atmospheric radiocarbon during the Younger Dryas. Geophys. Res. Lett. 35, https://doi.org/10.1029/2008GL034074 (2008).
Hughen, K. 14C activity and global carbon cycle changes over the past 50,000 years. Science 303, 202–207 (2004).
Barker, S. et al. Interhemispheric Atlantic seesaw response during the last deglaciation. Nature 457, 1097–1102 (2009).
Shakun, J. D. et al. Global warming preceded by increasing carbon dioxide concentrations during the last deglaciation. Nature 484, 49–54 (2012).
Muschitiello, F. & Wohlfarth, B. Time-transgressive environmental shifts across Northern Europe at the onset of the Younger Dryas. Quat. Sci. Rev. 109, 49–56 (2015).
Muschitiello, F. et al. Fennoscandian freshwater control on Greenland hydroclimate shifts at the onset of the Younger Dryas. Nat. Commun. 6, 8939 (2015).
Dokken, T. M., Nisancioglu, K. H., Li, C., Battisti, D. S. & Kissel, C. Dansgaard-Oeschger cycles: Interactions between ocean and sea ice intrinsic to the Nordic seas. Paleoceanography 28, 491–502 (2013).
Barker, S. et al. Icebergs not the trigger for North Atlantic cold events. Nature 520, 333–336 (2015).
Bakke, J. et al. Rapid oceanic and atmospheric changes during the Younger Dryas cold period. Nat. Geosci. 2, 202–205 (2009).
Lane, C. S., Brauer, A., Blockley, S. P. E. & Dulski, P. Volcanic ash reveals time-transgressive abrupt climate change during the Younger Dryas. Geology 41, 1251–1254 (2013).
Carlson, A. E. et al. Subtropical Atlantic salinity variability and Atlantic meridional circulation during the last deglaciation. Geology 36, 991–994 (2008).
Ng, H. C. et al. Coherent deglacial changes in western Atlantic Ocean circulation. Nat. Commun. 9, 2947 (2018).
Valley, S., Lynch-Stieglitz, J. & Marchitto, T. M. Timing of Deglacial AMOC Variability From a High-Resolution Seawater Cadmium Reconstruction. Paleoceanography 32, 1195–1203 (2017).
Lynch-Stieglitz, J. et al. Muted change in Atlantic overturning circulation over some glacial-aged Heinrich events. Nat. Geosci. 7, 144–150 (2014).
Yu, E.-F., Francois, R. & Bacon, M. P. Similar rates of modern and last-glacial ocean thermohaline circulation inferred from radiochemical data. Nature 379, 689–694 (1996).
Henry, L. G. et al. North Atlantic ocean circulation and abrupt climate change during the last glaciation. Science 353, 470–474 (2016).
Liu, Z. et al. Transient simulation of last deglaciation with a new mechanism for Bolling-Allerod warming. Science 325, 310–314 (2009).
Menviel, L., Timmermann, A., Friedrich, T. & England, M. H. Hindcasting the continuum of Dansgaard-Oeschger variability: Mechanisms, patterns and timing. Clim. Past. 10, 63–77 (2014).
Broecker, W. S., Bond, G., Klas, M., Bonani, G. & Wolfli, W. A salt oscillator in the glacial Atlantic? 1. Concept. Paleoceanogr. 5, 469–477 (1990).
Peltier, W. R. & Vettoretti, G. Dansgaard-Oeschger oscillations predicted in a comprehensive model of glacial climate: A 'kicked' salt oscillator in the Atlantic. Geophys. Res. Lett. 41, 7306–7313 (2014).
Renssen, H. et al. Multiple causes of the Younger Dryas cold period. Nat. Geosci. 8, 946–949 (2015).
Lehner, F., Born, A., Raible, C. C. & Stocker, T. F. Amplified inception of European little Ice Age by sea ice-ocean-atmosphere feedbacks. J. Clim. 26, 7586–7602 (2013).
Drijfhout, S., Gleeson, E., Dijkstra, H. A. & Livina, V. Spontaneous abrupt climate change due to an atmospheric blocking-sea-ice-ocean feedback in an unforced climate model simulation. Proc. Natl Acad. Sci. 110, 19713–19718 (2013).
Zhang, J. et al. Asynchronous warming and δ18O evolution of deep Atlantic water masses during the last deglaciation. Proc. Natl Acad. Sci. USA 114, 11075–11080 (2017).
Marcott, S. A. et al. Ice-shelf collapse from subsurface warming as a trigger for Heinrich events. Proc. Natl Acad. Sci. USA 108, 13415–13419 (2011).
Muschitiello, F. et al. Timing of the first drainage of the Baltic Ice Lake synchronous with the onset of Greenland Stadial 1. Boreas 45, 322–334 (2016).
Schenk, F. et al. Warm summers during the Younger Dryas cold reversal. Nat. Commun. 9, 1634 (2018).
Sherriff-Tadano, S., Abe-Ouchi, A., Yoshimori, M., Oka, A. & Chan, W. Le. Influence of glacial ice sheets on the Atlantic meridional overturning circulation through surface wind change. Clim. Dyn. 50, 2881–2903 (2018).
Ebbesen, H. & Hald, M. Unstable Younger Dryas climate in the northeast North Atlantic. Geology 32, 673–676 (2004).
Oksman, M. et al. Younger Dryas ice margin retreat triggered by ocean surface warming in central-eastern Baffin Bay. Nat. Commun. 8, 1017 (2017).
Cabedo-Sanz, P., Belt, S. T., Knies, J. & Husum, K. Identification of contrasting seasonal sea ice conditions during the Younger Dryas. Quat. Sci. Rev. 79, 74–86 (2013).
Pearce, C. et al. Ocean lead at the termination of the Younger Dryas cold spell. Nat. Commun. 4, 1664 (2013).
Caesar, L., Rahmstorf, S., Robinson, A., Feulner, G. & Saba, V. Observed fingerprint of a weakening Atlantic Ocean overturning circulation. Nature 556, 191 (2018).
Thornalley, D. J. R. et al. Anomalously weak Labrador Sea convection and Atlantic overturning during the past 150 years. Nature 556, 227 (2018).
Eldevik, T. et al. A brief history of climate e the northern seas from the last Glacial Maximum to global warming. Quat. Sci. Rev. 106, 225–246 (2014).
Wohlfarth, B. et al. Climate and environment in southwest Sweden 15.3 – 11.3 cal. ka BP. Boreas https://doi.org/10.1111/bor.12310 (2018).
Gustafsson, M., Rayner, D. & Chen, D. Extreme rainfall events in southern Sweden: Where does the moisture come from? Tellus, Ser. A Dyn. Meteorol. Oceanogr. 62, 605–616 (2010).
Rach, O., Brauer, A., Wilkes, H. & Sachse, D. Delayed hydrological response to Greenland cooling at the onset of the Younger Dryas in western Europe. Nat. Geosci. 7, 109–112 (2014).
Muschitiello, F., Andersson, A., Wohlfarth, B. & Smittenberg, R. H. The C20 highly branched isoprenoid biomarker - A new diatom-sourced proxy for summer trophic conditions? Org. Geochem. 81, 27–33 (2015).
Muschitiello, F. Deglacial impact of the Scandinavian Ice Sheet on the North Atlantic climate system. (Stockholm University, 2016).
Rasmussen, C. E. & Williams, C. K. I. Gaussian processes for machine learning. Int. J. Neural. Syst. 14, 69–106 (2004).
Abbott, P. M. & Davies, S. M. Volcanism and the Greenland ice-cores: The tephra record. Earth-Sci. Rev. 115, 173–191 (2012).
Lohne, Ø. S., Mangerud, J. & Birks, H. H. IntCal13 calibrated ages of the Vedde and Saksunarvatn ashes and the Younger Dryas boundaries from Kråkenes, western Norway. J. Quat. Sci. 29, 506–507 (2014).
Matthews, I. P. et al. New age estimates and climatostratigraphic correlations for the borrobol and penifiler tephras: Evidence from Abernethy Forest, Scotland. J. Quat. Sci. 26, 247–252 (2011).
MacLeod, A., Matthews, I. P., Lowe, J. J., Palmer, A. P. & Albert, P. G. A second tephra isochron for the Younger Dryas period in northern Europe: The Abernethy Tephra. Quat. Geochronol. 28, 1–11 (2015).
Turney, C. S. M. Extraction of rhyolitic component of Vedde microtephra from minerogenic lake sediments. J. Paleolimnol. 19, 199–206 (1998).
Mortensen, A. K., Bigler, M., Grönvold, K., Steffensen, J. P. & Johnsen, S. J. Volcanic ash layers from the last glacial termination in the NGRIP ice core. J. Quat. Sci. 20, 209–219 (2005).
Small, D. & Fabel, D. A Lateglacial 10Be production rate from glacial lake shorelines in Scotland. J. Quat. Sci. 30, 509–513 (2015).
Simstich, J., Sarnthein, M. & Erlenkeuser, H. Paired δ18O signals of Neogloboquadrina pachyderma (s) and Turborotalita quinqueloba show thermal stratification structure in Nordic Seas. Mar. Micropaleontol. 48, 107–125 (2003).
Wacker, L., Fülöp, R. H., Hajdas, I., Molnár, M. & Rethemeyer, J. A novel approach to process carbonate samples for radiocarbon measurements with helium carrier gas. in Nuclear Instruments and Methods in Physics. Res., Sect. B: Beam Interact. Mater. At. 294, 214–217 (2013).
Hajdas, I., Bonani, G., Herrgesell, S. & Millie, Z. 14C ages of Ostracode from Pleistocene lake of sediments of the Western Great Basin, USA - Results of progressive acid leaching. Radiocarbon 46, 189–200 (2004).
Stuiver, M. & Polach, H. A. Discussion: reporting of 14 C. data. Radiocarb. 19, 355–363 (1977).
Freeman, E., Skinner, L. C., Reimer, R., Scrivner, A. & Fallon, S. Graphitization of small carbonate samples for paleoceanographic research at the godwin radiocarbon laboratory, University of Cambridge. Radiocarbon 58, 89–97 (2016).
Skinner, L. C., Fallon, S., Waelbroeck, C., Michel, E. & Barker, S. Ventilation of the Deep Southern Ocean and Deglacial CO2 Rise. Science 328, 1147–1151 (2010).
Barker, S., Broecker, W., Clark, E. & Hajdas, I. Radiocarbon age offsets of foraminifera resulting from differential dissolution and fragmentation within the sedidmentary bioturbated zone. Paleoceanography 22, https://doi.org/10.1029/2006PA001354 (2007).
Broecker, W., Mix, A., Andree, M. & Oeschger, H. Radiocarbon measurements on coexisting benthic and planktic foraminifera shells: potential for reconstructing ocean ventilation times over the past 20 000 years. Nucl. Inst. Methods Phys. Res. B 5, 331–339 (1984).
Franke, J., Schulz, M., Paul, A. & Adkins, J. F. Assessing the ability of the 14C projection‐age method to constrain the circulation of the past in a 3‐D ocean model. Geochemistry, Geophys. Geosystems 9, https://doi.org/10.1029/2008GC001943 (2008).
Blackwell, P. G. & Buck, C. E. Estimating radiocarbon calibration curves. Bayesian Anal. 3, 225–248 (2008).
Heaton, T. J., Blackwell, P. G. & Buck, C. E. A Bayesian approach to the estimation of radiocarbon calibration curves: The IntCal09 methodology. Radiocarbon 51, 1151–1164 (2009).
Niu, M., Heaton, T. J., Blackwell, P. G. & Buck, C. E. The Bayesian approach to radiocarbon calibration curve estimation: the IntCal13, Marine13, and SHCal13 methodologies. Radiocarbon 55, 1905–1922 (2013).
Finkel, R. C. & Nishiizumi, K. Beryllium 10 concentrations in the Greenland Ice Sheet Project 2 ice core from 3-40 ka. J. Geophys. Res. Ocean. 102, 26699–26706 (1997).
Woodruff, T. E., Welten, K. C., Caffee, M. W. & Nishiizumi, K. Interlaboratory comparison of 10Be concentrations in two ice cores from Central West Antarctica. in Nuclear Instruments and Methods in Physics. Res., Sect. B: Beam Interact. Mater. At. 294, 77–80 (2013).
Nishiizumi, K. et al. Absolute calibration of 10Be AMS standards. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 258, 403–413 (2007).
Siegenthaler, U., Heimann, M. & Oeschger, H. 14C variations caused by changes in the global carbon cycle. Radiocarbon 22, 177–191 (1980).
Köhler, P., Muscheler, R. & Fischer, H. A model-based interpretation of low-frequency changes in the carbon cycle during the last 120,000 years and its implications for the reconstruction of atmospheric Δ14C. Geochemistry, Geophys. Geosystems 7, https://doi.org/10.1029/2005GC001228 (2006).
Ramsey, C. B., van der Plicht, J. & Weninger, B. 'Wiggle Matching' Radiocarbon Dates. Radiocarbon 43, 381–389 (2001).
Belt, S. T. et al. A novel chemical fossil of palaeo sea ice: IP25. Org. Geochem. 38, 16–27 (2007).
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).
Belt, S. T. et al. A reproducible method for the extraction, identification and quantification of the Arctic sea ice proxy IP25 from marine sediments. Anal. Methods 4, 705 (2012).
Prahl, F. G., Muehlhausen, L. A. & Zahnle, D. L. Further evaluation of long-chain alkenones as indicators of paleoceanographic conditions. Geochim. Cosmochim. Acta 52, 2303–2310 (1988).
Brassell, S. C. et al. Palaeoclimatic signals recognized by chemometric treatment of molecular stratigraphic data. Org. Geochem. 10, 649–660 (1986).
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).
Muggeo, V. M. R. Estimating regression models with unknown break-points. Stat. Med. 22, 3055–3071 (2003).
Johnsen, S. J. et al. High-resolution record of Northern Hemisphere climate extending into the last interglacial period. Nature 431, 147–151 (2004).
Article ADS PubMed CAS Google Scholar
Peltier, W. R., Argus, D. F. & Drummond, R. Space geodesy constrains ice age terminal deglaciation: The global ICE‐6G_C (VM5a) model. J. Geophys. Res. Solid Earth 120, 450–487 (2015).
Menviel, L., Timmermann, A., Timm, O. E. & Mouchet, A. Deconstructing the Last Glacial termination: The role of millennial and orbital-scale forcings. Quat. Sci. Rev. 30, 1155–1172 (2011).
Schmittner, A. & Lund, D. C. Early deglacial Atlantic overturning decline and its role in atmospheric CO2 rise inferred from carbon isotopes (δ13C). Clim. Past. 11, 135–152 (2015).
Dreger, D. Decadal to Centennial Scale Sediment Records of Ice Advance on the Barents Shelf and Meltwater Discharge Into the Northeastern Norwegian Sea Over the Last 40 Kyr. https://doi.org/10.2312/reports-ifg.1999.3 (1999).
Polyak, L. & Mikhailov, V. Post-glacial environments of the southeastern Barents Sea: foraminiferal evidence. Geol. Soc. Lond., Spec. Publ. 111, 323–337 (1996).
Eglinton, G. & Hamilton, R. J. Leaf Epicuticular Waxes. Sci. (80-.). 156, 1322–1335 (1967).
Cao, L., Fairbanks, R. G., Mortlock, R. A. & Risk, M. J. Radiocarbon reservoir age of high latitude North Atlantic surface water during the last deglacial. Quat. Sci. Rev. 26, 732–742 (2007).
Buizert, C. et al. Greenland temperature response to climate forcing during the last deglaciation. Science 345, 1177–1180 (2014).
Matsumoto, K. & Yokoyama, Y. Atmospheric Δ14C reduction in simulations of Atlantic overturning circulation shutdown. Glob. Biogeochem. Cycles 27, 296–304 (2013).
We are grateful to D.I. Blindheim for assistance with foraminifera sampling, H. Sadatzky for assistance with sampling for lipid biomarker analyses, and S. Rasmussen for helpful discussions. The research leading to these results has received funding from the European Research Council under the European Community's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement 610055 as part of the 'ice2ice' project. The 10Be measurements in the WAIS Divide core were funded by NSF grants ANT 0839042 (to M.W.C) and 0839137 (to K.C.W.). Additional funding by the Columbia Climate Center of Columbia University (to F.M and W.J.D.) and the Vetlesen Foundation (W.J.D) is gratefully acknowledged. This study is a contribution to the INTIMATE project.
Department of Geography, University of Cambridge, Cambridge, CB2 3EN, UK
Francesco Muschitiello
Lamont-Doherty Earth Observatory, Columbia University, Palisades, NY, 10964, USA
Francesco Muschitiello, William J. D'Andrea & Nicole deRoberts
NORCE Norwegian Research Centre and Bjerknes Centre for Climate Research, 5007, Bergen, Norway
Francesco Muschitiello, Margit H. Simon & Trond M. Dokken
College of Earth, Ocean, and Atmospheric Sciences, Oregon State University, Corvallis, OR, 97331-5503, USA
Andreas Schmittner
School of Mathematics and Statistics, University of Sheffield, Sheffield, S3 7RH, UK
Timothy J. Heaton
Department of Geology, College of William and Mary, Williamsburg, VA, 23187, USA
Nicholas L. Balascio
Department of Physics and Astronomy, Purdue University, West Lafayette, IN, 47907, USA
Marc W. Caffee & Thomas E. Woodruff
Department of Earth, Atmospheric, and Planetary Sciences, Purdue University, West Lafayette, IN, 47907, USA
Marc W. Caffee
Space Sciences Laboratory, University of California, Berkeley, CA, 94720, USA
Kees C. Welten
Godwin Laboratory for Palaeoclimate Research, Department of Earth Sciences, University of Cambridge, Cambridge, CB2 3EQ, UK
Luke C. Skinner
William J. D'Andrea
Nicole deRoberts
Thomas E. Woodruff
Margit H. Simon
Trond M. Dokken
F.M. conceived and designed the study, and wrote the first draft of the manuscript. W.J.D. contributed to the quantification and evaluation of biomarker data. A.S. designed and performed the climate model experiments. T.J.H. developed the Random Walk Model. N.L.B. performed the tephra analyses. N.D. together with F.M. performed the biomarker analyses. M.W.C., T.E.W., and K.W measured the WAIS Divide Ice Core 10Be data. L.C.S. provided complementary 14C data for core MD99-2284. M.H.S. provided insight into regional paleoenvironment. T.M.D provided the background data for core MD99-2284 and together with F.M and W.J.D. secured necessary financial support. All authors contributed to the interpretation of the results and editing of the manuscript.
Correspondence to Francesco Muschitiello.
Journal peer review information: Nature Communications thanks Lowell Stott and Anders Svensson for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Data 1
Muschitiello, F., D'Andrea, W.J., Schmittner, A. et al. Deep-water circulation changes lead North Atlantic climate during deglaciation. Nat Commun 10, 1272 (2019). https://doi.org/10.1038/s41467-019-09237-3
Central Mediterranean rainfall varied with high northern latitude temperatures during the last deglaciation
Andrea Columbu
Christoph Spötl
Jo De Waele
Communications Earth & Environment (2022)
Synchronous vegetation response to the last glacial-interglacial transition in northwest Europe
Stefan Engels
Christine S. Lane
Frederike Wagner-Cremer
Evolution of the North Atlantic Current and Barents Ice Sheet as revealed by grain size populations in the northern Norwegian Sea during the last 60 ka
Weiguo Wang
Mengwei Zhao
Acta Oceanologica Sinica (2021)
Onset and termination of Heinrich Stadial 4 and the underlying climate dynamics
Hai Cheng
Yao Xu
R. Lawrence Edwards
Eurasian Ice Sheet collapse was a major source of Meltwater Pulse 1A 14,600 years ago
Jo Brendryen
Haflidi Haflidason
Bjarte Hannisdal
Nature Geoscience (2020)


text:
Home / Historical / How to Ensnare a Highlander: The MacGregor Lairds, Book 2
How to Ensnare a Highlander: The MacGregor Lairds, Book 2
When Lady Elizabet Harding's family is targeted by the notorious Highland Highwayman, she finds herself wounded and an unwelcome guest of the charming outlaw who'd ruined her reputation earlier that day. Elizabet only wants to get even for her family's downfall, but the more time she spends with the scoundrel, the more complicated her feelings become.
Genre: Historical. Author: McLean Michelle.
By day, Laird John MacGregor graces the court of Charles II. By night, he exacts revenge on his enemies – one of whom is betrothed to the fiery heiress who has, quite literally, fallen into his arms. The daughter of one enemy and promised to another, Elizabet should be the last person John wants. Yet even as she exasperates him, she is all he has ever desired for himself. When her life is put in danger, he will do anything necessary to protect the women he never expected to love…even if it means they cannot be together.
Be the first to review "How to Ensnare a Highlander: The MacGregor Lairds, Book 2" Cancel reply
Valentine: The Brotherhood of Fallen Angels
Whispers of Shadow & Flame: Earthsinger Chronicles, Book 2
A Wicked Way to Win an Earl: The Sutherland Scandals, Book 1
Knight of Her Heart: Conquering the Heart Series	$29.95
The Virgin and the Rogue: The Rogue Files	$19.95
Last Gentleman Standing	$16.95

redpajama_set_name: RedPajamaCommonCrawl


The MAXFire® 100 Series gas igniter incorporates the "protected primary combustion zone" design in each of the 5 capacities available. The ignition process starts inside the guide tube away from the burner environment resulting in an extremely stable and reliable igniter. Capacities (nominal) range from 1 mmBtu/hr to 50 mmBtu/hr.
MAXFire® gas igniters provide the most reliable source of energy for main flame ignition, flame stabilization, and boiler warm-up.
The MAXFire® Series 100 gas igniters can be used as a Class 1, Class 2, or Class 3 igniter, depending on burner capacity. The MAXFire® igniter is easily installed as a retrofit into most burner types based on its size, capacity, and reliable operation.
MAXFire® gas igniters consist of four major components: a guide tube, a gas transport tube, a spark rod, and a stainless steel diffuser for pilot stabilization. The guide tube houses the spark rod and the gas transport tube, which are easily removed from the guide tube for maintenance purposes by releasing a "V" band retainer, requiring no tools. Cooling/combustion airflow is directed down the guide tube to support and maintain the primary combustion zone during operation. The controlled mixture of fuel and air in the primary combustion zone results in a highly reliable, stable igniter flame. While out of service, airflow cools the igniter and prevents debris from migrating into the end of the guide tube.
MAXFire® gas igniters come standard with Forney's High Energy Spark Igniter (HESI) as the spark source. The HESI produces a powerful 12-joule spark approximately three times per second at the primary combustion zone. The HESI spark tip is a nonfouling, surface-gap tip which is replaceable as a unit by a threaded connection. Optional equipment for the MAXFire® includes a High Tension Spark (HTSI) source in place of the HESI, mounting hardware, hoses, flame rod and an integral or externally mounted optical flame detector.
Capacity range from 0.3 to 50 MBtu/hr.
The nonfouling HESI comes standard. HTSI provides high performance at an economical price.


The way to design Home Design Maker? To remedy these inquiries, we have to primary should recognize with this. Home Design Maker can be an idea or perhaps suggestions for decorate your house or perhaps room so it seems to be practical in addition to pleasing to the eye. Building Home Design Maker is not complicated. You recently have to have a tiny ingenuity and also a powerful will, and you will work with your imagination right after finding photos available right here.
Clearly, Home Design Maker can be developed in addition to designed by your current loving. Often simply by using a smart style, along with modify the format with the bedroom or residence consequently be more successful to put together and check pretty.
Truly creations to switch one thing to generate that search stunning will not come to an end when you are developing Home Design Maker, particularly if the property or perhaps include your individual area. If you believe tired of your purchase, you can right away change it for a liking. Basically, creating Home Design Maker is almost endless. You cannot think that is simply on a the quality, quantity factor furthermore influences the design of the house or maybe area to produce the item appear more pleasing along with particular.
Beyond the quality in addition to quantity being a aspect with regard to designing a place as well as the house, you can find some other components. Decide what sort of design as well as strategy of the house whilst your room later on. And additional outside components such as photographs within this internet site, plus the concept of your designer residence and locations were being really motivated all of us obviously. The Home Design Maker will not likely work efficiently unless you focus on accuracy and reliability in lots of variables that were explained.
Technology likewise takes on a significant position inside creating Home Design Maker. The actual technology we can easily always be made it easier for both from the instrument and not making it easier for the perform creating a property as well as space.
What this content identified earlier mentioned can be a minor factors which can help anyone design and style as well as develop a house or even a place as occupancy warranted being filled. As well as it is possible to grab works of art or maybe tips associated with that which you view on this site both associated with Articles and Pictures Home Design Maker as a result build a amazing strategy regarding the idea of residence as well as place particular. You may also conserve photographs are generally fixed perfectly in most write-up from this web page being located in a laptop or computer, notebook or cell phone while your own stored data file to develop a home as well as bedroom in case whenever they want you will need. Might be valuable for all of us.
Charming Download DreamPlan Home Design Software .
Marvelous :D Online Home Design, 3d Home Design, .
Nice Home Design Maker Of Inspiring House Blueprint Generator Building Brilliant Floor Plan Minecraft .
Delightful Free Home Design Software Mac Modern Ideas ... .
House Floor Plans App Outstanding House Plan Builder Free Best Of Home Design Maker Elegant Outstanding .
Great More 6 Unique Online House Design Maker .
... Interesting 3d House Plan Maker Ideas Best Inspiration Home Amazing Design ... .
Unique Home Blue Prints And Home Blueprints Maker Home Design Blueprint Home Brilliant Home Design Blueprint .
Superior 3d Home Design Projects Inspiration Home Design Maker And More Bedroom D S More Home Design Single . 3d Home Design ... .
Easy Floor Plan Maker Beautiful Www Home Design Plan Fresh Open Floor House Plans Open Floor Plan .
Home Design Planner House Planner Bedroom Planner Wondrous Home Design Maker House Plan Generator Floor Creator . Home Design ... .
Lovely Online Floor Plan Maker Lovely Design Ideas Floor Planner Line For Modern Home Design .
Meagan Lane, Ritual Home Design, Furniture Maker, Furniture Maker San Diego, Craftsmanship .
Good Floor Plan Design .
Amazing Free Home Design Programs For Pc Proposal Maker Software Lovely Interior Design Software Cost Luxury Of .
Elegant Easy 3d Drawing Software Simple House Plan Software For Mac Best 98 3d Home Design Maker .
Perfect Home Design Maker Home Design Maker, Indem Sie Keine Mittel Gehen Aus Von Designs. .
Exceptional Sweet Home ... .
Captivating Easy 3d Drawing Software Simple House Plan Software For Mac Best 98 3d Home Design Maker .
Beautiful Online Floor Plan Maker Lovely Line Floor Plan Designer Line Home Design Plans Fresh Line Floor 6381 .
Awesome Unique Floor Plan Maker Online For Home Design Ideas Or Floor Plan Maker Online .
Wonderful ... Photos 3d House Floor Plan Maker Decor Design Crafty Simple Plans 12 ... .
High Quality Miraculous Square Rugs For Your Home Design 4x4 .
Simple Floor Plan Software Modern Architecture Wiring Office Plants Maker Creator Home Within 10 ... .
Easy Floor Plan Maker Best Of How To Design A House Plan Awesome Easy Home Design Beautiful 10 2345 .
Online House Design Free Single Family Home Design Online House Design Maker 3d .
3d Home Design Online House Plan Maker Online Inspirational House Plan Studio House Plan Maker Online .
Free Architectural Drawing Software Home Design Interior .
House Designs Exterior With House Plans Fresh House Plans With Interior Emergencymanagementsummit Of House Designs Exterior .
Home Design Maker Blueprint Home Design Awesome Modern House Blueprints Home Design Ideas Maker Images Houses . Home Design Maker ... .
3d House Plans Mac Unique 98 3d Home Design Maker Line Free Floor .
... House Plan Maker Fresh On Impressive Design Example Png Bn ... .
... Blueprint Maker Home Design Layout Plan Home Design Home Design Draw A Floor Plan Draw Floor ... .
Sketches Houses Plan Sketching Home Design Sketch Plans .
Roomsketcher Create Floor Plans And Home Designs Online Rh Roomsketcher Com Free Online 3d House Design Maker .
Design House Online 3d Free .
3d House Maker Home Design Also With A Design A House Also With A Home Interior . 3d House Maker House Design Maker Download Home ... .
... Plain Decoration Draw Your Own House Plans Software Draw Your Own House Plans Software Luxury 98 ... .
Captivating Home Design Planner House Plan Maker Floor Creator Pertaining To Best Planning Software Lovely Home .
Bedroom Layout Ideas For Rectangular Rooms Awesome Fresh Floor Plan Maker Concept Home Interior Design .
Exterior House Design Exterior Home Design Online House Map Elevation Exterior House Design Maker .
Coolest Home Design Maker For Ergonomic Home Remodeling 95 With Home Design Maker .
... Home Design Software Best Of House Floor Plan Maker New 3d Home Plans Home Plan Designer .
Home Design Blueprint Basic Floor Plan Best Of Maker House Plans Photos Simple Small Rectangular Fresh Modern Blueprints .
... Decor Ideas My House 3d Home Design Free Woodworking Timely Best Software The Fornting Designing Onlinent ... .
Home Design Maker House Plan Maker House Floor Plan Maker Home Design Software Free Remodelling .
Sketchup Home Design Best 3d Floor Plan Maker Lovely 3d Floor Plan Awesome Sketchup House Globosynovedades.com.co Valid Sketchup Home Design ... .
3d Home Plan Creator Elegant Home Design Maker Line Haus Plan Maker Home Design Und Dekor Ideen .
3d Floor Plan Maker Unique Free 3d Home Plans Fresh Home Design Plans Home Plan Designer .
Cute Delightful Latest Interior Design Magazine Zaila Mens Or Men Home Designs And Minimalist Environment Apartment_apartment .
Free And Online 3d Home Design Planner Homebyme Rh Home By Me .
3d Floor Plan Maker Inspirational Floor Plan Maker Best Draw A Floor Plan Inspirational Home Design .
House Plan Maker Home Floor Plan Creator Decorating Ideas Simple Home Design Floor Plans .
Az Home Plans Fresh House Floor Plan Maker Inspirational Az Home Plans Fresh Home Design Of .
Luxury House Blueprints Maker Or Dream House Maker House Blueprints Maker Blueprints Design Home Design Blueprint .
Article With Tag: House Floor Plans Maker .
Home Design Maker Free Floor Plan Maker Design House Floor Plans Creative Cheap Modern Home .
Five Various Ways To Do Free Online House Plans | Free Online House Plans .
Wooden Beams In House Baffling House Plans 3d S New 3d Home Plans Fresh Home Design . .
Home Plan Maker Online Unique What Is The Best Home Design Consumer Software Of Home Plan Maker Online .
3d Floor Plan 2D Floor Plan 3D Site Plan Design 3D Floor Plan Modeling Designer And .
1108 Sq Ft Single Floor Contemporary Home Design Free Online House Planner 3d Online House Design .
Floor Plan Maker Online Inspirational Home Design Making Line Of Floor Plan Maker Online Awesome Floor .
Floor Plan Maker With Measurement Elegant House Blueprint Generator Stylish Inspiration Home Design .
Autocad Home Design Drawings Free Download Free Restaurant Blueprint Maker Floor Plan Plans Template Bask .
Modern Kids Room 6 .
Floor Plans Building Blueprint Maker Floor Plans House Layout Design .
3d House Maker House Creator Floor Plan Perfect Builder Free Home Design Automatic Online Best D . 3d House Maker House Design ... .
House Pro Designs Ware Inside Maker Lots Floor Survival Idea .
3d Home Design Software Free Home Plan Design Software Download Awesome Floor Plan Maker Drawing Plan .
House Blueprints Maker Free Blueprints For A House Home Design Plans Best Reverse Living House .
Free Floor Plan Designer New On Cool Maker With 3d Home Plans Rectangular Room Elegant L 0c45d662c22f32dd .
3d Floor Plan Maker Awesome Floor Plan Cad Software Inspirational Home Design Plans Home Plan Of .
... Simple Floor Plan Maker New Simple Home Designs New Www Home Design Plan New Simple Floor ... .
Simple Home Floor Plans Inspirational Image Of Simple Modern House Floor Plans Simple Modern House .
25 Home Design 3d Maker .


Krogerus is on the shortlist to receive the Financial Times and Mergermarket European M&A Award 2013 in the category of Finland Legal Advisor of the Year.
Shortlisted firms represent the top legal advisors across Europe that excelled in deal-making over the past year.
The ceremony also acknowledges top deals and financial advisors, as well as private equity, public relations and accountancy firms. The winners will be announced at a ceremony in London on 12 December.


text:
Home/Dubstep
Win Tickets To See Youngsta & Truth @ Verboten!
London-based Youngsta and New Zealand's dubstep duo Truth are respected as figureheads of dubstep culture who continue to support the scene through their…
Dub Gabriel – Raggabass Resistance (Album)
Raggabass Resistance is the 4th album release from Ragga/Dubstep producer Dub Gabriel. Three years in the making, he has been…
Mala – Radio 1 Essential Mix
Mala is a producer, DJ, label owner and visionary best known for his contributions to the scene known as dubstep.…
Noah D on BBC 1xtra (Guest Mix)
San Francisco-based electronic music composer Noah D stopped by MistaJam's Daily Dose of SMOG on BBC 1xtra for a guest…
Access All Areas 003
The 3rd show from Access All Areas, hosted by Antix and Dramatik. This time around the show features some brand…
The first show from Access All Areas hosted by Antix & Dramatik. Brand new tunes from the likes of Netsky,…


text:
Readout of President Joe Biden's Call with President López Obrador of Mexico
President Joseph R. Biden, Jr. spoke today with President Andres Manuel López Obrador of Mexico to discuss joint migration management efforts, including actions to reduce the number of individuals who unlawfully cross the U.S.-Mexico border and to expand legal pathways as an alternative to irregular migration. The two leaders also reviewed security cooperation, including joint efforts to combat the trafficking of fentanyl and new tools under the Bipartisan Safer Communities Act to prosecute gun traffickers, as well as reaffirmed their commitment to North American economic cooperation.
Joseph R. Biden, Readout of President Joe Biden's Call with President López Obrador of Mexico Online by Gerhard Peters and John T. Woolley, The American Presidency Project https://www.presidency.ucsb.edu/node/358427


As of late last week, Bishop Sharma D. Lewis was still in her hospital in Atlanta, following an additional surgery stemming from her hip replacement in late December. She would like to thank you for all your prayers, notes and cards. Her orthopedic and infectious disease doctors are continuing their treatment. She will continue to keep us updated on her progress. Please pray for Bishop Lewis' continued healing.
The 2019 Annual Conference Housing letter has been posted on the Virginia Conference website at www.vaumc.org/AC2019. The direct link to open the letter is http://www.vaumc.org/ncfilerepository/AC2019/2019HousingLetter.pdf. This letter will help you with your lodging/housing for the upcoming conference. It includes a list of area hotels and motels that have special rates for the conference. All hotels/motels holding rooms for the conference will make the special rate rooms available to conference attendees starting at 10 a.m. on Monday, Feb. 11, 2019. NO RESERVATIONS WILL BE MADE AT THE CONFERENCE RATE PRIOR TO THAT TIME. You will receive a conference packet in May that will include information on the conference including schedule, Book of Reports, name tag, registration card, and other useful information. If you need information on the conference including banquet meals before your packet arrives, please refer to the conference website at www.vaumc.org/AC2019.
1. Please join in this disaster response with a financial gift to the Virginia Conference Disaster Response Fund which will direct contributions to the crisis in Danville. All contributions can be sent to the Conference Treasurer's office with the reference "Conference Disaster Response Fund #5037".
2. Consider organizing and taking an UMVIM team to work on the homes of those affected. Please contact the Rev. Lisa Nordan sailorschk@yahoo.com or 757-513-0496.
Virginia United Methodist Pensions, Inc. (VUMPI) wants to ensure that Virginia Conference's clergy and laity have clear and accurate information regarding any potential impacts to the denominational pension plans prior to the Special Session of General Conference at the end of the month. Executive Director John Fuller shared that while the organization has received few questions regarding the various plans going before the conference and their potential impact, the organization wanted to ensure clear information is shared. You can access Wespath's Frequently Asked Questions document addressing the potential pension impact associated with the plans being considered by the Commission on a Way Forward at the following link: https://www.wespath.org/WayForwardWespathFAQ/.
"Perhaps the most important point identified in Wespath's document is that no denominational pension benefits or Wespath retirement account balances will be reduced as a result of the adoption of any of the plans being considered by the Commission on a Way Forward," said Fuller. "Wespath has worked with the Commission on several provisions to protect clergy pension benefits; in particular, the Commission is recommending (in Petition 90017) an amendment to the current denominational pension plan. Under that amendment, a clergyperson who has accrued a pension benefit and leaves The United Methodist Church would have that benefit converted to an actuarially equivalent account balance that is transferred to the clergyperson's United Methodist Personal Investment Plan (UMPIP) account. This approach is similar to the treatment of pension benefits offered by secular employers upon termination of employment." You can contact Fuller at 804-521-1128 or johnfuller@vaumc.org.
Volunteer in Mission (VIM) teams from faith communities in the Virginia Conference of The United Methodist Church who are interested in serving in mission with the Methodist Church of Haiti have a unique opportunity to receive matching funds for work in disaster response or humanitarian projects. At the hyperlinks you can access general guidelines and one-page application to request the matching funds. If any further information is needed, please contact Glenn Rowley at glennrowley@vaumc.org or 804-521-1143 or Cindy Ingroff at the Office of Justice and Missional Excellence at cindyingroff@vaumc.org or 804-521-1153.
You may also find information regarding the Haiti Matching Funds on the conference website under the Center for Justice and Missional Excellence.
3. Schedule a visit to one of our campus ministries so a personal introduction can be made.
Names and contact information for all our campus ministry directors can be found at https://vaumc.org/page.aspx?pid=3421.
Church's trust clause: Critical amid struggle?


text:
News briefs: Adena Food Rx program extended through winter
Chillicothe Gazette
Adena Food Rx program extended through winter
CHILLICOTHE – Throughout the growing season, medical providers with Adena Health System have been prescribing fresh fruit and vegetables to patients in Ross County who face barriers when it comes to eating healthy. The program, known as Food Rx, helps to eliminate the cost barriers of healthy eating and drives patients to fresh food options while supporting local farmers.
Food Rx was initially adopted by Adena providers in 2017 and traditionally ran through the summer months. However, thanks to community support, the Chillicothe Farmers Market will now be operating during the winter within the former Mount Logan Elementary School building at 841 E. Main St., Chillicothe. With this news, Adena will continue offering its Food Rx program to patients through the winter months.
"The Food Rx program has been extremely popular among Adena's patient population," said Kim Jones, director of Adena Community Health Advancement. "The availability of fresh produce in the winter months will enable our patients to continue developing those healthy eating changes they have established through the summer, making it an all year habit."
Patients who qualify, and are "prescribed" a Food Rx receive a $15 coupon to use at a participating local farmer's market. These patients include those with chronic illnesses, such as diabetes and heart disease.
In addition, people who qualify for Supplemental Nutrition Assistance Program (SNAP) are encouraged to bring their Ohio Direction Card to the Chillicothe Farmers Market. SNAP recipients are eligible to scan at least $10 in Electronic Benefit Transfer (EBT) tokens, while also receiving $10 in Fruit and Veggies Perks tokens.
Patients in Ross County who are interested in the Food Rx program should ask their primary care provider for a coupon. Those who do not have a primary care provider should call 740-779-FIND (3463) to become established with a PCP today.
The Chillicothe Famers Market, located inside the former Mount Logan Elementary School, will be open from 9 a.m. to 12 p.m. on Saturdays, November through April.
Local business selected for BBC gaming program
CHILLICOTHE - Objective Reality Games, located in Zane Plaza, is the only business in Ohio to be selected for the BBC's release of a new interactive Doctor Who experience.
Developers of the feature length virtual reality game Doctor Who: The Edge of Time are creating a custom, location based entertainment version of the game specifically for VR arcades. The UK based digital entertainment studio Maze Theory, the show's producer BBC Studios and the game's publisher PlayStack, have partnered with leading VR licensing and management companies; Synthesis VR, SpringboardVR, Ctrl V and Viveport Arcade to make the game available in more than 100 arcades worldwide.
This will give all fans and VR enthusiasts the opportunity to visit breath-taking worlds and come face to face with classic Doctor Who monsters. The game is described as a high octane and will launch on the official Doctor Who Day - Nov. 23 - and features three levels modified from the game to give players a mind-blowing 20-minute experience.
Players can travel through the iconic Doctor Who title sequence and time vortex and explore the TARDIS.


Premium 303 is the most powerful cured ink remover manufactured. It is made for all screen print applications. The ultra fast-drying formula removes ink marks quickly and does not spread. This reduces chances of ink bleeding around the cleaned area. A new "frost-reducing" additive eliminates moisture rings caused by condensation. Premium 303 removes cured plastisol ink, water and vinyl based inks, ballpoint and marker pen inks, felt-tip marks, adhesives and more. Unlike other screen print ink removers, Premium 303 does not contain any "filler" petroleum distillates.