2.1. Earlier work

Previous snow indentation studies have been carried out at larger scales and under different conditions than those utilized here. Reference SeligmanSeligman (1936) documented indentation tests on seasonal snow, with an eye toward avalanche studies and modeling. Reference Yong and FukueYong and Fukue (1977) performed confined compression tests in natural snows, both fresh and aged, and in artificial snows that were ground from ice. Their samples were cylinders, 5 cm in diameter and 4.6 cm high. Reference Edens and BrownEdens and Brown (1991) investigated the correlation between microstructure and macrostructural strength by examining morphological changes under snow compression, such as coordination number, grain size, bond radius, neck length, pore size, free surface area and grains per volume. Reference Shoop, Alger and NewcombShoop and Alger (1998) carried out tests, indenting a 20 cm diameter plate into artificial snow, and compared results with a computer modeling study. Additional references on snow indentation are provided by Reference LeeLee (2009).

Other geometries have also been used for indentation testing. These have been based on the rammsonde cone penetrometer that is also used in geology and soils science (e.g. Reference GublerGubler, 1975; Reference Dowd and BrownDowd and Brown, 1986). The snow micropenetrometer (SMP) we used for snow characterization has also been used to measure snow strength (Reference Johnson and SchneebeliJohnson and Schneebeli, 1999; Reference Schneebeli, Pielmeier and JohnsonSchneebeli and others, 1999; Reference Marshall and JohnsonMarshall and Johnson, 2009). Others have modified the cone to have a round tip (Reference Floyer and JamiesonFloyer and Jamieson, 2010).

It has been suggested that there may be some similarities between the mechanical properties of snow and foams, since both are porous, random heterogeneous materials (Reference Kirchner, Michot, Narita and SuzukiKirchner and others, 2001; Reference Yuan, Lee and GuilkeyYuan and others, 2010). Thus, earlier snow studies and studies on metallic and polymeric foams (Reference Olurin, Fleck and AshbyOlurin and others, 2000; Reference Andrews, Gioux, Onck and GibsonAndrews and others, 2001; Reference Kirchner, Michot, Narita and SuzukiKirchner and others, 2001; Reference Ramamurty and KumaranRamamurty and Kumaran, 2004; Reference Lu, Shen, Hou, Ruan and OngLu and others, 2008; Reference Fores-Johnson and LiFlores-Johnson and Li, 2010) are examined here.

Snow has also been described as a bonded granular material (Reference McClung and SelvaduraiMcClung, 1981; Reference BrownBrown, 1994), especially at low, plastic deformation rates. Such deformation rates induce time-dependent ductile failure. Stereology theory is used to calculate the snow strength based on the microstructure, and the timescales are sufficiently long that the snow’s metamorphism must be accounted for during deformation. In our study, deformation rates took place in the brittle region. Reference Kinosita and OuraKinosita (1967) provides a discussion of the strain rates that yield ductile vs brittle failure in snow.

Several different mechanisms for snow deformation are compared and modeled by Reference MellorMellor (1975) and Reference Johnson and HopkinsJohnson and Hopkins (2005).

We apply mathematical models to the laboratory test data. Reference Olurin, Fleck and AshbyOlurin and others (2000) studied metallic foams under cylindrical indentation, to examine the potential for energy absorption under impact. Flat-bottomed, cylindrical punches were used, and the results from our snow testing, which used similar cylindrical pins, are compared against these results. Reference Andrews, Gioux, Onck and GibsonAndrews and others (2001) studied energy absorption by metallic foams, and performed indentation tests with various geometries. Size effects were examined, and relationships between first peak strengths and pin sizes were determined. The results from our snow tests are compared against these results, as well as results from snow models (Reference LeeLee, 2009).

Others who have examined the behavior of foams utilizing small-scale testing include Reference Onck, Andrews and GibsonOnck and others (2001), Reference Ramamurty and KumaranRamamurty and Kumaran (2004) and Reference Lu, Shen, Hou, Ruan and OngLu and others (2008). Reference Fores-Johnson and LiFlores-Johnson and Li (2010) examined the behavior of polymeric foams, and Reference Kirchner, Michot, Narita and SuzukiKirchner and others (2001) modeled snow as a foam of ice.

Our study builds upon the earlier work, with the following key features intended for application to ground vehicle traversal:

• Size effects on mesoscales are examined.

• High strain rates are used, so testing occurs fully in the brittle region.

• Natural snows are used.

• Snow selection is based on what is found most commonly in populated regions.

