A High-Fidelity Modelling Method for Mine Haul Truck Dumping Process

: Dumping is one of the main unit operations of mining. Notwithstanding a long history of using large rear dump trucks in mining, little knowledge exists on the cascading behavior of the run-of-mine material during and after dumping. In order to better investigate this behavior, a method for generating high ﬁdelity models (HFMs) of dump proﬁles was devised and investigated. This method involved using unmanned aerial vehicles with mounted cameras to generate photogrammetric models of dumps. Twenty-eight dump proﬁles were created from twenty-three drone ﬂights. Their characteristics were presented and summarized. Four types of dump proﬁles were observed to exist. Factors that inﬂuence the determination of these proﬁles include the location of the truck relative to the dump crest, the movement of the underlying dump material during the dumping process and the differences in the dump proﬁle prior to dumping. The HFMs created in this study could possibly be used for calibrating computer simulations of dumps to better match reality.


Introduction
Mine-to-mill optimization is a longstanding goal of the mining industry [1][2][3][4]. This approach focuses on optimizing the entire process around run-of-mine (ROM) material characteristics, rather than optimizing the unit operations to material characteristics individually. The primary motivation for this approach is ensuring that the priority for optimization is given to the most demanding and costly process (grinding) [5]. Despite its success, one difficulty in holistic mine-to-mill optimization stems from a lack of understanding around material behavior during and between unit operations [4].
Unit operations are defined within the context of mining as the basic steps used to produce mineral value from a deposit. They generally fall into either the category of rock fragmentation or materials handling [6]. Materials handling, for surface truck and shovel mines, comprises three steps, known as loading, hauling and dumping. Dumping commonly consists of the haul truck spotting itself into position and dropping ROM material from the back. This material forms a small heap if it is dumped on a flat surface, or cascades along the edge of a dump face if dumped over a developing dump, stockpile or dump/heap leach [7]. While there are minor differences in each of these earthworks, for the simplification of this article, they are referred to as rock piles [8].
Notwithstanding a long history of using large rear dump trucks in mining, little knowledge exists on the cascading behavior of ROM material during and after dumping. At least two schools of thought comprise the knowledge that does exist. The first school of thought studies ROM cascading behavior from an external perspective. The second seeks to understand the smorgasbord of characteristic properties that govern this behavior internally. Since this paper focuses primarily on the external viewpoint, only a sample of the literature pertaining to the internal philosophy is presented in order to call attention to some of its challenges. McLemore et al. [8] provide a robust review that includes a trove of references from each school of thought, should readers wish to dive deeper.
Work related to factors that influence ROM material behavior from an internal philosophy are numerous, and a complete review of their literature is not within the scope of this paper. A plethora of material property variables relate to these factors, which include shear strength, bulk density, particle size and shape distribution, cohesive properties, friction angle, moisture content, etc. [8,9]. Several issues confound the ability to fully isolate and understand these variables. First, the size, shape and variability of the ROM material are greatly dependent on blasting, which is not yet a fully understood process [10,11]. Second, the successful testing of material property variables depends heavily on sampling and statistical estimation [12]. Furthermore, dump faces are hazardous, and observation of the dumping process is challenging [13]. Additionally, dozers and other equipment handle the material subsequent to dumping, which compound with the issues mentioned previously [7]. Moreover, some complex engineering issues, such as reclamation and slope stability, are commonly considered to be linked to factors involving these material property variables in ways not fully understood, which increases the debate related to them [13][14][15]. Mines may also be under the assumption that understanding these factors must originate from a first-principles approach, rather than a data-driven statistical and mathematical approach [16].
Traditional viewpoints and operational strategies within the mining industry generally hinder innovation [17,18]. This might be particularly true for each of the three kinds of rock piles. Stockpiles have been seen as only useful to mines as a buffer against production variability [19][20][21]. Dump/heap leach design optimization involves many meticulous considerations that take precedence over understanding the minutia of the dumping process, which may be considered to be optimized under simple guidelines [22]. Waste dumps have been reserved for material below the economic cut-off grade and, as a result, have historically been given little attention beyond safety and risk management [23]. While these traditional viewpoints would need to change in order for rock pile innovation to be successful, they need only be challenged in order for research in these areas to be justified.
Literature and news articles provide both cases and causes for some changing perspectives that motivate and support the endeavor of this paper. First, even though the primary function of stockpiles remains the same [24], COVID-19 and global supply chain disruptions have placed increased demand on their use [25][26][27][28][29][30][31][32][33][34]. Increased knowledge of stockpile assets as the result of understanding the dumping process might afford mining companies a competitive edge in a post-pandemic world [35][36][37][38]. Second, while much of the seminal work in heap leach modelling involved small laboratory column studies operating at the particle scale [39][40][41], there has been a recent shift towards modelling bulk scale phenomena (inter-/intra-particle diffusion, liquid holdup and hysteresis, gas flow, etc. [42,43]) as well as understanding the factors more closely correlated with the dumping process (stratification, segregation, breakage induced by ore stacking, etc. [44][45][46]). It is true that a better understanding of the dumping process will likely not lead to changes in the construction of heap leaches. However, it may yield an understanding of the gestalts about the bulk phenomena in existing heap leaches, and aid in the development of models that increase their profitability [47]. Third, what was once considered to be below cut-off grade may later become mineable ore [48][49][50][51][52]. Therefore, a deeper understanding of the dumping process as it relates to waste dumps may prove valuable for operations, where old waste dumps become economically viable [53][54][55][56]. Finally, the perfunctory amount of documented work, aimed at capitalizing on the opportunities hypothesized previously, may have less to do with the merit of such hypotheses, and more to do with the mining sector's lag in innovation [17,18,57].
To improve the understanding of end-dumping and rock pile construction, we present a method for creating a high-fidelity model (HFM) of the dumping process through a digital transformation approach. Digital transformation [57,58], employs an external philosophy for understanding ROM behavior, and is a process whereby real-world assets and processes are digitally transformed in ways that add value for decision makers. Innovation [59] and emerging technologies [60] make digital transformation possible. Unmanned aerial vehicles (UAVs), also known as drones, are an emerging technology [61] that allows for the digital transformation of the dumping process through data collection. Rock pile faces span large areas, and the nature of the material requires multiple vantage points for sensors or traditional surveys to be effective. By using photogrammetry, UAVs are capable of capturing multiple angles and covering large difficult areas, such as a dump face [47,62]. Thus, an investigation into whether photogrammetry can create the HFM of an individual truck end dump is of interest to the issue of better understanding the dumping process and rock pile construction.
HFMs act as reference systems to reality [63]. These models are commonly used in engineering to calibrate simulation models [64] where there is a need to rapidly prototype many different design permutations [65], or where measuring the real data being modelled is not practical [66]. HFMs are modular elements, and HFM integration is a modular framework that allows each aspect of the full model to be worked on independently [67]. This modularity means that simulations to match HFMs can be worked on independently from the work used to improve how HFMs match with reality.
Digital transformation techniques are fairly new, and little work has been conducted to digitally transform the process of rock pile construction via haul truck end-dumping. Zahl et al. [9] accurately assert that the formation and shape of mine rock piles are based mainly on topography. While this work is foundational to the engineering of rock pile construction, it offers little for the purposes of rock pile digital transformation. Mclemore et al. [8] provide an extensive review on the construction, the factors influencing the shear strength of soil, characterization, the effects of weathering and the stability of rock piles throughout the world. Their review is more informative than that of Zahl et al.; however, it likewise does not consider a digital approach. Zhao [68] developed a real-time 3D modelling and mapping technique for the stockpiles formed by stacking/reclaiming machines for iron ore. Zhao's work is pioneering in the area of the digital transformation of small, intermediary stockpile construction, but does not consider the large rock fills made by haul truck enddumping. The authors of this present paper previously illustrated a method for modelling and mapping large heap-filled stockpiles using fleet management system (FMS) data [7]. However, while these data were amicable to modelling, no external data were available for the validation of the model. Zhang and Liu [47] employed UAV aerial photography and performed image analysis to investigate particle size distribution along the face of a dump leach. However, while they demonstrated the ability of UAV aerial photography to capture relevant data, they did not capture the volumes of individual dump profiles that could be used as a baseline to digitally transform the end-dumping rock pile construction process. Servin et al. [4] present a digital transformation technique for holistic mine-to-mill distributed particle simulation where the gaps between unit operations are simulated using data from control systems and sensors. While other frameworks for similar integration have been proposed in the past, the focus Servin et al. place on a unit operation-centered framework is consistent with the aims of this paper and the future direction of research in this area.
The method presented in this work demonstrates how to create 3D HFMs of ROM end-dumping from photogrammetry. The method presented is very similar to that of Zhang and Liu [47]. The resulting HFMs are akin to the model developed by Zhao [68], only they involve end-dumped rock piles as opposed to stacker-made intermediary stockpiles. These HFMs could be potentially suitable for the technique presented by Servin et al. [4] for holistic mine-to-mill optimization, which includes large end-dumped rock piles. They may also be used as a modular test bed for the calibration of future simulations of mine haul truck dumping activity. The end goal of the HFMs this method creates is the parameterization of dumping as a unit operation process. Essentially, while the foundational work of Zahl et al. [9] explains that topography is the main factor influencing rock pile shape, the method presented in this paper may be used to increase our understanding of the expected variability and parameters in such topography common to the end-dumping and rock pile construction process, thereby making these shapes more predictable and amicable to modelling.

