Regional variability of Antarctic sea ice extent

Interannual and seasonal variability of regional distribution of Antarctic sea ice extent is studied using monthly mean data on sea ice concentration in 1970-2012. The correlation is estimated between the variations in the area of floating ice in West and East Antarctica as well as in the Atlantic, Pacific, and Indian sectors of the Southern Ocean and the indices of atmospheric circuiation in the Southern Hemisphere.     

An assessment and interpretation of the observed warming of West Antarctica in the austral spring

We synthesize variability and trends in multiple analyses of Antarctic near-surface temperature representing several independent source datasets and spatially complete reconstructions, and place these into the broader context of the behavior of other components of the climate system during the past 30–50?years. Along with an annual-mean trend during the past 50?years of about 0.1°C/decade averaged over Antarctica, there is a distinct seasonality to the trends, with insignificant change (and even some cooling) in austral summer and autumn in East Antarctica, contrasting with warming in austral winter and spring. Apart from the Peninsula, the seasonal warming is largest and most significant in West Antarctica in the austral spring since the late 1970s. Concurrent trends in sea ice are independent evidence of the observed warming over West Antarctic, with the decrease in sea ice area in the Amundsen and Bellingshausen Seas congruent with at least 50% of the inland warming of West Antarctica. Trends in near surface winds and geopotential heights over the high-latitude South Pacific are consistent with a role for atmospheric forcing of the sea ice and air temperature anomalies. Most of the circulation trend projects onto the two Pacific South American (PSA) modes of atmospheric circulation variability, while the Southern Annular Mode lacks a positive trend in spring that would otherwise cause a cooling tendency. The largest circulation trend is associated with the PSA-1 mode, a wave-train extending from the tropics to the high Southern latitudes. The PSA-1 mode is significantly correlated with SSTs in the southwestern tropical and subtropical Pacific. The increased SSTs in this region, together with the observed increase in rainfall, suggest that anomalous deep convection has strengthened or increased the occurrence of the Rossby wave-train associated with PSA-1. This hypothesis is supported by results from two ensembles of SST-forced atmospheric general circulation model simulations. Finally, the implications of the seasonality, timing, and spatial patterns of Antarctic temperature trends with respect to interpreting the relative roles of stratospheric ozone depletion, SSTs and increased atmospheric concentrations of greenhouse gasses are discussed.     

Recent precipitation trends over the polar ice sheets

Summary Meteorological and glaciological analyses are integrated to examine the precipitation trends during the last three decades over the ice sheets covering Antarctica and Greenland. For Antarctica, the best data source is provided by glaciologically-measured trends of snow accumulation, and for limited sectors of East Antarctica consistency with precipitation amounts calculated from the atmospheric water balance equation is obtained. For Greenland, precipitation rates parameterized from atmospheric analyses yield the only comprehensive depiction. The precipitation rate over Antarctica appears to have increased by about 5% over a time period spanning the accumulation means for the 1955–65 to 1965–75 periods, while over Greenland it has decreased by about 15% since 1983 with a secondary increase over the southern part of the ice sheet starting in 1977. At the end of the 10-year overlapping period, the global sea-level impact of the precipitation changes over Antarctica dominates that for Greenland and yields a net ice-sheet precipitation contribution of roughly 0.02 mm yr^–1. These changes are likely due to marked variations in the cyclonic forcing affecting the ice sheets, but are only weakly reflected in the temperature regime, consistent with the episodic nature of cyclonic precipitation. These conclusions are not founded on high quality data bases. The importance of such changes for understanding global sea-level variations argues for a modest research effort to collect simultaneous meteorological and glaciological observations in order to describe and understand the current precipitation variations over both ice sheets. Some suggestions are offered for steps that could be taken.With 8 Figures     

A continuous simulation of global ice volume over the past 1 million years with 3-D ice-sheet models

Sea-level records show large glacial-interglacial changes over the past million years, which on these time scales are related to changes of ice volume on land. During the Pleistocene, sea-level changes induced by ice volume are largely caused by the waxing and waning of the large ice sheets in the Northern Hemisphere. However, the individual contributions of ice in the Northern and Southern Hemisphere are poorly constrained. In this study, for the first time a fully coupled system of four 3-D ice-sheet models is used, simulating glaciations on Eurasia, North America, Greenland and Antarctica. The ice-sheet models use a combination of the shallow ice and shelf approximations to determine sheet, shelf and sliding velocities. The framework consists of an inverse forward modelling approach to derive a self-consistent record of temperature and ice volume from deep-sea benthic δ¹⁸O data over the past 1 million years, a proxy for ice volume and temperature. It is shown that for both eustatic sea level and sea water δ¹⁸O changes, the Eurasian and North American ice sheets are responsible for the largest part of the variability. The combined contribution of the Antarctic and Greenland ice sheets is about 10 % for sea level and about 20 % for sea water δ¹⁸O during glacial maxima. However, changes in interglacials are mainly caused by melt of the Greenland and Antarctic ice sheets, with an average time lag of 4 kyr between melt and temperature. Furthermore, we have tested the separate response to changes in temperature and sea level for each ice sheet, indicating that ice volume can be significantly influenced by changes in eustatic sea level alone. Hence, showing the importance of a simultaneous simulation of all four ice sheets. This paper describes the first complete simulation of global ice-volume variations over the late Pleistocene with the possibility to model changes above and below present-day ice volume, constrained by observations of benthic δ¹⁸O proxy data.     

Albedo of coastal landfast sea ice in Prydz Bay,Antarctica: Observations and parameterization

The snow/sea-ice albedo was measured over coastal landfast sea ice in Prydz Bay, East Antarctica(off Zhongshan Station)during the austral spring and summer of 2010 and 2011. The variation of the observed albedo was a combination of a gradual seasonal transition from spring to summer and abrupt changes resulting from synoptic events, including snowfall, blowing snow, and overcast skies. The measured albedo ranged from 0.94 over thick fresh snow to 0.36 over melting sea ice. It was found that snow thickness was the most important factor influencing the albedo variation, while synoptic events and overcast skies could increase the albedo by about 0.18 and 0.06, respectively. The in-situ measured albedo and related physical parameters(e.g., snow thickness, ice thickness, surface temperature, and air temperature) were then used to evaluate four different snow/ice albedo parameterizations used in a variety of climate models. The parameterized albedos showed substantial discrepancies compared to the observed albedo, particularly during the summer melt period, even though more complex parameterizations yielded more realistic variations than simple ones. A modified parameterization was developed,which further considered synoptic events, cloud cover, and the local landfast sea-ice surface characteristics. The resulting parameterized albedo showed very good agreement with the observed albedo.     

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.     

Potential effect of nuclear war smokefall on sea ice

A large nuclear war could produce massive quantities of smoke from burning cities and industries. A portion of this smoke would fall out on Arctic sea ice, thus lowering its albedo and potentially increasing the solar energy absorbed by the ice and the snow that covers it. We use a one-dimensional thermodynamic sea ice model to examine the effect of smokefall on the seasonal variation of sea ice. In particular, we test the sensitivity of the model results to the time of year, duration, and latitude of smokefall.Sea ice thickness variations and the period of summer ice-free conditions are sensitive to the season of smokefall. The largest sea ice perturbations are generated by smokefall in spring. In this case the period of ice-free conditions during the summer can increase by 2 – 3.5 months between 67.5° N and 82.5° N. In any given season, the annual cycle of sea ice is not very sensitive to the duration of smokefall. The equilibrium annual cycle of sea ice variation is restored within a few years of smokefall when the smoke is flushed out of the ice/snow system.Since the sea ice model used here is not a comprehensive global climate model, it is difficult to predict the mid-latitude climate effects of the massive, but temporary, Arctic sea ice changes. However, our results suggest that future global climate model simulations of the effects of nuclear war smoke include interactive sea ice calculations.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Coupled ice‐ocean variability in the Greenland Sea

Abstract

A major surface feature of the Greenland Sea during winter is the frequent eastward extension of sea ice south of 75°N and an associated embayment to the north. These features are nominally connected with the East Greenland Current, and both the promontory and the embayment are readily apparent on climatic ice charts. However, there are significant changes in these features on time‐scales as short as a few days. Using a combination of satellite microwave images (SSM/I) of ice cover, meteorological data and in situ velocity, temperature and salinity records, we relate the ice distribution and its changes to the developing structure and circulation of the upper ocean during winter 1988–1989. Our measurements illustrate the preconditioning that leads to convective overturn, which in turn brings warmer water to the surface and results in the rapid disappearance of ice. In particular, the surface was cooled to the freezing point by early December and the salinity then increased through ice formation (about 0.016 m d^‐1) and brine rejection. Once the vertical density gradient was sufficiently eroded, a period of high heat flux (>300 W m^‐2) in late January provided enough buoyancy loss to convectively mix the upper water column to at least 200 m. We estimate vertical velocities at about 3 cm s^‐1 downward during the initial sinking. The deepening of the thermocline raised surface temperatures by over 1°C resulting in nearly 1.5 × 10⁵ km² of ice‐melt within two days. Average rates of ice retreat are about 11 km d^‐1 southwestward, generally consistent with a wind‐driven flow. Comparison of hydrographic surveys from before and after the overturning indicate the fresh water was advected out of the area, possibly to the south and east of our moorings.     

