Search

The non-local model of mixing based on internal wave breaking, IDEMIX, is implemented as an enhancement of a turbulent kinetic energy closure model in three non-eddy resolving general circulation ocean models that differ in the discretization and choice of computational grids. In IDEMIX internal wave energy is generated by an energy flux resulting from near-inertial waves induced by wind forcing at the surface, and at the bottom, by an energy flux that parameterizes the transfer of energy between baroclinic and barotropic tides. In all model simulations with IDEMIX, the mixing work is increased compared to the reference solutions without IDEMIX, reaching values in better agreement with finestructure observations. Furthermore, the horizontal structure of the mixing work is more realistic as a consequence of the heterogeneous forcing functions. All models with IDEMIX simulate deeper thermocline depths related to stronger shallow overturning cells in the Indo-Pacific. In the North Atlantic, deeper mixed layers in simulations with IDEMIX are associated with an increased Atlantic overturning circulation and an increase of northward heat transports toward more realistic values. The response of the deep Indo-Pacific overturning circulation and the weak bottom cell of the Atlantic to the inclusion of IDEMIX is incoherent between the models, suggesting that additional unidentified processes and numerical mixing may confound the analysis. Applying different tidal forcing functions leads to simulation differences that are small compared to differences between the different models or between simulations with IDEMIX and without IDEMIX.

Existing global mean surface temperature reconstructions for the Holocene lack high-frequency variability that is essential for contextualising recent trends and extremes in the Earth's climate system. Here, we isolate and recombine archive-specific climate signals to generate a frequency-optimised record of interannual to multi-millennial temperature changes for the past 12 000 years. Average temperatures before ∼8000 years BP and after ∼4000 years BP were 0.26 (±2.84) °C and 0.07 (±2.11) °C cooler than the long-term mean (0–12 000 years BP), while the Holocene Climate Optimum ∼7000–4000 years BP was 0.40 (±1.86) °C warmer. Biased towards Northern Hemisphere summer temperatures, our multi-proxy record captures the spectral properties of transient Earth system model simulations for the same spatial and season domain. The new frequency-optimised trajectory emphasises the importance and complex interplay of natural climate forcing factors throughout the Holocene, with an approximation of the full range of past temperature changes providing novel insights for policymakers addressing the risks of recent anthropogenic warming.

State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

The Atlantic Meridional Overturning Circulation (AMOC) is a key feature of the North Atlantic with global ocean impacts. The AMOC's response to past changes in forcings during the Holocene provides important context for the coming centuries. Here, we investigate AMOC trends using an emerging set of transient simulations using multiple global climate models for the past 6,000 years. Although some models show changes, no consistent trend in overall AMOC strength during the mid-to-late Holocene emerges from the ensemble. We interpret this result to suggest no overall change in AMOC, which fits with our assessment of available proxy reconstructions. The decadal variability of the AMOC does not change in ensemble during the mid- and late-Holocene. There are interesting AMOC changes seen in the early Holocene, but their nature depends a lot on which inputs are used to drive the experiment.

The Atlantic Meridional Overturning Circulation (AMOC) is a key feature of the North Atlantic with global ocean impacts. The AMOC's response to past changes in forcings during the Holocene provides important context for the coming centuries. Here, we investigate AMOC trends using an emerging set of transient simulations using multiple global climate models for the past 6,000 years. Although some models show changes, no consistent trend in overall AMOC strength during the mid-to-late Holocene emerges from the ensemble. We interpret this result to suggest no overall change in AMOC, which fits with our assessment of available proxy reconstructions. The decadal variability of the AMOC does not change in ensemble during the mid- and late-Holocene. There are interesting AMOC changes seen in the early Holocene, but their nature depends a lot on which inputs are used to drive the experiment.

State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

