A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

Paleoclimate data-model comparison and the role of climate forcings over the past 1500 years

The past 1500 years provide a valuable opportunity to study the response of the climate system to external forcings. However, the integration of paleoclimate proxies with climate modeling is critical to improving the understanding of climate dynamics. In this paper, a climate system model and proxy records are therefore used to study the role of natural and anthropogenic forcings in driving the global climate. The inverse and forward approaches to paleoclimate data–model comparison are applied, and sources of uncertainty are identified and discussed. In the first of two case studies, the climate model simulations are compared with multiproxy temperature reconstructions. Robust solar and volcanic signals are detected in Southern Hemisphere temperatures, with a possible volcanic signal detected in the Northern Hemisphere. The anthropogenic signal dominates during the industrial period. It is also found that seasonal and geographical biases may cause multiproxy reconstructions to overestimate the magnitude of the long-term preindustrial cooling trend. In the second case study, the model simulations are compared with a coral δ18O record from the central Pacific Ocean. It is found that greenhouse gases, solar irradiance, and volcanic eruptions all influence the mean state of the central Pacific, but there is no evidence that natural or anthropogenic forcings have any systematic impact on El Niño–Southern Oscillation. The proxy climate relationship is found to change over time, challenging the assumption of stationarity that underlies the interpretation of paleoclimate proxies. These case studies demonstrate the value of paleoclimate data–model comparison but also highlight the limitations of current techniques and demonstrate the need to develop alternative approaches.

Influence of the Central American Seaway and Drake Passage on ocean circulation and neodymium isotopes: A model study

The sensitivity of the neodymium isotopic composition (ϵNd) to tectonic rearrangements of seaways is investigated using an Earth System Model of Intermediate Complexity. The shoaling and closure of the Central American Seaway (CAS) is simulated, as well as the opening and deepening of Drake Passage (DP). Multiple series of equilibrium simulations with various intermediate depths are performed for both seaways, providing insight into ϵNd and circulation responses to progressive throughflow evolutions. Furthermore, the sensitivity of these responses to the Atlantic Meridional Overturning Circulation (AMOC) and the neodymium boundary source is examined. Modeled ϵNd changes are compared to sediment core and ferromanganese (Fe-Mn) crust data. The model results indicate that the North Atlantic ϵNd response to the CAS shoaling is highly dependent on the AMOC state, i.e., on the AMOC strength before the shoaling to shallow depths (preclosure). Three scenarios based on different AMOC forcings are discussed, of which the model-data agreement favors a shallow preclosure (Miocene) AMOC (∼6 Sv). The DP opening causes a rather complex circulation response, resulting in an initial South Atlantic ϵNd decrease preceding a larger increase. This feature may be specific to our model setup, which induces a vigorous CAS throughflow that is strongly anticorrelated to the DP throughflow. In freshwater experiments following the DP deepening, ODP Site 1090 is mainly influenced by AMOC and DP throughflow changes, while ODP Site 689 is more strongly influenced by Southern Ocean Meridional Overturning Circulation and CAS throughflow changes. The boundary source uncertainty is largest for shallow seaways and at shallow sites.

Earth system commitments due to delayed mitigation

As long as global CO₂ emissions continue to increase annually, long-term committed Earth system changes grow much faster than current observations. A novel metric linking this future growth to policy decisions today is the mitigation delay sensitivity (MDS), but MDS estimates for Earth system variables other than peak temperature (ΔT max) are missing. Using an Earth System Model of Intermediate Complexity, we show that the current emission increase rate causes a ΔT max increase roughly 3–7.5 times as fast as observed warming, and a millenial steric sea level rise (SSLR) 7–25 times as fast as observed SSLR, depending on the achievable rate of emission reductions after the peak of emissions. These ranges are only slightly affected by the uncertainty range in equilibrium climate sensitivity, which is included in the above values. The extent of ocean acidification at the end of the century is also strongly dependent on the starting time and rate of emission reductions. The preservable surface ocean area with sufficient aragonite supersaturation for coral reef growth is diminished globally at an MDS of roughly 25%–80% per decade. A near-complete loss of this area becomes unavoidable if mitigation is delayed for a few years to decades. Also with respect to aragonite, 12%–18% of the Southern Ocean surface become undersaturated per decade, if emission reductions are delayed beyond 2015–2040. We conclude that the consequences of delaying global emission reductions are much better captured if the MDS of relevant Earth system variables is communicated in addition to current trends and total projected future changes.

