Terrestrial biosphere changes over the last 120 kyr and their impact on ocean δ 13C

A new global synthesis and biomization of long (>40 kyr) pollen-data records is presented, and used with simulations from the HadCM3 and FAMOUS climate models to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Global modelled (BIOME4) biome distributions over time generally agree well with those inferred from pollen data. The two climate models show good agreement in global net primary productivity (NPP). NPP is strongly influenced by atmospheric carbon dioxide (CO2) concentrations through CO2 fertilization. The combined effects of modelled changes in vegetation and (via a simple model) soil carbon result in a global terrestrial carbon storage at the Last Glacial Maximum that is 210–470 Pg C less than in pre-industrial time. Without the contribution from exposed glacial continental shelves the reduction would be larger, 330–960 Pg C. Other intervals of low terrestrial carbon storage include stadial intervals at 108 and 85 ka BP, and between 60 and 65 ka BP during Marine Isotope Stage 4. Terrestrial carbon storage, determined by the balance of global NPP and decomposition, influences the stable carbon isotope composition (δ13C) of seawater because terrestrial organic carbon is depleted in 13C. Using a simple carbon-isotope mass balance equation we find agreement in trends between modelled ocean δ13C based on modelled land carbon storage, and palaeo-archives of ocean δ13C, confirming that terrestrial carbon storage variations may be important drivers of ocean δ13C changes.

Exploring the impact of CMIP5 model biases on the simulation of North Atlantic decadal variability

Instrumental observations, palaeo-proxies, and climate models suggest significant decadal variability within the North Atlantic subpolar gyre (NASPG). However, a poorly sampled observational record and a diversity of model behaviours mean that the precise nature and mechanisms of this variability are unclear. Here, we analyse an exceptionally large multi-model ensemble of 42 present-generation climate models to test whether NASPG mean state biases systematically affect the representation of decadal variability. Temperature and salinity biases in the Labrador Sea co-vary and influence whether density variability is controlled by temperature or salinity variations. Ocean horizontal resolution is a good predictor of the biases and the location of the dominant dynamical feedbacks within the NASPG. However, we find no link to the spectral characteristics of the variability. Our results suggest that the mean state and mechanisms of variability within the NASPG are not independent. This represents an important caveat for decadal predictions using anomaly-assimilation methods.

Serial clustering of extratropical cyclones in a multi-model ensemble of historical and future simulations

This study has investigated serial (temporal) clustering of extra-tropical cyclones simulated by 17 climate models that participated in CMIP5. Clustering was estimated by calculating the dispersion (ratio of variance to mean) of 30 December-February counts of Atlantic storm tracks passing nearby each grid point. Results from single historical simulations of 1975-2005 were compared to those from historical ERA40 reanalyses from 1958-2001 ERA40 and single future model projections of 2069-2099 under the RCP4.5 climate change scenario. Models were generally able to capture the broad features in reanalyses reported previously: underdispersion/regularity (i.e. variance less than mean) in the western core of the Atlantic storm track surrounded by overdispersion/clustering (i.e. variance greater than mean) to the north and south and over western Europe. Regression of counts onto North Atlantic Oscillation (NAO) indices revealed that much of the overdispersion in the historical reanalyses and model simulations can be accounted for by NAO variability. Future changes in dispersion were generally found to be small and not consistent across models. The overdispersion statistic, for any 30 year sample, is prone to large amounts of sampling uncertainty that obscures the climate change signal. For example, the projected increase in dispersion for storm counts near London in the CNRMCM5 model is 0.1 compared to a standard deviation of 0.25. Projected changes in the mean and variance of NAO are insufficient to create changes in overdispersion that are discernible above natural sampling variations.

Improved seasonal prediction of the hot summer of 2003 over Europe through better representation of uncertainty in the land surface

Methods to explicitly represent uncertainties in weather and climate models have been developed and refined over the past decade, and have reduced biases and improved forecast skill when implemented in the atmospheric component of models. These methods have not yet been applied to the land surface component of models. Since the land surface is strongly coupled to the atmospheric state at certain times and in certain places (such as the European summer of 2003), improvements in the representation of land surface uncertainty may potentially lead to improvements in atmospheric forecasts for such events. Here we analyse seasonal retrospective forecasts for 1981–2012 performed with the European Centre for Medium-Range Weather Forecasts’ (ECMWF) coupled ensemble forecast model. We consider two methods of incorporating uncertainty into the land surface model (H-TESSEL): stochastic perturbation of tendencies, and static perturbation of key soil parameters. We find that the perturbed parameter approach considerably improves the forecast of extreme air temperature for summer 2003, through better representation of negative soil moisture anomalies and upward sensible heat flux. Averaged across all the reforecasts the perturbed parameter experiment shows relatively little impact on the mean bias, suggesting perturbations of at least this magnitude can be applied to the land surface without any degradation of model climate. There is also little impact on skill averaged across all reforecasts and some evidence of overdispersion for soil moisture. The stochastic tendency experiments show a large overdispersion for the soil temperature fields, indicating that the perturbation here is too strong. There is also some indication that the forecast of the 2003 warm event is improved for the stochastic experiments, however the improvement is not as large as observed for the perturbed parameter experiment.