The land is a pivotal component of the Earth and its climate system since many processes of natural variations in the climate system, which affect the environment and human society, are governed by the land surface. Hence, a good representation of the thermal and hydrological states of the land surface in climate models is important to have a realistic simulation of the coupling between the atmosphere and the lito-biosphere. An influencing factor for improving the realism of the ground energy and water balance in climate models is the depth of the land zero-flux Bottom Boundary Condition Placement (BBCP). Despite recent improvements in modeling land surface processes in climate models, only limited attention has been directed toward the effect of the BBCP in Land Surface Models (LSMs) and its impact on the representation of terrestrial thermodynamics. Previous analytical and modeling studies suggest that the simulation of subsurface thermodynamics in current-generation climate models is not accurate due to the zero-heat-flux BBCP being imposed too close to the surface. An insufficiently deep land component in current-generation climate models compromises the simulation of the terrestrial thermal state and can influence land-atmosphere interactions. Further improvements in LSMs relate to the representation and sensitivity of coupling processes between the ground thermodynamic and hydrological regimes. As moisture is one of the main drivers of near-surface climate interactions, the hydro-thermodynamic coupling is crucial for studying the impacts of perturbations caused by human activity. Under climate change conditions, some areas and ecosystems are more vulnerable to a rapidly warming world than others..., El suelo es un componente fundamental de la Tierra y su sistema climático, ya que muchos procesos de variabilidad natural en el sistema climático, que afectan al medio ambiente y la sociedad humana, se rigen por lo que ocurre en la superficie terrestre. Por lo tanto, una buena representación de los estados térmicos e hidrológicos de la superficie terrestre en los modelos climáticos es importante para tener una simulación realista del acoplamiento entre la atmósfera y la lito-biosfera. Un factor que influye en la mejora del realismo de la energía del suelo y el equilibrio hídrico en los modelos climáticos es la profundidad de la posición de la condición de contorno del fondo (BBCP) de la tierra. Los estudios anteriores sugieren que la simulación dela termodinámica del subsuelo en los modelos climáticos de la generación actual no es precisa debido a que el BBCP de flujo de calor cero se impone demasiado cerca de la superficie. Una componente terrestre insuficientemente profunda en los modelos climáticos compromete la simulación del estado térmico terrestre y puede influir en las interacciones tierra-atmósfera. Otras mejoras en los LSM se relacionan con la representación y la sensibilidad de los procesos de acoplamiento entre los regímenes hidrológico y termodinámico del suelo. Dado que la humedad es uno delos principales impulsores de las interacciones climáticas cercanas a la superficie, el acoplamiento hidro-termodinámico es crucial para estudiar los impactos de las perturbaciones causadas por la actividad humana. En condiciones de cambio climático, algunos hábitats y ecosistemas son más vulnerables que otros en un mundo que se calienta rápidamente en un contexto de vertiginoso cambio global...

Numerical modeling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison are from the same position along the Earth's orbit. Most models use the classical calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, Last Interglacial, and pre-industrial climate conditions. We came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the classical seasonal cycle is applied. (d) Finally, month-length adjusted values f

Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi-centennial variability. However, all the simulations without this

Numerical modeling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison are from the same position along the Earth's orbit. Most models use the classical calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, Last Interglacial, and pre-industrial climate conditions. We came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the classical seasonal cycle is applied. (d) Finally, month-length adjusted values f

Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi-centennial variability. However, all the simulations without this

Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique oppor-tunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120-130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120-150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemi-sphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi -centennial variability. However, all the simulations without th

Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique oppor-tunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120-130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120-150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemi-sphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi -centennial variability. However, all the simulations without th

Northern Hemisphere (NH) climate has experienced various coherent wintertime multidecadal climate trends in stratosphere, troposphere, ocean, and cryosphere. However, the overall mechanistic framework linking these trends is not well established. Here we show, using long-term transient forced coupled climate simulation, that large parts of the coherent NH-multidecadal changes can be understood within a damped coupled stratosphere/troposphere/ocean-oscillation framework. Wave-induced downward propagating positive stratosphere/troposphere-coupled Northern Annular Mode (NAM) and associated stratospheric cooling initiate delayed thermohaline strengthening of Atlantic overturning circulation and extratropical Atlantic-gyres. These increase the poleward oceanic heat transport leading to Arctic sea-ice melting, Arctic warming amplification, and large-scale Atlantic warming, which in turn initiates wave-induced downward propagating negative NAM and stratospheric warming and therefore reverse the oscillation phase. This coupled variability improves the performance of statistical models, which project further weakening of North Atlantic Oscillation, North Atlantic cooling and hiatus in wintertime North Atlantic-Arctic sea-ice and global surface temperature just like the 1950s-1970s.

Numerical modeling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison are from the same position along the Earth's orbit. Most models use the classical calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, Last Interglacial, and pre-industrial climate conditions. We came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the classical seasonal cycle is applied. (d) Finally, month-length adjusted values f