• Flat-bottomed cylindrical indenters are used.

2.2. Indentation model

The indentation model developed by Reference LeeLee (2009) uses cohesion, internal friction angle and hydrostatic-stress/volumetric-plastic-strain relationships in a plasticity model to predict stress/displacement curves. The curves are classified into three zones, as shown in a theoretical curve in Figure 1.

Fig. 1. Theoretical stress/displacement relationship for indentation into snow.

First, the stress rises to a peak, before an immediate loss of strength. This first peak is called the first failure strength of the snow. After these first transient effects, the indentation proceeds through three zones discussed in Reference LeeLee’s (2009) model. Zone I exhibits a linear (but plastic) rise in stress as the indenter proceeds through the snow. During zone I, the indenter accumulates a densified zone beneath it called a pressure bulb. Reference Muro and O’BrienMuro and O’Brien (2004) cite several papers (e.g. Reference YosidaYosida and others, 1956; Reference Muro and YongMuro and Yong, 1980) that discuss this densification zone. Both an initial first peak and a densification zone are also mentioned for a round-tipped penetrometer (Reference Floyer and JamiesonFloyer and Jamieson, 2010) and for a vertically loaded plate (Reference Shoop, Alger and NewcombShoop and Alger, 1998).

At the transition to zone II, the pressure bulb has reached its maximum size and it proceeds to move through the snow together with the pin, usually with a slight increase in resistive force due to snow hardening. The average stress during zone II is henceforth called the ‘plateau strength’.

Zone III is the region where the pressure bulb reaches the bottom of the snow and further compacts. This only occurs during deep indentation, and is provided by the bottom of the snow container in the laboratory, or the frozen soil or street surface in a ‘real world’ application. This is called a densification (finite depth) zone.

2.3. Plateau strength analysis of foams

Reference Olurin, Fleck and AshbyOlurin and others (2000) performed similar indentation tests on open-cell aluminum foams and calculated the resistive strength during indentation, separated into two components, a crushing term and a tearing term. The crushing term relates to material being compacted beneath the pin, and is proportional to the leading-face surface area. The tearing term relates to damage caused to the snow by the circumferential edge of the pin as it creates new surface while driving into the snow. It is therefore proportional to the pin circumference. Olurin reasoned that

where F _PL is the plateau strength (units of force length rather than stress), r is the pin radius and σ _crush and γ _tear are the crushing and tearing terms respectively. Thus, a plot of F _PL/(πr) as a function of pin radius is expected to be linear, with the slope being the crush strength, and the intercept giving the tearing strength.

3.1. Collection and storage of snow

Virgin (undisturbed), metamorphically stable, dry snow was collected in and around Fairbanks in locations far from sources of contamination, at temperatures of −9 to −30°C. Typically, snow was gathered a week after the previous snowfall. This gave adequate time for sufficient metamorphism to provide roughly uniform grains of approximately spherical shape. Air temperature, snow surface temperature and snow pit-temperature were recorded at each snow collection site.

Two types of snow were collected. The first was fine-grained (average grain size <1 mm) snow such as occurs near the surface of a snowpack. This is snow that has fallen within the previous month or so. The second was coarse-grained (average grain size >1.5 mm) depth hoar such as occurs at depths greater than 20 cm under a surface cover and typically where there is ice over liquid water, such as on a river or a pond. This is snow that has fallen at least 3 months earlier, and has been kept at a thermal gradient by exposure to liquid water beneath the ice, and insulated above by fresher snow layers.

Samples were collected in chilled plastic containers and transported to laboratory freezers in a cooler. Laboratory freezers were maintained at −30°C, cold enough to slow metamorphism and sintering to a very low rate (after 6 months in storage, little sintering had taken place among the snow grains). Storage time ranged from <1 day to 6 months.

3.2. Snow characterization

Snow was selected and characterized both in situ and in the laboratory using multiple methods, to ensure that snow sample-to-sample variation was minimal. The first in situ test utilized a high-resolution SMP. Samples were then removed from the snow bank and measured, weighed and sieved on-site to obtain density and grain-size distribution. Specimens collected and brought to the laboratory were additionally characterized via optical inspection with a light microscope and CT (computed tomography) scanning using three-dimensional X-ray microtomography.

The SMP used was developed by Reference Johnson and SchneebeliJohnson and Schneebeli (1999) for use by avalanche researchers. In this study, it was used to check snow properties in situ so that appropriate snow samples could be selected. Output of the SMP is penetration force and texture index, which is an empirical number used to describe the snow grain geometry. It is defined as the ratio of mean grain size (mm) to density (kg m⁻³).

