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

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

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

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

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

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

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

Composite data set of last glacial Dansgaard/Oeschger events obtained from stable oxygen isotopes in speleothems

Millennial scale climate variations called Dansgaard-Oeschger cycles occurred frequently during the last glacial, with their central impact on climate in the North Atlantic region. These events are, for example, well captured by the stable oxygen isotope composition in continental ice from Greenland, but also in records from other regions. Recently, it has been shown that a water isotope enabled general circulation model is able to reproduce those millennial-scale oxygen isotope changes from Greenland (Sime et al., 2019). On a global scale, this isotope-enabled model has not been tested in its performance, as stable oxygen isotope records covering those millennial scale variability were so far missing or not systematically compiled.info:eu-repo/semantics/publishedVersio

CMIP6/PMIP4 simulations of the mid-Holocene and Last Interglacial using HadGEM3: comparison to the pre-industrial era, previous model versions and proxy data

Results for these simulations are assessed firstly against the same model's pre-industrial control simulation (a simulation comparison, to describe and understand the differences between the pre-industrial – PI – and the two palaeo simulations) and secondly against previous versions of the same model relative to newly available proxy data (a model–data comparison, to compare all available simulations from the same model with proxy data to assess any improvements due to model advances). The introduction of this newly available proxy data adds further novelty to this study. Globally, for metrics such as 1.5 m temperature and surface rainfall, whilst both the recent palaeoclimate simulations are mostly capturing the expected sign and, in some places, magnitude of change relative to the pre-industrial, this is geographically and seasonally dependent. Compared to newly available proxy data (including sea surface temperature – SST – and rainfall) and also incorporating data from previous versions of the model shows that the relative accuracy of the simulations appears to vary according to metric, proxy reconstruction used for comparison and geographical location. In some instances, such as mean rainfall in the mid-Holocene, there is a clear and linear improvement, relative to proxy data, from the oldest to the newest generation of the model. When zooming into northern Africa, a region known to be problematic for models in terms of rainfall enhancement, the behaviour of the West African monsoon in both recent palaeoclimate simulations is consistent with current understanding, suggesting a wetter monsoon during the mid-Holocene and (more so) the Last Interglacial, relative to the pre-industrial era. However, regarding the well-documented “Saharan greening” during the mid-Holocene, results here suggest that the most recent version of the UK's physical model is still unable to reproduce the increases suggested by proxy data, consistent with all other previous models to date

Interconnectivity between volume transports through Arctic straits

Arctic heat and freshwater budgets are highly sensitive to volume transports through the Arctic‐Subarctic straits. Here we study the interconnectivity of volume transports through Arctic straits in three models; two coupled global climate models, one with a third‐degree horizontal ocean resolution (HiGEM1.1) and one with a twelfth‐degree horizontal ocean resolution (HadGEM3), and one ocean‐only model with an idealized polar basin (tenth‐degree horizontal resolution). The two global climate models indicate that there is a strong anti‐correlation between the Bering Strait throughflow and the transport through the Nordic Seas, a second strong anti‐correlation between the transport through the Canadian Artic Archipelago (CAA) and the Nordic Seas transport, and a third strong anti‐correlation is found between the Fram Strait and the Barents Sea throughflows. We find that part of the strait correlations is due to the strait transports being coincidentally driven by large‐scale atmospheric forcing patterns. However, there is also a role for fast wave adjustments of some straits flows to perturbations in other straits since atmospheric forcing of individual strait flows alone cannot lead to near mass balance fortuitously every year. Idealized experiments with an ocean model (NEMO3.6) that investigate such causal strait relations suggest that perturbations in the Bering Strait are compensated preferentially in the Fram Strait due to the narrowness of the western Arctic shelf and the deeper depth of the Fram Strait