Observational challenges in evaluating climate models

Activities like the Coupled Model Intercomparison Project (CMIP) have revolutionized climate modelling in terms of our ability to compare models and to process information about climate projections and their uncertainties. The evaluation of models against observations is now considered a key component of multi-model studies. While there are a number of outstanding scientific issues surrounding model evaluation, notably the open question of how to link model performance to future projections, here we highlight a specific but growing problem in model evaluation - that of uncertainties in the observational data that are used to evaluate the models. We highlight the problem using an example obtained from studies of the South Asian Monsoon but we believe the problem is a generic one which arises in many different areas of climate model evaluation and which requires some attention by the community.

The depression of tropical snow lines at the Last Glacial Maximum: what can we learn from climate and model experiments?

Analyses of simulations of the last glacial maximum (LGM) made with 17 atmospheric general circulation models (AGCMs) participating in the Paleoclimate Modelling Intercomparison Project, and a high-resolution (T106) version of one of the models (CCSR1), show that changes in the elevation of tropical snowlines (as estimated by the depression of the maximum altitude of the 0 °C isotherm) are primarily controlled by changes in sea-surface temperatures (SSTs). The correlation between the two variables, averaged for the tropics as a whole, is 95%, and remains >80% even at a regional scale. The reduction of tropical SSTs at the LGM results in a drier atmosphere and hence steeper lapse rates. Changes in atmospheric circulation patterns, particularly the weakening of the Asian monsoon system and related atmospheric humidity changes, amplify the reduction in snowline elevation in the northern tropics. Colder conditions over the tropical oceans combined with a weakened Asian monsoon could produce snowline lowering of up to 1000 m in certain regions, comparable to the changes shown by observations. Nevertheless, such large changes are not typical of all regions of the tropics. Analysis of the higher resolution CCSR1 simulation shows that differences between the free atmospheric and along-slope lapse rate can be large, and may provide an additional factor to explain regional variations in observed snowline changes.

Global monsoons in the mid-Holocene and oceanic feedback

The response of the six major summer monsoon systems (the North American monsoon, the northern Africa monsoon, the Asia monsoon, the northern Australasian monsoon, the South America monsoon and the southern Africa monsoon) to mid-Holocene orbital forcing has been investigated using a coupled ocean–atmosphere general circulation model (FOAM), with the focus on the distinct roles of the direct insolation forcing and oceanic feedback. The simulation result is also found to compare well with the NCAR CSM. The direct effects of the change in insolation produce an enhancement of the Northern Hemisphere monsoons and a reduction of the Southern Hemisphere monsoons. Ocean feedbacks produce a further enhancement of the northern Africa monsoon and the North American monsoon. However, ocean feedbacks appear to weaken the Asia monsoon, although the overall effect (direct insolation forcing plus ocean feedback) remains a strengthened monsoon. The impact of ocean feedbacks on the South American and southern African monsoons is relatively small, and therefore these regions, especially the South America, experienced a reduced monsoon regime compared to present. However, there is a strong ocean feedback on the northern Australian monsoon that negates the direct effects of orbital changes and results in a strengthening of austral summer monsoon precipitation in this region. A new synthesis is made for mid-Holocene paleoenvironmental records and is compared with the model simulations. Overall, model simulations produce changes in regional climates that are generally consistent with paleoenvironmental observations.

Mid-Holocene climates of the Americas: a dynamical response to changed seasonality

Simulations of the climatic response to mid-Holocene (6 ka BP) orbital forcing with two coupled ocean–atmosphere models (FOAM and CSM) show enhancement of monsoonal precipitation in parts of the American Southwest, Central America and northernmost South America during Northern Hemisphere summer. The enhanced onshore flow that brings precipitation into Central America is caused by a northward displacement of the inter-tropical convergence zone, driven by cooling of the equatorial and warming of the northern subtropical and mid-latitude ocean. Ocean feedbacks also enhance precipitation over the American Southwest, although the increase in monsoon precipitation there is largely driven by increases in land-surface temperature. The northward shift in the equatorial precipitation band that causes enhanced precipitation in Central America and the American Southwest has a negative feedback effect on monsoonal precipitation in northern South America. The simulations demonstrate that mid-Holocene aridity in the mid-continent of North America is dynamically linked to the orbitally induced enhancement of the summer monsoon in the American Southwest, with a spatial structure (wet in the Southwest and dry in the mid-continent) similar to that found in strong monsoon years today. Changes in winter precipitation along the west coast of North America, in Central America and along the Gulf Coast, caused by southward-displacement of the westerly storm tracks, indicate that changes in the Northern Hemisphere winter monsoon also play a role in regional climate changes during the mid-Holocene. Although the simulations with FOAM and CSM differ in detail, the general mechanisms and patterns are common to both. The model results thus provide a coherent dynamical explanation for regional patterns of increased or decreased aridity shown by vegetation, lake status and aeolian data from the Americas

