The influence of continental ice,atmospheric CO2, and land albedo on the climate of the last glacial maximum

The contributions of expanded continental ice, reduced atmospheric CO₂, and changes in land albedo to the maintenance of the climate of the last glacial maximum (LGM) are examined. A series of experiments is performed using an atmosphere-mixed layer ocean model in which these changes in boundary conditions are incorporated either singly or in combination. The model used has been shown to produce a reasonably realistic simulation of the reduced temperature of the LGM (Manabe and Broccoli 1985b). By comparing the results from pairs of experiments, the effects of each of these environmental changes can be determined.Expanded continental ice and reduced atmospheric CO₂ are found to have a substantial impact on global mean temperature. The ice sheet effect is confined almost exclusively to the Northern Hemisphere, while lowered CO₂ cools both hemispheres. Changes in land albedo over ice-free areas have only a minor thermal effect on a global basis. The reduction of CO₂ content in the atmosphere is the primary contributor to the cooling of the Southern Hemisphere. The model sensitivity to both the ice sheet and CO₂ effects is characterized by a high latitude amplification and a late autumn and early winter maximum.Substantial changes in Northern Hemisphere tropospheric circulation are found in response to LGM boundary conditions during winter. An amplified flow pattern and enhanced westerlies occur in the vicinity of the North American and Eurasian ice sheets. These alterations of the tropospheric circulation are primarily the result of the ice sheet effect, with reduced CO₂ contributing only a slight amplification of the ice sheet-induced pattern.     

The effect of atmospheric CO2 and ice sheet topography on LGM climate

The role of reduced atmospheric CO₂ concentration and ice sheet topography plus its associated land albedo on the LGM climate is investigated using a coupled atmosphere-ocean-sea ice climate system model. The surface cooling induced by the reduced CO₂ concentration is larger than that by the ice sheet topography plus other factors by about 30% for the surface air temperature and by about 100% for the sea surface temperature. A large inter-hemispheric asymmetry in surface cooling with a larger cooling in the Northern Hemisphere is found for both cases. This asymmetric inter-hemispheric temperature response is consistent in the ice sheet topography case with earlier studies using an atmospheric model coupled with a mixed-layer ocean representation, but contrasts with these results in the reduced CO₂ case. The incorporation of ocean dynamics presumably leads to a larger snow and sea ice feedback as a result of the reduction in northward ocean heat transport, mainly as a consequence of the decrease in the North Atlantic overturning circulation by the substantial freshening of the North Atlantic convection regions. A reversed case is found in the Southern Ocean. Overall, the reduction in atmospheric CO₂ concentration accounts for about 60% of the total LGM climate change.     

Relative roles of dynamic and thermodynamic processes in causing positive and negative global mean SST trends during the past 100 years

Global mean surface temperature (GMST) during 1910–2012 experienced four alternated rapid warming and warming hiatus phases. Such a temporal variation is primarily determined by global mean sea surface temperature (SST) component. The relative roles of ocean dynamic and thermodynamic processes in causing such global mean SST variations are investigated, using two methods. The first method is ocean mixed layer heat budget analysis. The budget diagnosis result shows that the thermodynamic processes dominate in the rapid warming phases, while the ocean dynamics dominate during the hiatus phases. The second method relies on the diagnosis of a simple equilibrium state model. This model captures well the horizontal distribution of SST difference between two warmer and cooler equilibrium states during either the rapid warming or hiatus phases. It is found that the SST difference during the rapid warming phases is primarily controlled by the increase of downward longwave radiation as both column integrated water vapor and CO₂ increase during the phases. During the hiatus phases, the water vapor induced greenhouse effect offsets the CO₂ effect, and the SST cooling tendency is primarily determined by the ocean dynamics over the Southern Ocean and tropical Pacific. The SST pattern associated with the Interdecadal Pacific Oscillation (IPO) might be responsible for the remote and local ocean dynamic responses through induced wind change.     

Oceanic transport and storage of carbon emissions

