Geotechnical Properties Determination of Thickened Fluid Fine Tailings

Fluid fine tailings (FFT) comprising clayey-silt solids pose environmental and financial challenges. Currently, mining operators are depositing thickened FFT in deep-pits counting on self-weight consolidation to form stable ground. The motivation of this study is to model the long-term prospects of such deposits utilizing consolidation and direct shear strength measurements. The tests were conducted using scroll decanter centrifuge separated FFT sediment referred to as cake. Drained direct shear tests of the cake gave a linear Mohr–Coulomb failure envelope of 1.2 kPa cohesion intercept and 9° internal friction angle for normal stresses up to 1 MPa. Hydraulic conductivity of the cake was non-linear with normal stress decreasing to 1.7 × 10–11 m/s at 300 kPa. Consolidation results confirmed that the cake exhibits properties similar to those of active clay minerals. The cake compression index is governed by the same relationship as for active clays. The coefficient of consolidation for the cake was nearly constant and had a mean value of 2.25 × 10–3 m2/y, also similar to that of active clays. The void ratio—effective stress—hydraulic conductivity power law empirical relations were used to simulate settlement with a finite-strain model. Numerical results show that the top portion of an FFT deep-pit deposit remains in the liquid state while the lower portion whose maximum solids content converges to 74% is in plastic state. These mean that options that improve hydraulic conductivity and increase the shear strength of the thickened FFT should be integrated prior to final placement in order to create stable ground.


Introduction
In the extraction of minerals water is used in grinding, transporting, conditioning, cleaning, and separation of the valuable mineral from gangue. The process generates large amounts of mine waste tailings in the form of slurry. Considering the hot-water extraction of bitumen from surface-mined oil sands specifically, a 10% (mass basis) oil-containing ore generates a 55% solids discharge which is 82% sand and 17% fines (Chalaturnyk et al. 2002). Upon discharge, the entire coarse solids with some fraction of trapped fines settles out quickly while the rest of the fines flows into tailings ponds. There is in excess of 1.25 Bm 3 of fluid fine tailings (FFT) in containment ponds from this operation in Northern Alberta (AER 2019) alone. Other mineral resource extractions, such as phosphate and bauxite mining, equally produce large volumes of FFT that have similar characteristics as the FFT of oil sands. These FFT are discharged into containment ponds for water recovery and further processing for final disposal. Moreover, the ponds afford long term storage to stagger reclamation activities and time to apply emerging reclamation technologies. The FFT of tailings ponds eventually form a persistent suspension (referred to as mature fine tailings) whose solids mass relative to the wet mass plateaus to under 40% (w/w) in the case of oil sands extraction. These ponds pose major concerns as they require perpetual observation, construction, maintenance, risk of dam failure and increasing long-term liabilities. In response government departments place regulations to ensure safe operation and restoration to dry ground of tailings ponds at end of mine life (AER 2016). Reducing the water content-dewatering (thickening)-is a necessary step for restoration of the FFT ponds as a dry ground. All current commercial technologies for dewatering FFT start with the addition of flocculants and coagulants, which gather the suspended fine particles to form larger and denser aggregates, making them more amenable to subsequent processes.
The industry currently uses a few well-known tailings management methods for final tailings disposal after flocculation. Those methods are thin-lift drying, deep-pit deposition of flocculated FFT (also known as rim-ditch accelerated dewatering), and deposition of thickened tailings after separation using scroll decanter centrifugation (SDC) (BGC 2010). All these technologies in practice result in a thickened FFT whose dry solids content relative to the total wet mass converges to a maximum of about 60 solids% (w/ w). At a void ratio of about 1.6, such a deposit of fine solids lacks strength to form a stable ground. In this study we investigated whether or not such a deposit can in the long term become the anticipated reclamation-ready deposit by self-weight consolidation as a vast amount of it is being deposited in thickened form.
The purpose of the research described here was to determine (i) the geotechnical properties of thickened FFT to assess whether thickened FFT can form stable dry ground as a deposit, and (ii) conduct measurements to obtain constants for numerical modelling to predict settlement. Thickened FFT samples obtained using SDC separation, commonly referred to as cake, was used in this study. The cake by conventional categorization can be considered a thickened tailings sample (50-70% w/w). The initial solids content of the sample was 57% (w/w) with liquid and plastic limits of 61% and 16%, respectively and other properties are as given in Table 1. As a specimen the cake is no different from FFT sample that is dewatered to the same void ratio through other routes, such as thickening or consolidation in deep pits. The majority of studies on FFT consolidation were mainly conducted starting with high-void-ratio slurries (Been and Sills 1981;Pollock 1988;Jeeravipoolvarn et al. 2009). For this reason, there is a scarcity of data on the geotechnical properties of FFT with void ratios less than or equal to 2, particularly for oil sands FFT (McKenna et al. 2016;Yao 2016). Therefore, starting with cake sample extends the consolidation study to a stage beyond which this kind of tailings has so far been investigated.
Two types of studies were conducted to determine the geotechnical properties of the cake: consolidation and drained direct shear tests. Consolidation  (Cernica 1995;Skempton 1964). The two complementary studies provide information about the prospects of creating dry ground out of thickened oil sands FFT deposits. Although the origin of the cake sample is oil sands operations, the results are relevant to other mining tailings whose solids consist predominantly silts and clay fines.

