NSF NCAR Climate and Global Dynamics Laboratory

Abstract An overview of the Community Earth System Model Version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long preindustrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low‐top” (40 km, with limited chemistry) and “high‐top” (130 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low‐latitude precipitation and shortwave cloud forcing biases; better representation of the Madden‐Julian Oscillation; better El Niño‐Southern Oscillation‐related teleconnections; and a global land carbon accumulation trend that agrees well with observationally based estimates. Most tropospheric and surface features of the low‐ and high‐top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1–5.3 °C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2's transient climate response of 1.9–2.0 °C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and they represent CESM2's contributions to the Coupled Model Intercomparison Project Phase 6.

Abstract While internal climate variability is known to affect climate projections, its influence is often underappreciated and confused with model error. Why? In general, modeling centers contribute a small number of realizations to international climate model assessments [e.g., phase 5 of the Coupled Model Intercomparison Project (CMIP5)]. As a result, model error and internal climate variability are difficult, and at times impossible, to disentangle. In response, the Community Earth System Model (CESM) community designed the CESM Large Ensemble (CESM-LE) with the explicit goal of enabling assessment of climate change in the presence of internal climate variability. All CESM-LE simulations use a single CMIP5 model (CESM with the Community Atmosphere Model, version 5). The core simulations replay the twenty to twenty-first century (1920–2100) 30 times under historical and representative concentration pathway 8.5 external forcing with small initial condition differences. Two companion 1000+-yr-long preindustrial control simulations (fully coupled, prognostic atmosphere and land only) allow assessment of internal climate variability in the absence of climate change. Comprehensive outputs, including many daily fields, are available as single-variable time series on the Earth System Grid for anyone to use. Early results demonstrate the substantial influence of internal climate variability on twentieth- to twenty-first-century climate trajectories. Global warming hiatus decades occur, similar to those recently observed. Internal climate variability alone can produce projection spread comparable to that in CMIP5. Scientists and stakeholders can use CESM-LE outputs to help interpret the observational record, to understand projection spread and to plan for a range of possible futures influenced by both internal climate variability and forced climate change.

Global temperatures are known to have varied over the past 1500 years, but the spatial patterns have remained poorly defined. We used a global climate proxy network to reconstruct surface temperature patterns over this interval. The Medieval period is found to display warmth that matches or exceeds that of the past decade in some regions, but which falls well below recent levels globally. This period is marked by a tendency for La Niña-like conditions in the tropical Pacific. The coldest temperatures of the Little Ice Age are observed over the interval 1400 to 1700 C.E., with greatest cooling over the extratropical Northern Hemisphere continents. The patterns of temperature change imply dynamical responses of climate to natural radiative forcing changes involving El Niño and the North Atlantic Oscillation-Arctic Oscillation.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate – the “global carbon budget” – is important tobetter understand the global carbon cycle, support the development ofclimate policies, and project future climate change. Here we describe andsynthesize data sets and methodology to quantify the five major componentsof the global carbon budget and their uncertainties. Fossil CO2emissions (EFOS) are based on energy statistics and cement productiondata, while emissions from land-use change (ELUC), mainlydeforestation, are based on land use and land-use change data andbookkeeping models. Atmospheric CO2 concentration is measured directlyand its growth rate (GATM) is computed from the annual changes inconcentration. The ocean CO2 sink (SOCEAN) and terrestrialCO2 sink (SLAND) are estimated with global process modelsconstrained by observations. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the lastdecade available (2010–2019), EFOS was 9.6 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.4 ± 0.5 GtC yr−1 when the cement carbonation sink is included), andELUC was 1.6 ± 0.7 GtC yr−1. For the same decade, GATM was 5.1 ± 0.02 GtC yr−1 (2.4 ± 0.01 ppm yr−1), SOCEAN 2.5 ± 0.6 GtC yr−1, and SLAND 3.4 ± 0.9 GtC yr−1, with a budgetimbalance BIM of −0.1 GtC yr−1 indicating a near balance betweenestimated sources and sinks over the last decade. For the year 2019 alone, thegrowth in EFOS was only about 0.1 % with fossil emissions increasingto 9.9 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.7 ± 0.5 GtC yr−1 when cement carbonation sink is included), and ELUC was 1.8 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5 ± 0.9 GtC yr−1 (42.2 ± 3.3 GtCO2). Also for 2019, GATM was5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEANwas 2.6 ± 0.6 GtC yr−1, and SLAND was 3.1 ± 1.2 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminarydata for 2020, accounting for the COVID-19-induced changes in emissions,suggest a decrease in EFOS relative to 2019 of about −7 % (medianestimate) based on individual estimates from four studies of −6 %, −7 %,−7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in thecomponents of the global carbon budget are consistently estimated over theperiod 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for therepresentation of semi-decadal variability in CO2 fluxes. Comparison ofestimates from diverse approaches and observations shows (1) no consensusin the mean and trend in land-use change emissions over the last decade, (2)a persistent low agreement between the different methods on the magnitude ofthe land CO2 flux in the northern extra-tropics, and (3) an apparentdiscrepancy between the different methods for the ocean sink outside thetropics, particularly in the Southern Ocean. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set (Friedlingstein et al.,2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014,2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2020 (Friedlingstein et al., 2020).

Abstract The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time‐evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

From 1953 to 2006, Arctic sea ice extent at the end of the melt season in September has declined sharply. All models participating in the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4) show declining Arctic ice cover over this period. However, depending on the time window for analysis, none or very few individual model simulations show trends comparable to observations. If the multi‐model ensemble mean time series provides a true representation of forced change by greenhouse gas (GHG) loading, 33–38% of the observed September trend from 1953–2006 is externally forced, growing to 47–57% from 1979–2006. Given evidence that as a group, the models underestimate the GHG response, the externally forced component may be larger. While both observed and modeled Antarctic winter trends are small, comparisons for summer are confounded by generally poor model performance.

Uncertainty in future climate change presents a key challenge for adaptation planning. In this study, uncertainty arising from internal climate variability is investigated using a new 40-member ensemble conducted with the National Center for Atmospheric Research Community Climate System Model Version 3 (CCSM3) under the SRES A1B greenhouse gas and ozone recovery forcing scenarios during 2000–2060. The contribution of intrinsic atmospheric variability to the total uncertainty is further examined using a 10,000-year control integration of the atmospheric model component of CCSM3 under fixed boundary conditions. The global climate response is characterized in terms of air temperature, precipitation, and sea level pressure during winter and summer. The dominant source of uncertainty in the simulated climate response at middle and high latitudes is internal atmospheric variability associated with the annular modes of circulation variability. Coupled ocean-atmosphere variability plays a dominant role in the tropics, with attendant effects at higher latitudes via atmospheric teleconnections. Uncertainties in the forced response are generally larger for sea level pressure than precipitation, and smallest for air temperature. Accordingly, forced changes in air temperature can be detected earlier and with fewer ensemble members than those in atmospheric circulation and precipitation. Implications of the results for detection and attribution of observed climate change and for multi-model climate assessments are discussed. Internal variability is estimated to account for at least half of the inter-model spread in projected climate trends during 2005–2060 in the CMIP3 multi-model ensemble.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate is critical to better understand the global carboncycle, support the development of climate policies, and project futureclimate change. Here we describe and synthesize data sets and methodologies toquantify the five major components of the global carbon budget and theiruncertainties. Fossil CO2 emissions (EFOS) are based on energystatistics and cement production data, while emissions from land-use change(ELUC), mainly deforestation, are based on land use and land-use changedata and bookkeeping models. Atmospheric CO2 concentration is measureddirectly, and its growth rate (GATM) is computed from the annualchanges in concentration. The ocean CO2 sink (SOCEAN) is estimatedwith global ocean biogeochemistry models and observation-baseddata products. The terrestrial CO2 sink (SLAND) is estimated withdynamic global vegetation models. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, withfossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission(including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1(40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with aBIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low orsinks were too high). The global atmospheric CO2 concentration averaged over2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest anincrease in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %)globally and atmospheric CO2 concentration reaching 417.2 ppm, morethan 50 % above pre-industrial levels (around 278 ppm). Overall, the meanand trend in the components of the global carbon budget are consistentlyestimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadalvariability in CO2 fluxes. Comparison of estimates from multipleapproaches and observations shows (1) a persistent large uncertainty in theestimate of land-use change emissions, (2) a low agreement between thedifferent methods on the magnitude of the land CO2 flux in the northernextratropics, and (3) a discrepancy between the different methods on thestrength of the ocean sink over the last decade. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set. The data presented inthis work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biosphere– the “global carbon budget” – is important to better understand theglobal carbon cycle, support the development of climate policies, andproject future climate change. Here we describe data sets and methodology toquantify the five major components of the global carbon budget and theiruncertainties. Fossil CO2 emissions (EFF) are based on energystatistics and cement production data, while emissions from land use change(ELUC), mainly deforestation, are based on land use and land use changedata and bookkeeping models. Atmospheric CO2 concentration is measureddirectly and its growth rate (GATM) is computed from the annual changesin concentration. The ocean CO2 sink (SOCEAN) and terrestrialCO2 sink (SLAND) are estimated with global process modelsconstrained by observations. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the lastdecade available (2009–2018), EFF was 9.5±0.5 GtC yr−1,ELUC 1.5±0.7 GtC yr−1, GATM 4.9±0.02 GtC yr−1 (2.3±0.01 ppm yr−1), SOCEAN 2.5±0.6 GtC yr−1, and SLAND 3.2±0.6 GtC yr−1, with a budgetimbalance BIM of 0.4 GtC yr−1 indicating overestimated emissionsand/or underestimated sinks. For the year 2018 alone, the growth in EFF wasabout 2.1 % and fossil emissions increased to 10.0±0.5 GtC yr−1, reaching 10 GtC yr−1 for the first time in history,ELUC was 1.5±0.7 GtC yr−1, for total anthropogenicCO2 emissions of 11.5±0.9 GtC yr−1 (42.5±3.3 GtCO2). Also for 2018, GATM was 5.1±0.2 GtC yr−1 (2.4±0.1 ppm yr−1), SOCEAN was 2.6±0.6 GtC yr−1, and SLAND was 3.5±0.7 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in EFF of +0.6 % (range of−0.2 % to 1.5 %) based on national emissions projections for China, theUSA, the EU, and India and projections of gross domestic product correctedfor recent changes in the carbon intensity of the economy for the rest ofthe world. Overall, the mean and trend in the five components of the globalcarbon budget are consistently estimated over the period 1959–2018, butdiscrepancies of up to 1 GtC yr−1 persist for the representation ofsemi-decadal variability in CO2 fluxes. A detailed comparison amongindividual estimates and the introduction of a broad range of observationsshows (1) no consensus in the mean and trend in land use change emissionsover the last decade, (2) a persistent low agreement between the differentmethods on the magnitude of the land CO2 flux in the northernextra-tropics, and (3) an apparent underestimation of the CO2variability by ocean models outside the tropics. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set (Le Quéré etal., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated bythis work are available at https://doi.org/10.18160/gcp-2019 (Friedlingsteinet al., 2019).

Abstract. Accurate assessment of anthropogenic carbon dioxide(CO2) emissions and their redistribution among the atmosphere,ocean, and terrestrial biosphere – the “global carbon budget” – isimportant to better understand the global carbon cycle, support thedevelopment of climate policies, and project future climate change. Here wedescribe data sets and methodology to quantify the five major components ofthe global carbon budget and their uncertainties. Fossil CO2emissions (EFF) are based on energy statistics and cementproduction data, while emissions from land use and land-use change (ELUC),mainly deforestation, are based on land use and land-use change data andbookkeeping models. Atmospheric CO2 concentration is measureddirectly and its growth rate (GATM) is computed from the annualchanges in concentration. The ocean CO2 sink (SOCEAN)and terrestrial CO2 sink (SLAND) are estimated withglobal process models constrained by observations. The resulting carbonbudget imbalance (BIM), the difference between the estimatedtotal emissions and the estimated changes in the atmosphere, ocean, andterrestrial biosphere, is a measure of imperfect data and understanding ofthe contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2008–2017), EFF was9.4±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.7±0.02 GtC yr−1,SOCEAN 2.4±0.5 GtC yr−1, and SLAND 3.2±0.8 GtC yr−1, with a budget imbalance BIM of0.5 GtC yr−1 indicating overestimated emissions and/or underestimatedsinks. For the year 2017 alone, the growth in EFF was about 1.6 %and emissions increased to 9.9±0.5 GtC yr−1. Also for 2017,ELUC was 1.4±0.7 GtC yr−1, GATM was 4.6±0.2 GtC yr−1, SOCEAN was 2.5±0.5 GtC yr−1, and SLAND was 3.8±0.8 GtC yr−1,with a BIM of 0.3 GtC. The global atmosphericCO2 concentration reached 405.0±0.1 ppm averaged over 2017.For 2018, preliminary data for the first 6–9 months indicate a renewedgrowth in EFF of +2.7 % (range of 1.8 % to 3.7 %) basedon national emission projections for China, the US, the EU, and India andprojections of gross domestic product corrected for recent changes in thecarbon intensity of the economy for the rest of the world. The analysispresented here shows that the mean and trend in the five components of theglobal carbon budget are consistently estimated over the period of 1959–2017,but discrepancies of up to 1 GtC yr−1 persist for the representationof semi-decadal variability in CO2 fluxes. A detailed comparisonamong individual estimates and the introduction of a broad range ofobservations show (1) no consensus in the mean and trend in land-use changeemissions, (2) a persistent low agreement among the different methods onthe magnitude of the land CO2 flux in the northern extra-tropics,and (3) an apparent underestimation of the CO2 variability by oceanmodels, originating outside the tropics. This living data update documentschanges in the methods and data sets used in this new global carbon budgetand the progress in understanding the global carbon cycle compared withprevious publications of this data set (Le Quéré et al., 2018, 2016,2015a, b, 2014, 2013). All results presented here can be downloaded fromhttps://doi.org/10.18160/GCP-2018.