The density was measured at the field collection sites using a graduated cylinder of known mass, hanging from a spring scale. This gave a quick, rough estimate of density that was later compared with the density obtained from the CT scan images.

Grain-size distribution was also obtained by sieving, with sieves of 1.4, 1 and 0.425 mm. From this sieving, average grain sizes were estimated to be 1 mm for the fine-grained snow and 2.2 mm for the coarse-grained snow.

CT scanning took place in a cold room, and samples were scanned before and after testing at 23.7 μm resolution using parameters as recommended by Reference YuanYuan (2007).

The SMP readings, density measurements, weighing and sieving and CT scanning all established baseline microstructures for different natural snows in their natural state, using the International Classification for Seasonal Snow on the Ground (ICSSG; Reference FierzFierz and others, 2009).

3.3. Establishment of testing parameters

3.3.2. Pin selection

Cylindrical indenters were chosen because they are most typically used in vehicle-terrain studies. They also have several advantages that make them ideally suited to this study. First, to study size effects, it is necessary to have an indenter for which the size ratios among the different linear dimensions do change with changing indenter size. To illustrate, one of the pieces of data captured in this study was the initial failure strength of the snow, which is a function of initial geometry. The differences among initial geometry due to size are illustrated in Figure 2.

Fig. 2. Size effects of different indenter geometries.

Below the green line, a larger radius and a smaller radius cone present the same leading-face geometry to the test material (the snow). Therefore, a cone tip test cannot yield a size effect for initial failure strength. A change in cone angle would change the tip geometry, but any differences shown in such a study would be due to a change in angle as well as linear measurements.

For a blade, the linear dimensions change in one direction, but not in the other (the blade thickness). Therefore, a blade test yields a less pronounced size effect for initial failure strength. Other geometric cross sections will work as well as a circular cylinder, but a circular cylinder was chosen to simplify the calculations, and to avoid adding vertical edges.

3.3.3. Confirmation of indentation model using CT scans

The CT scan images confirmed that densification did occur during snow indentation. This confirms the theory described in Section 2.2 and further discussed in Section 4.1. Figure 3 shows top, front and side cross-sectional views of a column of snow after indentation. The post-processing of CT scans enables analysis of snow particles slice by slice. Figure 4 shows a lower slice of a snow sample before and after indentation with a 12.7 mm diameter pin.

Fig. 3. CT scans of indented snow (white is ice).

Fig. 4. Cross sections of CT scan of indented snow, 12.7 mm indenter (white is ice).

Using these images, the original density of this snow sample was found to be 186 kg m⁻³ (80% porosity), and the density of the region directly beneath the indenter was found to be 649 kg m⁻³ (29% porosity). The change in density indicates whether the pin is compacting the snow or drilling through the snow. The larger the pin, the more compaction took place. Figure 5 shows the comparison fora 6.35 mm pin.

Fig. 5. Cross sections of CT scan of indented snow, 6.35 mm indenter (white is ice).

Using these images, the original density of this sample was 262 kg m⁻³ (71% porosity), and the compressed density under the pin after indentation was 270 kg m⁻³ (70% porosity). As discussed in Section 4.1, less compaction took place with smaller indenter-size to average grain-size ratio. More compaction can be seen in Figure 4 than in Figure 5 (the higher percentage of white regions indicates a higher density of ice). Both Figures 4 and 5 are consistent with the theory, described in Section 2.3, of a densification zone beneath the pin.

The CT scanning also provides a criterion for snow container size during testing as follows. From the post-indentation scans, the horizontal slices beneath the indenter were analyzed. For a given snow slice, areal density was measured, of circles and annular regions of increasing size, centered on the spot directly beneath the middle of the indenter. When the density of an annular region was about the density of virgin snow, it was determined that the pressure bulb had been contained within the inner circle. After this test was performed with the largest pin sizes and with a safety factor of 1.5, it was determined that with the container diameter four times the pin diameter, containment effects were negligible.

3.4. Tribometer set-up

A CETR-UMT tribometer was used for indentation testing. This tribometer has high-resolution load cells and is easily adapted to individual tests.

Two different load cells were used for their different capacities: one had a range of 5–500 mN and a resolution of 0.05 mN; the other had a range of 2–200 N and resolution of 0.01 N. The test region was maintained at a temperature of −20°C using a low-temperature forced airstream, set to an output of −40°C. The reason for this temperature difference is that the enclosure has some inevitable leakage due to the necessity of openings for the various interfaces with the tribometer.

