Holocene lowering of the Laurentide ice sheet affects North Atlantic gyre circulation and climate

The Laurentide ice sheet, which covered Canada during glacial periods, had a major influence on atmospheric circulation and surface climate, but its role in climate during the early Holocene (9–7 ka), when it was thinner and confined around Hudson Bay, is unclear. It has been suggested that the demise of the ice sheet played a role in the 8.2 ka event (an abrupt 1–3 °C Northern Hemisphere cooling lasting ~ 160 years) through the influence of changing topography on atmospheric circulation. To test this hypothesis, and to investigate the broader implications of changing ice sheet topography for climate, we analyse a set of equilibrium climate simulations with ice sheet topographies taken at 500 year intervals from 9.5 to 8.0 ka. Between 9.5 and 8.0 ka, our simulations show a 2 °C cooling south of Iceland and a 1 °C warming between 40° and 50°N in the North Atlantic. These surface temperature changes are associated with a weakening of the subtropical and subpolar gyres caused by a decreasing wind stress curl over the mid-North Atlantic as the ice sheet lowers. The climate response is strongest during the period of peak ice volume change (9.5–8.5 ka), but becomes negligible after 8.5 ka. The climatic effects of the Laurentide ice sheet lowering during the Holocene are restricted to the North Atlantic sector. Thus, topographic forcing is unlikely to have played a major role in the 8.2 ka event and had only a small effect on Holocene climate change compared to the effects of changes in greenhouse gases, insolation and ice sheet meltwater

Simulating stable carbon isotopes in the ocean component of the FAMOUS general circulation model with MOSES1 (XOAVI)

Ocean circulation and the marine carbon cycle can be indirectly inferred from stable and radiogenic carbon isotope ratios (δ13C and Δ14C, respectively), measured directly in the water column, or recorded in geological archives such as sedimentary microfossils and corals. However, interpreting these records is non-trivial because they reflect a complex interplay between physical and biogeochemical processes. By directly simulating multiple isotopic tracer fields within numerical models, we can improve our understanding of the processes that control large-scale isotope distributions and interpolate the spatiotemporal gaps in both modern and palaeo datasets. We have added the stable isotope 13C to the ocean component of the FAMOUS coupled atmosphere–ocean general circulation model, which is a valuable tool for simulating complex feedbacks between different Earth system processes on decadal to multi-millennial timescales. We tested three different biological fractionation parameterisations to account for the uncertainty associated with equilibrium fractionation during photosynthesis and used sensitivity experiments to quantify the effects of fractionation during air–sea gas exchange and primary productivity on the simulated δ13CDIC distributions. Following a 10 000-year pre-industrial spin-up, we simulated the Suess effect (the isotopic imprint of anthropogenic fossil fuel burning) to assess the performance of the model in replicating modern observations. Our implementation captures the large-scale structure and range of δ13CDIC observations in the surface ocean, but the simulated values are too high at all depths, which we infer is due to biases in the biological pump. In the first instance, the new 13C tracer will therefore be useful for recalibrating both the physical and biogeochemical components of FAMOUS

Millennial‐scale climate oscillations triggered by deglacial meltwater discharge in last glacial maximum simulations

Our limited understanding of millennial-scale variability in the context of the last glacial period can be explained by the lack of a reliable modeling framework to study abrupt climate changes under realistic glacial backgrounds. In this article, we describe a new set of long-run Last Glacial Maximum experiments where such climate shifts were triggered by different snapshots of ice-sheet meltwater derived from the early stages of the last deglaciation. Depending on the location and the magnitude of the forcing, we observe three distinct dynamical regimes and highlight a subtle window of opportunity where the climate can sustain oscillations between cold and warm modes. We identify the Eurasian Arctic and Nordic Seas regions as being most sensitive to meltwater discharge in the context of switching to a cold mode, compared to freshwater fluxes from the Laurentide ice sheets. These cold climates follow a consistent pattern in temperature, sea ice, and convection, and are largely independent from freshwater release as a result of effective AMOC collapse. Warm modes, on the other hand, show more complexity in their response to the regional pattern of the meltwater input, and within them, we observe significant differences linked to the reorganization of deep water formation sites and the subpolar gyre. Broadly, the main characteristics of the oscillations, obtained under full-glacial conditions with ice-sheet reconstruction derived meltwater patterns, share similar characteristics with δ18O records of the last glacial period, although our experiment design prevents detailed conclusions from being drawn on whether these represent actual Dansgaard-Oeschger events

