Stable water isotopes in HadCM3: isotopic signature of El Nino-Southern Oscillation and the tropical amount effect

Stable water isotopes have been added to the full hydrological cycle of the Hadley Centre Climate model (HadCM3) coupled atmosphere-ocean GCM. Simulations of delta O-18 in precipitation and at the ocean surface compare well with observations for the present-day climate. The model has been used to investigate the isotopic anomalies associated with ENSO; it is found that the anomalous delta O-18 in precipitation is correlated with the anomalous precipitation amount in accordance with the "amount effect.'' The El Nino delta O-18 anomaly at the ocean surface is largest in coastal regions because of the mixing of ocean water and the more depleted runoff from the land surface. Coral delta O-18 anomalies were estimated, using an established empirical relationship, and generally reflect ocean surface delta O-18 anomalies in coastal regions and sea surface temperatures away from the coast. The spatial relationship between tropical precipitation and delta O-18 was investigated for the El Nino anomaly simulated by HadCM3. Weighting the El Nino precipitation anomaly by the precipitation amount at each grid box gave a large increase in the spatial correlation between tropical precipitation and delta O-18. This improvement was most apparent over land points and between 10 and 20 degrees of latitude

Machine dependence and reproducibility for coupled climate simulations: the HadGEM3-GC3.1 CMIP Preindustrial simulation

When the same weather or climate simulation is run on different high-performance computing (HPC) platforms, model outputs may not be identical for a given initial condition. While the role of HPC platforms in delivering better climate projections is to some extent discussed in the literature, attention is mainly focused on scalability and performance rather than on the impact of machine-dependent processes on the numerical solution. Here we investigate the behaviour of the Preindustrial (PI) simulation prepared by the UK Met Office for the forthcoming CMIP6 (Coupled Model Intercomparison Project Phase 6) under different computing environments. Discrepancies between the means of key climate variables were analysed at different timescales, from decadal to centennial. We found that for the two simulations to be statistically indistinguishable, a 200-year averaging period must be used for the analysis of the results. Thus, constant-forcing climate simulations using the HadGEM3-GC3.1 model are reproducible on different HPC platforms provided that a sufficiently long duration of simulation is used. In regions where El Niño–Southern Oscillation (ENSO) teleconnection patterns were detected, we found large sea surface temperature and sea ice concentration differences on centennial timescales. This indicates that a 100-year constant-forcing climate simulation may not be long enough to adequately capture the internal variability of the HadGEM3-GC3.1 model, despite this being the minimum simulation length recommended by CMIP6 protocols for many MIP (Model Intercomparison Project) experiments. On the basis of our findings, we recommend a minimum simulation length of 200 years whenever possible

Greenland deglaciation puzzles

About 23,000 years ago, the southern margins of the great Northern Hemisphere ice sheets across Europe and North America began to melt. The melt rate accelerated ∼20,000 years ago, and global sea level eventually rose by ∼130 m as meltwater flowed into the oceans. Ice cores from the Greenland and Antarctic ice sheets show the rise in atmospheric CO2 concentrations that accompanied this shift in global ice volume and climate. However, discrepancies in the temperature reconstructions from these cores have raised questions about the long-term relationship between atmospheric CO2 concentrations and Arctic temperature. On page 1177 of this issue, Buizert et al. (1) report temperature reconstructions from three locations on the Greenland ice sheet that directly address these problem

Sea ice feedbacks influence the isotopic signature of Greenland ice sheet elevation changes: Last interglacial HadCM3 simulations