Figure 6 shows an indentation test in progress (the enclosure was opened only for photography purposes; in actual testing, it remained closed). This example shows a 3.175 mm pin being indented into a small amount of snow. As discussed in Section 3.3.3, larger containers were used for the larger pin diameters, to eliminate any confinement effects.

Fig. 6. Indentation test set-up.

The sample holders were made of flat-bottomed plastic cups to minimize thermal contamination, and the pins were machined of aluminum, which has a good balance of light weight to minimize inertial effects, and high stiffness so that only the effects of the snow are captured during the indentation. In other words, the stiffness of the aluminum was considered to be infinite compared with the stiffness of the snow.

3.5. Test parameters

Tests were carried out with the parameters shown in Table 1. Some variations were made due to logistical constraints, and were not part of the variations under study.

Table 1. Indentation test parameters

Varying different parameters and repeating tests with ambiguous results, ∼400 indentation tests were performed in total. Of these, only those that gave clear results and the values needed for a given calculation (e.g. plateau strength for a plateau strength calculation) were used in the data processing.

The ductile-to-brittle transition normal strain rate for snow is ∼10⁻⁴ s⁻¹ (Reference Kirchner, Michot, Narita and SuzukiKirchner and others, 2001). Strain rates for this set of testing were selected to be well above this strain, but low enough to reduce inertial effects. Thus, all tests were well within the brittle range.

In total, 118 sets of tests were performed, with 3–6 runs of each test. Of these 600-odd data points, ∼100–200 were used for each data report.

4.1. Stress/displacement curves

The stress/displacement curves qualitatively reflected the model described in Section 2.2. However, not all three zones were present in all cases. The larger the indenter size relative to the average snow grain size, the more distinct the zones. At the limit of very small indenter-size to snow grain-size ratio, the first peak became unrepeatable as it depended too highly on initial placement of the pin relative to the snow grains. Also, for small indenters, the effect was more that of the pin drilling rather than of the cylinder compressing, and no pressure bulb accumulated, so that zone I was not represented. At the limit of deep snow, where the pressure bulb did not contact the bottom surface, zone III did not occur. At the limit of very shallow snow, zone II was not represented because the pressure bulb hit the bottom before it attained its maximum size.

Figure 7 shows the stress/displacement relationship for four 12.7 mm pin indentations into fine-grained snow. The horizontal axis is the pin displacement into the snow. The vertical axis is the stress, i.e. the resultant force detected by the load cell coupled to the pin, divided by the cross-sectional area of the pin.

Fig. 7. Stress/displacement curves for 12.7 mm pin indentations into fine snow.

Figure 8 shows the stress/displacement relationships for two different pin sizes and two different snow types.

Fig. 8. Stress/displacement curves comparing pin sizes and snow types.

As is discussed in Sections 4.2 and 4.3, the overall forces are higher for the smaller pin. It is also qualitatively observed that zones I and III do not occur for the 3.175 mm pins, since snow was not compacted in those tests. The curves show a transition from the initial peak directly to zone II. Moreover, we see that the 3.175 mm indentation into fine-grained snow yielded a smoother curve than the indentation into the coarse-grained snow, because the pin was large enough with respect to the grain size to read an average force over multiple snow grains. In contrast, the indentation into coarse-grained snow yielded a much less consistent force because the indenter size was on the same order of magnitude as the snow grain size. Thus the behavior of individual snow grains was detected.

Figure 9 shows the same 12.7mm indentation data as in Figure 8, but with the 3.175 mm data removed so that the detail of the larger pin indentation may be seen.

Fig. 9. Stress/displacement curves for 12.7 mm pin indentations, comparing two different snow types.

Here, where the indenter size is sufficiently large that the ratio of indenter size to snow grain size is much greater than 1 in both cases, both graphs look similar.

4.2. First peak strength

The size of the indenter and the average grain size of the snow were varied to determine their effects on first peak strength. It was found (Fig. 10) that the larger the indenter, the lower the first peak failure strength. On log_e scales, a linear relationship is evident so that a power law is apparent:

where σ _fs is the failure stress, γ is a dimensionless material property, A is the pin area and k is a constant.

Fig. 10. Example of first failure strength data, fine-grained snow.

Table 2 shows results obtained for γ against other parameters. Each test ID represents four to six indenter sizes, with three to six indentations repeated for each size.