fertilization effect. Six Earth system models have consistently projected continuous increases of VPD throughout the current century. Our results highlight that the impacts of VPD on vegetation growth should be adequately considered to assess ecosystem responses to future climate conditions.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate is critical to better understand the global carboncycle, support the development of climate policies, and project futureclimate change. Here we describe and synthesize datasets and methodology toquantify the five major components of the global carbon budget and theiruncertainties. Fossil CO2 emissions (EFOS) are based on energystatistics and cement production data, while emissions from land-use change(ELUC), mainly deforestation, are based on land use and land-use changedata and bookkeeping models. Atmospheric CO2 concentration is measureddirectly, and its growth rate (GATM) is computed from the annualchanges in concentration. The ocean CO2 sink (SOCEAN) is estimatedwith global ocean biogeochemistry models and observation-baseddata products. The terrestrial CO2 sink (SLAND) is estimated withdynamic global vegetation models. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the firsttime, an approach is shown to reconcile the difference in our ELUCestimate with the one from national greenhouse gas inventories, supportingthe assessment of collective countries' climate progress. For the year 2020, EFOS declined by 5.4 % relative to 2019, withfossil emissions at 9.5 ± 0.5 GtC yr−1 (9.3 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 0.9 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission of10.2 ± 0.8 GtC yr−1 (37.4 ± 2.9 GtCO2). Also, for2020, GATM was 5.0 ± 0.2 GtC yr−1 (2.4 ± 0.1 ppm yr−1), SOCEAN was 3.0 ± 0.4 GtC yr−1, and SLANDwas 2.9 ± 1 GtC yr−1, with a BIM of −0.8 GtC yr−1. Theglobal atmospheric CO2 concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in EFOSrelative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budgetare consistently estimated over the period 1959–2020, but discrepancies ofup to 1 GtC yr−1 persist for the representation of annual tosemi-decadal variability in CO2 fluxes. Comparison of estimates frommultiple approaches and observations shows (1) a persistent largeuncertainty in the estimate of land-use changes emissions, (2) a lowagreement between the different methods on the magnitude of the landCO2 flux in the northern extra-tropics, and (3) a discrepancy betweenthe different methods on the strength of the ocean sink over the lastdecade. This living data update documents changes in the methods and datasets used in this new global carbon budget and the progress in understandingof the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; LeQuéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). Thedata presented in this work are available at https://doi.org/10.18160/gcp-2021 (Friedlingstein et al., 2021).

