Modeling Dry-Snow Densification without Abrupt Transition

An empirical model for the densification of dry snow has been calibrated using strain-rate data from Pine Island Glacier basin, Antarctica. The model provides for a smooth transition between Stage 1 and Stage 2 densification, and leads to an analytical expression for density as a function of depth. It introduces two new parameters with a simple physical basis: transition density ρT and a scaling factor, M, which controls the extent of the transition zone. The standard (Herron and Langway) parameterization is used for strain rates away from the transition zone. Calibration, though tentative, produces best parameter values of ρT = 580 kg m−3 and M = 7 for the region. Using these values, the transition model produces better simulations of snow profiles from Pine Island Glacier basin than the well-established Herron and Langway and Ligtenberg models, both of which postulate abrupt transition. Simulation of density profiles from other sites using M = 7 produces the best values of ρT = 550 kg m−3 for a high accumulation site and 530 kg m−3 for a low accumulation site, suggesting that transition density may vary with climatic conditions. The variation of bubble close-off depth and depth-integrated porosity with mean annual accumulation predicted by the transition model is similar to that predicted by the Simonsen model tuned for Greenland.


Introduction
In the accumulation areas of ice sheets, ice caps, and glaciers, snow is deposited on the surface and, with time, becomes denser until it turns into ice.This process of densification proceeds at a rate that depends on climatic conditions; slowly in the cold, desert regions in the interior of the great polar ice sheets, and more rapidly in warmer regions with higher precipitation.The question of how to calculate this rate from given climatic information is an important aspect of many areas of glaciological research.A recent review [1] cites two problems for which good snow densification models are vital: (1) the determination of past climate from ice-core data and (2) the determination of ice-mass balance changes from altimetry data.In the first case, we wish to know the depth at which snow becomes so dense that bubbles of air are trapped in the ice matrix.This allows an appropriate age to be allocated to chemical constituents of the air.In the second case, we wish to know how changes in climate affect the density profile of snow in order to distinguish between elevation change produced by changes in density and that produced by mass-balance change.Although physics-based snow models, such as CROCUS [2] and SNOWPACK [3,4], were originally developed for hydrological and avalanche forecasting, they have been used successfully to model dry polar snow, but they require detailed initialization and input data that are not always available.There is still a need for simple empirical models that require only the mean annual temperature and accumulation rate as input.
An early suggestion made by Robin [5] has led to a family of empirical-densification models widely used by the glaciological community.He proposed that, in a natural dry-snow cover, the increase in overburden stress is linearly related to the proportional decrease in air volume in the mean density profile.In models based on the Robin hypothesis, strain rate is proportional to a simple function of density.However, from the analysis of density profiles from ice cores, it soon appeared that the proportionality constant changed with density.This problem was overcome by dividing the profile into layers, each with its own constant.Herron and Langway [6] proposed different densification rates for snow with a density of ρ < 550 kg m −3 and 550 kg m −3 ≤ ρ ≤ 800 kg m −3 .Other authors followed suit.This implies a sharp transition between "Stage 1 densification" for lower densities, in which grain-boundary sliding and grain growth are thought to be the dominant mechanisms, and "Stage 2 densification" for higher densities, in which sintering is the primary mechanism.The proposed transition density is the maximum packing density of uniform spheres of ice.However, in natural snow, with a range of grain sizes, high-resolution measurements of density profiles [7] and direct measurements of strain rate [8] suggest that transition density may vary between sites, and that the transition may be gradual.
A better approach to empirical modelling of dry-snow densification might be to accept that strain rate is dependent on a more complex function of density than that proposed by Robin.This would allow a gradual change in densification rate with depth, as the relative importance of grain-boundary sliding, grain growth, and sintering changes.In this paper, a new empirical model of snow densification is proposed, calibrated using field measurements of strain rate as a function of density and validated using high-resolution profiles of ice-core density.The model adds a term to existing models that removes the abrupt transition between Stage 1 and Stage 2.
Snow-densification models need to provide a constitutive law to describe the dependence of volumetric strain rate, ε, on the applied stress, σ, and on the strength of the snow, characterised by its bulk density, temperature and sometimes by other variables such as grain size, coordination number [9,10] and impurity concentration [11].The material-following densification rate is related to ε by the definition where t is time.In polar firn, there are density fluctuations on an annual or subannual scale about vertically smoothed density ρ 0 that monotonically increases with depth.In the upper layer of the snow, densification is enhanced by annual and subannual temperature variations, but the amplitude of these variations declines exponentially with depth until, at ≈ 10 m, temperature is close to mean annual temperature T m .The steady-state approximation applies below the surface layer for constant T m and constant accumulation rate ā.We begin with the hypothesis that ε is a function of ρ 0 and write where ρ i is the density of ice and c(ρ 0 ) is a site-specific function of density.This implies that the effects of microstructure and impurity concentration on snow strength are secondary compared to the effect of density.The Robin hypothesis states that c(ρ 0 ) is a constant and the dependence of strain rate on density is described solely by the function (ρ i − ρ 0 )/ρ 0 .Hidden within this formulation is an assumption that overburden stress σ, which increases with water-equivalent depth, is offset by a corresponding increase in snow strength because of some other factor that also increases with depth, such as time since deposition τ.
In the one-dimensional case when horizontal velocity divergence εH is negligible, ε = εzz , where We defined a co-ordinate system in which vertical distance z and velocity w are positive upward.Given constant accumulation rate ā = −ρ 0 w, Equations ( 2) and (3) give In the special case where c is a constant (the Robin hypothesis), its value can be estimated from a single steady-state profile ρ(z) by plotting the best-fit straight line through logarithmic density function ln(ρ/(ρ i − ρ)) as a function of z over a section of the profile, say from point z where the water-equivalent height is q to point z + Z where the water-equivalent height is q + Q. Writing the (negative) gradient of this line as −kρ i gives where k is a positive constant that we call the vertical densification constant since it is derived from vertical snow density profiles.For simplicity, we call c the density-corrected strain rate whether or not it it is constant.

