Mid-Holocene and Last Glacial Maximum climate simulations with the IPSL model—part I: comparing IPSL_CM5A to IPSL_CM4

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 CO₂ 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 × CO₂ 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.     

Past and future polar amplification of climate change: climate model intercomparisons and ice-core constraints

Climate model simulations available from the PMIP1, PMIP2 and CMIP (IPCC-AR4) intercomparison projects for past and future climate change simulations are examined in terms of polar temperature changes in comparison to global temperature changes and with respect to pre-industrial reference simulations. For the mid-Holocene (MH, 6,000 years ago), the models are forced by changes in the Earth’s orbital parameters. The MH PMIP1 atmosphere-only simulations conducted with sea surface temperatures fixed to modern conditions show no MH consistent response for the poles, whereas the new PMIP2 coupled atmosphere–ocean climate models systematically simulate a significant MH warming both for Greenland (but smaller than ice-core based estimates) and Antarctica (consistent with the range of ice-core based range). In both PMIP1 and PMIP2, the MH annual mean changes in global temperature are negligible, consistent with the MH orbital forcing. The simulated last glacial maximum (LGM, 21,000 years ago) to pre-industrial change in global mean temperature ranges between 3 and 7°C in PMIP1 and PMIP2 model runs, similar to the range of temperature change expected from a quadrupling of atmospheric CO₂ concentrations in the CMIP simulations. Both LGM and future climate simulations are associated with a polar amplification of climate change. The range of glacial polar amplification in Greenland is strongly dependent on the ice sheet elevation changes prescribed to the climate models. All PMIP2 simulations systematically underestimate the reconstructed glacial–interglacial Greenland temperature change, while some of the simulations do capture the reconstructed glacial–interglacial Antarctic temperature change. Uncertainties in the prescribed central ice cap elevation cannot account for the temperature change underestimation by climate models. The variety of climate model sensitivities enables the exploration of the relative changes in polar temperature with respect to changes in global temperatures. Simulated changes of polar temperatures are strongly related to changes in simulated global temperatures for both future and LGM climates, confirming that ice-core-based reconstructions provide quantitative insights on global climate changes. An erratum to this article can be found at     

Impact of a realistic river routing in coupled ocean–atmosphere simulations of the Last Glacial Maximum climate

The presence of large ice sheets over North America and North Europe at the Last Glacial Maximum (LGM) strongly impacted Northern hemisphere river pathways. Despite the fact that such changes may significantly alter the freshwater input to the ocean, modified surface hydrology has never been accounted for in coupled ocean–atmosphere general circulation model simulations of the LGM climate. To reconstruct the LGM river routing, we use the ICE-5G LGM topography. Because of the uncertainties in the extent of the Fennoscandian ice sheet in the Eastern part of the Kara Sea, we consider two more realistic river routing scenarios. The first scenario is characterised by the presence of an ice dammed lake south of the Fennoscandian ice sheet, and corresponds to the ICE-5G topography. This lake is fed by the Ob and Yenisei rivers. In the second scenario, both these rivers flow directly into the Arctic Ocean, which is more consistent with the latest QUEEN ice sheet margin reconstructions. We study the impact of these changes on the LGM climate as simulated by the IPSL_CM4 model and focus on the overturning thermohaline circulation. A comparison with a classical LGM simulation performed using the same model and modern river basins as designed in the PMIP2 exercise leads to the following conclusions: (1) The discharge into the North Atlantic Ocean is increased by 2,000 m³/s between 38° and 54°N in both simulations that contain LGM river routing, compared to the classical LGM experiment. (2) The ice dammed lake is shown to have a weak impact, relative to the classical simulation, both in terms of climate and ocean circulation. (3) In contrast, the North Atlantic deep convection and meridional overturning are weaker than during the classical LGM run if the Ob and Yenisei rivers flow directly into the Arctic Ocean. The total discharge into the Arctic Ocean is increased by 31,000 m³/s, relative to the classical LGM simulation. Consequentially, northward ocean heat transport is weaker, and sea ice more extensive, in better agreement with existing proxy data.     

ENSO at 6ka and 21ka from ocean–atmosphere coupled model simulations

We analyze how the characteristics of El Niño-Southern Oscillation (ENSO) are changed in coupled ocean–atmosphere simulations of the mid-Holocene (MH) and the Last Glacial Maximum (LGM) performed as part of the Paleoclimate Modeling Intercomparison Project phase 2 (PMIP2). Comparison of the model results with present day observations show that most of the models reproduce the large scale features of the tropical Pacific like the SST gradient, the mean SST and the mean seasonal cycles. All models simulate the ENSO variability, although with different skill. Our analyses show that several relationships between El Niño amplitude and the mean state across the different control simulations are still valid for simulations of the MH and the LGM. Results for the MH show a consistent El Niño amplitude decrease. It can be related to the large scale atmospheric circulation changes. While the Northern Hemisphere receives more insolation during the summer time, the Asian summer monsoon system is strengthened which leads to the enhancement of the Walker circulation. Easterlies prevailing over the central eastern Pacific induce an equatorial upwelling that damps the El Niño development. Results are less conclusive for 21ka. Large scale dynamic competes with changes in local heat fluxes, so that model shows a wide range of responses, as it is the case in future climate projections.     

A comparison of PMIP2 model simulations and the MARGO proxy reconstruction for tropical sea surface temperatures at last glacial maximum

Results from multiple model simulations are used to understand the tropical sea surface temperature (SST) response to the reduced greenhouse gas concentrations and large continental ice sheets of the last glacial maximum (LGM). We present LGM simulations from the Paleoclimate Modelling Intercomparison Project, Phase 2 (PMIP2) and compare these simulations to proxy data collated and harmonized within the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface Project (MARGO). Five atmosphere–ocean coupled climate models (AOGCMs) and one coupled model of intermediate complexity have PMIP2 ocean results available for LGM. The models give a range of tropical (defined for this paper as 15°S–15°N) SST cooling of 1.0–2.4°C, comparable to the MARGO estimate of annual cooling of 1.7 ± 1°C. The models simulate greater SST cooling in the tropical Atlantic than tropical Pacific, but interbasin and intrabasin variations of cooling are much smaller than those found in the MARGO reconstruction. The simulated tropical coolings are relatively insensitive to season, a feature also present in the MARGO transferred-based estimates calculated from planktonic foraminiferal assemblages for the Indian and Pacific Oceans. These assemblages indicate seasonality in cooling in the Atlantic basin, with greater cooling in northern summer than northern winter, not captured by the model simulations. Biases in the simulations of the tropical upwelling and thermocline found in the preindustrial control simulations remain for the LGM simulations and are partly responsible for the more homogeneous spatial and temporal LGM tropical cooling simulated by the models. The PMIP2 LGM simulations give estimates for the climate sensitivity parameter of 0.67°–0.83°C per Wm⁻², which translates to equilibrium climate sensitivity for doubling of atmospheric CO₂ of 2.6–3.1°C.     

Northern hemisphere storm tracks during the last glacial maximum in the PMIP2 ocean-atmosphere coupled models: energetic study,seasonal cycle,precipitation

Mid-latitude eddies are an important component of the climatic system due to their role in transporting heat, moisture and momentum from the tropics to the poles, and also for the precipitation associated with their fronts, especially in winter. We study northern hemisphere storm-tracks at the Last Glacial Maximum (LGM) and their influence on precipitation using ocean-atmosphere general circulation model (OAGCM) simulations from the second phase of the Paleoclimate Modelling Intercomparison Project (PMIP2). The difference with PMIP1 results in terms of sea-surface temperature forcing, fundamental for storm-track dynamics, is large, especially in the eastern North Atlantic where sea-ice extends less to the south in OAGCMs compared to atmospheric-only GCMs. Our analyses of the physics of the eddies are based on the equations of eddy energetics. All models simulate a consistent southeastward shift of the North Pacific storm-track in winter, related to a similar displacement of the jet stream, partly forced by the eddies themselves. Precipitation anomalies are consistent with storm-track changes, with a southeastward displacement of the North Pacific precipitation pattern. The common features of North Atlantic changes in the LGM simulations consist of a thinning of the storm-track in its western part and an amplification of synoptic activity to the southeast, in the region between the Azores Islands and the Iberian Peninsula, which reflects on precipitation. This southeastward extension is related to a similar displacement of the jet, partly forced by the eddies. In the western North Atlantic, the synoptic activity anomalies are at first order related to baroclinic generation term anomalies, but the mean-flow baroclinicity increase due to the presence of the Laurentide ice-sheet is partly balanced by a loss of eddy efficiency to convert energy from the mean flow. Moisture availability in this region is greatly reduced due to more advection of dry polar air by stationary waves, leading to less synoptic-scale latent heat release and hence less precipitation also. In terms of seasonality, the stormy season is shifted later in the year by a few days to a month depending on the season and the model considered. This shift does not directly reflect on the first-order seasonal cycle of precipitation, which also depends on other mechanisms, especially in summer.     