Material and Haulage Equipement
The material studied was low-grade run-of-mine (ROM) gold ore from a surface mining operation near Perth, WA. No additional characterization of the material is required for presenting the concept of this method, although future simulation and modelling work will require more detailed characterization in order to ensure that models are properly calibrated.
The ore was transported and dumped using CAT 793F haul trucks with a payload of approximately 231 metric tons, a struck capacity of 112.6-151 m 3 and an inside bed width of 7334 mm.

UAV Flights and Photogrammetry Methodology
Field work and UAV flights were conducted in coordination with the mine's survey crew and through the use of their equipment. A total of 23 flights were recorded over the course of five day-shifts. These flights were performed using a DJI Matrice M200 drone with a Zenmuse X5S camera (DJI-Innovations Company Limited, Shenzhen, China). Table 1 shows details about the camera and flight settings used in these flights. The methodology for the flights was as follows: 1. An initial flight was performed to create and update the local orthomosaic map (approximately 6 min); 2.
Ground control points (GCPs) were marked and surveyed (3 to 7 GCPs were used for each flight); 3.
A second flight was performed with an adjusted flight path to include ground control points (approximately 7 min); 4.
A Quick Map was created from the second adjusted flight path so that the area of interest (dump area) could be readied, and the working flight path programmed; 5.
Once steps 1 through 4 were completed, flights were repeated using the site scan Quick Fly software function in the same area along the same flight path to capture the dump face during different time intervals and thereby capture before and after photos of the dumping process during regular operation (60 m flight height, 6 min flight times).
Once the flights were completed, pictures were uploaded to the 3DR cloud processing system, and photogrammetric point clouds were generated from the flight photos using Pix4D software (mapper 4.1 version, Pix4D SA, Lausanne, Switzerland). Absolute and relative accuracy were calculated for each GCP in the X, Y and Z coordinate directions. These calculations were performed automatically as part of the 3DR cloud processing. Figure 1 shows a top view (XY plane) of an example flight plan with GCP locations and the corresponding image overlap map generated from the flight. All flight plans were similar to Figure 1, with the only changes being the location of different GCPs. The thick black curved line in Figure 1 represents the upper crest of the dumping area, and it is a close approximation to the actual dump crest. As demonstrated by Figure 1, the dump crest area for each flight was photographed with 5 or more images of overlap. Figure 1 shows a top view (XY plane) of an example flight plan with GCP locations and the corresponding image overlap map generated from the flight. All flight plans were similar to Figure 1, with the only changes being the location of different GCPs. The thick black curved line in Figure 1 represents the upper crest of the dumping area, and it is a close approximation to the actual dump crest. As demonstrated by Figure 1, the dump crest area for each flight was photographed with 5 or more images of overlap.

