A systematic approach to identify the sources of tropical SST errors in coupled models using the adjustment of initialised experiments

Understanding the sources of systematic errors in climate models is challenging because of coupled feedbacks and errors compensation. The developing seamless approach proposes that the identification and the correction of short term climate model errors have the potential to improve the modeled climate on longer time scales. In previous studies, initialised atmospheric simulations of a few days have been used to compare fast physics processes (convection, cloud processes) among models. The present study explores how initialised seasonal to decadal hindcasts (re-forecasts) relate transient week-to-month errors of the ocean and atmospheric components to the coupled model long-term pervasive SST errors. A protocol is designed to attribute the SST biases to the source processes. It includes five steps: (1) identify and describe biases in a coupled stabilized simulation, (2) determine the time scale of the advent of the bias and its propagation, (3) find the geographical origin of the bias, (4) evaluate the degree of coupling in the development of the bias, (5) find the field responsible for the bias. This strategy has been implemented with a set of experiments based on the initial adjustment of initialised simulations and exploring various degrees of coupling. In particular, hindcasts give the time scale of biases advent, regionally restored experiments show the geographical origin and ocean-only simulations isolate the field responsible for the bias and evaluate the degree of coupling in the bias development. This strategy is applied to four prominent SST biases of the IPSLCM5A-LR coupled model in the tropical Pacific, that are largely shared by other coupled models, including the Southeast Pacific warm bias and the equatorial cold tongue bias. Using the proposed protocol, we demonstrate that the East Pacific warm bias appears in a few months and is caused by a lack of upwelling due to too weak meridional coastal winds off Peru. The cold equatorial bias, which surprisingly takes 30 years to develop, is the result of an equatorward advection of midlatitude cold SST errors. Despite large development efforts, the current generation of coupled models shows only little improvement. The strategy proposed in this study is a further step to move from the current random ad hoc approach, to a bias-targeted, priority setting, systematic model development approach.

The role of dust in climate changes today, at the last glacial maximum and in the future.

Natural mineral aerosol (dust) is an active component of the climate system and plays multiple roles in mediating physical and biogeochemical exchanges between the atmosphere, land surface and ocean. Changes in the amount of dust in the atmosphere are caused both by changes in climate (precipitation, wind strength, regional moisture balance) and changes in the extent of dust sources caused by either anthropogenic or climatically induced changes in vegetation cover. Models of the global dust cycle take into account the physical controls on dust deflation from prescribed source areas (based largely on soil wetness and vegetation cover thresholds), dust transport within the atmospheric column, and dust deposition through sedimentation and scavenging by precipitation. These models successfully reproduce the first-order spatial and temporal patterns in atmospheric dust loading under modern conditions. Atmospheric dust loading was as much as an order-of-magnitude larger than today during the last glacial maximum (LGM). While the observed increase in emissions from northern Africa can be explained solely in terms of climate changes (colder, drier and windier glacial climates), increased emissions from other regions appear to have been largely a response to climatically induced changes in vegetation cover and hence in the extent of dust source areas. Model experiments suggest that the increased dust loading in tropical regions had an effect on radiative forcing comparable to that of low glacial CO2 levels. Changes in land-use are already increasing the dust loading of the atmosphere. However, simulations show that anthropogenically forced climate changes substantially reduce the extent and productivity of natural dust sources. Positive feedbacks initiated by a reduction of dust emissions from natural source areas on both radiative forcing and atmospheric CO2 could substantially mitigate the impacts of land-use changes, and need to be considered in climate change assessments.

How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets.

Global syntheses of palaeoenvironmental data are required to test climate models under conditions different from the present. Data sets for this purpose contain data from spatially extensive networks of sites. The data are either directly comparable to model output or readily interpretable in terms of modelled climate variables. Data sets must contain sufficient documentation to distinguish between raw (primary) and interpreted (secondary, tertiary) data, to evaluate the assumptions involved in interpretation of the data, to exercise quality control, and to select data appropriate for specific goals. Four data bases for the Late Quaternary, documenting changes in lake levels since 30 kyr BP (the Global Lake Status Data Base), vegetation distribution at 18 kyr and 6 kyr BP (BIOME 6000), aeolian accumulation rates during the last glacial-interglacial cycle (DIRTMAP), and tropical terrestrial climates at the Last Glacial Maximum (the LGM Tropical Terrestrial Data Synthesis) are summarised. Each has been used to evaluate simulations of Last Glacial Maximum (LGM: 21 calendar kyr BP) and/or mid-Holocene (6 cal. kyr BP) environments. Comparisons have demonstrated that changes in radiative forcing and orography due to orbital and ice-sheet variations explain the first-order, broad-scale (in space and time) features of global climate change since the LGM. However, atmospheric models forced by 6 cal. kyr BP orbital changes with unchanged surface conditions fail to capture quantitative aspects of the observed climate, including the greatly increased magnitude and northward shift of the African monsoon during the early to mid-Holocene. Similarly, comparisons with palaeoenvironmental datasets show that atmospheric models have underestimated the magnitude of cooling and drying of much of the land surface at the LGM. The inclusion of feedbacks due to changes in ocean- and land-surface conditions at both times, and atmospheric dust loading at the LGM, appears to be required in order to produce a better simulation of these past climates. The development of Earth system models incorporating the dynamic interactions among ocean, atmosphere, and vegetation is therefore mandated by Quaternary science results as well as climatological principles. For greatest scientific benefit, this development must be paralleled by continued advances in palaeodata analysis and synthesis, which in turn will help to define questions that call for new focused data collection efforts.