A coupled climate model simulation of the Last Glacial Maximum,Part 2: approach to equilibrium

The climate of the last glacial maximum (LGM) is simulated with a coupled climate model. The simulated climate undergoes a rapid adjustment during the first several decades after imposition of LGM boundary conditions, as described in Part 1, and then evolves toward equilibrium over 900 model years. The climate simulated by the coupled model at this period is compared with observationally-based LGM reconstructions and with LGM results obtained with an atmosphere-mixed layer (slab) ocean version of the model in order to investigate the role of ocean dynamics in the LGM climate. Global mean surface air temperature and sea surface temperature (SST) decrease by about 10 °C and 5.6 °C in the coupled model which includes ocean dynamics, compared to decreases of 6.3 and 3.8 °C in slab ocean case. The coupled model simulates a cooling of about 6.5 °C over the tropics, which is larger than that of the CLIMAP reconstruction (1.7 °C) and larger than that of the slab ocean simulation (3.3 °C), but which is in reasonable agreement with some recent proxy estimates. The ocean dynamics of the coupled model captures features found in the CLIMAP reconstructions such as a relative maximum of ocean cooling over the tropical Pacific associated with a mean La Niña-like response and lead to a more realistic SST pattern than in the slab model case. The reduction in global mean precipitation simulated in the coupled model is larger (15%) than that simulated with the slab ocean model (~10%) in conjunction with the enhanced cooling. Some regions, such as the USA and the Mediterranean region, experience increased precipitation in accord with proxy paleoclimate evidence. The overall much drier climate over the ocean leads to higher sea surface salinity (SSS) in most ocean basins except for the North Atlantic where SSS is considerably lower due to an increase in the supply of fresh water from the Mississippi and Amazon rivers and presumably a decrease in salt transport by the weakened North Atlantic overturning circulation. The North Atlantic overturning stream function weakens to less than half of the control run value. The overturning is limited to a shallower depth (less than 1000 m) and its outflow is confined to the Northern Hemisphere. In the Southern Ocean, convection is much stronger than in the control run leading to a stronger overturning stream function associated with enhanced Antarctic Bottom Water formation. As a result, Southern Ocean water masses fill the entire deep ocean. The Antarctic Circumpolar Current (ACC) transport through the Drake Passage increases by about 25%. The ACC transport, despite weaker zonal winds, is enhanced due to changes in bottom pressure torque. The weakening of the overturning circulation in the North Atlantic and the accompanying 30% decrease in the poleward ocean heat transport contrasts with the strengthening of the overturning circulation in the Southern Ocean and a 40% increase in heat transport. As a result, sea ice coverage and thickness are affected in opposite senses in the two hemispheres. The LGM climate simulated by the coupled model is in reasonable agreement with paleoclimate proxy evidence. The dynamical response of the ocean in the coupled model plays an important role in determining the simulated, and undoubtedly, the actual, LGM climate.     

我国对末次冰期冰盛期东亚区域气候模拟的研究

末次冰期冰盛期（或称为末次盛冰期）,为第四纪更新世最后一个冰期的鼎盛时期,该时期的气候与当代气候迥然不同。近年来,中国科研工作者已就末次冰期冰盛期东亚区域气候开展了一些数值模拟工作。结果表明：该时期中国大陆地表气温降低,中东部地区降水显著减少,东亚冬季风增强、夏季风显著减弱。在国际古气候模拟比较计划（PMIP）标准试验的基础上,进一步指出东亚植被的反馈作用、青藏高原可能冰川的反馈作用、以及西太平洋表面温度的作用能够引起额外的气候效应,可在一定程度上改进PMIP标准试验的模拟效果。     

ENSO teleconnections in projections of future climate in ECHAM5/MPI-OM

