Hostname: page-component-7c8c6479df-24hb2 Total loading time: 0 Render date: 2024-03-29T06:53:42.070Z Has data issue: false hasContentIssue false

Thinning of the Ice Sheet Estimated from Total Gas Content of Ice Cores in Mizuho Plateau, East Antarctica

Published online by Cambridge University Press:  20 January 2017

Takao Kameda
Affiliation:
Institute of Low Temperature Science, Hokkaido University,Sapporo 060, Japan
Masayoshi Nakawo
Affiliation:
Nagaoka Institute of Snow and Ice Studies, National Research Center for Disaster Prevention, STA, Nagaoka 940, Japan
Shinji Mae
Affiliation:
Department of Applied Physics, Faculty of Engineering, Hokkaido University, Sapporo 060, Japan
Okitsugu Watanabe
Affiliation:
National Institute of Polar Research, 9–10, Kaga 1-chome, Itabashi-ku, Tokyo 173, Japan
Renji Naruse
Affiliation:
Institute of Low Temperature Science, Hokkaido University,Sapporo 060, Japan
Rights & Permissions [Opens in a new window]

Abstract

A linear relation between total gas content in ice and the elevation of ice formation (i.e. pore close-off) was obtained from seven shallow ice cores in Mizuho Plateau, Antarctica. The derived relation was applied to the vertical profile of total gas content in a 700 m long ice core at Mizuho Station. A general trend of gradual increase in total gas content was observed from 600 to 200 m in depth, from which toward the upper layer it showed a steep increase. After eliminating the effect of down-slope flow of ice around Mizuho Station, it was estimated that the thickness of the ice sheet decreased by about 350 m at maximum during the last 2000 years. This tendency also appears in the δ18O profile of the same ice core.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1990

Introduction

Air bubbles in polar glacier ice are formed by the densification process of firn into ice. Composition of air in ice is essentially the same as that of the atmosphere at the time of bubble formations; thus the analyses of air in ice allow us to investigate the past atmospheric conditions (Reference Scholander, Hemmingsen, Coachman and NuttScholander and others, 1961; Reference Matsuo and MiyakeMatsuo and Miyake, 1966; Reference Stauffer and BernerStauffer and Berner, 1978; Oeschger and others, 1982; Reference Barnola, Raynaud, Korotkevich and LoriusBarnola and others, 1987; Reference Etheridge, Pearman and de SilvaEtheridge and others, 1988).

Total gas content is defined as a volume of air involved in unit mass of ice measured under conditions of standard temperature and pressure (0°C and 1 atmosphere). Knowing the temperature at a point where the air is trapped within the ice, we can estimate the air pressure at the point from total gas content in ice. Hence, the total gas content is considered as an indicator of elevation at bubble formation (i.e. elevation at pore close-off), as stated by Lorious and others (1968). Raynaud and Lebel (Reference Nishio, Mae, Ohmae, Nakawo, Takahashi and Kawada1979) primarily presented an empirical relation between the total gas content and the elevation of sampling site by using ice cores drilled at different altitudes in the Antarctic and Greenland ice sheets. This relationship has been applied to estimate past ice-sheet elevation (Raynaud and Whillans Reference Raynaud and Whillans1982; Reference Jenssen and RadokJenssen and Radok, 1982; Reference Herron and LangwayHerron and Langway, 1987).

The ice-core drilling project in Mizuho Plateau, Antarctica, has been carried out by the Japanese Antarctic Research Expeditions (JARE) since 1971. More than ten shallow ice cores were recovered in this area, and the JARE-25 accomplished the drilling of an ice core to a depth of 700.56 m at Mizuho Station in August 1984.

The present paper gives an empirical relation between total gas content and the elevation at pore close-off using shallow ice cores from Mizuho Plateau. The derived relation was applied to the vertical profile of total gas content in a 700 m deep ice core of Mizuho Station to estimate a past elevation of the ice sheet around the area.

Experimental Procedure

Total gas contents were measured by both “melting method” and “dry extraction method”. The procedure of melting method was described in detail by Langway (Reference Langway1958). The outline of this method is: air in an ice sample is extracted by melting the sample in a liquid and is introduced into a gas burette; then the total volume of air, and temperature and atmospheric pressure, are measured. This method is easy to perform, and is suited for relatively large samples from 20 to 200 g.