Time-Varying Conditions
Given measurements of density ρ 1 (q) for a material element at time t and ρ 2 (q) for the same element at time t + ∆t, we may calculate the average value of c(ρ) over time as In the particular case where c is a constant, the value of F (also constant) can be estimated by averaging over a section of the profile from q to q + Q where Q is several years of accumulation.Then Let −kρ i now be the gradient of the best straight line through a plot of ln(ρ 1 /(ρ i − ρ 1 )) as a function of z over the section of the profile from q to q + Q and let ā be the mean annual accumulation over this section.The indirect method of estimating strain rate from a single density profile is based on the assumption that: F ≈ − āk (8) by analogy with the steady-state case (Equation ( 5)).However, if F is not constant over the section of the profile from q to q + Q, this indirect method breaks down.This is because at the particular point q, where F(q) = F, it is not necessarily the case that the local gradient d ln(ρ 1 (q)/(ρ i − ρ 1 (q)))/dz = −kρ i .
It is important to remember that, in reality, the accumulation rate is not constant and density profiles can only give an approximate idea of the time-varying, possibly intermittent, densification process that has produced them.Nevertheless, most existing models depend on profile data for calibration since direct measurements of strain rate are so limited.Since the time between storms is a small percentage (say 1%) of the time taken to build up the profile, this should be a reasonable approach.
If F = F 0 in the Stage 1 region, and F 1 in the Stage 2 region, the Robin hypothesis implies and Other authors have suggested modifications to the Herron and Langway model.For example, Arthern and others [12] suggested that both k 0 and k 1 are independent of ā from consideration of the physics, and found a good fit to Antarctic strain-rate data from sites with ā = 0.13-1.04m w.e. a −1 using for Stage 1 and for Stage 2. Here, the Naborro-Herring creep (activation energy 60 kJ mol −1 ) is offset by grain growth (activation energy 42.4 kJ mol −1 ), so that, for T ≈ T m , the effective activation energy is −17.6 kJ mol −1 , somewhat higher than the Herron and Langway Stage 1 value and lower than the Stage 2 value.The Arthern model is attractive because it is firmly-based on physical principles.However, Ligtenberg and others [13,14] found it necessary to add empirical tuning parameters containing ā to the Arthern model in order to match profiles from Greenland and Antarctica over a wider range of accumulation.They suggest multiplying the r.h.s. of Equation ( 13) by the ratios MO = 1.435 − 0.151 ln 10 3 m w.e.−1 a ā (15) for Antarctic sites and MO = 1.042 − 0.09161 ln 10 3 m w.e.−1 a ā ( 16) for Greenland sites.Similarly the r.h.s. of Equation ( 14) is multiplied by MO = 2.366 − 0.293 ln 10 3 m w.e.−1 a ā .(17) for Antarctic sites and MO = 1.734 − 0.2039 ln 10 3 m w.e.−1 a ā .(18) for Greenland.
Figure 1 shows vertical densification rates k * 0 and k * 1 , implied by these equations, as a function of ā and T m .In both stages, the Ligtenberg values vary with ā because of factor MO, and with T m because the activation energy is not the same as the Herron and Langway value.There is a further variation with ā in Stage 2 because of factor ā−1/2 , which does not appear in the Arthern and Ligtenberg equations.In both Stage 1 and Stage 2, however, the Ligtenberg values lie mostly within the error bounds of the Herron and Langway values, especially taking into account the fact that, in nature, high accumulation is associated with higher temperatures and vice versa.This is reassuring given that the Herron and Langway model was tuned using a relatively limited set of density profiles.Both models are empirical; the difference is that the Ligtenberg model constrains the activation energy and proposes a more complex dependence on accumulation than the Herron and Langway model.Simonsen et al. [15] modified the Arthern model by multiplying the r.h.s of Equation ( 13) by a simple scaling factor f 0 , and the r.h.s. of Equation ( 14) by another scaling factor f 1 and the terms ā−1/2 and exp(−3.8kJ mol −1 /R T m ).This produces a Stage 2 expression that follows the Herron and Langway Stage 2 equation in its dependence on ā and T m .As with the Ligtenberg model, extra parameters are introduced to produce better simulations of field data, but the fundamental assumption of an abrupt transition between Stage 1 and Stage 2 densification is not queried.
In this paper, we tried a different approach.In Section 2, we introduce a new equation for strain rate that allows smooth transition from Stage 1 to Stage 2. We calibrate the parameters of this "transition model" using strain-rate data, making the assumption that the Herron and Langway densification rates for Stage 1 and Stage 2 are correct.We derived the best values of the two new parameters using data from Antarctica and used these values to simulate density-depth profiles from the same area.We then simulated density-depth profiles from other areas and commented on possible changes in parameter values.Finally, in Section 3, we demonstrate the performance of the model using two of the FirnMICE intercomparison experiments.