Changes in the Greenland ice sheet (GIS) affect global sea level. Greenland stable water isotope (δ18O) records from ice cores offer information on past changes in the surface of the GIS. Here, we use the isotope-enabled Hadley Centre Coupled Model version 3 (HadCM3) climate model to simulate a set of last interglacial (LIG) idealised GIS surface elevation change scenarios focusing on GIS ice core sites. We investigate how δ18O depends on the magnitude and sign of GIS elevation change and evaluate how the response is altered by sea ice changes. We find that modifying GIS elevation induces changes in Northern Hemisphere atmospheric circulation, sea ice and precipitation patterns. These climate feedbacks lead to ice-core-averaged isotopic lapse rates of 0.49 ‰ (100 m)−1 for the lowered GIS states and 0.29 ‰ (100 m)−1 for the enlarged GIS states. This is lower than the spatially derived Greenland lapse rates of 0.62–0.72 ‰ (100 m)−1. These results thus suggest non-linearities in the isotope–elevation relationship and have consequences for the interpretation of past elevation and climate changes across Greenland. In particular, our results suggest that winter sea ice changes may significantly influence isotope–elevation gradients: winter sea ice effect can decrease (increase) modelled core-averaged isotopic lapse rate values by about −19 % (and +28 %) for the lowered (enlarged) GIS states, respectively. The largest influence of sea ice on δ18O changes is found in coastal regions like the Camp Century site

Comparison of the oxygen isotope signatures in speleothem records and iHadCM3 model simulations for the last millennium

Improving the understanding of changes in the mean and variability of climate variables as well as their interrelation is crucial for reliable climate change projections. Comparisons between general circulation models and paleoclimate archives using indirect proxies for temperature or precipitation have been used to test and validate the capability of climate models to represent climate changes. The oxygen isotopic ratio δ18O, a proxy for many different climate variables, is routinely measured in speleothem samples at decadal or higher resolution, and single specimens can cover full glacial–interglacial cycles. The calcium carbonate cave deposits are precisely dateable and provide well preserved (semi-)continuous albeit multivariate climate signals in the lower and mid-latitudes, where the measured δ18O in the mineral does not directly represent temperature or precipitation. Therefore, speleothems represent suitable archives to assess climate model abilities to simulate climate variability beyond the timescales covered by meteorological observations (101–102 years). Here, we present three transient isotope-enabled simulations from the Hadley Center Climate Model version 3 (iHadCM3) covering the last millennium (850–1850 CE) and compare them to a large global dataset of speleothem δ18O records from the Speleothem Isotopes Synthesis and AnaLysis (SISAL) database version 2 (Comas-Bru et al., 2020b). We systematically evaluate offsets in mean and variance of simulated δ18O and test for the main climate drivers recorded in δ18O for individual records or regions. The time-mean spatial offsets between the simulated δ18O and the speleothem data are fairly small. However, using robust filters and spectral analysis, we show that the observed archive-based variability of δ18O is lower than simulated by iHadCM3 on decadal and higher on centennial timescales. Most of this difference can likely be attributed to the records' lower temporal resolution and averaging or smoothing processes affecting the δ18O signal, e.g., through soil water residence times. Using cross-correlation analyses at site level and modeled grid-box level, we find evidence for highly variable but generally low signal-to-noise ratios in the proxy data. This points to a high influence of cave-internal processes and regional climate particularities and could suggest low regional representativity of individual sites. Long-range strong positive correlations dominate the speleothem correlation network but are much weaker in the simulation. One reason for this could lie in a lack of long-term internal climate variability in these model simulations, which could be tested by repeating similar comparisons with other isotope-enabled climate models and paleoclimate databases

Last Interglacial Arctic sea ice as simulated by the latest generation of climate models

The 16 models that simulated the Last Interglacial climate as part of the CMIP6/PMIP4 exercise consistently produce a smaller Arctic summer sea-ice area compared to the pre-industrial period, but their reduction ranges widely (28–96% of the pre-industrial area). Causes for these differences need further investigation

Impact of abrupt sea ice loss on Greenland water isotopes during the last glacial period

