Abstract
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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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 sealevel 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, while the weakening of the AMOC causes redistribution of heat towards lower latitudes. 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.
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 CO₂. The output from a suite of integrations conducted with these models is freely available online (seehttp://nomads.gfdl.noaa.gov/).
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for the 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, subgrid 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 ecosystem dynamics and hydrology. Its horizontal resolution is approximately 200 km, and its vertical resolution ranges approximately from 70 m near the earth’s surface to 1 to 1.5 km near the tropopause and 3 to 4 km in much of the stratosphere. Most basic circulation features in AM3 are simulated as realistically, or more so, as in AM2. In particular, dry biases have been reduced over South America. In coupledmode, 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 remains problematic, as in AM2. The most intense 0.2% of precipitation rates occur less frequently in AM3 than observed. The last two decades of the twentieth century warm in CM3 by 0.32°C relative to 1881–1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of 0.56° and 0.52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol–cloud interactions, and its warming by the late twentieth century is somewhat less realistic than in CM2.1, which warmed 0.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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract
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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Because of its persistence on seasonal time scales, 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-yr control integration and a suite of initialized ensemble forecasts from a fully coupled global climate model. This 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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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 Geophysical Fluid Dynamics Laboratory models incorporating eddy‐resolving (0.1°), eddy‐permitting (0.25°), and eddy‐parameterizing (1°) 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 sea surface temperature (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 Geophysical Fluid Dynamics Laboratory models, but we do find substantial impact on precipitation changes. The eddy‐parameterizing model projects a banded structure common to the Coupled Model Intercomparison Project (Phase 5) 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.
Plain Language Summary
Global coupled atmosphere‐ocean models have not included the explicit simulation of ocean mesoscale eddies (turbulent ocean flow with a horizontal scale of less than 100 km) until the recent few years. As a result, current estimates of precipitation responses to greenhouse gases increase have predominantly relied on models in which ocean eddies have to be parameterized. Given the importance of eddies on ocean heat transport and the strong dependence of precipitation on ocean surface temperature, the lack of ocean eddy resolution may have influenced the projection of precipitation changes. In this study, we compare precipitation changes from CO2 forcing in models with and without explicit ocean eddy simulation. It is found that the incorporation of ocean eddy resolution can substantially influence the simulation of precipitation changes in eddy‐rich regions and the deep tropics. Some land regions are also affected, albeit to a lesser degree compared to the oceanic regions.
Key Points
The incorporation of ocean eddy resolution has substantial influences on regional precipitation changes
The largest impact of ocean eddy resolution on precipitation projections is found in eddy‐rich regions and the deep tropics
Land regions are generally not very sensitive to ocean eddy resolution
Abstract
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 CO
2
forcing has yet to be investigated. This study analyzes precipitation changes from a suite of Geophysical Fluid Dynamics Laboratory models incorporating eddy‐resolving (0.1°), eddy‐permitting (0.25°), and eddy‐parameterizing (1°) 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 sea surface temperature (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 Geophysical Fluid Dynamics Laboratory models, but we do find substantial impact on precipitation changes. The eddy‐parameterizing model projects a banded structure common to the Coupled Model Intercomparison Project (Phase 5) 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.
Plain Language Summary
Global coupled atmosphere‐ocean models have not included the explicit simulation of ocean mesoscale eddies (turbulent ocean flow with a horizontal scale of less than 100 km) until the recent few years. As a result, current estimates of precipitation responses to greenhouse gases increase have predominantly relied on models in which ocean eddies have to be parameterized. Given the importance of eddies on ocean heat transport and the strong dependence of precipitation on ocean surface temperature, the lack of ocean eddy resolution may have influenced the projection of precipitation changes. In this study, we compare precipitation changes from CO
2
forcing in models with and without explicit ocean eddy simulation. It is found that the incorporation of ocean eddy resolution can substantially influence the simulation of precipitation changes in eddy‐rich regions and the deep tropics. Some land regions are also affected, albeit to a lesser degree compared to the oceanic regions.
Key Points
The incorporation of ocean eddy resolution has substantial influences on regional precipitation changes
The largest impact of ocean eddy resolution on precipitation projections is found in eddy‐rich regions and the deep tropics
Land regions are generally not very sensitive to ocean eddy resolution
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 steady-state 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.
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, areeffectiveoptical 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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK