CO2 sublimation in Martian gullies: laboratory experiments at varied slope angle and regolith grain sizes

Abstract Martian gullies were initially hypothesized to be carved by liquid water, due to their resemblance to gullies on Earth. Recent observations have highlighted significant sediment transport events occurring in Martian gullies at times and places where CO2 ice should be actively sublimating. Here we explore the role of CO2 sublimation in mobilizing sediment through laboratory simulation. In our previous experimental work, we reported the first observations of sediment slope movement triggered by the sublimation of CO2 frost. We used a Mars regolith simulant near the angle of repose. The current study extends our previous work by including two additional substrates, fine and coarse sand, and by testing slope angles down to 10°. We find that the Mars regolith simulant is active down to 17°, the fine sand is active only near the angle of repose and the coarse sand shows negligible movement. Using an analytical model, we show that under Martian gravity motion should be possible at even lower slope angles. We conclude that these mass-wasting processes could be involved in shaping Martian gullies at the present day and intriguingly the newly reported CO2-creep process could provide an alternative explanation for putative solifluction lobes on Mars.

When initially observed by the Mars Global Surveyor (Fig. 1), Martian gullies were reasoned to have been formed by flowing water (Malin & Edgett 2000).With the observation of active modification and extension of gullies (Dundas et al. 2012), the apparent lack of liquid water under current Martian surface conditions leads to increased interest in alternative formative mechanisms.CO 2 is abundant on Mars, both as a gas (accounting for c. 95% of the atmosphere; Bargery et al. 2011) and as solid ice at the surface.Although concentrated in seasonal polar ice caps (Hess et al. 1979), CO 2 ice has been detected on the surface at latitudes down to 32.3°S (Vincendon 2015).Like Earth, the climate on Mars is seasonal, owing to the tilt of its axis of rotation relative to its orbital plane (25.2°v. 23.4°for Earth) driving winter expansion and summer retreat of the seasonal ice caps as the ice sublimates into the atmosphere (Hess et al. 1979).
This work explores the ability of seasonal CO 2 frost to trigger slope failures of unconsolidated regolith.Previous experiments by Sylvest et al. (2016) established the ability of CO 2 frost to trigger mass wasting on slopes under Martian atmospheric conditions.Those experiments, however, were all conducted on slopes of JSC Mars-1 regolith simulant at approximately the angle of repose.This steep slope angle was chosen as it was the most likely to produce observable initial disruptions in a relatively short period of time.Yet clearly, Martian gullies and other mass wasting forms on Mars occur over a range of slope angles: for instance, gullies commonly occur on slopes ranging from 10°to 30° (Dickson et al. 2007;Conway et al. 2015).JSC Mars-1 is a finegrained volcanic tephra ranging in size from 13 µm to 1 mm (Allen et al. 1998), yet on Mars there is an exceedingly wide range of grain size distributions to be found in the soils at the planet's surface, ranging from very fine sand to coarse sand, together with a fine dust component (McGlynn et al. 2011).
While there have been Mars-focused experimental studies undertaken on the role of slope on the morphological characteristics of slope disturbances under both Earth surface conditions (Coleman et al. 2009;Jouannic et al. 2015) and Martian conditions (Jouannic et al. 2015), no studies have systematically examined the role of slope over a range of angles known to support a variety of mass wasting forms, including gullies.Similarly, a few experimental studies have systematically investigated the role of grain size in influencing the movement of debris under Martian atmospheric conditions (Conway et al. 2011a).These studies, however, examined water flow rather than CO 2 frost sublimation.
The experiments reported in this paper explore the influence of sediment grain size and initial slope angle on the ability of sublimating CO 2 ice to trigger mass wasting under Martian conditions.

Methods
A total of 28 experimental runs, including two control runs, were performed in the Mars Simulation Chamber (MSC) at the Open University, Milton Keynes, UK (Fig. 2).This cylindrical vacuum chamber, 2 m long by 1 m in diameter, is capable of replicating Martian atmospheric temperatures and pressures.Three regolith simulants were used at several initial slope angles, summarized in Table 1.
Each simulant was first tested at the angle of repose (AOR), as this is the least stable initial slope angle possible, and therefore the most likely condition to support the triggering of mass wasting (Sylvest et al. 2016).Subsequent angles were then selected based on the observed results of each previous run for that simulant.
The apparatus and procedures used for this work were adapted from Sylvest et al. (2016).Each run required preparation of the initial slope, preparation of the chamber atmosphere, cooling of the slope, condensation of CO 2 frost and, finally, sublimation of the frost.The only modifications to the methods of Sylvest et al. (2016) were in the initial slope preparation.Below, we present the slope preparation procedure and briefly summarize the remaining procedures.
Slope preparation started with vacuum drying the regolith, followed by stirring to disaggregate any clumps and to help maintain a consistent initial bulk density and grain sorting between runs (Sylvest et al. 2016).The initial slope angle was then prepared by raking the regolith into a smooth, uniform slope with a plastic spatula.The angle was checked with a spirit-level protractor at three locations across the crest slope zone (Fig. 4a).The nominal slope angle for the run (Table 1) was based on the crest slope zone.The same c.30 cm long, 20 cm wide, 12 cm deep test section from Sylvest et al. (2016) was reused.
The MSC was initially evacuated and purged with N 2 to remove moisture and provide a noncondensing atmosphere for the regolith cooling procedure.The regolith was cooled by flowing liquid N 2 through the test section coils (Sylvest et al. 2016) until the maximum near-surface temperature was cold enough to ensure condensation (−120°C, Table A1) of the relatively warm CO 2 gas (c.−20°C) at reduced chamber pressure (c.350 mbar, Table A1).The chamber pressure typically reduced to c. 160 mbar during the cooling procedure, without pumping.
During the condensation procedure, CO 2 gas, cooled in a heat exchanger by the exhausted liquid N 2 (Fig. 2), was gently diffused (Figs 3d & 4b) over the regolith slope with the test section lid closed (Figs 3f & 4b).By keeping the lid closed, the CO 2 Fig. 2. The Large Mars Planetary Environmental Simulations Chamber (centre), vacuum pump (bottom), liquid nitrogen supply (left) and CO 2 gas cooler (front left of centre).frost was prevented from forming on the horizontal photogrammetric markers which were glued to aluminium mounting surfaces.
The chamber pressure was pumped down to a representative Mars surface pressure of 5-7 mbar (Table A1; Hess et al. 1980) during the final 20 min of the frost condensation procedure.
Finally, the test section lid was opened and the frosted regolith surface exposed to a simulated insolation of c. 350 W m −2 maximum heat flux (Table A1) while recording the surface with stereo videography for photogrammetric analyses (Fig. 4b; Sylvest et al. 2016).
For each run, initial and terminal slope angles were recorded manually.Temperatures were recorded continuously at several depths in two locations along the centreline of the slope (Fig. 4a) and are summarized in Table A2.Pressure within the MSC was logged manually throughout each run (Sylvest et al. 2016).

Regolith simulants
JSC Mars-1, fine sand and coarse sand were selected as regolith simulants to identify possible controls on slope stability related to CO 2 ice sublimation, primarily owing to grain size distributions, but also to regolith composition.JSC Mars-1, a volcanic tephra, is the finest grained of the three simulants, with a broad grain size distribution from 13 µm to 1 mm (Allen et al. 1998).The fine sand regolith simulant is a well-sorted, dry silica sand with particle sizes between 100 and 300 µm.Based on sieve analysis, the coarse sand is poorly sorted, with particle sizes between 3 mm and 125 µm.The hydraulic properties pertinent to slope stability and fluid dynamic analyses of all three simulants are listed in Table 1.Grain size distributions are presented in Appendix A (Figs A1 -A3).

