NCAR Community Climate System Model Research Articles

Impact of ocean model resolution on CCSM climate simulations

The current literature provides compelling evidence suggesting that an eddy-resolving (as opposed to eddy-permitting or eddy-parameterized) ocean component model will significantly impact the simulation of the large-scale climate, although this has not been fully tested to date in multi-decadal global coupled climate simulations. The purpose of this paper is to examine how resolved ocean fronts and eddies impact the simulation of large-scale climate. The model used for this study is the NCAR Community Climate System Model version 3.5 (CCSM3.5)—the forerunner to CCSM4. Two experiments are reported here. The control experiment is a 155-year present-day climate simulation using a 0.5° atmosphere component (zonal resolution 0.625 meridional resolution 0.5°; land surface component at the same resolution) coupled to ocean and sea-ice components with zonal resolution of 1.2° and meridional resolution varying from 0.27° at the equator to 0.54° in the mid-latitudes. The second simulation uses the same atmospheric and land-surface models coupled to eddy-resolving 0.1° ocean and sea-ice component models. The simulations are compared in terms of how the representation of smaller scale features in the time mean ocean circulation and ocean eddies impact the mean and variable climate. In terms of the global mean surface temperature, the enhanced ocean resolution leads to a ubiquitous surface warming with a global mean surface temperature increase of about 0.2 °C relative to the control. The warming is largest in the Arctic and regions of strong ocean fronts and ocean eddy activity (i.e., Southern Ocean, western boundary currents). The Arctic warming is associated with significant losses of sea-ice in the high-resolution simulation. The sea surface temperature gradients in the North Atlantic, in particular, are better resolved in the high-resolution model leading to significantly sharper temperature gradients and associated large-scale shifts in the rainfall. In the extra-tropics, the interannual temperature variability is increased with the resolved eddies, and a notable increases in the amplitude of the El Nino and the Southern Oscillation is also detected. Changes in global temperature anomaly teleconnections and local air-sea feedbacks are also documented and show large changes in ocean–atmosphere coupling. In particular, local air-sea feedbacks are significantly modified by the increased ocean resolution. In the high-resolution simulation in the extra-tropics there is compelling evidence of stronger forcing of the atmosphere by SST variability arising from ocean dynamics. This coupling is very weak or absent in the low-resolution model.

Extratropical modulation on Asian summer monsoon at precessional bands

In a 46‐kyrs' transient accelerated simulation with the NCAR Community Climate System Model version 3 (CCSM3), the increased July insolation at precessional bands causes a drastic summer warming in North Africa, the mid‐latitudes of Eurasia from surface to the upper troposphere, with a cooling troposphere in the Northwest Pacific (NWP), which intensifies the land‐ocean thermal contrast between Eurasia and its adjacent oceans. Associated with the tropospheric temperature variations are a lower‐tropospheric strengthened low‐pressure system from North Africa to the mid‐high latitudes of Eurasia and a lower‐tropospheric strengthened high‐pressure system over the NWP. Meanwhile, the strengthened southwesterly or southerly winds prevail over the Asian monsoon region, indicating a stronger Asian summer monsoon. Summer precipitation increases in tropical Africa, South Asia, and the mid‐high latitudes of East Asia, with a decreased precipitation over South China and the NWP. This response of atmospheric circulation and monsoon precipitation to the precessional changes is explained in the context of the extratropical Asian‐Pacific Oscillation.

Evaluation of WRF Mean and Extreme Precipitation over Spain: Present Climate (1970–99)

