Impact of dust aerosol on glacial-interglacial climate

The temperature anomaly and dust concentrations recorded from central Antarctic ice core records display a strong negative correlation. The dust concentration recorded from an ice core in central Antarctica is 50-70 times higher during glacial periods than interglacial periods. This study investigated the impact of dust aerosol on glacial-interglacial climate, using a zonal energy balance model and dust concentration data from an Antarctica ice core. Two important effects of dust, the direct radiative effect and dust-albedo feedback, were considered. On the one hand, the direct radiative effect of dust significantly cooled the climate during the glacial period, with cooling during the last glacial maximum being as much as 2.05℃ in Antarctica. On the other hand, dust deposition onto the ice decreased the surface albedo over Antarctica, leading to increased absorption of solar radiation, inducing a positive feedback that warmed the region by as much as about 0.9℃ during the glacial period. However, cooling by the direct dust effect was found to be the controlling effect for the glacial climate and may be the major influence on the strong negative correlation between temperature and dust concentration during glacial periods.     

The effect of snow on Antarctic sea ice simulations in a coupled atmosphere-sea ice model

 The effect of a snow cover on sea ice accretion and ablation is estimated based on the ‘zero-layer’ version sea ice model of Semtner, and is examined using a coupled atmosphere-sea ice model including feedbacks and ice dynamics effects. When snow is disregarded in the coupled model the averaged Antarctic sea ice becomes thicker. When only half of the snowfall predicted by the atmospheric model is allowed to land on the ice surface sea ice gets thicker in most of the Weddell and Ross Seas but thinner in East Antarctic in winter, with the average slightly thicker. When twice as much snowfall as predicted by the atmospheric model is assumed to land on the ice surface sea ice also gets much thicker due to the large increase of snow-ice formation. These results indicate the importance of the correct simulation of the snow cover over sea ice and snow-ice formation in the Antarctic. Our results also illustrate the complex feedback effects of the snow cover in global climate models. In this study we have also tested the use of a mean value of 0.16 Wm^-1 K^-1 instead of 0.31 for the thermal conductivity of snow in the coupled model, based on the most recent observations in the eastern Antarctic and Bellingshausen and Amundsen Seas, and have found that the sea ice distribution changes greatly, with the ice becoming much thinner by about 0.2 m in the Antarctic and about 0.4 m in the Arctic on average. This implies that the magnitude of the thermal conductivity of snow is of considerable importance for the simulation of the sea ice distribution. An appropriate value of the thermal conductivity of snow is as crucial as the depth of the snow layer and the snowfall rate in a sea ice model. The coupled climate models require accurate values of the effective thermal conductivity of snow from observations for validating the simulated sea ice distribution under the present climate conditions. Received: 20 November 1997/Accepted: 27 July 1998     

Comments on the origin of dust in East Antarctica for present and ice age conditions

We have studied the distribution of 327 clay mineral particles retrieved from four Antaretic ice smaples corresponding to present and Last Glacial Maximum (LGM) climate conditions. Illite, chlorite, smectite and kaolinite were identified in all samples. Focusing on kaolinite, because of its use as a possible tracer of low latitude soils, we find a significantly smaller amount for LGM samples while the dust concentration in snow during the LGM was about 30 times higher than for present climate conditions. This can be interpreted as change in the contribution of the Australian source with climate.A second approach was based on the modeling of the desert dust cycle using an Atmospheric General Circulation Model (AGCM) under both present-day and ice age conditions. Unlike mineralogical results, the model suggests the prevalence of the Australian dust source in the deposits over East Antarctica under both present-day and LGM climate conditions. However the model fails to reproduce the strong increase in dust deposits during the LGM. This discrepancy could be partly due to the lack of a higher latitude dust source in the model.The stronger dust input recorded in ice cores for the LGM could be related to an additional active high latitude source (possibly close to South America) overlapping the atmospheric background coming from low latitude areas.     

Ensemble of sea ice initial conditions for interannual climate predictions