Past and future evolution of Abies alba forests in Europe - comparison of a dynamic vegetation model with palaeo data and observations

Information on how species distributions and ecosystem services are impacted by anthropogenic climate change is important for adaptation planning. Palaeo data suggest that Abies alba formed forests under significantly warmer-than-present conditions in Europe and might be a native substitute for widespread drought-sensitive temperate and boreal tree species such as beech (Fagus sylvatica) and spruce (Picea abies) under future global warming conditions. Here, we combine pollen and macrofossil data, modern observations, and results from transient simulations with the LPX-Bern dynamic global vegetation model to assess past and future distributions of A. alba in Europe. LPX-Bern is forced with climate anomalies from a run over the past 21 000 years with the Community Earth System Model, modern climatology, and with 21st-century multimodel ensemble results for the high-emission RCP8.5 and the stringent mitigation RCP2.6 pathway. The simulated distribution for present climate encompasses the modern range of A. alba, with the model exceeding the present distribution in north-western and southern Europe. Mid-Holocene pollen data and model results agree for southern Europe, suggesting that at present, human impacts suppress the distribution in southern Europe. Pollen and model results both show range expansion starting during the Bølling–Allerød warm period, interrupted by the Younger Dryas cold, and resuming during the Holocene. The distribution of A. alba expands to the north-east in all future scenarios, whereas the potential (currently unrealized) range would be substantially reduced in southern Europe under RCP8.5. A. alba maintains its current range in central Europe despite competition by other thermophilous tree species. Our combined palaeoecological and model evidence suggest that A. alba may ensure important ecosystem services including stand and slope stability, infrastructure protection, and carbon sequestration under significantly warmer-than-present conditions in central Europe.

Variations of the Atlantic meridional overturning circulation in control and transient simulations of the last millennium

The variability of the Atlantic meridional overturing circulation (AMOC) strength is investigated in control experiments and in transient simulations of up to the last millennium using the low-resolution Community Climate System Model version 3. In the transient simulations the AMOC exhibits enhanced low-frequency variability that is mainly caused by infrequent transitions between two semi-stable circulation states which amount to a 10 percent change of the maximum overturning. One transition is also found in a control experiment, but the time-varying external forcing significantly increases the probability of the occurrence of such events though not having a direct, linear impact on the AMOC. The transition from a high to a low AMOC state starts with a reduction of the convection in the Labrador and Irminger Seas and goes along with a changed barotropic circulation of both gyres in the North Atlantic and a gradual strengthening of the convection in the Greenland-Iceland-Norwegian (GIN) Seas. In contrast, the transition from a weak to a strong overturning is induced by decreased mixing in the GIN Seas. As a consequence of the transition, regional sea surface temperature (SST) anomalies are found in the midlatitude North Atlantic and in the convection regions with an amplitude of up to 3 K. The atmospheric response to the SST forcing associated with the transition indicates a significant impact on the Scandinavian surface air temperature (SAT) in the order of 1 K. Thus, the changes of the ocean circulation make a major contribution to the Scandinavian SAT variability in the last millennium.

Comparison of climate field reconstruction techniques: Application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.

Atmospheric CO₂ response to volcanic eruptions: The role of ENSO, season, and variability

Tropical explosive volcanism is one of the most important natural factors that significantly impact the climate system and the carbon cycle on annual to multi-decadal time scales. The three largest explosive eruptions in the last 50�years�Agung, El Chichón, and Pinatubo�occurred in spring/summer in conjunction with El Niño events and left distinct negative signals in the observational temperature and CO2 records. However, confounding factors such as seasonal variability and El Niño-Southern Oscillation (ENSO) may obscure the forcing-response relationship. We determine for the first time the extent to which initial conditions, i.e., season and phase of the ENSO, and internal variability influence the coupled climate and carbon cycle response to volcanic forcing and how this affects estimates of the terrestrial and oceanic carbon sinks. Ensemble simulations with the Earth System Model (Climate System Model 1.4-carbon) predict that the atmospheric CO2 response is �60 larger when a volcanic eruption occurs during El Niño and in winter than during La Niña conditions. Our simulations suggest that the Pinatubo eruption contributed 11�±�6 to the 25�Pg terrestrial carbon sink inferred over the decade 1990�1999 and �2�±�1 to the 22�Pg oceanic carbon sink. In contrast to recent claims, trends in the airborne fraction of anthropogenic carbon cannot be detected when accounting for the decadal-scale influence of explosive volcanism and related uncertainties. Our results highlight the importance of considering the role of natural variability in the carbon cycle for interpretation of observations and for data-model intercomparison.

Allowable carbon emissions lowered by multiple climate targets

Climate targets are designed to inform policies that would limit the magnitude and impacts of climate change caused by anthropogenic emissions of greenhouse gases and other substances. The target that is currently recognized by most world governments1 places a limit of two degrees Celsius on the global mean warming since preindustrial times. This would require large sustained reductions in carbon dioxide emissions during the twenty-first century and beyond2, 3, 4. Such a global temperature target, however, is not sufficient to control many other quantities, such as transient sea level rise5, ocean acidification6, 7 and net primary production on land8, 9. Here, using an Earth system model of intermediate complexity (EMIC) in an observation-informed Bayesian approach, we show that allowable carbon emissions are substantially reduced when multiple climate targets are set. We take into account uncertainties in physical and carbon cycle model parameters, radiative efficiencies10, climate sensitivity11 and carbon cycle feedbacks12, 13 along with a large set of observational constraints. Within this framework, we explore a broad range of economically feasible greenhouse gas scenarios from the integrated assessment community14, 15, 16, 17 to determine the likelihood of meeting a combination of specific global and regional targets under various assumptions. For any given likelihood of meeting a set of such targets, the allowable cumulative emissions are greatly reduced from those inferred from the temperature target alone. Therefore, temperature targets alone are unable to comprehensively limit the risks from anthropogenic emissions.

Northern Hemisphere temperature reconstruction during the last millennium using multiple annual proxies

Previous studies have either exclusively used annual tree-ring data or have combined tree-ring series with other, lower temporal resolution proxy series. Both approaches can lead to significant uncertainties, as tree-rings may underestimate the amplitude of past temperature variations, and the validity of non-annual records cannot be clearly assessed. In this study, we assembled 45 published Northern Hemisphere (NH) temperature proxy records covering the past millennium, each of which satisfied 3 essential criteria: the series must be of annual resolution, span at least a thousand years, and represent an explicit temperature signal. Suitable climate archives included ice cores, varved lake sediments, tree-rings and speleothems. We reconstructed the average annual land temperature series for the NH over the last millennium by applying 3 different reconstruction techniques: (1) principal components (PC) plus second-order autoregressive model (AR2), (2) composite plus scale (CPS) and (3) regularized errors-in-variables approach (EIV). Our reconstruction is in excellent agreement with 6 climate model simulations (including the first 5 models derived from the fifth phase of the Coupled Model Intercomparison Project (CMIP5) and an earth system model of intermediate complexity (LOVECLIM), showing similar temperatures at multi-decadal timescales; however, all simulations appear to underestimate the temperature during the Medieval Warm Period (MWP). A comparison with other NH reconstructions shows that our results are consistent with earlier studies. These results indicate that well-validated annual proxy series should be used to minimize proxy-based artifacts, and that these proxy series contain sufficient information to reconstruct the low-frequency climate variability over the past millennium.

Ocean (De)oxygenation Across the Last Deglaciation: Insights for the Future

Anthropogenic warming is expected to drive oxygen out of the ocean as the water temperature rises and the rate of exchange between subsurface waters and the atmosphere slows due to enhanced upper ocean density stratification. Observations from recent decades are tantalizingly consistent with this prediction, though these changes remain subtle in the face of natural variability. Earth system model projections unanimously predict a long-term decrease in the global ocean oxygen inventory, but show regional discrepancies, particularly in the most oxygen-depleted waters, owing to the complex interplay between oxygen supply pathways and oxygen consumption. The geological record provides an orthogonal perspective, showing how the oceanic oxygen content varied in response to prior episodes of climate change. These past changes were much slower than the current, anthropogenic change, but can help to appraise sensitivities, and point toward potentially dominant mechanisms of change. Consistent with the model projections, marine sediments recorded an overall expansion of low-oxygen waters in the upper ocean as it warmed at the end of the last ice age. This expansion was not linearly related with temperature, though, but reached a deoxygenation extreme midway through the warming. Meanwhile, the deep ocean became better oxygenated, opposite the general expectation. These observations require that significant changes in apparent oxygen utilization occurred, suggesting that they will also be important in the future.

Multiyear predictability of tropical marine productivity

With the emergence of decadal predictability simulations, research toward forecasting variations of the climate system now covers a large range of timescales. However, assessment of the capacity to predict natural variations of relevant biogeochemical variables like carbon fluxes, pH, or marine primary productivity remains unexplored. Among these, the net primary productivity (NPP) is of particular relevance in a forecasting perspective. Indeed, in regions like the tropical Pacific (30°N–30°S), NPP exhibits natural fluctuations at interannual to decadal timescales that have large impacts on marine ecosystems and fisheries. Here, we investigate predictions of NPP variations over the last decades (i.e., from 1997 to 2011) with an Earth system model within the tropical Pacific. Results suggest a predictive skill for NPP of 3 y, which is higher than that of sea surface temperature (1 y). We attribute the higher predictability of NPP to the poleward advection of nutrient anomalies (nitrate and iron), which sustain fluctuations in phytoplankton productivity over several years. These results open previously unidentified perspectives to the development of science-based management approaches to marine resources relying on integrated physical-biogeochemical forecasting systems.

Synoptic drivers of 400 years of summer temperature and precipitation variability on Mt. Olympus, Greece

The Mediterranean region has been identified as a global warming hotspot, where future climate impacts are expected to have significant consequences on societal and ecosystem well-being. To put ongoing trends of summer climate into the context of past natural variability, we reconstructed climate from maximum latewood density (MXD) measurements of Pinus heldreichii (1521–2010) and latewood width (LWW) of Pinus nigra (1617–2010) on Mt. Olympus, Greece. Previous research in the northeastern Mediterranean has primarily focused on inter-annual variability, omitting any low-frequency trends. The present study utilizes methods capable of retaining climatically driven long-term behavior of tree growth. The LWW chronology corresponds closely to early summer moisture variability (May–July, r = 0.65, p < 0.001, 1950–2010), whereas the MXD-chronology relates mainly to late summer warmth (July–September, r = 0.64, p < 0.001; 1899–2010). The chronologies show opposing patterns of decadal variability over the twentieth century (r = −0.68, p < 0.001) and confirm the importance of the summer North Atlantic Oscillation (sNAO) for summer climate in the northeastern Mediterranean, with positive sNAO phases inducing cold anomalies and enhanced cloudiness and precipitation. The combined reconstructions document the late twentieth—early twenty-first century warming and drying trend, but indicate generally drier early summer and cooler late summer conditions in the period ~1700–1900 CE. Our findings suggest a potential decoupling between twentieth century atmospheric circulation patterns and pre-industrial climate variability. Furthermore, the range of natural climate variability stretches beyond summer moisture availabilityobserved in recent decades and thus lends credibility to the significant drying trends projected for this region in current Earth System Model simulations.

Impact of a Reduced Arctic Sea Ice Cover on Ocean and Atmospheric Properties

The Arctic sea ice cover declined over the last few decades and reached a record minimum in 2007, with a slight recovery thereafter. Inspired by this the authors investigate the response of atmospheric and oceanic properties to a 1-yr period of reduced sea ice cover. Two ensembles of equilibrium and transient simulations are produced with the Community Climate System Model. A sea ice change is induced through an albedo change of 1 yr. The sea ice area and thickness recover in both ensembles after 3 and 5 yr, respectively. The sea ice anomaly leads to changes in ocean temperature and salinity to a depth of about 200 m in the Arctic Basin. Further, the salinity and temperature changes in the surface layer trigger a “Great Salinity Anomaly” in the North Atlantic that takes roughly 8 yr to travel across the North Atlantic back to high latitudes. In the atmosphere the changes induced by the sea ice anomaly do not last as long as in the ocean. The response in the transient and equilibrium simulations, while similar overall, differs in specific regional and temporal details. The surface air temperature increases over the Arctic Basin and the anomaly extends through the whole atmospheric column, changing the geopotential height fields and thus the storm tracks. The patterns of warming and thus the position of the geopotential height changes vary in the two ensembles. While the equilibrium simulation shifts the storm tracks to the south over the eastern North Atlantic and Europe, the transient simulation shifts the storm tracks south over the western North Atlantic and North America. The authors propose that the overall reduction in sea ice cover is important for producing ocean anomalies; however, for atmospheric anomalies the regional location of the sea ice anomalies is more important. While observed trends in Arctic sea ice are large and exceed those simulated by comprehensive climate models, there is little evidence based on this particular model that the seasonal loss of sea ice (e.g., as occurred in 2007) would constitute a threshold after which the Arctic would exhibit nonlinear, irreversible, or strongly accelerated sea ice loss. Caution should be exerted when extrapolating short-term trends to future sea ice behavior.