Using a global carbon cycle model (GLOCO) that considers seven terrestrial biomes, surface and deep ocean layers based on the HILDA model and a single mixed atmosphere, we analyzed the response of atmospheric CO₂ concentration and oceanic DIC and DOC depth profiles to additions of carbon to the atmosphere and ocean. The rate of transport of carbon to the deepest oceanic layers is rather insensitive to the atmosphereic-ocean surface gas exchange coefficient over a wide range, hence discrepancies between researchers on the precise global average value of this coefficient do not significantly affect predictions of atmospheric response to anthropogenic inputs. Upwelling velocity, on the other hand, amplifies oceanic response by increasing primary production in the upper ocean layers, resulting in a larger flux into DOC and sediments and increased carbon storage; experiments to reduce the uncertainty in this parameter would be valuable.The location of the carbon addition, whether it is released in the atmosphere or in the middle of the oceanic thermocline, has a significant impact on the maximum atmospheric CO₂ concentration (pCO₂) subsequently reached, suggesting that oceanic burial of a significant fraction of carbon emissions (e.g. via clathrate hydrides) may be an important management option for limiting pCO₂ buildup. Our analysis indicates that the effectiveness of ocean burial decreases asymptotically below about 1000 m depth. With a constant emissions scenario (at 1990 levels), pCO₂ at year 2100 is reduced from 501 ppmv considering all emissions go to the atmosphere, to 422 ppmv with ocean burial at a depth of 1000 m of 50% of the fossil fuel emissions. An alternative scenario looks at stabilizing pCO₂ at 450 ppmv; with no ocean burial of fossil fuel emissions, the rate of emissions has to be cut drastically after the year 2010, whereas oceanic burial of 2 GtC/yr allows for a smoother transition to alternative energy sources.     

The role of temperature on treeline migration for an eastern African mountain during the Last Glacial Maximum

Paleo-data suggest that East African mountain treelines underwent an altitudinal shift during the Last Glacial Maximum (LGM). Understanding the ecological and physiological processes underlying treeline response to such past climate change will help to improve forecasts of treeline change under future global warming. In spite of significant improvements in paleoclimatic reconstruction, the climatic conditions explaining this migration are still debated and important factors such as atmospheric CO₂ concentration, the impact of lapse rate decreasing temperature along altitudinal gradients and rainfall modifications due to elevation have often been neglected or simplified. Here, we assess the effects of these different factors and estimate the influence of the most dominant factors controlling changes in past treeline position using a multi-proxy approach based on simulations from BIOME4, a coupled biogeography and biogeochemistry model, modified to account for the effect of elevation on vegetation, compared with pollen, and isotopic data. The results indicate a shift in mountain vegetation at the LGM was controlled by low pCO₂ and low temperatures promoting species morphologically and physiologically better adapted to LGM conditions than many trees composing the forest belt limit. Our estimate that the LGM climate was cooler than today’s by ?4.5 °C (range: ?4.3 to ?4.6 °C) at the upper limit of the treeline, whereas at 831 m it was cooler by ?1.4 °C (range: ?2.6 to ?0.6 °C), suggests that a possible lapse rate modification strongly constrained the upper limit of treeline, which may limit its potential extension under future global warming.     

Cloud processes associated with past and future climate changes

 To investigate the cloud response during cold and warm periods, we have performed simulations of the Last Glacial Maximum (LGM-21ky BP) and of double CO₂ concentration using the LMD AGCM model. We observe that the thermal characteristics of these two climates are opposite, but the cloud response is more complex and does not display the same symmetry When doubling the CO₂, the warming of the troposphere and the cooling of the stratosphere are clearly linked with a reduction in low-level clouds and an increase of high-level clouds associated with relative humidity changes. For the LGM, the cloud response is more complex. In the inter tropical region, we show that the Hadley cell is reinforced during LGM (+20%) whereas it is reduced (−10%) for the double CO₂ experiments. The most important feature is that we observe an enlarged Hadley cell for LGM climate which strongly modifies the atmospheric dynamics and water transport. For LGM conditions, the cloud response is then mostly driven by these dynamical changes at low latitudes though at high latitudes the thermal changes explain a large part of the cloud response. Two different versions of the model, using different parametrizations for the precipitation show that cloud feedbacks may act differently for cold and warm climates; and that the cloud response may be more complex that previously expected, but also indicate that the details of these effects are model dependent.     