The land is a pivotal component of the Earth and its climate system since many processes of natural variations in the climate system, which affect the environment and human society, are governed by the land surface. Hence, a good representation of the thermal and hydrological states of the land surface in climate models is important to have a realistic simulation of the coupling between the atmosphere and the lito-biosphere. An influencing factor for improving the realism of the ground energy and water balance in climate models is the depth of the land zero-flux Bottom Boundary Condition Placement (BBCP). Despite recent improvements in modeling land surface processes in climate models, only limited attention has been directed toward the effect of the BBCP in Land Surface Models (LSMs) and its impact on the representation of terrestrial thermodynamics. Previous analytical and modeling studies suggest that the simulation of subsurface thermodynamics in current-generation climate models is not accurate due to the zero-heat-flux BBCP being imposed too close to the surface. An insufficiently deep land component in current-generation climate models compromises the simulation of the terrestrial thermal state and can influence land-atmosphere interactions. Further improvements in LSMs relate to the representation and sensitivity of coupling processes between the ground thermodynamic and hydrological regimes. As moisture is one of the main drivers of near-surface climate interactions, the hydro-thermodynamic coupling is crucial for studying the impacts of perturbations caused by human activity. Under climate change conditions, some areas and ecosystems are more vulnerable to a rapidly warming world than others..., El suelo es un componente fundamental de la Tierra y su sistema climático, ya que muchos procesos de variabilidad natural en el sistema climático, que afectan al medio ambiente y la sociedad humana, se rigen por lo que ocurre en la superficie terrestre. Por lo tanto, una buena representación de los estados térmicos e hidrológicos de la superficie terrestre en los modelos climáticos es importante para tener una simulación realista del acoplamiento entre la atmósfera y la lito-biosfera. Un factor que influye en la mejora del realismo de la energía del suelo y el equilibrio hídrico en los modelos climáticos es la profundidad de la posición de la condición de contorno del fondo (BBCP) de la tierra. Los estudios anteriores sugieren que la simulación dela termodinámica del subsuelo en los modelos climáticos de la generación actual no es precisa debido a que el BBCP de flujo de calor cero se impone demasiado cerca de la superficie. Una componente terrestre insuficientemente profunda en los modelos climáticos compromete la simulación del estado térmico terrestre y puede influir en las interacciones tierra-atmósfera. Otras mejoras en los LSM se relacionan con la representación y la sensibilidad de los procesos de acoplamiento entre los regímenes hidrológico y termodinámico del suelo. Dado que la humedad es uno delos principales impulsores de las interacciones climáticas cercanas a la superficie, el acoplamiento hidro-termodinámico es crucial para estudiar los impactos de las perturbaciones causadas por la actividad humana. En condiciones de cambio climático, algunos hábitats y ecosistemas son más vulnerables que otros en un mundo que se calienta rápidamente en un contexto de vertiginoso cambio global...

Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique oppor-tunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120-130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120-150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemi-sphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi -centennial variability. However, all the simulations without th

Highlights: • Multi-centennial oscillation with 100–200 years periods is evident in proxy data and model simulations during the Holocene. • Multi-centennial oscillation is a global signal and is more significant in the Northern Hemisphere high latitudes. • None of the external forcings is found to be the sole driver of the multi-centennial variability. • It indicates the multi-centennial oscillation may be due to potential internal drivers and essential feedbacks. Abstract: Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combi

Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique oppor-tunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120-130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120-150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemi-sphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi -centennial variability. However, all the simulations without th

The land is a pivotal component of the Earth and its climate system since many processes of natural variations in the climate system, which affect the environment and human society, are governed by the land surface. Hence, a good representation of the thermal and hydrological states of the land surface in climate models is important to have a realistic simulation of the coupling between the atmosphere and the lito-biosphere. An influencing factor for improving the realism of the ground energy and water balance in climate models is the depth of the land zero-flux Bottom Boundary Condition Placement (BBCP). Despite recent improvements in modeling land surface processes in climate models, only limited attention has been directed toward the effect of the BBCP in Land Surface Models (LSMs) and its impact on the representation of terrestrial thermodynamics. Previous analytical and modeling studies suggest that the simulation of subsurface thermodynamics in current-generation climate models is not accurate due to the zero-heat-flux BBCP being imposed too close to the surface. An insufficiently deep land component in current-generation climate models compromises the simulation of the terrestrial thermal state and can influence land-atmosphere interactions. Further improvements in LSMs relate to the representation and sensitivity of coupling processes between the ground thermodynamic and hydrological regimes. As moisture is one of the main drivers of near-surface climate interactions, the hydro-thermodynamic coupling is crucial for studying the impacts of perturbations caused by human activity. Under climate change conditions, some areas and ecosystems are more vulnerable to a rapidly warming world than others..., El suelo es un componente fundamental de la Tierra y su sistema climático, ya que muchos procesos de variabilidad natural en el sistema climático, que afectan al medio ambiente y la sociedad humana, se rigen por lo que ocurre en la superficie terrestre. Por lo tanto, una buena representación de los estados térmicos e hidrológicos de la superficie terrestre en los modelos climáticos es importante para tener una simulación realista del acoplamiento entre la atmósfera y la lito-biosfera. Un factor que influye en la mejora del realismo de la energía del suelo y el equilibrio hídrico en los modelos climáticos es la profundidad de la posición de la condición de contorno del fondo (BBCP) de la tierra. Los estudios anteriores sugieren que la simulación dela termodinámica del subsuelo en los modelos climáticos de la generación actual no es precisa debido a que el BBCP de flujo de calor cero se impone demasiado cerca de la superficie. Una componente terrestre insuficientemente profunda en los modelos climáticos compromete la simulación del estado térmico terrestre y puede influir en las interacciones tierra-atmósfera. Otras mejoras en los LSM se relacionan con la representación y la sensibilidad de los procesos de acoplamiento entre los regímenes hidrológico y termodinámico del suelo. Dado que la humedad es uno delos principales impulsores de las interacciones climáticas cercanas a la superficie, el acoplamiento hidro-termodinámico es crucial para estudiar los impactos de las perturbaciones causadas por la actividad humana. En condiciones de cambio climático, algunos hábitats y ecosistemas son más vulnerables que otros en un mundo que se calienta rápidamente en un contexto de vertiginoso cambio global...