Table 2. Exponent from a range of tests

The average value of γ is 0.83 (with a coefficient of variation (COV) of 0.16) for fine-grained snow, and 0.85 (with a COV of 0.23) for coarse-grained snow. The average value for both snow types is 0.84, with a COV of 0.15.

The actual values of first peak strength had high test-totest variability (as high as an order of magnitude), even when the same snow sample was used and all other test conditions were ostensibly the same. The source of this variability is not understood, but it is specific to the first peak strength values over time. On the same day, first peak strength values were highly repeatable. Also, plateau strength values (Section 4.3) did not show the same day-to-day variability. However, all the different sets of data yielded the fairly consistent values of γ shown in Table 2.

The other size relationship explored was that of grain size. Identical tests were performed on both fine- and coarse-grained snow. The first peak strengths had the same mean values and same COV for both fine- and coarse-grained snow, although grain size did affect other snow properties as discussed in Section 4.1.

4.3. Plateau strength (zone II)

The mean plateau strengths for several different test sets are given in Table 3, in which ∼150 tests are represented, with each test yielding one data point. Typical results are shown in Figure 11a and b. Each pin diameter in each test set was tested four times, but not all gave a measurable plateau.

Table 3. Mean plateau strengths from a range of tests

Fig. 11. Relationship between plateau strength and pin diameter, for (a) fine grained snow and (b) coarse-grained snow (there is only one data point at 12.7 mm).

About three-quarters of the plateau strengths were determined ‘semi-automatically’, by criteria based on the specific test run. The plateau was defined by the rule: ‘when the moving average of X points changes by more than Y and remains stable for Z data points, then Z is called the plateau’. X, Yand Z were manually selected on a test-by-test basis, due to the large test-to-test variation. A quarter of the resulting plateau strengths came out completely inappropriately; plateau ranges were hand-selected using careful visual judgment.

The raw data for the plateau strengths shows a monotonically decreasing trend with increasing diameter. With such a large spread and with only four different diameters yielding plateau strengths, the mathematical form of the relationship is difficult to determine. However, once the data are analyzed for material properties as discussed in Section 2.3, the plot output from these calculations is much more useful. Representative results of plateau force, F _PL/(πr) (Section 2.3) for two different test sets are shown in Figure 12a and b. All test sets reflect multiple days of testing, often in different years, so the spread of data is quite large. However, this enabled more data points to be used, ensuring broad representation. Each data point represents one test.

Fig. 12. (a) Plateau force/πr as a function of pin radius at 2 mm s⁻¹ for (a) fine-grained and (b) coarse-grained snow.

Results from the four datasets that displayed a linear trend are shown in Table 4, in which ∼150 tests are represented, with each test yielding one data point. The apparent negative values of the shear term can be taken to be zero. This seems to indicate that the tearing term is indeed negligible.

Table 4. Compaction and shear terms from a range of tests

4.4. Absorption of energy

Total energy absorption for each indentation, from first contact with the snow to the completion of the test (in this case, 10 mm indentation), was calculated by taking the area under the force/displacement curve. The entire curve was used because energy absorption, unlike first peak strength or plateau strength, is a property that is cumulative with indentation depth. The deeper the snow, the more stable (with respect to depth) is the measured energy absorption density, since the transient initial effects from the first peak and zone I become much smaller than the effects from the plateau region of zone II.

To normalize against geometric factors and arrive at a dimensionless number, the volumetric absorbed energy density was calculated by dividing the absorbed energy by the volume of snow displaced and by Young’s modulus for ice, for which we used 9 GPa (Reference SchulsonSchulson, 1999). Figure 13a and b show results from all tests, of energy absorption density, with no outliers removed. Some smaller scatter was expected in these data due to the cumulative nature of the calculation of the energy absorption. Note that energy absorption data could be collected for the two largest pin sizes, 34 and 42 mm, yet data for these two pin sizes yielded neither a first peak, nor a plateau strength. Additionally, note that fine-grained snow showed greater scatter than coarse-grained snow. This indicates a greater variability of snow behavior in fine-grained snow than coarse-grained.

Fig. 13. Energy absorption density during pin indentation into (a) fine-grained and (b) coarse-grained snow.

It was found that the average energy absorption density was 3.4 × 10⁻⁷ for fine-grained snow, and 4.1 × 10⁻⁷ for coarse-grained snow. These preliminary results suggest that for cylindrical indentation in this size range, energy absorption density was constant with respect to diameter for fine-grained snow, and decreased with increasing diameter for coarse-grained snow. Note that this decrease was not linear, flattening out to suggest energy absorption independence of size at larger diameters.