Point Cloud Analysis
Maptek™ PointStudio (2021.1 version, Maptek/KRJA Systems Inc, Golden, Colorado, USA) was used to model and analyze the resulting point cloud data. These data were imported as individual files and converted into triangulation solids. Then, they were analyzed for overlapping areas where dumps had occurred. Polygons were created around the dump areas, and the two triangulations were processed into extracted solids. These extracted solids are considered to be 3D volumetric HFMs of individual haul truck dump profiles within the simulation and modelling context described in the introduction. The reason these software tools, as well as the drone equipment, were used is because they were what was available and provided by the mine. The results section details the HFMs created by this study. All HFM files are available at the link, https://zenodo.org/record/ 5789951#.YgQiqJbMKUl.

Classification Method
Little taxonomic terminology yet exists for the information presented in this article. Classification involved qualitative analysis by the authors on the HFMs, with an emphasis on attributes hypothesized to be of interest to characterizing the final position and geometry of the dump profile. From the qualitative analysis, the categories and characteristics of the dumps were created, and each dump was assigned to the matching category. Data on the resulting classification are presented along with the statistical information of each category. Tables 2 and 3 show the root mean square (RMS) error values and the ground sampling distance (GSD) of the models computed from the drone flight data, respectively. Errors between the modelled coordinates of a GCP and its known survey coordinates are frequently used to represent the accuracy of a photogrammetry model [69]. RMS errors represent the quadratic means of these errors in the X, Y and Z coordinate directions. Both absolute and relative RMS error values are given. Absolute error represents the error in the coordinate location of the GCPs to their actual location on earth. Relative error represents the error in coordinate location of the GCPs as they relate to each other. Relative error is of more interest to this study, since the integrity of the HFMs is reliant on achieving a low relative error in the photogrammetry models.