Direct and ozone-mediated forcing of the Southern Annular Mode by greenhouse gases

We assess the roles of long-lived greenhouse gases and ozone depletion in driving meridional surface pressure gradients in the southern extratropics; these gradients are a defining feature of the Southern Annular Mode. Stratospheric ozone depletion is thought to have caused a strengthening of this mode during summer, with increasing long-lived greenhouse gases playing a secondary role. Using a coupled atmosphere-ocean chemistry-climate model, we show that there is cancelation between the direct, radiative effect of increasing greenhouse gases by the also substantial indirect—chemical and dynamical—feedbacks that greenhouse gases have via their impact on ozone. This sensitivity of the mode to greenhouse gas-induced ozone changes suggests that a consistent implementation of ozone changes due to long-lived greenhouse gases in climate models benefits the simulation of this important aspect of Southern Hemisphere climate.

Improving a global model from the boundary layer: total turbulent energy and the neutral limit Prandtl number

Model intercomparisons have identified important deficits in the representation of the stable boundary layer by turbulence parametrizations used in current weather and climate models. However, detrimental impacts of more realistic schemes on the large-scale flow have hindered progress in this area. Here we implement a total turbulent energy scheme into the climate model ECHAM6. The total turbulent energy scheme considers the effects of Earth’s rotation and static stability on the turbulence length scale. In contrast to the previously used turbulence scheme, the TTE scheme also implicitly represents entrainment flux in a dry convective boundary layer. Reducing the previously exaggerated surface drag in stable boundary layers indeed causes an increase in southern hemispheric zonal winds and large-scale pressure gradients beyond observed values. These biases can be largely removed by increasing the parametrized orographic drag. Reducing the neutral limit turbulent Prandtl number warms and moistens low-latitude boundary layers and acts to reduce longstanding radiation biases in the stratocumulus regions, the Southern Ocean and the equatorial cold tongue that are common to many climate models.

Recent progress in understanding and projecting regional and global mean sea-level change

Considerable progress has been made in understanding the present and future regional and global sea level in the 2 years since the publication of the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change. Here, we evaluate how the new results affect the AR5’s assessment of (i) historical sea level rise, including attribution of that rise and implications for the sea level budget, (ii) projections of the components and of total global mean sea level (GMSL), and (iii) projections of regional variability and emergence of the anthropogenic signal. In each of these cases, new work largely provides additional evidence in support of the AR5 assessment, providing greater confidence in those findings. Recent analyses confirm the twentieth century sea level rise, with some analyses showing a slightly smaller rate before 1990 and some a slightly larger value than reported in the AR5. There is now more evidence of an acceleration in the rate of rise. Ongoing ocean heat uptake and associated thermal expansion have continued since 2000, and are consistent with ocean thermal expansion reported in the AR5. A significant amount of heat is being stored deeper in the water column, with a larger rate of heat uptake since 2000 compared to the previous decades and with the largest storage in the Southern Ocean. The first formal detection studies for ocean thermal expansion and glacier mass loss since the AR5 have confirmed the AR5 finding of a significant anthropogenic contribution to sea level rise over the last 50 years. New projections of glacier loss from two regions suggest smaller contributions to GMSL rise from these regions than in studies assessed by the AR5; additional regional studies are required to further assess whether there are broader implications of these results. Mass loss from the Greenland Ice Sheet, primarily as a result of increased surface melting, and from the Antarctic Ice Sheet, primarily as a result of increased ice discharge, has accelerated. The largest estimates of acceleration in mass loss from the two ice sheets for 2003–2013 equal or exceed the acceleration of GMSL rise calculated from the satellite altimeter sea level record over the longer period of 1993–2014. However, when increased mass gain in land water storage and parts of East Antarctica, and decreased mass loss from glaciers in Alaska and some other regions are taken into account, the net acceleration in the ocean mass gain is consistent with the satellite altimeter record. New studies suggest that a marine ice sheet instability (MISI) may have been initiated in parts of the West Antarctic Ice Sheet (WAIS), but that it will affect only a limited number of ice streams in the twenty-first century. New projections of mass loss from the Greenland and Antarctic Ice Sheets by 2100, including a contribution from parts of WAIS undergoing unstable retreat, suggest a contribution that falls largely within the likely range (i.e., two thirds probability) of the AR5. These new results increase confidence in the AR5 likely range, indicating that there is a greater probability that sea level rise by 2100 will lie in this range with a corresponding decrease in the likelihood of an additional contribution of several tens of centimeters above the likely range. In view of the comparatively limited state of knowledge and understanding of rapid ice sheet dynamics, we continue to think that it is not yet possible to make reliable quantitative estimates of future GMSL rise outside the likely range. Projections of twenty-first century GMSL rise published since the AR5 depend on results from expert elicitation, but we have low confidence in conclusions based on these approaches. New work on regional projections and emergence of the anthropogenic signal suggests that the two commonly predicted features of future regional sea level change (the increasing tilt across the Antarctic Circumpolar Current and the dipole in the North Atlantic) are related to regional changes in wind stress and surface heat flux. Moreover, it is expected that sea level change in response to anthropogenic forcing, particularly in regions of relatively low unforced variability such as the low-latitude Atlantic, will be detectable over most of the ocean by 2040. The east-west contrast of sea level trends in the Pacific observed since the early 1990s cannot be satisfactorily accounted for by climate models, nor yet definitively attributed either to unforced variability or forced climate change.