The teleconnections of the El Niño/Southern Oscillation (ENSO) in future climate projections are investigated using results of the coupled climate model ECHAM5/MPI-OM. For this, the IPCC SRES scenario A1B and a quadrupled CO₂ simulation are considered. It is found that changes of the mean state in the tropical Pacific are likely to condition ENSO teleconnections in the Pacific North America (PNA) region and in the North Atlantic European (NAE) region. With increasing greenhouse gas emissions the changes of the mean states in the tropical and sub-tropical Pacific are El Niño-like in this particular model. Sea surface temperatures in the tropical Pacific are increased predominantly in its eastern part and redistribute the precipitation further eastward. The dynamical response of the atmosphere is such that the equatorial east–west (Walker) circulation and the eastern Pacific inverse Hadley circulation are decreased. Over the subtropical East Pacific and North Atlantic the 200 hPa westerly wind is substantially increased. Composite maps of different climate parameters for positive and negative ENSO events are used to reveal changes of the ENSO teleconnections. Mean sea level pressure and upper tropospheric zonal winds indicate an eastward shift of the well-known teleconnection patterns in the PNA region and an increasing North Atlantic oscillation (NAO) like response over the NAE region. Surface temperature and precipitation underline this effect, particularly over the North Pacific and the central North Atlantic. Moreover, in the NAE region the 200 hPa westerly wind is increasingly related to the stationary wave activity. Here the stationary waves appear NAO-like.     

Tropical paleoclimates at the Last Glacial Maximum: comparison of Paleoclimate Modeling Intercomparison Project (PMIP) simulations and paleodata

 Seventeen simulations of the Last Glacial Maximum (LGM) climate have been performed using atmospheric general circulation models (AGCM) in the framework of the Paleoclimate Modeling Intercomparison Project (PMIP). These simulations use the boundary conditions for CO₂, insolation and ice-sheets; surface temperatures (SSTs) are either (a) prescribed using CLIMAP data set (eight models) or (b) computed by coupling the AGCM with a slab ocean (nine models). The present-day (PD) tropical climate is correctly depicted by all the models, except the coarser resolution models, and the simulated geographical distribution of annual mean temperature is in good agreement with climatology. Tropical cooling at the LGM is less than at middle and high latitudes, but greatly exceeds the PD temperature variability. The LGM simulations with prescribed SSTs underestimate the observed temperature changes except over equatorial Africa where the models produce a temperature decrease consistent with the data. Our results confirm previous analyses showing that CLIMAP (1981) SSTs only produce a weak terrestrial cooling. When SSTs are computed, the models depict a cooling over the Pacific and Indian oceans in contrast with CLIMAP and most models produce cooler temperatures over land. Moreover four of the nine simulations, produce a cooling in good agreement with terrestrial data. Two of these model results over ocean are consistent with new SST reconstructions whereas two models simulate a homogeneous cooling. Finally, the LGM aridity inferred for most of the tropics from the data, is globally reproduced by the models with a strong underestimation for models using computed SSTs. Received: 9 September 1998 / Accepted: 18 March 1999     

Sensitivity of the Northern Hemisphere climate of the Last Glacial Maximum to sea surface temperatures

The role of prescribing sea surface temperature in paleoclimate atmospheric simulations has been investigated by comparing Last Glacial Maximum AGCMs experiments using different SSTs data sets as well as coupled atmosphere/oceanic mixed layer models. Changes in the SSTs and sea-ice margin generate different patterns of zonal asymmetries in the atmospheric circulation that are responsible for reorganisation of heat and moisture transport, leading to important variations of Northern Hemisphere regional climates, particularly in winter. Additional sensitivity experiments have been carried out to isolate the individual role of North Pacific and North Atlantic SSTs anomalies. We found that changes in North Pacific SSTs have a much stronger impact over all the northern continental surfaces, including Europe and Siberia, than changes in the North Atlantic SSTs. As these SSTs anomalies are of the order of the typical errors generated by coupled ocean-atmosphere models, this suggests that these more complete models will likely still have problems in simulating the regional climate change at the LGM. Received: 11 October 1999 / Accepted: 9 June 2000     

第四次国际古气候模拟比较计划（PMIP4）概况与评述