Holocene vegetation dynamics of circum-Arctic permafrost peatlands

Vegetation shifts in circum-Arctic permafrost peatlands drive feedbacks with important consequences for peatland carbon budgets and the extent of permafrost thaw under changing climate. Recent shrub expansion across Arctic tundra environments has led to an increase in above-ground biomass, but the long-term spatiotemporal dynamics of shrub and tree growth in circum-Arctic peatlands remain unquantified. We investigate changes in peatland vegetation composition during the Holocene using previously-published plant macrofossil records from 76 sites across the circum-Arctic permafrost zone. In particular, we assess evidence for peatland shrubification at the continental scale. We identify increasing abundance of woody vegetation in circum-Arctic peatlands from ∼8000 years BP to present, coinciding with declining herbaceous vegetation and widespread Sphagnum expansion. Ecosystem shifts varied between regions and present-day permafrost zones, with late-Holocene shrubification most pronounced where permafrost coverage is presently discontinuous and sporadic. After ∼600 years BP, we find a proliferation of non-Sphagnum mosses in Fennoscandia and across the present-day continuous permafrost zone; and rapid expansion of Sphagnum in regions of discontinuous and isolated permafrost as expected following widespread fen-bog succession, which coincided with declining woody vegetation in eastern and western Canada. Since ∼200 years BP, both shrub expansion and decline were identified at different sites across the pan-Arctic, highlighting the complex ecological responses of circum-Arctic peatlands to post-industrial climate warming and permafrost degradation. Our results suggest that shrubification of circum-Arctic peatlands has primarily occurred alongside surface drying, resulting from Holocene climate shifts, autogenic peat accumulation, and permafrost aggradation. Future shrubification of circum-Arctic peatlands under 21st century climate change will likely be spatially heterogeneous, and be most prevalent where dry microforms persist

The problem of mutually unbiased bases in dimension 6

We outline a discretization approach to determine the maximal number of mutually unbiased bases in dimension 6. We describe the basic ideas and introduce the most important definitions to tackle this famous open problem which has been open for the last 10 years. Some preliminary results are also listed

Reconstructing Holocene temperatures in time and space using paleoclimate data assimilation

Paleoclimatic records provide valuable information about Holocene climate, revealing aspects of climate variability for a multitude of sites around the world. However, such data also possess limitations. Proxy networks are spatially uneven, seasonally biased, uncertain in time, and present a variety of challenges when used in concert to illustrate the complex variations of past climate. Paleoclimatic data assimilation provides one approach to reconstructing past climate that can account for the diverse nature of proxy records while maintaining the physics-based covariance structures simulated by climate models. Here, we use paleoclimate data assimilation to create a spatially complete reconstruction of temperature over the past 12 000 years using proxy data from the Temperature 12k database and output from transient climate model simulations. Following the last glacial period, the reconstruction shows Holocene temperatures warming to a peak near 6400 years ago followed by a slow cooling toward the present day, supporting a mid-Holocene which is at least as warm as the preindustrial. Sensitivity tests show that if proxies have an overlooked summer bias, some apparent mid-Holocene warmth could actually represent summer trends rather than annual mean trends. Regardless, the potential effects of proxy seasonal biases are insufficient to align the reconstructed global mean temperature with the warming trends seen in transient model simulations

Transient climate simulations of the deglaciation 21–9 thousand years before present; PMIP4 Core experiment design and boundary conditions