As a liquid, we employed a water- and air-saturated dimethyl silicon oil (1.5 × 10−6 m2 s−1), not conventionally-used kerosene, because of its low vapour pressure, slight solubility for water and pure composition. We also measured the dissolved-air volume in the melt water. This procedure enables us to reduce the experimental error to about 1% for 80 g ice samples (Kameda and others, 1989).

The dry extraction method procedure was described in detail by Kameda and others (1989) in detail. The outline of this is: air in an ice sample is extracted by crushing the sample in a vacuum vessel, then the pressure and the temperature of released air are measured. Corrections were made for the amount of vapour pressure in air and also for the loss of air due to cutting air bubbles on the sample surface. The whole procedure was carried out in a cold room (about −20°C). For small samples (namely less than 20 g), this method is more accurate than the melting method. Since there existed numerous fine cracks at intervals of 10 mm or less in the 700 m deep core of Mizuho Station (Mizuho core), only small pieces of ice samples of about 5 to 10 mm thick could be obtained below 120 m in depth and were subjected to measurements with the dry extraction method. Overall error for this method was estimated as less than 1% (Kameda and others, 1989).

Fig. 1. Results of total gas contents from ice core from station Gl5 measured by melting method (triangles) and dry extraction method (solid circles).

To examine systematic differences between the results by melting method and dry extraction method, total gas contents of ice core from station Gl5 were measured in detail from 89 to 100.9 m in depth, and are shown in Figure 1. Although obtained data scattered considerably, mainly due to the seasonal and year-to-year variations in total gas content, mean values of total gas contents by these two methods are 89.9 × 10−6 m3 kg−1 and 89.6 × 10−6 m3 kg−1, respectively. The difference between these was within the error of the two methods. Consequently it may be noted that there are no systematic differences between these two procedures, and a mean value of several samples should be used at each station.

Relation Between Total Gas Content and Elevation at Pore Close-Off in Mizuho Plateau

Total gas contents were measured for seven shallow ice cores from different altitudes in Mizuho Plateau, East Antarctica. A map showing the drilling sites is given in Figure 2; locations, surface elevations (Reference Nakawo, Narita and IsobeNakawo and others, 1984; Reference NishioNishio, 1984; Reference Ageta, Kikuchi, Kamiyama and OkuhiraAgeta and others, 1987) and mean annual air temperatures estimated from 10 m snow temperatures are listed in Table I.

Total volume of air bubbles within a unit mass of core sample was measured from the depth of 50 to 90 m at station G15. The variation in bubble volumes with the bulk density of firn or ice is shown in Figure 3. It is clearly seen that the bubble volume increases gradually with increase in density from 720 to about 830 kg m−3 from which it then decreases with density. The former increasing tendency is considered to result from increase in the number of air bubbles which had closed off, while the latter is caused by the compaction of air bubbles in ice. Hence, we regard the density of 830 kg m−3, where the volume of air bubbles indicated a maximum value, as the bulk density at pore close-off, and we assume it holds all over the Mizuho Plateau area. Elevation at pore close-off in each site was then determined by subtracting, from the surface elevation, the depth of the critical density estimated from a density profile of firn and ice, also shown in Table I. All these ice cores were recovered from 1979 to 1986.

Fig. 2. Locations of ice coring sites in Mizuho Plateau, East Antarctica.

To establish a relation between total gas content and elevation at pore close-off, samples subjected to measurements were those from 20 to 50 m below the transition from firn to ice. Such ice cores contained almost spherical air bubbles, which favoured corrections for air loss due to cutting the sample surface (Reference Nakawo and NaritaNakawo and Narita, 1985), and were also considered to be at almost the same temperature and pressure as those at the level of pore close-off. The numbers of measurements of total gas content for each site were between two at station S48 and 16 at station G15, as shown in Table I. At each site except S48, total gas contents were measured continuously along more than 0.25 m of ice core, which is equivalent to snow accumulation of 1–3 years.

