Seasonal to interannual climate predictability in mid and high northern latitudes in a global coupled model

The upper limit of climate predictability in mid and high northern latitudes on seasonal to interannual time scales is investigated by performing two perfect ensemble experiments with the global coupled atmosphere–ocean–sea ice model ECHAM5/MPI-OM. The ensembles consist of six members and are initialized in January and July from different years of the model’s 300-year control integration. The potential prognostic predictability is analyzed for a set of oceanic and atmospheric climate parameters. The predictability of the atmospheric circulation is small except for southeastern Europe, parts of North America and the North Pacific, where significant predictability occurs with a lead time of up to half a year. The predictability of 2 m air temperature shows a large land–sea contrast with highest predictabilities over the sub polar North Atlantic and North Pacific. A combination of relatively high persistence and advection of sea surface temperature anomalies into these areas leads to large predictability. Air temperature over Europe, parts of North America and Asia shows significant predictability of up to half a year in advance. Over the ice-covered Arctic, air temperature is not predictable at time scales exceeding 2 months. Sea ice thickness is highly predictable in the central Arctic mainly due to persistence and in the Labrador Sea due to dynamics. Surface salinity is highly predictable in the Arctic Ocean, northern North Atlantic and North Pacific for several years in advance. We compare the results to the predictability due to persistence and show the importance of dynamical processes for the predictability.     

Interannual predictability of Arctic sea ice in a global climate model: regional contrasts and temporal evolution

The predictability of the Arctic sea ice is investigated at the interannual time scale using decadal experiments performed within the framework of the fifth phase of the Coupled Model Intercomparison Project with the CNRM-CM5.1 coupled atmosphere–ocean global climate model. The predictability of summer Arctic sea ice extent is found to be weak and not to exceed 2 years. In contrast, robust prognostic potential predictability (PPP) up to several years is found for winter sea ice extent and volume. This predictability is regionally contrasted. The marginal seas in the Atlantic sector and the central Arctic show the highest potential predictability, while the marginal seas in the Pacific sector are barely predictable. The PPP is shown to decrease drastically in the more recent period. Regarding sea ice extent, this decrease is explained by a strong reduction of its natural variability in the Greenland–Iceland–Norwegian Seas due to the quasi-disappearance of the marginal ice zone in the center of the Greenland Sea. In contrast, the decrease of predictability of sea ice volume arises from the combined effect of a reduction of its natural variability and an increase in its chaotic nature. The latter is attributed to a thinning of sea ice cover over the whole Arctic, making it more sensitive to atmospheric fluctuations. In contrast to the PPP assessment, the prediction skill as measured by the anomaly correlation coefficient is found to be mostly due to external forcing. Yet, in agreement with the PPP assessment, a weak added value of the initialization is found in the Atlantic sector. Nevertheless, the trend-independent component of this skill is not statistically significant beyond the forecast range of 3 months. These contrasted findings regarding potential predictability and prediction skill arising from the initialization suggest that substantial improvements can be made in order to enhance the prediction skill.     

Sensitivity of open-water ice growth and ice concentration evolution in a coupled atmosphere-ocean-sea ice model

A coupled atmosphere-ocean-sea ice model is applied to investigate to what degree the area-thickness distribution of new ice formed in open water affects the ice and ocean properties. Two sensitivity experiments are performed which modify the horizontal-to-vertical aspect ratio of open-water ice growth. The resulting changes in the Arctic sea-ice concentration strongly affect the surface albedo, the ocean heat release to the atmosphere, and the sea-ice production. The changes are further amplified through a positive feedback mechanism among the Arctic sea ice, the Atlantic Meridional Overturning Circulation (AMOC), and the surface air temperature in the Arctic, as the Fram Strait sea ice import influences the freshwater budget in the North Atlantic Ocean. Anomalies in sea-ice transport lead to changes in sea surface properties of the North Atlantic and the strength of AMOC. For the Southern Ocean, the most pronounced change is a warming along the Antarctic Circumpolar Current (ACC), owing to the interhemispheric bipolar seasaw linked to AMOC weakening. Another insight of this study lies on the improvement of our climate model. The ocean component FESOM is a newly developed ocean-sea ice model with an unstructured mesh and multi-resolution. We find that the subpolar sea-ice boundary in the Northern Hemisphere can be improved by tuning the process of open-water ice growth, which strongly influences the sea ice concentration in the marginal ice zone, the North Atlantic circulation, salinity and Arctic sea ice volume. Since the distribution of new ice on open water relies on many uncertain parameters and the knowledge of the detailed processes is currently too crude, it is a challenge to implement the processes realistically into models. Based on our sensitivity experiments, we conclude a pronounced uncertainty related to open-water sea ice growth which could significantly affect the climate system sensitivity.     