Photogrammetric methods
The photogrammetric methods used for this study are an extension of those developed by Sylvest et al. (2016).As in that work, two identical Sony HDR-CX330 camcorders were used to simultaneously and stereoscopically record slope surface activity in high definition (Fig. 4b).Coded photogrammetric targets were precisely located at multiple levels and orientations on the cooling box (Fig. 3e).These targets provide external control for all aspects of the photogrammetric project, establishing a common three-dimensional frame of reference (and therefore scale) for all aspects of the photogrammetric models.The video recordings were then processed (Fig. 5) to produce quantitative data.Each pair of videos was synchronized, and then image pairs were captured at 10 min time steps using Sony Vegas Pro 13 (Sylvest et al. 2016).The image pairs were imported into Agisoft PhotoScan Pro 1.2.6 and a series of digital elevation models (DEMs) was generated, one DEM for each time step.
There are two primary steps involved with the DEM production.The first step involves applying a photogrammetric bundle adjustment to determine interior (i.e.focal length, principal point and lens distortions) and exterior (i.e.camera positions and orientations) camera parameters.This step will henceforth be referred to as image alignment.The second step is the extraction of a dense 3D surface using multi-view stereo techniques.It is important to note that the two primary camera positions provided full stereoscopic coverage for all coded targets and the experimental surface.However, typically more than two camera positions are necessary to photogrammetrically model the above-mentioned camera parameters.Therefore, in order to improve the overall geometry of the photogrammetric block, the two primary camera positions were supplemented with images recorded from various positions around the box before the start of cooling, thus providing a much larger, virtual set of cameras.The use of these additional camera positions and coded targets with known coordinates (i.e.control points) provides a strong photogrammetric network with a reprojection error of less than one pixel.
During the image alignment step in PhotoScan, the regolith surface was masked out of each image and a sparse set of 3D tie points was generated using the 'Align Photos' tool.This process involves automatically detecting and matching corresponding features across overlapping images.The coded targets were also detected by PhotoScan and the known coordinates for these targets were provided.In some instances, PhotoScan incorrectly measured and/or labelled the coded targets and manual adjustments of the marker(s) were necessary (Fig. 5).After necessary marker corrections were made, the photogrammetric block was then optimized using the bundle adjustment.Next, a dense 3D surface was extracted from the primary image pair and the resulting surface geometry compared with the known dimensions of the box.If the model was unsatisfactory, the marker locations for each coded target were  manually refined in each image, and additional tie points could be added, based on the nature of the error.This process was repeated (Fig. 5) until consecutive results were unimproved.
After the image alignment procedure, a sequence of DEMs (e.g.Fig. 7b) at 1 mm/pixel resolution were produced (one for each pair of images), now masking all but the slope inside the box.Likewise, a sequence of corresponding orthophotos at c. 0.29 mm/ pixel (e.g.Fig. 7a) was produced.Rasters of the elevation differences between consecutive DEMs were calculated in ESRI ArcMap 10.3.1.The region of interest (ROI) for each sequence of DEMs was determined within ArcMap by combining all the orthophotos for the sequence into a single shape that encapsulates the full extents of all the DEMs.
The difference rasters were multiplied by the DEM cell size (1 × 1 mm) to determine the erosional and depositional volumes.The volume of frost condensed on each slope was taken as the net deposition between the pre-condensation time step and the start of the sublimation procedure.The DEMs were also used to generate long profiles, used to characterize slope angles.The long profiles were generated by constructing a ROI 4 cm wide along the length of the slope in each DEM, and averaging the heights every 4 mm along the slope.Linear least squares fits were applied to the binned data in the three zones, crest, mid and base, as indicated in Figure 4a.
In order to reduce systematic volumetric noise artefacts in the DEMs, a new filtering procedure was added to the original methodology of Sylvest et al. (2016).For each difference raster, representing the change elevations between two DEMs, ArcMap was used to find contiguous volumes of transported regolith.A mask was generated by eliminating all but those volumes with areal extents below an interactively determined threshold, typically around 40 cm 2 .This volume filter mask was then applied to the original difference raster, thereby reducing the photogrammetric noise.The filtered results were then imported into MathWorks Matlab 2016a in order to calculate slope angles and changes of angle for each slope zone (Fig. 4b), as well as to tabulate and plot the data.

Control runs and error estimation
The list of potential sources of error in the photogrammetric results, and the difficulty in assessing the contribution of any one of these sources, renders a direct error computation impractical (Raack et al. 2017).In order to quantify the 'uncertainty envelope' on a case-by-case basis, we evaluated differences in DEMs for slope regions known (by observation) to be devoid of sediment movement (as detailed later in this section).Table A3 lists the potential sources of error we have identified with our methods.
Control runs were performed for JSC Mars-1 and fine sand.A control run was deemed unnecessary for coarse sand, given the lack of substantial sediment movement observed for those runs.No sediment transport was observed for either of the two control runs.The control runs were configured at an initial slope at the static AOR and followed the same protocol used for the data runs, except that no CO 2 was introduced into the chamber during the condensation phase of the experiment.(See Sylvest et al. 2016 for full details of the experimental procedures.) We did not use the control runs to calculate the errors on volumes and slopes, as potential sources of error, such as those listed in Table 4, varied between runs.For example, dense accumulations of surface frost are particularly difficult to model under the harsh lighting of the simulated insolation.The direction and intensity of the lamp (Fig. 4b) tends to eliminate shadows necessary to detect the surface features, which in turn are required for the dense surface reconstruction.Also, the high albedo of the frost tends to oversaturate the camera sensor, further obscuring features on the surface of the frost.Hence, during runs with more extensive surface frost coverage, noise in the surface model was more significant than for other runs.To allow for these variabilities, we developed a method to estimate errors using the difference DEMs for each run individually.Specifically, for each run, a small ROI (4.67 cm 2 ) was located on a portion of the slope where no movement was visually observed.The vertical noise was estimated as the mean change of elevation over the ROI, when comparing the first and last DEMs of the run.As in Sylvest et al. (2016), horizontal coordinate error estimates were based on the RMS errors for the coded target locations reported by PhotoScan, which were c. 1.0 mm for all runs.The equivalent vertical error for each of the runs, derived from the noise estimation procedure described above, ranged from 0.43 to 2.29 mm (Appendix A, Table A5).As in Sylvest et al. (2016), the standard accumulation of errors formulae were applied to these numbers in order to provide error estimates for each measured result.

Results
In the following sections, we describe the shared features of all the experiments, including how these volumes of sediment were transported, their magnitudes and the observed changes in slope angle.Then, we discuss the results for each sediment type separately.

Quantitative results
Substantial volumes of sediment were transported for JSC Mars-1 at initial slope angles from 17.5°t o the AOR (c.30°) and for fine sand at the AOR (Table 2).Between 1% and 14% of the total sediment volume was transported at initial target slope angles between 17.5°and the static angle of repose (AOR c. 30°) for JSC Mars-1, while between 2% and 4% of the total sediment volume was transported at the AOR for fine sand, also c. 30°(Table 3).Measured volumes of observed sediment transport for coarse sand at the AOR were less than 1% of the slope volume, well below the noise level of the photogrammetric results of 1.3% of the total slope volume (run 26, Table 5).Owing to the lack of sediment transport for coarse sand at the AOR, no additional runs were conducted with gentler slopes (Table 2).The maximum volumes eroded and deposited for all runs were 10.9% and 14.0% of the total slope volume, respectively, both recorded for run 24 (Table 3).These volumes correspond to maximum changes in slope elevation 47 mm for erosion, and 42 mm for deposition (Table 3).
Initial and final slope angles for each slope zone (Fig. 4a) are reported in Table 3, along with their associated angular changes.For the majority of the 25 runs with photogrammetric data, the angles of the crest and mid-slope zones changed very little (−1.1 ± 0.4°and −1.3 ± 0.4°on average, respectively).Only two JSC Mars-1 runs (runs 20 and 24, Table 3) underwent crest slope angle reductions greater than 8°.Base slope zone angles reduced by 3.6 ± 0.3°on average.No statistically significant correlations were found between initial slope angle or initial frost volume with elevation, volumetric or angular changes (see Appendix A for correlation results).