Transition Model
A simple model that removes the abrupt transition between Stage 1 and Stage 2 densification can be constructed using an activation function for c centered around transition density ρ T , which is likely to be close to 550 kg m −3 (activation functions are used in artificial neural networks and are not connected with the activation energy in an Arrhenius equation).Let where X = (ρ − ρ T )/M 1/2 is a scaled density variable describing the distance away from the transition.Away from the transition, AX 2 » 1 and c → (D −A −1/2 ) or (D + A −1/2 ).D and A can be chosen, so that (D − A −1/2 )= − āk 0 and (D + A −1/2 )= − āk 1 , so that c tends toward its Stage 1 and Stage 2 values.Parameter M controls the abruptness of the transition.This particular activation function was chosen because it leads to an analytic solution for the density profile.Other functions that produce a smooth transition from one value to another exist and produce much of the same effect.Note that the assumption is that k 1 < k 0 .For constant ā, using Equation ( 4) leads to where B = −ρ T /M 1/2 and C = (ρ i − ρ T )/M 1/2 .The solution for the integral is given in Appendix A. Equation ( 20) allows us to calculate the bubble close-off (BCO) depth, important for the determination of past climate from ice-core chemistry.The transition model also leads to an analytic solution for the depth-integrated porosity (DIP), important for the determination of ice-mass balance from altimetry data.
where X s is the surface value of X.The solution for the integral in this equation is given in Appendix B.
The BCO depth and DIP are key indicators used in the FirnMICE intercomparison of densification models [1], and are thus useful quantities for demonstrating the effect of using the transition model (Section 3).

Calibration and Validation
In order to calibrate and validate the new model, detailed measurements of (1) the strain rate as a function of density, or (2) density as a function of depth, are required.There are two sources of such data: repeated nondestructive density profiling in situ using a neutron-scattering technique and high-resolution density profiling of ice cores using gamma-ray attenuation.Over the period from 2004 to 2011, detailed density profiles, at ≈3 cm resolution, were measured in Greenland along the EGIG line [16][17][18] using a neutron probe.From these, strain-rate profiles over periods ranging from a few days to five years could be calculated at the same resolution.Similar repeated density profiles along the iSTAR traverse in the Pine Island Glacier basin, Antarctica, allowed the strain rate of over one year to be calculated [8].Gravimetric measurements of density are also available from ice cores taken at some of these sites.Höhold et al. [7] provide high-resolution (gamma-ray) core density profiles for a wide range of climatic conditions.
Figure 2 shows the mean annual temperature and annual accumulation rate for each site that could potentially provide calibration or validation data.Note that six iSTAR sites and two core density sites have accumulation rates outside the range used to calibrate the Herron and Langway model, and greater than those found along the EGIG line.We used strain-rate data from the Antarctic sites to calibrate the transition model and define its performance with respect to other models using density profiles from the same sites.We then used density profiles from different sites to (tentatively) validate the model.Mean annual temperature and accumulation at sites for which strain rate data are available (in Greenland [16] and Antarctica [8]), where high-resolution (gamma-ray) core-density profiles are available [7], and where the cores used to calibrate the Herron and Langway densificaton equation were collected [6].

