The authors present results for simulated climate and climate change from a newly developed high-resolution global climate model [Geophysical Fluid Dynamics Laboratory Climate Model version 2.5 (GFDL CM2.5)]. The GFDL CM2.5 has an atmospheric resolution of approximately 50 km in the horizontal, with 32 vertical levels. The horizontal resolution in the ocean ranges from 28 km in the tropics to 8 km at high latitudes, with 50 vertical levels. This resolution allows the explicit simulation of some mesoscale eddies in the ocean, particularly at lower latitudes. Analyses are presented based on the output of a 280-yr control simulation; also presented are results based on a 140-yr simulation in which atmospheric CO2 increases at 1% yr−1 until doubling after 70 yr. Results are compared to GFDL CM2.1, which has somewhat similar physics but a coarser resolution. The simulated climate in CM2.5 shows marked improvement over many regions, especially the tropics, including a reduction in the double ITCZ and an improved simulation of ENSO. Regional precipitation features are much improved. The Indian monsoon and Amazonian rainfall are also substantially more realistic in CM2.5. The response of CM2.5 to a doubling of atmospheric CO2 has many features in common with CM2.1, with some notable differences. For example, rainfall changes over the Mediterranean appear to be tightly linked to topography in CM2.5, in contrast to CM2.1 where the response is more spatially homogeneous. In addition, in CM2.5 the near-surface ocean warms substantially in the high latitudes of the Southern Ocean, in contrast to simulations using CM2.1.

The authors characterize impacts on heat in the ocean climate system from transient ocean mesoscale eddies. Their tool is a suite of centennial-scale 1990 radiatively forced numerical climate simulations from three GFDL coupled models comprising the Climate Model, version 2.0–Ocean (CM2-O), model suite. CM2-O models differ in their ocean resolution: CM2.6 uses a 0.1° ocean grid, CM2.5 uses an intermediate grid with 0.25° spacing, and CM2-1deg uses a nominal 1.0° grid. Analysis of the ocean heat budget reveals that mesoscale eddies act to transport heat upward in a manner that partially compensates (or offsets) for the downward heat transport from the time-mean currents. Stronger vertical eddy heat transport in CM2.6 relative to CM2.5 accounts for the significantly smaller temperature drift in CM2.6. The mesoscale eddy parameterization used in CM2-1deg also imparts an upward heat transport, yet it differs systematically from that found in CM2.6. This analysis points to the fundamental role that ocean mesoscale features play in transient ocean heat uptake. In general, the more accurate simulation found in CM2.6 provides an argument for either including a rich representation of the ocean mesoscale in model simulations of the mean and transient climate or for employing parameterizations that faithfully reflect the role of eddies in both lateral and vertical heat transport.

Tropical cyclones (TCs) are a hazard to life and property and a prominent element of the global climate system; therefore, understanding and predicting TC location, intensity, and frequency is of both societal and scientific significance. Methodologies exist to predict basinwide, seasonally aggregated TC activity months, seasons, and even years in advance. It is shown that a newly developed high-resolution global climate model can produce skillful forecasts of seasonal TC activity on spatial scales finer than basinwide, from months and seasons in advance of the TC season. The climate model used here is targeted at predicting regional climate and the statistics of weather extremes on seasonal to decadal time scales, and comprises high-resolution (50 km × 50 km) atmosphere and land components as well as more moderate-resolution (~100 km) sea ice and ocean components. The simulation of TC climatology and interannual variations in this climate model is substantially improved by correcting systematic ocean biases through “flux adjustment.” A suite of 12-month duration retrospective forecasts is performed over the 1981–2012 period, after initializing the climate model to observationally constrained conditions at the start of each forecast period, using both the standard and flux-adjusted versions of the model. The standard and flux-adjusted forecasts exhibit equivalent skill at predicting Northern Hemisphere TC season sea surface temperature, but the flux-adjusted model exhibits substantially improved basinwide and regional TC activity forecasts, highlighting the role of systematic biases in limiting the quality of TC forecasts. These results suggest that dynamical forecasts of seasonally aggregated regional TC activity months in advance are feasible.

We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution.