Polar climate studies are severely hampered by the sparseness of the sea ice observations. We aim at filling this critical gap by producing two 5-member sea ice historical simulations strongly constrained by ocean and atmosphere observational data and covering the 1958–2006 and 1979–2012 periods. This is the first multi-member sea ice reconstruction covering more than 50 years. The obtained sea ice conditions are in reasonable agreement with the few available observations. These best estimates of sea ice conditions serve subsequently as initial sea ice conditions for a set of 28 3-year-long retrospective climate predictions. We compare it to a set in which the sea ice initial conditions are taken from a single-member sea ice historical simulation constrained by atmosphere observations only. We find an improved skill in predicting the Arctic sea ice area and Arctic near surface temperature but a slightly degraded skill in predicting the Antarctic sea ice area. We also obtain a larger spread between the members for the sea ice variables, thus more representative of the forecast error.     

A simple model for air/snow fractionation of aerosol components over the Antarctic Peninsula

A model has been set up to investigate the wet and dry aerosol removal processes which occur in clean air over the Antarctic Peninsula. Input for the model was obtained from bulk chemical analysis and scanning electron microscopy of aerosol and snow samples collected simultaneously at remote sites around the Peninsula. The model predicts that sulphate and sea-salt aerosol will be removed mainly in-cloud by riming of falling snow and ice crystals. Crustal aerosol is principally removed by acting as nuclei for these crystals and by impaction on falling snow. For the largest, locally-generated aerosol dry deposition is indicated as the major removal process. These findings suggest a possible mechanism for the observed air/snow fractionation.     

Response of the Antarctic ice sheet to future greenhouse warming

Possible future changes in land ice volume are mentioned frequently as an important aspect of the greenhouse problem. This paper deals with the response of the Antarctic ice sheet and presents a tentative projection of changes in global sea level for the next few hundred years, due to changes in its surface mass balance. We imposed a temperature scenario, in which surface air temperature rises to 4.2° C in the year 2100 AD and is kept constant afterwards. As GCM studies seem to indicate a higher temperature increase in polar latitudes, the response to a more extreme scenario (warming doubled) has also been investigated. The mass balance model, driven by these temperature perturbations, consists of two parts: the accumulation rate is derived from present observed values and is consequently perturbed in proportion to the saturated vapour pressure at the temperature above the inversion layer. The ablation model is based on the degree-day method. It accounts for the daily temperature cycle, uses a different degree-day factor for snow and ice melting and treats refreezing of melt water in a simple way. According to this mass balance model, the amount of accumulation over the entire ice sheet is presently 24.06 × 10¹¹ m³ of ice, and no runoff takes place. A 1°C uniform warming is then calculated to increase the overall mass balance by an amount of 1.43 × 10¹¹ m³ of ice, corresponding to a lowering of global sea level with 0.36 mm/yr. A temperature increase of 5.3°C is needed for the increase in ablation to become more important than the increase in accumulation and the temperature would have to rise by as much as 11.4°C to produce a zero surface mass balance. Imposing the Bellagio-scenario and accumulating changes in mass balance forward in time (static response) would then lower global sea level by 9 cm by 2100 AD. In a subsequent run with a high-resolution 3-D thermomechanic model of the ice sheet, it turns out that the dynamic response of the ice sheet (as compared to the direct effect of the changes in surface mass balance) becomes significant after 100 years or so. Ice-discharge across the grounding-line increases, and eventually leads to grounding-line retreat. This is particularly evident in the extreme case scenario and is important along the Antarctic Peninsula and the overdeepened outlet glaciers along the East Antarctic coast. Grounding-line retreat in the Ross and Ronne-Filchner ice shelves, on the other hand, is small or absent.     

Impact of Greenland and Antarctic ice sheet interactions on climate sensitivity

