Hostname: page-component-586b7cd67f-rcrh6 Total loading time: 0 Render date: 2024-11-26T10:20:28.214Z Has data issue: false hasContentIssue false

Using upper-air conditions to estimate South Cascade Glacier (Washington, U.S.A.) summer balance

Published online by Cambridge University Press:  08 September 2017

Lowell A. Rasmussen
Affiliation:
Department of Earth and Space Sciences, University of Washington, Box 351310, Seattle, Washington 98195-1310, U.S.A. E-mail: [email protected]
Howard B. Conway
Affiliation:
Department of Earth and Space Sciences, University of Washington, Box 351310, Seattle, Washington 98195-1310, U.S.A. E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

A simple model uses once-daily meteorological values in the U.S. National Centers for Environmental Prediction and U.S. National Center for Atmospheric Research (NCEP–NCAR) re-analysis database to estimate summer balance of South Cascade Glacier, Washington, U.S.A., each year over 1959–99. The rms error, 0.30 m w.e. (r2 = 0.71), is comparable to measurement error. The model relates summer balance linearly to temperature T > 0°C at 2000 m andto snow flux at 1650 m, the altitudes in recent years of the equilibrium line and terminus. The snow flux is the product of the humidity and westerly wind component at 850 hPa when temperature T <+2°C at 1650 m. Temperatures are interpolated linearly between the 850 and 700 hPa levels. Both the positive 2000 m temperature and the snow flux are summed from 26 April to 4 October. When the summer estimates are combined with those from a winter balance model using the same database, the rms error in estimating net balance is 0.40 m (r2 = 0.81). The indicated sensitivities of balance to warming of 1°C are −0.51 m for summer and −0.24 m for winter. On the assumption that the total −0.75 m °C−1 sensitivity exists at all altitudes, a warming of only 0.7°C would be sufficient to overcome the 1986–98 average net balance +0.5 m at the top of the glacier.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2003

1. Introduction

Many previous studies of summer balance have used meteorological observations on or adjacent to the glacier itself undertaken specifically for those studies, or from an official weather station at some distance from the glacier, usually at lower altitude. Some apply rigorous physical principles to calculating ablation; others relate it to accumulated positive temperature during the ablation season. Recent examples of these approaches are described in Reference Braithwaite and ZhangBraithwaite and Zhang (2000), Reference OerlemansOerlemans (2000) and Reference Torinesi, Letréguilly and VallaTorinesi and others (2002), which give numerous references to earlier work. Reference Hodge, Trabant, Krimmel, Heinrichs, March and JosbergerHodge and others (1998) examined correlations of South Cascade Glacier (Washington, U.S.A.) seasonal mass balances over 1966–95 with four different measures of large-scale patterns of atmospheric circulation and of sea surface temperature: four measures correlated significantly with winter balance, but only one, sea surface temperature in the equatorial Pacific, correlated with summer balance (r 2 = 0.30).

Recent work has shown that upper-air data can provide a good estimate of conditions on glaciers. For example, temperature on Blue Glacier in the Olympic Mountains of Washington State (Fig. 1) is much better approximated by interpolating from a nearby radiosonde than by applying a lapse rate to measurements at a lowland station (Reference Conway, Rasmussen and HayesConway and others, 1995). A simple flux model that uses upper-air data estimates the winter balance of South Cascade Glacier with error 0.24 m (Reference Rasmussen and ConwayRasmussen and Conway, 2001).

Fig. 1. Northwestern Washington State. The U.S. National Weather Service radiosonde was at Tatoosh Island prior to 1 August 1966 and at Quillayute afterward.The heavy curve is the smoothed coastline.

The principal objective of this study is to investigate the use of upper-air variables routinely measured and archived in the U.S. National Centers for Environmental Prediction and U.S. National Center for Atmospheric Research (NCEP–NCAR) re-analysis database for estimating the summer balance of South Cascade Glacier. Another objective is to combine summer balance b s estimates with those of winter balance b w and estimate net balance from b n = b w + b s, in which b s is defined to be a negative quantity. A third is to demonstrate the effect of a hypothetical warming of the atmosphere on each of the balance components. Here, all balance values are in meters water equivalent, abbreviated to m.

Two summer balance models are presented. The first simply relates b s linearly to daily measurements of upper-air temperature at the equilibrium-line altitude (ELA) of the glacier. A second, which offers a slight improvement, also uses upper-air winds and humidity to estimate the effects of summer snowfall.

2. South Cascade Glacier