Changes of interannual NAO variability in response to greenhouse gases forcing

Observations show that there was change in interannual North Atlantic Oscillation (NAO) variability in the mid-1970s. This change was characterized by an eastward shift of the NAO action centres, a poleward shift of zonal wind anomalies and a downstream extension of climate anomalies associated with the NAO. The NAO interannual variability for the period after the mid-1970s has an annular mode structure that penetrates deeply into the stratosphere, indicating a strengthened relationship between the NAO and the Arctic Oscillation (AO) and strengthened stratosphere-troposphere coupling. In this study we have investigated possible causes of these changes in the NAO by carrying out experiments with an atmospheric GCM. The model is forced either by doubling CO₂, or increasing sea surface temperatures (SST), or both. In the case of SST forcing the SST anomaly is derived from a coupled model simulation forced by increasing CO₂. Results indicate that SST and CO₂ change both force a poleward and eastward shift in the pattern of interannual NAO variability and the associated poleward shift of zonal wind anomalies, similar to the observations. The effect of SST change can be understood in terms of mean changes in the troposphere. The direct effect of CO₂ change, in contrast, can not be understood in terms of mean changes in the troposphere. However, there is a significant response in the stratosphere, characterized by a strengthened climatological polar vortex with strongly enhanced interannual variability. In this case, the NAO interannual variability has a strong link with the variability over the North Pacific, as in the annular AO pattern, and is also strongly related to the stratospheric vortex, indicating strengthened stratosphere-troposphere coupling. The similarity of changes in many characteristics of NAO interannual variability between the model response to doubling CO₂ and those in observations in the mid-1970s implies that the increase of greenhouse gas concentration in the atmosphere, and the resulting changes in the stratosphere, might have played an important role in the multidecadal change of interannual NAO variability and its associated climate anomalies during the late twentieth century. The weak change in mean westerlies in the troposphere in response to CO₂ change implies that enhanced and eastward extended mid-latitude westerlies in the troposphere might not be a necessary condition for the poleward and eastward shift of the NAO action centres in the mid-1970s.     

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.     

Wind-driven changes in Southern Ocean residual circulation, ocean carbon reservoirs and atmospheric CO2

The effect of idealized wind-driven circulation changes in the Southern Ocean on atmospheric CO₂ and the ocean carbon inventory is investigated using a suite of coarse-resolution, global coupled ocean circulation and biogeochemistry experiments with parameterized eddy activity and only modest changes in surface buoyancy forcing, each experiment integrated for 5,000 years. A positive correlation is obtained between the meridional overturning or residual circulation in the Southern Ocean and atmospheric CO₂: stronger or northward-shifted westerly winds in the Southern Hemisphere result in increased residual circulation, greater upwelling of carbon-rich deep waters and oceanic outgassing, which increases atmospheric pCO₂ by ～20 μatm; weaker or southward-shifted winds lead to the opposing result. The ocean carbon inventory in our model varies through contrasting changes in the saturated, disequilibrium and biogenic (soft-tissue and carbonate) reservoirs, each varying by O(10–100) PgC, all of which contribute to the net anomaly in atmospheric CO₂. Increased residual overturning deepens the global pycnocline, warming the upper ocean and decreasing the saturated carbon reservoir. Increased upwelling of carbon- and nutrient-rich deep waters and inefficient biological activity results in subduction of unutilized nutrients into the ocean interior, decreasing the biogenic carbon reservoir of intermediate and mode waters ventilating the Northern Hemisphere, and making the disequilibrium carbon reservoir more positive in the mode waters due to the reduced residence time at the surface. Wind-induced changes in the model carbon inventory are dominated by the response of the global pycnocline, although there is an additional abyssal response when the peak westerly winds change their latitude, altering their proximity to Drake Passage and changing the depth extent of the southward return flow of the overturning: a northward shift of the westerly winds isolates dense isopycnals, allowing biogenic carbon to accumulate in the deep ocean of the Southern Hemisphere, while a southward shift shoals dense isopycnals that outcrop in the Southern Ocean and reduces the biogenic carbon store in the deep ocean.     