东西伯利亚——波弗特海海冰多年代际变率及其与大西洋多年代际振荡的联系

基于哈得来中心（Hadley Centre）逐月的海表温度、海冰密集度资料以及美国国家环境预报中心/国家大气研究中心（NCEP/NCAR）的大气环流再分析资料,分析了1950—2020年秋季（8—10月）东西伯利亚—波弗特海（East Siberian-Beaufort,EsCB）海冰年代际变化的时空特征,并阐述了大西洋多年代际振荡（Atlantic Multidecadal Oscillation,AMO）对EsCB海冰年代际变率的可能调制作用。结果表明,EsCB是秋季北极海冰年代际变化最主要的区域,该区海冰密集度年代际变率可占其异常总方差的40%以上。进一步研究发现,AMO对秋季EsCB海冰存在明显的调制作用,在AMO正位相,北大西洋正海温异常激发向极传播的大气罗斯贝波列,有利于北极中部出现高压异常,相应的大气绝热下沉运动使得对流层低层出现明显的升温,从而有利于EsCB海冰的融化。与此同时,地表升温和EsCB海冰消融会引起局地云量的增多、大气向下长波辐射增大,这反过来又使得地表气温升高,这种地表气温-云-长波辐射的正反馈过程有利于年代际海冰信号的长时间维持。耦合模式的北大西洋“起搏...     

Arctic sea ice and Eurasian climate: A review

The Arctic plays a fundamental role in the climate system and has shown significant climate change in recent decades,including the Arctic warming and decline of Arctic sea-ice extent and thickness. In contrast to the Arctic warming and reduction of Arctic sea ice, Europe, East Asia and North America have experienced anomalously cold conditions, with record snowfall during recent years. In this paper, we review current understanding of the sea-ice impacts on the Eurasian climate.Paleo, observational and modelling studies are covered to summarize several major themes, including: the variability of Arctic sea ice and its controls; the likely causes and apparent impacts of the Arctic sea-ice decline during the satellite era,as well as past and projected future impacts and trends; the links and feedback mechanisms between the Arctic sea ice and the Arctic Oscillation/North Atlantic Oscillation, the recent Eurasian cooling, winter atmospheric circulation, summer precipitation in East Asia, spring snowfall over Eurasia, East Asian winter monsoon, and midlatitude extreme weather; and the remote climate response(e.g., atmospheric circulation, air temperature) to changes in Arctic sea ice. We conclude with a brief summary and suggestions for future research.     

The effect of atmospheric CO2 and ice sheet topography on LGM climate

The role of reduced atmospheric CO₂ concentration and ice sheet topography plus its associated land albedo on the LGM climate is investigated using a coupled atmosphere-ocean-sea ice climate system model. The surface cooling induced by the reduced CO₂ concentration is larger than that by the ice sheet topography plus other factors by about 30% for the surface air temperature and by about 100% for the sea surface temperature. A large inter-hemispheric asymmetry in surface cooling with a larger cooling in the Northern Hemisphere is found for both cases. This asymmetric inter-hemispheric temperature response is consistent in the ice sheet topography case with earlier studies using an atmospheric model coupled with a mixed-layer ocean representation, but contrasts with these results in the reduced CO₂ case. The incorporation of ocean dynamics presumably leads to a larger snow and sea ice feedback as a result of the reduction in northward ocean heat transport, mainly as a consequence of the decrease in the North Atlantic overturning circulation by the substantial freshening of the North Atlantic convection regions. A reversed case is found in the Southern Ocean. Overall, the reduction in atmospheric CO₂ concentration accounts for about 60% of the total LGM climate change.     

