Publications

2010
We show that intermodel variations in the anthropogenically-forced evolution of September sea ice extent (SSIE) in the Arctic stem mainly from two factors: the baseline climatological sea ice thickness (SIT) distribution, and the local climate feedback parameter. The roles of these two factors evolve over the course of the twenty-first century. The SIT distribution is the most important factor in current trends and those of coming decades, accounting for roughly half the intermodel variations in SSIE trends. Then, its role progressively decreases, so that around the middle of the twenty-first century the local climate feedback parameter becomes the dominant factor. Through this analysis, we identify the investments in improved simulation of Arctic climate necessary to reduce uncertainties both in projections of sea ice loss over the coming decades and in the ultimate fate of the ice pack.
Boé, J, A Hall, F Colas, JC McWilliams, X Qu, J Kurian, and S Kapnick. 2010. “What shapes mesoscale wind anomalies in coastal upwelling zones?” Climate Dynamics 36: 2037–2049. Publisher's Version Abstract
Observational studies have shown that mesoscale variations in sea surface temperature may induce mesoscale variations in wind. In eastern subtropical upwelling regions such as the California coast, this mechanism could be of great importance for the mean state and variability of the climate system. In coastal regions orography also creates mesoscale variations in wind, and the orographic effect may extend more than 100 km offshore. The respective roles of SST/wind links and coastal orography in shaping mesoscale wind variations in nearshore regions is not clear. We address this question in the context of the California Upwelling System, using a high-resolution regional numerical modeling system coupling the WRF atmospheric model to the ROMS oceanic model, as well as additional uncoupled experiments to quantify and separate the effects of SST/wind links and coastal orography on mesoscale wind variations. After taking into account potential biases in the representation of the strength of SST/wind links by the model, our results suggest that the magnitude of mesoscale wind variations arising from the orographic effects is roughly twice that of wind variations associated with mesoscale SST anomalies. This indicates that even in this region where coastal orography is complex and leaves a strong imprint on coastal winds, the role of SST/winds links in shaping coastal circulation and climate cannot be neglected.
Qu, X, A Hall, and J Boé. 2010. “Why does the Antarctic Peninsula warm in climate simulations?” Climate Dynamics 38 (5–6): 913–927. Publisher's Version Abstract
The Antarctic Peninsula has warmed significantly since the 1950s. This pronounced and isolated warming trend is collectively captured by 29 twentieth-century climate hindcasts participating in the version 3 Coupled Model Intercomparison Project. To understand the factors driving warming trends in the hindcasts, we examine trends in Peninsula region’s atmospheric heat budget in every simulation. We find that atmospheric latent heat release increases in nearly all hindcasts. These increases are generally anthropogenic in origin, and account for about 60% of the ensemble-mean warming trend in the Peninsula. They are driven primarily by well-understood features of the anthropogenic intensification of global hydrological cycle. As sea surface temperature increases, moisture contained in atmospheric flows increases. When such flows are forced to ascend the Peninsula’s topography, enhanced local latent heat release results. The mechanism driving the warming of the Antarctic Peninsula is therefore clear in the models. Evidence for a similar mechanism operating in the real world is seen in the increasing snow accumulation rates inferred from ice cores drilled in the Peninsula. However, the relative importance of this mechanism and other processes previously identified as potentially causing the observed warming, such as the recent sea ice retreat in the Bellingshausen Sea, is difficult to assess. Thus the relevance of the simulated warming mechanism to the observed warming is unclear, in spite of its robustness in the models.
2009
Hughes, M, A Hall, and RG Fovell. 2009. “Blocking in areas of complex topography, and its influence on rainfall distribution.” Journal of the Atmospheric Sciences 66: 508–518. Publisher's Version Abstract

Using a 6-km-resolution regional climate simulation of Southern California, the effect of orographic blocking on the precipitation climatology is examined. To diagnose whether blocking occurs, precipitating hours are categorized by a bulk Froude number. The precipitation distribution becomes much more spatially homogeneous as the Froude number decreases, and an inspection of winds confirms that this results from the increasing prevalence of orographic blocking. Low Froude (Froude approximately less than 1), blocked cases account for a large fraction of climatological precipitation, particularly at the coastline where more than half is attributable to blocked cases. Thus, the climatological precipitation–slope relationship seen in observations and in the simulation is a hybrid of blocked and unblocked cases.

