Stochastic Final Pit Limits: An Efﬁcient Frontier Analysis under Geological Uncertainty in the Open-Pit Mining Industry

: In the context of planning the exploitation of an open-pit mine, the ﬁnal pit limit problem consists of ﬁnding the volume to be extracted so that it maximizes the total proﬁt of exploitation subject to overall slope angles to keep pit walls stable. To address this problem, the ore deposit is discretized as a block model, and efﬁcient algorithms are used to ﬁnd the optimal ﬁnal pit. However, this methodology assumes a deterministic scenario, i.e., it does not consider that information, such as ore grades, is subject to several sources of uncertainty. This paper presents a model based on stochastic programming, seeking a balance between conﬂicting objectives: on the one hand, it maximizes the expected value of the open-pit mining business and simultaneously minimizes the risk of losses, measured as conditional value at risk, associated with the uncertainty in the estimation of the mineral content found in the deposit, which is characterized by a set of conditional simulations. This allows generating a set of optimal solutions in the expected return vs. risk space, forming the Pareto front or efﬁcient frontier of ﬁnal pit alternatives under geological uncertainty. In addition, some criteria are proposed that can be used by the decision maker of the mining company to choose which ﬁnal pit best ﬁts the return/risk trade off according to its objectives. This methodology was applied on a real case study, making a comparison with other proposals in the literature. The results show that our proposal better manages the relationship in controlling the risk of suffering economic losses without renouncing high expected proﬁt.


Introduction
In simple terms, the final pit limit problem is defined as determining the ultimate mining limits of an ore deposit exploited by the open-pit mining method such that a maximum undiscounted profit is obtained from its extraction, respecting the precedence constraints given by overall slope angles that ensure stable walls of the open-pit mine.
The final pit limit plays a crucial role in long-term open-pit mine production planning. It approximates the placement, size, shape, and depth of the mine at the end of its exploitation, determining important aspects for mine operation such as the layout of access roads, ramps, waste dumps, stockpiles, processing, and other facilities, as well as in developing a production schedule. Besides, different studies such as feasibility analyses, assessment of capital exposure, and corporate risk can be considerably affected by the results of non-optimum final pit limit determination [1].
As a first step, the ore deposit is modelled as a regular array of blocks, known as the block model. This model is constructed by interpolating information available at drill hole samples obtained from the terrain. Several attributes are estimated, such as rock types, concentrations (grades) of relevant elements, and density, among others, using geostatistical estimation techniques [2]. Then, an economic value that represents its profit is computed for each block. This value depends on the estimated geological attributes and economic and operational parameters such as ore price, costs, and metallurgical recoveries. The result of this is called the economic block model and is the primary input for pit limit optimization [3].
The final pit decision is contemplated in the initial stages of strategic open-pit mine planning and can be solved efficiently by using Lerchs-Grossmann [4] or Pseudoflow [5,6] algorithms. However, as indicated before, it is based on a single, smooth representation of the deposit [7]; therefore, the in-situ geological uncertainty is not taken into account. This type of uncertainty comes from an incomplete knowledge of the ore deposit considering a reduced number of drill hole samples, restricted by exploration costs, limiting the estimation accuracy of the geological attributes. According to [8], the quality of the estimation depends mainly on the following: (i) the placement, quantity, and quality of the samples, (ii) the orebody classification, and (iii) the method used to generate the estimations.
As mine exploitation progresses, new sample information is generated, and the block model is updated. However, strategic mine planning decisions (such as mining sequence) have already been made, causing deviations in the planned production promises. As a result of the traditional approach, mining companies are often incapable of fulfilling production promises due to their inability to evaluate and minimize the consequences of geological uncertainty in the early stages of a mining project [9][10][11][12].
Conditional simulation is a geostatistical method for modelling geological uncertainty and a valuable tool for risk assessment in pit designs, generating equally probable scenarios (realizations) that represent the in-situ orebody variability [13][14][15]. There are many approaches for conditional simulation of continuous random variables, such as ore grades as used in this work. Some are based on a normal score transformation and the assumption of multi-Gaussianity. These include sequential Gaussian simulation [2,16,17], turning bands simulation [18] and matrix decomposition simulation [19]. Other methods rely on an approximation of the conditional distribution done by defining thresholds and creating indicator variables, which are characterized through sequential indicator simulation [20]. In recent years, more advanced methods that use multiple-point statistics have been developed [21,22] to simulate, accounting for patterns. Simulated annealing [23], direct sequential simulation [24], and p-field simulation [25] provide other approaches that do not require the same assumptions of the abovementioned methods.
Assuming there is a set of scenarios characterizing the geological uncertainty, for instance rock type and ore grades, the first challenge in mine planning is how to use this information to assess its impact on the plan. One option could be to obtain one final pit per scenario separately, and from them, decide of a single final pit, based on some well-defined criteria ( Figure 1). A completely different option is a model that does not consider the input data (set of scenarios) separately but can consider them simultaneously in one run and thus compute a single robust optimal final pit ( Figure 2). However, because the final pit is one of the most critical and higher impact decisions in a mining project, some questions arise: • What final pit should be chosen for a given confidence level in terms of economic value or metal content from all possible options? or, • How would the final pit contour change if the ore grade variability is higher or lower than expected?
Unfortunately, the traditional methodology to compute final pit limits cannot give answers to these critical questions. For this reason, it is crucial to generate robust final pit methodologies that consider geological uncertainty, especially at the beginning of a mining project, when a high level of uncertainty is present.

