Evaluation of biospheric components in Earth system models using modern and palaeo-observations: the state-of-the-art

Earth system models (ESMs) are increasing in complexity by incorporating more processes than their predecessors, making them potentially important tools for studying the evolution of climate and associated biogeochemical cycles. However, their coupled behaviour has only recently been examined in any detail, and has yielded a very wide range of outcomes. For example, coupled climate–carbon cycle models that represent land-use change simulate total land carbon stores at 2100 that vary by as much as 600 Pg C, given the same emissions scenario. This large uncertainty is associated with differences in how key processes are simulated in different models, and illustrates the necessity of determining which models are most realistic using rigorous methods of model evaluation. Here we assess the state-of-the-art in evaluation of ESMs, with a particular emphasis on the simulation of the carbon cycle and associated biospheric processes. We examine some of the new advances and remaining uncertainties relating to (i) modern and palaeodata and (ii) metrics for evaluation. We note that the practice of averaging results from many models is unreliable and no substitute for proper evaluation of individual models. We discuss a range of strategies, such as the inclusion of pre-calibration, combined process- and system-level evaluation, and the use of emergent constraints, that can contribute to the development of more robust evaluation schemes. An increasingly data-rich environment offers more opportunities for model evaluation, but also presents a challenge. Improved knowledge of data uncertainties is still necessary to move the field of ESM evaluation away from a "beauty contest" towards the development of useful constraints on model outcomes.

Global-scale climate impact functions: the relationship between climate forcing and impact

Although there is a strong policy interest in the impacts of climate change corresponding to different degrees of climate change, there is so far little consistent empirical evidence of the relationship between climate forcing and impact. This is because the vast majority of impact assessments use emissions-based scenarios with associated socio-economic assumptions, and it is not feasible to infer impacts at other temperature changes by interpolation. This paper presents an assessment of the global-scale impacts of climate change in 2050 corresponding to defined increases in global mean temperature, using spatially-explicit impacts models representing impacts in the water resources, river flooding, coastal, agriculture, ecosystem and built environment sectors. Pattern-scaling is used to construct climate scenarios associated with specific changes in global mean surface temperature, and a relationship between temperature and sea level used to construct sea level rise scenarios. Climate scenarios are constructed from 21 climate models to give an indication of the uncertainty between forcing and response. The analysis shows that there is considerable uncertainty in the impacts associated with a given increase in global mean temperature, due largely to uncertainty in the projected regional change in precipitation. This has important policy implications. There is evidence for some sectors of a non-linear relationship between global mean temperature change and impact, due to the changing relative importance of temperature and precipitation change. In the socio-economic sectors considered here, the relationships are reasonably consistent between socio-economic scenarios if impacts are expressed in proportional terms, but there can be large differences in absolute terms. There are a number of caveats with the approach, including the use of pattern-scaling to construct scenarios, the use of one impacts model per sector, and the sensitivity of the shape of the relationships between forcing and response to the definition of the impact indicator.

Costing conservation: an expert appraisal of the pollinator habitat benefits of England’s entry level stewardship

Pollination services provided by insects play a key role in English crop production and wider ecology. Despite growing evidence of the negative effects of habitat loss on pollinator populations, limited policy support is available to reverse this pressure. One measure that may provide beneficial habitat to pollinators is England’s entry level stewardship agri-environment scheme. This study uses a novel expert survey to develop weights for a range of models which adjust the balance of Entry Level Stewardship options within the current area of spending. The annual costs of establishing and maintaining these option compositions were estimated at £59.3–£12.4 M above current expenditure. Although this produced substantial reduction in private cost:benefit ratios, the benefits of the scheme to pollinator habitat rose by 7–140 %; significantly increasing the public cost:benefit ratio. This study demonstrates that the scheme has significant untapped potential to provide good quality habitat for pollinators across England, even within existing expenditure. The findings should open debate on the costs and benefits of specific entry level stewardship management options and how these can be enhanced to benefit both participants and biodiversity more equitably.

Humans use predictive kinematic models to calibrate visual cues to three-dimensional surface slant

When the sensory consequences of an action are systematically altered our brain can recalibrate the mappings between sensory cues and properties of our environment. This recalibration can be driven by both cue conflicts and altered sensory statistics, but neither mechanism offers a way for cues to be calibrated so they provide accurate information about the world, as sensory cues carry no information as to their own accuracy. Here, we explored whether sensory predictions based on internal physical models could be used to accurately calibrate visual cues to 3D surface slant. Human observers played a 3D kinematic game in which they adjusted the slant of a surface so that a moving ball would bounce off the surface and through a target hoop. In one group, the ball’s bounce was manipulated so that the surface behaved as if it had a different slant to that signaled by visual cues. With experience of this altered bounce, observers recalibrated their perception of slant so that it was more consistent with the assumed laws of kinematics and physical behavior of the surface. In another group, making the ball spin in a way that could physically explain its altered bounce eliminated this pattern of recalibration. Importantly, both groups adjusted their behavior in the kinematic game in the same way, experienced the same set of slants and were not presented with low-level cue conflicts that could drive the recalibration. We conclude that observers use predictive kinematic models to accurately calibrate visual cues to 3D properties of world.