The Community Land Model is the land component of the Community Climate System Model. Here, we describe a broad set of model improvements and additions that have been provided through the CLM development community to create CLM4. The model is extended with a carbon-nitrogen (CN) biogeochemical model that is prognostic with respect to vegetation, litter, and soil carbon and nitrogen states and vegetation phenology. An urban canyon model is added and a transient land cover and land use change (LCLUC) capability, including wood harvest, is introduced, enabling study of historic and future LCLUC on energy, water, momentum, carbon, and nitrogen fluxes. The hydrology scheme is modified with a revised numerical solution of the Richards equation and a revised ground evaporation parameterization that accounts for litter and within-canopy stability. The new snow model incorporates the SNow and Ice Aerosol Radiation model (SNICAR)- which includes aerosol deposition, grain-size dependent snow aging, and vertically-resolved snowpack heating – as well as new snow cover and snow burial fraction parameterizations. The thermal and hydrologic properties of organic soil are accounted for and the ground column is extended to,50-m depth. Several other minor modifications to the land surface types dataset, grass and

Abstract The Pacific decadal oscillation (PDO), the dominant year-round pattern of monthly North Pacific sea surface temperature (SST) variability, is an important target of ongoing research within the meteorological and climate dynamics communities and is central to the work of many geologists, ecologists, natural resource managers, and social scientists. Research over the last 15 years has led to an emerging consensus: the PDO is not a single phenomenon, but is instead the result of a combination of different physical processes, including both remote tropical forcing and local North Pacific atmosphere–ocean interactions, which operate on different time scales to drive similar PDO-like SST anomaly patterns. How these processes combine to generate the observed PDO evolution, including apparent regime shifts, is shown using simple autoregressive models of increasing spatial complexity. Simulations of recent climate in coupled GCMs are able to capture many aspects of the PDO, but do so based on a balance of processes often more independent of the tropics than is observed. Finally, it is suggested that the assessment of PDO-related regional climate impacts, reconstruction of PDO-related variability into the past with proxy records, and diagnosis of Pacific variability within coupled GCMs should all account for the effects of these different processes, which only partly represent the direct forcing of the atmosphere by North Pacific Ocean SSTs.