Mean value of total gas content R (m3 kg−1) obtained for each site is shown in Table I, and plotted against the elevation at pore close-off E (m) in Figure 4. Values of R show a tendency to decrease almost linearly with E, except for the data of S25 where elevation is lowest among the seven stations. Site S25 is characterized by relatively high accumulation rate, i.e. about 0.25 m a−1 in water equivalent, and high temperature, i.e. −17.96°C at 10 m in depth (personal communication from K. Satow), compared with other sites Mizuho Plateau. These regional characteristics of depositional regime should cause changes in the densification process of firn, so that the density at pore close-off in S25 may differ from that in inland areas. Hence, excluding the point of S25, the following empirical relation was obtained from 55 observations of total gas contents in Mizuho Plateau:

Table I. Locations and Elevations of Drilling Sites, and Mean Total Gas Contents in Ice Cores

Fig. 3. Total volume of air bubbles in a unit mass of firn or ice core versus the bulk density at station G15.

(1)

The result of Raynaud and Lebel (Reference Raynaud and Lebel1979) is also shown with open squares and a broken line in Figure 4. The difference between these two lines is considered due to the difference in climate of the regions; accumulation rate may be the key factor to it.

Fig. 4. Mean total gas content in ice core from 20 to 50 m below the firn–ice transition versus the elevation (a.s.l.) at pore close-off in Mizuho Plateau. Maximum and minimum values of total gas content are indicated by bars. Open squares and a broken line are the results of Raynaud and Lebel (Reference Raynaud and Lebel1979), in which Antarctic and Greenland ice cores are used.

Fig. 5. Vertical profile of total gas content and δ18O (Reference Higashi, Nakawo, Narita, Fujii, Nishio and WatanabeHigashi and others, 1988) in the 700 m deep ice core at Mizuho Station. A solid line of total gas content is fitted to the measured data (solid circles and triangles). The broken line indicates the profile of original elevation at which pore close-off took place, calculated under a steady-state assumption. Along the lower abiscissa the estimated elevation at pore close-off is shown.

Past Ice-Sheet Elevation Around Mizuho Station

The vertical profile of total gas content plotted against depth at Mizuho Station is shown in the left column of Figure 5 and Table II. Circles and triangles indicate values measured by dry extraction method and melting method, respectively. The two open circles in a shallow layer are only for reference, since the air bubbles were not spherical and corrections for surface effect were insufficient.

Values of total gas content scattered considerably. The width of scattering is about 10 × 10−6 m3 kg−1 at maximum, which is almost the same as a seasonal variation in total gas content reported in Raynaud and Lebel (Reference Raynaud and Lebel1979). As for a general trend of change in total gas content R, we can see a gradual increase in R from 600 to about 200 m in depth, followed by a steep increase toward the upper layer. By using a least squares regression method, a straight line was fitted to the data in the deep layer and a curve to the data in the shallow layer as shown in Figure 5.

The original elevation at which the deep ice of the Mizuho core was transformed from firn, that is the elevation at pore close-off, was estimated assuming that the relation in Equation (1) was valid for the past period. The elevations are shown on the lower abscissa in Figure 5. It is clear that the deeper ice should indicate the lower value of total gas content since the deeper one has originated at higher elevation of the ice sheet. To discuss the causes for gradual then sudden change in total gas contents, one has to evaluate the effect of down-slope movement of ice.

Calculations of particle path lines from points on the ice-sheet surface to the corresponding depths in the Mizuho core were made by using a three-dimensional steady-state flow model (Reference Takahashi and NakawoTakahashi and Nakawo, 1982). Major input data for calculation are the topography of top and bottom surfaces of the ice sheet (Reference Nakawo, Narita and IsobeNakawo and others, 1984), the width between adjacent flow lines (Reference Naruse and ShimizuNaruse and Shimizu, 1978), and the surface mass balance (Reference Nakawo, Narita and IsobeNakawo and others, 1984) in the up stream area of Mizuho Station. Estimated original elevations at which pores would have closed off under the present flow regime are shown by a broken line in Figure 5 for various depths in the Mizuho core. Also shown on the right ordinate of Figure 5 is an estimated time scale of ice in the Mizuho core determined by Nakawo and others (Reference Nakawo, Ohmae, Nishio and Kameda1989). Usually the age of air in bubbles is different from that of surrounding ice due to air mixing in the firn layer (Reference Schwander, Stauffer and SiggSchwander and others, 1988). For the Mizuho core, however, we may regard the time scale of ice as that of the air in the bubbles, because little mixing may have occurred due to the existence of ice crusts, including multi-layered ice crusts, as shown in the vertical profile of stratigraphy in the Mizuho core (Reference Narita, Watanabe, Satow and OkuhiraNarita and others, 1978).