Figure 1.
A first approach to compute final pit limit: (a) a set of realizations representing logical uncertainty are used as input individually, (b) the optimization computes a final pit realization, (c) a set of partial solutions are given, and (d) applying some criteria, a single limit is selected.

Figure 2.
A second approach to compute final pit limit: (a) a set of realizations represen geological uncertainty are used as input simultaneously, (b) the optimization computes solution for the entire set of realizations, and (c) a robust final pit limit is returned.
However, because the final pit is one of the most critical and higher impact d in a mining project, some questions arise: • What final pit should be chosen for a given confidence level in terms of ec value or metal content from all possible options? or, • How would the final pit contour change if the ore grade variability is higher o than expected? Unfortunately, the traditional methodology to compute final pit limits cann answers to these critical questions. For this reason, it is crucial to generate robust methodologies that consider geological uncertainty, especially at the beginning o  However, because the final pit is one of the most critical and higher impact d in a mining project, some questions arise: • What final pit should be chosen for a given confidence level in terms of ec value or metal content from all possible options? or, • How would the final pit contour change if the ore grade variability is higher o than expected? Unfortunately, the traditional methodology to compute final pit limits cann answers to these critical questions. For this reason, it is crucial to generate robust methodologies that consider geological uncertainty, especially at the beginning o One of the first efforts to include geological uncertainty into the final pit limit is by [26], which presented a hybrid approach within a set-theory framework. This work considered final pits for several simulated realizations of the orebody and then defined the so-called hybrid pit. Later, Refs. [10,27,28] compute pits based on the hybrid approach and propose some reliability indexes for a pit. Ref. [29] proposed to use one of the simulated pits as the best pit design, capturing the upside potential of the orebody and minimizing the potential downside risk. Ref. [30] presented a workflow that includes the computation of a series of nested pits, considering the simulations of the block model and distributions of economic and geotechnical parameters as inputs.
The above approaches can be used as a guide to assess the risk associated with geological uncertainty; however, their main drawback is that the optimization process does not incorporate risk control. The solution is only obtained after a series of final pit limit problems are solved separately, following the approach presented in Figure 1.
On the other hand, under the approach showed in Figure 2, we can identify the following contributions: Ref. [31] presented the advantages of using a stochastic approach based on expected economic values rather than expected ore grades, from a set of simulated realizations of ore grade. Ref. [32] proposed a constrained version of the problem, incorporating the risk management directly into the optimization model in a probabilistic framework. In the same line, refs. [33,34] incorporated grade uncertainty on a modified version of the model for the final pit limit problem, where extraction and processing decisions are taken separately, proposing a two-stage scheme. Ref. [35] used risk measures to incorporate geology and market uncertainty into the investment, pit design, and mining sequence decisions, but unfortunately, they do not perform the methodology on a real case study.
More recently, Ref. [36] develops a mean-variance criterion approach to finding nearoptimal final pit limits by using a heuristic procedure considering geological uncertainty, generating a set of solutions on the mean-variance efficient frontier. Then, based on some stochastic dominance rules, the authors eliminate sub-optimal solutions. Some drawbacks of this approach are (i) standard deviation as a risk measure does not consider skewness and (ii) the proposed heuristic algorithm cannot guarantee an optimal solution because it uses a greedy optimizer, and, besides, some parameters of the algorithm must be chosen manually. A second work, ref. [37], focuses on the theoretical problem and proposes properties that a risk-averse measure for final pit limit should have (i) nestedness of pits for different risk aversion levels and (ii) additive consistency, which states that extraction order should not be affected by the precedence of farther parts of the mine. They show that only an entropic risk measure complies with these properties and propose an approximation scheme based on a two-stage stochastic model. Finally, interesting research was developed by [38], who presented a comprehensive overview of the applications and trends of multi-criteria decision-making methods applied in mining and mineral processing problems.
In this paper, we propose a methodology based on a stochastic optimization model that maximizes expected profit while controlling maximum risk in terms of conditional value at risk in mining projects and apply it to a real case study. The orebody is modelled through several simulated realizations to incorporate explicit risk control in open-pit mine design and evaluation. As a result, the efficient frontier of final pit limit alternatives in the expected profit vs. risk environment is obtained. Additionally, we propose some criteria to choose pits from the efficient frontier depending on the decision maker's risk aversion.

Materials and Methods: Modelling, Notation, and Problem Statement
In this section, we present the main notation and formulation for the final pit limit problem. Section 2.1 presents the model for the deterministic case of the final pit, and Section 2.2 introduces basic concepts of conditional value at risk. Based on these preliminary concepts, Section 2.3 states the problem addressed in this work by providing a formulation to compute the final pit considering a representation of geological uncertainty explicitly.

Final Pit Limit: Deterministic Case
Let B be the block model and b a block. PREC b ⊂ B − {b} denotes the subset of blocks above block b that must be extracted before it, respecting maximum slope angles. An economic block value v b is precomputed for block b. By using integer programming, the decision variables are defined as Thus, the final pit limit is the solution to the following where Equation (2) represents the maximum undiscounted total profit of blocks within final pit limits. Equation (3) ensures that the shape of the ultimate pit respects the overall slope angle, and Equation (4) indicates the decision variables are binary. There exist fast algorithms to solve this problem [5,6,39]. As stated above, the geological uncertainty is not considered in this formulation. Therefore, this model does not capture the risk associated with a non-accurate ore resource estimation.

Background on Risk Measures: Conditional Value at Risk
In mining under geological uncertainty, the term risk refers to the potential economic losses caused by a misleading ore resource estimation. One tool used to reduce the risk of losses is the value at risk (VaR). Let δ ∈ (0, 1) be the risk level and f (x, y) a loss function for a decision vector x ∈ X and a random vector y ∈ R m . For a given x, f (x, y) is a random variable having a distribution into R induced by y. Ψ(x, ζ) is the probability that the loss function f (x, y) does not exceed a threshold value ζ, that is, Ψ(x, ζ) is the cumulative distribution for the loss associated to x, which depends on ζ considering x as fixed. Thus, However, this risk measure has some drawbacks, as reported in [40], such as lack of sub-additivity or lack of control by losses exceeding VaR. Thus, Rockafellar and Uryasev [41] propose a different measure of risk, namely the conditional value at risk (CVaR), which represents the conditional expected loss given that the loss exceeds VaR. For continuous distributions where p(y) is the density of y. To avert difficulties with VaR δ (x) into Equation (6), Rockafellar and Uryasev give an alternative representation that characterizes both VaR δ (x) and CVaR δ (x) in terms of F δ on X × R: quantifying the losses that exceed VaR, acting as an upper bound for it at the same confidence level δ. Besides, CVaR is convex and a coherent risk measure [42], allowing the implementation of optimization algorithms to determine it. The integral in Equation (7) can be approximated in several ways: for instance, by sampling the probability distribution of y according to its density p(y) from which a collection of R equally probable vectors y 1 , . . . , y R is obtained. Therefore, the corresponding approximation to F δ (x, ζ) is where [t] + = max{t, 0}. As shown in [41], an approximation of minimum CVaR is found by minimizing Equation (8). Note that the optimization simultaneously gives the optimal decision x * , VaR δ (x * ), and CVaR δ (x * ). Henceforth, when we refer to CVaR, we strictly mean the approximation to CVaR, since geological uncertainty is represented employing a finite number of scenarios. Finally, to avoid the non-linear function [ f (x, y r ) − ζ] + into Equation (8), we will use auxiliary variables z r , for all r ∈ R = {1, . . . , R}, Further details on minimization of CVaR can be found in [41]; for properties and strong and weak features of the VaR and CVaR risk measures, see [42]; and for illustrations with several examples, see [43]. Applications on optimization problems, including CVaR constraints, are addressed in [44].

Problem Statement: Risk-Based Final Pit Limit
Consider several realizations y 1 , . . . , y R of random variable y from resource models indexed by r ∈ R to represent the geological uncertainty. For example, y r may be interpreted as the ore grade or metal content according to realization r. In this case, An economic block model can be stated per realization, i.e., v br indicates the block value according to the metal content y br from realization r ∈ R and block b ∈ B. Contrary to the above, the traditional methodology proposes that a single economic value v b is obtained per block when the metal content y b is estimated by using estimation methods such as inverse distance, kriging, or e-type [2].
As before, the decision vector x = (x b ) b∈B corresponds to the selection of blocks within the final pit limit. Regarding the loss function for each realization r ∈ R, we propose to measure the losses (or gains) associated with deviations of economic values from each realization regarding an economic value provided according to one of the deterministic/traditional methods. Therefore, the loss function may be expressed as Finally, considering all the above notation, we present the stochastic optimization model that seeks the maximum expected profit under CVaR constraints.
Equation (12) shows the objective function, maximizing expected undiscounted profit of the final pit limit. Equation (13) establishes slope precedence among blocks, and Equation (14) limits the maximum risk µ ≥ 0 allowed in terms of CVaR by considering a confidence level δ ∈ (0, 1) along R geological scenarios. Equations (15) and (16) are imposed by auxiliary variables, and Equation (17)  Note that an optimal solution (x * , ζ * , z * ) to P µ for a given µ ≥ 0 determines the blocks inside the final pit limits, the expected undiscounted profit, both VaR (ζ * ) and CVaR, and how the impact of each scenario r into CVaR is distributed.
Varying the parameter value µ in the model P µ allows analysing the trade-off between the (maximum) expected profit of the final pit limit and the (minimum) risk of loss in terms of CVaR. A procedure for defining the set of parameters µ to be considered is as follows: Procedure 1.

•
Find the maximum risk level: to solve the maximum-profit problem, ignoring any constraints on the level of risk, thus obtaining a maximum risk level µ max as CVaR.

•
Find the minimum risk level: to solve a minimum-risk problem, ignoring any profit constraint, thus obtaining a minimum risk level µ min as CVaR, and • Partition the set of possible risk levels: to solve the problem P µ for a set of values Figure 3a illustrates the shape of the curve obtained by varying the parameter µ in the expected profit vs. risk plane. Because no final pit above the frontier is possible and alternatives below the frontier are sub-optimal, this curve defines the efficient frontier of final pit alternatives, whether with a higher expected profit keeping the risk bounded or one with lower risk but without reducing profit (see Figure 3b). This boundary allows the investor to make the best decision on which is the final pit of maximum profit, for a given risk level, or to minimize the risk when a particular return is expected.
Equation (12) shows the objective function, maximizing expected undiscounted profit of the final pit limit. Equation (13) establishes slope precedence among blocks, and Equation (14) limits the maximum risk ≥ 0 allowed in terms of CVaR by considering a confidence level ∈ (0, 1) along geological scenarios. Equations (15) and (16) are imposed by auxiliary variables, and Equation (17) denotes the nature of decision variables.
Note that an optimal solution ( * , * , * ) to ( ) for a given ≥ 0 determines the blocks inside the final pit limits, the expected undiscounted profit, both VaR ( * ) and CVaR, and how the impact of each scenario into CVaR is distributed.
Varying the parameter value in the model ( ) allows analysing the trade-off between the (maximum) expected profit of the final pit limit and the (minimum) risk of loss in terms of CVaR. A procedure for defining the set of parameters to be considered is as follows: Procedure 1.

•
Find the maximum risk level: to solve the maximum-profit problem, ignoring any constraints on the level of risk, thus obtaining a maximum risk level as CVaR.

•
Find the minimum risk level: to solve a minimum-risk problem, ignoring any profit constraint, thus obtaining a minimum risk level as CVaR, and • Partition the set of possible risk levels: to solve the problem ( ) for a set of values Figure 3a illustrates the shape of the curve obtained by varying the parameter in the expected profit vs. risk plane. Because no final pit above the frontier is possible and alternatives below the frontier are sub-optimal, this curve defines the efficient frontier of final pit alternatives, whether with a higher expected profit keeping the risk bounded or one with lower risk but without reducing profit (see Figure 3b). This boundary allows the investor to make the best decision on which is the final pit of maximum profit, for a given risk level, or to minimize the risk when a particular return is expected. Some of the criteria that a decision maker could use to choose a specific final pit limit from the frontier are the following: 1. (C1) The ideal point is defined as the unfeasible solution located above the frontier, where each criterion is optimized separately (maximum expected return and minimum risk). This criterion proposes to select the final pit that presents the minimum distance to the ideal point. Some distance measures available are Euclidean, normalized, and Manhattan, among others [45]. 2. (C2) To select the final pit that reports the most significant difference between expected return and CVaR. Some of the criteria that a decision maker could use to choose a specific final pit limit from the frontier are the following: 1.
(C1) The ideal point is defined as the unfeasible solution located above the frontier, where each criterion is optimized separately (maximum expected return and minimum risk). This criterion proposes to select the final pit that presents the minimum distance to the ideal point. Some distance measures available are Euclidean, normalized, and Manhattan, among others [45]. 2.
(C2) To select the final pit that reports the most significant difference between expected return and CVaR. 3.
(C3) A priori, it is possible that due to the company's internal policies, the decisionmaker may not be able to accept projects with a higher risk than a given value µ 0 . In this case, the criterion suggests selecting the pit alternative with a higher value and risk limited by µ 0 .

4.
(C4) Contrary to the previous case, the decision-maker may want to minimize the risk but ensure a minimum return V 0 . In this case, the criterion proposes to select the final pit with a minimum CVaR so that the expected return is greater than or equal to V 0 .
The interesting results in the efficient frontier are found in the p-range (points of the efficient frontier that are in the p% of a more significant expected profit zone). In general, the final pit limit decision will be taken within this range, looking for the alternative that maintains a high expected profit and reduces the associated risk. Finally, note that it is essential to avoid comparing VaR and CVaR values for the same level of confidence, because they indicate information from different parts of the loss distribution.

Implementation and Results
The methods presented in Section 2 are applied to a real porphyry copper deposit as a case study (Section 3.1). The experiments are as follows: (i) generating several optimum final pit alternatives (efficient frontier) for a family of parameters µ (Section 3.2), from a set of simulated realizations, and (ii) comparing in terms os expected profit vs. risk environment, the efficient frontier with the final pit obtained by using the traditional methodology (deterministic case) and other approaches as presented in Figure 1 (Section 4).
To solve the instances of the model P µ , GUROBI was used as the optimization solver, version 9.0.3. The code's execution was performed on an Intel Xeon CPU E5-2660 V3 machine with 10 cores at 2.6 GHz, 128 Gb and running Windows 10 Pro.

Dataset and Instances
The block model corresponds to a porphyry copper deposit (known as BMT, for confidentiality reasons). This block model has 407,179 blocks of 10 m × 10 m × 10 m and contains information about spatial coordinates, density, and 50 conditional simulations of copper grade generated by a sequential Gaussian simulation algorithm. In this method, the original samples are transformed to normal scores, and a standard Gaussian random variable is simulated in a spatial grid. Each node of this grid is visited sequentially in a random order, performing a simple kriging estimation using the previously simulated nodes and the normal scores of the true data to condition its value, and drawing a random number from a Gaussian distribution with mean and variance given by the kriging estimate and variance, respectively, to obtain the simulated value, which is then back-transformed to the original grade scale [46]. We assume that the 50 simulated realizations are sufficient to capture the actual variability of copper grades, and we use the E-Type model (average of the 50 realizations) as the standard representation of the metal content estimation considered in the deterministic approach. Figure 4 shows three different simulated realizations of the copper grade, in addition to the average grade for a plan-view at z = 2015 m. Figure 5a shows the average histogram of copper grades, including error bars along all sets of realizations. Finally, grade-tonnage curves are plotted to quantify the recoverable resources to different cut-off grades (Figure 5b). Both figures represent error bars with the 5th and 95th percentile per interval, showing low uncertainty in copper grades and ore tonnages.
Precedence arcs and economic block values are provided. In this case, slope precedence requirements respect a 45 • overall slope angle in pit walls, but the model can be implemented with other slope angles.    Precedence arcs and economic block values are provided. In this case, slope precedence requirements respect a 45° overall slope angle in pit walls, but the model can be implemented with other slope angles.

Results from Stochastic Model
In this section the results for stochastic final pit limit obtained from model ( ) are shown considering a = 95% confidence level in risk assessment through the CVaR. Applying Procedure 1, Table 1 shows the main numerical results obtained (columns indicate from left to right): the parameter that determines the maximum risk allowed according to Equation (18), VaR, CVaR, expected profit, and the distance (Euclidean metric) to the ideal point (DIP), which is presented as a decision criterion.

Results from Stochastic Model
In this section the results for stochastic final pit limit obtained from model P µ are shown considering a δ = 95% confidence level in risk assessment through the CVaR. Applying Procedure 1, Table 1 shows the main numerical results obtained (columns indicate from left to right): the parameter α that determines the maximum risk µ allowed according to Equation (18), VaR, CVaR, expected profit, and the distance (Euclidean metric) to the ideal point (DIP), which is presented as a decision criterion. All results are within an optimality gap of 1%. Note that the final pit decision yields an empty pit when looking for the minimum risk, i.e., it is equivalent to perform no mining operation. These results allow us to plot the pair (CVaR, expected profit) forming the efficient frontier, or Pareto front (Figure 6a): the points represent different final pit alternatives, each with its respective expected profit and risk. The ideal point (brown) is plotted as a reference according to criterion (C1).  The interesting results on the efficient frontier are found in the p-range (as explained above). Considering a 20-range, reaching a reduction of up to 20% of expected profit, the final pit alternatives from parameters µ = 200 [MUSD] to µ = 75 [MUSD] are inside this range. The limit case (µ = 75 [MUSD]) includes the final pit whose risk is reduced by approximately 63%, while expected profit drops almost 16%. Regarding the expected tonnages of ore and waste material for each final pit alternative on the efficient frontier, Figure 6 shows how expected ore and waste tonnages decrease as the parameter µ decreases: this should be interpreted as a renunciation of riskier areas within the deposit. Ore tonnage and average grade deviations are represented through error bars showing both the 5th and 95th percentiles along with all geological scenarios. Figure 7 shows the contribution of z r to CVaR through R. These contributions are the ones that generate the differences between columns VaR and CVaR in Table 1. Note that the way scenarios contribute to the CVaR is not intuitive because the decision is strongly dependent on the upper bound µ on Equations (14) and (15). In this case study, not all geological scenarios contribute to CVaR in the optimization process; only 11 out of 50 scenarios present total net losses, with realization indexed by r = 5 being the one that presents the highest loss. Other scenarios show losses as well, decreasing the total contribution (in MUSD) as the parameter µ decreases, but the final pit with parameter µ = 120 presents VaR = CVaR; therefore, ∑ r∈R z r = 0. Finally, Figure 8 shows plan-views and cross-sections for some final pits, particularly those obtained with parameters = 200, 125, 25.
[MUSD]: significant differences in size and shape among the obtained pits can be observed, especially at the bottom of the pits.  Finally, Figure 8 shows plan-views and cross-sections for some final pits, particularly those obtained with parameters = 200, 125, 25.
[MUSD]: significant differences in size and shape among the obtained pits can be observed, especially at the bottom of the pits.

