Hostname: page-component-cd9895bd7-8ctnn Total loading time: 0 Render date: 2024-12-22T08:58:46.008Z Has data issue: false hasContentIssue false

Statistical-Dynamical Model of Accumulation on the Greenland Ice Sheet

Published online by Cambridge University Press:  20 January 2017

Richard A. Keen*
Affiliation:
Cooperative Institute for Research in Environmental Sciences, Campus Box 449, University of Colorado, Boulder, Colorado 80309, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

A relatively simple, three-parameter model is used to simulate the annual precipitation (accumulation) distribution for the Greenland ice sheet and surrounding regions. The three parameters are (1) the flux of relative vorticity at the 500 mbar level (a measure of cyclonic activity), (2) atmospheric moisture content, and (3) surface terrain. The climatological (1946–79) precipitation distribution predicted by the model displays major features of the observed distribution derived from pit studies. However, the model suggests that, due to changes in storm tracks during this period of 33 a, accumulation distribution maps based on pit studies for varying periods of record may not be representative of a true mean for a uniform period of record. The model is then applied to reconstructed ice-age conditions. Compared to present conditions, accumulation reductions of 60¾ or more are indicated for much of the southern half of Greenland; only slight reductions are noted for northern Greenland.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1984

1 Introduction

As is the case for precipitation anywhere in the world, the spatial distribution of snow accumulation on the Greenland ice sheet may vary greatly from year to year. This variability is largely due to the changing tracks and intensities of cyclones affecting the region. Persistent anomalies of cyclone activity may affect the representativeness of accumulation distributions based on short-term stratigraphies; this problem is magnified when recent stratigraphies are applied to extremely anomalous epochs such as the last ice age. These problems are addressed in this study through development and application of a simple statistical-dynamical model that relates annual accumulation to an index of cyclone activity. The model is calibrated against recent observed annual accumulations, and is then used to reconstruct the accumulation distribution across Greenland during the ice age.

2 Model Description

The vorticity flux index

Surface cyclones are reflected in the middle and upper levels of the troposphere by the “baroclinic short waves”. These short waves can be identified on upper-level charts as regions of positive relative vorticity, and the magnitude of the vorticity gives a measure of the intensity of the system. The vorticity field at the 500 mbar level is particularly useful in forecasting and analysis, because it is near the level of non-divergence where the short waves move with the larger-scale (“long wave”) wind field (Reference FjørtoftFjørtoft 1952) and where vertical velocities (and therefore precipitation generation) are at a maximum (Reference HaltinerHaltiner 1971).

For this climatological analysis, short waves are treated as a flux quantity by multiplying the short-wave positive relative vorticity field by the long-wave vector wind field. The resulting quantity, called the vorticity-flux index , is defined by the equation

(1)

The integration is approximated by a summation of daily values computed from National Meteorological Center (NMC) gridded 500 mbar height fields. Details of the computational techniques may be found in Reference KeenKeen (1980).

The mean over the Greenland region for 34 years (1946–80) is shown in Figure 1. It can be seen that upper-level systems approach Greenland from the west and south-west, on the average, and that the incidence of storms is greatest along the south-east coast.

Fig. 1 27-year vector mean vorticity flux index for the Greenland region. The lengths of the arrows are proportional to the magnitude of the ; the horizontal scale line indicates a magnitude of 10−4 m sec−2. Bases of the arrows are located at the computational grid points.

VFI and precipitation

Accumulated precipitation at a point is related to the time and height (pressure)-integrated upward vertical velocity by the equation (Reference HaltinerHaltiner 1971)

(2)

qsis the saturation specific humidity and f(T) is a thermodynamic function of temperature, ð is an on-off function set equal to one when both ω is upward and q equals qs, and otherwise equals zero.

In the equivalent barotropic model of the atmosphere (Reference HaltinerHaltiner 1971) the vertical velocity ω is related to the relative vorticity ζ by the equation

(3)