Materials and Methods
The FFT used in this study was obtained from one of the major oil sands operators in northern Alberta. Samples were taken from a single batch of 10 m 3 of mature fine tailings (MFT) directly transported from the tailings pond for tailings management program. MFT refers to FFT that has undergone settlement for a few years and is a consistent gel-like sludge that shows no change in physical properties over time. The solids content of the MFT as received was 36% (w/w) (dry solids mass/total mass of sample).

Particle Size Determination
The special procedure developed for determination of FFT PSD begins with cleaning of the solids using Dean-Stark extraction (Dean and Stark 1920). The extraction removes residual hydrocarbon (bitumen) which interferes with the measurement. The dry clean mineral solids were then disaggregated in a mortar and pestle with care so as not to grind the particles. Prior to measurement the solids were further dispersed in a NaHCO 3 buffer at pH 9.6 aided by heating and 20 min of sonication in a bath. PSD were measured with a Sedigraph TM 5100 (MicroMeritics, Norcross, GA, USA). The analysis is based on the settling rate of particles in liquid as measured by the transmittance of a narrow beam of X-ray. Details of the method is found in ISO 13317-3. PSD of the FFT slurry is given in Fig. 1. PSD of the FFT from the oil sands operator from whose tailings pond the sample was obtained is consistent not only within the same batch but even between different batches as can be compared with that reported in Demoz (2015).

