Abstract

An attempt has been made to use paleoclimatic data from the last glacial maximum to evaluate the sensitivity of two versions of an atmosphere/mixed-layer ocean model. Each of these models has been used to study the C0₂-induced changes in climate. The models differ in their treatment of cloudiness, with one using a fixed cloud distribution and the other using a simple parameterization to predict clouds. The models also differ in the magnitude of their response to a doubling of atmospheric C0₂, with the variable cloud model being nearly twice as sensitive as the fixed cloud version. Given the distributions of continental ice sheets, surface albedo, and the reduced carbon dioxide concentration of the ice age, the climate of the last glacial maximum (LGM) is simulated by each model and compared with the corresponding simulation of the present climate. Both models generate differences in sea surface temperature and surface air temperature which compare favorably with estimates of the actual differences in temperature between the LGM and the present. However, it is difficult to determine which version of the model is more realistic in simulating the ice age climate for two reasons: 1) the differences between the two models are relatively small; and 2) there are substantial uncertainties in the pateoclimatic data. Neverthless, the similarity between the LGM simulations and the available paleoclimatic data suggests that the estimates of C0₂-induced climate change obtained from these models may not be too far from reality.

Abstract

The role of mountains in maintaining extensive midlatitude arid regions in the Northern Hemisphere was investigated using simulations from the GFDL Global Climate Model with and without orography. In the integration with mountains, dry climates were simulated over central Asia and the interior of North America, in good agreement with the observed climate. In contrast, moist climates were simulated in the same regions in the integration without mountains. During all season but summer, large amplitude stationary waves occur in response to the Tibetan Plateau and Rocky Mountains. The midlatitude dry regions are located upstream of the troughs of these waves, where general subsidence and relatively infrequent storm development occur and precipitation is thus inhibited. In summer, this mechanism contributes to the dryness of interior North America as a stationary wave trough remains east of the Rockies, but is not effective in Eurasia due to seasonal changes in the atmospheric circulation. The dryness of interior Eurasia in summer results, in part, from the south Asian monsoon circulation induced by the Tibetan Plateau. Its rising branch is centered above the southeastern Tibetan Plateau, and its salient features are a cyclonic flow at low levels (the “south Asian low”) and an anticyclonic flow in the upper troposphere. This circulation is associated with a northward displacement of the storm track and a flow of relatively dry, subsiding air across much of central Asia. In addition, land surface–atmosphere feedback contributes to the dryness of all midlatitude dry regions. Although the effect of this feedback is small in winter, it is responsible for more than half of the reduction in summer precipitation. Orography also substantially reduces the moisture transport across the continental interiors. The results from this experiment suggest that midlatitude dryness is largely due to the existence of orography. This is an alternative to the traditional explanation that distance from oceanic moisture sources, accentuated locally by the presence of mountain barriers upwind, is the major cause of midlatitude dry regions. Paleoclimatic evidence of less aridity during the late Tertiary, before substantial uplift of the Rocky Mountains and Tibetan Plateau is believed to have occurred, supports this possibility.

Abstract

High-impact extratropical cyclones (ETCs) cause considerable damage along the northeast coast of the United States through strong winds and inundation, but these relatively rare events are difficult to analyze owing to limited historical records. Using a 1505-yr simulation from the GFDL FLOR coupled model, statistical analyses of extreme events are performed including exceedance probability computations to compare estimates from shorter segments to estimates that could be obtained from a record of considerable length. The most extreme events possess characteristics including exceptionally low central pressure, hurricane-force winds, and a large surge potential, which would greatly impact nearby regions. Return level estimates of metrics of ETC intensity using shorter, historical-length segments of the FLOR simulation are underestimated compared to levels determined using the full simulation. This indicates that if the underlying distributions of observed ETC metrics are similar to those of the 1505-yr FLOR distributions, the actual frequency of extreme ETC events could also be underestimated. Comparisons between FLOR and reanalysis products suggest that not all features of simulated high-impact ETCs are representative of observations. Spatial track densities are similar, but FLOR exhibits a negative bias in central pressure and a positive bias in wind speed, particularly for more intense events. Although the existence of these model biases precludes the quantitative use of model-derived return statistics as a substitute for those derived from shorter observational records, this work suggests that statistics from future models of higher fidelity could be used to better constrain the probability of extreme ETC events and their impacts.

Abstract

The commonly held view of the conditions in the North Atlantic at the last glacial maximum, based on the interpretation of proxy records, is of large-scale cooling compared to today, limited deep convection, and extensive sea ice, all associated with a southward displaced and weakened overturning thermohaline circulation (THC) in the North Atlantic. Not all studies support that view; in particular, the “strength of the overturning circulation” is contentious and is a quantity that is difficult to determine even for the present day. Quasi-equilibrium simulations with coupled climate models forced by glacial boundary conditions have produced differing results, as have inferences made from proxy records. Most studies suggest the weaker circulation, some suggest little or no change, and a few suggest a stronger circulation.

Here results are presented from a three-dimensional climate model, the Hadley Centre Coupled Model version 3 (HadCM3), of the coupled atmosphere–ocean–sea ice system suggesting, in a qualitative sense, that these diverging views could all have occurred at different times during the last glacial period, with different modes existing at different times. One mode might have been characterized by an active THC associated with moderate temperatures in the North Atlantic and a modest expanse of sea ice. The other mode, perhaps forced by large inputs of meltwater from the continental ice sheets into the northern North Atlantic, might have been characterized by a sluggish THC associated with very cold conditions around the North Atlantic and a large areal cover of sea ice. The authors’ model simulation of such a mode, forced by a large input of freshwater, bears several of the characteristics of the Climate: Long-range Investigation, Mapping, and Prediction (CLIMAP) Project’s reconstruction of glacial sea surface temperature and sea ice extent.

Abstract

Climate models of varying complexity have been used for decades to investigate the impact of mountains on the atmosphere and surface climate. Here, the impact of removing the continental topography on the present-day ocean climate is investigated using three different climate models spanning multiple generations. An idealized study is performed where all present-day land surface topography is removed and the equilibrium change in the oceanic mean state with and without the mountains is studied. When the mountains are removed, changes found in all three models include a weakening of the Atlantic meridional overturning circulation and associated SST cooling in the subpolar North Atlantic. The SSTs also warm in all the models in the western North Pacific Ocean associated with a northward shift of the atmospheric jet and the Kuroshio. In the ocean interior, the magnitude of the temperature and salinity response to removing the mountains is relatively small and the sign and magnitude of the changes generally vary among the models. These different interior ocean responses are likely related to differences in the mean state of the control integrations due to differences in resolution and associated subgrid-scale mixing parameterizations. Compared to the results from 4xCO₂ simulations, the interior ocean temperature changes caused by mountain removal are relatively small; however, the oceanic circulation response and Northern Hemisphere near-surface temperature changes are of a similar magnitude to the response to such radiative forcing changes.

Abstract

Feedback analysis in climate models commonly involves decomposing any change in the system’s energy balance into radiative forcing terms due to prescribed changes, and response terms due to the radiative effects of changes in model variables such as temperature, water vapor, clouds, sea ice, and snow. The established “partial radiative perturbation” (PRP) method allows an accurate separation of these terms, but requires processing large volumes of model output with an offline version of the model’s radiation code. Here, we propose an “approximate PRP” (APRP) method for the shortwave that provides an accurate estimate of the radiative perturbation, but derived from a quite modest amount of monthly mean model output.

The APRP method is based on a simplified shortwave radiative model of the atmosphere, where surface absorption and atmospheric scattering and absorption are represented by means of three parameters that are diagnosed for overcast and clear-sky portions of each model grid cell. The accuracy of the method is gauged relative to full PRP calculations in two experiments: one in which carbon dioxide concentration is doubled and another in which conditions of the Last Glacial Maximum (LGM) are simulated. The approximate PRP method yields a shortwave cloud feedback accurate in the global mean to within 7%. Forcings and feedbacks due to surface albedo and noncloud atmospheric constituents are also well approximated with errors of order 5%–10%. Comparison of two different model simulations of the LGM shows that the regional and global differences in their ice sheet albedo forcing fields are clearly captured by the APRP method. Hence this method is an efficient and satisfactory tool for studying and intercomparing shortwave forcing and feedbacks in climate models.

Abstract

Eastern North America contains densely populated, highly developed areas, making winter storms with strong winds and high snowfall among the costliest storm types. For this reason, it is important to determine how the frequency of high-impact winter storms, specifically, those combining significant snowfall and winds, will change in this region under increasing greenhouse gas concentrations. This study uses a high-resolution coupled global climate model to simulate the changes in extreme winter conditions from the present climate to a future scenario with doubled CO₂ concentrations (2XC). In particular, this study focuses on changes in high-snowfall, extreme-wind (HSEW) events, which are defined as the occurrence of 2-day snowfall and high winds exceeding thresholds based on extreme values from the control simulation, where greenhouse gas concentrations remain fixed. Mean snowfall consistently decreases across the entire region, but extreme snowfall shows a more inconsistent pattern, with some areas experiencing increases in the frequency of extreme-snowfall events. Extreme-wind events show relatively small changes in frequency with 2XC, with the exception of high-elevation areas where there are large decreases in frequency. As a result of combined changes in wind and snowfall, HSEW events decrease in frequency in the 2XC simulation for much of eastern North America. Changes in the number of HSEW events in the 2XC environment are driven mainly by changes in the frequency of extreme-snowfall events, with most of the region experiencing decreases in event frequency, except for certain inland areas at higher latitudes.