derived originally by Reference Jenssen and RadokJenssen and Radok (1964). Overbars indicate vertical averaging of the quantities between the surface and the level of computation; p0 and ω0 are surface values. B(p) and C(p) are slowly varying functions of pressure describing the vertical wind shear; climatological values of these functions optimize to yield largest vertical velocities in the mid-troposphere (near 500 mbar), where this model works best. The two terms on the right-hand side of Equation (3) represent, respectively, the dynamic and orographic components of the vertical velocity. The dynamic component is primarily induced by the advection of vorticity by synoptic-scale systems. Over a time interval ∆t, the integrated amount of upward vertical motion of the air at a point is therefore determined mainly by the integrated positive vorticity advection (PVA). The magnitude of the is roughly proportional to the total PVA, if the short waves are sinusoidal and their aggregate mean length (and thus the scale of the operator) remains approximately constant. This assumption eliminates much of the “noise” that would result from taking an additional derivative of the gridded analyses. To the extent that the is representative of the horizontal flow direction during precipitation, the orographic component of ω may be parameterized by , where h is the surface terrain smoothed to the 400 km resolution of the NMC grid. Most orographic uplift occurs below the 500 mbar level, but it will be seen that this parameterization is effective in reproducing the orographic precipitation distribution across Greenland.

The concepts in Equations (1), (2) and (3) are combined to give a parameterization of the essential mechanisms of precipitation generation

(4)

In this model the mean saturation specific humidity qs is presumed to be proportional to the climatological mean specific humidity q. Because of the large seasonal variation of q (by a factor of three between summer and winter), vorticity fluxes are calculated separately for each of the four three-month seasons, and are used in Equation (4) with climatological values of q (derived from 500 mbar dewpoints given by Reference Crutcher and MeserveCrutcher and Meserve (1970)) for the respective season. Equation (4) is then summed over the four seasons to produce a modeled precipitation for the accumulation year from September through August.

The empirical coefficients A0 and A1 are derived from multiple linear regression of the modeled orographic and dynamic components on observed annual accumulations. For this purpose, accumulations for the years from 1946–47 through 1954–55 derived from pit studies in central and north-western Greenland (Reference BensonBenson 1962) are smoothed to the 400 km grid resolution of the model. In this manner three gridded accumulation points (see Fig.2) are available for the nine years of record; the correlation with the model yields r = 0.61 (n = 27).

Fig. 2 Contribution to mean annual precipitation predicted by the dynamic component of the accumulation model. Major contour interval: 20 cm water equivalent. A, B, and C locate the three grid points used to calibrate the model.

Comparison with Observations

The 33-year (1946–47 through 1978–79) mean annual precipitation distribution over the Greenland region predicted by the dynamic and orographic components of the model are shown in Figures 2 and 3, respectively. Between north-western and extreme southern Greenland the dynamic component can be seen to range from 15 to 200 cm. This large variation is due to the much greater storm activity (Fig.1) and higher moisture content of the air over south-east Greenland. The orographic component enhances the accumulation west of the ice divide and reduces accumulations to the east; combined with the dynamic component the total accumulation distribution (Fig.4) shows major observed features as the dry “snow shadow” region of north-eastern Greenland and the wet zone on the west slope of the ice sheet (Fig.5). Taking the observed mean annual accumulations at face value (but keeping in mind the caveat in the next paragraph), the model apparently overestimates accumulation over southern Greenland by a factor of two and underestimates the magnitude of the “snow shadow”, suggesting that the derived A0 is too large and A1 is too small. This problem may arise from the location of the three calibration grid points in areas of weak upslope near the crest of the ice sheet. It is also possible that clear-sky precipitation (e.g. Reference Radok, Lile and BusingerRadok and Lile 1977), which is a process not included in the model, may be a significant contributor to total accumulation over northern Greenland that leads to overestimation of A0. The calibration could be improved by incorporating observed annual accumulations from other regions of the ice sheet.

Fig. 3 Contribution to mean annual precipitation predicted by the orographic component of the accumulation model. Contour interval: 5 cm water equivalent; zero contour omitted for clarity.

Fig. 4 Mean annual precipitation distribution predicted by combined dynamic and orographic components of the accumulation model. Major contour interval: 20 cm water equivalent.

Fig. 5 Gridpolnt values of mean annual accumulation (cm water equivalent), spatially smoothed to the 400 km grid resolution from observations published in Reference BensonBenson (1962) and Reference Barry, Kiladis, Radok, Barry, Jenssen, Keen, Kiladis and MclnnesBarry and Kiladis (1982).