Highlights: • Multi-centennial oscillation with 100–200 years periods is evident in proxy data and model simulations during the Holocene. • Multi-centennial oscillation is a global signal and is more significant in the Northern Hemisphere high latitudes. • None of the external forcings is found to be the sole driver of the multi-centennial variability. • It indicates the multi-centennial oscillation may be due to potential internal drivers and essential feedbacks. Abstract: Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combi

Highlights: • Multi-centennial oscillation with 100–200 years periods is evident in proxy data and model simulations during the Holocene. • Multi-centennial oscillation is a global signal and is more significant in the Northern Hemisphere high latitudes. • None of the external forcings is found to be the sole driver of the multi-centennial variability. • It indicates the multi-centennial oscillation may be due to potential internal drivers and essential feedbacks. Abstract: Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combi

Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique oppor-tunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120-130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120-150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemi-sphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi -centennial variability. However, all the simulations without th

Northern Hemisphere (NH) climate has experienced various coherent wintertime multidecadal climate trends in stratosphere, troposphere, ocean, and cryosphere. However, the overall mechanistic framework linking these trends is not well established. Here we show, using long-term transient forced coupled climate simulation, that large parts of the coherent NH-multidecadal changes can be understood within a damped coupled stratosphere/troposphere/ocean-oscillation framework. Wave-induced downward propagating positive stratosphere/troposphere-coupled Northern Annular Mode (NAM) and associated stratospheric cooling initiate delayed thermohaline strengthening of Atlantic overturning circulation and extratropical Atlantic-gyres. These increase the poleward oceanic heat transport leading to Arctic sea-ice melting, Arctic warming amplification, and large-scale Atlantic warming, which in turn initiates wave-induced downward propagating negative NAM and stratospheric warming and therefore reverse the oscillation phase. This coupled variability improves the performance of statistical models, which project further weakening of North Atlantic Oscillation, North Atlantic cooling and hiatus in wintertime North Atlantic-Arctic sea-ice and global surface temperature just like the 1950s-1970s.

Numerical modeling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison are from the same position along the Earth's orbit. Most models use the classical calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, Last Interglacial, and pre-industrial climate conditions. We came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the classical seasonal cycle is applied. (d) Finally, month-length adjusted values f

Numerical modeling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison are from the same position along the Earth's orbit. Most models use the classical calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, Last Interglacial, and pre-industrial climate conditions. We came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the classical seasonal cycle is applied. (d) Finally, month-length adjusted values f

Northern Hemisphere (NH) climate has experienced various coherent wintertime multidecadal climate trends in stratosphere, troposphere, ocean, and cryosphere. However, the overall mechanistic framework linking these trends is not well established. Here we show, using long-term transient forced coupled climate simulation, that large parts of the coherent NH-multidecadal changes can be understood within a damped coupled stratosphere/troposphere/ocean-oscillation framework. Wave-induced downward propagating positive stratosphere/troposphere-coupled Northern Annular Mode (NAM) and associated stratospheric cooling initiate delayed thermohaline strengthening of Atlantic overturning circulation and extratropical Atlantic-gyres. These increase the poleward oceanic heat transport leading to Arctic sea-ice melting, Arctic warming amplification, and large-scale Atlantic warming, which in turn initiates wave-induced downward propagating negative NAM and stratospheric warming and therefore reverse the oscillation phase. This coupled variability improves the performance of statistical models, which project further weakening of North Atlantic Oscillation, North Atlantic cooling and hiatus in wintertime North Atlantic-Arctic sea-ice and global surface temperature just like the 1950s-1970s.

Lack of constraint on spatial and long-term temporal variability of the El Niño southern Oscillation (ENSO) and its sensitivity to external forcing limit our ability to evaluate climate models and ENSO future projections. Current knowledge of Holocene ENSO variability derived from paleoclimate reconstructions does not separate the role of insolation forcing from internal climate variability. Using an updated synthesis of coral and bivalve monthly resolved records, we build composite records of seasonality and interannual variability in four regions of the tropical Pacific: Eastern Pacific (EP), Central Pacific (CP), Western Pacific (WP) and South West Pacific (SWP). An analysis of the uncertainties due to the sampling of chaotic multidecadal to centennial variability by short records allows for an objective comparison with transient simulations (mid-Holocene to present) performed using four different Earth System models. Sea surface temperature and pseudo-δ18O are used in model-data comparisons to assess the potential influence of hydroclimate change on records. We confirm the significance of the Holocene ENSO minimum (HEM) 3-6ka compared to low frequency unforced modulation of ENSO, with a reduction of ENSO variance of ∼50 % in EP and ∼80 % in CP. The approach suggests that the increasing trend of ENSO since 6ka can be attributed to insolation, while models underestimate ENSO sensitivity to orbital forcing by a factor of 4.7 compared to data, even when accounting for the large multidecadal variability. Precession-induced change in seasonal temperature range is positively linked to ENSO variance in EP and to a lesser extent in other regions, in both models and observations. Our regional approach yields insights into the past spatial expression of ENSO across the tropical Pacific. In the SWP, today under the influence of the South Pacific Convergence Zone (SPCZ), interannual variability was increased by ∼200 % during the HEM, indicating that SPCZ variability is indep

The land is a pivotal component of the Earth and its climate system since many processes of natural variations in the climate system, which affect the environment and human society, are governed by the land surface. Hence, a good representation of the thermal and hydrological states of the land surface in climate models is important to have a realistic simulation of the coupling between the atmosphere and the lito-biosphere. An influencing factor for improving the realism of the ground energy and water balance in climate models is the depth of the land zero-flux Bottom Boundary Condition Placement (BBCP). Despite recent improvements in modeling land surface processes in climate models, only limited attention has been directed toward the effect of the BBCP in Land Surface Models (LSMs) and its impact on the representation of terrestrial thermodynamics. Previous analytical and modeling studies suggest that the simulation of subsurface thermodynamics in current-generation climate models is not accurate due to the zero-heat-flux BBCP being imposed too close to the surface. An insufficiently deep land component in current-generation climate models compromises the simulation of the terrestrial thermal state and can influence land-atmosphere interactions. Further improvements in LSMs relate to the representation and sensitivity of coupling processes between the ground thermodynamic and hydrological regimes. As moisture is one of the main drivers of near-surface climate interactions, the hydro-thermodynamic coupling is crucial for studying the impacts of perturbations caused by human activity. Under climate change conditions, some areas and ecosystems are more vulnerable to a rapidly warming world than others...

The transient climate response (TCR) is 20% higher in the Alfred Wegener Institute Climate Model (AWI-CM) compared to the Max Planck Institute Earth System Model (MPI-ESM) whereas the equilibrium climate sensitivity (ECS) is by up to 10% higher in AWI-CM. These results are largely independent of the two considered model resolutions for each model. The two coupled CMIP6 models share the same atmosphere-land component ECHAM6.3 developed at the Max Planck Institute for Meteorology (MPI-M). However, ECHAM6.3 is coupled to two different ocean models, namely the MPIOM sea ice-ocean model developed at MPI-M and the FESOM sea ice-ocean model developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). A reason for the different TCR is related to ocean heat uptake in response to greenhouse gas forcing. Specifically, AWI-CM simulations show stronger surface heating than MPI-ESM simulations while the latter accumulate more heat in the deeper ocean. The vertically integrated ocean heat content is increasing slower in AWI-CM model configurations compared to MPI-ESM model configurations in the high latitudes. Weaker vertical mixing in AWI-CM model configurations compared to MPI-ESM model configurations seems to be key for these differences. The strongest difference in vertical ocean mixing occurs inside the Weddell and Ross Gyres and the northern North Atlantic. Over the North Atlantic, these differences materialize in a lack of a warming hole in AWI-CM model configurations and the presence of a warming hole in MPI-ESM model configurations. All these differences occur largely independent of the considered model resolutions.

The land is a pivotal component of the Earth and its climate system since many processes of natural variations in the climate system, which affect the environment and human society, are governed by the land surface. Hence, a good representation of the thermal and hydrological states of the land surface in climate models is important to have a realistic simulation of the coupling between the atmosphere and the lito-biosphere. An influencing factor for improving the realism of the ground energy and water balance in climate models is the depth of the land zero-flux Bottom Boundary Condition Placement (BBCP). Despite recent improvements in modeling land surface processes in climate models, only limited attention has been directed toward the effect of the BBCP in Land Surface Models (LSMs) and its impact on the representation of terrestrial thermodynamics. Previous analytical and modeling studies suggest that the simulation of subsurface thermodynamics in current-generation climate models is not accurate due to the zero-heat-flux BBCP being imposed too close to the surface. An insufficiently deep land component in current-generation climate models compromises the simulation of the terrestrial thermal state and can influence land-atmosphere interactions. Further improvements in LSMs relate to the representation and sensitivity of coupling processes between the ground thermodynamic and hydrological regimes. As moisture is one of the main drivers of near-surface climate interactions, the hydro-thermodynamic coupling is crucial for studying the impacts of perturbations caused by human activity. Under climate change conditions, some areas and ecosystems are more vulnerable to a rapidly warming world than others...

The transient climate response (TCR) is 20% higher in the Alfred Wegener Institute Climate Model (AWI-CM) compared to the Max Planck Institute Earth System Model (MPI-ESM) whereas the equilibrium climate sensitivity (ECS) is by up to 10% higher in AWI-CM. These results are largely independent of the two considered model resolutions for each model. The two coupled CMIP6 models share the same atmosphere-land component ECHAM6.3 developed at the Max Planck Institute for Meteorology (MPI-M). However, ECHAM6.3 is coupled to two different ocean models, namely the MPIOM sea ice-ocean model developed at MPI-M and the FESOM sea ice-ocean model developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). A reason for the different TCR is related to ocean heat uptake in response to greenhouse gas forcing. Specifically, AWI-CM simulations show stronger surface heating than MPI-ESM simulations while the latter accumulate more heat in the deeper ocean. The vertically integrated ocean heat content is increasing slower in AWI-CM model configurations compared to MPI-ESM model configurations in the high latitudes. Weaker vertical mixing in AWI-CM model configurations compared to MPI-ESM model configurations seems to be key for these differences. The strongest difference in vertical ocean mixing occurs inside the Weddell and Ross Gyres and the northern North Atlantic. Over the North Atlantic, these differences materialize in a lack of a warming hole in AWI-CM model configurations and the presence of a warming hole in MPI-ESM model configurations. All these differences occur largely independent of the considered model resolutions.