Mechanisms of the negative shortwave cloud feedback in mid to high latitudes

Increases in cloud optical depth and liquid water path (LWP) are robust features of global warming model simulations in high latitudes, yielding a negative shortwave cloud feedback, but the mechanisms are still uncertain. We assess the importance of microphysical processes for the negative optical depth feedback by perturbing temperature in the microphysics schemes of two aquaplanet models, both of which have separate prognostic equations for liquid water and ice. We find that most of the LWP increase with warming is caused by a suppression of ice microphysical processes in mixed-phase clouds, resulting in reduced conversion efficiencies of liquid water to ice and precipitation. Perturbing the temperature-dependent phase partitioning of convective condensate also yields a small LWP increase. Together, the perturbations in large-scale microphysics and convective condensate partitioning explain more than two-thirds of the LWP response relative to a reference case with increased SSTs, and capture all of the vertical structure of the liquid water response. In support of these findings, we show the existence of a very robust positive relationship between monthly-mean LWP and temperature in CMIP5 models and observations in mixed-phase cloud regions only. In models, the historical LWP sensitivity to temperature is a good predictor of the forced global warming response poleward of about 45°, although models appear to overestimate the LWP response to warming compared to observations. We conclude that in climate models, the suppression of ice-phase microphysical processes that deplete cloud liquid water is a key driver of the LWP increase with warming and of the associated negative shortwave cloud feedback.

Revisiting Swiss temperature trends 1959-2008

Temperature is a key variable for monitoring global climate change. Here we perform a trend analysis of Swiss temperatures from 1959 to 2008, using a new 2 × 2 km gridded data-set based on carefully homogenised ground observations from MeteoSwiss. The aim of this study is twofold: first, to discuss the spatial and altitudinal temperature trend characteristics in detail, and second, to quantify the contribution of changes in atmospheric circulation and local effects to these trends. The seasonal trends are all positive and mostly significant with an annual average warming rate of 0.35 °C/decade (∼1.6 times the northern hemispheric warming rate), ranging from 0.17 in autumn to 0.48 °C/decade in summer. Altitude-dependent trends are found in autumn and early winter where the trends are stronger at low altitudes (<800 m asl), and in spring where slightly stronger trends are found at altitudes close to the snow line. Part of the trends can be explained by changes in atmospheric circulation, but with substantial differences from season to season. In winter, circulation effects account for more than half the trends, while this contribution is much smaller in other seasons. After removing the effect of circulation, the trends still show seasonal variations with higher values in spring and summer. The circulation-corrected trends are closer to the values simulated by a set of ENSEMBLES regional climate models, with the models still tending towards a trend underestimation in spring and summer. Our results suggest that both circulation changes and more local effects are important to explain part of recent warming in spring, summer, and autumn. Snow-albedo feedback effects could be responsible for the stronger spring trends at altitudes close to the snow line, but the overall effect is small. In autumn, the observed decrease in fog frequency might be a key process in explaining the stronger temperature trends at low altitudes.

