Statistical decadal predictions for sea surface temperatures: a benchmark for dynamical GCM predictions

Accurate decadal climate predictions could be used to inform adaptation actions to a changing climate. The skill of such predictions from initialised dynamical global climate models (GCMs) may be assessed by comparing with predictions from statistical models which are based solely on historical observations. This paper presents two benchmark statistical models for predicting both the radiatively forced trend and internal variability of annual mean sea surface temperatures (SSTs) on a decadal timescale based on the gridded observation data set HadISST. For both statistical models, the trend related to radiative forcing is modelled using a linear regression of SST time series at each grid box on the time series of equivalent global mean atmospheric CO2 concentration. The residual internal variability is then modelled by (1) a first-order autoregressive model (AR1) and (2) a constructed analogue model (CA). From the verification of 46 retrospective forecasts with start years from 1960 to 2005, the correlation coefficient for anomaly forecasts using trend with AR1 is greater than 0.7 over parts of extra-tropical North Atlantic, the Indian Ocean and western Pacific. This is primarily related to the prediction of the forced trend. More importantly, both CA and AR1 give skillful predictions of the internal variability of SSTs in the subpolar gyre region over the far North Atlantic for lead time of 2 to 5 years, with correlation coefficients greater than 0.5. For the subpolar gyre and parts of the South Atlantic, CA is superior to AR1 for lead time of 6 to 9 years. These statistical forecasts are also compared with ensemble mean retrospective forecasts by DePreSys, an initialised GCM. DePreSys is found to outperform the statistical models over large parts of North Atlantic for lead times of 2 to 5 years and 6 to 9 years, however trend with AR1 is generally superior to DePreSys in the North Atlantic Current region, while trend with CA is superior to DePreSys in parts of South Atlantic for lead time of 6 to 9 years. These findings encourage further development of benchmark statistical decadal prediction models, and methods to combine different predictions.

Drag produced by trapped lee waves and propagating mountain waves in a two-layer atmosphere

The surface drag force produced by trapped lee waves and upward propagating waves in non-hydrostatic stratified flow over a mountain ridge is explicitly calculated using linear theory for a two-layer atmosphere with piecewise-constant static stability and wind speed profiles. The behaviour of the drag normalized by its hydrostatic single-layer reference value is investigated as a function of the ratio of the Scorer parameters in the two layers l_2/l_1 and of the corresponding dimensionless interface height l_1 H, for selected values of the dimensionless ridge width l_1 a and ratio of wind speeds in the two layers. When l_2/l_1 → 1, the propagating wave drag approaches 1 in approximately hydrostatic conditions, and the trapped lee wave drag vanishes. As l_2/l_1 decreases, the propagating wave drag progressively displays an oscillatory behaviour with l_1 H, with maxima of increasing magnitude due to constructive interference of reflected waves in the lower layer. The trapped lee wave drag shows localized maxima associated with each resonant trapped lee wave mode, occurring for small l_2/l_1 and slightly higher values of l_1 H than the propagating wave drag maxima. As l1a decreases, i.e. the flow becomes more non-hydrostatic, the propagating wave drag decreases and the regions of non-zero trapped lee wave drag extend to higher l_2/l_1. These results are confirmed by numerical simulations for l_2/l_1 = 0.2. In parameter ranges of meteorological relevance, the trapped lee wave drag may have a magnitude comparable to that of propagating wave drag, and be larger than the reference single-layer drag. This may have implications for drag parametrization in global climate and weather-prediction models.

The coupled model predictability of the Western North Pacific summer monsoon with different leading times

Leading time length is an important issue for modeling seasonal forecasts. In this study, a comparison of the interannual predictability of the Western North Pacific (WNP) summer monsoon between different leading months was performed by using one-, four-, and seven-month lead retrospective forecasts (hindcasts) of four coupled models from Ensembles-Based Predictions of Climate Changes and Their Impacts (ENSEMBLES) for the period of 1960-2005. It is found that the WNP summer anomalies, including lower-tropospheric circulation and precipitation anomalies, can be well predicted for all these leading months. The accuracy of the four-month lead prediction is only slightly weaker than that of the one-month lead prediction, although the skill decreases with the increase of leading months.