UAV Flights and Field Work Results
Low ground sampling distance (GSD) values are required for accurate photogrammetry [69]. GSD is the physical distance represented between the centers of two adjacent pixels, and it can be estimated in advance of a flight based on the camera and flight details. With the drone flying at a 60 m height, a GSD of 1.64 cm/pixel is considered to be the lowest possible GSD value based on the camera specifications. Slight deviations due to perspective, the vibrations of the camera, blur, depth of field and other factors lead to variations in GSD at every point of a photogrammetric model. The GSD values calculated for the photogrammetry models of this study can be found in Table 3.
In Tables 2 and 3, the accuracy of the photogrammetry models is given. These accuracy values are important to consider, since they determine the granularity of the HFM models. The accuracy error of a given model will be at least the GSD value in each coordinate direction. Considering that the highest average GSD value of all flights was 2.52 cm/pixel, an accuracy error of at least 16 cm 3 (2.52 cm × 2.52 cm × 2.52 cm) to the volume of the HFM is expected to exist. Another way to confirm these accuracy errors is to multiply the RMS errors (Table 2) for each coordinate direction. The multiplication of the maximum relative errors for each coordinate direction, as shown in Table 2, gives 45.48 cm 3 , which is roughly 2.8 times the accuracy error obtained by cubing the GSD. This discrepancy is consistent with other photogrammetry models that have been correctly scaled and reconstructed, which typically contain accuracy errors between one to three times that of the GSD value. In summary, the accuracy of the volumes generated from the photogrammetry models used in this study can generally be considered to be accurate to within 50 cm 3 .