The temporal variability of calculated precipitation for 1946–47 through 1978–79 is shown in Figure 6 for the three grid points used in the calibration of the model; for comparison the observed values for 1946–47 through 1954–55 are also shown. The model tends to exaggerate the interannual precipitation variability, possibly because the variability of actual precipitation is reduced by heavier clear-sky precipitation during seasons with less cyclonic activity. An outstanding feature of these time series is the decrease in calculated accumulation at the northern grid point (A in Fig.2) after 1956, and a corresponding increase at the central Greenland grid point C. An implication of this feature is that there are long-term changes in the accumulation distribution across Greenland, due to shifts in storm tracks. Therefore, mean accumulation distribution maps derived from pit studies for varying periods of record may not be representative of a true mean for a uniform period of record.

Fig. 6 Comparison of model precipitation (solid line) with observed accumulation (dashed line) for the three grid points used to calibrate the model (see Fig.2).

The coarseness of the model grid is compatible with the spacing of the meteorological observing network upon which the vorticity flux calculations are based. Finer-resolution versions of the model, using smaller-scale orography, have been tested against unsmoothed pit accumulations and coastal station precipitation totals, but the correlations are generally not as close as for the larger-scale model. To some extent this reflects the absence of small (<400 km) scale features in the 500 mbar data set. However, correlations between pits and/or stations separated by only a few kilometers are often rather poor, and extremely local factors such as drifting and snow-gauge exposure may also play a role (Reference Barry, Kiladis, Radok, Barry, Jenssen, Keen, Kiladis and MclnnesBarry and Kiladis 1982).

3 An Application of the Accumulation Mddel To Ice-Age Conditions

The model described in the previous section can be applied to any epoch for which the input parameters, vorticity flux, atmospheric moisture, and surface terrain, are known or can be estimated. An epoch of great interest is the last glacial maximum (ca. 18 ka BP), during which rates and distribution of precipitation probably differed greatly from those of the present. The model has been applied to ice-age conditions, with the three input fields estimated as follows:

Vorticity flux

Large-scale circulation patterns during the last glacial maximum can be reconstructed approximately by presuming that surface temperature anomalies (derived by a variety of methods) extend upwards through the troposphere. Reference Lamb and WoodroffeLamb and Woodroffe (1970) have done such a reconstruction for the northern hemisphere 1 000 to 500 mbar thickness field, and through dynamic considerations estimated the surface pressure field for summers and winters during the Würm ice maximum (ca. 18 ka BP). Combining the surface pressure and thickness fields yields reconstructed 500 mbar height fields; the departure of the reconstructed annual mean height field from present climatology is shown in Figure 7. The major features of this departure field are the positive (+100 m) anomaly over southern Greenland and the larger (−400 m) negative anomaly around the Great Lakes. The selection of this long-wave circulation pattern for the ice-age reconstruction is supported by varied observational and theoretical evidence reviewed by Reference CrowleyCrowley (1984).

Fig. 7 Departure of the reconstructed annual mean 500 mbar height field for 18 ka BP from present climatology. Contour interval: 10 dm.

Vorticity fluxes cannot be computed directly from a time-averaged 500 mbar height field, and so are derived by an analog method. Based on the similarity between observed seasonal 500 mbar height departures and the ice-age reconstructions, the analog seasons chosen are the fall and winter of 1976 to 1977, and_ the spring and summer of 1958. The Greenland (positive) and Great Lakes (negative) anomaly centers in the . analog seasons are about a third the magnitude of those during the ice age (Fig.7). A first attempt at reconstructing vorticity fluxes, therefore, simply tripled the vorticity flux anomalies (i.e. the vector departure from the 34-year average) for each of the analog seasons. However, this created an improbable reversal of the direction of motion of the storm track off south-east Greenland. Multiplying the analog season vorticity flux anomalies by 1.7 yielded reconstructed ice-age vorticity flux fields (Fig.8) that are consistent with the reconstructed mean 500 mbar fields. The outstanding features of this reconstruction are a marked weakening of the storm track south and east of Greenland, a relative strengthening of the storm track over Baffin Island and southern Baffin Bay, and an absolute increase in this storm track over northern Baffin Bay, northern Greenland, and into the polar basin. Cyclone activity over the central Greenland ice sheet is reduced to nil.

Fig. 8 Reconstructed annual mean for 18 ka BP. Meaning of arrows as in Figure 1.

Atmospheric moisture content

Correlation of climatological mean dewpoints and thicknesses for various Greenland stations and seasons gives a dewpoint change of 3°C per 100 m. This value is applied to Lamb and Woodroffe’s ice-age thickness departure fields to give summer and winter ice-age dewpoint departures (spring and autumn are interpolated) from present climatology (e.g. Fig.9). The ice-age moisture content is computed from these adjusted dewpoints.

Fig. 9 Departure of reconstructed annual mean 500 mbar dewpoint temperatures for 18 ka BP from present climatology. Contour interval: 2°C.

Surface terrain

Surface terrain during the ice age (Fig.10) is assumed to be unchanged from the present for Greenland, Ellesmere Island, and Iceland. A continental ice sheet approximating that modeled by Reference Budd and SmithBudd and Smith (1981) for 18 ka BP is placed over the Canadian archipelago, with heights rising from sea level along the coast of Baffin Island to nearly 5 000 m at the south-west corner of the grid (northern Hudson Bay).

Fig. 10 Smoothed terrain used in the accumulation model; 1 km contour interval. Estimated conditions for 18 ka BP are shown; these differ from present terrain only in the south-western corner of the map.

Results

The reconstructed annual precipitation rate for the Greenland region is shown in Figure 11, and is shown as a percentage of the 1946–79 average model precipitation in Figure 12. Because of the calibration problems noted earlier, the accumulation percentages in Figure 12 should be considered more reliable than the calculated accumulations in Figure 11. Decreases of 60% during the ice age are widespread over the southern half of Greenland, and the 7 cm accumulation rate at the center of the ice sheet is but 30% of the present model value. In general, about a third of the decreases are due to reduced atmospheric moisture content; most of the rest is due to the weaker storm activity. The 70% decrease over south-eastern Baffin Island is largely due to orographic descent of storm systems from the continental ice sheet. While still somewhat below present values, ice-age accumulations over northern Greenland show an increase relative to the southern two-thirds of the island. In extreme northern Greenland the strengthened storm track just offsets the moisture decrease to produce a small area of slightly increased precipitation.

Fig. 11 Reconstructed annual precipitation distribution predicted by the accumulation model for 18 ka BP. Major contour interval: 20 cm water equivalent.

Fig. 12 Reconstructed annual precipitation for 18 ka BP, expressed as a percentage of 1946–79 average model precipitation for each grid point. Major contour interval: 20%.

4 Conclusion

A relatively simple, three-parameter model is used to simulate the precipitation (accumulation) rate distribution for the Greenland ice sheet and surrounding regions. The model is calibrated against time series of observed pit accumulations for three grid points in central and north-western Greenland. The climatological (1946–79) precipitation distribution predicted by the model displays major features of the observed distribution derived from pit studies. However, the model suggests that due to changes in storm tracks during the 33-year period, accumulation distribution maps based on pit studies for varying periods of record may not be representative of a true mean for a uniform period of record.

The model is then applied to reconstructed ice-age conditions. While the results of the ice-age model are difficult to verify, this particular application suggests the potential uses of this, or any similar, model. The model can be applied to any large ice sheet for any period of time for which the input parameters are known or can be estimated. Furthermore, the input parameters can be time-dependent, and the precipitation model may be combined with an ice-sheet model for an interactive atmospheric-cryospheric simulation of ice-sheet development.

This research was supported by the National Science Foundation and the National Oceanic and Atmospheric Administration; computing facilities were provided by the National Center for Atmospheric Research. The author is grateful to George Kiladis for his assistance in developing the accumulation model, and to Uwe Radok and Roger Barry for their helpful comments.

References