Simulated precipitation distributions are compared to those predicted by a simple linear model that includes only rainfall arising from direct forced topographic ascent. The agreement is nearly perfect for high Froude (Froude substantially larger than 1) cases but degrades dramatically as the index decreases; as blocking becomes more prevalent, the precipitation–slope relationship becomes continuously weaker than that predicted by the linear model. Because of its high fidelity during unblocked cases, it is surmised that blocking effects are the primary limitation preventing the linear model from accurately representing precipitation climatology and that the representation would be significantly improved during low Froude hours by the addition of a term to reduce the effective slope of the topography. These results suggest orographic blocking may substantially affect climatological precipitation distributions in similarly configured coastal areas.

Boé, J, A Hall, and X Qu. 2009. “Current GCMs' unrealistic negative feedback in the Arctic.” . Journal of Climate 22: 4682–4695. Publisher's Version Abstract
The large spread of the response to anthropogenic forcing simulated by state-of-the-art climate models in the Arctic is investigated. A feedback analysis framework specific to the Arctic is developed to address this issue. The feedback analysis shows that a large part of the spread of Arctic climate change is explained by the longwave feedback parameter. The large spread of the negative longwave feedback parameter is in turn mainly due to variations in temperature feedback. The vertical temperature structure of the atmosphere in the Arctic, characterized by a surface inversion during wintertime, exerts a strong control on the temperature feedback and consequently on simulated Arctic climate change. Most current climate models likely overestimate the climatological strength of the inversion, leading to excessive negative longwave feedback. The authors conclude that the models’ near-equilibrium response to anthropogenic forcing is generally too small.
Fletcher, C, P Kushner, A Hall, and X Qu. 2009. “Circulation responses to snow albedo feedback in climate change.” Geophysical Research Letters 36: L09702. Publisher's Version Abstract
Climate change is expected to cause a reduction in the spatial extent of snow cover on land. Recent work suggests that this will exert a local influence on the atmosphere and the hydrology of snow‐margin areas through the snow‐albedo feedback (SAF) mechanism. A significant fraction of variability among IPCC AR4 general circulation model (GCM) predictions for future summertime climate change over these areas is related to the models' representation of springtime SAF. In this study, we demonstrate a nonlocal influence of SAF on the summertime circulation in the extratropical Northern Hemisphere. Increased land surface warming in models with stronger SAF is associated with large‐scale sea‐level pressure anomalies over the northern oceans and a poleward intensified subtropical jet. We find that up to 25–30% and, on average, 5–10% of the inter‐model spread in projections of the circulation response to climate change is linearly related to SAF strength.
Derevianko, G, C Deutsch, and A Hall. 2009. “On the relationship between DMS and solar radiation.” Geophysical Research Letters 36: L17606. Publisher's Version Abstract
Biologically produced dimethylsulfide (DMS) is an important source of sulfur to the marine atmosphere that may affect cloud formation and properties. DMS is involved in a complex set of biochemical transformations and ecological exchanges so its global distribution is influenced by numerous factors, including oxidative stress from UV radiation. We re‐examine correlations between global surface DMS concentrations and mixed layer solar radiation dose (SRD), and find that SRD accounts for only a very small fraction (14%) of total variance in DMS measurements when using minimal aggregation methods. Moreover this relationship arises in part from the fact that when mixed layers deepen, both SRD and DMS decrease. When we control for this confounding effect, the correlation between DMS and SRD is reduced even further. These results indicate that factors other than solar irradiance play a leading role in determining global DMS emissions.
Boé, J, A Hall, and X Qu. 2009. “September sea-ice cover in the Arctic Ocean projected to vanish by 2100.” Nature Geoscience 2: 341–343. Publisher's Version Abstract
The Arctic climate is changing rapidly1. From 1979 to 2006, September sea-ice extent decreased by almost 25% or about 100,000 km2 per year (ref. 2). In September 2007, Arctic sea-ice extent reached its lowest level since satellite observations began3and in September 2008, sea-ice cover was still low. This development has raised concerns that the Arctic Ocean could be ice-free in late summer in only a few decades, with important economic and geopolitical implications. Unfortunately, most current climate models underestimate significantly the observed trend in Arctic sea-ice decline4, leading to doubts regarding their projections for the timing of ice-free conditions. Here we analyse the simulated trends in past sea-ice cover in 18 state-of-art-climate models and find a direct relationship between the simulated evolution of September sea-ice cover over the twenty-first century and the magnitude of past trends in sea-ice cover. Using this relationship together with observed trends, we project the evolution of September sea-ice cover over the twenty-first century. We find that under a scenario with medium future greenhouse-gas emissions, the Arctic Ocean will probably be ice-free in September before the end of the twenty-first century.
Boé, J, A Hall, and X Qu. 2009. “Deep ocean heat uptake as a major source of spread in transient climate change simulations.” Geophysical Research Letters 36: L22701. Publisher's Version Abstract

Two main mechanisms can potentially explain the spread in the magnitude of global warming simulated by climate models: deep ocean heat uptake and climate feedbacks. Here, we show that deep oceanic heat uptake is a major source of spread in simulations of 21st century climate change. Models with deeper baseline polar mixed layers are associated with larger deep ocean warming and smaller global surface warming. Based on this result, we set forth an observational constraint on polar vertical oceanic mixing. This constraint suggests that many models may overestimate the efficiency of polar oceanic mixing and therefore may underestimate future surface warming. Thus to reduce climate change uncertainties at time‐scales relevant for policy‐making, improved understanding and modelling of oceanic mixing at high latitudes is crucial.

Fernandes, R, H Zhao, X Wang, J Key, X Qu, and A Hall. 2009. “Controls on northern hemisphere snow albedo feedback quantified using satelllite Earth observations.” Geophysical Research Letters 36: L21702. Publisher's Version Abstract
Observation based estimates of controls on snow albedo feedback (SAF) are needed to constrain the snow and albedo parameterizations in general circulation model (GCM) projections of air temperature over the Northern Hemisphere (NH) landmass. The total April‐May NH SAF, corresponding to the sum of the effect of temperature on surface albedo over snow covered surfaces (‘metamorphism’) and over surfaces transitioning from snow covered to snow free conditions (‘snow cover’), is derived with daily NH snow cover and surface albedo products using Advanced Very High Resolution Radiometer Polar Pathfinder satellite data and surface air temperature from ERA40 reanalysis data between 1982–1999. Without using snow cover information, the estimated total SAF, for land surfaces north of 30°N, of −0.93 ± 0.06%K−1 was not significantly different (95% confidence) from estimates based on International Satellite Cloud Climatology Project surface albedo data. The SAF, constrained to only snow covered areas, grew to −1.06 ± 0.08%K−1 with similar magnitudes for the ‘snow cover’ and ‘metamorphosis’ components. The SAF pattern was significantly correlated with the ‘snow cover’ component pattern over both North America and Eurasia but only over Eurasia for the ‘metamorphosis’ component. However, in contrast to GCM model based diagnoses of SAF, the control on the ‘snow cover’ component related to the albedo contrast of snow covered and snow free surfaces was not strongly correlated to the total SAF.
2008
Hall, A, X Qu, and JD Neelin. 2008. “Improving predictions of summer climate change in the United States.” Geophysical Research Letters 35: L01702. Publisher's Version Abstract
Across vast, agriculturally intensive regions of the United States, the spread in predictions of summer temperature and soil moisture under global warming is curiously elevated in current climate models. Some models show modest warming of 2–3C° and little drying or slight moistening by the 22nd century, while at the other extreme are simulations with warming as large as 7–8C° and 20–40% reductions in soil moisture. We show this region of large spread arises from differences in simulations of snow albedo feedback. During winter and early spring, models with strong snow albedo feedback exhibit large reductions in snowpack and hence water storage. This water deficit persists in summer soil moisture, with reduced evapotranspiration yielding warmer temperatures. Comparison of simulated feedback strength to observations of the feedback from the current climate's seasonal cycle suggests the inter‐model differences are excessive. At the same time, the multi‐model mean feedback strength agrees reasonably well with the observed value. We estimate that if the next generation of models were brought into line with observations of snow albedo feedback, the unusually wide divergence in simulations of summer warming and drying over the US would shrink by roughly one third to one half.
2007
Qu, X, and A Hall. 2007. “What controls the strength of snow albedo feedback?” Journal of Climate 20: 3971–39. Publisher's Version Abstract
The strength of snow-albedo feedback (SAF) in transient climate change simulations of the Fourth Assessment of the Intergovernmental Panel on Climate Change is generally determined by the surface-albedo decrease associated with a loss of snow cover rather than the reduction in snow albedo due to snow metamorphosis in a warming climate. The large intermodel spread in SAF strength is likewise attributable mostly to the snow cover component. The spread in the strength of this component is in turn mostly attributable to a correspondingly large spread in mean effective snow albedo. Models with large effective snow albedos have a large surface-albedo contrast between snow-covered and snow-free regions and exhibit a correspondingly large surface-albedo decrease when snow cover decreases. Models without explicit treatment of the vegetation canopy in their surface-albedo calculations typically have high effective snow albedos and strong SAF, often stronger than observed. In models with explicit canopy treatment, completely snow-covered surfaces typically have lower albedos and the simulations have weaker SAF, generally weaker than observed. The authors speculate that in these models either snow albedos or canopy albedos when snow is present are too low, or vegetation shields snow-covered surfaces excessively. Detailed observations of surface albedo in a representative sampling of snow-covered surfaces would therefore be extremely useful in constraining these parameterizations and reducing SAF spread in the next generation of models.
Hughes, M, A Hall, and RG Fovell. 2007. “Dynamical controls on the diurnal cycle of temperature in complex topography.” Climate Dynamics 29: 277–292. Publisher's Version Abstract
We examine the climatological diurnal cycle of surface air temperature in a 6 km resolution atmospheric simulation of Southern California from 1995 to the present. We find its amplitude and phase both have significant geographical structure. This is most likely due to diurnally-varying flows back and forth across the coastline and elevation isolines resulting from the large daily warming and cooling over land. Because the region’s atmosphere is generally stably stratified, these flow patterns result in air of lower (higher) potential temperature being advected upslope (downslope) during daytime (nighttime). This suppresses the temperature diurnal cycle amplitude at mountaintops where diurnal flows converge (diverge) during the day (night). The nighttime land breeze also advects air of higher potential temperature downslope toward the coast. This raises minimum temperatures in land areas adjacent to the coast in a manner analogous to the daytime suppression of maximum temperature by the cool sea breeze in these same areas. Because stratification is greater in the coastal zone than in the desert interior, these thermal effects of the diurnal winds are not uniform, generating spatial structures in the phase and shape of the temperature diurnal cycle as well as its amplitude. We confirm that the simulated characteristics of the temperature diurnal cycle as well as those of the associated diurnal winds are also found in a network of 30 observation stations in the region. This gives confidence in the simulation’s realism and our study’s findings. Diurnal flows are probably mainly responsible for the geographical structures in the temperature diurnal cycle in other regions of significant topography and surface heterogeneity, their importance depending partly on the degree of atmospheric stratification.
Liou, KN, WL Lee, and A Hall. 2007. “Radiative transfer in mountains: Application to the Tibetan Plateau.” Geophysical Research Letters 34: L23809. Publisher's Version Abstract
We developed a 3D Monte Carlo photon tracing program for the transfer of radiation in inhomogeneous and irregular terrain to calculate broadband solar and thermal infrared fluxes. We selected an area of 100 × 100 km2 in the Tibetan Plateau centered at Lhasa city and used the albedo and surface temperature from MODIS/Terra for this study. We showed that anomalies of surface solar fluxes with reference to a flat surface can be as large as 600 W/m2, depending on time of day, mountain configuration, and albedo. Surface temperature is the dominating factor in determining anomalies of the surface infrared flux distribution relative to a flat surface with values as high as 70 W/m2 at cold mountain surfaces. The average surface solar flux over regional domains of 100 × 100 km2 and 50 × 50 km2 comprising intense topography can deviate from the smoothed surface conventionally assumed in climate models and GCMs by 10–50 W/m2.
2006
Differences in simulations of climate feedbacks are sources of significant divergence in climate models' temperature response to anthropogenic forcing. Snow albedo feedback is particularly critical for climate change prediction in heavily‐populated northern hemisphere land masses. Here we show its strength in current models exhibits a factor‐of‐three spread. These large intermodel variations in feedback strength in climate change are nearly perfectly correlated with comparably large intermodel variations in feedback strength in the context of the seasonal cycle. Moreover, the feedback strength in the real seasonal cycle can be measured and compared to simulated values. These mostly fall outside the range of the observed estimate, suggesting many models have an unrealistic snow albedo feedback in the seasonal cycle context. Because of the tight correlation between simulated feedback strength in the seasonal cycle and climate change, eliminating the model errors in the seasonal cycle will lead directly to a reduction in the spread of feedback strength in climate change. Though this comparison to observations may put the models in an unduly harsh light because of uncertainties in the observed estimate that are difficult to quantify, our results map out a clear strategy for targeted observation of the seasonal cycle to reduce divergence in simulations of climate sensitivity.
Chen, Y, A Hall, and KN Liou. 2006. “Application of three-dimensional solar radiative transfer to mountains.” Journal of Geophysical Research—Atmospheres 111: D21111. Publisher's Version Abstract
We developed a three‐dimensional radiative transfer model simulating solar fluxes over mountain surfaces precisely given distributions of atmospheric scatterers and absorbers. The model quantifies direct, diffuse, terrain‐reflected, and coupling (i.e., photons reflected and scattered more than once) fluxes. We applied it to a midlatitude mountainous surface to study these components' diurnal, seasonal, and geographical variability under clear skies. Domain‐averaged direct and diffuse fluxes together comprise over 96% of the flux year‐round, with diffuse fluxes' relative importance varying inversely with that of direct radiation. Direct fluxes generally account for at least 80% of the total. However, the domain‐averaged diffuse flux proportion increases to nearly 40% at high zenith angles, and approaches 100% when neighboring slopes obscure the surface from the Sun. Terrain‐reflected and coupling components each account for less than 1% throughout much of the year. However, together they comprise ∼3% when surface albedo increases during winter and are similarly nonnegligible in deep valleys all year. We also studied controls on geographical variations in flux components: The sky view factor, a conventional predictor of diffuse fluxes, is surprisingly weakly correlated with them, posing a parameterization challenge. Terrain‐reflected and coupling fluxes may be easier to parameterize given topography. Finally, we assessed shortwave errors in General Circulation Models with smoothed topography by comparing results with the mountainous surface to identical calculations for a flat surface with the same mean elevation. The differences range from 5 to 20 W/m2 and arise because the atmosphere absorbs a different amount of sunshine when underlying topography is smoothed.
Conil, S, and A Hall. 2006. “Local regimes of atmospheric variability: A case study of Southern California.” Journal of Climate 19: 4308–4325. Publisher's Version Abstract
The primary regimes of local atmospheric variability are examined in a 6-km regional atmospheric model of the southern third of California, an area of significant land surface heterogeneity, intense topography, and climate diversity. The model was forced by reanalysis boundary conditions over the period 1995–2003. The region is approximately the same size as a typical grid box of the current generation of general circulation models used for global climate prediction and reanalysis product generation, and so can be thought of as a laboratory for the study of climate at spatial scales smaller than those resolved by global simulations and reanalysis products. It is found that the simulated circulation during the October–March wet season, when variability is most significant, can be understood through an objective classification technique in terms of three wind regimes. The composite surface wind patterns associated with these regimes exhibit significant spatial structure within the model domain, consistent with the complex topography of the region. These regimes also correspond nearly perfectly with the simulation’s highly structured patterns of variability in hydrology and temperature, and therefore are the main contributors to the local climate variability. The regimes are approximately equally likely to occur regardless of the phase of the classical large-scale modes of atmospheric variability prevailing in the Pacific–North American sector. The high degree of spatial structure of the local regimes and their tightly associated climate impacts, as well as their ambiguous relationship with the primary modes of large-scale variability, demonstrate that the local perspective offered by the high-resolution model is necessary to understand and predict the climate variations of the region.
Bony, S, R Colman, V Kattsov, RP Allan, CS Bretherton, J-L Dufresne, A Hall, et al. 2006. “How well do we understand climate change feedback processes?” Journal of Climate 19: 3445–3482. Publisher's Version Abstract
Processes in the climate system that can either amplify or dampen the climate response to an external perturbation are referred to as climate feedbacks. Climate sensitivity estimates depend critically on radiative feedbacks associated with water vapor, lapse rate, clouds, snow, and sea ice, and global estimates of these feedbacks differ among general circulation models. By reviewing recent observational, numerical, and theoretical studies, this paper shows that there has been progress since the Third Assessment Report of the Intergovernmental Panel on Climate Change in (i) the understanding of the physical mechanisms involved in these feedbacks, (ii) the interpretation of intermodel differences in global estimates of these feedbacks, and (iii) the development of methodologies of evaluation of these feedbacks (or of some components) using observations. This suggests that continuing developments in climate feedback research will progressively help make it possible to constrain the GCMs’ range of climate feedbacks and climate sensitivity through an ensemble of diagnostics based on physical understanding and observations.
Qu, X, and A Hall. 2006. “Assessing snow albedo feedback in simulated climate change.” Journal of Climate 19: 2617–2630. Publisher's Version Abstract
In this paper, the two factors controlling Northern Hemisphere springtime snow albedo feedback in transient climate change are isolated and quantified based on scenario runs of 17 climate models used in the Intergovernmental Panel on Climate Change Fourth Assessment Report. The first factor is the dependence of planetary albedo on surface albedo, representing the atmosphere's attenuation effect on surface albedo anomalies. It is potentially a major source of divergence in simulations of snow albedo feedback because of large differences in simulated cloud fields in Northern Hemisphere land areas. To calculate the dependence, an analytical model governing planetary albedo was developed. Detailed validations of the analytical model for two of the simulations are shown, version 3 of the Community Climate System Model (CCSM3) and the Geophysical Fluid Dynamics Laboratory global coupled Climate Model 2.0 (CM2.0), demonstrating that it facilitates a highly accurate calculation of the dependence of planetary albedo on surface albedo given readily available simulation output. In all simulations it is found that surface albedo anomalies are attenuated by approximately half in Northern Hemisphere land areas as they are transformed into planetary albedo anomalies. The intermodel standard deviation in the dependence of planetary albedo on surface albedo is surprisingly small, less than 10% of the mean. Moreover, when an observational estimate of this factor is calculated by applying the same method to the satellite-based International Satellite Cloud Climatology Project (ISCCP) data, it is found that most simulations agree with ISCCP values to within about 10%, despite further disagreements between observed and simulated cloud fields. This suggests that even large relative errors in simulated cloud fields do not result in significant error in this factor, enhancing confidence in climate models. The second factor, related exclusively to surface processes, is the change in surface albedo associated with an anthropogenically induced temperature change in Northern Hemisphere land areas. It exhibits much more intermodel variability. The standard deviation is about ⅓ of the mean, with the largest value being approximately 3 times larger than the smallest. Therefore this factor is unquestionably the main source of the large divergence in simulations of snow albedo feedback. To reduce the divergence, attention should be focused on differing parameterizations of snow processes, rather than intermodel variations in the attenuation effect of the atmosphere on surface albedo anomalies.
2005
Hall, A, A Clement, DWJ Thompson, A Broccoli, and C Jackson. 2005. “The importance of atmospheric dynamics in the northern hemisphere wintertime climate response to changes in earth's orbit.” Journal of Climate 18: 1315–1325. Publisher's Version Abstract
Milankovitch proposed that variations in the earth’s orbit cause climate variability through a local thermodynamic response to changes in insolation. This hypothesis is tested by examining variability in an atmospheric general circulation model coupled to an ocean mixed layer model subjected to the orbital forcing of the past 165 000 yr. During Northern Hemisphere summer, the model’s response conforms to Milankovitch’s hypothesis, with high (low) insolation generating warm (cold) temperatures throughout the hemisphere. However, during Northern Hemisphere winter, the climate variations stemming from orbital forcing cannot be solely understood as a local thermodynamic response to radiation anomalies. Instead, orbital forcing perturbs the atmospheric circulation in a pattern bearing a striking resemblance to the northern annular mode, the primary mode of simulated and observed unforced atmospheric variability. The hypothesized reason for this similarity is that the circulation response to orbital forcing reflects the same dynamics generating unforced variability. These circulation anomalies are in turn responsible for significant fluctuations in other climate variables: Most of the simulated orbital signatures in wintertime surface air temperature over midlatitude continents are directly traceable not to local radiative forcing, but to orbital excitation of the northern annular mode. This has paleoclimate implications: during the point of the model integration corresponding to the last interglacial (Eemian) period, the orbital excitation of this mode generates a 1°–2°C warm surface air temperature anomaly over Europe, providing an explanation for the warm anomaly of comparable magnitude implied by the paleoclimate proxy record. The results imply that interpretations of the paleoclimate record must account for changes in surface temperature driven not only by changes in insolation, but also by perturbations in atmospheric dynamics.

Pages