Summertime land–sea thermal contrast and atmospheric circulation over East Asia in a warming climate—Part II: Importance of CO2-induced continental warming

In this the second of a two-part study, we examine the physical mechanisms responsible for the increasing contrast of the land–sea surface air temperature (SAT) in summertime over the Far East, as observed in recent decades and revealed in future climate projections obtained from a series of transient warming and sensitivity experiments conducted under the umbrella of the Coupled Model Intercomparison Project phase 5. On a global perspective, a strengthening of land–sea SAT contrast in the transient warming simulations of coupled atmosphere–ocean general circulation models is attributed to an increase in sea surface temperature (SST). However, in boreal summer, the strengthened contrast over the Far East is reproduced only by increasing atmospheric CO₂ concentration. In response to SST increase alone, the tropospheric warming over the interior of the mid- to high-latitude continents including Eurasia are weaker than those over the surrounding oceans, leading to a weakening of the land–sea SAT contrast over the Far East. Thus, the increasing contrast and associated change in atmospheric circulation over East Asia is explained by CO₂-induced continental warming. The degree of strengthening of the land–sea SAT contrast varies in different transient warming scenarios, but is reproduced through a combination of the CO₂-induced positive and SST-induced negative contributions to the land–sea contrast. These results imply that changes of climate patterns over the land–ocean boundary regions are sensitive to future scenarios of CO₂ concentration pathways including extreme cases.     

Influences of the NAO on the North Atlantic CO_2 Fluxes in Winter and Summer on the Interannual Scale

The differences in the influences of the North Atlantic Oscillation (NAO) on the air–sea CO_2 fluxes (f CO_2) in the North Atlantic (NA) between different seasons and between different regions are rarely fully investigated. We used observation-based data of f CO_2, surface-ocean CO_2partial pressure (p CO_(2sea)), wind speed and sea surface temperature(SST) to analyze the relationship between the NAO and f CO_2 of the subtropical and subpolar NA in winter and summer on the interannual time scale. Based on power spectrum estimation, there are significant interannual signs with a 2–6 year cycle in the NAO indexes and area-averaged f CO_2 anomalies in winter and summer from 1980 to 2015. Regression analysis with the 2–6 year filtered data shows that on the interannual scale the response of the f CO_2 anomalies to the NAO has an obvious meridional wave-train-like pattern in winter, but a zonal distribution in summer. This seasonal difference is because in winter the f CO_2anomalies are mainly controlled by the NAO-driven wind speed anomalies, which have a meridional distribution pattern, while in summer they are dominated by the NAO-driven SST anomalies, which show distinct zonal difference in the subtropical NA. In addition, in the same season, there are different factors controlling the variation of p CO_(2sea)in different regions. In summer, SST is important to the interannual variation of p CO_(2sea)in the subtropical NA, while some biogeochemical variables probably control the p CO_(2sea) variation in the subpolar NA.     

Strengthened African summer monsoon in the mid-Piacenzian