The drivers of projected North Atlantic sea level change

Sea level change predicted by the CMIP5 atmosphere–ocean general circulation models (AOGCMs) is not spatially homogeneous. In particular, the sea level change in the North Atlantic is usually characterised by a meridional dipole pattern with higher sea level rise north of 40°N and lower to the south. The spread among models is also high in that region. Here we evaluate the role of surface buoyancy fluxes by carrying out simulations with the FAMOUS low-resolution AOGCM forced by surface freshwater and heat flux changes from CO2-forced climate change experiments with CMIP5 AOGCMs, and by a standard idealised surface freshwater flux applied in the North Atlantic. Both kinds of buoyancy flux change lead to the formation of the sea level dipole pattern, although the effect of the heat flux has a greater magnitude, and is the main cause of the spread of results among the CMIP5 models. By using passive tracers in FAMOUS to distinguish between additional and redistributed buoyancy, we show that the enhanced sea level rise north of 40°N is mainly due to the direct steric effect (the reduction of sea water density) caused by adding heat or freshwater locally. The surface buoyancy forcing also causes a weakening of the Atlantic meridional overturning circulation, and the consequent reduction of the northward ocean heat transport imposes a negative tendency on sea level rise, producing the reduced rise south of 40°N. However, unlike previous authors, we find that this indirect effect of buoyancy forcing is generally less important than the direct one, except in a narrow band along the east coast of the US, where it plays a major role and leads to sea level rise, as found by previous authors.

Attribution of the spatial pattern of CO2-forced sea level change to ocean surface flux changes

Climate models taking part in the coupled model intercomparison project phase 5 (CMIP5) all predict a global mean sea level rise for the 21st century. Yet the sea level change is not spatially uniform and differs among models. Here we evaluate the role of air–sea fluxes of heat, water and momentum (windstress) to find the spatial pattern associated to each of them as well as the spread they can account for. Using one AOGCM to which we apply the surface flux changes from other AOGCMs, we show that the heat flux and windstress changes dominate both the pattern and the spread, but taking the freshwater flux into account as well yields a sea level change pattern in better agreement with the CMIP5 ensemble mean. Differences among the CMIP5 control ocean temperature fields have a smaller impact on the sea level change pattern.

New models of North West European Holocene palaeogeography and inundation

This paper presents new 500 year interval palaeogeographic models for Britain, Ireland and the North West French coast from 11000 cal. BP to present. These models are used to calculate the varying rates of inundation for different geographical zones over the study period. This allows for consideration of the differential impact that Holocene sea-level rise had across space and time, and on past societies. In turn, consideration of the limitations of the models helps to foreground profitable areas for future research.

Effect of uncertainty in surface mass balance--elevation feedback on projections of the future sea level contribution of the Greenland ice sheet

We apply a new parameterisation of the Greenland ice sheet (GrIS) feedback between surface mass balance (SMB: the sum of surface accumulation and surface ablation) and surface elevation in the MAR regional climate model (Edwards et al., 2014) to projections of future climate change using five ice sheet models (ISMs). The MAR (Modèle Atmosphérique Régional: Fettweis, 2007) climate projections are for 2000–2199, forced by the ECHAM5 and HadCM3 global climate models (GCMs) under the SRES A1B emissions scenario. The additional sea level contribution due to the SMB– elevation feedback averaged over five ISM projections for ECHAM5 and three for HadCM3 is 4.3% (best estimate; 95% credibility interval 1.8–6.9 %) at 2100, and 9.6% (best estimate; 95% credibility interval 3.6–16.0 %) at 2200. In all results the elevation feedback is significantly positive, amplifying the GrIS sea level contribution relative to the MAR projections in which the ice sheet topography is fixed: the lower bounds of our 95% credibility intervals (CIs) for sea level contributions are larger than the “no feedback” case for all ISMs and GCMs. Our method is novel in sea level projections because we propagate three types of modelling uncertainty – GCM and ISM structural uncertainties, and elevation feedback parameterisation uncertainty – along the causal chain, from SRES scenario to sea level, within a coherent experimental design and statistical framework. The relative contributions to uncertainty depend on the timescale of interest. At 2100, the GCM uncertainty is largest, but by 2200 both the ISM and parameterisation uncertainties are larger. We also perform a perturbed parameter ensemble with one ISM to estimate the shape of the projected sea level probability distribution; our results indicate that the probability density is slightly skewed towards higher sea level contributions.