Quantifying Arctic contributions to climate predictability in a regional coupled ocean-ice-atmosphere model

The relative importance of regional processes inside the Arctic climate system and the large scale atmospheric circulation for Arctic interannual climate variability has been estimated with the help of a regional Arctic coupled ocean-ice-atmosphere model. The study focuses on sea ice and surface climate during the 1980s and 1990s. Simulations agree reasonably well with observations. Correlations between the winter North Atlantic Oscillation index and the summer Arctic sea ice thickness and summer sea ice extent are found. Spread of sea ice extent within an ensemble of model runs can be associated with a surface pressure gradient between the Nordic Seas and the Kara Sea. Trends in the sea ice thickness field are widely significant and can formally be attributed to large scale forcing outside the Arctic model domain. Concerning predictability, results indicate that the variability generated by the external forcing is more important in most regions than the internally generated variability. However, both are in the same order of magnitude. Local areas such as the Northern Greenland coast together with Fram Straits and parts of the Greenland Sea show a strong importance of internally generated variability, which is associated with wind direction variability due to interaction with atmospheric dynamics on the Greenland ice sheet. High predictability of sea ice extent is supported by north-easterly winds from the Arctic Ocean to Scandinavia.     

Climate predictability on interannual to decadal time scales: the initial value problem

Any initial value forecast of climate will be subject to errors originating from poorly known initial conditions, model imperfections, and by "chaos" in the sense that, even if the initial conditions were perfectly known, infinitesimal errors can amplify and spoil the forecast at some lead time. Here the latter source of error is examined using a "perfect model" approach whereby small perturbations are made to a coupled atmosphere-ocean general circulation model and the spread of nearby model trajectories, on time and space scales appropriate to seasonal-decadal climate variability, is measured to assess the lead time at which the error saturates. The study therefore represents an estimate of the upper limit of the predictability of climate (appropriate to the initial value problem) given a perfect model and near perfect knowledge of the initial conditions. It is found that, on average, surface air temperature anomalies are potentially predictable on seasonal to interannual time scales in the tropical regions and are potentially predictable on decadal time scales over the ocean in the North Atlantic. For mid-latitude surface air temperature anomalies over land, model trajectories rapidly diverge and there is little sign of any potential predictability on time scales greater than a season or so. For mean sea level pressure anomalies, there is potential predictability on seasonal time scales in the tropics, and for some global scale annual-decadal anomalies, although not those associated with the North Atlantic Oscillation. For precipitation, the only potential for predictability is for seasonal time anomalies associated with the El-Niño Southern Oscillation. For the majority of the highly populated regions of the world, climate predictability on interannual to decadal time scales based in the initial value approach is likely to be severely limited by chaotic error growth. It is found however that there can be cases in which the potential predictability can be higher than average indicating that there is perhaps some utility in making initial value forecasts of climate in those regions which show low predictability on average.     

A predictability study of simulated North Atlantic multidecadal variability

 The North Atlantic is one of the few places on the globe where the atmosphere is linked to the deep ocean through air–sea interaction. While the internal variability of the atmosphere by itself is only predictable over a period of one to two weeks, climate variations are potentially predictable for much longer periods of months or even years because of coupling with the ocean. This work presents details from the first study to quantify the predictability for simulated multidecadal climate variability over the North Atlantic. The model used for this purpose is the GFDL coupled ocean-atmosphere climate model used extensively for studies of global warming and natural climate variability. This model contains fluctuations of the North Atlantic and high-latitude oceanic circulation with variability concentrated in the 40–60 year range. Oceanic predictability is quantified through analysis of the time-dependent behavior of large-scale empirical orthogonal function (EOF) patterns for the meridional stream function, dynamic topography, 170 m temperature, surface temperature and surface salinity. The results indicate that predictability in the North Atlantic depends on three main physical mechanisms. The first involves the oceanic deep convection in the subpolar region which acts to integrate atmospheric fluctuations, thus providing for a red noise oceanic response as elaborated by Hasselmann. The second involves the large-scale dynamics of the thermohaline circulation, which can cause the oceanic variations to have an oscillatory character on the multidecadal time scale. The third involves nonlocal effects on the North Atlantic arising from periodic anomalous fresh water transport advecting southward from the polar regions in the East Greenland Current. When the multidecadal oscillatory variations of the thermohaline circulation are active, the first and second EOF patterns for the North Atlantic dynamic topography have predictability time scales on the order of 10–20 y, whereas EOF-1 of SST has predictability time scales of 5–7 y. When the thermohaline variability has weak multidecadal power, the Hasselmann mechanism is dominant and the predictability is reduced by at least a factor of two. When the third mechanism is in an extreme phase, the North Atlantic dynamic topography patterns realize a 10–20 year predictability time scale. Additional analysis of SST in the Greenland Sea, in a region associated with the southward propagating fresh water anomalies, indicates the potential for decadal scale predictability for this high latitude region as well. The model calculations also allow insight into regional variations of predictability, which might be useful information for the design of a monitoring system for the North Atlantic. Predictability appears to break down most rapidly in regions of active convection in the high-latitude regions of the North Atlantic. Received: 28 October 1996 / Accepted: 21 March 1997     

