Because of a spring predictability barrier, the seasonal forecast skill of Arctic summer sea ice is limited by the availability of melt-season sea ice thickness (SIT) observations. The first year-round SIT observations, retrieved from CryoSat-2 from 2011 to 2020, are assimilated into the GFDL ocean–sea ice model. The model's SIT anomaly field is brought into significantly better agreement with the observations, particularly in the Central Arctic. Although the short observational period makes forecast assessment challenging, we find that the addition of May–August SIT assimilation improves September local sea ice concentration (SIC) and extent forecasts similarly to SIC-only assimilation. Although most regional forecasts are improved by SIT assimilation, the Chukchi Sea forecasts are degraded. This degradation is likely due to the introduction of negative correlations between September SIC and earlier SIT introduced by SIT assimilation, contrary to the increased correlations found in other regions.
We use two coupled climate models, GFDL-CM4 and GFDL-ESM4, to investigate the physical response of the Southern Ocean to changes in surface wind stress, Antarctic meltwater, and the combined forcing of the two in a pre-industrial control simulation. The meltwater cools the ocean surface in all regions except the Weddell Sea, where the wind stress warms the near-surface layer. The limited sensitivity of the Weddell Sea surface layer to the meltwater is due to the spatial distribution of the meltwater fluxes, regional bathymetry, and large-scale circulation patterns. The meltwater forcing dominates the Antarctic shelf response and the models yield strikingly different responses along West Antarctica. The disagreement is attributable to the mean-state representation and meltwater-driven acceleration of the Antarctic Slope Current (ASC). In CM4, the meltwater is efficiently trapped on the shelf by a well resolved, strong, and accelerating ASC which isolates the West Antarctic shelf from warm offshore waters, leading to strong subsurface cooling. In ESM4, a weaker and diffuse ASC allows more meltwater to escape to the open ocean, the West Antarctic shelf does not become isolated, and instead strong subsurface warming occurs. The CM4 results suggest a possible negative feedback mechanism that acts to limit future melting, while the ESM4 results suggest a possible positive feedback mechanism that acts to accelerate melt. Our results demonstrate the strong influence the ASC has on governing changes along the shelf, highlighting the importance of coupling interactive ice sheet models to ocean models that can resolve these dynamical processes.
Research over the past decade has demonstrated that dynamical forecast systems can skillfully predict pan-Arctic sea ice extent (SIE) on the seasonal time scale; however, there have been fewer assessments of prediction skill on user-relevant spatial scales. In this work, we evaluate regional Arctic SIE predictions made with the Forecast-Oriented Low Ocean Resolution (FLOR) and Seamless System for Prediction and Earth System Research (SPEAR_MED) dynamical seasonal forecast systems developed at the NOAA/Geophysical Fluid Dynamics Laboratory. Compared to FLOR, we find that the recently developed SPEAR_MED system displays improved skill in predicting regional detrended SIE anomalies, partially owing to improvements in sea ice concentration (SIC) and thickness (SIT) initial conditions. In both systems, winter SIE is skillfully predicted up to 11 months in advance, whereas summer minimum SIE predictions are limited by the Arctic spring predictability barrier, with typical skill horizons of roughly 4 months. We construct a parsimonious set of simple statistical prediction models to investigate the mechanisms of sea ice predictability in these systems. Three distinct predictability regimes are identified: a summer regime dominated by SIE and SIT anomaly persistence; a winter regime dominated by SIE and upper-ocean heat content (uOHC) anomaly persistence; and a combined regime in the Chukchi Sea, characterized by a trade-off between uOHC-based and SIT-based predictability that occurs as the sea ice edge position evolves seasonally. The combination of regional SIE, SIT, and uOHC predictors is able to reproduce the SIE skill of the dynamical models in nearly all regions, suggesting that these statistical predictors provide a stringent skill benchmark for assessing seasonal sea ice prediction systems.
The continuing decline of the summertime sea ice cover has reduced the sea ice path that must be traversed to Arctic destinations and through the Arctic between the Atlantic and Pacific Oceans, stimulating interest in trans–Arctic Ocean routes. Seasonal prediction of the sea ice cover along these routes could support the increasing summertime ship traffic taking advantage of recent low ice conditions. We introduce the minimum Arctic sea ice path (MIP) between Atlantic and Pacific Oceans as a shipping-relevant metric that is amenable to multidecadal hindcast evaluation. We show, using 1992–2017 retrospective predictions, that bias correction is necessary for the GFDL Seamless System for Prediction and Earth System Research (SPEAR) forecast system to improve upon damped persistence seasonal forecasts of summertime daily MIP between the Atlantic and Pacific Oceans both east and west of Greenland, corresponding roughly to the Northeast and Northwest Passages. Without bias correction, only the Northwest Passage MIP forecasts have lower error than a damped persistence forecast. Using the forecast ensemble spread to estimate a lower bound on forecast error, we find large opportunities for forecast error reduction, especially at lead times of less than 2 months. Most of the potential improvement remains after linear removal of climatological and trend biases, suggesting that significant error reduction might come from improved initialization and simulation of subannual variability. Using a different passive microwave sea ice dataset for calculating error than was used for data assimilation increases the raw forecast errors but not the trend anomaly forecast errors.
The current GFDL seasonal prediction system, the Seamless System for Prediction and Earth System Research (SPEAR), has shown skillful prediction of Arctic sea ice extent with atmosphere and ocean constrained by observations. In this study we present improvements in subseasonal and seasonal predictions of Arctic sea ice by directly assimilating sea ice observations. The sea ice initial conditions from a data assimilation (DA) system that assimilates satellite sea ice concentration (SIC) observations are used to produce a set of reforecast experiments (IceDA) starting from the first day of each month from 1992 to 2017. Our evaluation of daily sea ice extent prediction skill concludes that the SPEAR system generally outperforms the anomaly persistence forecast at lead times beyond 1 month. We primarily focus our analysis on daily gridcell-level sea ice fields. SIC DA improves prediction skill of SIC forecasts prominently in the June-, July-, August-, and September-initialized reforecasts. We evaluate two additional user-oriented metrics: the ice-free probability (IFP) and ice-free date (IFD). IFP is the probability of a grid cell experiencing ice-free conditions in a given year, and IFD is the first date on which a grid cell is ice free. A combined analysis of IFP and IFD demonstrates that the SPEAR model can make skillful predictions of local ice melt as early as May, with modest improvements from SIC DA.
Compared to the Arctic, seasonal predictions of Antarctic sea ice have received relatively little attention. In this work, we utilize three coupled dynamical prediction systems developed at the Geophysical Fluid Dynamics Laboratory to assess the seasonal prediction skill and predictability of Antarctic sea ice. These systems, based on the FLOR, SPEAR_LO, and SPEAR_MED dynamical models, differ in their coupled model components, initialization techniques, atmospheric resolution, and model biases. Using suites of retrospective initialized seasonal predictions spanning 1992–2018, we investigate the role of these factors in determining Antarctic sea ice prediction skill and examine the mechanisms of regional sea ice predictability. We find that each system is capable of skillfully predicting regional Antarctic sea ice extent (SIE) with skill that exceeds a persistence forecast. Winter SIE is skillfully predicted 11 months in advance in the Weddell, Amundsen/Bellingshausen, Indian, and west Pacific sectors, whereas winter skill is notably lower in the Ross sector. Zonally advected upper-ocean heat content anomalies are found to provide the crucial source of prediction skill for the winter sea ice edge position. The recently developed SPEAR systems are more skillful than FLOR for summer sea ice predictions, owing to improvements in sea ice concentration and sea ice thickness initialization. Summer Weddell SIE is skillfully predicted up to 9 months in advance in SPEAR_MED, due to the persistence and drift of initialized sea ice thickness anomalies from the previous winter. Overall, these results suggest a promising potential for providing operational Antarctic sea ice predictions on seasonal time scales.
Lockwood, Joseph W., Carolina O Dufour, Stephen M Griffies, and Michael Winton, April 2021: On the role of the Antarctic Slope Front on the occurrence of the Weddell Sea polynya under climate change. Journal of Climate, 34(7), DOI:10.1175/JCLI-D-20-0069.12529-2548. Abstract
This study investigates the occurrence of the Weddell Sea polynya under an idealized climate change scenario by evaluating simulations from climate models of different ocean resolutions. The GFDL-CM2.6 climate model, with roughly 3.8-km horizontal ocean grid spacing in the high latitudes, forms a Weddell Sea polynya at similar time and duration under idealized climate change forcing as under preindustrial forcing. In contrast, all convective models forming phase 5 of the Coupled Model Intercomparison Project (CMIP5) show either a cessation or a slowdown of Weddell Sea polynya events under climate warming. The representation of the Antarctic Slope Current and related Antarctic Slope Front is found to be key in explaining the differences between the two categories of models, with these features being more realistic in CM2.6 than in CMIP5. In CM2.6, the freshwater input driven by sea ice melt and enhanced runoff found under climate warming largely remains on the shelf region since the slope front restricts the lateral spread of the freshwater. In contrast, for most CMIP5 models, open-ocean stratification is enhanced by freshening since the absence of a slope front allows coastal freshwater anomalies to spread into the open ocean. This enhanced freshening contributes to the slowdown the occurrence of Weddell Sea polynyas. Hence, an improved representation of Weddell Sea shelf processes in current climate models is desirable to increase our ability to predict the fate of the Weddell Sea polynyas under climate change.
The current GFDL seasonal prediction system achieved retrospective sea ice extent (SIE) skill without direct sea ice data assimilation. Here we develop sea ice data assimilation, shown to be a key source of skill for seasonal sea ice predictions, in GFDL’s next-generation prediction system, the Seamless System for Prediction and Earth System Research (SPEAR). Satellite sea ice concentration (SIC) observations are assimilated into the GFDL Sea Ice Simulator version 2 (SIS2) using the ensemble adjustment Kalman filter (EAKF). Sea ice physics is perturbed to form an ensemble of ice–ocean members with atmospheric forcing from the JRA-55 reanalysis. Assimilation is performed every 5 days from 1982 to 2017 and the evaluation is conducted at pan-Arctic and regional scales over the same period. To mitigate an assimilation overshoot problem and improve the analysis, sea surface temperatures (SSTs) are restored to the daily Optimum Interpolation Sea Surface Temperature version 2 (OISSTv2). The combination of SIC assimilation and SST restoring reduces analysis errors to the observational error level (~10%) from up to 3 times larger than this (~30%) in the free-running model. Sensitivity experiments show that the choice of assimilation localization half-width (190 km) is near optimal and that SIC analysis errors can be further reduced slightly either by reducing the observational error or by increasing the assimilation frequency from every 5 days to daily. A lagged-correlation analysis suggests substantial prediction skill improvements from SIC initialization at lead times of less than 2 months.
The Southern Ocean south of 30° S represents only one-third of the total ocean area, yet absorbs half of the total ocean anthropogenic carbon and over two-thirds of ocean anthropogenic heat. In the past, the Southern Ocean has also been one of the most sparsely measured regions of the global ocean. Here we use pre-2005 ocean shipboard measurements alongside novel observations from autonomous floats with biogeochemical sensors to calculate changes in Southern Ocean temperature, salinity, pH and concentrations of nitrate, dissolved inorganic carbon and oxygen over two decades. We find local warming of over 3 °C, salinification of over 0.2 psu near the Antarctic coast, and isopycnals are found to deepen between 65° and 40° S. We find deoxygenation along the Antarctic coast, but reduced deoxygenation and nitrate concentrations where isopycnals deepen farther north. The forced response of the Earth system model ESM2M does not reproduce the observed patterns. Accounting for meltwater and poleward-intensifying winds in ESM2M improves reproduction of the observed large-scale changes, demonstrating the importance of recent changes in wind and meltwater. Future Southern Ocean biogeochemical changes are likely to be influenced by the relative strength of meltwater input and poleward-intensifying winds. The combined effect could lead to increased Southern Ocean deoxygenation and nutrient accumulation, starving the global ocean of nutrients sooner than otherwise expected.
The decline of Arctic sea‐ice extent has created a pressing need for accurate seasonal predictions of regional summer sea ice. Recent work has shown evidence for an Arctic sea ice spring predictability barrier, which may impose a sharp limit on regional forecasts initialized prior to spring. However, the physical mechanism for this barrier has remained elusive. In this work, we perform a daily sea‐ice mass (SIM) budget analysis in large ensemble experiments from two global climate models to investigate the mechanisms that underpin the spring predictability barrier. We find that predictability is limited in winter months by synoptically‐driven SIM export and negative feedbacks from sea‐ice growth. The spring barrier results from a sharp increase in predictability at melt onset, when ice‐albedo feedbacks act to enhance and persist the preexisting export‐generated mass anomaly. These results imply that ice‐thickness observations collected after melt onset are particularly critical for summer Arctic sea‐ice predictions.
We compare equilibrium climate sensitivity (ECS) estimates from pairs of long (≥800‐year) control and abruptly quadrupled CO2 simulations with shorter (150‐ and 300‐year) coupled atmosphere‐ocean simulations and slab ocean models (SOMs). Consistent with previous work, ECS estimates from shorter coupled simulations based on annual averages for years 1–150 underestimate those from SOM (−8% ± 13%) and long (−14% ± 8%) simulations. Analysis of only years 21–150 improved agreement with SOM (−2% ± 14%) and long (−8% ± 10%) estimates. Use of pentadal averages for years 51–150 results in improved agreement with long simulations (−4% ± 11%). While ECS estimates from current generation U.S. models based on SOM and coupled annual averages of years 1–150 range from 2.6°C to 5.3°C, estimates based longer simulations of the same models range from 3.2°C to 7.0°C. Such variations between methods argues for caution in comparison and interpretation of ECS estimates across models.
We describe the baseline coupled model configuration and simulation characteristics of GFDL's Earth System Model Version 4.1 (ESM4.1), which builds on component and coupled model developments at GFDL over 2013–2018 for coupled carbon‐chemistry‐climate simulation contributing to the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's CM4.0 development effort that focuses on ocean resolution for physical climate, ESM4.1 focuses on comprehensiveness of Earth system interactions. ESM4.1 features doubled horizontal resolution of both atmosphere (2° to 1°) and ocean (1° to 0.5°) relative to GFDL's previous‐generation coupled ESM2‐carbon and CM3‐chemistry models. ESM4.1 brings together key representational advances in CM4.0 dynamics and physics along with those in aerosols and their precursor emissions, land ecosystem vegetation and canopy competition, and multiday fire; ocean ecological and biogeochemical interactions, comprehensive land‐atmosphere‐ocean cycling of CO2, dust and iron, and interactive ocean‐atmosphere nitrogen cycling are described in detail across this volume of JAMES and presented here in terms of the overall coupling and resulting fidelity. ESM4.1 provides much improved fidelity in CO2 and chemistry over ESM2 and CM3, captures most of CM4.0's baseline simulations characteristics, and notably improves on CM4.0 in (1) Southern Ocean mode and intermediate water ventilation, (2) Southern Ocean aerosols, and (3) reduced spurious ocean heat uptake. ESM4.1 has reduced transient and equilibrium climate sensitivity compared to CM4.0. Fidelity concerns include (1) moderate degradation in sea surface temperature biases, (2) degradation in aerosols in some regions, and (3) strong centennial scale climate modulation by Southern Ocean convection.
Simulation of coupled carbon‐climate requires representation of ocean carbon cycling, but the computational burden of simulating the dozens of prognostic tracers in state‐of‐the‐art biogeochemistry ecosystem models can be prohibitive. We describe a six‐tracer biogeochemistry module of steady‐state phytoplankton and zooplankton dynamics in Biogeochemistry with Light, Iron, Nutrients and Gas (BLING version 2) with particular emphasis on enhancements relative to the previous version and evaluate its implementation in Geophysical Fluid Dynamics Laboratory's (GFDL) fourth‐generation climate model (CM4.0) with ¼° ocean. Major geographical and vertical patterns in chlorophyll, phosphorus, alkalinity, inorganic and organic carbon, and oxygen are well represented. Major biases in BLINGv2 include overly intensified production in high‐productivity regions at the expense of productivity in the oligotrophic oceans, overly zonal structure in tropical phosphorus, and intensified hypoxia in the eastern ocean basins as is typical in climate models. Overall, while BLINGv2 structural limitations prevent sophisticated application to plankton physiology, ecology, or biodiversity, its ability to represent major organic, inorganic, and solubility pumps makes it suitable for many coupled carbon‐climate and biogeochemistry studies including eddy interactions in the ocean interior. We further overview the biogeochemistry and circulation mechanisms that shape carbon uptake over the historical period. As an initial analysis of model historical and idealized response, we show that CM4.0 takes up slightly more anthropogenic carbon than previous models in part due to enhanced ventilation in the absence of an eddy parameterization. The CM4.0 biogeochemistry response to CO2 doubling highlights a mix of large declines and moderate increases consistent with previous models.
Recent laboratory and field studies point to an increase of sea salt aerosol (SSA) emissions with temperature, suggesting that SSA may lower climate sensitivity. We assess the impact of a strong (4.2 % K‐1) and weak (0.7% K‐1) temperature response of SSA emissions on the climate sensitivity of the coupled climate model CM4. We find that the stronger temperature dependence improves the simulation of marine aerosol optical depth sensitivity to temperature and lowers CM4 Transient Climate Response (‐0.12K) and Equilibrium Climate Sensitivity (‐0.5K). At CO2 doubling, the higher SSA emission sensitivity causes a negative radiative feedback (‐0.125 W m‐2 K‐1), which can only be partly explained by changes in the radiative effect of SSA (‐0.08 W m‐2 K‐1). Stronger radiative feedbacks are dominated by more negative low‐level clouds feedbacks in the Northern Hemisphere, which are partly offset by more positive feedbacks in the Southern Hemisphere associated with a weaker Atlantic Meridional Overturning Circulation.
GFDL's new CM4.0 climate model has high transient and equilibrium climate sensitivities near the middle of the upper half of CMIP5 models. The CMIP5 models have been criticized for excessive sensitivity based on observations of present‐day warming and heat uptake and estimates of radiative forcing. An ensemble of historical simulations with CM4.0 produces warming and heat uptake that are consistent with these observations under forcing that is at the middle of the assessed distribution. Energy budget‐based methods for estimating sensitivities based on these quantities underestimate CM4.0's sensitivities when applied to its historical simulations. However, we argue using a simple attribution procedure that CM4.0's warming evolution indicates excessive transient sensitivity to greenhouse gases. This excessive sensitivity is offset prior to recent decades by excessive response to aerosol and land use changes.
Storm surge and coastal flooding caused by tropical cyclones (hurricanes) and extratropical cyclones (nor'easters) pose a threat to communities along the Atlantic coast of the United States. Climate change and sea level rise are altering the statistics of these extreme events in a rather complex fashion. Here we use a fully-coupled global weather/climate modeling system (GFDL CM4) to study characteristics of extreme daily sea level (ESL) along the US Atlantic coast and their response to global warming. We find that under natural weather processes, the Gulf of Mexico coast is most vulnerable to storm surge and related ESL. New Orleans is a striking hotspot with the highest surge efficiency in response to storm winds. Under a 1% per year atmospheric CO2 increase on centennial time scales, the anthropogenic signal in ESL is robust along the US East Coast. It can emerge from the background variability as soon as in twenty years, or even before global sea level rise is taken into account. The regional dynamic sea level rise induced by the weakening of the Atlantic meridional overturning circulation facilitates this early emergence, especially during wintertime coastal flooding associated with nor’easters. Along the Gulf Coast, ESL is sensitive to the modification of hurricane characteristics under the CO2 forcing.
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 (CMIP6/OMIP). The ocean component of OM4 uses version 6 of the Modular Ocean Model (MOM6) and the sea‐ice component uses version 2 of the Sea Ice Simulator (SIS2), which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments (CORE) 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.
MOM6 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 mid‐depth 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.
Seasonal forecast systems can skillfully predict summer Arctic sea‐ice up to four months in advance. For some regions, however, there is a springtime predictability barrier that causes forecasts initialized prior to May to be less skillful. Since this barrier has only been documented in a few general circulation models (GCMs), we evaluate GCMs participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5). We first show sea‐ice volume skillfully predicts summer sea ice‐area (SIA) and has similar skill to a perfect model experiment. Given this result, we assess regional SIA predictability across CMIP5 and find a universal predictability barrier in late spring. For SIA at each summer target month in the marginal seas of the Arctic basin, a notable drop in prediction skill occurs from June to May in each GCM. This suggests summer sea‐ice forecasts initialized after June 1 will have better prediction skill than forecasts initialized before.
Seasonal predictions of Arctic sea ice on regional spatial scales are a pressing need for a broad group of stakeholders, however, most assessments of predictability and forecast skill to date have focused on pan-Arctic sea–ice extent (SIE). In this work, we present the first direct comparison of perfect model (PM) and operational (OP) seasonal prediction skill for regional Arctic SIE within a common dynamical prediction system. This assessment is based on two complementary suites of seasonal prediction ensemble experiments performed with a global coupled climate model. First, we present a suite of PM predictability experiments with start dates spanning the calendar year, which are used to quantify the potential regional SIE prediction skill of this system. Second, we assess the system’s OP prediction skill for detrended regional SIE using a suite of retrospective initialized seasonal forecasts spanning 1981–2016. In nearly all Arctic regions and for all target months, we find a substantial skill gap between PM and OP predictions of regional SIE. The PM experiments reveal that regional winter SIE is potentially predictable at lead times beyond 12 months, substantially longer than the skill of their OP counterparts. Both the OP and PM predictions display a spring prediction skill barrier for regional summer SIE forecasts, indicating a fundamental predictability limit for summer regional predictions. We find that a similar barrier exists for pan-Arctic sea–ice volume predictions, but is not present for predictions of pan-Arctic SIE. The skill gap identified in this work indicates a promising potential for future improvements in regional SIE predictions.
Dynamical prediction systems have shown potential to meet the emerging need for seasonal forecasts of regional Arctic sea ice. Observationally constrained initial conditions are a key source of skill for these predictions, but the direct influence of different observation types on prediction skill has not yet been systematically investigated. In this work, we perform a hierarchy of Observing System Experiments with a coupled global data assimilation and prediction system to assess the value of different classes of oceanic and atmospheric observations for seasonal sea-ice predictions in the Barents Sea. We find notable skill improvements due to the inclusion of both sea-surface temperature (SST) satellite observations and subsurface conductivity-temperature-depth (CTD) measurements. The SST data is found to provide the crucial source of interannual variability, whereas the CTD data primarily provide climatological and trend improvements. Analysis of the Barents Sea ocean heat budget suggests that ocean heat content anomalies in this region are driven by surface heat fluxes on seasonal timescales.
We describe GFDL'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 pre‐industrial 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 quasi‐periodicity 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.
Meltwater from the Antarctic Ice Sheet is projected to cause up to one metre of sea-level rise by 2100 under the highest greenhouse gas concentration trajectory (RCP8.5) considered by the Intergovernmental Panel on Climate Change (IPCC). However, the effects of meltwater from the ice sheets and ice shelves of Antarctica are not included in the widely used CMIP5 climate models, which introduces bias into IPCC climate projections. Here we assess a large ensemble simulation of the CMIP5 model ‘GFDL ESM2M’ that accounts for RCP8.5-projected Antarctic Ice Sheet meltwater. We find that, relative to the standard RCP8.5 scenario, accounting for meltwater delays the exceedance of the maximum global-mean atmospheric warming targets of 1.5 and 2 degrees Celsius by more than a decade, enhances drying of the Southern Hemisphere and reduces drying of the Northern Hemisphere, increases the formation of Antarctic sea ice (consistent with recent observations of increasing Antarctic sea-ice area) and warms the subsurface ocean around the Antarctic coast. Moreover, the meltwater-induced subsurface ocean warming could lead to further ice-sheet and ice-shelf melting through a positive feedback mechanism, highlighting the importance of including meltwater effects in simulations of future climate.
The past few years have seen a growing investment in the development of global eddy‐resolving ocean models, but the impact of incorporating such high ocean resolution on precipitation responses to CO2 forcing has yet to be investigated. This study analyzes precipitation changes from a suite of GFDL models incorporating eddy‐resolving (0.1o), eddy‐permitting (0.25o) and eddy‐parameterizing (1o) ocean models. The incorporation of eddy resolution does not challenge the large‐scale structure of precipitation changes but results in substantial regional differences, particularly over ocean. These oceanic differences are primarily driven by the pattern of SST changes with greater sensitivity in lower latitudes. The largest impact of ocean resolution on SST changes occurs in eddy rich regions (e.g., boundary currents and the Southern Ocean), where impact on precipitation changes is also found to various degrees. In the Gulf Stream region where previous studies found considerable impact of eddy resolution on the simulation of climatological precipitation, we do not find such impact from the GFDL models but we do find substantial impact on precipitation changes. The eddy‐parameterizing model projects a banded structure common to the CMIP5 models, whereas the higher‐resolution models project a poleward shift of precipitation maxima associated with an enhanced Gulf Stream warming. Over land, precipitation changes are generally not very sensitive to ocean resolution. In eastern North America adjacent to the Gulf Stream region, moderate differences are found between resolutions. We discuss the mechanisms of land differences, which arise through the simulation of both climatological SST and SST changes.
In this two-part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part I, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode – with prescribed sea surface temperatures (SSTs) and sea ice distribution – is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top-of-atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part II, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.
In Part II of this two-part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part I. Part II provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.
Previous studies found large biases between individual observational and model estimates of historical ocean anthropogenic carbon uptake. We show that the largest bias between the Coupled Model Intercomparison Project phase 5 (CMIP5) ensemble mean and between two observational estimates of ocean anthropogenic carbon is due to a difference in start date. After adjusting the CMIP5 and observational estimates to the 1791-1995 period, all three carbon uptake estimates agree to within 3 Pg of C, about 4% of the total. The CMIP5 ensemble mean spatial bias compared to the observations is generally smaller than the observational error, apart from a negative bias in the Southern Ocean, and a positive bias in the Southern Indian and Pacific Oceans compensating each other in the global mean. This dipole pattern is likely due to an equatorward and weak bias in the position of Southern Hemisphere westerlies and lack of mode and intermediate water ventilation.
Due to its persistence on seasonal timescales, Arctic sea-ice thickness (SIT) is a potential source of predictability for summer sea-ice extent (SIE). New satellite observations of SIT represent an opportunity to harness this potential predictability via improved thickness initialization in seasonal forecast systems. In this work, the evolution of Arctic sea-ice volume anomalies is studied using a 700-year control integration and a suite of initialized ensemble forecasts from a fully-coupled global climate model. Our analysis is focused on the September sea-ice zone, as this is the region where thickness anomalies have the potential to impact the SIE minimum. The primary finding of this paper is that, in addition to a general decay with time, sea-ice volume anomalies display a summer enhancement, in which anomalies tend to grow between the months of May and July. This summer enhancement is relatively symmetric for positive and negative volume anomalies and peaks in July regardless of the initial month. Analysis of the surface energy budget reveals that the summer volume anomaly enhancement is driven by a positive feedback between the SIT state and the surface albedo. The SIT state affects surface albedo through changes in the sea-ice concentration field, melt-onset date, snow coverage, and ice-thickness distribution, yielding an anomaly in the total absorbed shortwave radiation between May and August, which enhances the existing SIT anomaly. This phenomenon highlights the crucial importance of accurate SIT initialization and representation of ice-albedo feedback processes in seasonal forecast systems.
Recent Arctic sea ice seasonal prediction efforts and forecast skill assessments have primarily focused on pan-Arctic sea-ice extent (SIE). In this work, we move towards stakeholder-relevant spatial scales, investigating the regional forecast skill of Arctic sea ice in a Geophysical Fluid Dynamics Laboratory (GFDL) seasonal prediction system. Using a suite of retrospective initialized forecasts spanning 1981–2015 made with a coupled atmosphere-ocean-sea ice-land model, we show that predictions of detrended regional SIE are skillful at lead times up to 11 months. Regional prediction skill is highly region and target month dependent, and generically exceeds the skill of an anomaly persistence forecast. We show for the first time that initializing the ocean subsurface in a seasonal prediction system can yield significant regional skill for winter SIE. Similarly, as suggested by previous work, we find that sea-ice thickness initial conditions provide a crucial source of skill for regional summer SIE.
Dufour, Carolina O., Adele K Morrison, Stephen M Griffies, Ivy Frenger, Hannah Zanowski, and Michael Winton, October 2017: Preconditioning of the Weddell Sea polynya by the ocean mesoscale and dense water overflows. Journal of Climate, 30(19), DOI:10.1175/JCLI-D-16-0586.1. Abstract
TheWeddell Sea polynya is a large opening in the open-ocean sea ice cover associated with intense deep convection in the ocean. A necessary condition to form and maintain a polynya is the presence of a strong subsurface heat reservoir. This study investigates the processes that control the stratification and hence the build-up of the subsurface heat reservoir in theWeddell Sea. To do so, a climate model run for 200 years under preindustrial forcing with two eddying resolutions in the ocean (0.25° CM2.5 and 0.10° CM2.6) is investigated. Over the course of the simulation, CM2.6 develops two polynyas in the Weddell Sea, while CM2.5 exhibits quasi-continuous deep convection, but no polynyas, exemplifying that deep convection is not a sufficient condition for a polynya to occur. CM2.5 features a weaker subsurface heat reservoir than CM2.6 due to weak stratification associated with episodes of gravitational instability and enhanced vertical mixing of heat, resulting in an erosion of the reservoir. In contrast, in CM2.6, the water column is more stably stratified, allowing the subsurface heat reservoir to build up. The enhanced stratification in CM2.6 arises from its refined horizontal grid spacing and resolution of topography which allows, in particular, a better representation of the restratifying effect by transient mesoscale eddies and of the overflows of dense waters along the continental slope.
Goddard, P, Carolina O Dufour, Jianjun Yin, Stephen M Griffies, and Michael Winton, October 2017: CO2-Induced Ocean Warming of the Antarctic Continental Shelf in an Eddying Global Climate Model. Journal of Geophysical Research: Oceans, 122(10), DOI:10.1002/2017JC012849. Abstract
Ocean warming near the Antarctic ice shelves has critical implications for future ice sheet mass loss and global sea level rise. A global climate model with an eddying ocean is used to quantify the mechanisms contributing to ocean warming on the Antarctic continental shelf in an idealized 2xCO2 experiment. The results indicate that relatively large warm anomalies occur both in the upper 100 m and at depths above the shelf floor, which are controlled by different mechanisms. The near-surface ocean warming is primarily a response to enhanced onshore advective heat transport across the shelf break. The deep shelf warming is initiated by onshore intrusions of relatively warm Circumpolar Deep Water (CDW), in density classes that access the shelf, as well as the reduction of the vertical mixing of heat. CO2-induced shelf freshening influences both warming mechanisms. The shelf freshening slows vertical mixing by limiting gravitational instabilities and the upward diffusion of heat associated with CDW, resulting in the build-up of heat at depth. Meanwhile, freshening near the shelf break enhances the lateral density gradient of the Antarctic Slope Front (ASF) and disconnect isopycnals between the shelf and CDW, making cross-ASF heat exchange more difficult. However, at several locations along the ASF, the cross-ASF heat transport is less inhibited and heat can move onshore. Once onshore, lateral and vertical heat advection work to disperse the heat anomalies across the shelf region. Understanding the inhomogeneous Antarctic shelf warming will lead to better projections of future ice sheet mass loss.
There is large uncertainty in the simulation of transient climate sensitivity. This study aims to understand how such uncertainty is related to the simulation of the base climate by comparing two simulations with the same model but in which CO2 is increased from either a preindustrial (1860) or a present-day (1990) control simulation. This allows different base climate ocean circulations that are representative of those in current climate models to be imposed upon a single model. As a result, the model projects different transient climate sensitivities that are comparable to the multimodel spread. The greater warming in the 1990-start run occurs primarily at high latitudes and particularly over regions of oceanic convection. In the 1990-start run, ocean overturning circulations are initially weaker and weaken less from CO2 forcing. As a consequence, there are smaller reductions in the poleward ocean heat transport, leading to less tropical ocean heat storage and less moderated high-latitude surface warming. This process is evident in both hemispheres, with changes in the Atlantic meridional overturning circulation and the Antarctic Bottom Water formation dominating the warming differences in each hemisphere. The high-latitude warming in the 1990-start run is enhanced through albedo and cloud feedbacks, resulting in a smaller ocean heat uptake efficacy. The results highlight the importance of improving the base climate ocean circulation in order to provide a reasonable starting point for assessments of past climate change and the projection of future climate change.
Ice bridges are static structures composed of tightly packed sea ice that can form during the course of its
flow through a narrow strait. Despite their important role in local ecology and climate, the formation and
breakup of ice bridges is not well understood and has proved difficult to predict. Using long-wave
approximations and a continuum description of sea ice dynamics, we develop a one-dimensional theory for
the wind-driven formation of ice bridges in narrow straits, which is verified against direct numerical
simulations. We show that for a given wind stress and minimum and maximum channel widths, a steadystate
ice bridge can only form beyond a critical value of the thickness and the compactness of the ice field.
The theory also makes quantitative predictions for ice fluxes, which are particularly useful to estimate the
ice export associated with the breakup of ice bridges. We note that similar ideas are applicable to dense
granular flows in confined geometries.
Rallabandi, B, Z Zheng, Michael Winton, and H Stone, July 2017: Formation of sea ice bridges in narrow straits in response to wind and water stresses. Journal of Geophysical Research: Oceans, 122(7), DOI:10.1002/2017JC012822. Abstract
Ice bridges are rigid structures composed of sea ice that form seasonally in the many straits and channels of the Canadian Arctic Archipelago. Driven primarily by atmospheric stresses, these ice bridges are formed when sufficiently thick ice “jams” during the course of its flow between land masses, resulting in a region of stationary compacted ice that is separated from a region of flowing open water (a polynya) by a static arch. Using a continuum description of sea ice that is widely used in climate modeling, we present an asymptotic theory of the process of formation of such bridges in slender channels when the motion of the ice is driven by external wind and water stresses. We show that for an arbitrary channel shape, ice bridges can only form within a range of ice properties that is determined by the channel geometry and the external stress. We then compare the results of our theory with direct numerical simulations and observational evidence. Finally, we provide simple analytical expressions for the mean speed of the ice flow as a function of the channel shape, the properties of the ice, and the wind and water stresses along the channel.
Gregory, Jonathan M., N Bouttes-Mauhourat, Stephen M Griffies, Helmuth Haak, William J Hurlin, J H Jungclaus, M Kelley, W G Lee, J Marshall, Anastasia Romanou, Oleg A Saenko, Detlef Stammer, and Michael Winton, November 2016: The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) contribution to CMIP6: Investigation of sea-level and ocean climate change in response to CO2 forcing. Geoscientific Model Development, 9(11), DOI:10.5194/gmd-9-3993-2016. Abstract
The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) aims to investigate the spread in simulations of sea-level and ocean climate change in response to CO2 forcing by atmosphere-ocean general circulation models (AOGCMs). It is particularly motivated by the uncertainties in projections of ocean heat uptake, global-mean sea-level rise due to thermal expansion and the geographical patterns of sea-level change due to ocean density and circulation change. FAFMIP has three tier-1 experiments, in which prescribed surface flux perturbations of momentum, heat and freshwater respectively are applied to the ocean in separate AOGCM simulations. All other conditions are as in the pre-industrial control. The prescribed fields are typical of pattern and magnitude of changes in these fluxes projected by AOGCMs for doubled CO2 concentration. Five groups have tested the experimental design with existing AOGCMs. Their results show diversity in the pattern and magnitude of changes, with some common qualitative features. Heat and water flux perturbation cause the dipole in sea-level change in the North Atlantic, while momentum and heat flux perturbation cause the gradient across the Antarctic Circumpolar Current. The Atlantic Meridional Overturning Circulation (AMOC) declines in response to the heat flux perturbation, and there is a strong positive feedback on this effect due to the consequent cooling of sea surface temperature in the North Atlantic, which enhances the local heat input to the ocean. The momentum and water flux perturbations do not substantially affect the AMOC. Heat is taken up largely as a passive tracer in the Southern Ocean, which is the region of greatest heat input, but elsewhere heat is actively redistributed towards lower latitude. Future analysis of these and other phenomena with the wider range of CMIP6 FAFMIP AOGCMs will benefit from new diagnostics of temperature and salinity tendencies, which will enable investigation of the model spread in behaviour in terms of physical processes as formulated in the models.
Griffies, Stephen M., Gokhan Danabasoglu, Paul J Durack, Alistair Adcroft, V Balaji, C Böning, Eric P Chassignet, Enrique N Curchitser, Julie Deshayes, H Drange, Baylor Fox-Kemper, Peter J Gleckler, Jonathan M Gregory, Helmuth Haak, Robert Hallberg, Helene T Hewitt, David M Holland, Tatiana Ilyina, J H Jungclaus, Y Komuro, John P Krasting, William G Large, S J Marsland, S Masina, Trevor J McDougall, A J George Nurser, James C Orr, Anna Pirani, Fangli Qiao, Ronald J Stouffer, Karl E Taylor, Anne Marie Treguier, Hiroyuki Tsujino, P Uotila, M Valdivieso, Michael Winton, and Stephen G Yeager, September 2016: OMIP contribution to CMIP6: experimental and diagnostic protocol for the physical component of the Ocean Model Intercomparison Project. Geoscientific Model Development, 9(9), DOI:10.5194/gmd-9-3231-2016. Abstract
The Ocean Model Intercomparison Project (OMIP) aims to provide a framework for evaluating, understanding, and improving the ocean and sea-ice components of global climate and earth system models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses these aims in two complementary manners: (A) by providing an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing, (B) by providing a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) offering details for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows that of the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II have become the standard method to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP (Scenario MIP), as well as the ocean-sea ice OMIP simulations. The bulk of this paper offers scientific rationale for saving these diagnostics.
The Southern Ocean plays a dominant role in anthropogenic oceanic heat uptake. Strong northward transport of the heat content anomaly limits warming of the sea surface temperature in the uptake region and allows the heat uptake to be sustained. Using an eddy-rich global climate model, the processes controlling the northward transport and convergence of the heat anomaly in the mid-latitude Southern Ocean are investigated in an idealized 1% yr−1 increasing CO2 simulation. Heat budget analyses reveal that different processes dominate to the north and south of the main convergence region. The heat transport northward from the uptake region in the south is driven primarily by passive advection of the heat content anomaly by the existing time mean circulation, with a smaller 20% contribution from enhanced upwelling. The heat anomaly converges in the mid-latitude deep mixed layers, because there is not a corresponding increase in the mean heat transport out of the deep mixed layers northward into the mode waters. To the north of the deep mixed layers, eddy processes drive the warming and account for nearly 80% of the northward heat transport anomaly. The eddy transport mechanism results from a reduction in both the diffusive and advective southward eddy heat transports, driven by decreasing isopycnal slopes and decreasing along-isopycnal temperature gradients on the northern edge of the peak warming.
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.
We assess the uptake, transport and storage of oceanic anthropogenic carbon and
heat over the period 1861 to 2005 in a new set of coupled carbon-climate Earth
System models conducted for the fifth Coupled Model Intercomparison Project
(CMIP5), with a particular focus on the Southern Ocean. Simulations show the
Southern Ocean south of 30°S, occupying 30% of global surface ocean area, accounts
for 43 ± 3% (42 ± 5 Pg C) of anthropogenic CO2 and 75 ± 22% (23 ± 9 *1022J) of heat
uptake by the ocean over the historical period. Northward transport out of the Southern
Ocean is vigorous, reducing the storage to 33 ± 6 Pg anthropogenic carbon and 12 ± 7
*1022J heat in the region. The CMIP5 models as a class tend to underestimate the
observational-based global anthropogenic carbon storage, but simulate trends in global
ocean heat storage over the last fifty years within uncertainties of observation-based
estimates. CMIP5 models suggest global and Southern Ocean CO2 uptake have been
largely unaffected by recent climate variability and change. Anthropogenic carbon and
heat storage show a common broad-scale pattern of change, but ocean heat storage is
more structured than ocean carbon storage. Our results highlight the significance of
the Southern Ocean for the global climate and as the region where models differ the
most in representation of anthropogenic CO2 and in particular heat uptake.
We characterize impacts on heat in the ocean climate system from transient ocean mesoscale eddies. Our tool is a suite of centennial-scale 1990 radiatively forced numerical climate simulations from three GFDL coupled models comprising the 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 nominally 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.
Recent studies have suggested that global mean surface temperature would remain approximately constant on multi-century timescales after CO2 emissions are stopped. Here we use Earth system model simulations of such a stoppage to demonstrate that in some models, surface temperature may actually increase on multi-century timescales after an initial century-long decrease. This occurs in spite of a decline in radiative forcing that exceeds the decline in ocean heat uptake—a circumstance that would otherwise be expected to lead to a decline in global temperature. The reason is that the warming effect of decreasing ocean heat uptake together with feedback effects arising in response to the geographic structure of ocean heat uptake overcompensates the cooling effect of decreasing atmospheric CO2 on multi-century timescales. Our study also reveals that equilibrium climate sensitivity estimates based on a widely used method of regressing the Earth’s energy imbalance against surface temperature change are biased. Uncertainty in the magnitude of the feedback effects associated with the magnitude and geographic distribution of ocean heat uptake therefore contributes substantially to the uncertainty in allowable carbon emissions for a given multi-century warming target.
We provide an assessment of sea level simulated in a suite of global ocean-sea ice models using the interannual CORE atmospheric state to determine surface ocean boundary buoyancy and momentum fluxes. These CORE-II simulations are compared amongst themselves as well as to observation-based estimates. We focus on the final 15 years of the simulations (1993-2007), as this is a period where the CORE-II atmospheric state is well sampled, and it allows us to compare sea level related fields to both satellite and in situ analyses. The ensemble mean of the CORE-II simulations broadly agree with various global and regional observation-based analyses during this period, though with the global mean thermosteric sea level rise biased low relative to observation-based analyses. The simulations reveal a positive trend in dynamic sea level in the west Pacific and negative trend in the east, with this trend arising from wind shifts and regional changes in upper 700 m ocean heat content. The models also exhibit a thermosteric sea level rise in the subpolar North Atlantic associated with a transition around 1995/1996 of the North Atlantic Oscillation to its negative phase, and the advection of warm subtropical waters into the subpolar gyre. Sea level trends are predominantly associated with steric trends, with thermosteric effects generally far larger than halosteric effects, except in the Arctic and North Atlantic. There is a general anti-correlation between thermosteric and halosteric effects for much of the World Ocean, associated with density compensated changes.
We present seasonal predictions of Arctic sea ice extent (SIE) over the 1982-2013 period using two suites of retrospective forecasts initialized from a fully coupled ocean-atmosphere-sea ice assimilation system. High skill scores are found in predicting year-to-year fluctuations of Arctic SIE, with significant correlations up to 7 month ahead for September detrended anomalies. Predictions over the recent era, which coincides with an improved observational coverage, outperform the earlier period for most target months. We find, however, a degradation of skill in September during the last decade, a period of sea ice thinning in observations. The two prediction models, CM2.1 and FLOR, share very similar ocean and ice component and initialization but differ by their atmospheric component. FLOR has improved climatological atmospheric circulation and sea ice mean state but its skill is overall similar to CM2.1 for most seasons, which suggests a key role for initial conditions in predicting seasonal SIE fluctuations.
We investigate the influence of ocean component resolution on simulation of climate sensitivity using variants of the GFDL CM2.5 climate model incorporating eddy-resolving (1/10o) and eddy-parameterizing (1o) ocean resolutions. Two parameterization configurations of the coarse-resolution model are used yielding a three-model suite with significant variation in the transient climate response (TCR). The variation of TCR in this suite and in an enhanced group of 10 GFDL models is found to be strongly associated with the control climate Atlantic meridional overturning circulation (AMOC) magnitude and its decline under forcing. We find it is the AMOC behavior rather than resolution per se that accounts for most of the TCR differences. A smaller difference in TCR stems from the eddy-resolving model having more Southern Ocean surface warming than the coarse models.
Climate models simulate a wide range of climate changes at high northern latitudes in response to increased CO2. They also have substantial disagreement on projected changes of the Atlantic meridional overturning circulation (AMOC). Here we use two pairs of closely related climate models - each containing members with large and small AMOC declines - to explore the influence of AMOC decline on the high latitude response to increased CO2. The models with larger AMOC decline have less high latitude warming and sea ice decline than their small AMOC decline counterpart. By examining differences in the perturbation heat budget of the 40�90�N region, it is shown that AMOC decline diminishes the warming by weakening poleward ocean heat transport and increasing the ocean heat uptake. The cooling impact of this AMOC forced surface heat flux perturbation difference is enhanced by shortwave feedback and diminished by longwave feedback and atmospheric heat transport differences. The magnitude of the AMOC decline within model pairs is positively related to the magnitudes of control climate AMOC and Labrador Sea convection. Because the 40degree 90degree N region accounts for up to 40% of the simulated global ocean heat uptake over one hundred years, the process described here influences the global heat uptake efficiency.
We examine the influence of alternative ocean and atmosphere subcomponents on climate model simulation of transient sensitivities by comparing three GFDL climate models used for the CMIP5. The base model ESM2M is closely related to GFDL's CMIP3 climate model CM2.1, and makes use of a depth coordinate ocean component. The second model, ESM2G, is identical to ESM2M but makes use of an isopycnal coordinate ocean model. We compare the impact of this "ocean swap" with an "atmosphere swap" that produces the CM3 climate model by replacing the AM2 atmosphere with AM3 while retaining a depth coordinate ocean model. The atmosphere swap is found to have much larger influence on sensitivities of global surface temperature and Northern Hemisphere sea ice cover. The atmosphere swap also introduces a multi-decadal response timescale through its indirect influence on heat uptake. Despite significant differences in their interior ocean mean states, the ESM2M and ESM2G simulations of these metrics of climate change are very similar, except for an enhanced high latitude salinity response accompanied by temporarily advancing sea ice in ESM2G. In the ESM2G historical simulation this behavior results in the establishment of a strong halocline in the subpolar North Atlantic during the early 20th century and an associated cooling which are counter to observations in that region. The Atlantic meridional overturning declines comparably in all three models.
The influence of changing ocean currents on climate change is evaluated by comparing an earth system model’s response to increased CO2 with and without an ocean circulation response. Inhibiting the ocean circulation response, by specifying a seasonally-varying preindustrial climatology of currents, has a much larger influence on the heat storage pattern than on the carbon storage pattern. The heat storage pattern without circulation changes resembles carbon storage (either with or without circulation changes) more than it resembles the heat storage when currents are allowed to respond. This is shown to be due to the larger magnitude of the redistribution transport – the change in transport due to circulation anomalies acting on control climate gradients – for heat than for carbon. The net ocean heat and carbon uptake are slightly reduced when currents are allowed to respond. Hence, ocean circulation changes potentially act to warm the surface climate. However, the impact of the reduced carbon uptake on radiative forcing is estimated to be small while the redistribution heat transport shifts ocean heat uptake from low to high latitudes increasing its cooling power. Consequently, global surface warming is significantly reduced by circulation changes. Circulation changes also shift the pattern of warming from broad northern hemisphere amplification to a more structured pattern with reduced warming at subpolar latitudes in both hemispheres and enhanced warming near the equator.
The non-Gaussian probability distribution of sea-ice concentration makes difficulties for directly assimilating sea-ice observations into a climate model. Because of the strong impact of the atmospheric and oceanic forcing on the sea-ice state, any direct assimilation adjustment on sea-ice states is easily overridden by model physics.A new approach implements sea-ice data assimilation in enthalpy space where a sea-ice model represents a nonlinear function that transforms a positive-definite space into the sea-ice concentration subspace.Results from observation-assimilation experiments using a conceptual pycnocline prediction model that characterizes the influences of sea-ice on the decadal variability of the climate system show that the new scheme efficiently assimilates “sea-ice observations” into the model – while improving “sea-ice” variability itself, it consistently improves the estimates of all “climate” components.The resulted coupled initialization that is physically consistent among all coupled components significantly improves decadal-scale predictability of the coupled model.
We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory’s previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.
The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol-cloud interactions, chemistry-climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical-system component of earth-system models and models for decadal prediction in the near-term future, for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model.
Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud-droplet activation by aerosols, sub-grid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with eco-system dynamics and hydrology.
Most basic circulation features in AM3 are simulated as realistically, or more so, than in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks and the intensity distributions of precipitation remain problematic, as in AM2.
The last two decades of the 20th century warm in CM3 by .32°C relative to 1881-1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of .56°C and .52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol cloud interactions, and its warming by late 20th century is somewhat less realistic than in CM2.1, which warmed .66°C but did not include aerosol cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud-aerosol interactions to limit greenhouse gas warming in a way that is consistent with observed global temperature changes.
This paper documents time mean simulation characteristics from the ocean and sea ice components in a new coupled climate model developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The climate model, known as CM3, is formulated with effectively the same ocean and sea ice components as the earlier GFDL climate model, CM2.1, yet with extensive developments made to the atmosphere and land model components. Both CM2.1 and CM3 show stable mean climate indices, such as large scale circulation and sea surface temperatures (SSTs). There are notable improvements in the CM3 climate simulation relative to CM2.1, including a modified SST bias pattern and reduced biases in the Arctic sea ice cover. We anticipate SST differences between CM2.1 and CM3 in lower latitudes through analysis of the atmospheric fluxes at the ocean surface in corresponding Atmospheric Model Intercomparison Project (AMIP) simulations. In contrast, SST changes in the high latitudes are dominated by ocean and sea ice effects absent in AMIP simulations. The ocean interior simulation in CM3 is generally warmer than CM2.1, which adversely impacts the interior biases.
The sensitivity of Northern Hemisphere sea ice cover to global temperature change is examined in a group of climate models and in the satellite era observations. The models are found to have well defined, distinguishable sensitivities in climate change experiments. The satellite era observations show a larger sensitivity - a larger decline per degree warming - than any of the models. To evaluate the role of natural variability in this discrepancy, the sensitivity PDF is constructed based upon the observed trends and natural variability of multi-decadal ice cover and global temperature trends in a long control run of the GFDL CM2.1 climate model. This comparison shows that the model sensitivities range from about one to more than two pseudo-standard deviations of the variability smaller than observations indicate. The impact of natural Atlantic multi-decadal temperature trends (as simulated by the GFDL model) on the sensitivity distribution is examined and found to be minimal.
The fast and slow components of global warming in a comprehensive climate model are isolated by examining the response to an instantaneous return to pre-industrial forcing. The response is characterized by an initial fast exponential decay with an e-folding time smaller than 5 years, leaving behind a remnant that evolves more slowly. The slow component is estimated to be small at present, as measured by the global mean near-surface air temperature, and, in the model examined, grows to 0.4C by 2100 in the A1B SRES scenario and then to 1.4C by 2300 if one holds radiative forcing fixed after 2100. The dominance of the fast component at present is supported by examining the response to an instantaneous doubling of CO2 and by the excellent fit to the model's ensemble mean 20th century evolution with a simple one-box model with no long times scales.
We propose a modification to the standard forcing/feedback diagnostic energy balance model to account for 1) differences between effective and equilibrium climate sensitivities and 2) the variation of effective sensitivity over time in climate change experiments with coupled atmosphere-ocean climate models. In the spirit of Hansen et al (2005) we introduce an efficacy factor to the ocean heat uptake. Comparing the time-evolution of the surface warming in high and low efficacy models demonstrates the role of this efficacy in the transient response to CO2 forcing. Abrupt CO2 increase experiments show that the large efficacy of the Geophysical Fluid Dynamics Laboratory's CM2.1 model sets up in the first two decades following the increase in forcing. The use of an efficacy is necessary to fit this model's global mean temperature evolution in periods with both increasing and stable forcing. The inter-model correlation of transient climate response with ocean heat uptake efficacy is greater than its correlation with equilibrium climate sensitivity in an ensemble of climate models used for the 3rd and 4th IPCC assessments. When computed at the time of doubling in the standard experiment with 1%/yr increase in CO2, the efficacy is variable amongst the models but is generally greater than 1, averages between 1.3 and 1.4, and is as large as 1.75 in several models.
Coordinated Ocean-ice Reference Experiments (COREs) are presented as a tool to explore the behaviour of global ocean-ice models under forcing from a common atmospheric dataset. We highlight issues arising when designing coupled global ocean and sea ice experiments, such as difficulties formulating a consistent forcing methodology and experimental protocol. Particular focus is given to the hydrological forcing, the details of which are key to realizing simulations with stable meridional overturning circulations.
The atmospheric forcing from [Large, W., Yeager, S., 2004. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies. NCAR Technical Note: NCAR/TN-460+STR. CGD Division of the National Center for Atmospheric Research] was developed for coupled-ocean and sea ice models. We found it to be suitable for our purposes, even though its evaluation originally focussed more on the ocean than on the sea-ice. Simulations with this atmospheric forcing are presented from seven global ocean-ice models using the CORE-I design (repeating annual cycle of atmospheric forcing for 500 years). These simulations test the hypothesis that global ocean-ice models run under the same atmospheric state produce qualitatively similar simulations. The validity of this hypothesis is shown to depend on the chosen diagnostic. The CORE simulations provide feedback to the fidelity of the atmospheric forcing and model configuration, with identification of biases promoting avenues for forcing dataset and/or model development.
Winton, Michael, December 2008: Sea ice - Albedo feedback and nonlinear Arctic climate change In Arctic Sea Ice Decline, Washington, DC, American Geophysical Union, 111-131. Abstract PDF
The potential for sea ice-albedo feedback to give rise to nonlinear climate change in the Arctic Ocean defined as a nonlinear relationship between polar and global temperature change or, equivalently, a time-varying polar amplification is explored in IPCC AR4 climate models. Five models supplying SRES A1B ensembles for the 21st century are examined and very linear relationships are found between polar and global temperatures (indicating linear Arctic Ocean climate change), and between polar temperature and albedo (the potential source of nonlinearity). Two of the climate models have Arctic Ocean simulations that become annually sea ice-free under the stronger CO2 increase to quadrupling forcing. Both of these runs show increases in polar amplification at polar temperatures above -5°C and one exhibits heat budget changes that are consisten with the small ice cap instability of simple energy balance models. Both models show linear warming up to a polar temperature of -5° well above the disappearance of their September ice covers at about -9°C. Below -5°C, surface albedo decreases smoothly as reductions move, progressively, to earlier parts of the sunlit period. Atmospheric heat transport exerts a strong cooling effect during the transition to annually ice-free conditions. Specialized experiments with atmosphere and coupled models show that the main damping mechanism for sea ice region surface temperature is reduced upward heat flux through the adjacent ice-free oceans resulting in reduced atmospheric heat transport into the region.
We
present
a
mechanism
for
exchange
of
quantities
between
components
of
a
coupled
Earth
system
model,
where
each
component
is
independently
discretized.
The
exchange
grid
is
formed
by
overlaying
two
grids,
such
that
each
exchange
grid
cell
has
a
unique
parent
cell
on
each
of
its
antecedent
grids.
In
Earth
System
models
in
particular,
processes
occurring
near
component
surfaces
require
special
surface
boundary
layer
physical
processes
to
be
represented
on
the
exchange
grid.
The
exchange
grid
is
thus
more
than
just
a
stage
in
a
sequence
of
regrid-
ding
between
component
grids.
We
present
the
design
and
use
of
a
2-dimensional
exchange
grid
on
a
horizontal
planetary
surface
in
the
GFDL
Flexible
Modeling
System
(FMS),
highlighting
issues
of
parallelism
and
performance
The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved.
Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2° latitude × 2.5° longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1° in latitude and longitude, with meridional resolution equatorward of 30° becoming progressively finer, such that the meridional resolution is 1/3° at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments.
The control simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic.
Both models have been used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO2. The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/).
Manuscript received 8 December 2004, in final form 18 March 2005
The current generation of coupled climate models run at the Geophysical Fluid Dynamics Laboratory (GFDL) as part of the Climate Change Science Program contains ocean components that differ in almost every respect from those contained in previous generations of GFDL climate models. This paper summarizes the new physical features of the models and examines the simulations that they produce. Of the two new coupled climate model versions 2.1 (CM2.1) and 2.0 (CM2.0), the CM2.1 model represents a major improvement over CM2.0 in most of the major oceanic features examined, with strikingly lower drifts in hydrographic fields such as temperature and salinity, more realistic ventilation of the deep ocean, and currents that are closer to their observed values. Regional analysis of the differences between the models highlights the importance of wind stress in determining the circulation, particularly in the Southern Ocean. At present, major errors in both models are associated with Northern Hemisphere Mode Waters and outflows from overflows, particularly the Mediterranean Sea and Red Sea.
The climate response to idealized changes in the atmospheric CO2 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 CO2 concentration. When forced by a 1% per year CO2 increase, the surface air temperature difference around the time of CO2 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 CO2. By the time of CO2 doubling, the weakening in CM2.1 is larger than that found in CM2.0: 7 and 4 Sv (1 Sv 106 m3 s−1), 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).
A group of twelve IPCC fourth assessment report (AR4) climate models have Arctic (60N–90N) warmings that are, on average, 1.9 times greater than their global warmings at the time of CO2 doubling in 1%/year CO2 increase experiments. Forcings and feedbacks that impact the warming response are estimated for both Arctic and global regions based on standard model diagnostics. Fitting a zero-dimensional energy balance model to each region, an expression is derived that gives the Arctic amplification as a function of these forcings and feedbacks. Contributing to Arctic amplification are the Arctic-global differences in surface albedo feedback (SAF), longwave feedback and the net top-of-atmosphere flux forcing (the sum of the surface flux and the atmospheric heat transport convergence). The doubled CO2 forcing and non-SAF shortwave feedback oppose Arctic amplification. SAF is shown to be a contributing, but not a dominating, factor in the simulated Arctic amplification and its intermodel variation.
Two IPCC fourth assessment report climate models have Arctic Ocean simulations that become sea-ice-free year around in 1%/year CO2 increase to quadrupling experiments. These runs are examined for evidence of accelerated climate change associated with the removal of sea ice, particularly due to increasing surface albedo feedback. Both models become seasonally ice-free at an annual mean polar temperature of −9°C without registering much impact on the surface albedo feedback or disturbing the linear relationship between Arctic Ocean climate change and that of the surrounding region. When the polar temperature rises above −5°C, however, there is a sharp increase in the surface albedo feedback of one of the models, driving an abrupt elimination of Arctic ice and an increase in temperature above that expected from warming of the surrounding region. The transition to ice-free conditions is more linear in the other model, with ocean heat flux playing the primary driving role.
Winton, Michael, 2006: Surface Albedo Feedback Estimates for the AR4 Climate Models. Journal of Climate, 19(3), 359-365. Abstract PDF
A technique for estimating surface albedo feedback (SAF) from standard monthly mean climate model diagnostics is applied to the 1% yr−1 carbon dioxide (CO2)-increase transient climate change integrations of 12 Intergovernmental Panel on Climate Change (IPCC) fourth assessment report (AR4) climate models. Over the 80-yr runs, the models produce a mean SAF at the surface of 0.3 W m−2 K−1 with a standard deviation of 0.09 W m−2 K−1. Relative to 2 × CO2 equilibrium run estimates from an earlier group of models, both the mean SAF and the standard deviation are reduced. Three-quarters of the model mean SAF comes from the Northern Hemisphere in roughly equal parts from the land and ocean areas. The remainder is due to Southern Hemisphere ocean areas. The SAF differences between the models are shown to stem mainly from the sensitivity of the surface albedo to surface temperature rather from the impact of a given surface albedo change on the shortwave budget.
This paper summarizes the formulation of the ocean component to the Geophysical Fluid Dynamics Laboratory's (GFDL) climate model used for the 4th IPCC Assessment (AR4) of global climate change. In particular, it reviews the numerical schemes and physical parameterizations that make up an ocean climate model and how these schemes are pieced together for use in a state-of-the-art climate model. Features of the model described here include the following: (1) tripolar grid to resolve the Arctic Ocean without polar filtering, (2) partial bottom step representation of topography to better represent topographically influenced advective and wave processes, (3) more accurate equation of state, (4) three-dimensional flux limited tracer advection to reduce overshoots and undershoots, (5) incorporation of regional climatological variability in shortwave penetration, (6) neutral physics parameterization for representation of the pathways of tracer transport, (7) staggered time stepping for tracer conservation and numerical efficiency, (8) anisotropic horizontal viscosities for representation of equatorial currents, (9) parameterization of exchange with marginal seas, (10) incorporation of a free surface that accommodates a dynamic ice model and wave propagation, (11) transport of water across the ocean free surface to eliminate unphysical "virtual tracer flux" methods, (12) parameterization of tidal mixing on continental shelves. We also present preliminary analyses of two particularly important sensitivities isolated during the development process, namely the details of how parameterized subgridscale eddies transport momentum and tracers.
Observational and modelling evidence suggest that poleward ocean heat transport (OHT) can vary in response to both natural climate variability and greenhouse warming. Recent modelling studies have shown that increased OHT warms both the tropical and global mean climates. Using two different coupled climate models with mixed-layer oceans, with and without OHT, along with a coupled model with a fixed-current ocean component in which the currents are uniformly reduced and increased by 50%, an attempt is made to explain why this may happen.
OHT warms the global mean climate by 1 to 1.6 K in the atmospheric general circulation (AGCM) ML model and 3.5 K in the AGCM fixed current model. In each model the warming is attributed to an increase in atmospheric greenhouse trapping, primarily clear-sky greenhouse trapping, and a reduction in albedo. This occurs as OHT moistens the atmosphere, particularly at subtropical latitudes. This is not purely a thermodynamic response to the reduction in planetary albedo at these latitudes. It is a change in atmospheric circulation that both redistributes the water vapour and allows for a global atmospheric moistening—a positive 'dynamical' water vapour feedback. With increasing OHT the atmospheric water vapour content increases as atmospheric convection spreads out of the deep tropics. The global mean planetary albedo is decreased with increased OHT. This change is explained by a decrease in subtropical and mid-latitude low cloudiness, along with a reduction in high-latitude surface albedo due to decreased sea ice. The climate models with the mixed layer oceans underestimate both the subtropical low cloud cover and the high-latitude sea ice/surface albedo, and consequently have a smaller warming response to OHT.
Winton, Michael, 2005: Simple optical models for diagnosing surface-atmosphere shortwave interactions. Journal of Climate, 18(18), 3796-3805. Abstract PDF
A technique is developed for diagnosing effective surface and atmospheric optical properties from climate model shortwave flux diagnostics. These properties can be used to distinguish the contributions of surface and atmospheric optical property changes to shortwave flux changes at the surface and top of the atmosphere. In addition to the four standard shortwave flux diagnostics (upward, downward, surface, and top of atmosphere), the technique makes use of surface-down and top-up fluxes over a zero-albedo surface obtained from an auxiliary online shortwave calculation. The simple model optical properties, when constructed from the time-mean fluxes, are effective optical properties, useful for predicting the time-mean response to optical property changes. The technique is tested against auxiliary online shortwave calculations at four validation albedos and shown to predict the monthly mean surface absorption with an rms error of less than 2% over the globe. The reasons for the accuracy of the technique are explored. Less accurate techniques that make use of existing shortwave diagnostics are presented and compared.
Winton, Michael, 2003: On the climatic impact of ocean circulation. Journal of Climate, 16(17), 2875-2889. Abstract PDF
Integrations of coupled climate models with mixed-layer and fixed-current ocean components are used to explore the climatic response to varying magnitudes of ocean circulation. Four mixed-layer ocean experiments without ocean heat transports are performed using two different atmosphere-land components--the new GFDL AM2 and the GFDL Manabe Climate Model (MCM)--and two different sea ice components, one dynamic and one thermodynamic. Both experiments employing the dynamic sea ice component develop unstable growth of sea ice while the experiments with a thermodynamic sea ice component develop very large but stable ice covers. The global cooling ranges from modest to extreme in the four experiments.
Using the fixed-current climate model, a trio of 100-yr integrations are made with control currents from a GFDL R30 ocean simulation, same currents reduced by 50%, and same currents increased by 50%. This suite is performed with two coupled models again employing the two atmosphere-land components, AM2 and MCM, for a total of six experiments. Both models show a large sensitivity of the sea ice extent to the magnitude of currents with increased currents reducing the extent and warming the high latitudes. Low cloud cover also responds to circulation changes in both models but in the opposite sense. In the AM2-based model, low cloudiness decreases as ocean circulation increases, reinforcing the sea ice changes in reducing the planetary reflectivity, and warming the climate. This cloudiness change is associated with a reduction in lower-atmospheric stability over the ocean. Because the AM2-based model is able to simulate the observed seasonal low cloud-stability relationship and the changes in these quantities with altered ocean circulation are consistent with this relationship, the AM2 interpretation of the cloud changes is favored.
Cubasch, U, Gerald A Meehl, G J Boer, Ronald J Stouffer, Martin R Dix, A Noda, Catherine A Senior, S C B Raper, K S Yap, A Abe-Ouchi, S Brinkop, M Claussen, Matthew Collins, J Evans, I Fischer-Bruns, John C Fyfe, A Ganopolski, Jonathan M Gregory, Zeng-Zhen Hu, Fortunat Joos, Thomas R Knutson, Reto Knutti, Christopher Landsea, L Mearns, P C D Milly, J F B Mitchell, T Nozawa, H Paeth, J Räisänen, R Sausen, Steven J Smith, T F Stocker, Axel Timmermann, U Ulbrich, A J Weaver, J Wegner, P Whetton, T M L Wigley, Michael Winton, and F Zwiers, 2001: Projections of future climate change In Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge, UK, Cambridge University Press, 526-582.
Gregory, Jonathan M., J A Church, G J Boer, Keith W Dixon, G M Flato, D R Jackett, J A Lowe, S P O'Farrell, M M Rienecker, G Russell, Ronald J Stouffer, and Michael Winton, 2001: Comparison of results from several AOGCMs for global and regional sea-level change 1900-2100. Climate Dynamics, 18(3/4), 225-240. Abstract PDF
Sea-level rise is an important aspect of climate change because of its impact on society and ecosystems. Here we present an intercomparison of results from ten coupled atmosphere-ocean general circulation models (AOGCMs) for sea-level changes simulated for the twentieth century and projected to occur during the twenty first century in experiments following scenario 1892a for greenhouse gases and sulphate aerosols. The model results suggest that the rate of sea-level rise due to thermal expansion of sea water has increased during the twentieth century, but the small set of tide gauges with long records might not be adequate to detect this acceleration. The rate of sea-level rise due to thermal expansion continues to increase throughout the twenty first century, and the projected total is consequently larger than in the twentieth century; for 1990-2090 it amounts to 0.20-0.37 m. This wide range results from systematic uncertainty in modeling of climate change and of heat uptake by the ocean. The AOGCMs agree that sea-level rise is expected to be geographically non-uniform, with some regions experiencing as much as twice the global average, and others practically zero, but they do not agree about the geographical pattern. The lack of agreement indicates that we cannot currently have confidence in projections of local sea-level changes, and reveals a need for detailed analysis and intercomparison in order to understand and reduce the disagreements.
Winton, Michael, 2000: A reformulated three-layer sea ice model. Journal of Atmospheric and Oceanic Technology, 17(4), 525-531. Abstract PDF
A model is presented that provides an efficient approximation to sea ice thermodynamics for climate studies. Semtner's three-layer framework is used, but the brine content of the upper ice is represented with a variable heat capacity as is done in more physically based models. A noniterative fully implicit time-stepping scheme is used for calculation of ice temperature. The results of the new model are compared to those of Semtner's original model.
Winton, Michael, 1999: Polar water column stability. Journal of Physical Oceanography, 29(6), 1368-1371. Abstract PDF
An expression is derived for the surface salt input needed to induce complete convective overturning of a polar water column consisting of 1) a layer of sea ice, 2) a freezing temperature mixed layer, 3) a pycnocline with linearly varying temperature and salinity, and 4) deep water with fixed temperature and salinity. This quantity has been termed the bulk stability by Martinson. The bulk stability is found to consist of three components. The first two make up Martinson's salt deficit and are the salt input needed to increase the density of the mixed layer and the pycnocline layer to that of the deep water (the mixed layer stability and pycnocline layer stability, respectively). The third component is Martinson's thermal barrier: the potential for pycnocline heat to melt ice, reducing the surface salinity. It is found that when the pycnocline density gradient due to temperature offsets more than one half of that due to salinity, the pycnocline layer stability is negative. Consequently, it is possible for a stably stratified water column to have zero or negative bulk stability.
Martinson, D G., David S Battisti, R Bradley, J Cole, R Fine, M Ghil, Y Kushnir, Syukuro Manabe, M S McCartney, P McCormick, Michael J Prather, E Sarachik, P P Tans, L Thompson, and Michael Winton, 1998: Decade-to-Century-Scale Climate Variability and Climate Change: A Science Strategy, Washington, D.C.: National Research Council, 42 pp.
An important component of the ocean's thermohaline circulation is the sinking of dense water from continental shelves to abyssal depths. Such downslope flow is thought to be a consequence of bottom stress retarding the alongslope flow of density-driven plumes. In this paper the authors explore the potential for explicitly simulating the simple mechanism in z-coordinate models. A series of experiments are performed using a twin density-coordinate model simulation as a standard of comparison. The adiabatic nature of the experiments and the importance of bottom slope make it more likely that the density-coordinate model will faithfully reproduce the solution. The difficulty of maintaining the density signal as the plume descends the slope is found to be the main impediment to accurate simulation in the z-coordinate model. The results of process experiments suggest that the model solutions will converge when the z-coordinate model has sufficient vertical resolution to resolve the bottom viscous layer and horizontal grid spacing equal to its vertical grid spacing divided by the maximum slope. When this criterion is met it is shown that the z-coordinate model converges to an analytical solution for a simple two-dimensional flow.
Winton, Michael, 1997: The damping effect of bottom topography on internal decadal-scale oscillations of the thermohaline circulation. Journal of Physical Oceanography, 27(1), 203-208. Abstract PDF
By comparing the response of flat and bowl-shaped basins to fixed heat fluxes of various magnitudes, it is determined that coastal topography has a considerable damping influence upon internal decadal oscillations of the thermohaline circulation. It is proposed that this is because the adjustment of baroclinic currents to the no-normal-flow boundary condition at weakly stratified coasts is aided in the topography case by the generation of substantial barotropic flow.
Winton, Michael, 1997: The effect of cold climate upon North Atlantic Deep Water formation in a simple ocean-atmosphere model. Journal of Climate, 10(1), 37-51. Abstract PDF
The sensitivity of North Atlantic Deep Water formation to variations in mean surfae temperature is explored with a meridional-vertical plane ocean model coupled to an energy balance atmosphere. It is found that North Atlantic Deep Water formation is favored by a warm climate, while cold climates are more likely to produce Southern Ocean deep water or deep-decoupling oscillations (when the Southern sinking region is halocline covered). This behavior is traced to a cooling-induced convective instability near the North Atlantic sinking region, that is, to unstable horizontal spreading of a halocline that stratifies part of the region. Under the convective instability it is found that climate cooling is generally equivalent to increased freshwater forcing. This is because in a cold climate, high-latitude water masses approach the temperature of maximum density and the convection-driving, upward thermal buoyancy flux induced by surface cooling becomes insufficient to overcome the stratifying effect of surface freshening (a downward buoyancy flux). An extensive halocline is then formed and this halocline interferes with the heat loss necessary for the steady production of North Atlantic Deep Water.
Winton, Michael, 1996: The role of horizontal boundaries in parameter sensitivity and decadal-scale variability of coarse-resolution ocean general circulation models. Journal of Physical Oceanography, 26(3), 289-304. Abstract PDF
Coarse-resolution f-plane and beta-plane frictional geostrophic models are used to study the response to restored surface buoyancies and fixed surface buoyancy fluxes. With restored surface buoyancies it is found that the overturning and meridional buoyancy transport generally follow the scaling from thermal wind and vertical advective-diffusive balance. When maximum midbasin rather than maximum overall streamfunction is used as the metric, the sensitivity of overturning magnitude to vertical diffusivity agrees quite closely with that of a two-dimensional Rayleigh frictional model that follows an analogous scaling. This measure of meridional overturning, as well as the meridional buoyancy transport, also roughly follow the predicted f -1/3 scaling and are relatively insensitive to variations in beta and horizontal viscosity.
The sensitivity experiments indicate that the coarse-resolution model overturning and thermodynamic structure are well characterized by the adjustment of a uniformly rotating viscous fluid to boundaries parallel to the surface forcing gradient. This adjustment is examined in more detail with linear and nonlinear models. In stratified regions, thermal wind currents normal to the coast initiate the adjustment by forcing pycnocline depth anomalies in the coastal (horizontal) Ekman layer. These anomalies propagate around the coast in the Kelvin wave direction, setting up geostrophic currents parallel to the coast. In a model initialized with a high-latitude baroclinic jet, a warm (cold) boundary signal spreads around the poleward (equatorward) part of the basin, initiating geostrophic currents connecting the flow onto the coast in the east with that away from the coast in the west to form two gyres. The warm coastal signal propagates slowly along the weakly stratified polar wall. In steady circulations, strong damping inhibits warm signal propagation, and the warm water on the poleward part of the eastern coast is forced to downwell.
Self-sustaining decadal-scale oscillation is a robust feature of the models when forced with fixed buoyancy fluxes. This variability is inherently three-dimensional (it does not occur in a two-dimensional frictional model) and involves the periodic growth and decay of a baroclinic jet in the poleward eastern corner. Decay occurs when a jetlike disturbance, normal to the coast, propagates cyclonically around the basin replacing the cold water along the boundary with warm. Forcing of thermal wind currents normal to weakly stratified coasts and weak damping of the resulting propagating boundary disturbances are found to be conducive to oscillations.