Precipitation distributions for explicit versus parameterized convection in a large-domain high-resolution tropical case study

Global climate and weather models tend to produce rainfall that is too light and too regular over the tropical ocean. This is likely because of convective parametrizations, but the problem is not well understood. Here, distributions of precipitation rates are analyzed for high-resolution UK Met Office Unified Model simulations of a 10 day case study over a large tropical domain (∼20°S–20°N and 42°E–180°E). Simulations with 12 km grid length and parametrized convection have too many occurrences of light rain and too few of heavier rain when interpolated onto a 1° grid and compared with Tropical Rainfall Measuring Mission (TRMM) data. In fact, this version of the model appears to have a preferred scale of rainfall around 0.4 mm h−1 (10 mm day−1), unlike observations of tropical rainfall. On the other hand, 4 km grid length simulations with explicit convection produce distributions much more similar to TRMM observations. The apparent preferred scale at lighter rain rates seems to be a feature of the convective parametrization rather than the coarse resolution, as demonstrated by results from 12 km simulations with explicit convection and 40 km simulations with parametrized convection. In fact, coarser resolution models with explicit convection tend to have even more heavy rain than observed. Implications for models using convective parametrizations, including interactions of heating and moistening profiles with larger scales, are discussed. One important implication is that the explicit convection 4 km model has temperature and moisture tendencies that favour transitions in the convective regime. Also, the 12 km parametrized convection model produces a more stable temperature profile at its extreme high-precipitation range, which may reduce the chance of very heavy rainfall. Further study is needed to determine whether unrealistic precipitation distributions are due to some fundamental limitation of convective parametrizations or whether parametrizations can be improved, in order to better simulate these distributions.

An earth-system prediction initiative for the twenty-first century

The necessity and benefits for establishing the international Earth-system Prediction Initiative (EPI) are discussed by scientists associated with the World Meteorological Organization (WMO) World Weather Research Programme (WWRP), World Climate Research Programme (WCRP), International Geosphere–Biosphere Programme (IGBP), Global Climate Observing System (GCOS), and natural-hazards and socioeconomic communities. The proposed initiative will provide research and services to accelerate advances in weather, climate, and Earth system prediction and the use of this information by global societies. It will build upon the WMO, the Group on Earth Observations (GEO), the Global Earth Observation System of Systems (GEOSS) and the International Council for Science (ICSU) to coordinate the effort across the weather, climate, Earth system, natural-hazards, and socioeconomic disciplines. It will require (i) advanced high-performance computing facilities, supporting a worldwide network of research and operational modeling centers, and early warning systems; (ii) science, technology, and education projects to enhance knowledge, awareness, and utilization of weather, climate, environmental, and socioeconomic information; (iii) investments in maintaining existing and developing new observational capabilities; and (iv) infrastructure to transition achievements into operational products and services.

Sea ice rheology

The polar oceans of Earth are covered by sea ice. On timescales much greater than a day, the motion and deformation of the sea ice cover (i.e., its dynamics) are primarily determined by atmospheric and oceanic tractions on its upper and lower surfaces and by internal ice forces that arise within the ice cover owing to its deformation. This review discusses the relationship between the internal ice forces and the deformation of the ice cover, focusing on representations suitable for inclusion within global climate models. I first draw attention to theories that treat the sea ice cover as an isotropic continuum and then to the recent development of anisotropic models that deal with the presence of oriented weaknesses in the ice cover, known as leads.

Effects of soil conditions and drought on egg hatching and larval survival of the clover root weevil (Sitona lepidus)