We use the Earth system model of intermediate complexity LOVECLIM to show the effect of coupling interactive ice sheets on the climate sensitivity of the model on a millennial time scale. We compare the response to a 2×CO₂ warming scenario between fully coupled model versions including interactive Greenland and Antarctic ice sheet models and model versions with fixed ice sheets. For this purpose an ensemble of different parameter sets have been defined for LOVECLIM, covering a wide range of the model??s sensitivity to greenhouse warming, while still simulating the present-day climate and the climate evolution over the last millennium within observational uncertainties. Additional freshwater fluxes from the melting ice sheets have a mitigating effect on the model??s temperature response, leading to generally lower climate sensitivities of the fully coupled model versions. The mitigation is effectuated by changes in heat exchange within the ocean and at the sea?Cair interface, driven by freshening of the surface ocean and amplified by sea?Cice-related feedbacks. The strength of the effect depends on the response of the ice sheets to the warming and on the model??s climate sensitivity itself. The effect is relatively strong in model versions with higher climate sensitivity due to the relatively large polar amplification of LOVECLIM. With the ensemble approach in this study we cover a wide range of possible model responses.     

Atmospheric inversion strength over polar oceans in winter regulated by sea ice

Low-level temperature inversions are a common feature of the wintertime troposphere in the Arctic and Antarctic. Inversion strength plays an important role in regulating atmospheric processes including air pollution, ozone destruction, cloud formation, and negative longwave feedback mechanisms that shape polar climate response to anthropogenic forcing. The Atmospheric Infrared Sounder (AIRS) instrument provides reliable measures of spatial patterns in mean wintertime inversion strength when compared with available radiosonde observations and reanalysis products. Here, we examine the influence of sea ice concentration on inversion strength in the Arctic and Antarctic. Correlation of inversion strength with mean annual sea ice concentration, likely a surrogate for the effective thermal conductivity of the wintertime ice pack, yields strong, linear relationships in the Arctic (r?=?0.88) and Antarctic (r?=?0.86). We find a substantially greater (stronger) linear relationship between sea ice concentration and surface air temperature than with temperature at 850?hPa, lending credence to the idea that sea ice controls inversion strength through modulation of surface heat fluxes. As such, declines in sea ice in either hemisphere may imply weaker mean inversions in the future. Comparison of mean inversion strength in AIRS and global climate models (GCMs) suggests that many GCMs poorly characterize mean inversion strength at high latitudes.     

Climate modification by future ice sheet changes and consequences for ice sheet mass balance

The future evolution of global ice sheets under anthropogenic greenhouse forcing and its impact on the climate system, including the regional climate of the ice sheets, are investigated with a comprehensive earth system model consisting of a coupled Atmosphere–Ocean General Circulation Model, a dynamic vegetation model and an ice sheet model. The simulated control climate is realistic enough to permit a direct coupling of the atmosphere and ice sheet components, avoiding the use of anomaly coupling, which represents a strong improvement with respect to previous modelling studies. Glacier ablation is calculated with an energy-balance scheme, a more physical approach than the commonly used degree-day method. Modifications of glacier mask, topographic height and freshwater fluxes by the ice sheets influence the atmosphere and ocean via dynamical and thermodynamical processes. Several simulations under idealized scenarios of greenhouse forcing have been performed, where the atmospheric carbon dioxide stabilizes at two and four times pre-industrial levels. The evolution of the climate system and the ice sheets in the simulations with interactive ice sheets is compared with the simulations with passively coupled ice sheets. For a four-times CO₂ scenario forcing, a faster decay rate of the Greenland ice sheet is found in the non-interactive case, where melting rates are higher. This is caused by overestimation of the increase in near-surface temperature that follows the reduction in topographic height. In areas close to retreating margins, melting rates are stronger in the interactive case, due to changes in local albedo. Our results call for careful consideration of the feedbacks operating between ice sheets and climate after substantial decay of the ice sheets.     

A thermomechanical model of ice flow in West Antarctica

 This study uses a three-dimensional thermo-mechanical model to investigate the internal flow dynamics of the West Antarctic Ice Sheet (WAIS). The model allows ice thickness, flow and temperature to interact freely. Its domain is prescribed as that of the present-day grounded WAIS. Realistic present-day climatic and topographical boundary conditions are employed. The analysis of a series of experiments pays particular attention to the location and dynamics of concentrations of ice flow (ice streams). Underlying topographic troughs are crucial in determining the strength and location of these concentrations of flow. The flow pattern generated by subglacial troughs is made more distinct by the inclusion of ice flow/temperature coupling. The inclusion of sliding leads to the generation of limit cycles in the ice flow. They are concentrated around the present-day ice streams B and C of the Siple Coast and have a period of 5 to 10 ky. There appears to be competition between several preferred ice flow pathways in this area. The two end members of the flow regime are a strong ice stream C with a weakened ice stream A/B complex, and strong ice streams A and B with a dormant ice stream C. Ice streams appear to require ice discharges above a certain threshold in order to maintain frictional heat generation and fast flow. Individual ice streams can therefore interact through changes in catchment-area size: a reduction in catchment area reduces the volume of ice entering a stream and can cause stagnation as the amount of frictional heating falls. Received: 22 July 1997/Accepted: 27 July 1998     

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.     

Soluble and insoluble impurities along the 950 m deep Vostok ice core (Antarctica) — Climatic implications

Simultaneous measurements of soluble and insoluble impurities were made on the 950 m deep Vostok (78°30′S, 106°54′E, 3420 m a.s.l.) ice core, spanning roughly 50000 yr, using various analytical techniques. We observed higher continental (×37) and marine (×5.1) inputs during the last glacial age than during the Holocene stage. A study of microparticle compositions and of volcanic indicators (Zn, H₂SO₄), shows that the high observed crustal input is not due to enhanced volcanism, but is rather of continental eolian origin. For the first time, the ionic balance along a deep ice core is established, mainly used in discussing the evolution of the Cl to Na ratio over central East Antarctica with changing climatic conditions: the presence of relatively high amounts of Na₂SO₄ in the marine aerosol at the Vostok site during the Holocene is demonstrated. Comparison with the Dome C (74°39′S, 124°10′E, 3040 m a.s.l.) results confirms the chronology of the major events: (i) maximum terrestrial input around the last glacial maximum (～18 ka BP); (ii) end of the high continental flux over Antarctica near 13 ka BP; (iii) marine input varying in an opposing manner to isotopic fluctuations with rather high concentrations beginning to decrease when isotopic values increase and reaching Holocene values at the end of the transition between cold and warmer climate conditions. Detailed comparison with results provided by deep ice cores from other sites which are probably more influenced by oceanic air masses seems to indicate that most of the aerosol reaching central East Antarctica travel over large distance probably at rather high altitude through the troposphere. We can consider that central East Antarctica is well representative of the upper part of the troposphere (higher than i.e., 3000 m) and should, therefore, provide valuable data for global and Antarctic paleoclimatological models.     

Simulation of sea ice in FGOALS-g2: Climatology and late 20th century changes

Sea ice is an important component in the Earth’s climate system. Coupled climate system models are indispensable tools for the study of sea ice, its internal processes, interaction with other components, and projection of future changes. This paper evaluates the simulation of sea ice by the Flexible Global Ocean-Atmosphere-Land System model Grid-point Version 2 (FGOALS-g2), in the fifth phase of the Coupled Model Inter-comparison Project (CMIP5), with a focus on historical experiments and late 20th century simulation. Through analysis, we find that FGOALS-g2 produces reasonable Arctic and Antarctic sea ice climatology and variability. Sea ice spatial distribution and seasonal change characteristics are well captured. The decrease of Arctic sea ice extent in the late 20th century is reproduced in simulations, although the decrease trend is lower compared with observations. Simulated Antarctic sea ice shows a reasonable distribution and seasonal cycle with high accordance to the amplitude of winter-summer changes. Large improvement is achieved as compared with FGOALS-g1.0 in CMIP3. Diagnosis of atmospheric and oceanic forcing on sea ice reveals several shortcomings and major aspects to improve upon in the future: (1) ocean model improvements to remove the artificial island at the North Pole; (2) higher resolution of the atmosphere model for better simulation of important features such as, among others, the Icelandic Low and westerly wind over the Southern Ocean; and (3) ocean model improvements to accurately receive freshwater input from land, and higher resolution for resolving major water channels in the Canadian Arctic Archipelago.     