Discussion
Once the efficient frontier is determined, it is interesting to know its relative position regarding several final pit decisions obtained from other strategies. In this section we present a comparison of our approach based on the maximization of expected profit but controlling risk, with other approaches: (i) E-Type, as the only deterministic representation of the deposit; (ii) Expected-profit, as shown in [31]; (iii) Best-simulation as pit design, similar to [29]; and (iv) Hybrid-pit, as presented in [26]. All these results are located in the feasible region within the profit-risk plane. For comparison purposes, the 'best' final pit from the efficient frontier will be chosen according to the criterion (C1) from Section 2.3,

Discussion
Once the efficient frontier is determined, it is interesting to know its relative position regarding several final pit decisions obtained from other strategies. In this section we present a comparison of our approach based on the maximization of expected profit but controlling risk, with other approaches: (i) E-Type, as the only deterministic representation of the deposit; (ii) Expected-profit, as shown in [31]; (iii) Best-simulation as pit design, similar to [29]; and (iv) Hybrid-pit, as presented in [26]. All these results are located in the feasible region within the profit-risk plane. For comparison purposes, the 'best' final pit from the efficient frontier will be chosen according to the criterion (C1) from Section 2.3, that is, by using parameter µ = 125 [MUSD]. However, the same analysis can be generalized by considering other points on the efficient frontier.
In the cases of (i) E-Type, and (ii) Expected-profit, an optimization program (P) defined by Equations (2)-(4) (differentiated by the economic block models) is used to compute a single final pit limit. To find the position of each final pit in the expected profit vs. risk plane, we must solve the program P µ max through Equations (12)- (17), but fixing the variables x b to the values obtained according to the predefined final pit. Therefore, the model P µ max only solves for the variable z r . With this, CVaR and expected profit values can be computed for each final pit. In the cases of (iii) Best-simulation and (iv) Hybrid-pit, 50 different results were obtained, but the chosen ones are those with minimum values according to criterion (C1) after applying model P µ max . Figure 9 shows the risk and profit for all methodologies and the reference ideal point. As expected, the final pit from [31] Table 2 shows a summary of the comparison: CVaR, expected profit and distance to ideal point (DIP, criterion C1), and their respective relative variations (RV, in %) for each alternative compared to the stochastic final pit obtained with µ = 125 [MUSD], which minimizes DIP, as shown in Table 1. In this case study, there is up to 32% improvement from our stochastic final pit compared to the deterministic one.
Thus, the stochastic final pits in the efficient frontier offer the best trade-off between lower risk in terms of CVaR and higher expected profit when compared to other alternatives, especially the deterministic one. Therefore, applying the proposed model based on stochastic programming reduces the risk and maximizes the expected profit. However, the quantification of this better performance depends on the grade distribution, the representation of geological uncertainty, and the parameters used in each case study. The efficient frontier's knowledge using geological scenarios allows for making better decisions in strategic open-pit mine production planning regarding expected profit, risk, or both.

Conclusions and Future Work
This paper develops a stochastic model that allows generating the efficient frontier of final pit limit alternatives under geological uncertainty in the expected profit vs. risk context. The variability of the deposit is modelled using several conditional simulations, and the risk of losses is measured in terms of conditional value at risk. The proposed methodology's main advantage is the generation of an optimal policy to manage the tradeoff between expected profit and risk, represented by the efficient frontier, allowing the evaluator to make the best decision according to the mining company's interests and its aversion to risk.
To test the performance of the proposed methodology to determine final pits, the efficient frontier obtained was compared with different approaches available in the literature. The criterion used was to minimize the distance to the ideal point (DIP, the unfeasible final pit with maximum expected profit and minimum risk along the efficient frontier). Our approach shows better results in controlling the risk of suffering economic losses without renouncing high expected profit.
Future research work includes developing new efficient algorithms for finding optimal solutions, because computing time will be a limitation when finding solutions of P µ for a model including millions of blocks and hundreds of conditional simulations representing the geological uncertainty. As an alternative approach, instead of computing the entire efficient frontier of final pits, one may optimize the criterion directly (a priori methods) and compute one optimal solution according to the decision-maker preferences [47,48].
Finally, it would be interesting to explore other criteria to compute optimal solutions or select them from the efficient frontier. This paper proposed some guidelines, and the stochastic dominance ideas developed by [36] are another interesting starting point to continue researching on this topic.

Data Availability Statement:
The data that support the findings of this study are available on reasonable request from the corresponding author. The data are not publicly available due to confidentiality issues.