古气候模拟比较计划（PMIP）是古气候数值模拟领域一项重大的国际合作研究计划,其主旨是为古气候模拟和模拟结果评估提供一个协调机制,理解过去气候变化的物理机制和气候反馈的重要作用,为未来气候预估提供科学依据。同时,通过对比分析验证模式的模拟性能,探索其不确定性,促进耦合气候系统模式的发展。PMIP目前进行到第四阶段（PMIP4)。PMIP4进一步加强了与第六次国际耦合模式比较计划（CMIP6）的协作,选取了5组共同关注的PMIP4-CMIP6古气候模拟试验（中全新世、末次盛冰期、过去千年、末次间冰期和上新世暖期),考察气候系统对不同气候背景的综合响应。除此以外,PMIP4还设计了众多敏感性试验研究不同外强迫因子的影响。PMIP4模拟试验不仅为古气候研究提供大量的模拟数据,还将服务于CMIP6及其他众多模式比较计划。     

The Southern Annular Mode (SAM) in PMIP2 simulations of the last glacial maximum

The increasing trend of the Southern Annular Mode （SAM） in recent decades has influenced climate change in the Southem Hemisphere （SH）.How the SAM will respond increased greenhouse gas concentrations in the future remains uncertain.Understanding the variability of the SAM in the past under a colder climate such as during the Last Glacial Maximum （LGM） might provide some understanding of the response of the SAM under a future warmer climate.We analyzed the changes in the SAM during the LGM in comparison to pre-industrial （PI） simulations using five coupled ocean-atmosphere models （CCSM,FGOALS,IPSL,MIROC,HadCM） from the second phase of the Paleoclimate Modelling Intercomparison Project （PMIP2）.In CCSM,MIROC,IPSL,and FGOALS,the variability of the simulated SAM appears to be reduced in the LGM compared to the PI simulations,with a decrease in the standard deviation of the SAM index.Overall,four out of the five models suggest a weaker SAM amplitude in the LGM consistent with a weaker SH polar vortex and westerly winds found in some proxy records and model analyses.The weakening of the SAM in the LGM was associated with an increase in the vertical propagation of Rossby waves in southern high latitudes.     

Mid-Holocene tropical Pacific climate state,annual cycle,and ENSO in PMIP2 and PMIP3

Using the Paleoclimate Modeling Inter-comparison Project Phase 2 and 3 (PMIP2 and PMIP3), we investigated the tropical Pacific climate state, annual cycle, and El Niño-Southern Oscillation (ENSO) during the mid-Holocene period (6,000 years before present; 6 ka run). When the 6 ka run was compared to the control run (0 ka run), the reduced sea surface temperature (SST) and the reduced precipitation due to the basin-wide cooling, and the intensified cross-equatorial surface winds due to the hemispheric discrepancy of the surface cooling over the tropical Pacific were commonly observed in both the PMIP2 and PMIP3, but changes were more dominant in the PMIP3. The annual cycle of SST was weaker over the equatorial eastern Pacific, because of the orbital forcing change and the deepening mixed layer, while it was stronger over the equatorial western pacific in both the PMIP2 and PMIP3. The stronger annual cycle of the equatorial western Pacific SST was accompanied by the intensified annual cycle of the zonal surface wind, which dominated in the PMIP3 in particular. The ENSO activity in the 6 ka run was significantly suppressed in the PMIP2, but marginally reduced in the PMIP3. In general, the weakened air-sea coupling associated with basin-wide cooling, reduced precipitation, and a hemispheric contrast in the climate state led to the suppression of ENSO activity, and the weakening of the annual cycle over the tropical eastern Pacific might lead to the intensification of ENSO through the frequency entrainment. Therefore, the two opposite effects are slightly compensated for by each other, which results in a small reduction in the ENSO activity during the 6 ka in the PMIP3. On the whole, in PMIP2/PMIP3, the variability of canonical (or conventional) El Niño tends to be reduced during 6 ka, while that of CP/Modoki El Niño tends to be intensified.     

北大西洋年际变率的海气耦合模式模拟II：热带太平洋强迫

利用一个全球海气耦合模式(BCM),结合观测资料,讨论了热带太平洋强迫对北大西洋年际气候变率的影响。研究表明,BCM能够相对合理地模拟赤道太平洋的年际变率模态及相应的海温距平型和大气遥相关型,尽管其准3年的振荡周期过于规则。来自数值模式和观测上的证据都表明,北大西洋冬季海温的主导性变率模态,即自北而南出现的“- -”的海温距平型,受到来自热带太平洋强迫的显著影响,其正位相与赤道中东太平洋冷事件相对应。换言之,赤道太平洋暖事件的发生,在太平洋-北美沿岸激发出PNA遥相关型,进而通过在北大西洋产生类似NAO负位相的气压距平型,削弱本来与NAO正位相直接联系的三核型海温距平。北大西洋三核型海温距平对热带太平洋强迫的响应,要滞后2—3个月的时间。     