Sea-ice effects on climate model sensitivity and low frequency variability

A change in a sea-ice parameter in a global coupled climate model results in a reduction in amplitude (of about 60%) and a shortening of the predominant period of decadal low frequency variability in the time series of globally averaged surface air temperature. These changes are global in extent and also are reflected in time series of area-averaged SSTs in the equatorial eastern Pacific Ocean, the principal components of the first EOFs of global surface air temperature and sea level pressure, Asian monsoon precipitations and other quantities. Coupled ocean-atmosphere-sea ice processes acting on a global scale are modified to produce these changes. Global climate sensitivity is reduced when ice albedo feedback is weakened due to the change in sea ice that makes it more difficult to melt. The changes in the amplitude and time scale of the low frequency variability in the model are traced to changes in the base state of the climate simulations as affected by modifications associated with the changes in sea ice. Making sea ice more difficult to melt results in increased sea-ice area, colder high latitudes, increased meridional surface temperature gradients, and, to a first order, stronger surface winds in most regions which strengthen near-surface currents, particularly in the Northern Hemisphere, and decreases the advection time scale in the upper ocean gyres. Additionally, in the North Atlantic there is enhanced meridional overturning due to increased density mainly in the Greenland Sea region. This also contributes to an intensified North Atlantic gyre. The changes in base state due to the sea ice change result in a more predominant decadal time scale of near 14 years and significantly reduced contributions from lower frequencies in the range of 15–40 year periods. Received: 11 December 1998 / Accepted: 4 October 1999     

Linking Arctic sea‐ice and atmospheric circulation anomalies on interannual and decadal timescales

Abstract

The relationship between Arctic sea‐ice concentration anomalies, particularly those associated with the “Great Salinity Anomaly” of 1968–1982, and atmospheric circulation anomalies north of 45°N is investigated. Empirical orthogonal function (EOF) analyses are performed on winter Arctic ice concentration from 1954 to 1990, sea level pressure and 500‐hPa heights from 1947 to 1994, and 850‐hPa temperatures from 1963 to 1994. Variability on both interannual and decadal timescales is apparent in the time series of the leading winter EOFs of all variables. The first EOF of winter sea‐ice concentration was found to characterize the patterns of ice variability associated with the Great Salinity Anomaly in the northern North Atlantic from 1968–82. Spatial maps of temporal correlation coefficients between the time series of the first EOF of winter sea‐ice concentration and the winter atmospheric anomaly fields are calculated at lags of 0 and ±7 year. Maximum correlations were found to exist when the time‐series of this ice EOF 1 leads the atmospheric anomaly fields by one year. A particularly interesting result is the connection between the presence of ice anomalies in the Greenland and Barents Seas and subsequent pressure anomalies of the same sign over the Irminger Basin and the Canadian Arctic. The main emphasis of the paper is to identify connections between Arctic sea‐ice and atmospheric circulation anomalies at interannual time‐scales.     

Arctic sea ice and the Madden–Julian Oscillation (MJO)