Climate models without pre-industrial volcanic forcing underestimate historical ocean thermal expansion

Episodic explosive volcanic eruptions are a natural part of the climate system but are often omitted from atmosphere-ocean general circulation model (AOGCM) preindustrial spin-up and control experiments. This omission imposes a negative bias on ocean heat uptake in simulations of the historical period. In models of a range of complexity, we find that global-mean sea level rise due to thermal expansion during the last ∼ 150 years is consequently underestimated by 5–30 mm, which is a substantial proportion of the model mean of 50 mm in Coupled Model Intercomparison Project Phase 3 AOGCMs with anthropogenic forcing only, and is therefore important in accounting for 20th century sea level rise. We test and recommend a procedure for removing the bias.

Twentieth-century global-mean sea-level rise: is the whole greater than the sum of the parts?

Confidence in projections of global-mean sea level rise (GMSLR) depends on an ability to account for GMSLR during the twentieth century. There are contributions from ocean thermal expansion, mass loss from glaciers and ice sheets, groundwater extraction, and reservoir impoundment. Progress has been made toward solving the “enigma” of twentieth-century GMSLR, which is that the observed GMSLR has previously been found to exceed the sum of estimated contributions, especially for the earlier decades. The authors propose the following: thermal expansion simulated by climate models may previously have been underestimated because of their not including volcanic forcing in their control state; the rate of glacier mass loss was larger than previously estimated and was not smaller in the first half than in the second half of the century; the Greenland ice sheet could have made a positive contribution throughout the century; and groundwater depletion and reservoir impoundment, which are of opposite sign, may have been approximately equal in magnitude. It is possible to reconstruct the time series of GMSLR from the quantified contributions, apart from a constant residual term, which is small enough to be explained as a long-term contribution from the Antarctic ice sheet. The reconstructions account for the observation that the rate of GMSLR was not much larger during the last 50 years than during the twentieth century as a whole, despite the increasing anthropogenic forcing. Semiempirical methods for projecting GMSLR depend on the existence of a relationship between global climate change and the rate of GMSLR, but the implication of the authors' closure of the budget is that such a relationship is weak or absent during the twentieth century.

Causes of plasma flow bursts and dayside auroral transients: An evaluation of two models Invoking reconnection pulses and changes in the Y component of the magnetosheath field

Longitudinal flow bursts observed by the European Incoherent Scatter (EISCAT) radar, in association with dayside auroral transients observed from Svalbard, have been interpreted as resulting from pulses of enhanced reconnection at the dayside magnetopause. However, an alternative model has recently been proposed for a steady rate of magnetopause reconnection, in which the bursts of longitudinal flow are due to increases in the field line curvature force, associated with the By component of the magnetosheath field. We here evaluate these two models, using observations on January 20, 1990, by EISCAT and a 630-nm all-sky camera at Ny Ålesund. For both models, we predict the behavior of both the dayside flows and the 630-nm emissions on newly opened field lines. It is shown that the signatures of steady reconnection and magnetosheath By changes could possibly resemble the observed 630-nm auroral events, but only for certain locations of the observing site, relative to the ionospheric projection of the reconnection X line: however, in such cases, the flow bursts would be seen between the 630-nm transients and not within them. On the other hand, the model of reconnection rate pulses predicts that the flows will be enhanced within each 630-nm transient auroral event. The observations on January 20, 1990, are shown to be consistent with the model of enhanced reconnection rate pulses over a background level and inconsistent with the effects of periodic enhancements of the magnitude of the magnetosheath By component. We estimate that the reconnection rate within the pulses would have to be at least an order of magnitude larger than the background level between the pulses.

Analysis of the regional pattern of sea level change due to ocean dynamics and density changes for 1993-2099 in observations and CMIP5 AOGCMs