Abstract The ability of the Weather Research and Forecasting model (WRF) to simulate precipitation over Spain is evaluated from a climatological point of view. The complex topography and the large rainfall variability make the Iberian Peninsula a particularly interesting region and permit assessment of model performance under very demanding conditions. Three high-resolution (10 km) simulations over the Iberian Peninsula have been completed spanning a 30-yr period (1970–99) and driven by different datasets: the 40-yr European Centre for Medium-Range Weather Forecasts Re-Analysis (ERA-40) as “perfect boundary conditions” and two general circulation models (GCMs), the Max Planck Institute ECHAM5 model (ECHAM5/MPI) and the NCAR Community Climate System Model, version 3 (CCSM3). The daily precipitation observational grid Spain02 is employed to evaluate the model at varying time scales. Not only are the long-term means (annual, seasonal, and monthly) examined but also the high-order statistics (extreme events). The WRF provides valuable information on precipitation at high resolution and enhances local spatial distribution due to orographic features. Although substantial errors are still observed in terms of monthly precipitation, especially during the spring, the model is largely able to capture the various precipitation regimes. The major benefits of using WRF are related to the spatial distribution of rainfall and the simulation of extreme events, two facets of climate that can be barely explored with GCMs. This study shows that WRF can be a useful tool for generating high-resolution climate information for Spanish precipitation at spatial and temporal scales that are crucial for both the environment and human life.

Exploring Mechanisms of Variability and Predictability of Atlantic Meridional Overturning Circulation in Two Coupled Climate Models

Abstract Multidecadal variability in the Atlantic meridional overturning circulation (AMOC) of the ocean is diagnosed in the NCAR Community Climate System Model, version 3 (CCSM3), and the GFDL Coupled Model (CM2.1). Common diagnostic approaches are applied to draw out similarities and differences between the two models. An index of AMOC variability is defined, and the manner in which key variables covary with it is determined. In both models the following is found. (i) AMOC variability is associated with upper-ocean (top 1 km) density anomalies (dominated by temperature) on the western margin of the basin in the region of the Mann eddy with a period of about 20 years. These anomalies modulate the trajectory and strength of the North Atlantic Current. The importance of the western margin is a direct consequence of the thermal wind relation and is independent of the mechanisms that create those density anomalies. (ii) Density anomalies in this key region are part of a larger-scale pattern that propagates around the subpolar gyre and acts as a “pacemaker” of AMOC variability. (iii) The observed variability is consistent with the primary driving mechanism being stochastic wind curl forcing, with Labrador Sea convection playing a secondary role. Also, “toy models” of delayed oscillator form are fitted to power spectra of key variables and are used to infer “quality factors” (Q-factors), which characterize the bandwidth relative to the center frequency and hence AMOC predictability horizons. The two models studied here have Q-factors of around 2, suggesting that prediction is possible out to about two cycles, which is likely larger than the real AMOC.

ENSO and Pacific Decadal Variability in the Community Climate System Model Version 4

Abstract This study presents an overview of the El Niño–Southern Oscillation (ENSO) phenomenon and Pacific decadal variability (PDV) simulated in a multicentury preindustrial control integration of the NCAR Community Climate System Model version 4 (CCSM4) at nominal 1° latitude–longitude resolution. Several aspects of ENSO are improved in CCSM4 compared to its predecessor CCSM3, including the lengthened period (3–6 yr), the larger range of amplitude and frequency of events, and the longer duration of La Niña compared to El Niño. However, the overall magnitude of ENSO in CCSM4 is overestimated by ~30%. The simulated ENSO exhibits characteristics consistent with the delayed/recharge oscillator paradigm, including correspondence between the lengthened period and increased latitudinal width of the anomalous equatorial zonal wind stress. Global seasonal atmospheric teleconnections with accompanying impacts on precipitation and temperature are generally well simulated, although the wintertime deepening of the Aleutian low erroneously persists into spring. The vertical structure of the upper-ocean temperature response to ENSO in the north and south Pacific displays a realistic seasonal evolution, with notable asymmetries between warm and cold events. The model shows evidence of atmospheric circulation precursors over the North Pacific associated with the “seasonal footprinting mechanism,” similar to observations. Simulated PDV exhibits a significant spectral peak around 15 yr, with generally realistic spatial pattern and magnitude. However, PDV linkages between the tropics and extratropics are weaker than observed.

Improved Madden–Julian Oscillations with Improved Physics: The Impact of Modified Convection Parameterizations