Arctic sea ice responds to atmospheric forcing in primarily a top-down manner, whereby near-surface air circulation and temperature govern motion, formation, melting, and accretion. As a result, concentrations of sea ice vary with phases of many of the major modes of atmospheric variability, including the North Atlantic Oscillation, the Arctic Oscillation, and the El Niño-Southern Oscillation. However, until this present study, variability of sea ice by phase of the leading mode of atmospheric intraseasonal variability, the Madden–Julian Oscillation (MJO), which has been found to modify Arctic circulation and temperature, remained largely unstudied. Anomalies in daily change in sea ice concentration were isolated for all phases of the real-time multivariate MJO index during both summer (May–July) and winter (November–January) months. The three principal findings of the current study were as follows. (1) The MJO projects onto the Arctic atmosphere, as evidenced by statistically significant wavy patterns and consistent anomaly sign changes in composites of surface and mid-tropospheric atmospheric fields. (2) The MJO modulates Arctic sea ice in both summer and winter seasons, with the region of greatest variability shifting with the migration of the ice margin poleward (equatorward) during the summer (winter) period. Active regions of coherent ice concentration variability were identified in the Atlantic sector on days when the MJO was in phases 4 and 7 and the Pacific sector on days when the MJO was in phases 2 and 6, all supported by corresponding anomalies in surface wind and temperature. During July, similar variability in sea ice concentration was found in the North Atlantic sector during MJO phases 2 and 6 and Siberian sector during MJO phases 1 and 5, also supported by corresponding anomalies in surface wind. (3) The MJO modulates Arctic sea ice regionally, often resulting in dipole-shaped patterns of variability between anomaly centers. These results provide an important first look at intraseasonal variability of sea ice in the Arctic.     

BCC-CSM2-MR模式对北极海冰和气候的模拟及预估

基于第六次耦合模式比较计划（CMIP6）,使用新一代全球模式BCC-CSM2-MR的历史试验和未来共享社会经济路径（SSPs）数据,依据Hadley中心的海表面温度和海冰密集度数据及NCEP/NCAR I再分析资料,评估了BCC-CSM2-MR模式对北极海冰及北极气候的模拟能力,并对未来变化进行了预估。结果表明：BCC-CSM2-MR模式可以较好再现北极海冰密集度、近地层大气平均温度和海表温度的多年平均空间分布特征。但模式对北极局地大气平均温度模拟存在一定偏差,可能在一定程度上导致相应地区海冰的模拟存在差异。21世纪,北极海冰范围持续减少,9月减少趋势显著,3月减少趋势相对较弱。3月北极大部地区表现为一致的增温,仅在北大西洋局部出现一定程度的降温,9月北极大气增温幅度弱于3月。与地表平均温度不同,3月和9月的北极大部地区海表温度均出现增加,且9月海表温度的增幅大于3月,仅拉布拉多海海温出现下降。     

Time-dependent response of a zonally averaged ocean–atmosphere–sea ice model to Milankovitch forcing

An ocean–atmosphere–sea ice model is developed to explore the time-dependent response of climate to Milankovitch forcing for the time interval 5–3 Myr BP. The ocean component is a zonally averaged model of the circulation in five basins (Arctic, Atlantic, Indian, Pacific, and Southern Oceans). The atmospheric component is a one-dimensional (latitudinal) energy balance model, and the sea-ice component is a thermodynamic model. Two numerical experiments are conducted. The first experiment does not include sea ice and the Arctic Ocean; the second experiment does. Results from the two experiments are used to investigate (1) the response of annual mean surface air and ocean temperatures to Milankovitch forcing, and (2) the role of sea ice in this response. In both experiments, the response of air temperature is dominated by obliquity cycles at most latitudes. On the other hand, the response of ocean temperature varies with latitude and depth. Deep water formed between 45°N and 65°N in the Atlantic Ocean mainly responds to precession. In contrast, deep water formed south of 60°S responds to obliquity when sea ice is not included. Sea ice acts as a time-integrator of summer insolation changes such that annual mean sea-ice conditions mainly respond to obliquity. Thus, in the presence of sea ice, air temperature changes over the sea ice are amplified, and temperature changes in deep water of southern origin are suppressed since water below sea ice is kept near the freezing point.     