Point Cloud Analyses and Results
Of the 23 flights flown, 14 solids were extracted, from which 29 dump activities were identified and 28 were considered useable HFMs of the dumping process. Of these 28 HFMs, four types of dump profiles were determined to exist. The four dump types are named after their shape, as follows: Oval-type dump profiles are the most commonly occurring type of dump profile. These dump profiles are characterized by their oval shape when viewed from a vantage point normal to the dump face. These dumps have narrow ends at the crest and toe of the dump face, and a maximum width midway through the dump face. Comet dump profiles are characterized by a large volume near the base of the dump and a narrow trail of material extending upwards along the dump face. Rectangular dump profiles cover the entire dump face (or a large portion of it) evenly to a uniform width. Sloughed heap profiles occur when a portion of the material is not dumped over the edge of the berm, but rather on the floor of the upper level of the dump. This causes a portion of the material to bunch near the dump berm, and a part of it to slough over the edge of the dump. Table 4 shows example grayscale renderings of what is observed in PointStudio for each of the four types of dumping profiles determined by a qualitative analysis of the dump HFMs. These example grayscale renderings are showcased via both front and side views. The front view shows the perspective facing the dump and in-line/parallel to the dump, with the dumping process proceeding from top to bottom. The side view represents the perspective perpendicular to the front view, with the dumping direction proceeding from top right to bottom left. Measurements of these profiles were performed using software functions found within PointStudio.
The HFM solids created in PointStudio were classified using the following variables: volume (in cubic meters), maximum height, length, width and thickness (in meters) and angle (in degrees) from bottom left to top right. Volume was determined by querying the properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section.
In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all dumps were along the edge of a 30 m high dump crest. Table 5 reveals that large volume does not always indicate large maximum width, height or length. High angles do not indicate large volume, but seem to increase the maximum width, height or length. In the case of sloughed heaps, the low angle is due to the fact that it is not measuring to the top of the berm, but to the extent of the dump profile, which occurs on the floor of the upper level of the dump area. properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum lengt were considered to be the vertical and horizontal legs of the right triangle formed b connecting the bottom left and top right points of the HFM solid. Maximum width wa considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between th two farthest horizontal points of the HFM solid from the side-view perspective. Angl was considered to be the inner angle of the right triangle formed by connecting th bottom-left and top-right points of the HFM solid. These measurements are tabulated i Table 5. An additional visualization of the HFMs data containing all of the extracted solid is available in the data cache associated with this article. The link accessed to this data ca be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 show information from all of the dump profiles. With the exception of dumps 1, 2 and 3, a properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum lengt were considered to be the vertical and horizontal legs of the right triangle formed b connecting the bottom left and top right points of the HFM solid. Maximum width wa considered to be the distance between the two farthest horizontal points of the HFM soli from the front-view perspective. Thickness was considered to be the distance between th two farthest horizontal points of the HFM solid from the side-view perspective. Angl was considered to be the inner angle of the right triangle formed by connecting th bottom-left and top-right points of the HFM solid. These measurements are tabulated i Table 5. An additional visualization of the HFMs data containing all of the extracted solid is available in the data cache associated with this article. The link accessed to this data ca be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 show information from all of the dump profiles. With the exception of dumps 1, 2 and 3, a properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all properties of the HFM solids after their creation. Maximum height and maximum lengt were considered to be the vertical and horizontal legs of the right triangle formed b connecting the bottom left and top right points of the HFM solid. Maximum width wa considered to be the distance between the two farthest horizontal points of the HFM soli from the front-view perspective. Thickness was considered to be the distance between th two farthest horizontal points of the HFM solid from the side-view perspective. Angl was considered to be the inner angle of the right triangle formed by connecting th bottom-left and top-right points of the HFM solid. These measurements are tabulated i Table 5. An additional visualization of the HFMs data containing all of the extracted solid is available in the data cache associated with this article. The link accessed to this data ca be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 show information from all of the dump profiles. With the exception of dumps 1, 2 and 3, a properties of the HFM solids after their creation. Maximum height and maximum length were considered to be the vertical and horizontal legs of the right triangle formed by connecting the bottom left and top right points of the HFM solid. Maximum width was considered to be the distance between the two farthest horizontal points of the HFM solid from the front-view perspective. Thickness was considered to be the distance between the two farthest horizontal points of the HFM solid from the side-view perspective. Angle was considered to be the inner angle of the right triangle formed by connecting the bottom-left and top-right points of the HFM solid. These measurements are tabulated in Table 5. An additional visualization of the HFMs data containing all of the extracted solids is available in the data cache associated with this article. The link accessed to this data can be found in the "Supplementary Materials" section. In total, 28 dumps were profiled and classified via qualitative analysis. Table 5 shows information from all of the dump profiles. With the exception of dumps 1, 2 and 3, all 10 m 10 m 10 m In order to investigate the statistical differences among the classification data of each type, box plots of their data for each variable are shown in Figures 2-7. In the box plots, the black dot represents the mean average value within the data, the line represents the median value, the box edges represent 50% of the data between the first and third quartiles and the lines above and below the box represent the maximum and minimum values of the data.
Mining 2022, 2, FOR PEER REVIEW 10 dump types seem to involve the movement of additional material contained in the dump face.          Figure 5 shows the same general differences between dump types as Figure 4. This is likely due to the fact that length is a measurement of how much of the horizontal dimension of the dump is covered by the dump profile. Length values are lower for sloughed heap-and comet-type dump profiles compared to the other dump types. This may be related to the same reason why height values for these types are also lower.   Figure 5 shows the same general differences between dump types as Figure 4. This is likely due to the fact that length is a measurement of how much of the horizontal dimension of the dump is covered by the dump profile. Length values are lower for sloughed heap-and comet-type dump profiles compared to the other dump types. This may be related to the same reason why height values for these types are also lower.  As shown in Figure 7, the comet and sloughed heap dump profiles typically have a higher thickness than the oval and rectangular ones. Assuming that the dumps are of similar total volume, this is to be expected, since oval and rectangular dump profiles typically have higher height and length. Therefore, the dump material for oval and rectangular dumps is spread thinner across more surface area, which leads to less thickness. Comet dump profiles may also interact with loose material on the dump surface, and the resulting solid might include some of that material in the thickness, as it might with the total volume.   As shown in Figure 7, the comet and sloughed heap dump profiles typically have a higher thickness than the oval and rectangular ones. Assuming that the dumps are of similar total volume, this is to be expected, since oval and rectangular dump profiles typically have higher height and length. Therefore, the dump material for oval and rectangular dumps is spread thinner across more surface area, which leads to less thickness. Comet dump profiles may also interact with loose material on the dump surface, and the resulting solid might include some of that material in the thickness, as it might with the total volume.   Figure 2 shows that the comet and rectangular dump types have a slightly higher average volume compared with the oval and sloughed heap types. Volume is influenced by how much material was loaded into the truck before dumping, and also by how much material in the dump face moved at the time of dumping. This increase in the average volume of material for the comet and rectangular dumps may be due to the fact that both dump types seem to involve the movement of additional material contained in the dump face. Figure 3 shows a lot of similarity in the average angles of the comet, oval and rectangular dump types. Sloughed heap angles are lower because the material does not typically extend the full length of the dump face. Where the material of the sloughed heap extends along the dump face, the angle of the material matches the angles for the other dump types. Figure 4 shows that the oval and rectangular dumps have the highest dump height. Dump height is a measurement of how much of the vertical dimension of the dump is covered by the dump profile. Sloughed heap and comet height values are lower because the material does not typically extend the full length of the dump face. Figure 5 shows the same general differences between dump types as Figure 4. This is likely due to the fact that length is a measurement of how much of the horizontal dimension of the dump is covered by the dump profile. Length values are lower for sloughed heapand comet-type dump profiles compared to the other dump types. This may be related to the same reason why height values for these types are also lower. Figure 6 shows that comet dump profiles tend to be the widest. This may be the result of additional material from the dump face aggregating with the dump mass as it cascades, resulting in an increase in width. Rectangular dump profiles tend to have the lowest width values, which is interesting because they typically have the highest volume values, and this might be explained by a low amount of frictional resistance on the dump face compared to the cohesion of the dump mass.
As shown in Figure 7, the comet and sloughed heap dump profiles typically have a higher thickness than the oval and rectangular ones. Assuming that the dumps are of similar total volume, this is to be expected, since oval and rectangular dump profiles typically have higher height and length. Therefore, the dump material for oval and rectangular dumps is spread thinner across more surface area, which leads to less thickness. Comet dump profiles may also interact with loose material on the dump surface, and the resulting solid might include some of that material in the thickness, as it might with the total volume.