South Cascade Glacier (48.36° N, 121.06° W) is a small, north-facing, temperate glacier on the crest of the North Cascade Range of northwestern Washington, 250 km from the Pacific Ocean (Fig. 1). It has an altitude range 1630–2130 m and an area of 2.0 km2. Most precipitation on the glacier falls as snow, with ≈3 m accumulating each winter at the head of the glacier. Annual ablation at the terminus usually exceeds 8 m.

The glacier is not in equilibrium with the recent climate, having retreated and lost much mass over the past century (Reference MillerMiller, 1969; Reference KrimmelKrimmel, 2000). A glacier topography in equilibrium with the 1986–98 mass balance would occupy only about 30% of its 1999 extent (Reference Rasmussen and ConwayRasmussen and Conway, 2001). Mean 1986–98 summer balance varied from about −2.5 m at the head of the glacier to −7.0 m at the terminus (Reference Rasmussen and ConwayRasmussen and Conway, 2001, fig. 4). Over that period the ELA was 1990 m.

3. Data Sources

The U.S. Geological Survey established a research program at South Cascade Glacier in the mid-1950s. Values of winter balance and net balance for each year, both averaged over that year’s glacier topography, have been determined annually since 1958 (Reference KrimmelKrimmel, 2001). Here, revisions of previously published 1986–91 values (Reference KrimmelKrimmel, 2000) are used.

Winter balance at South Cascade Glacier is measured near the end of the accumulation season by probing at numerous points spanning the altitude range of the glacier and converting the snow depth to its w.e. amount by using an altitude variation of density obtained from direct measurements at usually one or two altitudes. Net balance is measured at the end of the ablation season by using the previous summer surface as the reference. In early years net balance was measured at several stakes in holes drilled into the ice, afterward at a smaller number of stakes in positions learned in the early years to be strategic for determining the altitude variation over the whole glacier. In recent years three or four ablation stakes were installed and snow depth was probed at 20–50 sites (Reference KrimmelKrimmel, 2000). Summer balance is calculated as the difference between the net balance and winter balance.

The U.S.National Weather Service operated a radiosonde station with twice-daily soundings from 1948 until 1 August 1966 atTatoosh Island, 275 km west of South Cascade Glacier, and since then at Quillayute, about 50 km south and 10 km east of Tatoosh Island (Fig. 1). Soundings measure temperature, humidity and wind at many levels in the atmosphere. The NCEP–NCAR re-analysis database consists of radiosonde observations, along with surface and satellite observations, systematically merged into a set of internally consistent meteorological values, with 6 hour resolution since 1948, interpolated at several vertical levels to a regular grid spanning the entire Earth (Reference KalnayKalnay and others,1996).

The re-analysis data used here were obtained from NCAR in April 2001. Data include temperature, winds and humidity measured at 1200 UTC (Universal Time Coordinated; 0400 local time) over 1959–99, at 850 hPa (≈1450 m) and 700 hPa (≈3000 m), at a gridpoint at 47.5° N, 125° W, which is 310 km west-southwest of the glacier (Fig. 1). Data from that gridpoint were acquired because it is nearest Quillayute, the radiosonde station from which data were used in earlier work on a closely related subject (Reference Rasmussen, Conway and HayesRasmussen and others, 2001). Horizontal variations of upper-air temperature were scaled from an upper-air atlas (United States Navy, 1966), which showed that the temperature gradient swings from south-southwest in winter to east-southeast in summer and doubles in magnitude. Projection of the gradient on the distance between the glacier and the re-analysis gridpoint yields an adjustment ΔT, which is assumed to be constant over the altitude range of the glacier. Monthly mean 1959–99 temperatures at the gridpoint interpolated at 2000 m between the 850 and 700 hPa levels, along with those estimated for the glacier, are shown in Figure 2.

Fig. 2. Seasonal variation of temperature at South Cascade Glacier and at re-analysis gridpoint 47.5° N, 125° W. The value at the gridpoint is the 1959–99 mean of the temperature interpolated at 2000 m between the 850 and 700 hPa levels in the NCEP–NCAR re-analysis database. Temperature at the glacier is calculated by applying to it the horizontal gradient of mean monthly 850 hPa temperature (United States Navy, 1966).

4. Regional Climatology

Northwestern Washington has mild, wet winters and mild, dry summers. The four wettest months, November–February, account for about 57% of the annual precipitation, and the four driest, June–September, for only about 11%. Nearly all precipitation in this region is produced by extratropical cyclones from the Pacific Ocean, often occurring in rapid succession in winter and only occasionally in summer. Convective storms are uncommon, even in summer. The north–south-trending Cascade Range is an effective climate barrier that delineates the humid marine region to the west from the arid continental region to the east.

5. Estimating South Cascade Glacier Summer Balance

