Tri-Criterion Model for Constructing Low-Carbon Mutual Fund Portfolios: A Preference-Based Multi-Objective Genetic Algorithm Approach

Sustainable finance, which integrates environmental, social and governance criteria on financial decisions rests on the fact that money should be used for good purposes. Thus, the financial sector is also expected to play a more important role to decarbonise the global economy. To align financial flows with a pathway towards a low-carbon economy, investors should be able to integrate into their financial decisions additional criteria beyond return and risk to manage climate risk. We propose a tri-criterion portfolio selection model to extend the classical Markowitz’s mean-variance approach to include investor’s preferences on the portfolio carbon risk exposure as an additional criterion. To approximate the 3D Pareto front we apply an efficient multi-objective genetic algorithm called ev-MOGA which is based on the concept of ε-dominance. Furthermore, we introduce a-posteriori approach to incorporate the investor’s preferences into the solution process regarding their climate-change related preferences measured by the carbon risk exposure and their loss-adverse attitude. We test the performance of the proposed algorithm in a cross-section of European socially responsible investments open-end funds to assess the extent to which climate-related risk could be embedded in the portfolio according to the investor’s preferences.


Introduction
Climate change will pose a challenge for the financial sector seeking a balance between purely financial goals-looking for high returns-and sustainability making a positive impact on the environment and society. Since 2015, by adopting the Paris Agreement on climate-change and the UN 2030 for Sustainable Development, there has been a clear commitment, especially in the European Union, to align financial flows with a pathway towards a low-carbon, more resource-efficient and sustainable economy. In 2018, the EU launched and Action Plan to set out a strategy for sustainable finance, that is, "the process of taking due account of environmental, social and governance (ESG) considerations in investment decision-making, leading to increased investments in longer-term and sustainable activities" [1]. As stated in this report, to date, environmental and climate risks had not been appropriately considered by the financial sector, which is why if the EU wants to reorient capital flows to a more sustainable economy, environmental and social goals will have to be included in the financial decision-making. To this end, the Markets in Financial Instruments Directive (MIFID II) and the Insurance Distribution Directive (IDD) provides that investment firms and insurance distributors should ask their clients' investments objectives as regard sustainability and take their preferences into account when providing financial advice. This implies that investors should be able to integrate into their financial decisions additional criteria beyond return and risk and then to extend the classical bi-criterion portfolio selection problem based on Markowitz's mean-variance approach [2] by adding one more criterion.
In the literature, tri-criterion portfolio selection problems have been addressed by several authors making use of multicriteria decision-making problems (MCDM). While, the first attempts to deal with three-criterion portfolio selection problems were made by means of exact methods, several heuristic approaches such as multi-objective evolutionary algorithms (MOEAs) have been employed to solve multi-objective portfolio optimization problems. A better understanding of recent MCDM approaches to deal with tri-criterion portfolio selection problems by reviewing the literature on exact methods and multi-objective genetic algorithms techniques is undertaken in Section 2. From this review, we find there is a clear gap in the literature in the use of MOEAs to asses and manage the impact of carbon risk related preferences in the portfolio selection problem.
In line with the goals of the Paris Agreement, the financial flows should be consistent with a pathway towards low greenhouse gas emissions. In recent years, in response to an increasing climate-conscious financial products demand, Morningstar, the most important information provider in the mutual fund industry, introduced the low carbon designation (LCD) eco-label [3]. This new mutual fund eco-label helps investors to easily recognize which mutual funds are aligned with the transition to a low-carbon economy [4]. The LCD is composed of two indices, the Carbon Risk score and the Fossil Fuel involvement. In our research, we only consider the fund-level Carbon Risk score, which is obtained by weighting the firm-level exposure and management of material carbon issues. As Krueger et al. [5] highlighted, institutional investors increasingly address climate-change related risks and they are also viewed as catalysing driving firms to meet the reduced emission target. Thus, the mutual fund industry, and in particular institutional investors is an ideal setting to test our proposal.
The main contributions of this paper are summarized as follows: (i) We propose an approach based on a recent multi-objective genetic algorithm called ev-MOGA [6] to assess and manage the impact of climate-change related risk including as a third objective a new measure of carbon risk based on the LCD benchmark provided by Morningstar. The proposed ev-MOGA ensure the convergence of uniformly distributed solutions with no conditions related to the type of the objective function.
(ii) Once we obtain the approximated 3D non-dominated surface for the tri-criterion multi-objective problem, we introduce a-posteriori preference handling approach including the investor's preferences about carbon risk and risk aversion to come up with a single solution. Thus, according to the investor's preferences, we provide a bounded ε-Pareto front with information about minimum carbon risk, minimum risk aversion, and maximum return. With this tool, the investors could address more effectively with tri-criterion problems in the sense that they can easily visualize and find a desirable trade-off among conflicting objectives.
The rest of the paper is structured as follows. A review of the literature about the tri-criterion portfolio selection problem addressed either by exact or heuristic methodologies is presented in Section 2. In Section 3, the proposed tri-criterion genetic multi-objective evolutionary algorithm for constructing low carbon portfolios is formulated. In addition, we provide a-posteriori approach to integrate carbon risk related preferences in financial decision-making. In Section 4, we present the numerical results obtained by the application of the ev-MOGA using different investor profiles for a data set of European socially responsible investments (SRI) open-end funds. The discussion of the results is presented in Section 5. Finally, Section 6 concludes the paper.