Sediment transport types
Sediment transport was observed in 20 of the 26 experimental runs conducted.Four types of sediment transport were identified in various combinations for each of the runs in Table 4.The most readily observed type is discrete granular flow.These flows are up to several centimetres wide, and can extend as far as the base of the slope, where they form depositional fans and lobes.They generally display well-defined areal extents, and occur on all but the gentlest slopes.When occurring concurrently, discrete flows frequently coalesce into larger composite flows, occasionally encompassing the full breadth and majority of the length of the slope.These flows can initiate on slope regions with no topographically induced instabilities, but only initiate on regions devoid of visible surface frost.Discrete granular flows were present in 19 of the 20 runs in which sediment transport was observed and recorded.
Second, in terms of volume transport, creep movements are more difficult to observe, generally requiring accelerated video playback to see their subtler movements.However, the areal extents of some creep movements were large enough to suggest substantial sediment transport volumes.Fourteen of the 26 runs were affected by creep (Table 4).Creep was not confined to any particular slope region; however, it was frequently present at the boundaries of retreating surface ice accumulations.These movements presented as narrow (<1 cm wide) bands of sediment, adjacent to and following the trailing edge of the retreating ice.
Sixteen of the 26 runs displayed obvious signs of gas entrainment of regolith particles.Although a common phenomenon, entrainment did not appear to be an effective sediment transport process on its own.Entrainment was most obvious at the toe of the slope, where sediment was deposited over surface frost.In these instances, the dusty surface gave the impression of boiling, owing to the vigorous sublimation of the shallowly buried frost.Evidence of escaping gas was also frequently observed where the thermocouple stations came close to, or in some cases became exposed at, the slope surface.These circular areas of activity were typically smaller and less active than the active areas at the slope base.For runs with JSC Mars-1 at initial slope angles less than the angle of repose, fine dust deposits on the upslope box boundaries suggest that dust particles were entrained in the flow of CO 2 gas escaping from under the regolith along the warming box edge.The fourth sediment transport type was the rolling and tumbling of individual, isolated grains; henceforth simply referred to as tumbling.Although tumbling grains of frost were common for many of the runs, they had no apparent effect on the underlying regolith; nor were they observed to trigger other modes of transport.Tumbling regolith particles were only observed in the two coarse sand runs (Table 4).This may be due in part to the much smaller, more uniform grain sizes of JSC Mars-1 and fine sand (Table 1).As for creep, accelerated video playback speeds facilitate observation of these small (1-2 particle diameters) movements which occurred over the entire extent of the slope.As for tumbling ice particles, tumbling sediment neither moved substantial volumes of material, nor did it trigger other sediment movements

JSC Mars-1 sediment movement types
JSC Mars-1 was the most active of the three regolith simulants.Discrete granular flows were detected for all but the one run with the gentlest initial slope angle (11.4 ± 0.1°for run 9 in Table 4).Creep was detected in 10 of the 16 runs with JSC Mars-1; and gas entrainment was observed in all but one run (run 12 in Table 4).Correlating sediment movements (Table 4) with the initial slope angle of the corresponding slope zone (Fig. 4a), discrete flow initiation was identified on slopes between 12.2°and 39.5°, while creep initiated on slopes between 11.4°and 33.9°(Table 3).The observed behaviour of discrete granular flows revealed a dependence on initial slope angle.For JSC Mars-1, run 9 (Table 4), with an initial slope angle of c. 10°, was the only run that did not display any discrete flow events.Flows became both larger and more numerous with increasing steeper slope angles.The maximum width of flows increased with initial slope angle from c. 1 cm at 15°slopes, and reaching the full slope width (c.20 cm) for initial slopes of 20°.All runs starting at the AOR displayed discrete flows up to the full width of the slope, with runouts from crest to the bottom limit of the box.Maximum runouts started at less than 2 cm for 15°slopes and increasing to 6 cm at 17.5°initial slopes.For initial slopes below 20°, flows primarily took the form of slumps, leaving an upslope recess and a downslope berm of displaced sediment just below the recess.At steeper angles, flows appeared to skim, rather than slide over, the downslope surface, particularly flows starting at the AOR.For these steepest runs, flows frequently initiated at the crest, leaving a sharp break of slope, which receded with continued flow activity.
As seen in Table 4, creep was detected in runs covering the full span of tested initial slope angles (10°to c. 30°).Conspicuously, creep was not detected for any of the runs starting at 20°.With the exception for run 9 (Table 4), creep was not detected in the crest zone (Fig. 4a), and in that instance, was only present at the edges, adjacent to the box.Gas entrainment activity was relatively independent of slope angle, except where it modified deposits from previously transported sediment.

JSC Mars-1 morphology observations
Slope failures on JSC Mars-1 slopes produced a variety of morphological changes on the surface.Discrete slope failures leave clearly defined areas of erosion and deposition, with distinct boundaries (Fig. 6).While the steepness of crest regions near the top of the slope did generally decrease, sharp breaks in slope were maintained through crest retreat, observed in several runs.Scarp-like features were frequently observed on eroded slope faces (Fig. 6a, b).For larger volumes, discrete slope failures, eroded material was carried to the bottom of the slope, colliding with the end of the test section with sufficient energy to create large clouds of dust, briefly obscuring the entire slope.Centimetre-scale puffs of dust, entrained by jets of escaping subsurface gas, were also observed.The timing of these events is difficult to constrain, given their brief duration and small scale, combined with near-vertical camera angles.Broadly, they occurred within the periods of gas entrainment activity listed in Table 4, beginning as early as the onset of sublimation up to halfway (48%) through the period, and ending from 32% to 98% through the sublimation period (Table 4).Jets adjacent to the box sides left remnant fans of dust just above the slope on the box sides (Fig. 6c).Escaping gas also left behind relatively large areas of pitting with millimetre-scale pores, principally at the foot of the slope where surface frost was buried by deposition (Fig. 6e), and also centred above one or both thermocouple trees (Figs 4a & 6d).
Looking at the difference raster on the right in Figure 7b, the greatest amount of erosion was on the right-hand side of the slope, at and below the original crestline, visible as the darker orange region.Deposition is deepest at the base slope zone (dark blue) and extends into the mid-slope zone, relatively symmetrically along the centreline.The difference raster on the left in Figure 8b shows the accumulation of frost at t = 0 min, the start of sublimation.However, over the majority of the slope, the thickness of the condensed frost is on the order of the photogrammetric noise.Hence, the pattern of frost, clearly visible in the corresponding orthophoto at t = 0 min (Fig. 7a), is not readily visible in the difference raster.The crest of the slope remains well defined and retreats, as seen in the orthophotos (Fig. 7a) and the long profiles (Fig. 7c).A small, arcuate alcove is clearly visible in the centre of the crest, along with lobate debris aprons downslope at t = 80 min.(Fig. 7a, b).
By way of comparison, Figure 8a-c illustrates the slope evolution for run 14 (JSC Mars-1 at 17.5°), the gentlest slope tested with substantial regolith movement (Table 2).In Figure 8a & c, we see that, unlike run 20 in Figure 7a-c, the initial slope is relatively smooth, with no break of slope at the crest, and curvature at the toe of the slope is confined to the bottom c. 30 mm (Fig. 8c), v. c. 100 mm for run 20 (Fig. 7c).The colour classifications representing slope elevation changes are repeated from Figure 8b to illustrate the relative magnitudes of slope modifications between the two runs.By comparison, all of the changes in run 14 are near the error estimates for elevation change of ±0.69 mm (Table 3).Although not readily apparent in Figure 8a, visual observation of the video recordings confirms the location and approximate magnitude of the long profiles in Figure 8c.Where essentially the entire slope was modified in the steeper run 20, only the lower half of run 14 underwent substantial alteration.
Quantitative results for JSC Mars-1 JSC Mars-1 was the most active of the three regolith types tested, undergoing the largest changes in slope elevations, volumetric displacement and the largest slope angle changes for all three slope zones (base, mid-slope and crest; Table 3).The maximum elevation change associated with erosion was 47 mm, with a mean of 2.3 mm; and the maximum change owing to deposition was 42 mm, with a mean of 3.0 mm (Table 3).The limit of detection for changes in slope height (ΔZ ) is 1.4 mm, based on the noise estimates described under Controls and Error Estimation.The maximum erosional volume was 374 cm 3 , with a mean of 60 cm 3 , and the maximum depositional volume was 482 cm 3 , with a mean of 82 cm 3 (Table 3).The mean decrease in base slope zone angle was 4.9°, with a maximum of 15°.Both the mid-slope and crest zones had a mean decrease in slope of 1.6°, with maximum reductions of 13°a nd 15°, respectively (Table 3).