Table II. Numerical Values of Total Gas Content of Mizuho 700 m Ice Core

In Figure 5, the solid line is smaller than the broken one in the layer deeper than 200 m, by 7.8 × l0−6 m3 kg−1. These values are equivalent to about 350 m in elevation difference, from Equation (1).

Values of total gas contents below 120 m in depth, however, might be smaller than the original values due to the gas release through micro-cracks which cannot be detected in sample preparations, since these samples were selected between visible cracks in the ice core. Hence the difference of about 350 m as shown in Figure 5 may be a maximum value which can be estimated from the result of total gas content profile.

Consequently, if we assume that the present relationship of air temperature and pressure with elevation would be valid through the last few thousand years, thickness of the ice sheet decreased by about 350 m at maximum from about 2000 years BP.

In order to compare the results of total gas content with that of δ18O, the vertical profile of δ18O along the same ice core (Higashi and others, 1989) is shown in the right column of Figure 5. A thin zig-zag line represents every 5 m value and a solid line shows the running mean over ten successive values.

The general trend of increase in δ18O with decrease in depth from 700 to about 200 m is due mostly to the down-slope movement of ice from the high land. However, the remarkable increase in δ18O (about 2.8‰) is recognized from 200 to 100 m in depths, where total gas content data shows also sudden increase. To discuss the cause of this remarkable increase of δ18O, we have to evaluate the effect of down-slope movement of ice as the same way of total gas content analyses.

From calculations of particle path lines described previously, it was estimated that the original elevation of 200 m-deep ice was about 70 m higher than that of 100 m-deep ice. Using the gradient of δ18O of deposited snow against 10 m snow temperature, 0.91‰/°C (by personal communication from Y. Fujii), and the mean lapse rate of 10 m snow temperature in Mizuho Plateau (−1.25°C/100 m) derived by Satow and Kikuchi (1989), it was estimated that the difference of 70 m in surface elevation results in a difference of 0.8‰ in δ18O between 200 and 100 m-deep ice.

On the other hand, to estimate a relative warming at Mizuho Station caused by 350 m thinning of the ice sheet, an aerological temperature gradient above the inversion layer at Mizuho Station (−0.5°C/100 m; figure 3.8(a) in National Institute of Polar Research, 1988) can be used. It is, therefore, estimated that the 350 m decrease in ice-sheet elevation may caused by a relative warming about 1.8°C of mean annual temperature. This values also corresponds with an increase of about 1.6‰ in δ18O.

Consequently, most of the increase in δ18O from 200 to 100 m (2.8‰) was explained by relative warming due to thinning of the ice sheet (1.6‰), and by the down-slope movement of the ice (0.8‰). The remainder (about 0.4‰) of δ18O may correspond to climatic warming around Mizuho Station. However, the contribution of this climatic warming (about 0.5°C) may be slightly underestimated since 350 m thinning of the ice sheet is possibly an overestimate, for reason described previously.

Concluding Remarks

The vertical profile of total gas content of the Mizuho 700 m ice core showed a sudden increase from 200 to 100 m in depth. On the basis of a relation between total gas content and elevation at pore close-off derived at Mizuho Plateau, and also of the ice-sheet flow calculations under steady-state assumptions, this increase was explained by a 350 m decrease of the ice-sheet elevation at maximum from 2000 years B.P. The thinning tendency was recognized also in the δ18O profile of the same ice core from 200 to 100 m in depth. It was found that this increase in δ18O was mainly attributed to a relative warming caused by the 350 m decrease in the ice-sheet elevation.