The effects of aggressive mitigation on steric sea level rise and sea ice changes

With an increasing political focus on limiting global warming to less than 2 °C above pre-industrial levels it is vital to understand the consequences of these targets on key parts of the climate system. Here, we focus on changes in sea level and sea ice, comparing twenty-first century projections with increased greenhouse gas concentrations (using the mid-range IPCC A1B emissions scenario) with those under a mitigation scenario with large reductions in emissions (the E1 scenario). At the end of the twenty-first century, the global mean steric sea level rise is reduced by about a third in the mitigation scenario compared with the A1B scenario. Changes in surface air temperature are found to be poorly correlated with steric sea level changes. While the projected decreases in sea ice extent during the first half of the twenty-first century are independent of the season or scenario, especially in the Arctic, the seasonal cycle of sea ice extent is amplified. By the end of the century the Arctic becomes sea ice free in September in the A1B scenario in most models. In the mitigation scenario the ice does not disappear in the majority of models, but is reduced by 42 % of the present September extent. Results for Antarctic sea ice changes reveal large initial biases in the models and a significant correlation between projected changes and the initial extent. This latter result highlights the necessity for further refinements in Antarctic sea ice modelling for more reliable projections of future sea ice.     

Recent southern Indian Ocean climate variability inferred from a Law Dome ice core: new insights for the interpretation of coastal Antarctic isotopic records

Stable isotopes in water have been measured along a very high accumulation ice core from Law Dome on the east Antarctic coast. These enable a detailed comparison of the isotopic records over sixty years (1934–1992) with local (Antarctic station data) and remote meteorological observations (atmospheric reanalyses and sea-surface temperature estimates) on a seasonal to inter-annual time scale. Using both observations and isotopic atmospheric general circulation model (GCM) results, we quantify the relationships between stable isotopes (d ¹⁸O, dD and deuterium excess; d = dD –8 × d ¹⁸O) with site and source temperature at seasonal and decadal time scales, showing the large imprint of source conditions on Law Dome isotopes. These calibrations provide new insights for the quantitative interpretation of temporal isotopic fluctuations from coastal Antarctic ice cores. An abrupt change in the local meridional atmospheric circulation is clearly identified from Law Dome deuterium excess during the 1970s and analysed using GCM simulations.     

The demise of Arctic sea ice during stabilisation at high greenhouse gas concentrations

A climate model experiment was conducted using the HadCM3 climate model and a scenario in which the atmospheric CO₂ concentration was increased over 70 years from pre-industrial concentrations to 4 times this level and then stabilised for more than a 1,000 years. During the period of stabilisation the global atmospheric surface temperatures continued to rise as the deep oceans adjusted towards a new equilibrium. However, even after 1,000 years this new equilibrium had not been reached. During the first 600 years, Arctic and Antarctic winter sea ice thickness and area covered declined with a significant impact on the global radiation budget. After this period the area of the Arctic covered by sea ice entered a 150 years period during which time it underwent a series of oscillations. Following the oscillation the centre of the Arctic basin became ice free throughout the year. A sensitivity experiment demonstrates that although the sea ice extent can be greatly reduced through the artificial heating of the mixed layer, prior to the onset of the oscillatory phase the ice recovers over 15 years. Understanding the causes of this oscillatory phase may elucidate the mechanisms of variability in the Arctic in the present climate and in future policy relevant scenarios. We have investigated the atmospheric and oceanic forcing on the ice during the oscillatory phase, and find that the behaviour is linked to a redistribution of Arctic Ocean heat stores.     