Future changes of wind energy potentials over Europe in a large CMIP5 multi-model ensemble

A statistical-dynamical downscaling method is used to estimate future changes of wind energy output (Eout) of a benchmark wind turbine across Europe at the regional scale. With this aim, 22 global climate models (GCMs) of the Coupled Model Intercomparison Project Phase 5 (CMIP5) ensemble are considered. The downscaling method uses circulation weather types and regional climate modelling with the COSMO-CLM model. Future projections are computed for two time periods (2021–2060 and 2061–2100) following two scenarios (RCP4.5 and RCP8.5). The CMIP5 ensemble mean response reveals a more likely than not increase of mean annual Eout over Northern and Central Europe and a likely decrease over Southern Europe. There is some uncertainty with respect to the magnitude and the sign of the changes. Higher robustness in future changes is observed for specific seasons. Except from the Mediterranean area, an ensemble mean increase of Eout is simulated for winter and a decreasing for the summer season, resulting in a strong increase of the intra-annual variability for most of Europe. The latter is, in particular, probable during the second half of the 21st century under the RCP8.5 scenario. In general, signals are stronger for 2061–2100 compared to 2021–2060 and for RCP8.5 compared to RCP4.5. Regarding changes of the inter-annual variability of Eout for Central Europe, the future projections strongly vary between individual models and also between future periods and scenarios within single models. This study showed for an ensemble of 22 CMIP5 models that changes in the wind energy potentials over Europe may take place in future decades. However, due to the uncertainties detected in this research, further investigations with multi-model ensembles are needed to provide a better quantification and understanding of the future changes.

Clouds and the atmospheric circulation response to warming

We study the effect of clouds on the atmospheric circulation response to CO2 quadrupling in an aquaplanet model with a slab-ocean lower boundary. The cloud effect is isolated by locking the clouds to either the control or 4xCO2 state in the shortwave (SW) or longwave (LW) radiation schemes. In our model, cloud-radiative changes explain more than half of the total poleward expansion of the Hadley cells, midlatitude jets, and storm tracks under CO2 quadrupling, even though they cause only one-fourth of the total global-mean surface warming. The effect of clouds on circulation results mainly from the SW cloud-radiative changes, which strongly enhance the Equator-to-pole temperature gradient at all levels in the troposphere, favoring stronger and poleward-shifted midlatitude eddies. By contrast, quadrupling CO2 while holding the clouds fixed causes strong polar amplification and weakened midlatitude baroclinicity at lower levels, yielding only a small poleward expansion of the circulation. Our results show that (a) the atmospheric circulation responds sensitively to cloud-driven changes in meridional and vertical temperature distribution, and (b) the spatial structure of cloud feedbacks likely plays a dominant role in the circulation response to greenhouse gas forcing. While the magnitude and spatial structure of the cloud feedback are expected to be highly model-dependent, an analysis of 4xCO2 simulations of CMIP5 models shows that the SW cloud feedback likely forces a poleward expansion of the tropospheric circulation in most climate models.

Improved Arctic sea-ice thickness projections using bias corrected CMIP5 simulations

Projections of Arctic sea ice thickness (SIT) have the potential to inform stakeholders about accessibility to the region, but are currently rather uncertain. The latest suite of CMIP5 Global Climate Models (GCMs) produce a wide range of simulated SIT in the historical period (1979–2014) and exhibit various biases when compared with the Pan-Arctic Ice Ocean Modelling and Assimilation System (PIOMAS) sea ice reanalysis. We present a new method to constrain such GCM simulations of SIT via a statistical bias correction technique. The bias correction successfully constrains the spatial SIT distribution and temporal variability in the CMIP5 projections whilst retaining the climatic fluctuations from individual ensemble members. The bias correction acts to reduce the spread in projections of SIT and reveals the significant contributions of climate internal variability in the first half of the century and of scenario uncertainty from mid-century onwards. The projected date of ice-free conditions in the Arctic under the RCP8.5 high emission scenario occurs in the 2050s, which is a decade earlier than without the bias correction, with potentially significant implications for stakeholders in the Arctic such as the shipping industry. The bias correction methodology developed could be similarly applied to other variables to reduce spread in climate projections more generally.

Current model capabilities for simulating black carbon and sulfate concentrations in the Arctic atmosphere: a multi-model evaluation using a comprehensive measurement data set