Fine sand sediment movement types
Fine sand transport was only detected in the three runs starting at the AOR (Table 4).For the first run (run 6 in Table 4), sediment movement was almost imperceptible.The cumulative magnitude of the by guest on April 23, 2018 http://sp.lyellcollection.org/Downloaded from minute creep observed along the retreating boundaries of surface frost was only revealed in the photogrammetric results (Table 3).For the second run (run 7 in Table 4), slow, steady creep eventually led to a large, slope-wide discrete flow.This was the only run with fine sand that exhibited gas entrainment, which in this case was observed as bubbling of the sediment deposited atop surface frost at the toe of the slope.This event which lasted c. 139 s, was an immediate consequence of a discrete flow, which over-topped an accumulation of surface frost at the toe of the slope, just over halfway through the c. 2 h run (Table 4).For the third run (run 8 in Table 4), a single large discrete flow was accompanied by creep, which was most evident in the reshaping of the depositional fan generated by the initial, discrete flow event.For these three runs, the initial slope angle was photogrammetrically measured between 27.6°and 30.5°(Table 3).Discrete flows and creep were correlated with photogrammetrically determined slope angles, each between 29°and 35°( Tables 3 & 4).

Fine sand morphology observations
Morphological changes of fine sand slopes were restricted to erosional lowering of steeper slope segments and development of depositional fans.Run 8, the most active of the fine sand runs, underwent the greatest slope angle reduction for all three slope zones (Fig. 4a and Table 3).As illustrated in Figure 9c, roughly the top third of the long profile was eroded, while the balance of the slope accumulated the associated deposition.Both the pre-frost and t = 100 min profiles highlight the somewhat featureless character of the sand slope.The t = 100 min orthophoto in Figure 9a does reveal somewhat extensive depositional features over much of the slope.The t = 100 min profile, recorded along the slope centreline, fails to capture the full depth of erosion visible at the left-hand side of the crest in the difference raster (Fig. 9b).Noting that interference of the box at the toe of the slope limits interpretation of runout features, depositional material did, in some instances, accumulate along the upslope boundary of dense surface frost without reaching the box boundary.Subsequent sublimation of this frost left steepened lobate termini.
In comparison with the JSC Mars-1 runs, fine sand slopes remained relatively smooth, with no sharp breaks or scarps.Fine sand did, in contrast, develop more complex depositional features, fans frequently colliding and coalescing (Fig. 9a).Entrainment features were much less evident than for JSC Mars-1, only one entrainment event being observed.3).(c) Dust fans along back wall of box are the result of centimetre-scale jets of escaping CO 2 gas carrying entrained dust (from run 18 Table 3).(d) Pitting around the mid-slope thermocouple tree (upper-most thermocouple has been exposed at the surface), caused by escaping CO 2 gas (from run 15 Table 3).(e) Pitting at the base slope zone from run 24 (Table 3).

Fine sand quantitative results
For the three active fine sand runs (all at the angle of repose, Table 2), the maximum change in slope elevation associated with erosion was 12 mm, with a mean of 9.0 mm, compared with a maximum of 47 mm and mean of 10 mm for JSC Mars-1.The maximum change associated with deposition was 23 mm, with a mean of 15 mm, compared with a maximum of 42 mm and mean of 12 mm for JSC Mars-1 (Table 3).The fine sand slopes underwent 68.7-88.8cm 3 of erosion, with a mean of 77.8 cm 3 .Deposition was between 1.55 and 132.0 cm 3 , with a mean of 59.8 cm 3 (Table 3).The mean change of angle for the base slope zone was −3.9°.The midslope and crest zones were essentially unchanged, the measured values (−1.3°and 0.5°respectively) being well within the estimated error of ±1.4°( Table 3).The largest slope angle changes were 4.6°, −3.4°and −5.75°for the crest, mid-slope and base regions respectivelyall of these values coming from run 8 (Table 3).

Coarse sand sediment movement types
Coarse sand was the least active of the three regolith simulants tested.Only two runs, both at the AOR were conducted with coarse sand, as no appreciable sediment transport was observed.Poor camera placement prevented development of a usable photogrammetric model for run 25; however, minute surface alterations were observed, visually.In that run, a single, shallow (<2 mm deep) discrete sediment flow (c. 10 mm wide, 20 mm long) was observed, apparently driven downslope by a similarly small mass of surface frost which broke away from the side of the test section.Creep can be detected along the retreating edges of surface frost in the video recordings.Larger sediment grains can also be seen to roll or tumble in the video recordings.While this activity is more easily observed than the creep, it does not appear to cause any substantial slope movement.
A single discrete slope failure was observed in run 26, over halfway through the sublimation process; in this case without the influence of frost falling from the box side.The extent of the sediment transport was comparable with that in run 25, and had a similarly limited effect on the surface morphology.As for run 25, rolling and tumbling of larger sediment grains was active throughout the duration of the run, widely distributed across the entire slope.Based on photogrammetric results, the discrete regolith flow was detected between 30°and 32°, and creep was detected at 30°.