Pollen-based biome reconstructions for China at 0 and 6000 years

Biomization provides an objective and robust method of assigning pollen spectra to biomes so that pollen data can be mapped and compared directly with the output of biomgeographic models. We have tested the applicability of this procedure, originally developed for Europe, to assign modern surface samples from China to biomes. The procedure successfully delineated the major vegetation types of China. When the same procedure was applied to fossil pollen samples for 6000 years ago, the reconstructions showed systematic differences from present, consistent with previous interpretations of vegetation changes since the mid-Holocene. In eastern China, the forest zones were systematically shifted northwards, such that cool mixed forests displaced taiga in northeastern China, while broad-leaved evergreen forest extended c. 300 km and temperate deciduous forestc. 500–600 km beyond their present northern limits. In northwestern China, the area of desert and steppe vegetation was reduced compared to present. On the Tibetan Plateau, forest vegetation extended to higher elevations than today and the area of tundra was reduced. These shifts in biome distributions imply significant changes in climate since 6000 years ago that can be interpreted qualitatively as a response to orbital forcing and its secondary effects on the Asian monsoon.

Climatology of Upper Tropospheric–Lower Stratospheric (UTLS) Jets and Tropopauses in MERRA

A global climatology (1979–2012) from the Modern-Era Retrospective Analysis for Research and Applications (MERRA) shows distributions and seasonal evolution of upper tropospheric jets and their relationships to the stratospheric subvortex and multiple tropopauses. The overall climatological patterns of upper tropospheric jets confirm those seen in previous studies, indicating accurate representation of jet stream dynamics in MERRA. The analysis shows a Northern Hemisphere (NH) upper tropospheric jet stretching nearly zonally from the mid-Atlantic across Africa and Asia. In winter–spring, this jet splits over the eastern Pacific, merges again over eastern North America, and then shifts poleward over the North Atlantic. The jets associated with tropical circulations are also captured, with upper tropospheric westerlies demarking cyclonic flow downstream from the Australian and Asian monsoon anticyclones and associated easterly jets. Multiple tropopauses associated with the thermal tropopause “break” commonly extend poleward from the subtropical upper tropospheric jet. In Southern Hemisphere (SH) summer, the tropopause break, along with a poleward-stretching secondary tropopause, often occurs across the tropical westerly jet downstream of the Australian monsoon region. SH high-latitude multiple tropopauses, nearly ubiquitous in June–July, are associated with the unique polar winter thermal structure. High-latitude multiple tropopauses in NH fall–winter are, however, sometimes associated with poleward-shifted upper tropospheric jets. The SH subvortex jet extends down near the level of the subtropical jet core in winter and spring. Most SH subvortex jets merge with an upper tropospheric jet between May and December; although much less persistent than in the SH, merged NH subvortex jets are common between November and April.

Late Quaternary lake-level record from northern Eurasia

Lake records from northern Eurasia show regionally coherent patterns of changes during the late Quaternary. Lakes peripheral to the Scandinavian ice sheet were lower than those today but lakes in the Mediterranean zone were high at the glacial maximum, reflecting the dominance of glacial anticyclonic conditions in northern Europe and a southward shift of the Westerlies. The influence of the glacial anticyclonic circulation attenuated through the late glacial period, and the Westerlies gradually shifted northward, such that drier conditions south of the ice sheet were confined to a progressively narrower zone and the Mediterranean became drier. The early Holocene shows a gradual shift to conditions wetter than present in central Asia, associated with the expanded Asian monsoon, and in the Mediterranean, in response to local, monsoon-type circulation. There is no evidence of mid-continental aridity in northern Eurasia during the mid-Holocene. In contrast, the circum-Baltic region was drier, reflecting the increased incidence of blocking anticyclones centered on Scandinavia in summer. There is a gradual transition to modern conditions after ca. 5000 yr B.P. Although these broad-scale patterns are interrupted by shorter term fluctuations, the long-term trends in lake behavior show a clear response to changes in insolation and glaciation.

Mid-holocene and last glacial maximum climate simulations with the IPSL model-part I: comparing IPSLCM5A to IPSLCM4