Using model results from the first phase of the Pliocene Model Intercomparison Project (PlioMIP) and four experiments with CAM4, the intensified African summer monsoon (ASM) in the mid-Piacenzian and corresponding mechanisms are analyzed. The results from PlioMIP show that the ASM intensified and summer precipitation increased in North Africa during the mid-Piacenzian, which can be explained by the increased net energy in the atmospheric column above North Africa. Further experiments with CAM4 indicated that the combined changes in the mid-Piacenzian of atmospheric CO₂ concentration and SST, as well as the vegetation change, could have substantially increased the net energy in the atmospheric column over North Africa and further intensified the ASM. The experiments also demonstrated that topography change had a weak effect. Overall, the combined changes of atmospheric CO₂ concentration and SST were the most important factor that brought about the intensified ASM in the mid-Piacenzian.     

Factors influencing anthropogenic carbon dioxide uptake in the North Atlantic in models of the ocean carbon cycle

The uptake and storage of anthropogenic carbon in the North Atlantic is investigated using different configurations of ocean general circulation/carbon cycle models. We investigate how different representations of the ocean physics in the models, which represent the range of models currently in use, affect the evolution of CO₂ uptake in the North Atlantic. The buffer effect of the ocean carbon system would be expected to reduce ocean CO₂ uptake as the ocean absorbs increasing amounts of CO₂. We find that the strength of the buffer effect is very dependent on the model ocean state, as it affects both the magnitude and timing of the changes in uptake. The timescale over which uptake of CO₂ in the North Atlantic drops to below preindustrial levels is particularly sensitive to the ocean state which sets the degree of buffering; it is less sensitive to the choice of atmospheric CO₂ forcing scenario. Neglecting physical climate change effects, North Atlantic CO₂ uptake drops below preindustrial levels between 50 and 300 years after stabilisation of atmospheric CO₂ in different model configurations. Storage of anthropogenic carbon in the North Atlantic varies much less among the different model configurations, as differences in ocean transport of dissolved inorganic carbon and uptake of CO₂ compensate each other. This supports the idea that measured inventories of anthropogenic carbon in the real ocean cannot be used to constrain the surface uptake. Including physical climate change effects reduces anthropogenic CO₂ uptake and storage in the North Atlantic further, due to the combined effects of surface warming, increased freshwater input, and a slowdown of the meridional overturning circulation. The timescale over which North Atlantic CO₂ uptake drops to below preindustrial levels is reduced by about one-third, leading to an estimate of this timescale for the real world of about 50 years after the stabilisation of atmospheric CO₂. In the climate change experiment, a shallowing of the mixed layer depths in the North Atlantic results in a significant reduction in primary production, reducing the potential role for biology in drawing down anthropogenic CO₂.     

Transport and storage of CO2 in the ocean ——an inorganic ocean-circulation carbon cycle model

Inorganic carbon in the ocean is modelled as a passive tracer advected by a three-dimensional current field computed from a dynamical global ocean circulation model. The carbon exchange between the ocean and atmosphere is determined directly from the (temperature-dependent) chemical interaction rates in the mixed layer, using a standard CO₂ flux relation at the air-sea interface. The carbon cycle is closed by coupling the ocean to a one-layer, horizontally diffusive atmosphere. Biological sources and sinks are not included. In this form the ocean carbon model contains essentially no free tuning parameters. The model may be regarded as a reference for interpreting numerical experiments with extended versions of the model including biological processes in the ocean (Bacastow R and Maier-Reimer E in prep.) and on land (Esser G et al in prep.). Qualitatively, the model reproduces the principal features of the observed CO₂ distribution bution in the surface ocean. However, the amplitudes of surface pCO₂ are underestimated in upwelling regions by a factor of the order of 1.5 due to the missing biological pump. The model without biota may, nevertheless, be applied to compute the storage capacity of the ocean to first order for anthropogenic CO₂ emissions. In the linear regime, the response of the model may be represented by an impulse response function which can be approximated by a superposition of exponentials with different amplitudes and time constants. This provides a simple reference for comparison with box models. The largest-amplitude (0.35) exponential has a time constant of 300 years. The effective storage capacity of the oceans is strongly dependent on the time history of the anthropogenic input, as found also in earlier box model studies.     

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     