Greenland ice cores provide excellent evidence of past abrupt climate changes. However, there is no universally accepted theory of how and why these Dansgaard–Oeschger (DO) events occur. Several mechanisms have been proposed to explain DO events, including sea ice, ice shelf buildup, ice sheets, atmospheric circulation, and meltwater changes. DO event temperature reconstructions depend on the stable water isotope (δ18O) and nitrogen isotope measurements from Greenland ice cores: interpretation of these measurements holds the key to understanding the nature of DO events. Here, we demonstrate the primary importance of sea ice as a control on Greenland ice core δ18O: 95% of the variability in δ18O in southern Greenland is explained by DO event sea ice changes. Our suite of DO events, simulated using a general circulation model, accurately captures the amplitude of δ18O enrichment during the abrupt DO event onsets. Simulated geographical variability is broadly consistent with available ice core evidence. We find an hitherto unknown sensitivity of the δ18O paleothermometer to the magnitude of DO event temperature increase: the change in δ18O per Kelvin temperature increase reduces with DO event amplitude. We show that this effect is controlled by precipitation seasonality

Model-data comparison of Antarctic winter sea-ice extent and Southern Ocean sea-surface temperatures during Marine Isotope Stage 5e

Marine Isotope Stage (MIS) 5e (130-116 ka) represents a ‘laboratory’ for evaluating climate model performance under warmer-than-present conditions. Climate model simulations for MIS 5e have previously failed to produce Southern Ocean sea-surface temperatures (SST) and sea-ice extent reconstructed from marine sediment core proxy records. Here we compare state of the art HadGEM3 and HadCM3 simulations of Peak MIS 5e Southern Ocean summer SST and September sea-ice concentrations with the latest marine sediment core proxy data. The model outputs and proxy records show the least consistency in the regions located near the present-day Southern Ocean gyre boundaries, implying the possibility that model simulations are currently unable to fully realise changes in gyre extent and position during MIS 5e. Including Heinrich 11 meltwater forcing in Peak MIS 5e climate simulations improves the likeness to proxy data but it is clear that longer (3-4 ka) run times are required to fully test the consistency between models and data

Sea ice led to poleward-shifted winds at the Last Glacial Maximum: the influence of state dependency on CMIP5 and PMIP3 models

Latitudinal shifts in the Southern Ocean westerly wind jet could drive changes in the glacial to interglacial ocean CO2 inventory. However, whilst CMIP5 model results feature consistent future-warming jet shifts, there is considerable disagreement in deglacial-warming jet shifts. We find here that the dependence of pre-industrial (PI) to Last Glacial Maximum (LGM) jet shifts on PI jet position, or state dependency, explains less of the shifts in jet simulated by the models for the LGM compared with future-warming scenarios. State dependence is also weaker for intensity changes, compared to latitudinal shifts in the jet. Winter sea ice was considerably more extensive during the LGM. Changes in surface heat fluxes, due to this sea ice change, probably had a large impact on the jet. Models that both simulate realistically large expansions in sea ice and feature PI jets which are south of 50° S show an increase in wind speed around 55° S and can show a poleward shift in the jet between the PI and the LGM. However, models with the PI jet positioned equatorwards of around 47° S do not show this response: the sea ice edge is too far from the jet for it to respond. In models with accurately positioned PI jets, a +1° difference in the latitude of the sea ice edge tends to be associated with a −0.85° shift in the 850 hPa jet. However, it seems that around 5° of expansion of LGM sea ice is necessary to hold the jet in its PI position. Since the Gersonde et al. (2005) data support an expansion of more than 5°, this result suggests that a slight poleward shift and intensification was the most likely jet change between the PI and the LGM. Without the effect of sea ice, models simulate poleward-shifted westerlies in warming climates and equatorward-shifted westerlies in colder climates. However, the feedback of sea ice counters and reverses the equatorward trend in cooler climates so that the LGM winds were more likely to have also been shifted slightly poleward

Global reorganization of atmospheric circulation during Dansgaard-Oschger cycles

Ice core records from Greenland provide evidence for multiple abrupt warming events recurring at millennial time scales during the last glacial interval. Although climate transitions strongly resembling these Dansgaard-Oeschger (DO) transitions have been identified in several speleothem records, our understanding of the climate and ecosystem impacts of the Greenland warming events in lower latitudes remains incomplete.info:eu-repo/semantics/publishedVersio