The climates of the mid-Holocene (MH), 6,000 years ago, and of the Last Glacial Maximum (LGM), 21,000 years ago, have extensively been simulated, in particular in the framework of the Palaeoclimate Modelling Intercomparion Project. These periods are well documented by paleo-records, which can be used for evaluating model results for climates different from the present one. Here, we present new simulations of the MH and the LGM climates obtained with the IPSL_CM5A model and compare them to our previous results obtained with the IPSL_CM4 model. Compared to IPSL_CM4, IPSL_CM5A includes two new features: the interactive representation of the plant phenology and marine biogeochemistry. But one of the most important differences between these models is the latitudinal resolution and vertical domain of their atmospheric component, which have been improved in IPSL_CM5A and results in a better representation of the mid-latitude jet-streams. The Asian monsoon’s representation is also substantially improved. The global average mean annual temperature simulated for the pre-industrial (PI) period is colder in IPSL_CM5A than in IPSL_CM4 but their climate sensitivity to a CO2 doubling is similar. Here we show that these differences in the simulated PI climate have an impact on the simulated MH and LGM climatic anomalies. The larger cooling response to LGM boundary conditions in IPSL_CM5A appears to be mainly due to differences between the PMIP3 and PMIP2 boundary conditions, as shown by a short wave radiative forcing/feedback analysis based on a simplified perturbation method. It is found that the sensitivity computed from the LGM climate is lower than that computed from 2 × CO2 simulations, confirming previous studies based on different models. For the MH, the Asian monsoon, stronger in the IPSL_CM5A PI simulation, is also more sensitive to the insolation changes. The African monsoon is also further amplified in IPSL_CM5A due to the impact of the interactive phenology. Finally the changes in variability for both models and for MH and LGM are presented taking the example of the El-Niño Southern Oscillation (ENSO), which is very different in the PI simulations. ENSO variability is damped in both model versions at the MH, whereas inconsistent responses are found between the two versions for the LGM. Part 2 of this paper examines whether these differences between IPSL_CM4 and IPSL_CM5A can be distinguished when comparing those results to palaeo-climatic reconstructions and investigates new approaches for model-data comparisons made possible by the inclusion of new components in IPSL_CM5A.

Late Quaternary environmental change in the interior South American tropics: new insight from leaf wax stable isotopes

Stable isotope analysis of leaf waxes in a sediment core from Laguna La Gaiba, a shallow lake located at the Bolivian margin of the Pantanal wetlands, provides new perspective on vegetation and climate change in the lowland interior tropics of South America over the past 40,000 years. The carbon isotopic compositions (δ13C) of long-chain n-alkanes reveal large shifts between C3-and C4-dominated vegetation communities since the last glacial period, consistent with landscape reconstructions generated with pollen data from the same sediment core. Leaf wax δ13C values during the last glacial period reflect an open landscape composed of C4grasses and C3herbs from 41–20ka. A peak in C4abundance during the Last Glacial Maximum (LGM, ∼21ka) suggests drier or more seasonal conditions relative to the earlier glacial period, while the development of a C3-dominated forest community after 20 ka points to increased humidity during the last deglaciation. Within the Holocene, large changes in the abundance of C4 vegetation indicate a transition from drier or more seasonal conditions during the early/mid-Holocene to wetter conditions in the late Holocene coincident with increasing austral summer insolation. Strong negative correlations between leaf wax δ13C and δD values over the entire record indicate that the majority of variability in leaf wax δD at this site can be explained by variability in the magnitude of biosynthetic fractionation by different vegetation types rather than changes in meteoric water δD signatures. However, positive δD deviations from the observed δ13C–δD trends are consistent with more enriched source water and drier or more seasonal conditions during the early/mid-Holocene and LGM. Overall, our record adds to evidence of varying influence of glacial boundary conditions and orbital forcing on South American Summer Monsoon precipitation in different regions of the South American tropics. Moreover, the relationships between leaf wax stable isotopes and pollen data observed at this site underscore the complementary nature of pollen and leaf wax δ13C data for reconstructing past vegetation changes and the potentially large effects of such changes on leaf wax δD signatures.

Seasonality of westerly moisture transport in the East Asian summer monsoon and its implications for interpreting precipitation δ18O