Modelling Discussion
The HFMs presented in this paper offer a novel look at individual dumps. These HFMs can be used as a basis to calibrate future simulation models of dumps via the V-model for calibration described by Quist [63], Hofmann [66] and others. Figure 8 illustrates the V-model for simulation calibration and validation, which occurs over three levels. At the bottom level, validation is carried out by calibrating the individual parameters of particles in a laboratory setting. The second level involves calibration based on the aggregated behavior of the particles through multiple flow regimes. The third level compares simulated outcomes with real industrial-scale operations.
OR PEER REVIEW 13 4. Discussion

Modelling Discussion
The HFMs presented in this paper offer a novel look at individual dumps. These HFMs can be used as a basis to calibrate future simulation models of dumps via the Vmodel for calibration described by Quist [63], Hofmann [66] and others. Figure 8 illustrates the V-model for simulation calibration and validation, which occurs over three levels. At the bottom level, validation is carried out by calibrating the individual parameters of particles in a laboratory setting. The second level involves calibration based on the aggregated behavior of the particles through multiple flow regimes. The third level compares simulated outcomes with real industrial-scale operations. To the knowledge of the authors, no work has been conducted to validate dumping behavior at a real mine against simulated outcomes; however, much work has been performed to calibrate and validate ROM material simulations at lower levels of the V-model, particularly for mineral processing applications [70]. In instances such as dumping, where reality cannot be used to calibrate models, calibration involves adjusting parameters to fit HFMs instead. There is no mention of a method for creating HFMs of mine haul truck dumping activity in the literature. The method used for this study is easily achievable for any mine operation where UAVs and photogrammetry are presently in use. To the knowledge of the authors, no work has been conducted to validate dumping behavior at a real mine against simulated outcomes; however, much work has been performed to calibrate and validate ROM material simulations at lower levels of the V-model, particularly for mineral processing applications [70]. In instances such as dumping, where reality cannot be used to calibrate models, calibration involves adjusting parameters to fit HFMs instead. There is no mention of a method for creating HFMs of mine haul truck dumping activity in the literature. The method used for this study is easily achievable for any mine operation where UAVs and photogrammetry are presently in use.