We conducted the first synchronously coupled atmosphere-ocean general circulation model simulation from the Last Glacial Maximum to the Bølling-Allerød (BA) warming. Our model reproduces several major features of the deglacial climate evolution, suggesting a good agreement in climate sensitivity between the model and observations. In particular, our model simulates the abrupt BA warming as a transient response of the Atlantic meridional overturning circulation (AMOC) to a sudden termination of freshwater discharge to the North Atlantic before the BA. In contrast to previous mechanisms that invoke AMOC multiple equilibrium and Southern Hemisphere climate forcing, we propose that the BA transition is caused by the superposition of climatic responses to the transient CO(2) forcing, the AMOC recovery from Heinrich Event 1, and an AMOC overshoot.

A consistent weakening and poleward expansion of the Hadley circulation is diagnosed in the climate change simulations of the IPCC AR4 project. Associated with this widening is a poleward expansion of the subtropical dry zone. Simple scaling analysis supports the notion that the poleward extent of the Hadley cell is set by the location where the thermally driven jet first becomes baroclinically unstable. The expansion of the Hadley cell is caused by an increase in the subtropical static stability, which pushes poleward the baroclinic instability zone and hence the outer boundary of the Hadley cell.

Abstract. Rises in surface and lower troposphere air temperatures through the 21st century are projected to be especially pronounced over the Arctic Ocean during the cold season. This Arctic amplification is largely driven by loss of the sea ice cover, allowing for strong heat transfers from the ocean to the atmosphere. Consistent with observed reductions in sea ice extent, fields from both the NCEP/NCAR and JRA-25 reanalyses point to emergence of surface-based Arctic amplification in the last decade.