Literature Review
The idea of determining the Pareto efficient frontier in portfolio selection from a mean-variance (M-V) optimization was initially conceived by Markowitz [2]. The essence of the M-V model is that risk is the investor's main concern and he/she tries to minimize risk for the desired level of expected returns. Over the years, the Markowitz's model has been extended either through more complex risk measures or through additional constraints, and in recent years through the possibility to include additional objectives. In this context, two main approaches to deal with the extended portfolio optimization problem can be found: (i) Exact methods or (ii) Heuristic methodologies. In what follows, we provide an integrated review of the theoretical and practical works that have used exact methods or heuristic approaches for extended M-V portfolio selection.

M-V Extended Approaches by Exact Methods
Since the early 1970s, several authors have attempted to expand the classical bi-criteria portfolio selection model beyond the expected return and variance with exact methods. Three main groups of studies can be identified by dealing with this problem. The first group of authors have expanded the Markowitz's model by introducing additional constraints such as cardinality, round lots or buy-in threshold [7][8][9][10]. Alternative risk measures such as down-side risk measures or CVaR have been proposed in the second group of studies [11][12][13]. A literature review on risk measures in terms of computational comparison is conducted in Mansini et al. [14].
Not until the 20th century was the idea of additional objectives was further boosted by the third group of studies. A tri-criterion non dominated surface can be found in [15][16][17] using a constrained linear program (QCLP) approach by solving a quad-lin-lin optimization problem where the third objective is linear. By defining several measures of liquidity in Lo et al. [18] a three-dimensional mean-variance-liquidity frontier is constructed. A general framework for computing the non-dominated surface in tri-criterion portfolio selection that extends the Markowitz's portfolio selection approach to an additional linear criterion (dividends, liquidity or sustainability) is addressed in Hirschberger et al. [15]. By solving a quad-lin-lin program, they provide an exact method for computing the non-dominated surface that can outperform standard portfolio strategies for multicriteria decision makers. An empirical application where the third criterion is sustainability is developed to illustrate how to compose the non-dominated surface.
In Utz et al. [16] sustainability is included as the third criterion to obtain the variance-expected return-sustainability efficient frontier in order to explain how the sustainable mutual fund industry can increase its levels of sustainability. The tri-criterion non-dominated surface is computed through the Quadratic Constrained Linear Program (QCLP) approach, and from the experimental results, it can be concluded that there was room to expand the sustainability levels without hampering the levels of risk and return.
However, the existing proposals based on exact procedures to solve tri-criteria portfolio selection problems have limited capabilities when the third objective is non-linear. In such cases, heuristic techniques have been recently applied to solve multi-objective problems and to provide fair approximations of the Pareto front.