Soil-dwelling insect herbivores are significant pests in many managed ecosystems. Because eggs and larvae are difficult to observe, mathematical models have been developed to predict life-cycle events occurring in the soil. To date, these models have incorporated very little empirical information about how soil and drought conditions interact to shape these processes. This study investigated how soil temperature (10, 15, 20 and 25 °C), water content (0.02 (air dried), 0.10 and 0.25 g g−1) and pH (5, 7 and 9) interactively affected egg hatching and early larval lifespan of the clover root weevil (Sitona lepidus Gyllenhal, Coleoptera: Curculionidae). Eggs developed over 3.5 times faster at 25 °C compared with 10 °C (hatching after 40.1 and 11.5 days, respectively). The effect of drought on S. lepidus eggs was investigated by exposing eggs to drought conditions before wetting the soil (2–12 days later) at four temperatures. No eggs hatched in dry soil, suggesting that S. lepidus eggs require water to remain viable. Eggs hatched significantly sooner in slightly acidic soil (pH 5) compared with soils with higher pH values. There was also a significant interaction between soil temperature, pH and soil water content. Egg viability was significantly reduced by exposure to drought. When exposed to 2–6 days of drought, egg viability was 80–100% at all temperatures but fell to 50% after 12 days exposure at 10 °C and did not hatch at all at 20 °C and above. Drought exposure also increased hatching time of viable eggs. The effects of soil conditions on unfed larvae were less influential, except for soil temperature which significantly reduced larval longevity by 57% when reared at 25 °C compared with 10 °C (4.1 and 9.7 days, respectively). The effects of soil conditions on S. lepidus eggs and larvae are discussed in the context of global climate change and how such empirically based information could be useful for refining existing mathematical models of these processes.

Ozone database in support of CMIP5 simulations: results and corresponding radiative forcing