Radiative forcing and response of a GCM to ice age boundary conditions: cloud feedback and climate sensitivity

 A general circulation model is used to examine the effects of reduced atmospheric CO₂, insolation changes and an updated reconstruction of the continental ice sheets at the Last Glacial Maximum (LGM). A set of experiments is performed to estimate the radiative forcing from each of the boundary conditions. These calculations are used to estimate a total radiative forcing for the climate of the LGM. The response of the general circulation model to the forcing from each of the changed boundary conditions is then investigated. About two-thirds of the simulated glacial cooling is due to the presence of the continental ice sheets. The effect of the cloud feedback is substantially modified where there are large changes to surface albedo. Finally, the climate sensitivity is estimated based on the global mean LGM radiative forcing and temperature response, and is compared to the climate sensitivity calculated from equilibrium experiments with atmospheric CO₂ doubled from present day concentration. The calculations here using the model and palaeodata support a climate sensitivity of about 1 Wm^-2 K^-1 which is within the conventional range. Received: 8 February 1997 / Accepted: 4 June 1997     

Mechanisms for projected future changes in south Asian monsoon precipitation

Results are first presented from an analysis of a global coupled climate model regarding changes in future mean and variability of south Asian monsoon precipitation due to increased atmospheric CO₂ for doubled (2 × CO₂) and quadrupled (4 × CO₂) present-day amounts. Results from the coupled model show that, in agreement with previous studies, mean area-averaged south Asian monsoon precipitation increases with greater CO₂ concentrations, as does the interannual variability. Mechanisms producing these changes are then examined in a series of AMIP2-style sensitivity experiments using the atmospheric model (taken from the coupled model) run with specified SSTs. Three sets of ensemble experiments are run with SST anomalies superimposed on the AMIP2 SSTs from 1979–97: (1) anomalously warm Indian Ocean SSTs, (2) anomalously warm Pacific Ocean SSTs, and (3) anomalously warm Indian and Pacific Ocean SSTs. Results from these experiments show that the greater mean monsoon precipitation is due to increased moisture source from the warmer Indian Ocean. Increased south Asian monsoon interannual variability is primarily due to warmer Pacific Ocean SSTs with enhanced evaporation variability, with the warmer Indian Ocean SSTs a contributing but secondary factor. That is, for a given interannual tropical Pacific SST fluctuation with warmer mean SSTs in the future climate, there is enhanced evaporation and precipitation variability that is communicated via the Walker Circulation in the atmosphere to the south Asian monsoon to increase interannual precipitation variability there. This enhanced monsoon variability occurs even with no change in interannual SST variability in the tropical Pacific.     

Effects of tropical instability wave (TIW)-induced surface wind feedback in the tropical Pacific Ocean

Tropical instability waves (TIWs) arise from oceanic instability in the eastern tropical Pacific and Atlantic Oceans, having a clear atmospheric signature that results in coupled atmosphere–ocean interactions at TIW scales. In this study, the extent to which TIW-induced surface wind feedback influences the ocean is examined using an ocean general circulation model (OGCM). The TIW-induced wind stress (τ_TIW) part is diagnostically determined using an empirical τ_TIW model from sea surface temperature (SST) fields simulated in the OGCM. The interactively represented TIW wind tends to reduce TIW activity in the ocean and influence the mean state, with largest impacts during TIW active periods in fall and winter. In December, the interactive τ_TIW forcing induces a surface cooling (an order of ?0.1 to ?0.3 °C), an increased heat flux into the ocean, a shallower mixed layer and a weakening of the South Equatorial Current in the eastern equatorial Pacific. Additionally, the TIW wind effect yields a pronounced latitudinal asymmetry of sea level field across the equator, and a change to upper thermal structure, characterized by a surface cooling and a warming below in the thermocline, leading to a decreased temperature gradient between the mixed layer and the thermocline. Processes responsible for the τ_TIW–induced cooling effects are analyzed. Vertical mixing and meridional advection are the two terms in the SST budget that are dominantly affected by the TIW wind feedback: the cooling effect from the vertical mixing on SST is enhanced, with the maximum induced cooling in winter; the warming effect from the meridional advection is reduced in July–October, but enhanced in November–December. Additional experiments are performed to separate the relative roles the affected surface momentum and heat fluxes play in the cooling effect on SST. This ocean-only modeling work indicates that the effect of TIW-induced wind feedback is small but not negligible, and may need to be adequately taken into account in large-scale climate modeling.     