Comparative inverse analysis of satellite (MODIS) and ground (PM10) observations to estimate dust emissions in East Asia

Soil dust aerosol is the largest contributor to aerosol mass concentrations in the troposphere and has considerable effects on air quality and climate. Arid and semi-arid areas of East Asia are one of the important dust source regions thus it is crucial to understand dust mobilization and accurately estimate dust emissions in East Asia. However, present dust models still contain large uncertainties with dust emissions that remain a significant contributor to the overall uncertainties in the model. In this study, we attempt to reduce these uncertainties by using an inverse modeling technique and obtain optimized dust emissions. We use Moderate Resolution Imaging Spectrometer (MODIS) aerosol optical depths (AODs) and groundbased mass concentrations of particles less than 10 μm in aerodynamic diameter (PM₁₀) observations over East Asia in May 2007. The MODIS AODs are validated with AErosol RObotic NETwork (AERONET) AODs. The inversion uses the maximum a posteriori method and the GEOS-Chem chemical transport model (CTM) as a forward model. The model error is large over dust source regions including the Gobi Desert and Mongolia. We find that inverse modeling analyses from the MODIS and PM₁₀ observations consistently result in decrease of dust emissions over Mongolia and the Gobi Desert. Whereas over the Taklamakan Desert and Manchuria, the inverse modeling analyses from both observations yield contrast results such as increase of dust sources using MODIS AODs, while decrease of those using PM₁₀ observations. We discuss some limitations of both observations to obtain the optimized dust emissions and suggest several strategies for the improvement of dust emission estimates in the model.     

An investigation of the small ice cap instability in the Southern Hemisphere with a coupled atmosphere-sea ice-ocean-terrestrial ice model

 A simple climate model has been developed to investigate the existence of the small ice cap instability in the Southern Hemisphere. The model consists of four coupled components: an atmospheric energy balance model, a thermodynamic snow-sea ice model, an oceanic mixed layer model and a terrestrial ice model. Results from a series of experiments involving different degrees of coupling in the model show that the instability appears only in those cases when an explicit representation of the Antarctic ice sheet is not included in the model. In order to determine which physical processes in the ice sheet model lead to a stabilization of the system we have conducted several sensitivity experiments in each of which a given ice sheet process has been removed from the control formulation of the model. Results from these experiments suggest that the feedback between the elevation of the ice sheet and the snow accumulation-ice ablation balance is responsible for the disappearance of the small ice cap instability in our simulation. In the model, the mass balance of the ice sheet depends on the air temperature at sea level corrected for altitude and it is, therefore, a function of surface elevation. This altitude-mass balance feedback effectively decouples the location of the ice edge from any specific sea level isotherm, thus decreasing the model sensitivity to the albedo-temperature feedback, which is responsible for the appearance of the instability. It is also shown that the elevation-radiative cooling feedback tends to stabilize the ice sheet, although its effect does not seem to be strong enough to remove the instability. Another interesting result is that for those simulations which include the terrestrial ice model with elevation-dependent surface mass balance, hysteresis is exhibited, where for a given level of external forcing, two stable solutions with different, non-zero ice-sheet volume and area and different air and ocean temperature fields occur. However, no unstable transition between the two solutions is ever observed. Our results suggest that the small ice cap instability mechanism could be unsuitable for explaining the inception of glaciation in Antarctica. Received: 14 April 1997 / Accepted: 22 October 1997     

过去0.8 Ma期间的地表温度数值模拟研究

使用一个改进的二维能量平衡模式模拟了过去0.8 Ma冰期-间冰期旋回期间北半球各纬度带的地表温度,并以65°N的地表温度为代表与南极冰芯记录进行了比较.通过敏感性试验,分析了日射量、温室气体、沙尘气溶胶强迫和水汽反馈的辐射-气候效应.结果显示,日射量变化不足以解释冰期-间冰期旋回期间北半球的地表温度变化,大气温室气体(...