A continuous tropospheric and stratospheric vertically resolved ozone time series, from 1850 to 2099, has been generated to be used as forcing in global climate models that do not include interactive chemistry. A multiple linear regression analysis of SAGE I+II satellite observations and polar ozonesonde measurements is used for the stratospheric zonal mean dataset during the well-observed period from 1979 to 2009. In addition to terms describing the mean annual cycle, the regression includes terms representing equivalent effective stratospheric chlorine (EESC) and the 11-yr solar cycle variability. The EESC regression fit coefficients, together with pre-1979 EESC values, are used to extrapolate the stratospheric ozone time series backward to 1850. While a similar procedure could be used to extrapolate into the future, coupled chemistry climate model (CCM) simulations indicate that future stratospheric ozone abundances are likely to be significantly affected by climate change, and capturing such effects through a regression model approach is not feasible. Therefore, the stratospheric ozone dataset is extended into the future (merged in 2009) with multimodel mean projections from 13 CCMs that performed a simulation until 2099 under the SRES (Special Report on Emission Scenarios) A1B greenhouse gas scenario and the A1 adjusted halogen scenario in the second round of the Chemistry-Climate Model Validation (CCMVal-2) Activity. The stratospheric zonal mean ozone time series is merged with a three-dimensional tropospheric data set extracted from simulations of the past by two CCMs (CAM3.5 and GISSPUCCINI)and of the future by one CCM (CAM3.5). The future tropospheric ozone time series continues the historical CAM3.5 simulation until 2099 following the four different Representative Concentration Pathways (RCPs). Generally good agreement is found between the historical segment of the ozone database and satellite observations, although it should be noted that total column ozone is overestimated in the southern polar latitudes during spring and tropospheric column ozone is slightly underestimated. Vertical profiles of tropospheric ozone are broadly consistent with ozonesondes and in-situ measurements, with some deviations in regions of biomass burning. The tropospheric ozone radiative forcing (RF) from the 1850s to the 2000s is 0.23Wm−2, lower than previous results. The lower value is mainly due to (i) a smaller increase in biomass burning emissions; (ii) a larger influence of stratospheric ozone depletion on upper tropospheric ozone at high southern latitudes; and possibly (iii) a larger influence of clouds (which act to reduce the net forcing) compared to previous radiative forcing calculations. Over the same period, decreases in stratospheric ozone, mainly at high latitudes, produce a RF of −0.08Wm−2, which is more negative than the central Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) value of −0.05Wm−2, but which is within the stated range of −0.15 to +0.05Wm−2. The more negative value is explained by the fact that the regression model simulates significant ozone depletion prior to 1979, in line with the increase in EESC and as confirmed by CCMs, while the AR4 assumed no change in stratospheric RF prior to 1979. A negative RF of similar magnitude persists into the future, although its location shifts from high latitudes to the tropics. This shift is due to increases in polar stratospheric ozone, but decreases in tropical lower stratospheric ozone, related to a strengthening of the Brewer-Dobson circulation, particularly through the latter half of the 21st century. Differences in trends in tropospheric ozone among the four RCPs are mainly driven by different methane concentrations, resulting in a range of tropospheric ozone RFs between 0.4 and 0.1Wm−2 by 2100. The ozone dataset described here has been released for the Coupled Model Intercomparison Project (CMIP5) model simulations in netCDF Climate and Forecast (CF) Metadata Convention at the PCMDI website (http://cmip-pcmdi.llnl.gov/).

Southern annular mode dynamics in observations and models, Part II: Eddy feedbacks

Many global climate models (GCMs) have trouble simulating Southern Annular Mode (SAM) variability correctly, particularly in the Southern Hemisphere summer season where it tends to be too persistent. In this two part study, a suite of experiments with the Canadian Middle Atmosphere Model (CMAM) is analyzed to improve our understanding of the dynamics of SAM variability and its deficiencies in GCMs. Here, an examination of the eddy-mean flow feedbacks is presented by quantification of the feedback strength as a function of zonal scale and season using a new methodology that accounts for intraseasonal forcing of the SAM. In the observed atmosphere, in the summer season, a strong negative feedback by planetary scale waves, in particular zonal wavenumber 3, is found in a localized region in the south west Pacific. It cancels a large proportion of the positive feedback by synoptic and smaller scale eddies in the zonal mean, resulting in a very weak overall eddy feedback on the SAM. CMAM is deficient in this negative feedback by planetary scale waves, making a substantial contribution to its bias in summertime SAM persistence. Furthermore, this bias is not alleviated by artificially improving the climatological circulation, suggesting that climatological circulation biases are not the cause of the planetary wave feedback deficiency in the model. Analysis of the summertime eddy feedbacks in the CMIP-5 models confirms that this is indeed a common problem among GCMs, suggesting that understanding this planetary wave feedback and the reason for its deficiency in GCMs is key to improving the fidelity of simulated SAM variability in the summer season.

Influences of increasing temperature on Indian wheat: quantifying limits to predictability

As climate changes, temperatures will play an increasing role in determining crop yield. Both climate model error and lack of constrained physiological thresholds limit the predictability of yield. We used a perturbed-parameter climate model ensemble with two methods of bias-correction as input to a regional-scale wheat simulation model over India to examine future yields. This model configuration accounted for uncertainty in climate, planting date, optimization, temperature-induced changes in development rate and reproduction. It also accounts for lethal temperatures, which have been somewhat neglected to date. Using uncertainty decomposition, we found that fractional uncertainty due to temperature-driven processes in the crop model was on average larger than climate model uncertainty (0.56 versus 0.44), and that the crop model uncertainty is dominated by crop development. Simulations with the raw compared to the bias-corrected climate data did not agree on the impact on future wheat yield, nor its geographical distribution. However the method of bias-correction was not an important source of uncertainty. We conclude that bias-correction of climate model data and improved constraints on especially crop development are critical for robust impact predictions.

A methodology for in-situ and remote sensing of microphysical and radiative properties of contrails as they evolve into cirrus

Contrails and especially their evolution into cirrus-like clouds are thought to have very important effects on local and global radiation budgets, though are generally not well represented in global climate models. Lack of contrail parameterisations is due to the limited availability of in situ contrail measurements which are difficult to obtain. Here we present a methodology for successful sampling and interpretation of contrail microphysical and radiative data using both in situ and remote sensing instrumentation on board the FAAM BAe146 UK research aircraft as part of the COntrails Spreading Into Cirrus (COSIC) study.

Exploring socioeconomic impacts of forest based mitigation projects: Lessons from Brazil and Bolivia

This paper aims to contribute new insights globally and regionally on how carbon forest mitigation contributes to sustainable development in South America. Carbon finance has emerged as a potential policy option to tackling global climate change, degradation of forests, and social development in poor countries. This paper focuses on evaluating the socioeconomic impacts of a set of forest based mitigation pilot projects that emerged under the United Nations Framework Convention on Climate Change. The paper reviews research conducted in 2001–2002, drawing from empirical data from four pilot projects, derived from qualitative stakeholder interviews, and complemented by policy documents and literature. Of the four projects studied three are located in frontier areas, where there are considerable pressures for conversion of standing forest to agriculture. In this sense, forest mitigation projects have a substantial role to play in the region. Findings suggest however, that all four projects have experienced cumbersome implementation processes specifically, due to weak social objectives, poor communication, as well as time constraints. In three out of four cases, stakeholders highlighted limited local acceptance at the implementation stages. In the light of these findings, we discuss opportunities for implementation of future forest based mitigation projects in the land use sector.

Tackling drought stress: RECEPTOR-LIKE KINASES present new approaches

Global climate change and a growing population require tackling the reduction in arable land and improving biomass production and seed yield per area under varying conditions. One of these conditions is suboptimal water availability. Here, we review some of the classical approaches to dealing with plant response to drought stress and we evaluate how research on RECEPTOR-LIKE KINASES (RLKs) can contribute to improving plant performance under drought stress. RLKs are considered as key regulators of plant architecture and growth behavior, but they also function in defense and stress responses. The available literature and analyses of available transcript profiling data indeed suggest that RLKs can play an important role in optimizing plant responses to drought stress. In addition, RLK pathways are ideal targets for nontransgenic approaches, such as synthetic molecules, providing a novel strategy to manipulate their activity and supporting translational studies from model species, such as Arabidopsis thaliana, to economically useful crops.

Systematic study of the impact of fresh water fluxes on the glacial carbon cycle

During glacial periods, atmospheric CO2 concentration increases and decreases by around 15 ppm. At the same time, the climate changes gradually in Antarctica. Such climate changes can be simulated in models when the AMOC (Atlantic Meridional Oceanic Circulation) is weakened by adding fresh water to the North Atlantic. The impact on the carbon cycle is less straightforward, and previous studies give opposite results. Because the models and the fresh water fluxes were different in these studies, it prevents any direct comparison and hinders finding whether the discrepancies arise from using different models or different fresh water fluxes. In this study we use the CLIMBER-2 coupled climate carbon model to explore the impact of different fresh water fluxes. In both preindustrial and glacial states, the addition of fresh water and the resulting slow-down of the AMOC lead to an uptake of carbon by the ocean and a release by the terrestrial biosphere. The duration, shape and amplitude of the fresh water flux all have an impact on the change of atmospheric CO2 because they modulate the change of the AMOC. The maximum CO2 change linearly depends on the time integral of the AMOC change. The different duration, amplitude, and shape of the fresh water flux cannot explain the opposite evolution of ocean and vegetation carbon inventory in different models. The different CO2 evolution thus depends on the AMOC response to the addition of fresh water and the resulting climatic change, which are both model dependent. In CLIMBER-2, the rise of CO2 recorded in ice cores during abrupt events can be simulated under glacial conditions, especially when the sinking of brines in the Southern Ocean is taken into account. The addition of fresh water in the Southern Hemisphere leads to a decline of CO2, contrary to the addition of fresh water in the Northern Hemisphere.

The reversibility of sea level rise

During the last century, global climate has been warming, and projections indicate that such a warming is likely to continue over coming decades. Most of the extra heat is stored in the ocean, resulting in thermal expansion of seawater and global mean sea level rise. Previous studies have shown that after CO2 emissions cease or CO2 concentration is stabilized, global mean surface air temperature stabilizes or decreases slowly, but sea level continues to rise. Using idealized CO2 scenario simulations with a hierarchy of models including an AOGCM and a step-response model, the authors show how the evolution of thermal expansion can be interpreted in terms of the climate energy balance and the vertical profile of ocean warming. Whereas surface temperature depends on cumulative CO2 emissions, sea level rise due to thermal expansion depends on the time profile of emissions. Sea level rise is smaller for later emissions, implying that targets to limit sea level rise would need to refer to the rate of emissions, not only to the time integral. Thermal expansion is in principle reversible, but to halt or reverse it quickly requires the radiative forcing to be reduced substantially, which is possible on centennial time scales only by geoengineering. If it could be done, the results indicate that heat would leave the ocean more readily than it entered, but even if thermal expansion were returned to zero, the geographical pattern of sea level would be altered. Therefore, despite any aggressive CO2 mitigation, regional sea level change is inevitable.