Abstract Two modifications are made to the deep convection parameterization in the NCAR Community Climate System Model, version 3 (CCSM3): a dilute plume approximation and an implementation of the convective momentum transport (CMT). These changes lead to significant improvement in the simulated Madden–Julian oscillations (MJOs). With the dilute plume approximation, temperature and convective heating perturbations become more positively correlated. Consequently, more available potential energy is generated and the intraseasonal variability (ISV) becomes stronger. The organization of ISV is also improved, which is manifest in coherent structures between different MJO phases and an improved simulation of the eastward propagation of MJOs with a reasonable eastward speed. The improved propagation can be attributed to a better simulation of the low-level zonal winds due to the inclusion of CMT. The authors posit that the large-scale zonal winds are akin to a selective conveyor belt that facilitates the organization of ISVs into highly coherent structures, which are important features of observed MJOs. The conclusions are supported by two supplementary experiments, which include the dilute plume approximation and CMT separately.

Probabilistic Projections of Climate Change over China under the SRES A1B Scenario Using 28 AOGCMs

Probabilistic projection of climate change consists of formulating the climate change information in a probabilistic manner at either global or regional scale. This can produce useful results for studies of the impact of climate change impact and change mitigation. In the present study, a simple yet effective approach is proposed with the purpose of producing probabilistic results of climate change over China for the middle and end of the twenty-first century under the Special Report on Emissions Scenarios A1B (SRES A1B) emission scenario. Data from 28 coupled atmosphere–ocean general circulation models (AOGCMs) are used. The methodology consists of ranking the 28 models, based on their ability to simulate climate over China in terms of two model evaluation metrics. Different weights were then given to the models according to their performances in present-day climate. Results of the evaluation for the current climate show that five models that have relatively higher resolutions—namely, the Istituto Nazionale di Geofisica e Vulcanologia ECHAM4 (INGV ECHAM4), the third climate configuration of the Met Office Unified Model (UKMO HadCM3), the CSIRO Mark version 3.5 (Mk3.5), the NCAR Community Climate System Model, version 3 (CCSM3), and the Model for Interdisciplinary Research on Climate 3.2, high-resolution version [MIROC3.2 (hires)]—perform better than others over China. Their corresponding weights (normalized to 1) are 0.289, 0.096, 0.058, 0.048, and 0.044, respectively. Under the A1B scenario, surface air temperature is projected to increase significantly for both the middle and end of the twenty-first century, with larger magnitude over the north and in winter. There are also significant increases in rainfall in the twenty-first century under the A1B scenario, especially for the period 2070–99. As far as the interannual variability is concerned, the most striking feature is that there are high probabilities for the future intensification of interannual variability of precipitation over most of China in both winter and summer. For instance, over the Yangtze–Huai River basin (28°–35°N, 105°–120°E), there is a 60% probability of increased interannual standard deviation of precipitation by 20% in summer, which is much higher than that of the mean precipitation. In general there are small differences between weighted and unweighted projections, but the uncertainties in the projected changes are reduced to some extent after weighting.

Evaluation of the HOMME Dynamical Core in the Aquaplanet Configuration of NCAR CAM4: Rainfall

Abstract The NCAR Community Climate System Model, version 4 (CCSM4), includes a new dynamical core option based on NCAR’s High-Order Method Modeling Environment (HOMME). HOMME is a petascale-capable high-order element-based conservative dynamical core developed on the cubed-sphere grid. Initial simulations have been completed in an aquaplanet configuration of the Community Atmosphere Model, version 4 (CAM4), the atmospheric component of CCSM4. The authors examined the results of this simulation and assessed its fidelity in simulating rainfall, which is one of the most important components of the earth’s climate system. For this they compared the results from two other dynamical cores of CAM4: the finite volume (FV) and Eulerian (EUL). Instantaneous features of rainfall in HOMME are similar to FV and EUL. Similar to EUL and FV, HOMME simulates a single-peak intertropical convergence zone (ITCZ) over the equator. The strength of the ITCZ is found to be almost the same in HOMME and EUL but more in FV. It is observed that in HOMME and EUL, there is higher surface evaporation, which supplies more moisture to the deep tropics and gives more rainfall over the ITCZ. The altitude of maximum precipitation is found to be at almost the same level in all three dynamical cores. The eastward propagation of rainfall bands is organized and more prominent in FV and HOMME than in EUL. The phase speed of the eastward propagation in HOMME is found to be higher than in FV. The results show that, in general, the rainfall simulated by HOMME falls in a regime between that of FV and EUL. Hence, they conclude that the key aspects of rainfall simulation with HOMME falls into an acceptable range, as compared to the existing dynamical cores used in the model.