Abstract Spatial variations in sea surface temperature (SST) and rainfall changes over the tropics are investigated based on ensemble simulations for the first half of the twenty-first century under the greenhouse gas (GHG) emission scenario A1B with coupled ocean–atmosphere general circulation models of the Geophysical Fluid Dynamics Laboratory (GFDL) and National Center for Atmospheric Research (NCAR). Despite a GHG increase that is nearly uniform in space, pronounced patterns emerge in both SST and precipitation. Regional differences in SST warming can be as large as the tropical-mean warming. Specifically, the tropical Pacific warming features a conspicuous maximum along the equator and a minimum in the southeast subtropics. The former is associated with westerly wind anomalies whereas the latter is linked to intensified southeast trade winds, suggestive of wind–evaporation–SST feedback. There is a tendency for a greater warming in the northern subtropics than in the southern subtropics in accordance with asymmetries in trade wind changes. Over the equatorial Indian Ocean, surface wind anomalies are easterly, the thermocline shoals, and the warming is reduced in the east, indicative of Bjerknes feedback. In the midlatitudes, ocean circulation changes generate narrow banded structures in SST warming. The warming is negatively correlated with wind speed change over the tropics and positively correlated with ocean heat transport change in the northern extratropics. A diagnostic method based on the ocean mixed layer heat budget is developed to investigate mechanisms for SST pattern formation. Tropical precipitation changes are positively correlated with spatial deviations of SST warming from the tropical mean. In particular, the equatorial maximum in SST warming over the Pacific anchors a band of pronounced rainfall increase. The gross moist instability follows closely relative SST change as equatorial wave adjustments flatten upper-tropospheric warming. The comparison with atmospheric simulations in response to a spatially uniform SST warming illustrates the importance of SST patterns for rainfall change, an effect overlooked in current discussion of precipitation response to global warming. Implications for the global and regional response of tropical cyclones are discussed.

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

A consistent poleward and upward shift and intensification of the storm tracks is found in an ensemble of 21st century climate simulations performed by 15 coupled climate models. The shift of the storm tracks is accompanied by a poleward shift and upward expansion of the midlatitude baroclinic regions associated with enhanced warming in the tropical upper troposphere and increased tropopause height. The poleward shift in baroclinicity is augmented in the Southern Hemisphere and partially offset in the Northern Hemisphere by changes in the surface meridional temperature gradient. The poleward shift of the storm tracks also tends to be accompanied by poleward shifts in surface wind stress and precipitation, and a shift towards the high index state of the annular modes. These results highlight the integral role that the storm tracks play in the climate system, and the importance of understanding how and why they will change in the future.