Moeas and the Extended M-V Portfolio Optimization Problem
The increasing complexity of financial decision-making problems has led researchers to apply heuristic procedures inspired by biological processes such as multi-objective evolutionary algorithms (MOEAs). Suggested at the beginning of the 1990s, MOEAs have been applied in several fields including finance, and in particular to solve the portfolio selection problem [19,20]. These techniques provide satisfactory approximations of the efficient frontier even when the problem involves non-convexity, discontinuity or non-integer variables. In Arnone et al. [21], a MOEA was proposed for the first time for optimal portfolio selection by using lower partial moments as a measure of risk. The first attempts to propose MOEAs as an extension of the M-V model aimed at considering additional constraints such as, cardinality, lower and upper bounds, transaction costs, transaction round lots, non-negativity constraints or sector capitalization constraints [10,[22][23][24][25][26][27][28][29][30]. A review of the state of the art of MOEAs in portfolio selection can be found in Metaxiotis and Liagkouras [31]. In Silva et al. [32] a unified method with the ability of tackling several constraints (floor and ceiling, cardinality, round-lot and pre-assignment) is presented.
Another group of researchers have also tried to propose alternative risk measures to variance, the most popular being: semi-variance, value at risk (VaR) and conditional value at risk (CVaR), the lower partial moments (LPM), the Expected Shortfall, the Skewness, and Risk parity [33][34][35][36][37]. Other approaches used MOEAs with copula functions to consider the dependence structure of asset returns, see for example Babaei et al. [38]. Some scholars have been focused on preference-based algorithms for portfolio selection to direct the search of the non-dominated solutions towards de region of interest [39].
Regarding the number of objectives, while the two-objective case is the most widely used among the authors, the tri-objective problem has risen in popularity in the last few years. A tri-objective optimization problem is proposed in Anagnostopoulos and Mamanis [40] to find the trade-off between risk, return and the number of securities in the portfolio. In this paper, the authors compare three evolutionary multi-objective optimization techniques for finding the best trade-off between risk, return and the cardinality of the portfolio. A multi-swarm multi-objective optimizer based on p-optimality criteria called p-MSMOEAs is proposed with three-objectives in Hu et al. [41]. In Rangel-González et al. [42], a fuzzy multi-objective particle swarm optimization (FOMOPSO) algorithm that implements a fuzzy controller is implemented to solve a three-objective portfolio optimization problem.
A recent approach, based on the elitist evolutionary multi-objective optimization algorithm called ev-MOGA [6], is adopted in Garcia-Bernabeu et al. [43] to derive the non-dominated mean-variance-sustainability surface. While in the previous study the third objective was an overall score of sustainability including an integrated measure of ESG issues, we now focus on climate-change related objectives such as to reduce low carbon emissions. To the best of our knowledge, this is the first time that a MOEA has been used to assess and manage the impact of climate-related risks in the portfolio selection problem as a response to the upcoming more stringent climate policies.

The Tri-Criterion Multi-Objective Approach by Ev-Moga to Manage Carbon Risk Exposure
During the last two decades, MOEAs for portfolio management have attracted scholars and practitioners attention as stated in Section 2.2. Next, some previous notions on multi-objective optimization and genetic multi-objective optimization techniques are provided.

Background on Multi-Objective Optimization and Ev-Moga
Multi-objective optimization is an important subclass of multiple criteria decision making techniques involving more than one objective function to be optimized simultaneously. Since the conflict degree between the objectives makes it impossible to find a feasible solution that simultaneously optimizes all the objective functions, there is a set of Pareto optimal solutions denoted as Pareto front at which none of the objectives can be improved without deteriorating at least one of the others. In general, a multi-objective optimization problem is stated as follows: where the vector w " rω 1 , ω 2 , . . . , ω n s T is a n-parameter set included in the decision space S, and f i pwq : R n Ñ R, i " 1, . . . , m, are the objectives to be minimized at the same time.
In recent years MOEAs have been widely accepted as useful tools for solving real-world multi-objective problems. Within MOEAs several powerful stochastic search techniques that mimic Darwinian principles of natural selection are included. In this study, we focus on the ev-MOGA algorithm proposed in Herrero [6], which combines the concept of Pareto optimality and ε-dominance due to Laumanns et al. [44], thus providing an approximated ε-Pareto set. Definition 1. Dominance: Let w 1 , w 2 P R n be two feasible solutions, and let fpw 1 q, fpw 2 q P R m be their image solutions in the objective space. Then, assuming that the objective functions have to be minimized, w 1 is said to dominate w 2 , denoted as fpw 1 q ă fpw 2 q, iff: Definition 2. Pareto front: Let Ω Ď R n be a set of vectors of feasible solutions with fpΩq as their image solutions. Then the Pareto front fpΩ P q is defined as follows: fpΩ P q contains all vectors fpw u q P fpΩq that are not dominated by any vector fpw v q P fpΩq, i.e., From this definition, it is easily to deduce that any vector fpw v q P fpΩq z fpΩ P q is dominated by at least one vector fpw u q P fpΩ P q, i.e., @ fpw v q P fpΩq z fpΩ P q : D fpw u q P fpΩ P q | fpw u q ă fpw v q Definition 3. ε-dominance: Let w 1 , w 2 P R n be two feasible solutions, an let fpw 1 q, fpw 2 q P R m be their image solutions in the objective space. Then, assuming that the objective functions have to be minimized, w 1 is said to ε-dominate w 2 for some ε ą 0, denoted as fpw 1 q ă ε fpw 2 q, iff: Definition 4. ε-approximate Pareto front: Let Ω Ď R n be a set of feasible solution vectors with fpΩq as their image solutions. Then, fpΩP q is called a ε-approximate Pareto front if any vector fpw v q P fpΩq z fpΩP q is ε-dominated by at least one vector fpw u q P fpΩP q , i.e., The set of all ε-approximate Pareto fronts of fpΩq is denoted as the ε-Pareto front fpΩ P q.
The most outstanding feature of this algorithm is that the optimal solutions are distributed uniformly across the ε-Pareto front. To this end, the ε-Pareto front is split into a fixed number of boxes forming a grid, so that the algorithm ensures that just one solution is stored by one box. The size of the boxes is determined by the value of ε i , which is calculated as follows: where, fi and f i˚c orrespond to the maximum and minimum value of the objective function f i , and n box is the number of boxes. In addition, ev-MOGA is able to adjust the width of ε i dynamically and prevent solutions belonging to the extremes of the front from being lost. For solving the ev-MOGA, the main population P k whose size is N indP explores the searching space S defined by the multi-objective problem during a number k of iterations. In the archive population A k the ε i -nondominated solutions are stored. Then, at the end of the iteration process, A k is an ε-approximate Pareto front fpΩP q, with pn box q m´1 as the maximum number of solutions, where m is the number of objectives. Furthermore, in the case that more than one ε-nondominated solution is detected, the solution that prevails in A k will be the one that is closest to the center of the box. Next, the new individuals obtained by crossover or mutation with probability of crossing/mutation P c{m are included in the auxiliary population GA k . Before running the algorithm, the following parameters should be defined by the analyst: • N indP " Size of the main population. • N indGA " Size of the auxiliary population. • k max " Maximum algorithm iterations. • P c{m " Probability of crossing/mutation. • n box " Number of boxes.
The main advantage of ev-MOGA is that they generate good approximations of a well-distributed Pareto front in a single run and within limited computational time. The original ev-MOGA algorithm is available at Matlab Central [45]: ev-MOGA in Matlab Central.

The Ev-Moga Tri-Criterion Portfolio Selection
In this study, beyond risk and return, we wish to consider an additional objective that minimizes the carbon risk exposure of a portfolio. Then, by introducing a third objective into the portfolio optimization model the efficient frontier becomes a surface in the three-dimensional space. The tri-criterion portfolio selection problem where the objectives are the risk of the portfolio, the returns, and the portfolio carbon risk exposure can be mathematically formulated as follows: where N denotes the available assets, µ i is the expected return of asset i (i " 1, 2, . . . , N), σ ij is the covariance between asset i and j. In addition, c i is the carbon risk score and ω i denotes the proportion of asset i in the portfolio. Algorithm 1 details the ev-MOGA steps to obtain the well-distributed ε-approximate Pareto front.
2: Initialize the population of candidate solutions P 0 and set A 0 " H 3: Conduct the multi-objective evaluation of portfolios from P 0 using Equations (7)- (10) 4: Detect the ε-nondominated portfolios from P 0 and store in the archive A 0 5: while k ď k max do 6: Generate the auxiliary population GA k from the main population P k and the archive population A k following this procedure: 7: for j Ð 1, N indGA {2 do 8: Randomly select two portfolios X P and X A from P k and A k , respectively 9: Generate a random number u P r0, 1s 10: If u ą P c{m , X P and X A are crossed over by means of the extended linear recombination technique, generating two new portfolios for GA k

11:
If u ď P c{m , X P and X A are mutated using random mutation with Gaussian distribution and then included in GA k

12:
end for 13: Evaluate population GA k using the tri-criterion multi-objective portfolio model defined by Equations (7)-(10). 14: Check which portfolios in GA k must be included in A k`1 on the basis of their location in the objective space. A k`1 will contain all the portfolios from A k that are not ε-dominated by elements of GA k , and all the portfolios from GA k which are not ε-dominated by elements of A k

15:
Update population P k`1 with portfolios from GA k . Every portfolio X GA from GA k is compared with a portfolio X P that is randomly selected from the portfolios in P k . X GA will replace X P in P k`1 if it dominates X P . Otherwise, X P will not be replaced 16: k Ð k`1 17: end while

Defining A-Posteriori Preferences for Each Investor Profile
With the previous multi-objective optimization design a vast region of the tri-objective whole Pareto front is generated. Even though it is true that the non-dominated surface allows us to better understand the trade-off between the three objectives, this solution does not provide a useful tool from the user's perspective. To come up with a single solution, we assume that the decision-maker is available to take part in the solution process. According to [46,47] the articulation of preferences may be done either before (a-priori), during (progressive), or after (a-posteriori) the optimization process. In what follows, we assume that once the investor has seen an overview of the Pareto optimal solutions, they takes part in the final solution. Thus, we propose an a-posteriori approach.
The analyst supporting a-posteriori methodology has to inform the decision-maker either providing a list of solutions or to provide a visualization of the Pareto front [48]. In a tri-objective case, two main approaches have been used to visualize the Pareto frontier: (i) three-dimensional graph, and (ii) decision maps. However, a new graphical visualization called Level Diagram is proposed in Blasco et al. [49] to represent n-dimensional Pareto fronts. The Level Diagrams tool also allows the incorporation of decision-maker preferences, and it offers an excellent tool to help in the decision-making process.
In our proposal, information on carbon risk and risk aversion related preferences is given by the investor, which determine the upper bounds for the ε-Pareto front. Let us denote the vector for the preferences about green investments defined by the carbon risk score objective function in Equation (9) as P g , and the preferences for the loss aversion attitude defined by Equation (7) as P r .
Concerning the carbon risk related preferences, we consider three types of green investor profile. They are defined as follows: 1.
Strong green investor. This profile is defined by low level of allowed carbon risk p s g .

2.
Moderate green investor. This profile is defined by medium level of allowed carbon risk p m g .

3.
Weak green investor. This profile is defined by high level of allowed carbon risk p w g .
Thus, the reference vector for the green investor could be stated as follows: Concerning the investor's loss aversion attitude, we consider three types of risk aversion investor profile. They are defined as follows: 1.
Conservative investor. This profile is characterized by investing in lower-risk securities, namely, a high loss aversion attitude p c r . 2.
Cautious investor. This profile is defined by medium risk tolerance, and consequently a moderate loss aversion attitude p k r . 3.
Aggressive investor. It includes investors that actively seek stocks with higher risk, but a chance for higher reward, that is a low loss aversion attitude p a r .
Thus, the reference vector regarding the risk aversion could be stated as follows: " p c r , p k r , p a r ı (12) Figure 1a shows a three-dimensional representation of the ε-approximate Pareto front generated by ev-MOGA for the tri-criterion problem defined in Equations (7)- (10). The point's color is faded from gray to blue: the more blue the color is, the more close to the global optimal point we are, which is colored in red. The optimal point is the one with the minimum normalized distance to the ideal. As may be seen, ev-MOGA provides uniformly distributed efficient portfolios along the ε-Pareto front.
The impact of considering the investor's preferences (11) and (12) into the decision-making process is examined in Figure 1b. The new colored area, faded from blue to green, defines the investor's region of interest throughout a bounded ε-Pareto front according to the previously defined preferences. Furthermore, in red we have colored the portfolios that can be proposed to the investor to maximize the returns (max ret), minimize the risk (min risk) or minimize the emission risk (min carb). The red point located in the green area is associated with the optimal portfolio that simultaneously achieves the three objectives within the investor's region of interest.

Empirical Application
We use a set of monthly returns on 22 institutional SRI European open-end funds offered in Spain for the period 2009-2019. The empirical information includes the time series of 120 monthly returns and the carbon risk indices. As a previous step the expected return vector µ " pµ 1 , . . . , µ 22 q T and the covariance matrix Σ " rσ ij s, i, j " 1, . . . , 22 are computed. For the carbon risk score c i , we use the Morningstar R Portfolio Carbon Risk Score, which indicates the risk that companies face from the transition to a low-carbon economy. In this set, scores c i range from 0 to 10, where lower scores are better, indicating lower carbon risk levels. All the numerical information to be used on this opportunity set comes from Morningstar database. Table 1 shows the parameter setting applied to the ev-MOGA algorithm. The size of the main population is N indP " 10 4 , while the population of the archive A k is N indGA " 500. For the probability of crossing/mutation we select P c{m " 0.2. Finally, the space of each objective function has been divided in 300 boxes. Table 1. Parameter setting of the ev-MOGA.

Parameter Value
Size of the main population N indP " 10 4 Size of the auxiliary population N indGA " 500 Maximum algorithm iterations k max " 10 5 Probability of crossing/mutation P c{m " 0.20

Number of boxes n box " 300
Maximum portfolio weights ω max " 0.20 The proposed approach was coded in Matlab R R2020a, and it was executed on an Intel i7 3.20 GHz, 20 GB (2667 MHz DDR4) of RAM. As a result, the ev-MOGA algorithm provided an ε-approximate Pareto front of N p " 7361 solutions stored in two data arrays: first, an pN pˆ3 q array containing the values of the three objectives; second, an pN pˆ2 2q array with the corresponding portfolio weights, which have been rounded to three decimal places. The computation of the ε-approximate Pareto front required a time of 34 h.
Next, in the a-posteriori approach we obtain the bounded ε-approximate Pareto front according to the following carbon risk and risk aversion related preferences:

•
Concerning the carbon risk related preferences, the green investors are classified into three profiles according to expression (11) by using the 25th percentile for the Strong green investor p s g " P 25 , the 55th percentile for the Moderate green investor p m g " P 55 , and the 75th percentile for the Weak green investor p w g " P 75 . Thus, the reference vector for the green investor is obtained by applying these percentiles to the column of Carbon risk at the pN pˆ3 q array containing the ε-approximate Pareto front: P g " rP 25 , P 55 , P 75 s " r2.803, 3.001, 3.135s • Considering the investor's loss aversion attitude, the investors are classified into three profiles according to expression (12) by establishing the 50th percentile for a Conservative investor p c r " P 50 , the 75th percentile for a cautious investor p k r " P 75 , and the 100th percentile for an Aggressive investor p a r " P 100 . Thus, the reference vector regarding the risk aversion is obtained by applying these percentiles to the column of Risk at the pN pˆ3 q array containing the ε-approximate Pareto front: P r " rP 50 , P 75 , P 100 s " r9.575, 10.095, 11.633s From Tables 2-4, a comparison of efficient portfolios for Weak, Moderate and Strong green investors is made in terms of different loss aversion attitude. To this end, for each profile, we display a numerical description of the portfolio composition (F i " 100¨ω i ) and the objective function values attained by the portfolios. Figures 2-4 show the 3D representation of the ε-approximate Pareto front, thus providing the non-dominated mean-variance-emission surface for the three types of Green investor profile and each level of loss aversion. Notice that, as the level of loss aversion attitude decreases, the Green Investor non-dominated surface (colored in blue) grows.

Discussion
The results are displayed in Tables 2-4, for each green investor profile. Each table contains three groups of rows corresponding to the risk aversion profile, and two groups of columns with the information about the funds' portfolio composition (F i " 100¨ω i ) and the objective function values attained by the portfolios [Risk, Ret., Carb.]. Given a green investor profile and a risk aversion profile, we consider four strategies as follows: 1.
In the first row, we highlight in bold optimal funds allocation for the nearest solution in the ε-approximate Pareto front to the ideal values of return, risk, and carbon risk, when the investor is seeking the optimization of the three objectives simultaneously [opt].

2.
In the second row, we display the optimal funds' allocation when the investor is focused on achieving the minimum risk [min risk].

3.
In the third row, we report the optimal funds' allocation when the investor tries to minimize the carbon risk exposure of the portfolio [min carb].

4.
In the fourth row, we display the optimal funds allocation when the investor seeks to maximize the return of the portfolio [max ret].
As a complement to the information supplied by above mentioned tables, we add Figures 2-4. Thus, each figure is associate to a green investor profile and includes three representations of the ε-approximate Pareto front for the corresponding risk aversion profile: (a) conservative, (b) cautious, and (c) aggressive.
Let us see, for example, the case of a Weak green investor with a conservative attitude toward risk, which information is provided in Table 2 and Figure 2a. If the investor is seeking the optimization of the three objectives simultaneously, the optimal portfolio allocation is given by F 3 , F 10 , F 11 , F 12 , F 14 , F 16 , and F 21 [opt]. We can also view the 3D non-dominated surface, Figure 2a, in which the whole ε-approximate Pareto front is colored in grey, and the bounded investor's region of interest is faded from blue to green. While the nearest solution to the ideal values for the three objectives lies at the center of the bounded region of interest [opt], the corner solutions indicate the objective values involving minimum risk [min risk], minimum carbon risk [min carb], and maximum return [max ret]. A red dot marks these optimal values. Note that the bounded region of interest increases as the investor's risk aversion decreases (see Figure 2b,c). We can observe that the more aggressive attitude toward risk, the more returns and the lower carbon risk values are reached.
When analyzing the Moderate and Strong green investor results, displayed in Table 3- Figure 3 and Table 4- Figure 4, respectively, we would like to remark the same behavior as explained for the Weak green investor. Notice that, as the level of loss aversion attitude decreases, the bounded ε-approximate Pareto front (colored in blue-green) grows.
Comparing data through Tables 2-4, we stress at this point that the lower level of returns corresponds to a Strong green investor regardless of the investor's risk aversion preferences, as shown in Figure 4. As expected, the minimum region of interest corresponds to a Strong green investor with a conservative attitude toward risk, as shown in Figure 4a. We also would like to remark that a Moderate green investor with an Aggressive attitude toward risk can achieve a high level of return (see Figure 3c). Thus, we conclude that green investors have a leeway to decrease the emission risk of the portfolio at no cost to risk and returns.

Conclusions
In this paper, we have proposed a new application of the ev-MOGA algorithm to handle a tri-criterion portfolio optimization problem in which the third criterion is the carbon risk score of the portfolio. Moreover, we have incorporated the investor's preferences regarding the risk emissions and the loss aversion attitude into the solution process by defining different investor profiles. This allows us to propose a solution to the investor in terms of their carbon risk and risk aversion related preferences.
Given the urgency around climate change, investors are becoming increasingly aware of the need to make the transition to a low carbon economy and to address climate-change related risks. To the best of our knowledge, this is the first time that a MOEA based on the concept of ε-dominance called ev-MOGA has been used to assess and manage the impact of climate-related risks in the portfolio selection problem. The ev-MOGA allows us to derive a 3D-Pareto front in a well-distributed manner with limited memory resources. Next, we have introduced a-posteriori approach to include the investors' preferences about green investments and risk aversion to understand the trade-off between the three objectives. So, by considering different investor profiles, we can provide a more approximated solution according to their preferences. In applying this approach to a sample of institutional SRI European open-end funds, we obtain a global ε-approximate Pareto front including a vast region of solutions. The a-posteriori analysis shows the bounded ε-approximate Pareto front for each investor profile and, we find that the more aggressive attitude towards risk, the more returns and the lower the carbon risk values, namely, aggressive investors looking for high returns are allowed to invest in funds with a lower level of carbon risk scores. Moreover, we can also conclude that green investors have a leeway to decrease the carbon risk of the portfolio at no cost to risk and returns.
Because of the possibility to obtain an efficient frontier in three dimensions while including the choices on risk and climate-change related risks, we believe this is a useful tool for investors, especially for those who are willing to rebalance their portfolios towards more climate-conscious firms.