Coarse sand morphological and quantitative results
Only superficial morphological changes were observed for both coarse sand runs.Both runs displayed the tumbling of individual sediment grains.These grains were visually estimated to comprise <5% of all surface grains, distributed widely and uniformly across the entire slope area.These sub-centimetre-scale movements produced no discernible surface features on either of the two coarse sand slopes.Creep was visible across the mid and crest slope zones for run 25, with somewhat more evident movement along the box boundaries.Creep was less widespread in run 26 than in run 25, only apparent along the box edges and in association with the retreating edge of surface frost accumulations.For both runs, creep-induced movements were too small to produce photogrammetrically detectable surface alterations.
Coarse sand, initially at the nominal angle of repose, showed no substantial changes of elevation: on average <1 mm for the entire slope (Table 3), which is less than the estimated vertical photogrammetric noise of 1.2 mm (Table 3).The measured erosion for run 26 was 10.0 cm 3 (0.3% of the total regolith volume), while deposition was 32.2 cm 3 (0.9% of the total regolith volume) and the volume of frost was 64.8 cm 3 (Table 3).To put these volumes in context, the erosion and deposition were 0.3% and 0.9% of the total regolith volume, respectively.Remaining surface frost, visible in the final image pair for run 26, is included in the reported deposition volume.The mean change of slope angle was −3.21 ± 0.73°for the base slope zone, −0.6 ± 0.66°f or the mid-slope zone and −3.4 ± 0.37°for the crest zone (Table 3).Sylvest et al. (2016) proposed that the observed mass-wasting behaviour of slopes of JSC Mars-1, initially at the angle of repose, was triggered by the sublimation of CO 2 frost condensed within the regolith pore space.They hypothesized that the rapid production of gas produced via sublimation caused the pore pressure within the sediment to increase and therefore initiated failure.In this work, we have tested two new substrate types and different initial slope angles, and we found:

Discussion
(1) volumes of sediment moved by JSC Mars-1 remaining similar down to slope angles of 20°(and are slightly reduced at 17°) and volumes of sediment moving in fine sand experiments of the same magnitude, but only near the angle of repose (little movement was detected at other initial angles) and only very limited sediment transport occurring for coarse sand near the angle of repose; (2) four different types of sediment movement, discrete flow, creep, gas entrainment and grain tumbling (creep was not reported by Sylvest et al. 2016).
In the following discussion, we integrate our experimental results into a discussion of the physics of these movements and then discuss their applicability to Mars.Our discussion focuses on the discrete and creep flows, as these cause the most sediment transport.

Mechanism and physics
In order to assess whether the mass wasting we observe in our experiments is caused by the reduction in friction angle of the sediment owing to the gas flow reducing the intergranular pressure, we consider a simple one-dimensional analytical model.The model solves continuity equations for energy and CO 2 , and balances downslope gravitational attraction against Coulomb friction.Suppose that the bed comprises sand with bulk density ρ s , thermal conductivity κ, permeability k, specific heat capacity c and some initial amount of CO 2 ice with density ρ i .We assume that at t = 0 the bed is all at the sublimation temperature T s , and then a radiant heat flux Q is applied to the surface.We measure distance downwards normal to the surface using the coordinate x, and define the point X(t) as the boundary between CO 2 ice and pure sand.Initially X(0) = 0, but over time, as heat is conducted into the bed, this point will move downwards (increasing x).We assume that the heat flux is due only to heat conduction between the sand grains; thus the temperature T(t, x) satisfies the equation where the subscripts t and x are partial derivatives with respect to time and position.The boundary conditions are representing the heat flux at the surface, the temperature being at the sublimation temperature at the CO 2 ice interface.The movement of the point X(t) is given by a differential equation that describes the heat flux driving the sublimation of the CO 2 ice: where e is the enthalpy of sublimation.We assume that temperature variations do not significantly affect any material properties, including the carbon dioxide gas density ρ g .The behaviour of the solution is best understood by defining the following time, length and temperatures scales, t * = kr 2 i e 2 /r s cQ 2 , x * = (kr i e/r s cQ) and T * = (r i e/r s c).The model is then non-dimensionalized by writing and where f and g are non-dimensional functions.The resulting non-dimensional system and an approximate analytic solution are described in Appendix B. A key point is that, for time less than the order of t*, most of the heat is going into sublimating ice, and the frost front advances linearly.Conversely, for time greater than t*, the energy balance changes and most of the heat goes into warming ice-free sand, the flux decaying as 1/ t √ .The largest gas flux, and hence highest grain mobility, occurs for short values of time.A power series solution for f and g can be developed in nondimensionalized time, which gives the frost point advance as Now, the mass production rate of CO 2 is r i Ẋ, and we assume that the gas density ρ g is constant and equal to the value at the sublimation temperature T s and background pressure p 0 .The volume flux of CO 2 for x < X(t) is therefore where the subscript s denotes the derivative.Within the bed, Darcy's law gives the stress on the sand grains as S = (n/k)q(t).If we define the nondimensional Darcy stress, S * = (n/k)(r i /r g )(x * /t * ) (1/r s g) = (nQ/ker g r s g), then S = ρ s gS*f s , S being the reduction in normal stress due to the gas flow.
As the maximum shear stress that can be supported by the grains is proportional to the normal stress between the grains, this reduction increases the probability of failure.The Coulomb failure criterion is independent of depth, and after dividing by xρ s g, at any point x < X(t), it is where μ is the coefficient of friction and θ is the slope angle.Thus, it can be seen that S* is the key nondimensional group that determines whether the bed is likely to mobilize, and that this does not depend on the concentration of ice.The ice concentration will, however, determine how long mobilization will occur, owing to the dependence on t*.For short times (t < t*), f s = 1; hence, the CO 2 flux can initially be calculated by assuming that all the radiant heat is subliming CO 2 ice.Solving equation ( 8) for θ, using the parameters for the experimental conditions, gives the new, reduced internal friction angle for the slope under those conditions.This can be solved exactly, but an approximate formula, for small q, is more convenient: The input experimental parameters are listed in Table B1.If the model provides an accurate description of the physics, then we would expect that the reduction in friction angle corresponds to the initial slopes at which we observe movement in our experiments.Note however, that there is considerable (as much as 5°) stochastic variation in the failure angle of granular materials.
The results of applying our analytical model to our experimental data are presented in Table 5.Because coarse sand shows no significant transport, even at the angle of repose, this substrate is expected to have the lowest reduction in friction angle.Conversely, because JSC Mars-1 transports substantial volumes of sediment on slopes down to 17°( Table 3), it is expected to undergo the highest reduction in friction angle.Fine sand should therefore be intermediate, between the other two sediments, as it only shows substantial volume transport near the angle of repose.The different substrates do have the expected relative reductions in friction angle predicted by the model (Table 5).The reduction in friction angles for coarse and fine sand under terrestrial gravity (0.5°and 1.9°respectively, Table 5) match our experimental results, whereby no appreciable movement is seen for coarse sand, and appreciable movement is only observed near the angle of repose for fine sand.We note that the finer grain fraction in the coarse sand did show some signs of mobilization; however, we infer that the presence of larger grains both impeded mass movement and also increased the substrate permeability, preventing the build-up of pressure within the pore space.For JSC Mars-1, we get a slightly larger change of internal friction angle than expected, of 17.6°.Based on our experimental results, this would imply activity in the 15°experiments, which was not observed.However, unlike fine and coarse sand, we did not make our own independent measurements of the permeability for the JSC Mars-1 material (Sizemore & Mellon 2008), and slight variations can change the outcome of the calculation.Equally, our measurements indicated that the bulk density of the JSC Mars-1 could vary by ±0.07 g cm −3 , even with similar preparation, which could also contribute to this discrepancy.The volume flux velocity for our experiments was based on the rate of sublimation dictated by the supplied radiant heat and thermodynamic properties of the CO 2 .As these parameters were held constant for all runs, the volume flux velocity was also constant, at 0.012 m s −1 (Table 5).
The scale time, t*, is the cross-over time, when radiant heat input transitions from primarily sublimating CO 2 frost to primarily heating the regolith.In practice, t* is neither readily observed nor measured, but we might expect it to have the same order of magnitude as the duration of the phase   5).In contrast to fine sand, the behaviour of JSC Mars-1 is different (Fig. 10b).In this case, all six temperature traces are essentially constant at the onset of the sublimation process, followed by a steady increase.None of the JSC Mars-1 traces displays a plateau suggesting active frost sublimation, hence we cannot compare these data with the sublimation times of the fine sand runs, nor with the predicted t* values from our analytical results.
An additional complication for understanding the experiments is the complicated geometry and initial conditions.As shown in Figure 10, the initial temperature profile is far from uniform, and the entire bed is not at the sublimation temperature as was assumed in the analytical model.To investigate this further, we numerically solved the one-dimensional heat conduction equation with sublimation.We assumed zero heat flux on the lower boundary, and the same radiant heat flux used in the analytical model on the upper boundary (350 W m −2 , Table A1).The frost load was estimated to be 20 kg m −3 .We attempted a simple match of the initial temperature profile to the thermocouple readings.The results can be seen in Figure 10, along with the experimental traces.We see reasonable agreement between the experimental and numerical traces for JSC Mars-1 (Fig. 10b, d).The traces for fine sand (Fig. 10a, c) do not match, which suggests that, although the theory is reasonable, the uncertainties in the initial and boundary conditions substantially affect the energy budget, and thus the sublimation rate.
These null results highlight that the predictive power of these physical models is limited by the quality of the measured parameters, particularly the regolith properties.Thermal conductivity, for example, can vary by an order of magnitude depending on how the regolith was handled during preparation of the slope.JSC Mars-1 provides particularly difficult challenges for modelling, owing to the increased influence of inter-particle interactions compared with even fine sand.However, even with these caveats, we are able to obtain results which capture certain aspects of our experimental results, and that provide a better understanding of the basic physical mechanisms driving the sediment transport in our experiments.This then allows us to transfer our results to Mars, which we present in the following section.