Regional Arctic sea ice variations as predictor for winter climate conditions

Seasonal prediction skill of winter mid and high northern latitudes climate from sea ice variations in eight different Arctic regions is analyzed using detrended ERA-interim data and satellite sea ice data for the period 1980–2013. We find significant correlations between ice areas in both September and November and winter sea level pressure, air temperature and precipitation. The prediction skill is improved when using November sea ice conditions as predictor compared to September. This is particularly true for predicting winter NAO-like patterns and blocking situations in the Euro-Atlantic area. We find that sea ice variations in Barents Sea seem to be most important for the sign of the following winter NAO—negative after low ice—but amplitude and extension of the patterns are modulated by Greenland and Labrador Seas ice areas. November ice variability in the Greenland Sea provides the best prediction skill for central and western European temperature and ice variations in the Laptev/East Siberian Seas have the largest impact on the blocking number in the Euro-Atlantic region. Over North America, prediction skill is largest using September ice areas from the Pacific Arctic sector as predictor. Composite analyses of high and low regional autumn ice conditions reveal that the atmospheric response is not entirely linear suggesting changing predictive skill dependent on sign and amplitude of the anomaly. The results confirm the importance of realistic sea ice initial conditions for seasonal forecasts. However, correlations do seldom exceed 0.6 indicating that Arctic sea ice variations can only explain a part of winter climate variations in northern mid and high latitudes.     

Sensitivity of last glacial maximum climate to sea ice conditions in the Nordic Seas

Published reconstructions of last glacial maximum (LGM) sea surface temperatures and sea ice extent differ significantly. We here test the sensitivity of simulated North Atlantic climates to two different reconstructions by using these reconstructions as boundary conditions for model experiments. An atmospheric general circulation model has been used to perform two simulations of the (LGM) and a modern-day control simulation. Standard (CLIMAP) reconstructions of sea ice and sea surface temperatures have been used for the first simulation, and a set of new reconstructions in the Nordic Seas/Northern Atlantic have been used for the second experiment. The new reconstruction is based on 158 core samples, and represents ice-free conditions during summer in the Nordic Seas, with accordingly warmer sea surface temperatures and less extensive sea ice during winter as well. The simulated glacial climate is globally 5.7 K colder than modern day, with the largest changes at mid and high latitudes. Due to more intense Hadley circulation, the precipitation at lower latitudes has increased in the simulations of the LGM. Relative to the simulation with the standard CLIMAP reconstructions, reduction of the sea ice in the North Atlantic gives positive local responses in temperature, precipitation and reduction of the sea level pressure. Only very weak signatures of the wintertime Icelandic Low occur when the standard CLIMAP sea surface temperature reconstruction is used as the lower boundary condition in LGM. With reduced sea ice conditions in the Nordic Seas, the Icelandic Low becomes more intense and closer to its present structure. This indicates that thermal forcing is an important factor in determining the strength and position of the Icelandic Low. The Arctic Oscillation is the most dominant large scale variability feature on the Northern Hemisphere in modern day winter climate. In the simulation of the LGM with extensive sea ice this pattern is significantly changed and represents no systematic large scale variability over the North Atlantic. Reduction of the North Atlantic sea ice extent leads to stronger variability in monthly mean sea level pressure in winter. The synoptic variability appears at a lower level in the simulation when standard reconstructions of the sea surface in the LGM are used. A closer inspection of storm tracks in this model experiment shows that that the synoptic lows follow a narrow band along the ice edge during winter. The trajectories of synoptic lows are not constrained to the sea ice edge to the same degree when the sea ice extent is reduced. Seasonally open waters in the Nordic Seas in the new reconstruction apparently act as a moisture source, consistent with the current understanding of the rapid growth of the Fennoscandian and Barents Ice Sheets, during the LGM. The signal from the intensified thermal forcing in the North Atlantic in Boreal winter is carried zonally by upper tropospheric waves, and thus generates non-local responses to the changed sea ice cover.     

The effect of Antarctic sea ice on the Southern Hemisphere atmosphere during the southern summer