Methylene Blue (MB) Test
MB spot test for the determination of clay content including those of the American Society for Testing and Materials (ASTM) are based on the cationic MB adsorption and the first stable appearance (5 min) of a blue halo color on filter paper at titration equivalence point (ASTM 1984;Chiappone et al. 2004). To ensure complete interaction of MB with the FFT solids the standard methods were modified with additional cleaning and dispersion to handle oil sands samples (Currie et al. 2014;Young and Sethi 1981). MB titration of the solids after Dean-Stark extraction gave 12.73 mequiv per 100 g solids, which translates to a calculated clay content of 91 mass%. In contrast, MB titration of the as-received FFT slurry titration gave a clay content of 58 mass% which is more representative of the actual tailings. The difference in clay content between the two methods is known to arise from the effect of sample preparation on surface area. Thickened samples were prepared by solid-liquid separation of the FFT using SDC. The FFT was first diluted to 23% solids (mass of solids/total sample mass) in a holding tank and flocculated using inline static mixers before the SDC inlet. The SDC was operated so that [ 99% of the FFT solids was captured as cake. The cake sample was kept in sealed 20 L pails. The full details of the FFT continuous SDC separation process is found in Demoz (2018).

Atterberg Limit Measurements
The plastic limit was measured using the soil-thread method in accordance with ASTM Standard D4318-10 (ASTM 2007). The liquid limit, W L , was measured following the ASTM cone penetration method ASTM D3441-16 (ASTM 2016).

Direct Shear Strength Measurement
The shear strength of specimens was measured using the direct shear test method ASTM D3080-98 (ASTM 1998). A 100 mm circular shearbox was filled with the saturated pasty cake specimen (57.5% w/w solids) directly using a spatula. The samples solids contents varied by \ 0.5% w/w from sample to sample. The pasty cake was easily worked by hand and also using circular aluminum pusher to completely fill the shearbox. Drainage was facilitated by placing the cake between two porous stones: one underneath the stainless steel vertical-force loading pad and the other above the bottom stainless steel retaining grid of the shearbox. To maintain the saturation of the sample, deionized water was added around the edges of the stainless steel loading pad as needed and there was process water in the shearbox carriage all the time. A high-vacuum grease (Dow Corning Corp. Midland, MI, USA) was lightly applied around the edges of the shearbox to reduce edge-to-edge friction and to seal any gaps around the edges of the shearbox.
Drained direct shear test was carried in two phases. The first phase was consolidation under the test normal stress. Multiple smaller load steps had to be used to reach the test normal stress because otherwise the FFT would exude around the edges of the shearbox. The consolidated specimen was thereafter sheared for 10 mm horizontal displacement at a slow rate of 0.025 mm/min while the test normal stress was maintained constant. This shear rate was assessed to be sufficient for excess pore pressure that could develop during shearing to drain. The direct shear test instrument used was a Humboldt E5706 instrument (Humboldt Mfg. Co., Elgin, IL, USA). This equipment is instrumented for automatic logging of vertical and horizontal displacement using linear variable displacement transducers (LVDTs). Both horizontal and vertical loads were measured by 10 kN maximum load cells. Normal load on sample was applied using pneumatic method. A geared stepper motor applied horizontal deformation.

Compressibility Measurement
The compressibility behaviour of the cake was measured using the Humboldt E5706 instrument, too, but operated as a consolidometer (oedometer). Cake sample preparation and normal stress application for the compression tests were the same as for the direct shear tests. The normal stress was increased sequentially following ASTM D2435-11 (ASTM 2011). Monitoring of the outputs of the Humboldt E5706 LVDTs and control of the instrument were all done by a computer using NEXT TM application software from Humboldt Mfg. Co.

Hydraulic Conductivity Measurement
Hydraulic conductivity of the cake was measured in a large-strain consolidation (LSC) cell. Small constant heads provided low-pressure gradients for hydraulic conductivity measurements. The particular cell design and its use for permeability measurements are well referenced in many publications of the University of Alberta, Department of Civil and Environmental Engineering (e.g., Suthaker and Scott 1994;Scott et al. 2008). The constant head permeability measurement method was selected because it allows for extremely small head drop tests.

Direct shear strength
Most reports of the shear strength of flocculated fine tailings are from undrained strengths obtained using field vane or cone penetration testing (Elias and Beier 2017;Fisseha 2020;Sobkowicz 2013;Wilson et al. 2018). Any drained shear strength measurement via direct shear test of oil sands tailings available in the literature thus far has been on samples that have a sizeable fraction of coarse particles as well, which is not representative of the large inventory of FFT being held in tailings ponds (Fisseha 2020;Kouakou 2014).
Direct shear testing has long been used to estimate drained shear strength parameters of soil for use in the analysis of slope-stability, retaining-wall, and bearingcapacity problems (Cernica 1995;Skempton 1964). In direct shear testing the sample is sheared in a prescribed plane, and although this may be an issue for soils that have developed structure, it is less relevant for uniform, short-lived specimens like FFT cake. The grain sizes of the cake (see Fig. 1) are small, and shear can develop within the cake without any influence from the edges of the shearbox (Palmeira and Milligan 1989).
Drained shear strength of clays and silts depends on factors such as the mineralogical composition of the grains, grain size, grain shape, surface texture, grading, moisture content, normal stress and shear rate. The PSD of the FFT and the Atterberg limits of the cake are shown in Fig. 1 and Table 1, respectively. The PSD shows that the FFT solid particles can be considered 100% fines. The D 10 value for the cake was below the detection limit of the measurement method used, which is 0.3 lm. The coefficients of curvature and uniformity for the FFT have no relevance as the fines content is way more than 10% (Bowles 1984).
The mean particle size of the FFT is 1.37 lm (Fig. 1). As a clayey material with substantive affinity for water, the long times needed for full consolidation were maintained before the shear strengths were measured. By contrast, just a brief consolidation time is adequate for clean coarse sand, which drains readily, as noted by Kovacs (1997). A soil whose volume remains constant under a constant state of stress can be said to be fully consolidated. When the displacement stops to change, it indicates that the applied stress is born by the grains in the sample and the pore pressure has fully dissipated. For saturated soil the volume change is due to water escape from the voids. In the absence of pore pressure measurement options, the cessation of vertical displacement of specimens was used as the condition of full consolidation in this study. Total pore pressure dissipation is virtually reached when the specimen height ceases to show any change with time. All specimens were consolidated for comparable time periods as that presented in Fig. 4 before starting the direct shear measurements.
The realistic and sustained state of thickened FFT discharge sites is that the deposits remain saturated and open to draining. The direct shear tests were therefore conducted on saturated and drained specimens so as to correspond to real conditions. Figure 2a presents drained direct shear stress versus horizontal displacement at different normal stresses up to that higher than the self-weight of a deep-pit deposit can exert. The shear stress-displacement plot does not display peaks with six of the eight test normal stresses and the shear stress increased with shear displacement without approaching a constant. This is to be expected as peak stress and dilatancy tend to be observed with overconsolidated or dense sandy samples both of which a cake is not. Figure 2b illustrates vertical contraction with horizontal displacement of cake. There were slight contractions which ranged from 1.8% to 4.4% relative to horizontal displacement. In Fig. 2 Curves of shear strength (a) and vertical displacement (b) versus horizontal displacement for cake in a direct shear test these kinds of large displacement tests up to 2% of contraction are artifacts of measurements as reported by Bareither et al. (2008). The very small contractions of the cake indicate that the shear zone remained localized in the plane near or at the interface between the two shearbox halves.
Shear failure is frequently associated with the point in the test where the maximum shear stress is reached. The drained direct shear test results show that a shear displacement of 1.2 mm was sufficient to mobilize the peak shear stress; beyond that point the shear stress kept the same plateauing trend. There was also a lag of about 0.5 mm at the start between shear displacement and shear stress. The lag was not specific to the specimen but was equally observed in tests conducted without sample done to determine the shearbox edgeto-edge friction indicating that the lag was due the presence of lateral slack in the shearbox carriage system. Where there is no clearly defined failure point, a select strain can be used as the failure point. In this study the shear stress at 3 mm displacement is reported as the failure shear stress. The uncertainty in shear strength due to this interpretation is minimal because the increase of the measured shear stress with displacement is small. Failure in a soil occurs when the stresses on any plane are such that where s 0 is the effective shear stress at failure; c 0 is the effective (drained) cohesion intercept; / 0 is the effective internal friction angle; and r 0 is the effective normal stress at failure. The line plotted for shear stress dependence on effective stress forms the Mohr-Coulomb failure envelope, as shown in Fig. 3. The low effective stress plot intercept is very close to the shear strength of fresh cake that is commonly observed using vane shear tools (Sobkowicz 2013). Some portion of the shear strength is expected to be generated by the friction between the specimen and the edges of the shearbox. To correct for this, shear tests were conducted without any sample but with a thin film of the high-vacuum grease and cake mix applied on the edges of the shearbox. A linear dependence of shear stress and normal stress was observed for the edge effect, and this correlation was used to correct the edge friction effect. The experimental drained shear strength of the cake after correction of the edge friction effect is as presented in Fig. 3.
Generally, the shear strength of freshly produced cakes is around 1.2 kPa at most (Sobkowicz 2013). A constraint of a 1.2 kPa intercept was therefore introduced in deriving the Mohr-Coulomb envelope (Merrien-Soukatchoff and Omraci 2004). The generalized Mohr-Coulomb failure model after correcting for the edge-to-edge friction effect over the test stress range yielded an internal friction angle of 9°. The Mohr-Coulomb failure envelope for the consolidated, drained cake (thickened FFT) tested in this study followed the empirical relationship given below: The intercept, s 0 and r 0 of Eq. 2 are all in kPa. Both the effective cohesion strength and the effective internal friction angle of cake are small. It has been known that the failure envelope at low confinements can differ from that obtained over relatively large confinements (Nafisi et al. 2020). Even nonlinear failure envelopes have been proposed to fit data that could not be expressed with linear failure envelope expressions (Sharma et al. 2011). A bilinear failure envelope can be observed when the confinement stress is limited to two ranges; one for the low confinement range up to100 kPa and an upper confinement range from 200 kPa. A higher internal friction angle of 17°was obtained in the lower confinement range as shown in Fig. 3 inset. The bilinear low confinement fit is relevant for comparison of low confinement stress Inset shows low confinement stress versus shear stress for a bilinear fit. The error bars are standard deviations from measurements of up to three freshly prepared samples measurements but not for practical considerations as high confinements would be required for stable ground formation. From MB titration and PSD, we determined that the cake is largely made of clay minerals and silt. Its lack of intercalating and cementing minerals explains the small cohesive shear strength of cakes. The cake can be considered a remoulded sample, where no extended structures exist between the fine particles. The strength of a remoulded soil sample is lower than that of a soil in the undisturbed state. The more finer grained particles there are in a soil, the smaller the inclination of the failure envelope, as is also confirmed in this study (Mitchell and Soga 2005). Hu et al. (2017) have reported the strength indices of iron and copper fine tailings whose D 50% were 30 and 60 lm, respectively. The strength indices of both the iron and copper fine tailings were more than three times that of the oil sands fine indices and all three in reverse order with respect to their particle sizes.
The FFT index properties show that the solids have high activity and that their plasticity index (PI) is very high (PI [ 40) (Das 1998;Rankka et al. 2004). The index properties indicate that the FFT are composed largely of swelling clay particles. Higher clay contents in a soil cause higher plasticity, greater shrinkage and swell potential, lower hydraulic conductivity, higher compressibility, and lower internal angle of friction, as were observed for the cake. Each of these properties is a drawback for the creation of a stable ground from FFT (Mitchell and Soga 2005).
The deposit needs to be firm to stiff in consistency (i.e., at least 50 kPa shear strength) for it to support small earth-working equipment (McKenna et al. 2016). An effective normal stress of 308 kPa is needed to get a 50 kPa drained shear strength according to Eq. (2). Including the difficulty of laying overburden on high void ratio deposit, the high normal stress threshold indicates that the upper part of a thickened FFT deposit does not have the strength to support earth-working equipment.

Cake compressibility
The shear strength laboratory analysis above does not give an understanding of the period to attain full consolidation. The vertical deformation of the cake sample from the consolidometer tests is presented in Fig. 4. The time to reach full consolidation is shorter the higher the applied normal stress. For instance, the normalized vertical displacement changes after 500 min of consolidation were 60, 89, and 91% of their final settlement at normal stresses of 6.25, 100, and 750 kPa, respectively in this test. The consolidation times at the three stresses in increasing order, for instance, were 73, 41, and 32 h, respectively. Therefore, the consolidation periods were longer at lower loads so as to ensure completeness of primary consolidation.
The void ratio, e-log r 0 relationship for ultra-soft soil is nonlinear with a large compression index and no gain in strength (Jeeravipoolvarn et al. 2009). The elog r 0 relationship for the cake was linear, similar to that for soils that naturally have a far smaller compressibility by comparison. The compression index, C c , of the cake for void ratio from 1.6 to 0.57 was constant and equal to 0.45. This is practically the same as that from the well-known C c relationship to liquid limit, C c = 0.009*(W L -10), for normally consolidated clays, which for the cake yields C c-= 0.46 (Das 1998;Terzaghi et al 1996). The cake C c is in the middle of the range for silty clays, which is 0.36-0.54, and double that of sandy clays (Whitlow 2001).
A compressibility curve in general may consist of three distinctive regions (Terzaghi et al. 1996;Mitchell and Soga 2005;Khan and Azam 2016). It may start with a horizontal initial portion that pertains to apparent pre-consolidation compression, followed by a linear portion for primary consolidation, and finally by a flatter curve corresponding to secondary compression that is associated with creep settlement. The cake compression curve exhibited neither the initial horizontal portion nor the post primary compression curve in the studied effective stress range. The absence of an initial flat curve is to be expected, as this was a freshly produced material with no consolidation history. In general, fine tails are underconsolidated and therefore should follow a virgin compression line (Suthaker and Scott 1994). The secondary compression, stage was not observed at the highest applied axial loads, either. Figure 4 illustrates the settlement of the cake with time at increasing normal loads. In order to express its compressibility, the void ratio at the end of each consolidation stage was calculated using Eq. 3. The vertical displacement data and the void ratio obtained from moisture content measurement after the sample was removed from the consolidometer were used to determine the void ratio and the corresponding effective stress dependence as shown in Fig. 5.
where e o is the void ratio and H o is the corresponding sample height, De and DH represent changes of these parameters at completion of each consolidation step. The void ratio calculated using Eq. 3 dependence with effective stress is shown in Fig. 5. The best fit curve is also included as part of finding the constitutive equation for effective stress-void ratio relationship. This equation is used to calculate the settlement of cake using the finite strain consolidation model later on. The relationship between void ratio and stress for highly compressible material and soils is often expressed by a power function (Carrier et al. 1983;Pollock 1988). The e-r 0 power fit relation below has a best least squares fit correlation factor of 0.98.
where r 0 is effective stress in Pa.

Coefficient of consolidation and hydraulic conductivity
The coefficient of volume compressibility, m v , and the coefficient of consolidation, C v , dependence on e and r 0 are plotted in Fig. 6. C v is regularly obtained using Terzaghi's infinitesimal-strain consolidation theory (Das 1998;Whitlow 2001). The theory is inapplicable for thickened FFT (cake) because three of the conditions of the infinitesimal theory: specifically very small settlement, linear compressibility and liner permeability were all not observed with this sample.
There are other methods besides Terzaghi's consolidation-theory for determining C v (Islam et al 2020). In this study C v were determined using Eq. 6. This approach was preferred because it fully relies on the measured parameters of the sample itself. where the coefficient of compressibility, a v , = de/dr 0 ; k c is the hydraulic conductivity constant (m/s); and c w is the unit weight of water (9.81 kN/m 3 ). All parameters for eq. [6] were obtained empirically from consolidometer and LSC test results. Cake hydraulic conductivities were measured at a fewer number of effective stresses and up to a maximum effective stress of 300 kPa. In order to increase the reliability of C v and m v of Fig. 6 only empirically obtained values were utilized in their calculations. The mean of the cake C v is 2.2 9 10 -3 ± 8.5 9 10 -4 m 2 /y, which is in the range reported for active clays (Lambe and Whitman 1979;Mitchell and Soga 2005). As a comparison, C v of iron, copper and coal fine tailings all are more than two orders of magnitude higher (Hu et al 2017;Islam et al 2020). The value for m v decreased with increasing r 0 and increased with increasing e, as expected. For cake, as for tails with higher water content, an increase in effective stress causes the hydraulic conductivity to decrease due to a significant decrease in void ratio (Pollock 1988;Khan and Azam 2016). The value of m v varied from 0.03 to 0.0002 m 2 /kN, similar to the range measured for underwater dredged fine slurries (Imai et al. 1984).
The cake hydraulic conductivity, k c , was measured under load using the constant head method (Suthaker and Scott 1994;Scott et al. 2008). Figure 7 shows the k c of cake at different void ratios. Values of k c between 10 -7 m/s -1 and 10 -9 m/s are understood to be poor, whereas deposits with values below that are considered practically impervious (Terzaghi et al. 1996).
The hydraulic conductivity of cake decreased by a factor of about 60 for a decrease in void ratio of only 1. The k ce relationship for the cake has the typical power function: where A is a constant with units of velocity (m/s) and B is dimensionless constant. The values of A and B for the cake are 1.23*10 -10 m/s and 3.40, respectively. These are of similar magnitude to thickened FFT reported by Owolagba and Azzam (2015) as well as for phosphatic clays (Abu-Hejleh et al. 1996;Khan and Azam 2016). Finally having determined the constitutive relations between e and k c and e and r 0 of the cake finite-strain consolidation theory was applied to simulate the settlement of its deposit over long times.

1D modelling
The question that this study examines is evaluate whether thickened FFT and cake could form timely geotechnically stable deposits. FSConsol, a finitestrain consolidation modelling software (GWP Geo Software Inc., Edmonton, AB) that is widely used in industry, was used for the analyses in this study (Tito 2015). The modelling software is based on the finitestrain settlement solutions given by Gibson et al. (1967Gibson et al. ( , 1981. The application provides the sought outputs of solids content profile with depth and the deep-pit deposit settlement as a function of time (Tito 2015). In addition to the constitutive equations (Eqs. 4 and 7), the inputs of the program include the initial solids content, pit capacity, and rate of filling. To complete the inputs to the model, a 30 m deep pit is set to fill with cake in 8 years and specified to consolidate in two ways; by upward drainage only and drainage through both upward and downward flows. The deposit filling conditions and the constitutive relations used in the simulation are given in  Figure 8 presents the solid content profile at the same four periods cited above. The increase in solids content begins from the bottom, which is to be expected as it experiences the higher compression.
The modelling also shows that the dewatering converges towards a maximum solids content of 74%. The model results show that even after hundreds of years of consolidation the top third of the deposit is in the liquid state, absent the impact of environmental effects such as freeze-thaw and desiccation. As well, the portion of the deposit that consolidates to the maximum solids content is in the plastic state closer to the liquid state as its liquidity index of 0.42 indicates. Realistically, the self-weight consolidated cake or thickened (dewatered) FFT deposit will remain saturated and of such a high moisture content that it cannot by itself be relied on as a route toward reclamation in the form of a dry ground. Similarly, thickening of FFT using filtration gives a deposit no different in properties than that of SDC separation. Backfilling in all its variations is an operation that cannot be exercised with surface mines like the oil sands (Yilmaz et al. 2012). These suggest that co-disposal of the thickened tailings with mine coarse tailings is the only realistic alternative if FFT ponds are to be reclaimed as dryland. The advantages of the co-disposal processing are likely to be realized if the coarse tailings is added into thickened (plastic consistency) FFT rather than the high void ratio slurry (viscous) FFT (Wilson et al. 2018).

Summary and conclusion
The large and still increasing inventories of FFT from global mining has serious engineering, environmental and economic costs. The geotechnical properties of FFT from the oil sands extraction were determined to evaluate whether such fine tails can develop into dryland through in-pit consolidation. Direct shear strength for confinement up to 1 MPA gave a linear Mohr-Coulomb failure envelope with very small cohesive intercept and small internal friction angle indicating that even thickened deposits of FFT will remain metastable and liquifiable. Thickened FFT hydraulic conductivity was very low and similar to that of active clays. The e-r 0 -k c empirical power law curves constitutive constants were used as model inputs to predict tailings settlement using a finite-strain 1D model software. The model output indicated that consolidation by selfweight would take an extremely long time to create reclamation-ready deposits. The implication of this study is that the integration of options in the operation which increase shear strength and hydraulic conductivity are a necessity for stable dryland formation. Codisposal with available coarse tailings is the realistic route to improve both the shear strength and hydraulic conductivity. Work is underway to determine analytical relations between hydraulic conductivity and sand-to-fine ratio for co-disposal schemes.
Acknowledgements Funding for this work was provided by the Natural Resources Canada Program of Energy Research and Development. Sincere thanks to colleagues in the Environmental Impacts team, in particular Dr. K. Kasperski for team leadership, and J. Elias and C. McMullen for help with the LSC measurements.
Funding For this work was provided by the Natural Resources Canada Program for Energy Research and Development.

Declarations
Conflict of interest The author has no relevant financial or non-financial interests to disclose.
Code availability FSConsol (GWP Geo Software Inc. Edmonton, AB) is a commercially available software.
Data availability All data used is included in the manuscript and can be provided upon request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.