Atmospheric circulation anomalies due to 115 kyr BP climate forcing dominated by changes in the North Pacific Ocean

Climate at the time of inception of the Laurentide Ice Sheet (LIS) at ~115 kyr BP is simulated with the fully coupled NCAR Community Climate System Model (CCSM3) and compared to a simulated preindustrial climate (circa 1870) in order to better understand land surface and atmospheric responses to orbital and greenhouse cooling at inception. The interaction between obliquity and eccentricity produces maximum decrease in TOA insolation in JJA over the Arctic but increases occur over the tropics in DJF. The land surface response is dominated by widespread summer cooling in the Northern Hemisphere (NH), increases in snowfall, and decreases in melt rates and total precipitation. CCSM3 responds to the climate forcing at 115 kyr BP by producing incipient glaciation in the areas of LIS nucleation. We find that the inception of the LIS could have occurred with atmospheric circulation patterns that differ little from the present. The location of the troughs/ridges, mean flow over the Canadian Arctic and dominant modes of the atmospheric circulation are all very similar to the present. Larger changes in mean sea level pressure occur upstream of the inception region in the North Pacific Ocean and downstream in Western Europe. In the North Pacific region, the 115 kyr BP anomalies weaken both the Pacific high and Aleutian low making NH summers look more like the PREIND winters and vice versa. The occurrence of cold JJA anomalies at 115 kyr BP favors outbreaks of cold air not in the winter as in contemporary climates but during the summer instead and reinforces the cooling from orbital and GHG reductions. Increased poleward eddy transport of heat and moisture characterizes the atmospheric response in addition to reduced total cloud cover in the Arctic.

Seasonal Variations of Climate Feedbacks in the NCAR CCSM3

Abstract This study investigates the annual cycle of radiative contributions to global climate feedbacks. A partial radiative perturbation (PRP) technique is used to diagnose monthly radiative perturbations at the top of atmosphere (TOA) due to CO2 forcing; surface temperature response; and water vapor, cloud, lapse rate, and surface albedo feedbacks using NCAR Community Climate System Model, version 3 (CCSM3) output from a Special Report on Emissions Scenarios (SRES) A1B emissions-scenario-forced climate simulation. The seasonal global mean longwave TOA radiative feedback was found to be minimal. However, the global mean shortwave (SW) TOA cloud and surface albedo radiative perturbations exhibit large seasonality. The largest contributions to the negative SW cloud feedback occur during summer in each hemisphere, marking the largest differences with previous results. Results suggest that intermodel spread in climate sensitivity may occur, partially from cloud and surface albedo feedback seasonality differences. Further, links between the climate feedback and surface temperature response seasonality are investigated, showing a strong relationship between the seasonal climate feedback distribution and the seasonal surface temperature response.

The Greenland Ice Sheet Response to Transient Climate Change

Abstract This study applies a multiphase, multiple-rheology, scalable, and extensible geofluid model to the Greenland Ice Sheet (GrIS). The model is driven by monthly atmospheric forcing from global climate model simulations. Novel features of the model, referred to as the scalable and extensible geofluid modeling system (SEGMENT-Ice), include using the full Navier–Stokes equations to account for nonlocal dynamic balance and its influence on ice flow, and a granular sliding layer between the bottom ice layer and the lithosphere layer to provide a mechanism for possible large-scale surges in a warmer future climate (granular basal layer is for certain specific regions, though). Monthly climate of SEGMENT-Ice allows an investigation of detailed features such as seasonal melt area extent (SME) over Greenland. The model reproduced reasonably well the annual maximum SME and total ice mass lost rate when compared observations from the Special Sensing Microwave Imager (SSM/I) and Gravity Recovery and Climate Experiment (GRACE) over the past few decades. The SEGMENT-Ice simulations are driven by projections from two relatively high-resolution climate models, the NCAR Community Climate System Model, version 3 (CCSM3) and the Model for Interdisciplinary Research on Climate 3.2, high-resolution version [MIROC3.2(hires)], under a realistic twenty-first-century greenhouse gas emission scenario. They suggest that the surface flow would be enhanced over the entire GrIS owing to a reduction of ice viscosity as the temperature increases, despite the small change in the ice surface topography over the interior of Greenland. With increased surface flow speed, strain heating induces more rapid heating in the ice at levels deeper than due to diffusion alone. Basal sliding, especially for granular sediments, provides an efficient mechanism for fast-glacier acceleration and enhanced mass loss. This mechanism, absent from other models, provides a rapid dynamic response to climate change. Net mass loss estimates from the new model should reach ~220 km3 yr−1 by 2100, significantly higher than estimates by the Intergovernmental Panel on Climate Change (IPCC) Assessment Report 4 (AR4) of ~50–100 km3 yr−1. By 2100, the perennial frozen surface area decreases up to ~60%, to ~7 × 105 km2, indicating a massive expansion of the ablation zone. Ice mass change patterns, particularly along the periphery, are very similar between the two climate models.

Geographical Distribution of Climate Feedbacks in the NCAR CCSM3.0

Abstract This study performs offline, partial radiative perturbation calculations to determine the geographical distributions of climate feedbacks contributing to the top-of-atmosphere (TOA) radiative energy budget. These radiative perturbations are diagnosed using monthly mean model output from the NCAR Community Climate System Model version 3 (CCSM3.0) forced with the Special Report Emissions Scenario (SRES) A1B emission scenario. The Monte Carlo Independent Column Approximation (MCICA) technique with a maximum–random overlap rule is used to sample monthly mean cloud frequency profiles to perform the radiative transfer calculations. It is shown that the MCICA technique provides a good estimate of all feedback sensitivity parameters. The radiative perturbation results are used to investigate the spatial variability of model feedbacks showing that the shortwave cloud and lapse rate feedbacks exhibit the most and second most spatial variability, respectively. It has been shown that the model surface temperature response is highly correlated with the change in the TOA net flux, and that the latter is largely determined by the total feedback spatial pattern rather than the external forcing. It is shown by representing the change in the TOA net flux as a linear combination of individual feedback radiative perturbations that the lapse rate explains the most spatial variance of the surface temperature response. Feedback spatial patterns are correlated with the model response and other feedback spatial patterns to investigate these relationships. The results indicate that the model convective response is strongly correlated with cloud and water vapor feedbacks, but the lapse rate feedback geographic distribution is strongly correlated with the climatological distribution of convection. The implication for the water vapor–lapse rate anticorrelation is discussed.

CLARREO shortwave observing system simulation experiments of the twenty-first century: Simulator design and implementation