The last deglaciation, which marked the transition between the last glacial and present interglacial periods, was punctuated by a series of rapid (centennial and decadal) climate changes. Numerical climate models are useful for investigating mechanisms that underpin the events, especially now that some of the complex models can be run for multiple millennia. We have set up a Paleoclimate Modelling Intercomparison Project (PMIP) working group to coordinate efforts to run transient simulations of the last deglaciation, and to facilitate the dissemination of expertise between modellers and those engaged with reconstructing the climate of the last 21 thousand years. Here, we present the design of a coordinated Core simulation over the period 21–9 thousand years before present (ka) with time varying orbital forcing, greenhouse gases, ice sheets, and other geographical changes. A choice of two ice sheet reconstructions is given, but no ice sheet or iceberg meltwater should be prescribed in the Core simulation. Additional focussed simulations will also be coordinated on an ad-hoc basis by the working group, for example to investigate the effect of ice sheet and iceberg meltwater, and the uncertainty in other forcings. Some of these focussed simulations will focus on shorter durations around specific events to allow the more computationally expensive models to take part

Simulating marine neodymium isotope distributions using Nd v1.0 coupled to the ocean component of the FAMOUS–MOSES1 climate model: sensitivities to reversible scavenging efficiency and benthic source distributions

The neodymium (Nd) isotopic composition of seawater is a widely used ocean circulation tracer. However, uncertainty in quantifying the global ocean Nd budget, particularly constraining elusive non-conservative processes, remains a major challenge. A substantial increase in modern seawater Nd measurements from the GEOTRACES programme, coupled with recent hypotheses that a seafloor-wide benthic Nd flux to the ocean may govern global Nd isotope distributions (εNd), presents an opportunity to develop a new scheme specifically designed to test these paradigms. Here, we present the implementation of Nd isotopes (143Nd and 144Nd) into the ocean component of the FAMOUS coupled atmosphere–ocean general circulation model (Nd v1.0), a tool which can be widely used for simulating complex feedbacks between different Earth system processes on decadal to multi-millennial timescales. Using an equilibrium pre-industrial simulation tuned to represent the large-scale Atlantic Ocean circulation, we perform a series of sensitivity tests evaluating the new Nd isotope scheme. We investigate how Nd source and sink and cycling parameters govern global marine εNd distributions and provide an updated compilation of 6048 Nd concentrations and 3278 εNd measurements to assess model performance. Our findings support the notions that reversible scavenging is a key process for enhancing the Atlantic–Pacific basinal εNd gradient and is capable of driving the observed increase in Nd concentration along the global circulation pathway. A benthic flux represents a major source of Nd to the deep ocean. However, model–data disparities in the North Pacific highlight that under a uniform benthic flux, the source of εNd from seafloor sediments is too non-radiogenic in our model to be able to accurately represent seawater measurements. Additionally, model–data mismatch in the northern North Atlantic alludes to the possibility of preferential contributions from “reactive” non-radiogenic detrital sediments. The new Nd isotope scheme forms an excellent tool for exploring global marine Nd cycling and the interplay between climatic and oceanographic conditions under both modern and palaeoceanographic contexts

Latest Miocene restriction of the Mediterranean Outflow Water:a perspective from the Gulf of Cádiz

The Mediterranean-Atlantic water mass exchange provides the ideal setting for deciphering the role of gateway evolution in ocean circulation. However, the dynamics of Mediterranean Outflow Water (MOW) during the closure of the Late Miocene Mediterranean-Atlantic gateways are poorly understood. Here, we define the sedimentary evolution of Neogene basins from the Gulf of Cádiz to the West Iberian margin to investigate MOW circulation during the latest Miocene. Seismic interpretation highlights a middle to upper Messinian seismic unit of transparent facies, whose base predates the onset of the Messinian salinity crisis (MSC). Its facies and distribution imply a predominantly hemipelagic environment along the Atlantic margins, suggesting an absence or intermittence of MOW preceding evaporite precipitation in the Mediterranean, simultaneous to progressive gateway restriction. The removal of MOW from the Mediterranean-Atlantic water mass exchange reorganized the Atlantic water masses and is correlated to a severe weakening of the Atlantic Meridional Overturning Circulation (AMOC) and a period of further cooling in the North Atlantic during the latest Miocene