Search

The existence and magnitude of the recently suggested global warming hiatus, or slowdown, have been strongly debated1–3. Although various physical processes4–8 have been examined to elucidate this phenomenon, the accuracy and completeness of observational data that comprise global average surface air temperature (SAT) datasets is a concern9,10. In particular, these datasets lack either complete geographic coverage or in situ observations over the Arctic, owing to the sparse observational network in this area 9 . As a consequence, the contribution of Arctic warming to global SAT changes may have been underestimated, leading to an uncertainty in the hiatus debate. Here, we constructed a new Arctic SAT dataset using the most recently updated global SATs 2 and a drifting buoys based Arctic SAT dataset 11 through employing the ‘data interpolating empirical orthogonal functions’ method 12 . Our estimate of global SAT rate of increase is around 0.112 °C per decade, instead of 0.05 °C per decade from IPCC AR5 1 , for 1998–2012. Analysis of this dataset shows that the amplified Arctic warming over the past decade has significantly contributed to a continual global warming trend, rather than a hiatus or slowdown. The Arctic is under-represented in surface temperature datasets and this could affect estimates of global warming. A new dataset with greater coverage of the Arctic shows a higher warming rate of 0.112 °C per decade compared to 0.005 °C from IPCC AR5.

Decadal prediction experiments of Beijing Climate Center climate system model version 1.1 (BCC-CSM1.1) participated in Coupled Model Intercomparison Project Phase 5 (CMIP5) had poor skill in extratropics of the North Atlantic, the initialization of which was done by relaxing modeled ocean temperature to the Simple Ocean Data Assimilation (SODA) reanalysis data. This study aims to improve the prediction skill of this model by using the assimilation technique in the initialization. New ocean data are firstly generated by assimilating the sea surface temperature (SST) of the Hadley Centre Sea Ice and Sea Surface Temperature (HadISST) dataset to the ocean model of BCC-CSM1.1 via Ensemble Optimum Interpolation (EnOI). Then a suite of decadal re-forecasts launched annually over the period 1961-2005 is carried out with simulated ocean temperature restored to the assimilated ocean data. Comparisons between the re-forecasts and previous CMIP5 forecasts show that the re-forecasts are more skillful in mid-to-high latitude SST of the North Atlantic. Improved prediction skill is also found for the Atlantic multi-decadal oscillation (AMO), which is consistent with the better skill of Atlantic meridional overturning circulation (AMOC) predicted by the re-forecasts. We conclude that the EnOI assimilation generates better ocean data than the SODA reanalysis for initializing decadal climate prediction of BCC-CSM1.1 model.

To highlight the compatibility of climate model simulation and proxy reconstruction at different timescales, a timescale separation merging method combining proxy records and climate model simulations is presented. Annual mean surface temperature anomalies for the last millennium (851–2005 AD) at various scales over the land of the Northern Hemisphere were reconstructed with 2° × 2° spatial resolution, using an optimal interpolation (OI) algorithm. All target series were decomposed using an ensemble empirical mode decomposition method followed by power spectral analysis. Four typical components were obtained at inter-annual, decadal, multidecadal, and centennial timescales. A total of 323 temperature-sensitive proxy chronologies were incorporated after screening for each component. By scaling the proxy components using variance matching and applying a localized OI algorithm to all four components point by point, we obtained merged surface temperatures. Independent validation indicates that the most significant improvement was for components at the inter-annual scale, but this became less evident with increasing timescales. In mid-latitude land areas, 10–30% of grids were significantly corrected at the inter-annual scale. By assimilating the proxy records, the merged results reduced the gap in response to volcanic forcing between a pure reconstruction and simulation. Difficulty remained in verifying the centennial information and quantifying corresponding uncertainties, so additional effort should be devoted to this aspect in future research.

Sea ice plays an important role in the air–ice–ocean interaction, but it is often represented simply in many regional atmospheric models. The Noah sea ice scheme, which is the only option in the current Weather Research and Forecasting (WRF) model (version 3.6.1), has a problem of energy imbalance due to its simplification in snow processes and lack of ablation and accretion processes in ice. Validated against the Surface Heat Budget of the Arctic Ocean (SHEBA) in situ observations, Noah underestimates the sea ice temperature which can reach −10 °C in winter. Sensitivity tests show that this bias is mainly attributed to the simulation within the ice when a time-dependent ice thickness is specified. Compared with the Noah sea ice model, the high-resolution thermodynamic snow and ice model (HIGHTSI) uses more realistic thermodynamics for snow and ice. Most importantly, HIGHTSI includes the ablation and accretion processes of sea ice and uses an interpolation method which can ensure the heat conservation during its integration. These allow the HIGHTSI to better resolve the energy balance in the sea ice, and the bias in sea ice temperature is reduced considerably. When HIGHTSI is coupled with the WRF model, the simulation of sea ice temperature by the original Polar WRF is greatly improved. Considering the bias with reference to SHEBA observations, WRF-HIGHTSI improves the simulation of surface temperature, 2 m air temperature and surface upward long-wave radiation flux in winter by 6, 5 °C and 20 W m−2, respectively. A discussion on the impact of specifying sea ice thickness in the WRF model is presented. Consistent with previous research, prescribing the sea ice thickness with observational information results in the best simulation among the available methods. If no observational information is available, we present a new method in which the sea ice thickness is initialized from empirical estimation and its further change is predicted by a complex thermodynamic sea ice model. The ice thickness simulated by this method depends much on the quality of the initial guess of the ice thickness and the role of the ice dynamic processes.

Regional climate models often lack detailed description of ice sheet surface and, as a result, are limited in their capability to provide useful information for Antarctic climate research and field campaigns. In this study, an upgraded scheme of surface physics for Antarctic ice sheet (IST) is developed to improve the surface temperature simulations in Antarctica. Through stand-alone simulations, IST shows advantages over the Noah glacial module, a commonly utilized scheme in the widely used Weather Research and Forecasting (WRF) model. These improvements are mainly attributed to the incorporation of detailed snow physics and optimized surface layer parameterization, which results in better simulations of both the surface albedo in summer and the turbulent sensible heat flux in winter. When coupled with IST instead of Noah, WRF models show improved simulation of surface temperatures throughout the year. The bias and root-mean-square-error of annual mean surface temperatures are reduced from 5.7 and 6.0 to 0.2 and 2.7 K.

Simulation and projection of the characteristics of heat waves over China were investigated using 12 CMIP5 global climate models and the CN05.1 observational gridded dataset. Four heat wave indices (heat wave frequency, longest heat wave duration, heat wave days, and high temperature days) were adopted in the analysis. Evaluations of the 12 CMIP5 models and their ensemble indicated that the multi-model ensemble could capture the spatiotemporal characteristics of heat wave variation over China. The inter-decadal variations of heat waves during 1961–2005 can be well simulated by multi-model ensemble. Based on model projections, the features of heat waves over China for eight different global warming targets (1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0 °C) were explored. The results showed that the frequency and intensity of heat waves would increase more dramatically as the global mean temperature rise attained higher warming targets. Under the RCP8.5 scenario, the four China-averaged heat wave indices would increase from about 1.0 times/year, 2.5, 5.4, and 13.8 days/year to about 3.2 times/year, 14.0, 32.0, and 31.9 days/year for 1.5 and 5.0 °C warming targets, respectively. Those regions that suffer severe heat waves in the base climate would experience the heat waves with greater frequency and severity following global temperature rise. It is also noteworthy that the areas in which a greater number of severe heat waves occur displayed considerable expansion. Moreover, the model uncertainties exhibit a gradual enhancement with projected time extending from 2006 to 2099.

A new dataset of surface temperature over North America has been constructed by merging climate model results and empirical tree-ring data through the application of an optimal interpolation algorithm. Errors of both the Community Climate System Model version 4 (CCSM4) simulation and the tree-ring reconstruction were considered to optimize the combination of the two elements. Variance matching was used to reconstruct the surface temperature series. The model simulation provided the background field, and the error covariance matrix was estimated statistically using samples from the simulation results with a running 31-year window for each grid. Thus, the merging process could continue with a time-varying gain matrix. This merging method (MM) was tested using two types of experiment, and the results indicated that the standard deviation of errors was about 0.4 °C lower than the tree-ring reconstructions and about 0.5 °C lower than the model simulation. Because of internal variabilities and uncertainties in the external forcing data, the simulated decadal warm–cool periods were readjusted by the MM such that the decadal variability was more reliable (e.g., the 1940–1960s cooling). During the two centuries (1601–1800 AD) of the preindustrial period, the MM results revealed a compromised spatial pattern of the linear trend of surface temperature, which is in accordance with the phase transition of the Pacific decadal oscillation and Atlantic multidecadal oscillation. Compared with pure CCSM4 simulations, it was demonstrated that the MM brought a significant improvement to the decadal variability of the gridded temperature via the merging of temperature-sensitive tree-ring records.

In this paper, we analyze the two dominant modes of winter surface air temperature (SAT) variations over China and their relationships with large-scale circulation anomalies. We then examine the fidelities of 20 individual models participating in the Coupled Model Inter-comparison Project Phase 5 in reproducing these two perspectives. Results showed that the winter SAT variations over China are dominated by two modes, a homogeneous warming pattern and a tripole pattern with warm departure in Northwest and Northeast China and cold departure in central and southern China. Consistent with the previous studies which documented the variations of the two modes are associated with the Siberian high and Arctic Oscillation (AO) anomalies, respectively, it is newly found that the variation of Empirical Orthogonal Function 2 (EOF2) mode is associated with the Northwest Pacific south–north dipole sea surface temperature anomaly in addition to the AO anomaly. Through comparisons with the observations, we identified that eight models outperform the others in simulating the two dominant modes and their relationships with large-scale circulation anomalies. These high-performing models were then selected to project future winter SAT changes over China under the Representative Concentration Pathway 4.5 (RCP4.5) scenario. Based on the multi-model ensemble mean, a nationwide warming was projected relative to the present climatology (1970–1999), with the largest increase in the Tibetan Plateau of 1.45 ± 0.62 °C by the period 2010–2039 and 2.87 ± 0.82 °C by the period 2050–2079; followed by Northeast China, Northwest China, North China, East China, Southwest China, and, finally, Southeast China.

Large spread appears in the projection of air-sea CO2 fluxes using the latest simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5). Here, two methods are applied to narrow this spread in 13 CMIP5 models. One method involves model selection based on the ability of models to reproduce the observed air-sea CO2 fluxes from 1980 to 2005. The other method involves constrained estimation based on the strong relationship between the historical and future air-sea CO2 fluxes. The estimated spread of the projected air-sea CO2 fluxes is effectively reduced by using these two approaches. These two approaches also show great agreement in the global ocean and three regional oceans of the equatorial Pacific Ocean, the North Atlantic Ocean and the Southern Ocean, including the average state and evolution characteristics. Based on the projections of the two approaches, the global ocean carbon uptake will increase in the first half of the 21st century then remain relatively stable and is projected to be 3.68–4.57 PgC/yr at the end of 21st century. The projections indicate that the increase in the CO2 uptake by the oceans will cease at the year of approximately 2070.

In the version of this Letter originally published, the increments on the y axis of Fig. 3 were incorrectly labelled as ‘0.0; 0.2; 0.2; 0.3’; they should have read ‘0.0; 0.1; 0.2; 0.3’. This has now been corrected in all versions of the Letter.

Large-scale climate history of the past millennium reconstructed solely from tree-ring data is prone to underestimate the amplitude of low-frequency variability. In this paper, we aimed at solving this problem by utilizing a novel method termed "MDVM", which was a combination of the ensemble empirical mode decomposition (EEMD) and variance matching techniques. We compiled a set of 211 tree-ring records from the extratropical Northern Hemisphere (30-90°N) in an effort to develop a new reconstruction of the annual mean temperature by the MDVM method. Among these dataset, a number of 126 records were screened out to reconstruct temperature variability longer than decadal scale for the period 850-2000 AD. The MDVM reconstruction depicted significant low-frequency variability in the past millennium with evident Medieval Warm Period (MWP) over the interval 950-1150 AD and pronounced Little Ice Age (LIA) cumulating in 1450-1850 AD. In the context of 1150-year reconstruction, the accelerating warming in 20th century was likely unprecedented, and the coldest decades appeared in the 1640s, 1600s and 1580s, whereas the warmest decades occurred in the 1990s, 1940s and 1930s. Additionally, the MDVM reconstruction covaried broadly with changes in natural radiative forcing, and especially showed distinct footprints of multiple volcanic eruptions in the last millennium. Comparisons of our results with previous reconstructions and model simulations showed the efficiency of the MDVM method on capturing low-frequency variability, particularly much colder signals of the LIA relative to the reference period. Our results demonstrated that the MDVM method has advantages in studying large-scale and low-frequency climate signals using pure tree-ring data.

Based on the historical and future outputs of 17 coupled model intercomparison project phase 5 (CMIP5) models, simulation of the precipitation extremes in China was evaluated under baseline climate condition compared to a gridded daily observation dataset CN05.1. The variations in precipitation extremes for eight global warming targets were also projected. The 17 individual models and the multi-model ensemble accurately reproduced the spatial distribution of precipitation extremes, although they were limited in their ability to capture the temporal characteristics. A notable dry bias existed in Southeast China, while a wet bias was present in North and Northwest China. The precipitation extremes in China were projected to be more frequent and more intense as global temperature rise reached the 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0 °C warming targets. The projected percentage changes in the annual number of days with precipitation >50 mm (R50) and total precipitation during days in which the daily precipitation exceeds the 99th percentile (R99p) are projected to increase by 25.81 and 69.14 % relative to the baseline climate for a 1.5 °C warming target, and by 95.52 and 162.00 % for a 4.0 °C warming target, respectively. As the global mean temperature rise increased from 1.5 to 5 °C, the subregions considerably affected by the East Asian summer monsoon (e.g., Southwest China, South China, and the Yangtze-Huai River Valley) were projected to experience a more dramatic increase in extreme precipitation events, in both number of days and intensity, while North and Northwest China were projected to suffer from relatively slight increases. The model uncertainties in the projected precipitation extremes in China by 17 CMIP5 models increase as global temperature rise increases.