[1] Projected changes in the Earth system will likely be manifested in changes in reflected solar radiation. This paper introduces an operational Observational System Simulation Experiment (OSSE) to calculate the signals of future climate forcings and feedbacks in top-of-atmosphere reflectance spectra. The OSSE combines simulations from the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report for the NCAR Community Climate System Model (CCSM) with the MODTRAN radiative transfer code to calculate reflectance spectra for simulations of current and future climatic conditions over the 21st century. The OSSE produces narrowband reflectances and broadband fluxes, the latter of which have been extensively validated against archived CCSM results. The shortwave reflectance spectra contain atmospheric features including signals from water vapor, liquid and ice clouds, and aerosols. The spectra are also strongly influenced by the surface bidirectional reflectance properties of predicted snow and sea ice and the climatological seasonal cycles of vegetation. By comparing and contrasting simulated reflectance spectra based on emissions scenarios with increasing projected and fixed present-day greenhouse gas and aerosol concentrations, we find that prescribed forcings from increases in anthropogenic sulfate and carbonaceous aerosols are detectable and are spatially confined to lower latitudes. Also, changes in the intertropical convergence zone and poleward shifts in the subsidence zones and the storm tracks are all detectable along with large changes in snow cover and sea ice fraction. These findings suggest that the proposed NASA Climate Absolute Radiance and Refractivity Observatory (CLARREO) mission to measure shortwave reflectance spectra may help elucidate climate forcings, responses, and feedbacks.

Erroneous Model Field Representations in Multiple Pseudoproxy Studies: Corrections and Implications*

Abstract Pseudoproxy experiments evaluate statistical methods used to reconstruct climate fields from paleoclimatic proxies during the Common Era. These experiments typically employ output from millennial simulations by general circulation models (GCMs). It is demonstrated that multiple published pseudoproxy studies have used erroneously processed GCM surface temperature fields: the NCAR Community Climate System Model (CCSM), version 1.4, field was incorrectly oriented geographically and the GKSS ECHO-g FOR1 field was corrupted by a hemispheric-scale smoothing in the Western Hemisphere. These problems are not associated with the original model simulations; they instead arose because of incorrect processing of the model data for the pseudoproxy experiments. The consequences of these problems are evaluated for the studies in which the incorrect fields were used. Some quantitative results are invalidated by the findings: these include all experiments that used the corrupted ECHO-g field and those aspects of previous CCSM experiments that focused on Niño-3 reconstructions. Other results derived from the CCSM field can be reinterpreted based on the information provided herein and their qualitative characteristics remain similar.

Convection Parameterization, Tropical Pacific Double ITCZ, and Upper-Ocean Biases in the NCAR CCSM3. Part II: Coupled Feedback and the Role of Ocean Heat Transport

Abstract This study investigates the coupled atmosphere–ocean feedback and the role of ocean dynamic heat transport in the formation of double ITCZ over the tropical Pacific in the NCAR Community Climate System Model, version 3 (CCSM3) and its alleviation when a revised Zhang–McFarlane (ZM) convection scheme is used. A hierarchy of coupling strategy is employed for this purpose. A slab ocean model is coupled with the atmospheric component of the Community Atmosphere Model, version 3 (CAM3) to investigate the local feedback between the atmosphere and the ocean. It is shown that the net surface energy flux differences in the southern ITCZ region between the revised and original ZM scheme seen in the stand-alone CAM3 simulations can cool the SST by up to 1.5°C. However, the simulated SST distribution is very sensitive to the prescribed ocean heat transport required in the slab ocean model. To understand the role of ocean heat transport, the fully coupled CCSM3 model is used. The analysis of CCSM3 simulations shows that the altered ocean dynamic heat transport when the revised ZM scheme is used is largely responsible for the reduction of SST bias in the southern ITCZ region, although surface energy flux also helps to cool the SST in the first few months of the year in seasonal variation. The results, together with those from Part I, suggest that the unrealistic simulation of convection over the southern ITCZ region in the standard CCSM3 leads to the double-ITCZ bias through complex coupled interactions between atmospheric convection, surface winds, latent heat flux, cloud radiative forcing, SST, and upper-ocean circulations. The mitigation of the double-ITCZ bias using the revised ZM scheme is achieved by altering this chain of interactions.

Influence of ENSO and the Atlantic Multidecadal Oscillation on Drought over the United States