The concentrations of sulfate, black carbon (BC) and other aerosols in the Arctic are characterized by high values in late winter and spring (so-called Arctic Haze) and low values in summer. Models have long been struggling to capture this seasonality and especially the high concentrations associated with Arctic Haze. In this study, we evaluate sulfate and BC concentrations from eleven different models driven with the same emission inventory against a comprehensive pan-Arctic measurement data set over a time period of 2 years (2008–2009). The set of models consisted of one Lagrangian particle dispersion model, four chemistry transport models (CTMs), one atmospheric chemistry-weather forecast model and five chemistry climate models (CCMs), of which two were nudged to meteorological analyses and three were running freely. The measurement data set consisted of surface measurements of equivalent BC (eBC) from five stations (Alert, Barrow, Pallas, Tiksi and Zeppelin), elemental carbon (EC) from Station Nord and Alert and aircraft measurements of refractory BC (rBC) from six different campaigns. We find that the models generally captured the measured eBC or rBC and sulfate concentrations quite well, compared to previous comparisons. However, the aerosol seasonality at the surface is still too weak in most models. Concentrations of eBC and sulfate averaged over three surface sites are underestimated in winter/spring in all but one model (model means for January–March underestimated by 59 and 37 % for BC and sulfate, respectively), whereas concentrations in summer are overestimated in the model mean (by 88 and 44 % for July–September), but with overestimates as well as underestimates present in individual models. The most pronounced eBC underestimates, not included in the above multi-site average, are found for the station Tiksi in Siberia where the measured annual mean eBC concentration is 3 times higher than the average annual mean for all other stations. This suggests an underestimate of BC sources in Russia in the emission inventory used. Based on the campaign data, biomass burning was identified as another cause of the modeling problems. For sulfate, very large differences were found in the model ensemble, with an apparent anti-correlation between modeled surface concentrations and total atmospheric columns. There is a strong correlation between observed sulfate and eBC concentrations with consistent sulfate/eBC slopes found for all Arctic stations, indicating that the sources contributing to sulfate and BC are similar throughout the Arctic and that the aerosols are internally mixed and undergo similar removal. However, only three models reproduced this finding, whereas sulfate and BC are weakly correlated in the other models. Overall, no class of models (e.g., CTMs, CCMs) performed better than the others and differences are independent of model resolution.

Local adaptation in migrated interior Douglas-fir seedlings is mediated by ectomycorrhizas and other soil factors

Separating edaphic impacts on tree distributions from those of climate and geography is notoriously difficult. Aboveground and belowground factors play important roles, and determining their relative contribution to tree success will greatly assist in refining predictive models and forestry strategies in a changing climate. In a common glasshouse, seedlings of interior Douglas-fir (Pseudotsuga menziesii var. glauca) from multiple populations were grown in multiple forest soils. Fungicide was applied to half of the seedlings to separate soil fungal and nonfungal impacts on seedling performance. Soils of varying geographic and climatic distance from seed origin were compared, using a transfer function approach. Seedling height and biomass were optimized following seed transfer into drier soils, whereas survival was optimized when elevation transfer was minimised. Fungicide application reduced ectomycorrhizal root colonization by c. 50%, with treated seedlings exhibiting greater survival but reduced biomass. Local adaptation of Douglas-fir populations to soils was mediated by soil fungi to some extent in 56% of soil origin by response variable combinations. Mediation by edaphic factors in general occurred in 81% of combinations. Soil biota, hitherto unaccounted for in climate models, interacts with biogeography to influence plant ranges in a changing climate.

Links between tropical Pacific seasonal, interannual and orbital variability during the Holocene

The El Niño/Southern Oscillation (ENSO) is the leading mode of interannual climate variability. However, it is unclear how ENSO has responded to external forcing, particularly orbitally induced changes in the amplitude of the seasonal cycle during the Holocene. Here we present a reconstruction of seasonal and interannual surface conditions in the tropical Pacific Ocean from a network of high-resolution coral and mollusc records that span discrete intervals of the Holocene. We identify several intervals of reduced variance in the 2 to 7 yr ENSO band that are not in phase with orbital changes in equatorial insolation, with a notable 64% reduction between 5,000 and 3,000 years ago. We compare the reconstructed ENSO variance and seasonal cycle with that simulated by nine climate models that include orbital forcing, and find that the models do not capture the timing or amplitude of ENSO variability, nor the mid-Holocene increase in seasonality seen in the observations; moreover, a simulated inverse relationship between the amplitude of the seasonal cycle and ENSO-related variance in sea surface temperatures is not found in our reconstructions. We conclude that the tropical Pacific climate is highly variable and subject to millennial scale quiescent periods. These periods harbour no simple link to orbital forcing, and are not adequately simulated by the current generation of models.