Abstract

Using simulations performed with 24 coupled atmosphere–ocean global climate models from phase 5 of the Coupled Model Intercomparison Project (CMIP5), projections of Northern Hemisphere daily snowfall events under the RCP8.5 emissions scenario are analyzed for the periods of 2021–50 and 2071–2100 and compared to the historical period of 1971–2000. The overall frequency of daily snowfall events is simulated to decrease across much of the Northern Hemisphere, except at the highest latitudes such as northern Canada, northern Siberia, and Greenland. Seasonal redistributions of daily snowfall event frequency and average daily snowfall are also projected to occur in some regions. For example, large portions of the Northern Hemisphere, including much of Canada, Tibet, northern Scandinavia, northern Siberia, and Greenland, are projected to experience increases in average daily snowfall and event frequency in midwinter. But in warmer months, the regions with increased snowfall become fewer in number and are limited to northern Canada, northern Siberia, and Greenland. These simulations also show changes in the frequency distribution of daily snowfall event intensity, including an increase in heavier snowfall events even in some regions where the overall snowfall decreases. The projected changes in daily snowfall event frequency exhibit some dependence on the temperature biases of the individual models in certain regions and times of the year, with colder models typically toward the positive end of the distribution of event frequency changes and warmer models toward the negative end, particularly in regions near the transition zone between increasing and decreasing snowfall.

Abstract

The response of the equatorial Pacific Ocean’s seasonal cycle to orbital forcing is explored using idealized simulations with a coupled atmosphere–ocean GCM in which eccentricity, obliquity, and the longitude of perihelion are altered while other boundary conditions are maintained at preindustrial levels. The importance of ocean dynamics in the climate response is investigated using additional simulations with a slab ocean version of the model. Precession is found to substantially influence the equatorial Pacific seasonal cycle through both thermodynamic and dynamic mechanisms, while changes in obliquity have only a small effect. In the precession experiments, western equatorial Pacific SSTs respond in a direct thermodynamic manner to changes in insolation, while the eastern equatorial Pacific is first affected by the propagation of thermocline temperature anomalies from the west. These thermocline signals result from zonal wind anomalies associated with changes in the strength of subtropical anticyclones and shifts in the regions of convection in the western equatorial Pacific. The redistribution of heat from these thermocline signals, aided by the direct thermodynamic effect of insolation anomalies, results in large changes to the strength and timing of the eastern equatorial Pacific seasonal cycle. A comparison of 10 CMIP5 mid-Holocene experiments, in which the primary forcing is due to precession, shows that this response is relatively robust across models. Because equatorial Pacific SST anomalies have local climate impacts as well as nonlocal impacts through teleconnections, these results may be important to understanding paleoclimate variations both inside and outside of the tropical Pacific.

Abstract

The climate response to idealized changes in the atmospheric CO₂ concentration by the new GFDL climate model (CM2) is documented. This new model is very different from earlier GFDL models in its parameterizations of subgrid-scale physical processes, numerical algorithms, and resolution. The model was constructed to be useful for both seasonal-to-interannual predictions and climate change research. Unlike previous versions of the global coupled GFDL climate models, CM2 does not use flux adjustments to maintain a stable control climate. Results from two model versions, Climate Model versions 2.0 (CM2.0) and 2.1 (CM2.1), are presented.

Two atmosphere–mixed layer ocean or slab models, Slab Model versions 2.0 (SM2.0) and 2.1 (SM2.1), are constructed corresponding to CM2.0 and CM2.1. Using the SM2 models to estimate the climate sensitivity, it is found that the equilibrium globally averaged surface air temperature increases 2.9 (SM2.0) and 3.4 K (SM2.1) for a doubling of the atmospheric CO₂ concentration. When forced by a 1% per year CO₂ increase, the surface air temperature difference around the time of CO₂ doubling [transient climate response (TCR)] is about 1.6 K for both coupled model versions (CM2.0 and CM2.1). The simulated warming is near the median of the responses documented for the climate models used in the 2001 Intergovernmental Panel on Climate Change (IPCC) Working Group I Third Assessment Report (TAR).

The thermohaline circulation (THC) weakened in response to increasing atmospheric CO₂. By the time of CO₂ doubling, the weakening in CM2.1 is larger than that found in CM2.0: 7 and 4 Sv (1 Sv ≡ 10⁶ m³ s⁻¹), respectively. However, the THC in the control integration of CM2.1 is stronger than in CM2.0, so that the percentage change in the THC between the two versions is more similar. The average THC change for the models presented in the TAR is about 3 or 4 Sv; however, the range across the model results is very large, varying from a slight increase (+2 Sv) to a large decrease (−10 Sv).