5.1. An ablation formulation, the T model

Of the upper-air variables examined, temperature is the one that correlates best with summer balance (Table 1). In a simple linear relation, balance in year j is estimated from

(1)

in which the constant term γ 1 presumably accounts for radiative melt in the absence (T ≤ 0) of turbulent exchange. Here is the temperature interpolated at 2000 m (the ELA in recent years) between the 850 and 700 hPa levels at the re-analysis gridpoint, adjusted to the location of the glacier by the horizontal gradient (Fig. 2) and averaged from day t 0 through day t 1. Model parameters (t 0, t 1, α 1 and γ 1) were determined empirically to minimize the rms error ∊s. Over the 41 years j = 1959,…,1999, the rms of the residuals is 0.33 m.Where the standard deviation of b s,j is σ = 0.55 m, the coefficient of determination r 2 = 0.65 is calculated from

(2)

Optimum start and end dates are t 0 = 116 (26 April) and t 1 = 277 (4 October), with α 1 = −0.50 and γ 1 = −0.67. Over 1959–99, and . Because α < 0, the correlation between and b s is negative, r = −0.81.

Table 1. Correlations over 1959–99

Forming a regression by summing T ELAT′ on only those days when it exceeds the threshold T ELA > T′, such as in positive-degree-day models that use T′ = 0°C, degraded results (Fig. 3d). Because summer balance at South Cascade Glacier consists of a large negative ablation component and a small positive accumulation component, the inclusion of negative temperatures in forming presumably accounts for summer snowfall events. In addition, T > 0 can sometimes occur low on the glacier, producing melt there, when T ELA < 0.

Fig. 3. Sensitivity of summer balance models to values of model parameters: (a) beginning date t0 for temperature averaging; (b) endingdate t1 for temperatureaveraging; (c) critical direction ϕ′ for calculating f (Equations (3) and (4)); (d) threshold temperature T′ above which to sum TELA in Equation (1). Aseach parameter(a–c)intheT; f model(Equation (5)) is varied, the others are held at their optimum values: t0 = 116 (26 April), t1 = 277 (4 October) and ϕ′ = 271°.The T model (Equation (1)) differed only in having optimum t1 = 279; itsoptimum T′ beingbelow the lowest TELA in the period of record indicates all TELA should be used in formin .

5.2. An ablation–accumulation formulation, the T, f model

Snowfall can be approximated explicitly by taking it to be proportional to the snow flux f when T TERM at the glacier terminus, about 1650 m, is less than +2°C.The snow flux is given by

(3)

in which 0 ≤ RH ≤ 1 is the 850 hPa relative humidity, and U is the component of the wind in the empirically determined critical direction ϕ′.That is,

(4)

where ϕ 850 is the direction and the speed in m s−1 of the 850 hPa wind.The wind, RH, and T TERM are at the re-analysis gridpoint, with T TERM being interpolated between the 850 and 700 hPa levels and adjusted to the glacier by the difference ΔT (Fig. 2). The optimum ϕ′ = 271° is adopted from the results of running a winter balance model (Reference Rasmussen and ConwayRasmussen and Conway, 2001) with 1200 UTC re-analysis data.Tests of the model using either RH or specific humidity (Reference Rasmussen, Conway and HayesRasmussen and others, 2001, table 4) indicate that RH gives better results, presumably because precipitation is a saturation phenomenon. Section 6 of that study also showed that the simple flux model gave better results than a two-dimensional (3) mesoscale model. Another test showed that at five stations in the central Cascades of Washington over a 143 day period during the 1996/97 winter, it produced substantially better results than an advanced mesoscale precipitation model (Reference Hayes, Rasmussen and ConwayHayes and others, 2002, table 6).

The two-term model for summer balance, with both T ELA and f averaged over days 116–277 (the optimum dates from the one-term model)

(5)

has ∊s = 0.30 m and r 2 = 0.71. Introducing f to account explicitly for accumulation relieves T ELA of needing to do it, and results improve whe is formed from the sum of only T ELA > 0, which is represented by . The coefficients, again determined empirically to minimize the error, are α 2 = −0.39, β 2 = 1.27, γ 2 = −0.48. Over 1959–99, and with ϕ′ = 271°.Sensitivity of results to model parameters t 0, t 1, ϕ′ and T′ (Fig. 3) is determined by varying each one experimentally while holding the others at their optimum values.The parameters are strongly orthogonal, withthe apparent optimum value of one being markedly insensitive to the values adopted for the others.

Improvement of the T, f model over the T model (Equation (1)) is significant at the 0.01 level according to the statistical F test. The 41 observed b s,j and calculated are compared in Figure 4.

Fig. 4. Results from summer balance T, f model (Equation (5)) over 1959–99, using t0 = 116 ( 26 April), t1 = 277 (4 October) and ϕ′ = 271°. Residuals for 1959–78 (solid circles) have mean value −0.07 m and for 1979–99 (open circles) have mean value +0.07 m.

5.3. Correlation of summer balance with meteorological variables

Although both turbulent exchange and the downward longwave flux are increasing functions of T ELA, the models do not explicitly account for variations in solar radiation. High T ELA, however, is usually accompanied by clear weather with greater atmospheric transmittance (low RH), which is indicated by the strong negative correlation between T ELA and RH (Table 1) and the positive correlation with the number of dayswith low RH (below 0.65) and withthe incidence of dry days at Diablo Dam, which is 40 km north of the glacier (Fig. 1). The 0.65 RH threshold was also used to distinguish between wet and dry days in modeling precipitation in Colorado (Reference RheaRhea, 1978) and in Washington (Reference Rasmussen, Conway and HayesRasmussen and others, 2001). In addition, the season of high T ELA is July and August (Fig. 2), when clear weather is most likely. Incoming solar radiation increases with high transmittance associated with clear weather, and by July and August much of the glacier surface has changed from high-albedo snow to low-albedo ice, thus permitting more of the received radiation to be absorbed.

Previous calculations (Reference Rasmussen and ConwayRasmussen and Conway, 2001) show that about 9% of the total accumulation occurs during summer. The summer snow flux increases with the number of cold days in summer and also during cold summers, which is reflected by the positive correlation between and the number of T ELA < 0 days and the negative correlation between and (Table 1),which was calculated by averaging all T ELA < 0 over the summer season. Summer balance becomes less negative a increases, and that occurs with more negative temperatures.

Correlations of b s with other upper-air variables are also shown inTable 1. It correlates positively with 850 hPa relative humidity RH and negatively with the number of days when RH is low. Correlation of b s with the number of low-RH days in spring (r = −0.60) is stronger than with the number in late summer (r = −0.43); albedo is lower in late summer, but sun angles are higher in spring. Correlations of b s with and are strongly negative.

5.4. Distribution of errors

There are several sources of error. The models only crudely represent the relevant physical processes and coarsely sample atmospheric conditions. They use data at only one geographic point, from only two vertical levels and only once a day. The model parameters, however, appear to be little affected by sampling error as shown by the split-sample tests detailed in the Appendix.

Determining summer balance from glacier observations is also subject to several sources of error. It is not measured directly, but is the difference between measurements of net balance and winter balance, which are both susceptible to sampling error in the horizontal. Reference BraithwaiteBraithwaite (1986) found a variance of 0.25 m in net balance measured over 1979–84 at three closely spaced stakes in the ablation zone of a glacier in southwest Greenland. Measuring winter balance is subject to the possibility of not probing exactly to the previous summer surface and of using imprecise density values to convert the thickness to a water equivalent amount. Undetected ablation after the last measurements of the season introduces error in determining net balance b n. The rms of errors that Reference KrimmelKrimmel (2000, p.77) estimated over 1986–91 are δ w = 0.24 m for winter balance and δ n = 0.23 m for net balance. On the assumption that the errors are uncorrelated, because the two balance components are independently measured, the error in summer balance b s is δ s = 0.33 m. That is, (Reference BevingtonBevington, 1969, p.64).

5.5. Effect of changing topography

There is an incoherence between the glacier variable b s and the upper-air variable and in Equations (1) and (5). Conditions aloft at the re-analysis gridpoint are unaffected by changes in the topography of South Cascade Glacier, but published glacier-total balance values describe the effect of climate on a varying glacier topography.Values of b n and b w were determined year by year over a glacier topography that was steadily losing area at low altitude, where both components are more negative than at higher altitude. The effect of loss of area in the region of strongly negative balance was offset slightly by the effect of thinning else-where, which lowered the surface along the vertical gradient of balance ∂b/∂z > 0.

Published values of b n and b w for the early years are less positive than they would have been if the b n(z) and b w(z) from which they were calculated had been integrated over the mean 1959–99 topography, and for the later years are more positive than they would have been if they had been integrated over the mean topography. In a definitive analysis of these effects, Reference Elsberg, Harrison, Echelmeyer and KrimmelElsberg and others (2001) found that the cumulative net balance reported over 1970–97 (−15 m) was 16% more positive than if the b n curves had all been integrated over the 1970 topography. If the constant-topography values differed from the published b n by an increment that grew steadily by a constant amount Δb each year as the glacier shrank, the values for 1970, 1971,…, 1997 would be

from which subtracting b 70 + b 71 + … + b 97 = −15 gives Δb = +0.006 m each year. The effect of changing topography is reflected by the distribution of residuals (Fig. 4), with the mean value of over 1959–78 being −0.07 m and over 1979–99 being +0.07 m. That is, in the early years, when the glacier had more area at low altitude, the published b s values were more negative on average than the model indicates; in the later years, when it had less area there, they were more positive. Most of the scatter between measured and estimated values, however, is due to imperfections in the model, error in meteorological data and error in the mass-balance measurements.

6. Estimating Net Balance

Net balance b n can be estimated by adding from the T, f model (Equation (5)) to estimates of from the winter balance model (Reference Rasmussen and ConwayRasmussen and Conway, 2001). When the models use data from the same re-analysis gridpoint over 1959–99, ∊s = 0.30 m and ∊w = 0.29 m, which yields ∊n = 0.40 m. This is less than the error if the residuals were uncorrelated (0.42) because of a slight negative correlation (r = −0.07) between ∊s and ∊w. Combining the ∊n = 0.40 m and the standard deviation of the net balance σ n = 0.91 in Equation (2) yields r 2 = 0.81. The two balance components b s and b w are positively correlated (r = +0.21), which is thought to result from high albedo persisting longer into the ablation season following heavy winter accumulation, rather than from correlation between winter and summer meteorological conditions.

7. Trends

Discontinuities in South Cascade Glacier summer balance b s, temperature at 2000 m at the gridpoint, and snow flux calculated at 1650 m (Equation (3)) all occur between 1986 and 1987 (Fig. 5). Reference Hodge, Trabant, Krimmel, Heinrichs, March and JosbergerHodge and others (1998) showed the 1986/87 discontinuity in b s to be significant at the 0.01 level according to four statistical tests. Summer snow flux decreased 33% from 0.42 to 0.28 m s−1, due in almost equal parts to drying and warming; that is, a decrease in the total amount of moisture reaching the glacier and in the part of it occurring below +2°C at 1650 m.This trend is also evident in the mean May–September precipitation at Diablo Dam, which decreased 17% from 340 mm over 1959–86 to 280 mm over 1987–99.

Fig. 5. Variations of summer balance, temperature and snow flux over 1959–99. (a) South Cascade Glacier summer balance bs; (b) average temperatur ; (c) average snow flu at 1650 m. Averaging period for T and f is 26 April–4 October. Time of discontinuity between the two stages of the best-fitting piecewise-constant function is determined empirically, the best-fitting constant in a stage being the mean of the values in that stage. For each series the discontinuity is between 1986 and 1987, with the magnitude of the discontinuity shown in the figure.

Trends are expressed here in terms of a step discontinuity superposed on a time series assumed to be stationary both in the stage before the discontinuity and in the stage after it. The best-fitting constant in a stage is the mean of the values in that stage, and the time of discontinuity between one stage and the next is determined empirically to give the best fit overall. Time series approximated by two-stage piecewise-constant functions over 1958–99 can be interpreted as a section of five-stage functions fit to several climatological variables in the North Pacific from 1870 to 1990 (Reference MinobeMinobe, 1997, fig. 1c). Linear fits to the series in Figure 5 have slopes indicating annual rates of change −0.011 m for summer balance, +0.032°C for 2000 m temperature and −0.005 m s−1 for snow flux; in all three cases the piecewise-constant fits explain more of the variance in the time series than the linear fits do.

Sensitivity to temperature is given by the T model (Equation (1)) as ∂b s /∂T = α 1 = −0:50 m °C−1. The T, f model gives a similar result; upon differentiating Equation (5), assuming the perturbation is the same at all altitudes,

(6)

The derivative is determined by calculating f when setting the rain–snow discriminator to +1°C instead of +2°C while holding all other model parameters constant. This summer sensitivity compares with 0.3–0.5 m °C−1 shown by the crosses in figure 3 of Reference Oerlemans and FortuinOerlemans and Fortuin (1992) for 12 selected glaciers.

The summer sensitivity consists of −0.39 m °C−1 due to enhanced ablation, and −0.12 m °C−1 due to reduced accumulation (Equation (6)). Combining these with −0.24 m °C−1 calculated for winter, which is mainly accumulation, shows that the sensitivities of annual accumulation (−0.36 m °C−1) and annual ablation (−0.39 m °C−1) are nearly equal. A slightly lower sensitivity (−0.30 m °C−1) of annual accumulation was estimated in section 4 of Reference Rasmussen, Conway and HayesRasmussen and others (2000) at Blue Glacier, which is 200 km to the west-southwest with about one-third greater accumulation.

Uniform warming of 1°C through the year at all altitudes at South Cascade Glacier would make b w less positive by 0.24 m and b s more negative by 0.51 m, for a total decrease in b n of 0.75 m. The 0.75 m °C−1 sensitivity nearly coincides with the 0.78 m °C−1 calculated for South Cascade Glacier by R. J. Braithwaite (personal communication, 2003), who also found this weak variation with altitude at South Cascade Glacier. The weak variation may not exist at glaciers with large vertical ranges.

Vertical profiles b n(z) for each year 1986–98 have two significant features.The profile for each year is given in Reference KrimmelKrimmel (2000, Reference Krimmel2001), or in reports for individual years cited in Reference KrimmelKrimmel (2001). The13 year mean value of b n at the top of the glacier was only ≈+0.5 m, and in four years (1986, 1987, 1992, 1998) b n(z) was negative over the entire altitude range of the glacier, as it was again in 2001 (Reference KrimmelKrimmel, 2002). Thus, in a climate similar to 1986–98 except for being warmer by ≈ 0.7°C, the mean b n(z) would be negative over the entire altitude range because 0.7°C × 0.75 m °C−1 > 0.5 m. This inference assumes future b n(z) would differ from the average profile over 1986–98 by a constant amount Δb at all z, which is consistent with the assumption of weak variation with altitude of the sensitivity to temperature perturbations. It was shown in section 4.2 of Reference Rasmussen and ConwayRasmussen and Conway (2001) to distinguish the profiles for individual years better than adjusting the profiles by an altitude increment; that is, using b n(z + Δz). The glacier would probably then shrink to a narrow strip along the shaded western edge of the upper part of the glacier, which is the only place where 2000/01 accumulation survived after the 2001 ablation season (Reference KrimmelKrimmel, 2002, fig.14). Profiles b(z) are width-averaged at each z, which is the reason b n can be positive in some places even when b n(z) is negative over the entire altitude range.

8. Discussion and Summary

The most significant result of the analysis is discovering how valuable the upper-air data are. The summer balance model has been tested and applied to Wolverine and Gulkana Glaciers, Alaska, U.S.A., over 1966–99 and Place Glacier, western Canada, over 1965–95. The resulting r 2 were, respectively, 0.55, 0.65 and 0.53, compared with 0.71 at South Cascade Glacier found in this study. It is remarkable that such a simple model does so well. Upper-air conditions must have a profound and direct effect on glacier mass balance; for both summer and winter mass balance, model errors are comparable to measurement errors. Success of the models depends on optimizing their parameters from balance measurements and using a temporal resolution that is sufficient to capture the non-linearities in relations between meteorological conditions and accumulation and ablation.

Reference Oerlemans and ReichertOerlemans and Reichert (2000) used three different types of data to model the net balance for Nigardsbreen, in a maritime setting in Norway. Results using data from a distant low-land station (Bergen) gave r = 0.65; monthly mean data from three nearby weather stations gave r = 0.76; re-analysis data from the European Centre for Medium-Range Weather Forecasts gave r = 0.85, which is slightly lower than r = 0.90 found for South Cascade Glacier in this study.

Meteorological data cannot substitute for glacier measurements, because a balance record is needed to calibrate the model. Nevertheless in the absence of a glacier record, the model can be used to give a qualitative indication of relative balance changes. Taking the critical direction ϕ′ to be in the onshore direction, and interpolating temperature at the altitudes of the equilibrium line and the terminus will probably give results similar to those produced by the optimum values for those model parameters.

The two models could be used to extend the South Cascade Glacier mass-balance series back to 1948, when regular upper-air measurements began, provided major data inhomogeneities in early years can be mitigated. For example, it is not clear whether the apparent jump of ≈1.5°C in the upper-air temperature record in the late 1950s is real or whether it is an artifact of sensor changes. The International Geophysical Year (1957/58) inspired increases in the quantity and quality of radiosonde measurements, and observation times changed on 1 June 1957 from 0300 and 1500 UTC, to 0000 and 1200 UTC (Reference AngellAngell, 1988; Reference KistlerKistler and others, 2001). The history of making and recording radiosonde temperature observations is detailed by Reference GaffenGaffen (1994), and a brief summary of inhomogeneities in the humidity record is given in Reference Rasmussen, Conway and HayesRasmussen and others (2001, section 4).

The models could also be used to predict how mass balance might change under possible future climate scenarios, such as from global circulation models. In the absence of change to either precipitation or atmospheric transmittance, a 1°C warming would make b s more negative by 0.51 m and b w less positive by 0.24 m. The net balance of South Cascade Glacier would decrease by 0.75 m, which would overcome the 1986–98 average net balance of +0.5 m at the top of the glacier.

Acknowledgements

This work was funded from the U.S. National Science Foundation grant OPP-0240861 and NASA grant 1215388. It was possible because of the accomplishment of U.S. Geological Survey researchers in operating the South Cascade Glacier program, principally M. F. Meier, A. S. Post, W.V. Tangborn and R. M. Krimmel, and because of Krimmel’s superb documentation of measurements made. Comments by R. Braithwaite, by an anonymous reviewer and by scientific editor R. Naruse have all produced many improvements in the paper. We thank E. Recker for extracting the re-analysis data from NCAR.

Appendix

When the T model (Equation (1)) was calibrated on only even-numbered years, the optimum averaging interval [t 0, t 1] was days 117–283; it had ∊s = 0.35 and r 2 = 0.60 when applied to the entire period 1959–99.When it was calibrated on only odd-numbered years over the same period, it was days 116–272; it had ∊s = 0.34 and r 2 = 0.62. When it was calibrated on all years, [t 0 , t 1] was days 116–277; over the early period 1959–86 it had ∊s = 0.28 and r 2 = 0.66 and over 1987–99 it had ∊s = 0.38 and r 2 = 0.56.

In the T, f model (Equation (5)), the averaging interval [t 0, t 1] was fixed at the optimum from the T model, and the critical direction ϕ′ for computing the flux f was fixed at the optimum 271° from the winter balance model. When calibrated on only even-numbered years, it had ∊s = 0.30 and r 2 = 0.70 when applied to the entire period 1959–99.When calibrated on only odd-numbered years over the same period, it had ∊s = 0.31 and r 2 = 0.69. When it was calibrated on all years, over 1959–86 it had ∊s = 0.27 and r 2 = 0.69 and over 1987–99 it had ∊s = 0.35 and r 2 = 0.64.

References

Angell, J. K. 1988.Variations and trends in tropospheric and stratospheric global temperatures, 1958–87. J. Climate, 1(12), 12961313.Google Scholar
Bevington, P. R. 1969. Data reduction and error analysis for the physical sciences. NewYork, McGraw-Hill.Google Scholar
Braithwaite, R. J. 1986. Assessment of mass-balance variations within a sparse stake network, Qamanârssûp sermia, West Greenland. J. Glaciol., 32(110), 5053.Google Scholar
Braithwaite, R. J. and Zhang, Y.. 2000. Sensitivity of mass balance of five Swiss glaciers to temperature changes assessed by tuning a degree-day model. J. Glaciol., 46(152), 714.Google Scholar
Conway, H., Rasmussen, L. A. and Hayes, P.. 1995. On the use of radiosondes to model glacier ablation. Ann. Glaciol., 21, 245250.CrossRefGoogle Scholar
Elsberg, D. H., Harrison, W. D., Echelmeyer, K. A. and Krimmel, R. M.. 2001. Quantifying the effects of climate and surface change on glacier mass balance. J. Glaciol., 47(159), 649658.CrossRefGoogle Scholar
Gaffen, D. J. 1994. Temporal inhomogeneities in radiosonde temperature records. J. Geophys. Res., 99(D2), 36673676.Google Scholar
Hayes, P. S., Rasmussen, L. A. and Conway, H.. 2002. Estimating precipitation in the central Cascades of Washington. J. Hydrometeorol., 3(3), 335346.Google Scholar
Hodge, S. M., Trabant, D. C., Krimmel, R. M., Heinrichs, T. A., March, R. S. and Josberger, E. G.. 1998. Climate variations and changes in mass of three glaciers in western North America. J. Climate, 11(9), 21612179.Google Scholar
Kalnay, E. and 21 others. 1996. The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc., 77(3), 437471.Google Scholar
Kistler, R. and 12 others. 2001. The NCEP/NCAR 50-year reanalysis: monthly means CD-ROM and documentation. Bull. Am. Meteorol. Soc., 82(2), 427467.Google Scholar
Krimmel, R. M. 2000.Water, ice, meteorological and speed measurements at South Cascade Glacier, Washington, 1986–1991 balance years. U.S. Geol. Surv.Water-Resour. Invest. Rep. 00-4006.Google Scholar
Krimmel, R. M. 2001. Water, ice, meteorological, and speed measurements at South Cascade Glacier, Washington, 1999 balance year. U.S. Geol. Surv. Water-Resour. Invest. Rep. 00-4625.Google Scholar
Krimmel, R. M. 2002. Water, ice and meteorological measurements at South Cascade Glacier, Washington, 2000–01 balance years. U.S. Geol. Surv.Water-Resour. Invest. Rep. 02-4165.Google Scholar
Miller, C. D. 1969. Chronology of Neoglacial moraines in the Dome Peak area, North Cascade Range, Washington. Arct. Alp. Res., 1(1), 4966.Google Scholar
Minobe, S. 1997. A 50–70 year climatic oscillation over the North Pacific and North America. Geophys. Res. Lett., 24(6), 683686.Google Scholar
Oerlemans, J. 2000. Analysis of a 3 year meteorological record from the ablation zone of Morteratschgletscher, Switzerland: energy and mass balance. J. Glaciol., 46(155), 571579.Google Scholar
Oerlemans, J. and Fortuin, J. P. F.. 1992. Sensitivity of glaciers and small ice caps to greenhouse warming. Science, 258(5079), 115117.Google Scholar
Oerlemans, J. and Reichert, B. K.. 2000. Relating glacier mass balance to meteorological data by using a seasonal sensitivity characteristic. J. Glaciol., 46(152), 16.Google Scholar
Rasmussen, L. A. and Conway, H.. 2001. Estimating South Cascade Glacier (Washington, U.S.A.) mass balance from a distant radiosonde and comparison with Blue Glacier. J. Glaciol., 47(159), 579588.Google Scholar
Rasmussen, L. A., Conway, H. and Hayes, P. S.. 2000. The accumulation regime of Blue Glacier, U.S.A., 1914–96. J. Glaciol., 46(153), 326334.Google Scholar
Rasmussen, L. A., Conway, H. and Hayes, P. S.. 2001. Estimating Olympic Peninsula precipitation from upper-air wind and humidity. J. Geophys. Res., 106(D2), 14931501.Google Scholar
Rhea, J. O. 1978. Orographic precipitation model for hydrometeorological use. (Ph.D. thesis, Colorado State University.)Google Scholar
Torinesi, O., Letréguilly, A. and Valla, F.. 2002. A century reconstruction of the mass balance of Glacier de Sarennes, French Alps. J. Glaciol., 48(160), 142148.Google Scholar
United States Navy. 1966. Selected level temperatures and dew points for the Northern Hemisphere. Washington, DC, U.S. Navy. Chief of Naval Operations. (Document NAVAIR 50-1C-52.)Google Scholar
Figure 0

Fig. 1. Northwestern Washington State. The U.S. National Weather Service radiosonde was at Tatoosh Island prior to 1 August 1966 and at Quillayute afterward.The heavy curve is the smoothed coastline.

Figure 1

Fig. 2. Seasonal variation of temperature at South Cascade Glacier and at re-analysis gridpoint 47.5° N, 125° W. The value at the gridpoint is the 1959–99 mean of the temperature interpolated at 2000 m between the 850 and 700 hPa levels in the NCEP–NCAR re-analysis database. Temperature at the glacier is calculated by applying to it the horizontal gradient of mean monthly 850 hPa temperature (United States Navy, 1966).

Figure 2

Table 1. Correlations over 1959–99

Figure 3

Fig. 3. Sensitivity of summer balance models to values of model parameters: (a) beginning date t0 for temperature averaging; (b) endingdate t1 for temperatureaveraging; (c) critical direction ϕ′ for calculating f (Equations (3) and (4)); (d) threshold temperature T′ above which to sum TELA in Equation (1). Aseach parameter(a–c)intheT; f model(Equation (5)) is varied, the others are held at their optimum values: t0 = 116 (26 April), t1 = 277 (4 October) and ϕ′ = 271°.The T model (Equation (1)) differed only in having optimum t1 = 279; itsoptimum T′ beingbelow the lowest TELA in the period of record indicates all TELA should be used in formin .

Figure 4

Fig. 4. Results from summer balance T, f model (Equation (5)) over 1959–99, using t0 = 116 ( 26 April), t1 = 277 (4 October) and ϕ′ = 271°. Residuals for 1959–78 (solid circles) have mean value −0.07 m and for 1979–99 (open circles) have mean value +0.07 m.

Figure 5

Fig. 5. Variations of summer balance, temperature and snow flux over 1959–99. (a) South Cascade Glacier summer balance bs; (b) average temperatur ; (c) average snow flu at 1650 m. Averaging period for T and f is 26 April–4 October. Time of discontinuity between the two stages of the best-fitting piecewise-constant function is determined empirically, the best-fitting constant in a stage being the mean of the values in that stage. For each series the discontinuity is between 1986 and 1987, with the magnitude of the discontinuity shown in the figure.