Application to Mars
One of the underlying assumptions in the analytical models outlined above was that CO 2 ice is emplaced below or within the sediment, in order for sediment transport to occur.This reflects a key observation from both the experiments reported here and in Sylvest et al. (2016): sediment transport is only triggered where CO 2 frost is in the subsurface, and not where it has accumulated on the surface.Reporting of CO 2 ice on Mars is limited to surface observations of frost and/or slab ice (e.g.Gardin et al. 2010;Appéré et al. 2011;Piqueux et al. 2015).However, when surface CO 2 is not present, trapping of CO 2 ice within the regolith should be possible, as has been previously formulated for the emplacement of water ice in the subsurface (e.g.Mellon et al. 1993).Assuming Mars' regolith is dry, the atmosphere can diffuse into the subsurface pore space.The surface of Mars undergoes a diurnal and annual temperature cycle, and therefore the physical processes described in 'Mechanism and physics' section could take place repeatedly on time and depth scales concurrent with diurnal, annual and secular temperature variations.As the surface temperature variation propagates into the regolith, it experiences a lagging and damping effect (for illustration, see Hagermann 2005).This means that there are periods in the subsurface thermal cycle when the temperature wave decreases to below the sublimation point of CO 2 such that the CO 2 sublimation horizon moves upwards.As atmospheric CO 2 diffuses into the soil, subsurface CO 2 deposition occurs.At low pressures, this is best described as molecules following the vapour pressure gradient until they are deposited as ice at depth.Our experiments reveal that the nature of the CO 2 frost, and exactly where it ends up on or in the slope, should control the type and quantity of sediment transport owing to subsequent sublimation of the frost.
Although beyond the scope of the current study, temperature gradients through the depth of the slope, during the condensation phase of the experiments, are expected to control the depth of subsurface frost formation, and possibly the density of both subsurface and surface frost.Our analytical models demonstrate that these are key factors in determining whether sediment transport will occur.
Given these arguments, we believe that the emplacement of CO 2 into the regolith pore-space should be possible on Mars, and therefore, we have applied our simple analytical model above using Martian gravity instead of terrestrial (Table 4).This reveals that the reduction in friction angle should be even greater on Mars for any given substrate, compared with the reductions observed on Earth.Hence, these results should be applicable to a larger range of slopes for any given substrate than expected, based on our laboratory results without scaling.The calculations under Martian gravity suggest that movements in fine sand should be observable on slopes as much as 5°lower than the nominal angle of repose.
In summary, our results predict that mass wasting (creep and discrete granular flows) can be triggered by the sublimation of CO 2 frost on Mars, where temperatures in the near-surface regolith dip below the condensation temperature of CO 2 , and where CO 2 surface frost or ice is not present at the surface when the sediment transport occurs.Our experiments do not allow us to directly compare morphologies observed on the Martian surface with morphologies observed in our experiments owing to the difference in scale.However, from our results we can infer likely locations on Mars where this process could occur by examining regions which have similar slopes and grainsizes to those found to be active in our experiments.Many of the present-day sediment transport events associated with Martian gullies fall into this category.Notable exceptions include: (1) polar pit gullies at 68°S, where Raack et al. (2015) noted that the recent dark flows occurred in the spring, when the surface is still covered with CO 2 ice; (2) activity in large apron gullies on sand dunes between 40°S and 60°S, which occur in winter, when CO 2 frost is still present (Diniega et al. 2010;Pasquon et al., this volume, in review); and (3) mass wasting events on north polar dunes, that are active in the mid-winter, under the CO 2 ice slab (which although not strictly gullies, could represent processes active in gullies; Diniega et al. 2017).These three exceptions represent a small proportion of known active gullies.The remaining active gully sites comprise 'classic', mid-latitude gullies and linear dune gullies (Auld & Dixon 2016;Conway et al. 2017).As reported in Vincendon (2015) and Dundas et al. (2017), activity in mid-latitude classic gullies is limited to periods when thin (microns to centimetres) and patchy CO 2 or H 2 O frost is present during, or just prior to the noted activity.So far, c. 67 such sites have been catalogued as active (Dundas et al. 2017), which is c. 18% of monitored sites in the southern hemisphere, but activity is sporadic, rather than annual.As reported in Reiss & Jaumann (2003), Reiss et al. (2010), Pasquon et al. (2016;this volume, in review and Jouannic et al. (this volume, in review), activity in linear dune gullies happens when the last CO 2 ice is disappearing from the pole-facing crest of the dune.Linear dune gullies are particularly active, with most of the 33 sites showing annual changes.Hence, in terms of timing and frost observations for both active 'classic' gullies and active linear gullies, the mechanism we have observed in the laboratory could be at work.
We re-emphasize here that our experiments cannot tell us what role sublimation of subsurface CO 2 frost is playing in forming the morphologies of these gully types, and we note that it may only be a secondary process (see further discussion on this point below).However, we think it should be considered among the candidates for morphological changes in these gullies for the following reasons.The substrate type is better known for linear dune gullies than for classic gullies.Sand dunes on Mars have been investigated in situ by rovers, the most recent study revealing grains ranging between 50 and 350 µm with a mean size of 113 µm (Ewing et al. 2017).Activity in linear gullies is only found where the crest of the dune is at 20°or higher, but the changes themselves occur on slopes down to 5-10° (Pasquon et al. 2016;Jouannic et al., this volume, in review).Both the substrate type and slope angles used in our experiments are consistent with these data, hence our results would predict movement in this context if CO 2 is condensed in to the subsurface.For classic gullies, we must rely on orbital observations, with thermal inertia measurements suggesting that gullies reside in materials classed as unconsolidated 'sandto pebble-sized grains' (Reiss et al. 2009;Harrison et al. 2014).Source areas for recent motions in classic gullies are usually diffuse, suggesting remobilization of a loose surface cover (Dundas et al. 2017) and source areas of classic gullies range upwards from 20°( median 25°; Conway et al. 2015).Results from our laboratory work and application of our analytical model reveal that these grain sizes and slope angles should be able to sustain the CO 2 sublimationtriggered failures that we have investigated in our experiments.
The range of sediments that could be mobilized on Mars could be wider than that encompassed by the fine sand and JSC Mars1 sediments used in our laboratory work.Our results indicate that the fine component of the sediment plays an important role in triggering activity at angles lower than the angle of repose.On Mars, dust (generally accepted as grains with >30 µm diameter) is abundant on the surface and in the atmosphere (e.g.Christensen 1986;Tomasko et al. 1999).However, further experiments would be required to assess how important this factor might be in triggering movement.
We also consider it possible that the mechanism of CO 2 sublimation-triggered failures we observed in the laboratory could also have occurred under climate conditions different from those observed on Mars today, which may explain some of the sediment transport that has contributed to gully formation on timescales of millions of years (Reiss et al. 2004;Schon et al. 2009;de Haas et al. 2015a).However, further modelling work to understand the plausible temporal and spatial extent of this process at the present day would be needed in order to confidently extrapolate this process into the past, an endeavour beyond the scope of this present work.
Our experiments investigated the triggering of failures of unconsolidated materials, and we observed two models of subsequent transport, granular flow and creep.A subtype of recent gully activity termed 'bright flows' has been found to have a morphology consistent with a classic granular flow (Pelletier et al. 2008;Kolb et al. 2010).However, other features of recent (and past) activity, including the transport of metre-scale boulders, the formation of levees and lobate termini (e.g.Dundas et al. 2010Dundas et al. , 2014Dundas et al. , 2017;;Johnsson et al. 2014;de Haas et al. 2015b) require some viscosity or fluidization of the flow beyond that of a simple granular flow.Further, in linear gullies the new morphologies are complex, encompassing albedo changes, formation of pits, channels with levees and complex distributary networks on relatively low slopes (Diniega et al. 2013;Pasquon et al. 2016;this volume, in review;McKeown et al. 2017;Jouannic et al., this volume, in review).As discussed by Stewart & Nimmo (2002), CO 2 gas should dissipate too quickly to enable durable fluidization of a sublimationtriggered flow.However, the fluidization of CO 2 sublimation-triggered granular flows has not yet been studied in the laboratory and should therefore be a focus of future work in order to substantiate these calculations.Our observations that CO 2 sublimation can generate a creep-like movement in unconsolidated sediments is of particular relevance to Mars, because it could provide an explanation for lobate features (Gallagher & Balme 2011;Gallagher et al. 2011;Johnsson et al. 2012;Balme et al. 2013;Soare et al. 2016), which are often associated with Martian gullies.The closest terrestrial analogue for these features is solifluction lobes, which are uniquely associated with creep generated by freeze-thaw cycling of water in the ground, so our results provide a possible alternative that needs to be explored further.
Our laboratory work results only pertain to sediment transport of unconsolidated materials.The present-day observations of activity in classic gullies seem to only encompass the transport of unconsolidated sediments, as outlined above.However, these landforms, which have been dated to millions of years (Reiss et al. 2004;Schon et al. 2009;Johnsson et al. 2014), are incised into consolidated materials including both the ice-rich Latitude Dependent Mantle and bedrock (e.g.Dickson et al. 2015;de Haas et al. 2017).Erosion of the unconsolidated sublimation lag believed to be on top of the LDM should engender loss of the interstitial ice by sublimation, rendering further sediment available for transport (e.g.Pilorget & Forget 2016).However, for gully alcoves cut into bedrock, it remains an open question as to how material is weathered to render it transportable.It should be acknowledged that whether this is a primary feature of the gully-forming process is under debate (cf.de Haas et al. 2015a;Dickson et al. 2015).The relationship between this weathering and the action of CO 2 condensation-sublimation cycles is an area for future work.

Conclusions
We have experimentally investigated the effect of sublimating CO 2 on the downslope mass wasting of sediment under Martian atmospheric conditions.We tested three substrate types, fine sand (mean diameter 168 µm), coarse sand (mean diameter 594 µm) and a Mars regolith simulant (JSC Mars-1) over slope angles ranging from 10°to near the angle of repose, adding to the work of Sylvest et al. (2016), which only considered JSC Mars-1 near the angle of repose.We observed four principal movement types: discrete flows, creep, gas entrainment and grain tumbling.Of these, creep and grain tumbling were not reported in Sylvest et al. (2016).The observed sediment movement types were influenced both by initial slope angle and by the nature of the regolith.
We found that significant volumes of sediment were only transported by the discrete flows and creep movements.In fine sand, these processes were only active at slope angles near the angle of repose and the volumes transported were of the same order as those transported in experiments using JSC Mars-1 at the angle of repose (82 cm 3 mean erosion for a surface area of 473 cm 2 ).For JSC Mars-1, these processes continued to transport equivalent volumes of sediment down to 20°, slightly less at 17°and negligible amounts at slope angles <17°.In our previous work we hypothesized that mass wasting was triggered by a lowering of the static friction angle by gas escaping through the substrate from sublimation.We tested this hypothesis by constructing an analytical model describing the physics of this process and found that when applied to our experimental parameters this model successfully predicts the activity observed in our experiments.With this validated model, we were able to predict that, under Martian gravity for equivalent sediment types, mass wasting could be triggered at even lower slope angles (movement should be possible for coarse sand near angle of repose, 25°for fine sand and on any slope for JSC Mars-1).This model also reveals that the reduction in permeability in the JSC Mars-1 owing to the presence of fines could be the key parameter for explaining this substrate's enhanced activity range compared with the two sands.Our results suggest that the absolute amount of CO 2 in the subsurface may control the amount of sediment moved.Further the vertical temperature profiles from our experiments reveal that CO 2 ice emplacement in the subsurface, in terms of vertical distribution and density, is a key parameter to elucidate in order to better understand the limits of sediment transport by CO 2 sublimation.
On Mars, we find that the CO 2 sublimationtriggered mass movements observed in our experiments could be applicable for explaining some of the movements seen in present-day mid-latitude gullies and linear dune gullies.Specifically, our experiments and analytical model reveal that the grain sizes and slopes should be compatible with this type of motion.Some of the features, including movements of metre-scale boulders, levees and lobate termini, are features that would require further experimentation to determine if CO 2 sublimation can explain them.Finally, we present the first observations of a creep-like motion caused by CO 2 sublimation and this could provide a viable alternative to water-ice freeze-thaw cycles to explain the origin of lobate features often found in association with Martian gullies.
Funding was provided for MES by the Open University, Europlanet, the Royal Astronomical Society and the University of Arkansas Graduate School.Europlanet 2020 R.I. has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement no.654208.MRP was funded under STFC grant ST/P001262/1 and from the European Union's Horizon 2020 research and innovation programme under grant agreement no.633127 (UPWARDS).JNM was partly supported by NASA grant NNX14AO21G.AH was supported by the UKSA, grant number ST/J005304/1.We thank W. Johnston for his guidance and support with instrumentation and T. Ringrose for his collaboration and guidance in the laboratory.
The final equation comes from requiring that the mean error in the conduction equation is zero, or that the total heat input matches the sublimation energy and heat increase in the sand.That is This equation simplifies to an ODE in f(s) This equation cannot be solved exactly.It is like an Abel equation, but with a quartic, rather than a cubic, form.It can easily be solved numerically, or approximated numerically.We are primarily concerned with the derivative, f′, since this sets the gas flux.This is well approximated by and has the exact large and small s behaviour.The flux rate falls to half its initial value when s ≈ 2.5.With the numbers in Table B1 (Mars conditions), we get a 5.1°reduction of bed friction angle, θ, for fine sand and a 50.3°reduction for JSC Mars-1.The unphysically large reduction for JSC Mars-1 (greater than the angle of repose) implies that the sublimed CO 2 gas is capable of levitating grains of this regolith at any angle, which is consistent with the observed entrainment behaviour described above.The much smaller reduction of bed friction angle for fine sand is also consistent with our results, these slopes only failing at angles near the static angle of repose.

Fig. 1 .
Fig. 1.A HiRISE gully monitoring image of a series of typical Mars gullies with source alcoves converging downslope (a) into a transport channel (b), which opens onto extensive depositional fans (c).Image: ESP_048424_1105_red -NASA/JPL/University of Arizona.

Fig. 3 .
Fig. 3.A view of an experiment inside the Mars chamber, at the start of sublimation.(a) The copper-coil regolith cooling box, externally insulated with open-celled foam and Mylar.(b) Insulated liquid nitrogen supply and exhaust lines.(c) Thermocouple wires were routed along the box wall and bottom to reach the thermocouple trees (Fig. 4a) with minimal influence on the regolith.(d) The CO 2 inlet diffuser reduced the likelihood of disturbing the regolith surface during the condensation procedure.(e) Coded photogrammetric targets were affixed at multiple heights and orientations.(f ) The box lid and actuator were used to isolate the regolith from thermal radiation during the cooling and condensation procedures, exposing the surface only for the sublimation procedure.(g) The regolith simulant for this run was JSC Mars-1.Parallel, long-slope marks were left from the protractor used during slope preparation.

Fig. 4 .
Fig. 4. Mars Chamber schematics.(a) Cut-away of copper cooling box with nominal dimensions.The two thermocouple trees were positioned along the centreline of the box.In all but the first six runs, an additional thermocouple (not shown) was positioned near the surface at the toe of the slope, c. 2 cm from the box wall.The 'Crest', 'Mid-slope' and 'Base' slope zones are also indicated.(b) Interior of the Mars Chamber, illustrating the relative positions of the cooling box, video cameras and heat lamp.The cooling lid remained closed throughout the cooling and condensation procedures, and then opened for sublimation.

Fig. 5 .
Fig. 5. Photogrammetric pipeline.The procedure for extracting data from the stereo video recordings made use of three principal software packages.Sony Vegas Pro was used to synchronize and extract coordinated image pairs for each DEM.Photoscan Pro was used to develop the DEMs based on the image pairs and physically measured target locations.Additionally, known hardware dimensions and measured slope angles were used to evaluate the accuracy of each DEM.The DEMs were then processed with ArcMap to derive volumetric and angular changes between DEMs.Finally, MatLab was used to compute slope angles and prepare the data for interpretation and presentation.

Fig. 6 .
Fig. 6.Resultant slope morphologies.(a) Oblique view of the crest slope zone of Run 20 (Table3) displaying scarp and ridge morphology at the end of the experimental run.The clearly defined areas of failure, delimited by fresh scarps faces at the slope crest, parallel ridges to either side and debris aprons, are typical of the observed discrete sediment flows.(b) Overhead view of a small scarp which formed c. 2 cm below the slope crest from (run 23 Table3).(c) Dust fans along back wall of box are the result of centimetre-scale jets of escaping CO 2 gas carrying entrained dust (from run 18 Table3).(d) Pitting around the mid-slope thermocouple tree (upper-most thermocouple has been exposed at the surface), caused by escaping CO 2 gas (from run 15 Table3).(e) Pitting at the base slope zone from run 24 (Table3).

Fig. 7 .
Fig. 7. Slope evolution for run 20 (JSC Mars-1 at AOR): (a) orthophotos showing visible appearance of the surface prior to frost condensation, at t = 0 (the start of sublimation), t = 80 min (the end of sublimation).At t = 0, surface frost is clearly visible as a white beard on the lower half of the slope, and as a smaller, dense accumulation at the back edge of the slope, adjacent to the box.(b) Difference rasters illustrating topographic changes between the pre-frost slope and the slope at the start and end of sublimation.Red represents erosion and blue represents deposition.Both the orthophotos and the difference rasters are oriented with the slope crest at the top of each image.(c) Evolution of the topographic long profiles for the slope pre-frost, and at the start and end of sublimation.The slope zones used to measure the slope angles are indicated in Figure 4a.

Fig. 8 .
Fig. 8. Slope evolution for run 14 (JSC Mars-1 at 17.5°): (a) orthophotos showing visible appearance of the surface prior to frost condensation, at t = 0 min (the start of sublimation), t = 90 min (the end of sublimation).At t = 0 min, a white beard of surface frost is clearly visible on the base slope zone, a thin (1-2 cm) concentration of frost encircling the balance of the slope along the box edges.(b) Difference rasters illustrating topographic changes between the pre-frost slope and the slope at the start and end of sublimation.The colour classifications representing elevation changes are the same as those in Figure 7.The singular classification illustrates that the scale of slope modifications was much smaller than that of run 20, in Figure 7. (c) Evolution of the topographic long profiles for the slope pre-frost, and at the start and end of sublimation.The slope zones used to measure the slope angles are indicated in Figure 4a.

Fig. 9 .
Fig. 9. Slope evolution for run 8 (fine sand at AOR): (a) orthophotos showing visible appearance of the surface prior to frost condensation, at t = 0 (the start of sublimation), t = 80 min (the end of sublimation).At t = 0, surface frost is clearly visible at the toe of the slope, and faintly visible on the lower half of the slope and along the edges of the box.(b) Difference rasters illustrating topographic changes between the pre-frost slope and the slope at the start and end of sublimation.Red represents erosion and blue represents deposition.Both the orthophotos and the difference rasters are oriented with the slope crest at the top of each image.(c) Evolution of the topographic long profiles for the slope pre-frost, and at the start and end of sublimation.The slope zones used to measure the slope angles are indicated in Figure 4a.
change.Example experimental temperature traces over the duration of sublimation for fine sand and JSC Mars-1, both initially near the angle of repose, are presented in Figure10a& b.The longer scale time (t*) output from our analytical model for JSC Mars-1 might be expected to be evidenced by a delayed start of regolith temperature increase, relative to fine sand.Neither of these predictions is borne out in the temperature data.In Figure10afor fine sand, initially, the deepest parts of the slope begin to warm with the cessation of active cooling and the upper regions continue to cool.Each trace, roughly in order of increasing depth, warms to the sublimation temperature (c.−120°C), where it remains constant during sublimation of the condensed frost before continuing to warm.The phase change dominates the temperature curves for hours, v. the <2 min suggested by t* in our analytical model (Table

Fig. 10 .
Fig. 10.Experimental and modelled temperature traces during sublimation.(a) Run 8 (fine sand at AOR).(b) Run 23 (JSC Mars-1 at AOR).(c) and (d) Simulations of the one-dimensional heat conduction equation with sublimation for fine sand and JSC Mars-1, respectively.Thermocouple locations are indicated in Figure 4a.Vertical positions are in centimetres above the bottom of the box.Solid lines correspond to mid-slope thermocouples; dashed rear.Blue traces correspond to the deepest locations; red traces are the closest to the slope surface and orange are in-between.

Table 2 .
Summary of regolith movement v. initial slope angle

Table 3 .
Slope angle, elevation and volume changes n.d., No data.

Table 4 .
Observed sediment transport types

Table 5 .
Analytical results for internal friction angle calculations