5.1. First peak strength

First peak snow strength was highly dependent upon surface conditions and showed high test-to-test variation. Even when known factors were carefully controlled, a different day would yield different first peak strengths, even though the spread of first peak data for a given day was very small. This suggested variations were due not to the equipment or operator, but to intrinsic changes that occur in the surface properties of the snow with storage time. The variability between tests was much smaller for plateau strengths, suggesting that changes that led to the high variability in initial failure strength did not affect snow behavior beyond the first failure. Once the pressure bulb had accumulated, the bulb itself exhibited volumetric properties that were consistent from test to test. By the time the test had reached zone II, the transient effects of the initial conditions had effectively disappeared and the plateau strength better reflected the intrinsic properties of the snow.

The value of the exponent in the power law relationship between first peak strength and pin area was consistent from snow sample to snow sample, so although the absolute first peak strength was highly sample-dependent, the magnitude of the size effect was consistent.

The average value of the power law exponent for both snow types is 0.84. As a comparison, the exponent for brittle foams was estimated to be 0.5 (Reference Gibson and MFGibson and Ashby, 1999), lower than that for the snow tested, indicating that the size effect displayed by first peak strength in snow is greater than that for brittle foams.

5.2. Plateau strength

The plateau strengths used to calculate the crush strength of snow were repeatable and are believed to reflect fundamental properties of the bulk snow. The plateau strength is also of more use to engineers and scientists, as it is reasonable to expect a tire, ski, sled runner, animal paw or hoof to compress the snow a small amount before coming to a stop after a certain compression amount. What limits this compression amount is of more interest than the initial collapse strength.

Reference Olurin, Fleck and AshbyOlurin’s (2000) values for σ _crush for aluminum foam range from 1.4 to 2.8 MPa. Data for the bulk material properties of Alporas, the trade name for the aluminum foam used in Olurin’s studies, could not be found, but the ultimate strength of a few randomly selected brittle aluminum alloys on www.matweb.com is in the 500 MPa range. The ultimate brittle failure strength of water ice, in comparison, is ∼10 MPa (Reference GoldGold, 1977). The ratio of ultimate failure strength of bulk material to crush strength of associated foam is two orders of magnitude for aluminum/aluminum foam and four orders of magnitude for ice/snow. Although a direct comparison cannot be made because the exact material properties of bulk Alporas are unknown, it seems reasonable that the ratio of snow strength to ice strength be lower than the ratio of metal foam strength to bulk metal strength, because snow is formed from ice crystals sitting in contact, with perhaps some, but certainly no guaranteed, sintered bonds between individual crystals. Alporas foam, in contrast, is formed by injecting air bubbles into a molten metal matrix, then cooling it (Reference Miyoshi, Itoh, Akiyama and KitaharaMiyoshi and others, 2000). This is counter to ice, and, in theory, the metal matrix could have the same crystal structure as the bulk metal. It is also of interest to note that for the snow density range studied here, the ratio of snow tensile strength to that of ice can be very small as summarized by Reference PetrovicPetrovic (2003).

Despite the higher test-to-test repeatability and fundamental nature of the plateau strength, it still exhibits considerable scatter, This is to be expected from a natural, non-engineered material. Smaller scatter has been obtained for artificial snows (Reference Shoop, Alger and NewcombShoop and Alger, 1998). Differences in mechanical performance among natural, sorted and artificial snows are discussed by Reference Yong and FukueYong and Fukue (1977). The literature on synthetic foams also exhibits very little scatter in engineered materials (Reference Andrews, Gioux, Onck and GibsonAndrews and others, 2001; Reference Ramamurty and KumaranRamamurty and Kumaran, 2004; Reference Lu, Shen, Hou, Ruan and OngLu and others, 2008; Reference Fores-Johnson and LiFlores-Johnson and Li, 2010).

5.3. Energy absorption

Coarse-grained snow exhibited decreasing energy absorption ability with increasing indenter diameter. This decrease was not linear, flattening out to suggest energy absorption independence of size at larger diameters. Fine-grained snow showed much greater scatter, suggesting there was greater variability of absorption capabilities with different samples than for coarse-grained snow. This could be because snow collection of fine-grained samples took place over many months (November–March of each boreal season), with much more potential for varying snow conditions, while collection of coarse-grained samples only occurred in March, allowing less variation.