Barry, R G, Kiladis, G N 1982 Climatic characteristics of Greenland. In Radok, U, Barry, R G, Jenssen, D, Keen, R A, Kiladis, G N, Mclnnes, B (eds) Climatic and physical characteristics of the Greenland ice sheet. Boulder, CO, University of Colorado. Cooperative Institute for Research in Environmental Sciences: 733 Google Scholar
Benson, C S 1962 Stratigraphie studies in the snow and firn of the Greenland ice sheet. SIPRE Research Report 70 Google Scholar
Budd, W F, Smith, I N 1981 The growth and retreat of ice sheets in response to orbital radiation changes. International Association of Hydrological Sciences Publication 131 (Symposium at Canberra 1979 – Sea Level, Ice and Climatic Change): 369409 Google Scholar
Crowley, T J 1984 Atmospheric circulation patterns during glacial inception: a possible candidate. Quaternary Research 21(1): 105110 CrossRefGoogle Scholar
Crutcher, H L, Meserve, J M 1970 Selected level heights, temperatures, and dew points for the northern hemisphere. Washington, DC, US Naval Weather Service Command (NAVAIR Report 50–1C–52)Google Scholar
Fjørtoft, R 1952 On a numerical method of integrating the barotropic vorticity equation. Tellus 4(3): 179194 Google Scholar
Haltiner, G 1971 Numerical leather prediction. New York, John Wiley Google Scholar
Jenssen, D, Radok, U 1964 Diabatic heating and cooling in the equivalent-barotropic atmosphere. nature 202: 11041106 Google Scholar
Keen, R A 1980 Temperature and circulation anomalies in the eastern Canadian Arctic, 1946–1976. university of Colorado. Institute of Arctic and Alpine Research. Occasional Paper 34 Google Scholar
Lamb, H H, Woodroffe, A 1970 Atmospheric circulation during the last ice age. Quaternary Research (1): 2958 CrossRefGoogle Scholar
Radok, U, Lile, R C 1977 A year of snow accumulation at Plateau station. In Businger, J A (ed) Mateoro-logical studies at Plateau station, Antarctica. Washington, DC, American Geophysical Union (Antarctic Research Series 25): 1726 CrossRefGoogle Scholar
Figure 0

Fig. 1 27-year vector mean vorticity flux index for the Greenland region. The lengths of the arrows are proportional to the magnitude of the ; the horizontal scale line indicates a magnitude of 10−4 m sec−2. Bases of the arrows are located at the computational grid points.

Figure 1

Fig. 2 Contribution to mean annual precipitation predicted by the dynamic component of the accumulation model. Major contour interval: 20 cm water equivalent. A, B, and C locate the three grid points used to calibrate the model.

Figure 2

Fig. 3 Contribution to mean annual precipitation predicted by the orographic component of the accumulation model. Contour interval: 5 cm water equivalent; zero contour omitted for clarity.

Figure 3

Fig. 4 Mean annual precipitation distribution predicted by combined dynamic and orographic components of the accumulation model. Major contour interval: 20 cm water equivalent.

Figure 4

Fig. 5 Gridpolnt values of mean annual accumulation (cm water equivalent), spatially smoothed to the 400 km grid resolution from observations published in Benson (1962) and Barry and Kiladis (1982).

Figure 5

Fig. 6 Comparison of model precipitation (solid line) with observed accumulation (dashed line) for the three grid points used to calibrate the model (see Fig.2).

Figure 6

Fig. 7 Departure of the reconstructed annual mean 500 mbar height field for 18 ka BP from present climatology. Contour interval: 10 dm.

Figure 7

Fig. 8 Reconstructed annual mean for 18 ka BP. Meaning of arrows as in Figure 1.

Figure 8

Fig. 9 Departure of reconstructed annual mean 500 mbar dewpoint temperatures for 18 ka BP from present climatology. Contour interval: 2°C.

Figure 9

Fig. 10 Smoothed terrain used in the accumulation model; 1 km contour interval. Estimated conditions for 18 ka BP are shown; these differ from present terrain only in the south-western corner of the map.

Figure 10

Fig. 11 Reconstructed annual precipitation distribution predicted by the accumulation model for 18 ka BP. Major contour interval: 20 cm water equivalent.

Figure 11

Fig. 12 Reconstructed annual precipitation for 18 ka BP, expressed as a percentage of 1946–79 average model precipitation for each grid point. Major contour interval: 20%.