Inherent sea ice predictability in the rapidly changing Arctic environment of the Community Climate System Model, version 3

Seasonal predictions of Arctic sea ice have typically been based on statistical regression models or on results from ensemble ice model forecasts driven by historical atmospheric forcing. However, in the rapidly changing Arctic environment, the predictability characteristics of summer ice cover could undergo important transformations. Here global coupled climate model simulations are used to assess the inherent predictability of Arctic sea ice conditions on seasonal to interannual timescales within the Community Climate System Model, version 3. The role of preconditioning of the ice cover versus intrinsic variations in determining sea ice conditions is examined using ensemble experiments initialized in January with identical ice?Cocean?Cterrestrial conditions. Assessing the divergence among the ensemble members reveals that sea ice area exhibits potential predictability during the first summer and for winter conditions after a year. The ice area exhibits little potential predictability during the spring transition season. Comparing experiments initialized with different mean ice conditions indicates that ice area in a thicker sea ice regime generally exhibits higher potential predictability for a longer period of time. In a thinner sea ice regime, winter ice conditions provide little ice area predictive capability after approximately 1?year. In all regimes, ice thickness has high potential predictability for at least 2?years.     

Variability of Fram Strait sea ice export: causes, impacts and feedbacks in a coupled climate model

Analyses of a 500-year control integration of the global coupled atmosphere–sea ice–ocean model ECHAM5.0/MPI-OM show a high variability in the ice export through Fram Strait on interannual to decadal timescales. This variability is mainly determined by variations in the sea level pressure gradient across Fram Strait and thus geostrophic wind stress. Ice thickness anomalies, formed at the Siberian coast and in the Chukchi Sea, propagate across the Arctic to Fram Strait and contribute to the variability of the ice export on a timescale of about 9 years. Large anomalies of the ice export through Fram Strait cause fresh water signals, which reach the Labrador Sea after 1–2 years and lead to significant changes in the deep convection. The associated anomalies in ice cover and ocean heat release have a significant impact on air temperature in the Labrador Sea and on the large-scale atmospheric circulation. This affects the sea ice transport and distribution in the Arctic again. Sensitivity studies, simulating the effect of large ice exports through Fram Strait, show that the isolated effect of a prescribed ice/fresh water anomaly is very important for the climate variability in the Labrador Sea. Thus, the ice export through Fram Strait can be used for predictability of Labrador Sea climate up to 2 years in advance.     

A possible feedback mechanism involving the Arctic freshwater, the Arctic sea ice, and the North Atlantic Drift

Model studies point to enhanced warming and to increased freshwater fluxes to high northern latitudes in response to global warming. In order to address possible feedbacks in the ice-ocean system in response to such changes, the combined effect of increased freshwater input to the Arctic Ocean and Arctic warming--the latter manifested as a gradual melting of the Arctic sea ice--is examined using a 3-D isopycnic coordinate ocean general circulation model. A suite of three idealized experiments is carried out: one control integration, one integration with a doubling of the modern Arctic river runoff, and a third more extreme case, where the river runoff is five times the modern value. In the two freshwater cases, the sea ice thickness is reduced by 1.5-2 m in the central Arctic Ocean over a 50-year period. The modelled ocean response is qualitatively the same for both perturbation experiments: freshwater propagates into the Atlantic Ocean and the Nordic Seas, leading to an initial weakening of the North Atlantic Drift.Furthermore, changes in the geostrophic currents in the central Arctic and melting of the Arctic sea ice lead to an intensified Beaufort Gyre, which in turn increases the southward volume transport through the Canadian Archipelago. To compensate for this southward transport of mass, more warm and saline Atlantic water is carried northward with the North Atlantic Drift. It is found that the increased transport of salt into the northern North Atlantic and the Nordic Seas tends to counteract the impact of the increased freshwater originating from the Arctic, leading to a stabilization of the North Atlantic Drift.     

Is There a Link between Arctic Sea Ice Loss and Increasing Frequency of Extremely Cold Winters in Eurasia and North America? Synthesis of Current Research