Simulations of Common Era climate evolution coordinated by PMIP's "Past2K" working group together with multiproxy reconstructions from the PAGES 2k Network provide pivotal understanding for the evolution of the modern climate system and for expected changes in the near future.

© Author(s) 2020. We gratefully acknowledge the IlModels (CGL2014-59644-R) and GreatModelS (RTI2018-102305-B-C21) projects. We also thank the CESM1(CAM5) Last Millennium Ensemble Community Project and supercomputing resources provided by the NSF, CISL and Yellowstone. This research has been supported by the Spanish Ministry of Economy, Industry and Competitiveness (grant no. BES-2015-075019)., Borehole-based reconstruction is a wellestablished technique to recover information of the past climate variability based on two main hypotheses: (1) past ground surface temperature (GST) histories can be recovered from borehole temperature profiles (BTPs); (2) the past GST evolution is coupled to surface air temperature (SAT) changes, and thus, past SAT changes can be recovered from BTPs. Compared to some of the last millennium (LM) proxy-based reconstructions, previous studies based on the borehole technique indicate a larger temperature increase during the last few centuries. The nature of these differences has fostered the assessment of this reconstruction technique in search of potential causes of bias. Here, we expand previous works to explore potential methodological and physical biases using pseudo-proxy experiments with the Community Earth System Model Last Millennium Ensemble (CESM-LME). A heat-conduction forward model driven by simulated surface temperature is used to generate synthetic BTPs that are then inverted using singular value decomposition. This procedure is applied to the set of simulations that incorporates all of the LM external forcing factors as well as those that consider the concentration of the green house gases (GHGs) and the land use land cover (LULC) changes forcings separately. The results indicate that methodological issues may impact the representation of the simulated GST at different spatial scales, with the temporal logging of the BTPs as the main sampling issue that may lead to an underestimation of the simulated GST 20th-century trends. Our analysis also shows that in the surrogate reality of the CESM-LME the GST does not fully capture the SAT warming during the industrial period, and thus, there may be a further underestimation of the past SAT changes due to physical processes. Globally, this effect is mainly influenced by the GHG forcing, whereas regionally, LULC changes and other forcings factors also contribute. These findi, Ministerio de Economía y Competitividad (MINECO), Fac. de Ciencias Físicas, TRUE, pub

© Author(s) 2020. We gratefully acknowledge the IlModels (CGL2014-59644-R) and GreatModelS (RTI2018-102305-B-C21) projects. We also thank the CESM1(CAM5) Last Millennium Ensemble Community Project and supercomputing resources provided by the NSF, CISL and Yellowstone. This research has been supported by the Spanish Ministry of Economy, Industry and Competitiveness (grant no. BES-2015-075019)., Borehole-based reconstruction is a wellestablished technique to recover information of the past climate variability based on two main hypotheses: (1) past ground surface temperature (GST) histories can be recovered from borehole temperature profiles (BTPs); (2) the past GST evolution is coupled to surface air temperature (SAT) changes, and thus, past SAT changes can be recovered from BTPs. Compared to some of the last millennium (LM) proxy-based reconstructions, previous studies based on the borehole technique indicate a larger temperature increase during the last few centuries. The nature of these differences has fostered the assessment of this reconstruction technique in search of potential causes of bias. Here, we expand previous works to explore potential methodological and physical biases using pseudo-proxy experiments with the Community Earth System Model Last Millennium Ensemble (CESM-LME). A heat-conduction forward model driven by simulated surface temperature is used to generate synthetic BTPs that are then inverted using singular value decomposition. This procedure is applied to the set of simulations that incorporates all of the LM external forcing factors as well as those that consider the concentration of the green house gases (GHGs) and the land use land cover (LULC) changes forcings separately. The results indicate that methodological issues may impact the representation of the simulated GST at different spatial scales, with the temporal logging of the BTPs as the main sampling issue that may lead to an underestimation of the simulated GST 20th-century trends. Our analysis also shows that in the surrogate reality of the CESM-LME the GST does not fully capture the SAT warming during the industrial period, and thus, there may be a further underestimation of the past SAT changes due to physical processes. Globally, this effect is mainly influenced by the GHG forcing, whereas regionally, LULC changes and other forcings factors also contribute. These findi, Ministerio de Economía y Competitividad (MINECO), Fac. de Ciencias Físicas, TRUE, pub

The North Atlantic ocean warming hole has been linked to reduced tropical heat import. Model simulations show an anthropogenically forced increased heat export poleward from the region, by overturning and gyre circulation, and shortwave cloud feedback control the warming hole formation and growth. Despite global warming, a region in the North Atlantic ocean has been observed to cool, a phenomenon known as the warming hole. Its emergence has been linked to a slowdown of the Atlantic meridional overturning circulation, which leads to a reduced ocean heat transport into the warming hole region. Here we show that, in addition to the reduced low-latitude heat import, increased ocean heat transport out of the region into higher latitudes and a shortwave cloud feedback dominate the formation and temporal evolution of the warming hole under greenhouse gas forcing. In climate model simulations of the historical period, the low-latitude Atlantic meridional overturning circulation decline does not emerge from natural variability, whereas the accelerating heat transport to higher latitudes is clearly attributable to anthropogenic forcing. Both the overturning and the gyre circulation contribute to the increased high-latitude ocean heat transport, and therefore are critical to understand the past and future evolutions of the warming hole.

The North Atlantic ocean warming hole has been linked to reduced tropical heat import. Model simulations show an anthropogenically forced increased heat export poleward from the region, by overturning and gyre circulation, and shortwave cloud feedback control the warming hole formation and growth. Despite global warming, a region in the North Atlantic ocean has been observed to cool, a phenomenon known as the warming hole. Its emergence has been linked to a slowdown of the Atlantic meridional overturning circulation, which leads to a reduced ocean heat transport into the warming hole region. Here we show that, in addition to the reduced low-latitude heat import, increased ocean heat transport out of the region into higher latitudes and a shortwave cloud feedback dominate the formation and temporal evolution of the warming hole under greenhouse gas forcing. In climate model simulations of the historical period, the low-latitude Atlantic meridional overturning circulation decline does not emerge from natural variability, whereas the accelerating heat transport to higher latitudes is clearly attributable to anthropogenic forcing. Both the overturning and the gyre circulation contribute to the increased high-latitude ocean heat transport, and therefore are critical to understand the past and future evolutions of the warming hole.

The North Atlantic ocean warming hole has been linked to reduced tropical heat import. Model simulations show an anthropogenically forced increased heat export poleward from the region, by overturning and gyre circulation, and shortwave cloud feedback control the warming hole formation and growth. Despite global warming, a region in the North Atlantic ocean has been observed to cool, a phenomenon known as the warming hole. Its emergence has been linked to a slowdown of the Atlantic meridional overturning circulation, which leads to a reduced ocean heat transport into the warming hole region. Here we show that, in addition to the reduced low-latitude heat import, increased ocean heat transport out of the region into higher latitudes and a shortwave cloud feedback dominate the formation and temporal evolution of the warming hole under greenhouse gas forcing. In climate model simulations of the historical period, the low-latitude Atlantic meridional overturning circulation decline does not emerge from natural variability, whereas the accelerating heat transport to higher latitudes is clearly attributable to anthropogenic forcing. Both the overturning and the gyre circulation contribute to the increased high-latitude ocean heat transport, and therefore are critical to understand the past and future evolutions of the warming hole.

The North Atlantic ocean warming hole has been linked to reduced tropical heat import. Model simulations show an anthropogenically forced increased heat export poleward from the region, by overturning and gyre circulation, and shortwave cloud feedback control the warming hole formation and growth. Despite global warming, a region in the North Atlantic ocean has been observed to cool, a phenomenon known as the warming hole. Its emergence has been linked to a slowdown of the Atlantic meridional overturning circulation, which leads to a reduced ocean heat transport into the warming hole region. Here we show that, in addition to the reduced low-latitude heat import, increased ocean heat transport out of the region into higher latitudes and a shortwave cloud feedback dominate the formation and temporal evolution of the warming hole under greenhouse gas forcing. In climate model simulations of the historical period, the low-latitude Atlantic meridional overturning circulation decline does not emerge from natural variability, whereas the accelerating heat transport to higher latitudes is clearly attributable to anthropogenic forcing. Both the overturning and the gyre circulation contribute to the increased high-latitude ocean heat transport, and therefore are critical to understand the past and future evolutions of the warming hole.

Reconstructions of the global mean annual temperature evolution during the Holocene yield conflicting results. One temperature reconstruction shows global cooling during the late Holocene. The other reconstruction reveals global warming. Here we show that both a global warming mode and a cooling mode emerge when performing a spatio-temporal analysis of annual temperature variability during the Holocene using data from a transient climate model simulation. The warming mode is most pronounced in the tropics. The simulated cooling mode is determined by changes in the seasonal cycle of Arctic sea-ice that are forced by orbital variations and volcanic eruptions. The warming mode dominates in the mid-Holocene, whereas the cooling mode takes over in the late Holocene. The weighted sum of the two modes yields the simulated global temperature trend evolution. Our findings have strong implications for the interpretation of proxy data and the selection of proxy locations to compute global mean temperatures.

Reconstructions of the global mean annual temperature evolution during the Holocene yield conflicting results. One temperature reconstruction shows global cooling during the late Holocene. The other reconstruction reveals global warming. Here we show that both a global warming mode and a cooling mode emerge when performing a spatio-temporal analysis of annual temperature variability during the Holocene using data from a transient climate model simulation. The warming mode is most pronounced in the tropics. The simulated cooling mode is determined by changes in the seasonal cycle of Arctic sea-ice that are forced by orbital variations and volcanic eruptions. The warming mode dominates in the mid-Holocene, whereas the cooling mode takes over in the late Holocene. The weighted sum of the two modes yields the simulated global temperature trend evolution. Our findings have strong implications for the interpretation of proxy data and the selection of proxy locations to compute global mean temperatures.

Borehole-based reconstruction is a well-established technique to recover information of the past climate variability based on two main hypotheses: (1) past ground surface temperature (GST) histories can be recovered from borehole temperature profiles (BTPs); (2) the past GST evolution is coupled to surface air temperature (SAT) changes, and thus, past SAT changes can be recovered from BTPs. Compared to some of the last millennium (LM) proxy-based reconstructions, previous studies based on the borehole technique indicate a larger temperature increase during the last few centuries. The nature of these differences has fostered the assessment of this reconstruction technique in search of potential causes of bias. Here, we expand previous works to explore potential methodological and physical biases using pseudo-proxy experiments with the Community Earth System Model Last Millennium Ensemble (CESM-LME). A heat-conduction forward model driven by simulated surface temperature is used to generate synthetic BTPs that are then inverted using singular value decomposition. This procedure is applied to the set of simulations that incorporates all of the LM external forcing factors as well as those that consider the concentration of the green house gases (GHGs) and the land use land cover (LULC) changes forcings separately. The results indicate that methodological issues may impact the representation of the simulated GST at different spatial scales, with the temporal logging of the BTPs as the main sampling issue that may lead to an underestimation of the simulated GST 20th-century trends. Our analysis also shows that in the surrogate reality of the CESM-LME the GST does not fully capture the SAT warming during the industrial period, and thus, there may be a further underestimation of the past SAT changes due to physical processes. Globally, this effect is mainly influenced by the GHG forcing, whereas regionally, LULC changes and other forcings factors also contribute. These find

The North Atlantic ocean warming hole has been linked to reduced tropical heat import. Model simulations show an anthropogenically forced increased heat export poleward from the region, by overturning and gyre circulation, and shortwave cloud feedback control the warming hole formation and growth. Despite global warming, a region in the North Atlantic ocean has been observed to cool, a phenomenon known as the warming hole. Its emergence has been linked to a slowdown of the Atlantic meridional overturning circulation, which leads to a reduced ocean heat transport into the warming hole region. Here we show that, in addition to the reduced low-latitude heat import, increased ocean heat transport out of the region into higher latitudes and a shortwave cloud feedback dominate the formation and temporal evolution of the warming hole under greenhouse gas forcing. In climate model simulations of the historical period, the low-latitude Atlantic meridional overturning circulation decline does not emerge from natural variability, whereas the accelerating heat transport to higher latitudes is clearly attributable to anthropogenic forcing. Both the overturning and the gyre circulation contribute to the increased high-latitude ocean heat transport, and therefore are critical to understand the past and future evolutions of the warming hole.

The North Atlantic ocean warming hole has been linked to reduced tropical heat import. Model simulations show an anthropogenically forced increased heat export poleward from the region, by overturning and gyre circulation, and shortwave cloud feedback control the warming hole formation and growth. Despite global warming, a region in the North Atlantic ocean has been observed to cool, a phenomenon known as the warming hole. Its emergence has been linked to a slowdown of the Atlantic meridional overturning circulation, which leads to a reduced ocean heat transport into the warming hole region. Here we show that, in addition to the reduced low-latitude heat import, increased ocean heat transport out of the region into higher latitudes and a shortwave cloud feedback dominate the formation and temporal evolution of the warming hole under greenhouse gas forcing. In climate model simulations of the historical period, the low-latitude Atlantic meridional overturning circulation decline does not emerge from natural variability, whereas the accelerating heat transport to higher latitudes is clearly attributable to anthropogenic forcing. Both the overturning and the gyre circulation contribute to the increased high-latitude ocean heat transport, and therefore are critical to understand the past and future evolutions of the warming hole.

Reconstructions of the global mean annual temperature evolution during the Holocene yield conflicting results. One temperature reconstruction shows global cooling during the late Holocene. The other reconstruction reveals global warming. Here we show that both a global warming mode and a cooling mode emerge when performing a spatio-temporal analysis of annual temperature variability during the Holocene using data from a transient climate model simulation. The warming mode is most pronounced in the tropics. The simulated cooling mode is determined by changes in the seasonal cycle of Arctic sea-ice that are forced by orbital variations and volcanic eruptions. The warming mode dominates in the mid-Holocene, whereas the cooling mode takes over in the late Holocene. The weighted sum of the two modes yields the simulated global temperature trend evolution. Our findings have strong implications for the interpretation of proxy data and the selection of proxy locations to compute global mean temperatures.