Abstract Composites based on observations and model outputs from the Climate Variability and Predictability (CLIVAR) drought experiments were used to examine the impact of El Niño–Southern Oscillation (ENSO) and the Atlantic multidecadal oscillation (AMO) on drought over the United States. Because drought implies persistent dryness, the 6-month standardized precipitation index, standardized runoff index, and soil moisture anomalies are used to represent drought. The experiments were performed by forcing an AGCM with prescribed sea surface temperature anomalies (SSTAs) superimposed on the monthly mean SST climatology. Four model outputs from the NCEP Global Forecast System (GFS), NASA’s Seasonal-to-Interannual Prediction Project, version 1 (NSIPP1), GFDL’s global atmospheric model, version 2.1 (AM2.1), and the Lamont-Doherty Earth Observatory (LDEO)/NCAR Community Climate System Model, version 3 (CCM3) were analyzed in this study. Each run lasts from 36 to 51 yr. The impact of ENSO on drought over the United States is concentrated over the Southwest, the Great Plains, and the lower Colorado River basin, with cold (warm) ENSO events favoring drought (wet spells). Over the East Coast and the Southeast, the impact of ENSO is small because the precipitation responses to ENSO are opposite in sign for winter and summer. For these areas, a prolonged ENSO does not always favor either drought or wet spells. The direct influence of the AMO on drought is small. The major influence of the AMO is to modulate the impact of ENSO on drought. The influence is large when the SSTAs in the tropical Pacific and in the North Atlantic are opposite in phase. A cold (warm) event in a positive (negative) AMO phase amplifies the impact of the cold (warm) ENSO on drought. The ENSO influence on drought is much weaker when the SSTAs in the tropical Pacific and in the North Atlantic are in phase.

Convection Parameterization, Tropical Pacific Double ITCZ, and Upper-Ocean Biases in the NCAR CCSM3. Part I: Climatology and Atmospheric Feedback

Abstract The role of convection parameterization in the formation of double ITCZ and associated upper-ocean biases in the NCAR Community Climate System Model, version 3 (CCSM3) is investigated by comparing the simulations using the original and revised Zhang–McFarlane (ZM) convection schemes. Ten-year model climatologies show that the simulation with the original ZM scheme produces a typical double ITCZ bias, whereas all biases related to the spurious double ITCZ and overly strong cold tongue in precipitation, sea surface temperature (SST), wind stress, ocean thermocline, upper-ocean currents, temperature, and salinity are dramatically reduced when the revised ZM scheme is used. These results demonstrate that convection parameterization plays a critical role in the formation of double ITCZ bias in the CCSM3. To understand the physical mechanisms through which the modifications of the convection scheme in the atmospheric model alleviate the double ITCZ bias in the CCSM3, the authors investigate the impacts of convection schemes on the atmospheric forcing and feedback in the uncoupled Community Atmospheric Model, version 3 (CAM3). It is shown that the CAM3 simulation with the original ZM scheme also produces a signature of double ITCZ bias in precipitation, whereas the simulation with the revised ZM scheme does not. Diagnostic analyses have identified three factors on the atmospheric side (i.e., the sensitivity of convection to SST, the convection–shortwave flux–SST feedback, and the convection–wind–evaporation–SST feedback) that may contribute to the differences in the coupled simulations.

The sensitivity of the climate response to the magnitude and location of freshwater forcing: last glacial maximum experiments

Proxy records indicate that the locations and magnitudes of freshwater forcing to the Atlantic Ocean basin as iceberg discharges into the high-latitude North Atlantic, Laurentide meltwater input to the Gulf of Mexico, or meltwater diversion to the North Atlantic via the St. Lawrence River and other eastern outlets may have influenced the North Atlantic thermohaline circulation and global climate. We have performed Last Glacial Maximum (LGM) simulations with the NCAR Community Climate System Model (CCSM3) in which the magnitude of the freshwater forcing has been varied from 0.1 to 1 Sv and inserted either into the subpolar North Atlantic Ocean or the Gulf of Mexico. In these glacial freshening experiments, the less dense freshwater provides a lid on the ocean water below, suppressing ocean convection and interaction with the atmosphere above and reducing the Atlantic Meridional Overturning Circulation (AMOC). This is the case whether the freshwater is added directly to the area of convection south of Greenland or transported there by the subtropical and subpolar gyres when added to the Gulf of Mexico. The AMOC reduction is less for the smaller freshwater forcings, but is not linear with the size of the freshwater perturbation. The recovery of the AMOC from a “slow” state is ∼200 years for the 0.1 Sv experiment and ∼500 years for the 1 Sv experiment. For glacial climates, with large Northern Hemisphere ice sheets and reduced greenhouse gases, the cold subpolar North Atlantic is primed to respond rapidly and dramatically to freshwater that is either directly dumped into this region or after being advected from the Gulf of Mexico. Greenland temperatures cool by 6–8 °C in all the experiments, with little sensitivity to the magnitude, location or duration of the freshwater forcing, but exhibiting large seasonality. Sea ice is important for explaining the responses. The Northern Hemisphere high latitudes are slow to recover. Antarctica and the Southern Ocean show a bipolar response, with warming and reduced sea ice. This warming continues after the cessation of the freshwater forcing and shows a dependence on the duration of the freshwater forcing. Equatorward of the expanded sea ice, the simulated temperature and salinity anomalies are sensitive to the amount of colder and fresher waters that are advected out of the subpolar North Atlantic. In the tropical Atlantic, the recovery of the Intertropical Convergence Zone (ITCZ) from its more southerly position during the freshwater forcing is much more rapid than the recovery of the AMOC, and is more related to the recovery of low-latitude surface temperatures than Greenland temperature or sea ice. These results have implications for using proxy records as indirect measures of the AMOC.

Climate response to large, high‐latitude and low‐latitude volcanic eruptions in the Community Climate System Model

Explosive volcanism is known to be a leading natural cause of climate change. The second half of the 13th century was likely the most volcanically perturbed half‐century of the last 2000 years, although none of the major 13th century eruptions have been clearly attributed to specific volcanoes. This period was in general a time of transition from the relatively warm Medieval period to the colder Little Ice Age, but available proxy records are insufficient on their own to clearly assess whether this transition is associated with volcanism. This context motivates our investigation of the climate system sensitivity to high‐ and low‐latitude volcanism using the fully coupled NCAR Community Climate System Model (CCSM3). We evaluate two sets of ensemble simulations, each containing four volcanic pulses, with the first set representing them as a sequence of tropical eruptions and the second representing eruptions occurring in the mid‐high latitudes of both the Northern and Southern hemispheres. The short‐term, direct radiative impacts of tropical and high‐latitude eruptions include significant cooling over the continents in summer and cooling over regions of increased sea‐ice concentration in Northern Hemisphere (NH) winter. A main dynamical impact of moderate tropical eruptions is a winter warming pattern across northern Eurasia. Furthermore, both ensembles show significant reductions in global precipitation, especially in the summer monsoon regions. The most important long‐term impact is the cooling of the high‐latitude NH produced by multiple tropical eruptions, suggesting that positive feedbacks associated with ice and snow cover could lead to long‐term climate cooling in the Arctic.

Early Eocene Arctic climate sensitivity to pCO2 and basin geography

We present results from new early Eocene (∼55–45 Ma) climate modeling experiments with the NCAR Community Climate System Model. These experiments test the sensitivity of climate to a large increase in atmospheric greenhouse gases, such as may have occurred at the Paleocene‐Eocene (P‐E) boundary (∼55.5 Ma), and also allow us to explore the role of Arctic basin configuration on climate. Experiments were run with pCO2 at 560 and 2240 ppm, and a third experiment, at 2240 ppm, incorporates a passage to a neighboring ocean to explore the potential effect of the ocean on Arctic warming, were the Arctic not isolated. Quadrupling pCO2 warms the Arctic by ∼8°C in the annual average, doubles atmospheric moisture content in this region and eliminates Arctic sea ice, consistent with proxy estimates of warming at the P‐E boundary. Opening the Arctic Ocean warms mean annual sea surface temperature by an additional ∼4°C.