Studies dealing with impact of the Arctic warming and related sea ice decline on the Northern Hemisphere atmospheric circulation are considered. The causes of occurrence of extremely cold winters over the mid-latitude continents observed in the recent decades against the warming background are discussed. Several conceptions are outlined which explain potential reasons for occurrence of this phenomenon. The paper discusses impacts of the Arctic sea ice loss on the large-scale atmospheric circulation, oscillations of planetary waves. It also discusses issues related to sea ice changes in the Barents and Kara seas and their link to the frequency of extremely cold winters observed in Eurasia and North America, the contribution of internal atmospheric variability to the increasing frequency of cold weather, and the role of the Atlantic Multidecadal Oscillation in the Arctic sea ice reduction.     

Centennial-scale climate change from decadally-paced explosive volcanism: a coupled sea ice-ocean mechanism

Northern Hemisphere summer cooling through the Holocene is largely driven by the steady decrease in summer insolation tied to the precession of the equinoxes. However, centennial-scale climate departures, such as the Little Ice Age, must be caused by other forcings, most likely explosive volcanism and changes in solar irradiance. Stratospheric volcanic aerosols have the stronger forcing, but their short residence time likely precludes a lasting climate impact from a single eruption. Decadally paced explosive volcanism may produce a greater climate impact because the long response time of ocean surface waters allows for a cumulative decrease in sea-surface temperatures that exceeds that of any single eruption. Here we use a global climate model to evaluate the potential long-term climate impacts from four decadally paced large tropical eruptions. Direct forcing results in a rapid expansion of Arctic Ocean sea ice that persists throughout the eruption period. The expanded sea ice increases the flux of sea ice exported to the northern North Atlantic long enough that it reduces the convective warming of surface waters in the subpolar North Atlantic. In two of our four simulations the cooler surface waters being advected into the Arctic Ocean reduced the rate of basal sea-ice melt in the Atlantic sector of the Arctic Ocean, allowing sea ice to remain in an expanded state for?>?100 model years after volcanic aerosols were removed from the stratosphere. In these simulations the coupled sea ice-ocean mechanism maintains the strong positive feedbacks of an expanded Arctic Ocean sea ice cover, allowing the initial cooling related to the direct effect of volcanic aerosols to be perpetuated, potentially resulting in a centennial-scale or longer change of state in Arctic climate. The fact that the sea ice-ocean mechanism was not established in two of our four simulations suggests that a long-term sea ice response to volcanic forcing is sensitive to the stability of the seawater column, wind, and ocean currents in the North Atlantic during the eruptions.     

Decadal‐to‐interdecadal fluctuations of arctic sea‐ice cover and the atmospheric circulation during 1954–1994

Abstract

The relationship between the Arctic and subarctic sea‐ice concentration (SIC) anomalies, particularly those associated with the decadal‐scale Greenland and Labrador Seas “Ice and Salinity Anomalies (ISAs) “, and the overlying atmospheric circulation fluctuations is investigated using the singular value decomposition (SVD) and composite map analysis methods. The data analyzed are monthly SIC and sea level pressure (SLP) anomalies, which cover the northern hemisphere poleward of 45°N and extend over the 41‐year period 1954–1994.

The SVD₁ (first) mode of the coupled variability, which accounts for 57% of the square covariance, is for the most part an atmosphere‐to‐ice forcing mode characterized by the decadal timescale. The aforementioned ISA anomalies are clearly captured by this mode whose SIC anomalies are dominated by a strong dipole across Greenland. However, as part of the same mode, there is also a weaker SIC dipole in the northern North Pacific which has opposite‐signed anomalies in the Sea of Okhotsk and the Bering Sea. It is also shown that there exists a significant negative correlation between the decadal SIC variability in the Greenland‐Barents Seas region associated with this mode and the North Atlantic Oscillation, whose spectrum also exhibits a quasi‐decadal signal.

The SVD₂ mode accounts for 12% of the square covariance and shows no evidence of a dominant forcing field of either SIC or SLP. This SVD mode exhibits very low frequency (interdecadal) variability, and its co‐variability is mainly concentrated in the northern North Pacific. It appears to be a high‐latitude extension of the recently investigated interdecadal North Pacific Oscillation. The spatial structure of the second mode complements the case of the first SVD mode whose co‐variability mainly occurs in the northern North Atlantic.