Strain-Rate Profiles
At most, iSTAR and EGIG line sites neutron probe profiles were obtained at around 13 m depth.At this depth at sites along the EGIG line, density is still below 550 kg m −3 , so the data can only be used to determine Stage 1 densification rates, and hence the combination of parameters (D − A −1/2 ).Even for iSTAR sites, with relatively high surface density, the transition to Stage 2 is not necessarily completed at 13 m depth, so the strain-rate data cannot be used to derive the parameters A and D separately a priori.However, if Stages 1 and 2 densification rates from an existing model are accepted, then at least some of the neutron-probe data can be used to derive suitable values for M and ρ T .For simplicity, the Herron and Langway expressions for k 0 and k 1 are used (Equations ( 9) and ( 10)) to calculate A and D.
As an example of the data available from the iSTAR traverse, Figure 3 shows density and volumetric strain rate ε as a function of q for iSTAR Site 21 where ā = 0.75 m w.e. a −1 and T m = −22.3• C. The variations in density about the fitted curve show the annual layering.The transition density of ≈550 kg m −3 is reached around 6 m below the surface at q ≈ −3.5 m w.e., suggesting that the transition between Stage 1 and Stage 2 densification should be observable at this site.The volumetric strain rate decreases in magnitude with depth, with enhanced values near the surface where summer warming increases the densification rate.The other iSTAR sites showed a similar pattern of behavior.
Figure 4 shows the average density-corrected volumetric strain rate F measured at iSTAR Site 21 as a function of density and the density-corrected volumetric strain rate c predicted by the transition model using best-fit values M = 2.8 and ρ T = 590 kg m −3 .The curve for c tends to the values expected from the Herron and Langway model away from the transition density.The values of F are very noisy but it is possible to see that the Herron and Langway densification rates for Stages 1 and 2 are reasonable, provided that a transition region is allowed.Best-fit values for all sites, taken individually, show wide variation in M, but the values of ρ T cluster around 580 kg m −3 .Fixing this parameter and optimizing for M produces 22 values with a mean of 17 ± 8. Excluding three values more than three standard deviations from this mean reduces it to 7 ± 2. Figure 5 shows the effect of using M = 7 and ρ T = 580 kg m −3 at Site 21.

Density Profiles
Proper calibration of the transition model requires more extensive strain-rate data, ideally from a site with low surface density, so that there is a clear Stage 1 region unaffected by the surface temperature.In order to observe the full transition to Stage 2, data needed to be collected from depths of the order of 20 m.This poses a methodological challenge for neutron-probe measurements.Nevertheless, it was still interesting to investigate the performance of the model as it stands.With the assumption of a constant accumulation rate, this can be done using profiles of ln(ρ/(ρ i − ρ)) with depth, provided that they are sufficiently detailed.Figure 6 shows that the transition model does indeed produce a better fit to the neutron probe data at Site 21 than the Herron and Langway model, and that the gradient of a straight line fitted to all the data below the Herron and Langway transition point is greater than the true Stage 2 value to which it tends.Morris et al. [8] defined transitional densification rates, k t > k 1 , for densification below ρ = 550 kg m −3 , but using the new transitional model makes these unnecessary.As an aside, we note that in the upper 4 m of the profile the densification rate was enhanced by summer warming.Both models would be able to simulate this if snow temperature T rather than the mean annual temperature T m were used to calculate k 0 in Equation ( 9).In order to quantify the improvement in fit, we define a cost function where values of z are calculated at intervals of 5 kg m −3 from ρ = 500 kg m −3 .For the neutron probe density profiles, ρ(z), the sum goes to ρ = 600 kg m −3 and N = 20.The observed values are taken from the polynomial fit ρ 0 (z).Thus, the cost function describes the goodness of fit of a given model to the smoothed density data in the deeper part of the profile, where the effect of surface-temperature variations is not significant.Table 1 shows the cost functions calculated using neutron-probe data from 22 iSTAR sites with mean annual accumulation ranging between 0.23 and 0.80 m w.e. a −1 .At all sites, the addition of the transition model improves the performance of the Herron and Langway model, i.e., reduces the cost function.At all sites except Site 15, which had the highest accumulation, the Ligtenberg model was also an improvement on the Herron and Langway model.However, it is also the case that at all sites except Sites 10 and 11, which had the lowest accumulation, the transition model performed better than the Ligtenberg model.
At some iSTAR sites, ice cores were collected and gravimetric density measurements made to a depth of 50 m.These may be used to extend the neutron-probe data and test the model over a wider range of density.Figure 7 shows a better fit to the neutron probe and core data at Site 4, again obtained by using the transition model rather than the Herron and Langway alone.At this site, ā = 0.58 m w.e. a −1 and T m = −23.6 • C.
Table 2 shows the cost functions calculated for ρ = 500-800 kg m −3 using a polynomial fit to the core densities to give observed values of z.In this case, N = 60.As before, at all sites, the addition of the transition model improved the performance of the Herron and Langway model, i.e., reduced the cost function.At all sites, including Site 15 this time, the Ligtenberg model was also an improvement on the Herron and Langway model.However, it was again the case that at all sites except Site 10, the transition model performed better than the Ligtenberg model.Although the cost functions in Table 2 place more emphasis on higher densities and less on the transition region than those in Table 1, the conclusions we drew were the same.Independent high-resolution density data are available from the gamma-ray profiling of ice cores, and provide a test of the transition model in regions outside the Pine Island Glacier basin for which it was optimized.Figure 8 shows data from an ice core (B39) collected in Dronning Maud Land, Antarctica [19], at a site where ā = 0.77 m w.e. a −1 and T m = −17.9• C. The data are available at 1 mm resolution but, in the figure, the data have been smoothed by a 31-point running mean (≈3 cm) for clarity.The model lines terminate at the BCO horizon (the point when density reaches 815 kg m −3 [1]) below which they are not expected to be valid.The cost functions for the Herron and Langway and Ligtenberg models for the density range 500-800 kg m −3 are Ψ = 0.163 and Ψ = 0.156, respectively.Using the transition model with ρ T = 580 kg m −3 and M = 7 produces a cost function of Ψ = 0.128.That is, both the Ligtenberg and the transition model improve on the Herron and Langway fit, with the transition model being better for this high-accumulation site.If ρ T is reduced to 550 kg m −3 an improved fit with Ψ = 0.059 is obtained with the transition model.Note that none of the models predicts the slight increase in densification rate just below the BCO horizon associated with an increase in the density variability.For this, a more complex model including the effect of microstructure and impurities on snow strength would be required [11].In their analysis of high-resolution (gamma-ray) density profiles, Hörhold et al. [7] noted a weak transition in the slope of the density-depth profiles at densities between 550 and 580 kg m −3 for high accumulation sites (such as B39), whereas lower accumulation sites (such as B26) showed this transition at much lower densities, below 500 kg m −3 .Figure 9 shows that a slightly improved fit to the data was obtained with the transition model using ρ T = 530 kg m −3 for core B26 collected in North Greenland [20] at a site with T m = −30.6 • C and ā = 0.18 m w.e. a −1 .The Herron and Langway cost function Ψ = 0.087 is reduced to Ψ = 0.085.Using ρ T = 580 kg m −3 increased the cost function to Ψ = 0.118.However, even with this transition density, the model performed better that the Ligtenberg model with Ψ = 0.170.It is clear that the transition model could be improved by allowing ρ T to vary with T m or ā.When more strain-rate data are collected, or more high-resolution (gamma-ray) density profiles are released, this can be investigated.

Discussion
The results in the previous section show that the transition model is capable of simulating density profiles through the transition between Stage 1 and Stage 2 more accurately than models with an abrupt transition.This is useful for glaciological problems that require knowledge of transition densities, such as the determination of annual accumulation series from airborne radar data or estimation of past climate from near-surface temperature profiles.More importantly, we note that density profiles are a valuable indirect source of information on strain rates in dry polar snow.By fitting Equation ( 20) to a density profile we can deduce the parameters of the constitutive law for strain rate, Equation (19).Because of the limitations of the available data, in this paper we have chosen to fix two model parameters (D and A) using the Herron and Langway model, and vary the two remaining parameters (ρ T and M).In the future, when more extensive high-resolution measurements of ice-core density are available, one might optimise for all parameters and possibly improve on the Herron and Langway empirical expressions for Stage 1 and Stage 2 densification.
If Stage 2 densification rates are less than Stage 1 rates, the effect of adding a transition zone centered on density ρ T = 580 kg m −3 greater than the usual transition density of 550 kg m −3 will be to increase predicted Stage 2 densities.This would be the case whatever the shape of the activation function and is because densification proceeds at a higher rate over the transition zone before it settles to the Stage 2 value.How important this "over-shoot" is depends on the difference between k 0 and k 1 .If the Herron and Langway values are used, it is clear from Equations ( 9) and ( 10) that this difference increases with increasing ā.The consequence is that the transition model predicts BCO depths that are less than those predicted by the Herron and Langway model, by an amount that increases with ā.For the same reason, the depth-integrated porosity (DIP) predicted by the transition model is less than that predicted by the Herron and Langway model.There is some evidence [7] that, at very low accumulation rates, transition density may be less than 550 kg m −3 , in which case the BCO depth and DIP would be greater than the Herron and Langway values.This point needs to be explored when more high-resolution field data become available.
Figure 10 shows BCO as a function of accumulation (in this case given in m ice a −1 ) for T m = −30 • C, ρ S = 360 kg m −3 and can be compared to Figure 2e in the FirnMICE comparison paper [1].Using the transition model produces BCO depths that are similar to those predicted by the Simonsen model (using f 0 = 0.8 and f 1 = 1.25), but greater than those predicted by the Ligtenberg model for Antarctic sites.From Figure 1 we see that, for ā = 0.1-0.3m w.e. a −1 , the Ligtenberg Stage 2 vertical densification rates are greater than the Herron and Langway rates, so that the difference between the curves in Figures 10 and 11 is as expected.Interestingly, the two tuning parameters added to the Herron and Langway model by Simonsen have much the same effect as the two parameters used in the transition model.However, before taking this comparison too far, we need to remember that the Ligtenberg and Simonsen models were designed for use with different Regional Climate Models (RACMO and HIRHAM5), and the values of ā and T m derived from these RCMs are not necessarily the same as the values that would be measured in situ.

Conclusions
The transition model provides an alternative way of adapting the Herron and Langway model for density-corrected volumetric strain rates c 0 and c 1 to produce a better fit to field data, using the same number of extra parameters as the Simonsen model, and fewer parameters than the Ligtenberg model.Of course, it would be possible to combine the transition model with these other models for c 0 and c 1 , although this would produce a less parsimonious model.The attraction of the transition model is that it is based on a simple physical concept that is already well-established in the glaciological literature, that the transition between grain-boundary sliding and sintering as mechanisms for snow densification is a gradual one, and, hence, the two extra parameters, ρ T and M, which describe the location and width of the transition, have a simple physical basis.

Figure 1 .
Figure 1.Vertical densification rates for (a) Stage 1 and (b) Stage 2 densification.Black lines show the Herron and Langway values (solid lines) with their uncertainties (dashed lines).Colored lines show the Ligtenberg values for mean annual temperatures of 230 (blue), 240 (cyan), 250 (green), 260 (magenta), and 270 K (red), with solid lines for Antarctic sites and dashed lines for Greenland sites.

Figure 2 .
Figure 2.Mean annual temperature and accumulation at sites for which strain rate data are available (in Greenland[16] and Antarctica[8]), where high-resolution (gamma-ray) core-density profiles are available[7], and where the cores used to calibrate the Herron and Langway densificaton equation were collected[6].

1 )Figure 5 .
Figure 5. Measured values of F as a function of density for iSTAR Site 21 (red dots), values of c for Stage 1 and Stage 2 from the Herron and Langway model (cyan lines), and c as a function of density from the transition model with ρ T = 580 kg m −3 , M = 7 (black line).

Figure 10 .
Figure 10.Variation of steady-state BCO depth with ā using the Herron and Langway (brown line), Ligtenberg (green line), Simonsen (pink line), and transition (blue line) models.Simonsen and Ligtenberg results are taken from [1].

Figure 11 .
Figure 11.Variation of steady-state DIP with ā using the Herron and Langway (brown line), Ligtenberg (green line), Simonsen (pink line) and transition (blue line) models.Simonsen and Ligtenberg results are taken from [1].

Table 1 .
Cost functions over the range of 500-600 kg m −3 for the Herron and Langway, Ligtenberg, and transition models for sites in the Pine Island Glacier basin.

Table 2 .
Cost functions over the range 500-800 kg m −3 for the Herron and Langway, Ligtenberg, and transition models for sites in the Pine Island Glacier basin.