Predictions of twenty-first century sea level change show strong regional variation. Regional sea level change observed by satellite altimetry since 1993 is also not spatially homogenous. By comparison with historical and pre-industrial control simulations using the atmosphere–ocean general circulation models (AOGCMs) of the CMIP5 project, we conclude that the observed pattern is generally dominated by unforced (internal generated) variability, although some regions, especially in the Southern Ocean, may already show an externally forced response. Simulated unforced variability cannot explain the observed trends in the tropical Pacific, but we suggest that this is due to inadequate simulation of variability by CMIP5 AOGCMs, rather than evidence of anthropogenic change. We apply the method of pattern scaling to projections of sea level change and show that it gives accurate estimates of future local sea level change in response to anthropogenic forcing as simulated by the AOGCMs under RCP scenarios, implying that the pattern will remain stable in future decades. We note, however, that use of a single integration to evaluate the performance of the pattern-scaling method tends to exaggerate its accuracy. We find that ocean volume mean temperature is generally a better predictor than global mean surface temperature of the magnitude of sea level change, and that the pattern is very similar under the different RCPs for a given model. We determine that the forced signal will be detectable above the noise of unforced internal variability within the next decade globally and may already be detectable in the tropical Atlantic.

An assessment of the impact of local processes on dust lifting in martian climate models

Simulation of the lifting of dust from the planetary surface is of substantially greater importance on Mars than on Earth, due to the fundamental role that atmospheric dust plays in the former’s climate, yet the dust emission parameterisations used to date in martian global climate models (MGCMs) lag, understandably, behind their terrestrial counterparts in terms of sophistication. Recent developments in estimating surface roughness length over all martian terrains and in modelling atmospheric circulations at regional to local scales (less than O(100 km)) presents an opportunity to formulate an improved wind stress lifting parameterisation. We have upgraded the conventional scheme by including the spatially varying roughness length in the lifting parameterisation in a fully consistent manner (thereby correcting a possible underestimation of the true threshold level for wind stress lifting), and used a modification to account for deviations from neutral stability in the surface layer. Following these improvements, it is found that wind speeds at typical MGCM resolution never reach the lifting threshold at most gridpoints: winds fall particularly short in the southern midlatitudes, where mean roughness is large. Sub-grid scale variability, manifested in both the near-surface wind field and the surface roughness, is then considered, and is found to be a crucial means of bridging the gap between model winds and thresholds. Both forms of small-scale variability contribute to the formation of dust emission ‘hotspots’: areas within the model gridbox with particularly favourable conditions for lifting, namely a smooth surface combined with strong near-surface gusts. Such small-scale emission could in fact be particularly influential on Mars, due both to the intense positive radiative feedbacks that can drive storm growth and a strong hysteresis effect on saltation. By modelling this variability, dust lifting is predicted at the locations at which dust storms are frequently observed, including the flushing storm sources of Chryse and Utopia, and southern midlatitude areas from which larger storms tend to initiate, such as Hellas and Solis Planum. The seasonal cycle of emission, which includes a double-peaked structure in northern autumn and winter, also appears realistic. Significant increases to lifting rates are produced for any sensible choices of parameters controlling the sub-grid distributions used, but results are sensitive to the smallest scale of variability considered, which high-resolution modelling suggests should be O(1 km) or less. Use of such models in future will permit the use of a diagnosed (rather than prescribed) variable gustiness intensity, which should further enhance dust lifting in the southern hemisphere in particular.

Metabolic theory predicts whole-ecosystem properties

Understanding the effects of individual organisms on material cycles and energy fluxes within ecosystems is central to predicting the impacts of human-caused changes on climate, land use, and biodiversity. Here we present a theory that integrates metabolic (organism-based bottom-up) and systems (ecosystem-based top-down) approaches to characterize how the metabolism of individuals affects the flows and stores of materials and energy in ecosystems. The theory predicts how the average residence time of carbon molecules, total system throughflow (TST), and amount of recycling vary with the body size and temperature of the organisms and with trophic organization. We evaluate the theory by comparing theoretical predictions with outputs of numerical models designed to simulate diverse ecosystem types and with empirical data for real ecosystems. Although residence times within different ecosystems vary by orders of magnitude—from weeks in warm pelagic oceans with minute phytoplankton producers to centuries in cold forests with large tree producers—as predicted, all ecosystems fall along a single line: residence time increases linearly with slope = 1.0 with the ratio of whole-ecosystem biomass to primary productivity (B/P). TST was affected predominantly by primary productivity and recycling by the transfer of energy from microbial decomposers to animal consumers. The theory provides a robust basis for estimating the flux and storage of energy, carbon, and other materials in terrestrial, marine, and freshwater ecosystems and for quantifying the roles of different kinds of organisms and environments at scales from local ecosystems to the biosphere.

Comparative analysis of intestinal tract models

The human gut is a complex ecosystem occupied by a diverse microbial community. Modulation of this microbiota impacts health and disease. The definitive way to investigate the impact of dietary intervention on the gut microbiota is a human trial. However, human trials are expensive and can be difficult to control; thus, initial screening is desirable. Utilization of a range of in vitro and in vivo models means that useful information can be gathered prior to the necessity for human intervention. This review discusses the benefits and limitations of these approaches.