Dust sources and deposition during the last glacial maximum and current climate: a comparison of model results with paleodata from ice cores and marine sediments

Mineral dust aerosols in the atmosphere have the potential to affect the global climate by influencing the radiative balance of the atmosphere and the supply of micronutrients to the ocean. Ice and marine sediment cores indicate that dust deposition from the atmosphere was at some locations 2–20 times greater during glacial periods, raising the possibility that mineral aerosols might have contributed to climate change on glacial-interglacial time scales. To address this question, we have used linked terrestrial biosphere, dust source, and atmospheric transport models to simulate the dust cycle in the atmosphere for current and last glacial maximum (LGM) climates. We obtain a 2.5-fold higher dust loading in the entire atmosphere and a twenty-fold higher loading in high latitudes, in LGM relative to present. Comparisons to a compilation of atmospheric dust deposition flux estimates for LGM and present in marine sediment and ice cores show that the simulated flux ratios are broadly in agreement with observations; differences suggest where further improvements in the simple dust model could be made. The simulated increase in high-latitude dustiness depends on the expansion of unvegetated areas, especially in the high latitudes and in central Asia, caused by a combination of increased aridity and low atmospheric [CO2]. The existence of these dust source areas at the LGM is supported by pollen data and loess distribution in the northern continents. These results point to a role for vegetation feedbacks, including climate effects and physiological effects of low [CO2], in modulating the atmospheric distribution of dust.

Climate-driven spatial mismatches between British orchards and their pollinators: increased risks of pollination deficits

Understanding how climate change can affect crop-pollinator systems helps predict potential geographical mismatches between a crop and its pollinators, and therefore identify areas vulnerable to loss of pollination services. We examined the distribution of orchard species (apples, pears, plums and other top fruits) and their pollinators in Great Britain, for present and future climatic conditions projected for 2050 under the SRES A1B Emissions Scenario. We used a relative index of pollinator availability as a proxy for pollination service. At present there is a large spatial overlap between orchards and their pollinators, but predictions for 2050 revealed that the most suitable areas for orchards corresponded to low pollinator availability. However, we found that pollinator availability may persist in areas currently used for fruit production, but which are predicted to provide sub-optimal environmental suitability for orchard species in the future. Our results may be used to identify mitigation options to safeguard orchard production against the risk of pollination failure in Great Britain over the next 50 years; for instance choosing fruit tree varieties that are adapted to future climatic conditions, or boosting wild pollinators through improving landscape resources. Our approach can be readily applied to other regions and crop systems, and expanded to include different climatic scenarios.

Mid-Holocene land-surface conditions in northern Africa and the Arabian peninsula: a data set for the analysis of biogeophysical feedbacks in the climate system

Large changes in the extent of northern subtropical arid regions during the Holocene are attributed to orbitally forced variations in monsoon strength and have been implicated in the regulation of atmospheric trace gas concentrations on millenial timescales. Models that omit biogeophysical feedback, however, are unable to account for the full magnitude of African monsoon amplification and extension during the early to middle Holocene (˜9500–5000 years B.P.). A data set describing land-surface conditions 6000 years B.P. on a 1° × 1° grid across northern Africa and the Arabian Peninsula has been prepared from published maps and other sources of palaeoenvironmental data, with the primary aim of providing a realistic lower boundary condition for atmospheric general circulation model experiments similar to those performed in the Palaeoclimate Modelling Intercomparison Project. The data set includes information on the percentage of each grid cell occupied by specific vegetation types (steppe, savanna, xerophytic woods/scrub, tropical deciduous forest, and tropical montane evergreen forest), open water (lakes), and wetlands, plus information on the flow direction of major drainage channels for use in large-scale palaeohydrological modeling.

Arctic warming, atmospheric blocking and cold European winters in CMIP5 models

Amplified Arctic warming is expected to have a significant longterm influence on the midlatitude atmospheric circulation by the latter half of the 21st century. Potential influences of recent and near future Arctic changes on shorter timescales are much less clear, despite having received much recent attention in the literature. In this letter, climate models from the recent CMIP5 experiment are analysed for evidence of an influence of Arctic temperatures on midlatitude blocking and cold European winters in particular. The focus is on the variability of these features in detrended data and, in contrast to other studies, limited evidence of an influence is found. The occurrence of cold European winters is found to be largely independent of the temperature variability in the key Barents–Kara Sea region. Positive correlations of the Barents–Kara temperatures with Eurasian blocking are found in some models, but significant correlations are limited.

Predictability of the western North Pacific summer climate associated with different ENSO phases by ENSEMBLES multi-model seasonal forecasts

Predictability of the western North Pacific (WNP) summer climate associated with different El Niño–Southern Oscillation (ENSO) phases is investigated in this study based on the 1-month lead retrospective forecasts of five state-of-the-art coupled models from ENSEMBLES. During the period from 1960 to 2005, the models well capture the WNP summer climate anomalies during most of years in different ENSO phases except the La Niña decaying summers. In the El Niño developing, El Niño decaying and La Niña developing summers, the prediction skills are high for the WNP summer monsoon index (WNPMI), with the prediction correlation larger than 0.7. The high prediction skills of the lower-tropospheric circulation during these phases are found mainly over the tropical western Pacific Ocean, South China Sea and subtropical WNP. These good predictions correspond well to their close teleconnection with ENSO and the high prediction skills of tropical SSTs. By contrast, for the La Niña decaying summers, the prediction skills are considerably low with the prediction correlation for the WNPMI near to zero and low prediction skills around the Philippines and subtropical WNP. These poor predictions relate to the weak summer anomalies of the WNPMI during the La Niña decaying years and no significant connections between the WNP lower-tropospheric circulation anomalies and the SSTs over the tropical central and eastern Pacific Ocean in observations. However, the models tend to predict an apparent anomalous cyclone over the WNP during the La Niña decaying years, indicating a linearity of the circulation response over WNP in the models prediction in comparison with that during the El Niño decaying years which differs from observations. In addition, the models show considerable capability in describing the WNP summer anomalies during the ENSO neutral summers. These anomalies are related to the positive feedback between the WNP lower-tropospheric circulation and the local SSTs. The models can capture this positive feedback but with some uncertainties from different ensemble members during the ENSO neutral summers.

Numerical analyses of Runge–Kutta implicit–explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models

With the prospect of exascale computing, computational methods requiring only local data become especially attractive. Consequently, the typical domain decomposition of atmospheric models means horizontally-explicit vertically-implicit (HEVI) time-stepping schemes warrant further attention. In this analysis, Runge-Kutta implicit-explicit schemes from the literature are analysed for their stability and accuracy using a von Neumann stability analysis of two linear systems. Attention is paid to the numerical phase to indicate the behaviour of phase and group velocities. Where the analysis is tractable, analytically derived expressions are considered. For more complicated cases, amplification factors have been numerically generated and the associated amplitudes and phase diagnosed. Analysis of a system describing acoustic waves has necessitated attributing the three resultant eigenvalues to the three physical modes of the system. To do so, a series of algorithms has been devised to track the eigenvalues across the frequency space. The result enables analysis of whether the schemes exactly preserve the non-divergent mode; and whether there is evidence of spurious reversal in the direction of group velocities or asymmetry in the damping for the pair of acoustic modes. Frequency ranges that span next-generation high-resolution weather models to coarse-resolution climate models are considered; and a comparison is made of errors accumulated from multiple stability-constrained shorter time-steps from the HEVI scheme with a single integration from a fully implicit scheme over the same time interval. Two schemes, “Trap2(2,3,2)” and “UJ3(1,3,2)”, both already used in atmospheric models, are identified as offering consistently good stability and representation of phase across all the analyses. Furthermore, according to a simple measure of computational cost, “Trap2(2,3,2)” is the least expensive.

Species distribution models for crop pollination: a modelling framework applied to Great Britain

Insect pollination benefits over three quarters of the world's major crops. There is growing concern that observed declines in pollinators may impact on production and revenues from animal pollinated crops. Knowing the distribution of pollinators is therefore crucial for estimating their availability to pollinate crops; however, in general, we have an incomplete knowledge of where these pollinators occur. We propose a method to predict geographical patterns of pollination service to crops, novel in two elements: the use of pollinator records rather than expert knowledge to predict pollinator occurrence, and the inclusion of the managed pollinator supply. We integrated a maximum entropy species distribution model (SDM) with an existing pollination service model (PSM) to derive the availability of pollinators for crop pollination. We used nation-wide records of wild and managed pollinators (honey bees) as well as agricultural data from Great Britain. We first calibrated the SDM on a representative sample of bee and hoverfly crop pollinator species, evaluating the effects of different settings on model performance and on its capacity to identify the most important predictors. The importance of the different predictors was better resolved by SDM derived from simpler functions, with consistent results for bees and hoverflies. We then used the species distributions from the calibrated model to predict pollination service of wild and managed pollinators, using field beans as a test case. The PSM allowed us to spatially characterize the contribution of wild and managed pollinators and also identify areas potentially vulnerable to low pollination service provision, which can help direct local scale interventions. This approach can be extended to investigate geographical mismatches between crop pollination demand and the availability of pollinators, resulting from environmental change or policy scenarios.

Can a climate model reproduce extreme regional precipitation events over England and Wales?

The ability of the HiGEM climate model to represent high-impact, regional, precipitation events is investigated in two ways. The first focusses on a case study of extreme regional accumulation of precipitation during the passage of a summer extra-tropical cyclone across southern England on 20 July 2007 that resulted in a national flooding emergency. The climate model is compared with a global Numerical Weather Prediction (NWP) model and higher resolution, nested limited area models. While the climate model does not simulate the timing and location of the cyclone and associated precipitation as accurately as the NWP simulations, the total accumulated precipitation in all models is similar to the rain gauge estimate across England and Wales. The regional accumulation over the event is insensitive to horizontal resolution for grid spacings ranging from 90km to 4km. Secondly, the free-running climate model reproduces the statistical distribution of daily precipitation accumulations observed in the England-Wales precipitation record. The model distribution diverges increasingly from the record for longer accumulation periods with a consistent under-representation of more intense multi-day accumulations. This may indicate a lack of low-frequency variability associated with weather regime persistence. Despite this, the overall seasonal and annual precipitation totals from the model are still comparable to those from ERA-Interim.

An interdisciplinary modelling approach asessing the cost-effectiveness of agri-environmental measures on reducing nutrient concentration to WFD thresholds under climate change: the case of the Louros catchment

The application of the Water Framework Directive (WFD) in the European Union (EU) targets certain threshold levels for the concentration of various nutrients, nitrogen and phosphorous being the most important. In the EU, agri-environmental measures constitute a significant component of Pillar 2—Rural Development Policies in both financial and regulatory terms. Environmental measures also are linked to Pillar 1 payments through cross-compliance and the greening proposals. This paper drawing from work carried out in the REFRESH FP7 project aims to show how an INtegrated CAtchment model of plant/soil system dynamics and instream biogeochemical and hydrological dynamics can be used to assess the cost-effectiveness of agri-environmental measures in relation to nutrient concentration targets set by the WFD, especially in the presence of important habitats. We present the procedures (methodological steps, challenges and problems) for assessing the cost-effectiveness of agri-environmental measures at the baseline situation, and climate and land use change scenarios. Furthermore, we present results of an application of this methodology to the Louros watershed in Greece and discuss the likely uses and future extensions of the modelling approach. Finally, we attempt to reveal the importance of this methodology for designing and incorporating alternative environmental practices in Pillar 1 and 2 measures.

Global radiative and climate effect of the water vapour continuum at visible and near-infrared wavelengths

Recent laboratory measurements show that absorption by the water vapour continuum in near-infrared windows may be about an order of magnitude higher than assumed in many radiation codes. The radiative impact of the continuum at visible and near-infrared wavelengths is examined for the present day and for a possible future warmer climate (with a global-mean total column water increase of 33%). The calculations use a continuum model frequently used in climate models (‘CKD’) and a continuum model where absorption is enhanced at wavelengths greater than 1 µm based on recent measurements (‘CAVIAR’). The continuum predominantly changes the partitioning between solar radiation absorbed by the surface and the atmosphere; changes in top-of-atmosphere net irradiances are smaller. The global-mean clear-sky atmospheric absorption is enhanced by 1.5 W m−2 (about 2%) and 2.8 W m−2 (about 3.5%) for CKD and CAVIAR respectively, relative to a hypothetical no-continuum case, with all-sky enhancements about 80% of these values. The continuum is, in relative terms, more important for radiation budget changes between the present day and a possible future climate. Relative to the no-continuum case, the increase in global-mean clear-sky absorption is 8% higher using CKD and almost 20% higher using CAVIAR; all-sky enhancements are about half these values. The effect of the continuum is estimated for the solar component of the water vapour feedback, the reduction in downward surface irradiance and precipitation change in a warmer world. For CKD and CAVIAR respectively, and relative to the no-continuum case, the solar component of the water vapour feedback is enhanced by about 4 and 9%, the change in clear-sky downward surface irradiance is 7 and 18% more negative, and the global-mean precipitation response decreases by 1 and 4%. There is a continued need for improved continuum measurements, especially at atmospheric temperatures and at wavelengths below 2 µm.

Disparity-driven vs blur-driven models of accommodation and convergence in binocular vision and intermittent strabismus

Background. Current models of concomitant, intermittent strabismus, heterophoria, convergence and accommodation anomalies are either theoretically complex or incomplete. We propose an alternative and more practical way to conceptualize clinical patterns. Methods. In each of three hypothetical scenarios (normal; high AC/A and low CA/C ratios; low AC/A and high CA/C ratios) there can be a disparity-biased or blur-biased “style”, despite identical ratios. We calculated a disparity bias index (DBI) to reflect these biases. We suggest how clinical patterns fit these scenarios and provide early objective data from small illustrative clinical groups. Results. Normal adults and children showed disparity bias (adult DBI 0.43 (95%CI 0.50-0.36), child DBI 0.20 (95%CI 0.31-0.07) (p=0.001). Accommodative esotropes showed less disparity-bias (DBI 0.03). In the high AC/A and low CA/C scenario, early presbyopes had mean DBI of 0.17 (95%CI 0.28-0.06), compared to DBI of -0.31 in convergence excess esotropes. In the low AC/A and high CA/C scenario near exotropes had mean DBI of 0.27, while we predict that non-strabismic, non-amblyopic hyperopes with good vision without spectacles will show lower DBIs. Disparity bias ranged between 1.25 and -1.67. Conclusions. Establishing disparity or blur bias, together with knowing whether convergence to target demand exceeds accommodation or vice versa explains clinical patterns more effectively than AC/A and CA/C ratios alone. Excessive bias or inflexibility in near-cue use increases risk of clinical problems. We suggest clinicians look carefully at details of accommodation and convergence changes induced by lenses, dissociation and prisms and use these to plan treatment in relation to the model.

Land use and climate change impacts on soil organic carbon stocks in semi-arid Spain

Purpose The sensitivity of soil organic carbon to global change drivers, according to the depth profile, is receiving increasing attention because of its importance in the global carbon cycle and its potential feedback to climate change. A better knowledge of the vertical distribution of SOC and its controlling factors—the aim of this study—will help scientists predict the consequences of global change. Materials and methods The study area was the Murcia Province (S.E. Spain) under semiarid Mediterranean conditions. The database used consists of 312 soil profiles collected in a systematic grid, each 12 km2 covering a total area of 11,004 km2. Statistical analysis to study the relationships between SOC concentration and control factors in different soil use scenarios was conducted at fixed depths of 0–20, 20–40, 40–60, and 60–100 cm. Results and discussion SOC concentration in the top 40 cm ranged between 6.1 and 31.5 g kg−1, with significant differences according to land use, soil type and lithology, while below this depth, no differences were observed (SOC concentration 2.1–6.8 g kg−1). The ANOVA showed that land use was the most important factor controlling SOC concentration in the 0–40 cm depth. Significant differences were found in the relative importance of environmental and textural factors according to land use and soil depth. In forestland, mean annual precipitation and texture were the main predictors of SOC, while in cropland and shrubland, the main predictors were mean annual temperature and lithology. Total SOC stored in the top 1 m in the region was about 79 Tg with a low mean density of 7.18 kg Cm−3. The vertical distribution of SOC was shallower in forestland and deeper in cropland. A reduction in rainfall would lead to SOC decrease in forestland and shrubland, and an increase of mean annual temperature would adversely affect SOC in croplands and shrubland. With increasing depth, the relative importance of climatic factors decreases and texture becomes more important in controlling SOC in all land uses. Conclusions Due to climate change, impacts will be much greater in surface SOC, the strategies for C sequestration should be focused on subsoil sequestration, which was hindered in forestland due to bedrock limitations to soil depth. In these conditions, sequestration in cropland through appropriate management practices is recommended.