An overview of the performance of CMIP6 models in the tropical Atlantic: mean state,variability, and remote impacts

General circulation models of the Coupled Model Intercomparison Project Phase 6 (CMIP6) are examined with respect to their ability to simulate the mean state and variability of the tropical Atlantic and its linkage to the tropical Pacific. While, on average, mean state biases have improved little, relative to the previous intercomparison (CMIP5), there are now a few models with very small biases. In particular the equatorial Atlantic warm SST and westerly wind biases are mostly eliminated in these models. Furthermore, interannual variability in the equatorial and subtropical Atlantic is quite realistic in a number of CMIP6 models, which suggests that they should be useful tools for understanding and predicting variability patterns. The evolution of equatorial Atlantic biases follows the same pattern as in previous model generations, with westerly wind biases during boreal spring preceding warm sea-surface temperature (SST) biases in the east during boreal summer. A substantial portion of the westerly wind bias exists already in atmosphere-only simulations forced with observed SST, suggesting an atmospheric origin. While variability is relatively realistic in many models, SSTs seem less responsive to wind forcing than observed, both on the equator and in the subtropics, possibly due to an excessively deep mixed layer originating in the oceanic component. Thus models with realistic SST amplitude tend to have excessive wind amplitude. The models with the smallest mean state biases all have relatively high resolution but there are also a few low-resolution models that perform similarly well, indicating that resolution is not the only way toward reducing tropical Atlantic biases. The results also show a relatively weak link between mean state biases and the quality of the simulated variability. The linkage to the tropical Pacific shows a wide range of behaviors across models, indicating the need for further model improvement.

     

Origin of regional climate differences: role of boundary conditions and model formulation in two GCMs

Model differences in projections of extratropical regional climate change due to increasing greenhouse gases are investigated using two atmospheric general circulation models (AGCMs): ECHAM4 (Max Planck Institute, version 4) and CCM3 (National Center for Atmospheric Research Community Climate Model version 3). Sea-surface temperature (SST) fields calculated from observations and coupled versions of the two models are used to force each AGCM in experiments based on time-slice methodology. Results from the forced AGCMs are then compared to coupled model results from the Coupled Model Intercomparison Project 2 (CMIP2) database. The time-slice methodology is verified by showing that the response of each model to doubled CO₂ and SST forcing from the CMIP2 experiments is consistent with the results of the coupled GCMs. The differences in the responses of the models are attributed to (1) the different tropical SST warmings in the coupled simulations and (2) the different atmospheric model responses to the same tropical SST warmings. Both are found to have important contributions to differences in implied Northern Hemisphere (NH) winter extratropical regional 500 mb height and tropical precipitation climate changes. Forced teleconnection patterns from tropical SST differences are primarily responsible for sensitivity differences in the extratropical North Pacific, but have relatively little impact on the North Atlantic. There are also significant differences in the extratropical response of the models to the same tropical SST anomalies due to differences in numerical and physical parameterizations. Differences due to parameterizations dominate in the North Atlantic. Differences in the control climates of the two coupled models from the current climate, in particular for the coupled model containing CCM3, are also demonstrated to be important in leading to differences in extratropical regional sensitivity.     

Synoptic-scale perturbations in AGCM simulations of the present and Last Glacial Maximum climates

 The conditions of development of mid-latitude depressions (synoptic eddies) in the winter Northern Hemisphere mid-latitudes at the Last Glacial Maximum (LGM, 21 000 years ago) are very different from the present ones: this period is characterised by a general cooling of the extra-tropics, with massive ice sheets over the Northern Hemisphere continents and sea-ice extending very far south over the North Atlantic. The present work uses regression analysis to study the characteristics of the synoptic eddies in present-day and LGM climate simulations by the Atmospheric General Circulation Model (AGCM) of the UK Universities' Global Atmospheric Programme (UGAMP). In the LGM experiment, the structure of the Pacific eddies is similar to the present-day (PD) situation, but they are weaker. On the other hand, the Atlantic eddies show an increased zonal wavelength and a much shallower structure in the temperature and vertical wind perturbations. To understand the changes of these characteristics from present-day to LGM, we compare them to those computed for the most unstable modes of the corresponding mean flows, determined using a dry primitive equation model. A normal-mode stability analysis is carried both on zonally symmetric and asymmetric flows for each of the Northern Hemisphere storm-tracks. The changes in the most unstable normal modes found by both these analyses give a good account of changes in the structure of the perturbations as retrieved from the AGCM, suggesting that changes in the mean state (especially the temperature gradient) is the main driver of these changes. However in the case of the present-day Atlantic storm-track, the growth rate of these modes is found to be very low compared to the other cases. A complementary analysis evaluates the importance of non-modal growth, in the form of downstream development of perturbations, for each of the storm-tracks. This type of growth is found to be especially important in the case of the present-day Atlantic storm-track. Received: 29 September 1999 / Accepted: 17 November 1999     

Weak ENSO Asymmetry Due to Weak Nonlinear Air–Sea Interaction in CMIP5 Climate Models

State-of-the-art climate models have long-standing intrinsic biases that limit their simulation and projection capabilities.Significantly weak ENSO asymmetry and weakly nonlinear air–sea interaction over the tropical Pacific was found in CMIP5(Coupled Model Intercomparison Project, Phase 5) climate models compared with observation. The results suggest that a weak nonlinear air–sea interaction may play a role in the weak ENSO asymmetry. Moreover, a weak nonlinearity in air–sea interaction in the models may be associated with the biases in the mean climate—the cold biases in the equatorial central Pacific. The excessive cold tongue bias pushes the deep convection far west to the western Pacific warm pool region and suppresses its development in the central equatorial Pacific. The deep convection has difficulties in further moving to the eastern equatorial Pacific, especially during extreme El Ni o events, which confines the westerly wind anomaly to the western Pacific. This weakens the eastern Pacific El Ni o events, especially the extreme El Ni o events, and thus leads to the weakened ENSO asymmetry in climate models. An accurate mean state structure(especially a realistic cold tongue and deep convection) is critical to reproducing ENSO events in climate models. Our evaluation also revealed that ENSO statistics in CMIP5 climate models are slightly improved compared with those of CMIP3. The weak ENSO asymmetry in CMIP5 is closer to the observation. It is more evident in CMIP5 that strong ENSO activities are usually accompanied by strong ENSO asymmetry, and the diversity of ENSO amplitude is reduced.     

Structure and variances of equatorial zonal circulation in a multimodel ensemble

The structure and variance of the equatorial zonal circulation, as characterized by the atmospheric mass flux in the equatorial zonal plane, is examined and inter-compared in simulations from 9 CMIP3 coupled climate models with multiple ensemble members and the NCEP-NCAR and ERA-40 reanalyses. The climate model simulations analyzed here include twentieth century (20C3M) and twenty-first century (SRES A1B) simulations. We evaluate the 20C3M modeled zonal circulations by comparing them with those in the reanalyses. We then examine the variability of the circulation, its changes with global warming, and the associated thermodynamic maintenance. The tropical zonal circulation involves three major components situated over the Pacific, Indian, and Atlantic oceans. The three cells are supported by the corresponding diabatic heating extending deeply throughout the troposphere, with heating centers apparent in the mid-troposphere. Seasonal features appear in the zonal circulation, including variations in its intensity and longitudinal migration. Most models, and hence the multi-model mean, represent the annual and seasonal features of the circulation and the associated heating reasonably well. The multi-model mean reproduces the observed climatology better than any individual model, as indicated by the spatial pattern correlation and mean square difference of the mass flux and the diabatic heating compared to the reanalysis based values. Projected changes in the zonal circulation under A1B forcing are dominated by mass flux changes over the Pacific and Indian oceans. An eastward shift of the Pacific Walker circulation is clearly evident with global warming, with anomalous rising motion apparent over the equatorial central Pacific and anomalous sinking motions in the west and east, which favors an overall strengthening of the Walker circulation. The zonal circulation weakens and shifts westwards over the Indian Ocean under external forcing, whereas it strengthens and shifts slightly westwards over the Atlantic Ocean. The forced circulation changes are associated with broad SST and atmospheric diabatic heating changes in the tropics. Linear trends of these forced circulation changes, as characterized by regional spatial maximum amplitudes of mass fluxes and their longitudes over the three oceans, are statistically significant at the 5?% level for 2000–2099 for the multi-model mean. However, wide differences of the trends are apparent across the models, because of both deficiencies in the simulation of the circulations in different models and the high internal variability of the circulations.