The result of a thinning ice sheet from 2000 years B.P. may confirm that the thinning phenomenon, by a rate from 0.5 m a−1 to about 1 m a−1, directly measured with a short time interval in Mizuho Plateau (Reference NaruseNaruse, 1979; Nishio and others, 1989), would be a long term trend of ice-sheet variation.

Although these considerations may be somewhat speculative at this stage, the results given in this paper should provide an insight to ice-sheet variations, until further analyses have been carried out.

Acknowledgements

This work is a contribution from the Glaciological Research Program in east Dronning Maud Land, Antarctica. We thank the members of JARE-20, 21, 24, 25 and 26, in particular those who made enormous efforts to recover ice cores. We are indebted also to Professor A. Higashi of International Christian University for his generous supervision of the Glaciological Research Program and to Dr Y. Fujii of the National Institute of Polar Research for discussion about δ18O data in Mizuho Plateau.

References

Ageta, Y, Kikuchi, T, Kamiyama, K and Okuhira, F 1987. Glaciological research program in east Queen Maud Land, East Antarctica. Part 5, 1985. Tokyo, National Institute of Polar Research. (JARE Data Reports 125 (Glaciology 14).)Google Scholar
Barnola, J.M, Raynaud, D, Korotkevich, Y.S and Lorius, C 1987. Vostok ice core provides 160,000-year record of atmospheric CO2. Nature, 329(6138), 408414.CrossRefGoogle Scholar
Etheridge, D.M, Pearman, G.I and de Silva, F 1988. Atmospheric trace-gas variations as revealed by air trapped in an ice core from Law Dome, Antarctica. Ann. Glaciol., 10, 2833.CrossRefGoogle Scholar
Herron, S.L and Langway, C.C jr. 1987. Derivation of paleoelevations from total air content of two deep Greenland ice cores. International Association of Hydrological Sciences Publication 170 (Symposium at Vancouver 1987 — The Physical Basis of Ice Sheet Modelling), 283295.Google Scholar
Higashi, Α, Nakawo, M, Narita, H, Fujii, Y, Nishio, F and Watanabe, O 1988. Preliminary results of analyses of 700 m ice cores retrieved at Mizuho Station, Antarctica. Ann. Glaciol., 10, 5256.CrossRefGoogle Scholar
Jenssen, D and Radok, U 1982. On the joint interpretation of total gas contents and stable isotope ratios in ice cores. Ann. Glaciol., 3, 152155.CrossRefGoogle Scholar
Kameda, T, Nakawo, M, Nagoshi, M and Mae, S In press. Measurements of total gas content of an ice core from Mizuho Station, Antarctica. Proceedings of the NIPR Symposium on Polar Meteorology and Glaciology 3.Google Scholar
Langway, C.C jr. 1958. Bubble pressures in Greenland glacier ice. International Association of Scientific Hydrology Publication 47 (Symposium at Chamonix 1958 — Physics of the Movement of the Ice), 336349.Google Scholar
Lorius, C, Raynaud, D and Dolle, L 1968. Densité de la glace et étude des gaz en profondeur dans un glacier antarctique. Tellus, 20(3), 449460.Google Scholar
Matsuo, S and Miyake, Y 1966. Gas composition in ice samples from Antarctica. J. Geophys. Res., 71(22), 52355241.CrossRefGoogle Scholar
Nakawo, M and Narita, H 1985. Density profile of a 413.5 m deep fresh core recovered at Mizuho Station, East Antarctica. Mem. Natl. Inst. Polar Res. Spec. Issue 39, 141156.Google Scholar
Nakawo, M, Narita, H and Isobe, T 1984. Glaciological research program in east Queen Maud Land, East Antarctica. Part 2. 1983. Tokyo, National Institute of Polar Research. (JARE Data Reports 96 (Glaciology 11).)Google Scholar
Nakawo, M, Ohmae, H, Nishio, F and Kameda, T 1989. Dating the Mizuho 700-m core from core ice fabric data. Proceedings of the NIPR Symposium on Polar Meteorology and Glaciology 2, 105110.Google Scholar
Narita, H, Watanabe, O, Satow, K and Okuhira, F 1978. Compiled stratigraphic data from cores drilled at Mizuho Station. Mem. Natl. Inst. Polar Res. Spec. Issue 10, 132135.Google Scholar
Naruse, R 1979. Thinning of the ice sheet in Mizuho Plateau, East Antarctica. J. Glaciol., 24(90), 4552.CrossRefGoogle Scholar
Naruse, R and Shimizu, H 1978. Flow line of the ice sheet over Mizuho Plateau. Mem. Natl. Inst. Polar Res. Spec. Issue 7, 227234.Google Scholar
National Institute of Polar Research. 1988. In Kawaguchi, S and Kanzawa, H eds. Science in Antarctica. 3. Meteorology. Tokyo, National Institute of Polar Research. [In Japanese.]Google Scholar
Nishio, F ed. 1984. Glaciological research program in east Queen Maud Land, East Antarctica. Part 1, 1982–1983. Tokyo, National Institute of Polar Research. (JARE Data Reports 94 (Glaciology 10).)Google Scholar
Nishio, F, Mae, S, Ohmae, H, Nakawo, M, Takahashi, S and Kawada, Κ 1984. Dynamical behaviour of the ice sheet on Mizuho Plateau, East Antarctica. (Abstract.) Ann. Glaciol., 12, 205.CrossRefGoogle Scholar
Oeschger, Η, Stauffer, B, Neftel, A, Schwander, J and Zumbrunn, R 1982. Atmospheric CO2 content in the past deduced from ice-core analyses. Ann. Glaciol., 3 227232.CrossRefGoogle Scholar
Raynaud, D and Lebel, B 1979. Total gas content and surface elevation of polar ice sheet. Nature, 281(5729), 289291.CrossRefGoogle Scholar
Raynaud, D and Whillans, I.M 1982. Air content of the Byrd core and past changes in the West Antarctic ice sheet. Ann. Glaciol., 3, 269273.CrossRefGoogle Scholar
Satow, K and Kikuchi, T In press. The 10 m snow temperature in the ice sheet. In Higashi, A ed. Antarctica: East Queen Maud Land. Enderby Land; Glaciological Folio. Tokyo. National Institute of Polar Research.Google Scholar
Scholander, P.F, Hemmingsen, E.A, Coachman, L.K and Nutt, D.C 1961. Composition of gas bubbles in Greenland icebergs. J. Glaciol., 3(29), 813822.CrossRefGoogle Scholar
Schwander, J, Stauffer, B and Sigg, A 1988. Air mixing in firn and the age of the air at pore close-off. Ann. Glaciol., 10, 141145.CrossRefGoogle Scholar
Stauffer, B and Berner, W 1978. CO2 in natural ice. J. Glaciol., 21(85), 291300.CrossRefGoogle Scholar
Takahashi, S and Nakawo, M 1982. An estimation of the steady flow field along a flow line of Shirase Glacier, Antarctica. Seppyo, 44(4), 189195. [In Japanese.]Google Scholar
Figure 0