Seasonally asymmetric transition of the Asian monsoon in response to ice age boundary conditions

Modulation of a monsoon under glacial forcing is examined using an atmosphere?Cocean coupled general circulation model (AOGCM) following the specifications established by Paleoclimate Modelling Intercomparison Project phase 2 (PMIP2) to understand the air?Csea?Cland interaction under different climate forcing. Several sensitivity experiments are performed in response to individual changes in the continental ice sheet, orbital parameters, and sea surface temperature (SST) in the Last Glacial Maximum (LGM: 21?ka) to evaluate the driving mechanisms for the anomalous seasonal evolution of the monsoon. Comparison of the model results in the LGM with the pre-industrial (PI) simulation shows that the Arabian Sea and Bay of Bengal are characterized by enhancement of pre-monsoon convection despite a drop in the SST encompassing the globe, while the rainfall is considerably suppressed in the subsequent monsoon period. In the LGM winter relative to the PI, anomalies in the meridional temperature gradient (MTG) between the Asian continents minus the tropical oceans become positive and are consistent with the intensified pre-monsoon circulation. The enhanced MTG anomalies can be explained by a decrease in the condensation heating relevant to the suppressed tropical convection as well as positive insolation anomalies in the higher latitude, showing an opposing view to a warmer future climate. It is also evident that a latitudinal gradient in the SST across the equator plays an important role in the enhancement of pre-monsoon rainfall. As for the summer, the sensitivity experiments imply that two ice sheets over the northern hemisphere cools the air temperature over the Asian continent, which is consistent with the reduction of MTG involved in the attenuated monsoon. The surplus pre-monsoon convection causes a decrease in the SST through increased heat loss from the ocean surface; in other words, negative ocean feedback is also responsible for the subsequent weakening of summer convection.     

Towards quantifying uncertainty in transient climate change

Ensembles of coupled atmosphere–ocean global circulation model simulations are required to make probabilistic predictions of future climate change. “Perturbed physics” ensembles provide a new approach in which modelling uncertainties are sampled systematically by perturbing uncertain parameters. The aim is to provide a basis for probabilistic predictions in which the impact of prior assumptions and observational constraints can be clearly distinguished. Here we report on the first perturbed physics coupled atmosphere–ocean model ensemble in which poorly constrained atmosphere, land and sea-ice component parameters are varied in the third version of the Hadley Centre model (the variation of ocean parameters will be the subject of future study). Flux adjustments are employed, both to reduce regional sea surface temperature (SST) and salinity biases and also to admit the use of combinations of model parameter values which give non-zero values for the global radiation balance. This improves the extent to which the ensemble provides a credible basis for the quantification of uncertainties in climate change, especially at a regional level. However, this particular implementation of flux-adjustments leads to a weakening of the Atlantic overturning circulation, resulting in the development of biases in SST and sea ice in the North Atlantic and Arctic Oceans. Nevertheless, model versions are produced which are of similar quality to the unperturbed and un-flux-adjusted version. The ensemble is used to simulate pre-industrial conditions and a simple scenario of a 1% per year compounded increase in CO₂. The range of transient climate response (the 20 year averaged global warming at the time of CO₂ doubling) is 1.5–2.6°C, similar to that found in multi-model studies. Measures of global and large scale climate change from the coupled models show simple relationships with associated measures computed from atmosphere-mixed-layer-ocean climate change experiments, suggesting that recent advances in computing the probability density function of climate change under equilibrium conditions using the perturbed physics approach may be extended to the transient case.