We describe the Geophysical Fluid Dynamics Laboratory's CM4.0 physical climate model, with emphasis on those aspects that may be of particular importance to users of this model and its simulations. The model is built with the AM4.0/LM4.0 atmosphere/land model and OM4.0 ocean model. Topics include the rationale for key choices made in the model formulation, the stability as well as drift of the preindustrial control simulation, and comparison of key aspects of the historical simulations with observations from recent decades. Notable achievements include the relatively small biases in seasonal spatial patterns of top‐of‐atmosphere fluxes, surface temperature, and precipitation; reduced double Intertropical Convergence Zone bias; dramatically improved representation of ocean boundary currents; a high‐quality simulation of climatological Arctic sea ice extent and its recent decline; and excellent simulation of the El Niño‐Southern Oscillation spectrum and structure. Areas of concern include inadequate deep convection in the Nordic Seas; an inaccurate Antarctic sea ice simulation; precipitation and wind composites still affected by the equatorial cold tongue bias; muted variability in the Atlantic Meridional Overturning Circulation; strong 100 year quasiperiodicity in Southern Ocean ventilation; and a lack of historical warming before 1990 and too rapid warming thereafter due to high climate sensitivity and strong aerosol forcing, in contrast to the observational record. Overall, CM4.0 scores very well in its fidelity against observations compared to the Coupled Model Intercomparison Project Phase 5 generation in terms of both mean state and modes of variability and should prove a valuable new addition for analysis across a broad array of applications.

A new high-resolution Geophysical Fluid Dynamics Laboratory (GFDL) coupled model [the High-Resolution Forecast-Oriented Low Ocean Resolution (FLOR) model (HiFLOR)] has been developed and used to investigate potential skill in simulation and prediction of tropical cyclone (TC) activity. HiFLOR comprises high-resolution (~25-km mesh) atmosphere and land components and a more moderate-resolution (~100-km mesh) sea ice and ocean component. HiFLOR was developed from FLOR by decreasing the horizontal grid spacing of the atmospheric component from 50 to 25 km, while leaving most of the subgrid-scale physical parameterizations unchanged. Compared with FLOR, HiFLOR yields a more realistic simulation of the structure, global distribution, and seasonal and interannual variations of TCs, as well as a comparable simulation of storm-induced cold wakes and TC-genesis modulation induced by the Madden–Julian oscillation (MJO). Moreover, HiFLOR is able to simulate and predict extremely intense TCs (Saffir–Simpson hurricane categories 4 and 5) and their interannual variations, which represents the first time a global coupled model has been able to simulate such extremely intense TCs in a multicentury simulation, sea surface temperature restoring simulations, and retrospective seasonal predictions.

The Intergovernmental Panel on Climate Change (IPCC) fifth assessment of projected global and regional ocean temperature change is based on global climate models that have coarse (∼100 km) ocean and atmosphere resolutions. In the Northwest Atlantic, the ensemble of global climate models has a warm bias in sea surface temperature due to a misrepresentation of the Gulf Stream position; thus, existing climate change projections are based on unrealistic regional ocean circulation. Here we compare simulations and an atmospheric CO2 doubling response from four global climate models of varying ocean and atmosphere resolution. We find that the highest resolution climate model (∼10 km ocean, ∼50 km atmosphere) resolves Northwest Atlantic circulation and water mass distribution most accurately. The CO2 doubling response from this model shows that upper‐ocean (0–300 m) temperature in the Northwest Atlantic Shelf warms at a rate nearly twice as fast as the coarser models and nearly three times faster than the global average. This enhanced warming is accompanied by an increase in salinity due to a change in water mass distribution that is related to a retreat of the Labrador Current and a northerly shift of the Gulf Stream. Both observations and the climate model demonstrate a robust relationship between a weakening Atlantic Meridional Overturning Circulation (AMOC) and an increase in the proportion of Warm‐Temperate Slope Water entering the Northwest Atlantic Shelf. Therefore, prior climate change projections for the Northwest Atlantic may be far too conservative. These results point to the need to improve simulations of basin and regional‐scale ocean circulation.

Responses of tropical cyclones (TCs) to CO 2 doubling are explored using coupled global climate models (GCMs) with increasingly refined atmospheric/land horizontal grids (~ 200 km, ~ 50 km and ~ 25 km). The three models exhibit similar changes in background climate fields thought to regulate TC activity, such as relative sea surface temperature (SST), potential intensity, and wind shear. However, global TC frequency decreases substantially in the 50 km model, while the 25 km model shows no significant change. The ~ 25 km model also has a substantial and spatially-ubiquitous increase of Category 3-4-5 hurricanes. Idealized perturbation experiments are performed to understand the TC response. Each model's transient fullycoupled 2 × CO 2 TC activity response is largely recovered by "time-slice" experiments using time-invariant SST perturbations added to each model's own SST climatology. The TC response to SST forcing depends on each model's background climatological SST biases: removing these biases leads to a global TC intensity increase in the ~ 50 km model, and a global TC frequency increase in the ~ 25 km model, in response to CO 2 -induced warming patterns and CO 2 doubling. Isolated CO 2 doubling leads to a significant TC frequency decrease, while isolated uniform SST warming leads to a significant global TC frequency increase; the ~ 25 km model has a greater tendency for frequency increase. Global TC frequency responds to both (1) changes in TC "seeds", which increase due to warming (more so in the ~ 25 km model) and decrease due to higher CO 2 concentrations, and (2) less efficient development of these"seeds" into TCs, largely due to the nonlinear relation between temperature and saturation specific humidity.