Fig. 1. Results of total gas contents from ice core from station Gl5 measured by melting method (triangles) and dry extraction method (solid circles).

Figure 1

Fig. 2. Locations of ice coring sites in Mizuho Plateau, East Antarctica.

Figure 2

Table I. Locations and Elevations of Drilling Sites, and Mean Total Gas Contents in Ice Cores

Figure 3

Fig. 3. Total volume of air bubbles in a unit mass of firn or ice core versus the bulk density at station G15.

Figure 4

Fig. 4. Mean total gas content in ice core from 20 to 50 m below the firn–ice transition versus the elevation (a.s.l.) at pore close-off in Mizuho Plateau. Maximum and minimum values of total gas content are indicated by bars. Open squares and a broken line are the results of Raynaud and Lebel (1979), in which Antarctic and Greenland ice cores are used.

Figure 5

Fig. 5. Vertical profile of total gas content and δ18O (Higashi and others, 1988) in the 700 m deep ice core at Mizuho Station. A solid line of total gas content is fitted to the measured data (solid circles and triangles). The broken line indicates the profile of original elevation at which pore close-off took place, calculated under a steady-state assumption. Along the lower abiscissa the estimated elevation at pore close-off is shown.

Figure 6

Table II. Numerical Values of Total Gas Content of Mizuho 700 m Ice Core