Dump Profile Discussion
Even though this study found only four types of dump profiles, there may exist additional types. The factors that determine which type a dump will be categorized into are many. These factors likely fall into several categories, including the way the material is loaded into the truck, how the truck dumps the material, how the pre-existing dump face interacts with the dumped load and how the material behaves on its own. None of these factors were investigated in this paper How the truck is loaded directly influences the volume of the resulting dump profile. It may also influence the determination of the dump profile type, forming either comet, oval or rectangular profiles. For example, if the truck is unevenly loaded with more material in the back than towards the front, this may cause the resulting dump profile to favor a comet format. However, the exact interplay between how the trucks were loaded and the resulting dump profiles remains unclear, and no information on truck loading was gathered during this study.
How the truck dumps the material clearly influences whether or not the dump profile becomes a sloughed heap or one of the other types. This is because if the truck dumps far from the crest of the dump face, it will create a sloughed heap. When the truck dumps against the crest of the dump face, the type of the resulting dump profile is either comet, oval or rectangular; however, it is unclear which one it will become from this information alone.

Terminology Discussion
The lack of terminology around dump profile behavior is one challenge to improving our understanding. The authors present new terminology here for the four types of dump shapes, but much more terminology is likely needed to fully describe and characterize dump profiles. This is made true by the fact that many dump profiles display characteristics of more than one type and, therefore, additional description may be required to fully classify them. For example, a sloughed heap may slough into a comet shape at the bottom of the dump, or an oval may have an extremely long rectangular section. There also needs to be reflection around how much effort should be spent creating terminology, since it is well known that the shape of mine rock piles is mainly based on topography [9].

Practicality Discussion
Admittedly, some findings in this study may not offer much practical applicability, especially considering the status quo of the mining industry. It is unknown whether the HFMs presented in this study represent the global characteristics of the dumping process across all mines, or whether they are limited to the mine used for the study. However, the method used for data capture was seamlessly incorporated into routine operation, and was practical from the standpoint of simplicity and ease of realization. Many mines are capable of measuring their own dump characteristics and creating HFMs for themselves. It is entirely probable that many other shapes of fallow land exist, and their discovery and classification can increase our understanding of the cascading process of ROM material.
Additionally, as has been previously stated, the findings of this study are practical for the purpose of verifying the accuracy of the simulation modelling of individual truck end dumps. Without an HFM to verify simulations against, there would be less confidence in the accuracy and relevance of the simulation. It is computationally intensive to simulate an entire dump area. Simulating individual dumps may allow for larger areas to be simulated with less computational power through the use of pseudo-particles [4].
Dumps generally conform to whatever shape the local topography provides, but knowledge of the common shapes provided via a routine dumping process for a particular ore at a given mine has historically been mostly speculatory. Speculation of this kind causes engineers and mine planners to place an unknown amount of dependance on operators that work with the material constantly to ensure the process is conforming to plan. In future autonomous mining scenarios, this dependance will not be allowed. Thus, digitally transforming the dumping process is hypothesized to support continuity in the transition towards autonomous mining. Without operators at the helm, there is little knowledge to ensure that dumping is occurring correctly. Additionally, without a basic reference, there would be no way to determine the performance quality of the autonomous equipment. The fact that autonomous equipment will be covered in sensors to map the work area continuously will mean very little, unless it is known what conformity should actually looks like.
As another hypothesis, there might be a relationship of practical synergy between predictive modelling and UAV surveys. UAV surveys have become increasingly ubiquitous at mines, and can accomplish tasks to a level of quality unachievable by traditional survey crews within a greatly reduced timeframe. They could be made to occur at such frequent intervals as to overwhelm mining engineers and long-range planners. The HFMs described in this article demonstrate the parameterization of the dumping process, which could facilitate the training and validation of predictive models to help automate design conformity to UAV surveys of rock piles, thereby decreasing the cognitive load placed on the domain experts.

Conclusions
HFMs of dump profiles for 28 dumps were created. These HFMs show characteristic behaviors classifiable into four types, named comet, oval, rectangular and sloughed heap. These classifications may make it easier to examine the dumping process as a whole. While more terminology and modelling will be required to gain a complete understanding of the dumping process, the HFMs examined here provide a basis to articulate new terminology and calibrate new modelling.
Since further investigation is required, some recommendations for future study in this area include: