Search

The anthropogenically intensified greenhouse effect has caused a radiative imbalance at the top of the atmosphere during the industrial period. This, in turn, has led to an energy surplus in various components of the Earth system, with the ocean storing the largest part. The land contribution ranks second with the latest observational estimates based on borehole temperature profiles, which quantify the terrestrial energy surplus to be 6 % in the last 5 decades, whereas studies based on state-of-the-art climate models scale it down to 2 %. This underestimation stems from land surface models (LSMs) having a subsurface that is too shallow, which severely constrains the land heat uptake simulated by Earth system models (ESMs). A forced simulation of the last 2000 years with the Max Planck Institute ESM (MPI-ESM) using a deep LSM captures 4 times more heat than the standard shallow MPI-ESM simulations in the historical period, well above the estimates provided by other ESMs. However, deepening the LSM does not remarkably affect the simulated surface temperature. It is shown that the heat stored during the historical period by an ESM using a deep LSM component can be accurately estimated by considering the surface temperatures simulated by the ESM using a shallow LSM and propagating them with a standalone forward model. This result is used to derive estimates of land heat uptake using all available observational datasets, reanalysis products, and state-of-the-art ESM experiments. This approach yields values of 10.5–16.0 ZJ for 1971–2018, which are 12 %–42 % smaller than the latest borehole-based estimates (18.2 ZJ).

Understanding the differences between regional simulations of land–atmosphere interactions and near-surface conditions is crucial for a more reliable representation of past and future climate. Here, we explore the effect of changes in the model's horizontal resolution on the simulated energy balance at the surface and near-surface conditions using the Weather Research and Forecasting (WRF) model. To this aim, an ensemble of 12 simulations using three different horizontal resolutions (25, 50 and 100 km) and four different land surface model (LSM) configurations over North America from 1980 to 2013 is developed. Our results show that finer resolutions lead to higher surface net shortwave radiation and maximum temperatures at mid and high latitudes. At low latitudes over coastal areas, an increase in resolution leads to lower values of sensible heat flux and higher values of latent heat flux, as well as lower values of surface temperatures and higher values of precipitation, and soil moisture in summer. The use of finer resolutions leads then to an increase in summer values of latent heat flux and convective and non-convective precipitation and soil moisture at low latitudes. The effect of the LSM choice is larger than the effect of horizontal resolution on the near-surface temperature conditions. By contrast, the effect of the LSM choice on the simulation of precipitation is weaker than the effect of horizontal resolution, showing larger differences among LSM simulations in summer and over regions with high latent heat flux. Comparison between observations and the simulation of daily maximum and minimum temperatures and accumulated precipitation indicates that the CLM4 LSM yields the lowest biases in maximum and minimum mean temperatures but the highest biases in extreme temperatures. Increasing horizontal resolution leads to larger biases in accumulated precipitation over all regions particularly in summer. The reasons behind this are related to the partition between convective and non-convective precipitation, specially noticeable over western USA.

Energy exchanges among climate subsystems are of critical importance to determine the climate sensitivity of the Earth's system to greenhouse gases, to quantify the magnitude and evolution of the Earth's energy imbalance, and to project the evolution of future climate. Thus, ascertaining the magnitude of and change in the Earth's energy partition within climate subsystems has become urgent in recent years. Here, we provide new global estimates of changes in ground surface temperature, ground surface heat flux, and continental heat storage derived from geothermal data using an expanded database and new techniques. Results reveal markedly higher changes in ground heat flux and heat storage within the continental subsurface than previously reported, with land temperature changes of 1 K and continental heat gains of around 12 ZJ during the last part of the 20th century relative to preindustrial times. Half of the heat gain by the continental subsurface since 1960 has occurred in the last 20 years.

The representation and projection of extreme temperature and precipitation events in regional and global climate models are of major importance for the study of climate change impacts. However, state-of-the-art global and regional climate model simulations yield a broad inter-model range of intensity, duration and frequency of these extremes. Here, we present a modeling experiment using the Weather Research and Forecasting (WRF) model to determine the influence of the land surface model (LSM) component on uncertainties associated with extreme events. First, we analyze land–atmosphere interactions within four simulations performed by the WRF model from 1980 to 2012 over North America, using three different LSMs. Results show LSM-dependent differences at regional scales in the frequency of occurrence of events when surface conditions are altered by atmospheric forcing or land processes. The inter-model range of extreme statistics across the WRF simulations is large, particularly for indices related to the intensity and duration of temperature and precipitation extremes. Our results show that the WRF simulation of the climatology of heat extremes can be 5 ∘C warmer and 6 d longer depending on the employed LSM component, and similarly for cold extremes and heavy precipitation events. Areas showing large uncertainty in WRF-simulated extreme events are also identified in a model ensemble from three different regional climate model (RCM) simulations participating in the Coordinated Regional Climate Downscaling Experiment (CORDEX) project, revealing the implications of these results for other model ensembles. Thus, studies based on multi-model ensembles and reanalyses should include a variety of LSM configurations to account for the uncertainty arising from this model component or to test the performance of the selected LSM component before running the whole simulation. This study illustrates the importance of the LSM choice in climate simulations, supporting the development of new modeling studies using different LSM components to understand inter-model differences in simulating extreme temperature and precipitation events, which in turn will help to reduce uncertainties in climate model projections.

This is the second of two papers that document the creation of the New European Wind Atlas (NEWA). In Part 1, we described the sensitivity experiments and accompanying evaluation done to arrive at the final mesoscale model setup used to produce the mesoscale wind atlas. In this paper, Part 2, we document how we made the final wind atlas product, covering both the production of the mesoscale climatology generated with the Weather Research and Forecasting (WRF) model and the microscale climatology generated with the Wind Atlas Analysis and Applications Program (WAsP). The paper includes a detailed description of the technical and practical aspects that went into running the mesoscale simulations and the downscaling using WAsP. We show the main results from the final wind atlas and present a comprehensive evaluation of each component of the NEWA model chain using observations from a large set of tall masts located all over Europe. The added value of the WRF and WAsP downscaling of wind climatologies is evaluated relative to the performance of the driving ERA5 reanalysis and shows that the WRF downscaling reduces the mean wind speed bias and spread relative to that of ERA5 from -1.50±1.30 to 0.02±0.78 m s−1. The WAsP downscaling has an added positive impact relative to that of the WRF model in simple terrain. In complex terrain, where the assumptions of the linearized flow model break down, both the mean bias and spread in wind speed are worse than those from the raw mesoscale results.

This is the first of two papers that document the creation of the New European Wind Atlas (NEWA). It describes the sensitivity analysis and evaluation procedures that formed the basis for choosing the final setup of the mesoscale model simulations of the wind atlas. The suitable combination of model setup and parameterizations, bound by practical constraints, was found for simulating the climatology of the wind field at turbine-relevant heights with the Weather Research and Forecasting (WRF) model. Initial WRF model sensitivity experiments compared the wind climate generated by using two commonly used planetary boundary layer schemes and were carried out over several regions in Europe. They confirmed that the most significant differences in annual mean wind speed at 100 m a.g.l. (above ground level) mostly coincide with areas of high surface roughness length and not with the location of the domains or maximum wind speed. Then an ensemble of more than 50 simulations with different setups for a single year was carried out for one domain covering northern Europe for which tall mast observations were available. We varied many different parameters across the simulations, e.g. model version, forcing data, various physical parameterizations, and the size of the model domain. These simulations showed that although virtually every parameter change affects the results in some way, significant changes in the wind climate in the boundary layer are mostly due to using different physical parameterizations, especially the planetary boundary layer scheme, the representation of the land surface, and the prescribed surface roughness length. Also, the setup of the simulations, such as the integration length and the domain size, can considerably influence the results. We assessed the degree of similarity between winds simulated by the WRF ensemble members and the observations using a suite of metrics, including the Earth Mover's Distance (EMD), a statistic that measures the distance between two probability distributions. The EMD was used to diagnose the performance of each ensemble member using the full wind speed and direction distribution, which is essential for wind resource assessment. We identified the most realistic ensemble members to determine the most suitable configuration to be used in the final production run, which is fully described and evaluated in the second part of this study (Dörenkämper et al., 2020).

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 findings suggest that despite the larger temperature increase suggested by the borehole estimations during the last few centuries of the LM relative to some other proxy reconstructions, both the methodological and physical biases would result in a underestimation of the 20th-century warming.

Past climate variations may be uncovered via reconstruction methods that use proxy data as predictors. Among them, borehole reconstruction is a well-established technique to recover the long-term past surface air temperature (SAT) evolution. It is based on the assumption that SAT changes are strongly coupled to ground surface temperature (GST) changes and transferred to the subsurface by thermal conduction. We evaluate the SAT–GST coupling during the last millennium (LM) using simulations from the Community Earth System Model LM Ensemble (CESM-LME). The validity of such a premise is explored by analyzing the structure of the SAT–GST covariance during the LM and also by investigating the evolution of the long-term SAT–GST relationship. The multiple and single-forcing simulations in the CESM-LME are used to analyze the SAT–GST relationship within different regions and spatial scales and to derive the influence of the different forcing factors on producing feedback mechanisms that alter the energy balance at the surface. The results indicate that SAT–GST coupling is strong at global and above multi-decadal timescales in CESM-LME, although a relatively small variation in the long-term SAT–GST relationship is also represented. However, at a global scale such variation does not significantly impact the SAT–GST coupling, at local to regional scales this relationship experiences considerable long-term changes mostly after the end of the 19th century. Land use land cover changes are the main driver for locally and regionally decoupling SAT and GST, as they modify the land surface properties such as albedo, surface roughness and hydrology, which in turn modifies the energy fluxes at the surface. Snow cover feedbacks due to the influence of other external forcing are also important for corrupting the long-term SAT–GST coupling. Our findings suggest that such local and regional SAT–GST decoupling processes may represent a source of bias for SAT reconstructions from borehole measurement, since the thermal signature imprinted in the subsurface over the affected regions is not fully representative of the long-term SAT variations.

Understanding natural climate variability and its driving factors is crucial to assessing future climate change. Therefore, comparing proxy-based climate reconstructions with forcing factors as well as comparing these with paleoclimate model simulations is key to gaining insights into the relative roles of internal versus forced variability. A review of the state of modelling of the climate of the last millennium prior to the CMIP5–PMIP3 (Coupled Model Intercomparison Project Phase 5–Paleoclimate Modelling Intercomparison Project Phase 3) coordinated effort is presented and compared to the available temperature reconstructions. Simulations and reconstructions broadly agree on reproducing the major temperature changes and suggest an overall linear response to external forcing on multidecadal or longer timescales. Internal variability is found to have an important influence at hemispheric and global scales. The spatial distribution of simulated temperature changes during the transition from the Medieval Climate Anomaly to the Little Ice Age disagrees with that found in the reconstructions. Thus, either internal variability is a possible major player in shaping temperature changes through the millennium or the model simulations have problems realistically representing the response pattern to external forcing. A last millennium transient climate response (LMTCR) is defined to provide a quantitative framework for analysing the consistency between simulated and reconstructed climate. Beyond an overall agreement between simulated and reconstructed LMTCR ranges, this analysis is able to single out specific discrepancies between some reconstructions and the ensemble of simulations. The disagreement is found in the cases where the reconstructions show reduced covariability with external forcings or when they present high rates of temperature change.

The spatial context is critical when assessing present-day climate anomalies, attributing them to potential forcings and making statements regarding their frequency and severity in a long-term perspective. Recent international initiatives have expanded the number of high-quality proxy-records and developed new statistical reconstruction methods. These advances allow more rigorous regional past temperature reconstructions and, in turn, the possibility of evaluating climate models on policy-relevant, spatio-temporal scales. Here we provide a new proxy-based, annually-resolved, spatial reconstruction of the European summer (June–August) temperature fields back to 755 CE based on Bayesian hierarchical modelling (BHM), together with estimates of the European mean temperature variation since 138 BCE based on BHM and composite-plus-scaling (CPS). Our reconstructions compare well with independent instrumental and proxy-based temperature estimates, but suggest a larger amplitude in summer temperature variability than previously reported. Both CPS and BHM reconstructions indicate that the mean 20th century European summer temperature was not significantly different from some earlier centuries, including the 1st, 2nd, 8th and 10th centuries CE. The 1st century (in BHM also the 10th century) may even have been slightly warmer than the 20th century, but the difference is not statistically significant. Comparing each 50 yr period with the 1951–2000 period reveals a similar pattern. Recent summers, however, have been unusually warm in the context of the last two millennia and there are no 30 yr periods in either reconstruction that exceed the mean average European summer temperature of the last 3 decades (1986–2015 CE). A comparison with an ensemble of climate model simulations suggests that the reconstructed European summer temperature variability over the period 850–2000 CE reflects changes in both internal variability and external forcing on multi-decadal time-scales. For pan-European temperatures we find slightly better agreement between the reconstruction and the model simulations with high-end estimates for total solar irradiance. Temperature differences between the medieval period, the recent period and the Little Ice Age are larger in the reconstructions than the simulations. This may indicate inflated variability of the reconstructions, a lack of sensitivity and processes to changes in external forcing on the simulated European climate and/or an underestimation of internal variability on centennial and longer time scales.

East Asia has experienced strong warming since the 1960s accompanied by an increased frequency of heat waves and shrinking glaciers over the Tibetan Plateau and the Tien Shan. Here, we place the recent warmth in a long-term perspective by presenting a new spatially resolved warm-season (May-September) temperature reconstruction for the period 1-2000 CE using 59 multiproxy records from a wide range of East Asian regions. Our Bayesian Hierarchical Model (BHM) based reconstructions generally agree with earlier shorter regional temperature reconstructions but are more stable due to additional temperature sensitive proxies. We find a rather warm period during the first two centuries CE, followed by a multi-century long cooling period and again a warm interval covering the 900-1200 CE period (Medieval Climate Anomaly, MCA). The interval from 1450 to 1850 CE (Little Ice Age, LIA) was characterized by cooler conditions and the last 150 years are characterized by a continuous warming until recent times. Our results also suggest that the 1990s were likely the warmest decade in at least 1200 years. The comparison between an ensemble of climate model simulations and our summer reconstructions since 850 CE shows good agreement and an important role of internal variability and external forcing on multi-decadal time-scales.