This study examines the influence of Antarctic sea ice distribution on the large scale circulation of the Southern Hemisphere using a fully coupled GCM where the sea ice submodel is replaced by a climatology of observed extremes in sea ice concentration. Three 150-year simulations were completed for maximum, minimum and average sea ice concentrations and the results for the austral summer (January?CMarch) were compared using the surface temperatures forced by the sea ice distributions as a filter for creating the composite differences. The results indicate that in the austral summer the polar cell expands (contracts) under minimum (maximum) sea ice conditions with corresponding shifts in the midlatitude Ferrell cell. We suggest that this response occurs because sea ice lies in the margin between the polar and midlatitude cells. The polarity of the Southern Hemisphere Annular (SAM) mode is also influenced such that when sea ice is at a minimum (maximum) the polarity of the SAM tends to be negative (positive).     

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.     

Tropical cooling and the isotopic composition of precipitation in general circulation model simulations of the ice age climate

 We test the climate effects of changes in the tropical ocean by imposing three different patterns of tropical SSTs in ice age general circulation model simulations that include water source tracers and water isotope tracers. The continental air temperature and hydrological cycle response in these simulations is substantial and should be directly comparable to the paleoclimatic record. With tropical cooling imposed, there is a strong temperature response in mid- to high-latitudes resulting from changes in sea ice and disturbance of the planetary waves; the results suggest that tropical/subtropical ocean cooling leads to significant dynamical and radiative feedbacks that might amplify ice age cycles. The isotopes in precipitation generally follow the temperature response at higher latitudes, but regional δ¹⁸O/air temperature scaling factors differ greatly among the experiments. In low-latitudes, continental surface temperatures decrease congruently with the adjacent SSTs in the cooling experiments. Assuming CLIMAP SSTs, ¹⁸O/¹⁶O ratios in low-latitude precipitation show no change from modern values. However, the experiments with additional cooling of SSTs produce much lower tropical continental δ¹⁸O values, and these low values result primarily from an enhanced recycling of continental moisture (as marine evaporation is reduced). The water isotopes are especially sensitive to continental aridity, suggesting that they represent an effective tracer of the extent of tropical cooling and drying. Only one of the tropical cooling simulations produces generalized low-latitude aridity. These results demonstrate that the geographic pattern of cooling is most critical for promoting much drier continents, and they underscore the need for accurate reconstructions of SST gradients in the ice age ocean. Received: 26 July 1999 / Accepted: 10 July 2000     

Effective sea ice area based on a thickness threshold

As Arctic sea ice declines in response to climate change, a shift from thick multiyear ice to a thinner ice cover is occurring. With this transition, ice thicknesses approach a threshold below which ice no longer insulates the atmosphere from oceanic surface fluxes. While this is well known, there are no estimates of the magnitude of this threshold, nor of the proportion of sea ice area that is below this threshold as ice thins. We determine this threshold by simulating the atmospheric response to varying thicknesses, ranging from 0.0 to 2.0 m and determine that threshold to be 0.40–0.50 m. The resulting “effective” ice area is 4–14% lower than reported total ice area, as 0.39–0.97 × 10⁶ km² of the total ice area falls below the threshold throughout the twentieth century, including during notable ice minima. The atmosphere above large non-insulating ice-covered regions is susceptible to more than 2 °C of warming despite ice presence. Observed mean Arctic Ocean ice thickness is projected to fall below this threshold as early as the mid-2020s. Studies on ocean–atmosphere interactions in relation to sea ice area should focus on this insulating sea ice area, where ice is at least 0.40–0.50 m thick, and treat ice regions below 0.40–0.50 m thickness with caution.

     

Stable isotopes in surface snow along a traverse route from Zhongshan station to Dome A, East Antarctica

Samples of surface snow were collected for stable isotope analysis along the traverse route from Zhongshan to Dome A (East Antarctica) from Dec 28th, 2007 to Feb. 8th, 2008. The local relationship between δD and surface temperature is established to be 6.4 ± 0.2 ‰ per °C, very similar to the average for East Antarctic. The deuterium excess shows a pattern of high values over Antarctica, particularly at Dome A. We compare our data with an atmospheric general circulation model which includes stable water isotopes (ECHAM5-wiso). The model simulation captures the right levels of δD, but overestimates δ¹⁸O. This study provides support for the ongoing deep ice core project at Dome A.     

The initiation of ice sheet growth,Milankovitch solar radiation variations,and the 100 ky ice age cycle

Much work has gone into deciphering the causes of the large scale glacial/interglacial variations in the climate system over the last 900 000 years. While variations on the 41 thousand year (ky) and 23 ky time scales seem to be linearly linked to the variations in the distribution of solar radiation at the top of the atmosphere, Milankovitch solar radiation variations, the causes of the dominant 100 ky cycle in the geologic record are still unknown. One of the aspects of this cycle that is not well understood is how large scale ice sheet growth is initiated. Here we describe the mechanisms by which large scale ice sheet growth may have been initiated by the changes in the seasonal and latitudinal distribution of solar radiation over the past 160 ky. This is done through the use of a coupled energy balance climate-thermodynamic sea ice model that includes a hydrologic cycle which computes precipitation, and a land surface energy balance which determines the net accumulation of snow and ice. Results indicate that the initiation of ice sheet growth is possible during times of extremely low summer solstice solar radiation as a result of a large decrease in ablation during the critical melt season.     

Coupled climate modelling of ocean circulation changes during ice age inception

Freshening of high latitude surface waters can change the large-scale oceanic transport of heat and salt. Consequently, atmospheric and sea ice perturbations over the deep water production sites excite a large-scale response establishing an oceanic "teleconnection" with time scales of years to centuries. To study these feedbacks, a coupled atmosphere-ocean-sea ice model consisting of a two dimensional atmospheric energy and moisture balance model (EMBM) coupled to a thermodynamic sea ice model and an ocean general circulation model is utilised. The coupled model reproduces many aspects of the present oceanic circulation. We also investigate the climate impact of changes in fresh water balance during an ice age initiation. In this experiment part of the precipitation over continents is stored within continental ice sheets. During the buildup of ice sheets the oceanic stratification in the North Atlantic is weakened by a reduced continental run-off leading to an enhanced thermohaline circulation. Under these conditions salinity is redistributed such that deep water is more saline than under present conditions. Once the ice sheets built up, we simulate an ice age climate without net fresh water storage on the continents. In this case the coupled model reproduces the shallow and weak overturning cell, an ice edge advance insulating the upper ocean, and many other aspects of the glacial circulation.     

Climatic variations of ice conditions in different regions of the Barents Sea

Presented are the results of studying the regional peculiarities of climatic variations of spatiotemporal distribution of ice in the Barents Sea water area in 1977?C2010. Demonstrated is the dynamics of the interannual and seasonal variability of main elements of the ice regime (ice cover area, ice edge position, and ice period duration). Revealed are the common features and differences in the ice conditions in the water areas under study. It has been found that there is a significant feedback between the specific ice coverage in different areas of the sea. The climatic variations of the total ice coverage of the Barents Sea for the period of 1960?C2010 are analyzed using the electronic database on the Barents Sea ice coverage. It can be supposed that the current warm phase of climatic variations in the Barents Sea is coming to the end.     

Influence of coupling on atmosphere, sea ice and ocean regional models in the Ross Sea sector, Antarctica

Air–sea ice–ocean interactions in the Ross Sea sector form dense waters that feed the global thermohaline circulation. In this paper, we develop the new limited-area ocean–sea ice–atmosphere coupled model TANGO to simulate the Ross Sea sector. TANGO is built up by coupling the atmospheric limited-area model MAR to a regional configuration of the ocean–sea ice model NEMO. A method is then developed to identify the mechanisms by which local coupling affects the simulations. TANGO is shown to simulate realistic sea ice properties and atmospheric surface temperatures. These skills are mostly related to the skills of the stand alone atmospheric and oceanic models used to build TANGO. Nonetheless, air temperatures over ocean and winter sea ice thickness are found to be slightly improved in coupled simulations as compared to standard stand alone ones. Local atmosphere ocean feedbacks over the open ocean are found to significantly influence ocean temperature and salinity. In a stand alone ocean configuration, the dry and cold air produces an ocean cooling through sensible and latent heat loss. In a coupled configuration, the atmosphere is in turn moistened and warmed by the ocean; sensible and latent heat loss is therefore reduced as compared to the stand alone simulations. The atmosphere is found to be less sensitive to local feedbacks than the ocean. Effects of local feedbacks are increased in the coastal area because of the presence of sea ice. It is suggested that slow heat conduction within sea ice could amplify the feedbacks. These local feedbacks result in less sea ice production in polynyas in coupled mode, with a subsequent reduction in deep water formation.     

Ice and climate modeling: An editorial essay

The growth and decay of ice sheets are driven by forces affecting the seasonal cycles of snowfall and snowmelt. The external forces are likely to be variations in the earth's orbit which cause differences in the solar radiation received. Radiational control of snowmelt is modulated by the seasonal cycles of snow albedo and cloud cover. The effects of orbital changes can be magnified by feedbacks involving atmospheric CO₂ content, ocean temperatures and desert areas. Climate modeling of the causes of the Pleistocene ice ages involves modeling the interactions of all components of the climate system; snow, sea ice, glacier ice, the ocean, the atmosphere, and the solid earth. Such modeling is also necessary for interpreting oxygen isotope records from ice and ocean as paleoclimatic evidence.