East Asian summer monsoon (EASM) rainfall impacts the world's most populous regions. Accurate EASM rainfall prediction necessitates robust paleoclimate reconstructions from proxy data and quantitative linkage to modern climatic conditions. Many precisely dated oxygen isotope records from Chinese stalagmites have been interpreted as directly reflecting past EASM rainfall amount variability, but recent research suggests that such records instead integrate multiple hydroclimatic processes. Using a Lagrangian precipitation moisture source diagnostic, we demonstrate that EASM rainfall is primarily derived from the Indian Ocean. Conversely, Pacific Ocean moisture export peaks during winter, and the moisture uptake area does not differ significantly between summer and winter and is thus a minor contributor to monsoonal precipitation. Our results are substantiated by an accurate reproduction of summer and winter spatial rainfall distributions across China. We also correlate modern EASM rainfall oxygen isotope ratios with instrumental rainfall amount and our moisture source data. This analysis reveals that the strength of the source effect is geographically variable, and differences in atmospheric moisture transport may significantly impact the isotopic signature of EASM rainfall at the Hulu, Dongge, and Wanxiang Cave sites. These results improve our ability to isolate the rainfall amount signal in paleomonsoon reconstructions and indicate that precipitation across central and eastern China will directly respond to variability in Indian Ocean moisture supply.

The mean evolution and variability of the Asian summer monsoon: Comparison of ECMWF and NCEP-NCAR reanalyses

The behavior of the Asian summer monsoon is documented and compared using the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA) and the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) Reanalysis. In terms of seasonal mean climatologies the results suggest that, in several respects, the ERA is superior to the NCEP-NCAR Reanalysis. The overall better simulation of the precipitation and hence the diabatic heating field over the monsoon domain in ERA means that the analyzed circulation is probably nearer reality. In terms of interannual variability, inconsistencies in the definition of weak and strong monsoon years based on typical monsoon indices such as All-India Rainfall (AIR) anomalies and the large-scale wind shear based dynamical monsoon index (DMI) still exist. Two dominant modes of interannual variability have been identified that together explain nearly 50% of the variance. Individually, they have many features in common with the composite flow patterns associated with weak and strong monsoons, when defined in terms of regional AIR anomalies and the large-scale DMI. The reanalyses also show a common dominant mode of intraseasonal variability that describes the latitudinal displacement of the tropical convergence zone from its oceanic-to-continental regime and essentially captures the low-frequency active/break cycles of the monsoon. The relationship between interannual and intraseasonal variability has been investigated by considering the probability density function (PDF) of the principal component of the dominant intraseasonal mode. Based on the DMI, there is an indication that in years with a weaker monsoon circulation, the PDF is skewed toward negative values (i,e., break conditions). Similarly, the PDFs for El Nino and La Nina years suggest that El Nino predisposes the system to more break spells, although the sample size may limit the statistical significance of the results.

The effect of regional changes in anthropogenic aerosols on rainfall of the East Asian summer monsoon

The response of East Asian Summer Monsoon (EASM) precipitation to long term changes in regional anthropogenic aerosols (sulphate and black carbon) is explored in an atmospheric general circulation model, the atmospheric component of the UK High-Resolution Global Environment Model v1.2 (HiGAM). Separately, sulphur dioxide (SO2) and black carbon (BC) emissions in 1950 and 2000 over East Asia are used to drive model simulations, while emissions are kept constant at year 2000 level outside this region. The response of the EASM is examined by comparing simulations driven by aerosol emissions representative of 1950 and 2000. The aerosol radiative effects are also determined using an off-line radiative transfer model. During June, July and August, the EASM was not significantly changed as either SO2 or BC emissions increased from 1950 to 2000 levels. However, in September, precipitation is significantly decreased by 26.4% for sulphate aerosol and 14.6% for black carbon when emissions are at the 2000 level. Over 80% of the decrease is attributed to changes in convective precipitation. The cooler land surface temperature over China in September (0.8 °C for sulphate and 0.5 °C for black carbon) due to increased aerosols reduces the surface thermal contrast that supports the EASM circulation. However, mechanisms causing the surface temperature decrease in September are different between sulphate and BC experiments. In the sulphate experiment, the sulphate direct and the 1st indirect radiative effects contribute to the surface cooling. In the BC experiment, the BC direct effect is the main driver of the surface cooling, however, a decrease in low cloud cover due to the increased heating by BC absorption partially counteracts the direct effect. This results in a weaker land surface temperature response to BC changes than to sulphate changes. The resulting precipitation response is also weaker, and the responses of the monsoon circulation are different for sulphate and black carbon experiments. This study demonstrates a mechanism that links regional aerosol emission changes to the precipitation changes of the EASM, and it could be applied to help understand the future changes in EASM precipitation in CMIP5 simulations.

An investigation of QBO signals in the east Asian and Indian monsoon in GCM experiments

Observations have shown that the monsoon is a highly variable phenomenon of the tropical troposphere, which exhibits significant variance in the temporal range of two to three years. The reason for this specific interannual variability has not yet been identified unequivocally. Observational analyses have also shown that EI Niño indices or western Pacific SSTs exhibit some power in the two to three year period range and therefore it was suggested that an ocean-atmosphere interaction could excite and support such a cycle. Similar mechanisms include land-surface-atmosphere interaction as a possible driving mechanism. A rather different explanation could be provided by a forcing mechanism based on the quasi-biennial oscillation of the zonal wind in the lower equatorial stratosphere (QBO). The QBO is a phenomenon driven by equatorial waves with periods of some days which are excited in the troposphere. Provided that the monsoon circulation reacts to the modulation of tropopause conditions as forced by the QBO, this could explain monsoon variability in the quasi-biennial window. The possibility of a QBO-driven monsoon variability is investigated in this study in a number of general circulation model experiments where the QBO is assimilated to externally controlled phase states. These experiments show that the boreal summer monsoon is significantly influenced by the QBO. A QBO westerly phase implies less precipitation in the western Pacific, but more in India, in agreement with observations. The austral summer monsoon is exposed to similar but weaker mechanisms and the precipitation does not change significantly.

Preferred response of the East Asian summer monsoon to local and non-local anthropogenic sulphur dioxide emissions

In this study, the atmospheric component of a state-of-the-art climate model (HadGEM2-ES) that includes earth system components such as interactive chemistry and eight species of tropospheric aerosols considering aerosol direct, indirect, and semi-direct effects, has been used to investigate the impacts of local and non-local emissions of anthropogenic sulphur dioxide on the East Asian summer monsoon (EASM). The study focuses on the fast responses (including land surface feedbacks, but without sea surface temperature feedbacks) to sudden changes in emissions from Asia and Europe. The initial responses, over days 1–40, to Asian and European emissions show large differences. The response to Asian emissions involves a direct impact on the sulphate burden over Asia, with immediate consequences for the shortwave energy budget through aerosol–radiation and aerosol–cloud interactions. These changes lead to cooling of East Asia and a weakening of the EASM. In contrast, European emissions have no significant impact on the sulphate burden over Asia, but they induce mid-tropospheric cooling and drying over the European sector. Subsequently, however, this cold and dry anomaly is advected into Asia, where it induces atmospheric and surface feedbacks over Asia and the Western North Pacific (WNP), which also weaken the EASM. In spite of very different perturbations to the local aerosol burden in response to Asian and European sulphur dioxide emissions, the large scale pattern of changes in land–sea thermal contrast, atmospheric circulation and local precipitation over East Asia from days 40 onward exhibits similar structures, indicating a preferred response, and suggesting that emissions from both regions likely contributed to the observed weakening of the EASM. Cooling and drying of the troposphere over Asia, together with warming and moistening over the WNP, reduces the land–sea thermal contrast between the Asian continent and surrounding oceans. This leads to high sea level pressure (SLP) anomalies over Asia and low SLP anomalies over the WNP, associated with a weakened EASM. In response to emissions from both regions warming and moistening over the WNP plays an important role and determines the time scale of the response.

The effect of doubled CO2 and model basic state biases on the monsoon-ENSO system. I: Mean response and interannual variability

The impact of doubled CO2 concentration on the Asian summer monsoon is studied using a coupled ocean-atmosphere model. Both the mean seasonal precipitation and interannual monsoon variability are found to increase in the future climate scenario presented. Systematic biases in current climate simulations of the coupled system prevent accurate representation of the monsoon-ENSO teleconnection, of prime importance for seasonal prediction and for determining monsoon interannual variability. By applying seasonally varying heat flux adjustments to the tropical Pacific and Indian Ocean surface in the future climate simulation, some assessment can be made of the impact of systematic model biases on future climate predictions. In simulations where the flux adjustments are implemented, the response to climate change is magnified, with the suggestion that systematic biases may be masking the true impact of increased greenhouse gas forcing. The teleconnection between ENSO and the Asian summer monsoon remains robust in the future climate, although the Indo-Pacific takes on more of a biennial character for long periods of the flux-adjusted simulation. Assessing the teleconnection across interdecadal timescales shows wide variations in its amplitude, despite the absence of external forcing. This suggests that recent changes in the observed record cannot be distinguished from internal variations and as such are not necessarily related to climate change.