Eco-Driving in Railway Lines Considering the Uncertainty Associated with Climatological Conditions

: Eco-driving is a keystone in energy reduction in railways and a fundamental tool to contribute to the Sustainable Development Goals in the transport sector. However, its results in real applications are subject to uncertainties such as climatological factors that are not considered in the train driving optimisation. This paper aims to develop an eco-driving model to design efﬁcient driving commands considering the uncertainty of climatological conditions. This uncertainty in temperature, pressure, and wind is modelled by means of fuzzy numbers, and the optimisation problem is solved using a Genetic Algorithm with fuzzy parameters making use of an accurate railway simulator. It has been applied to a realistic Spanish high-speed railway scenario, proving the importance of considering the uncertainty of climatological parameters to adapt driving commands to them. The results obtained show that the energy savings expected without considering climatological factors account for 29.76%, but if they are considered, savings can rise up to 34.7% in summer conditions. With the proposed model, a variation in energy of 5.32% is obtained when summer and winter scenarios are compared while punctuality constraints are fulﬁled. In conclusion, the model allows the operator to estimate better energy by obtaining optimised driving adapted to the climate.


Introduction
The current global context of sustainable development urges all countries to mitigate the consequences of climate change. It is also a time when the current trend of societies in developed countries is to concentrate the population in huge metropolia. The railway sector plays a vital role in transporting people between these large metropolitan areas. Its immediate future lies in making the train more attractive to potential users. In this regard, high-speed railways (HSR) has been an advance, providing a fast and comfortable transport mode for long-distance services. HSR is considered an efficient transport mode and is part of the strategy to avoid the pollution produced by the massive use of private vehicles. Europe is leading among the developed countries in managing decarbonisation in all sectors with imminent effect. In this way, the Sustainable Development Goals (SDGs) are the main framework in which measures have to be taken. Therefore, the most immediate action to carry through by society, administrations, and companies is to develop technology by making it more efficient in terms of emissions of greenhouse gases.
This rational use of energy is particularly relevant for the transport sector because of the massive amount of energy required in developed countries [1,2], especially when the trend is towards unification into an electrified system. In this manner, railway transport is the most efficient in energy consumption and pollution emission. Railway operations account for a low percentage of energy consumption, but, even so, it is a large consumer and, therefore, its efficiency is susceptible to being improved [3]. Railway operators, manual driving. Direct search algorithms have been applied to minimise energy traction in [47] while Brute Force algorithm is used in [48] to minimise energy consumption in substations. A Monte Carlo Simulation model was proposed to determine the system energy flow considering regenerative braking trains [49]. Differential Evolution has been presented in [50] to optimise railway operation. Optimisation by minimising energy traction is carried out by Simulated Annealing in [51] and searching optimal coasting points in [52]. On the other hand, multi-objective models have been developed using numerical methods. Among them, Indicator Based Evolutionary Algorithm (IBEA) was proposed to optimise both the running time and energy consumption [53]. Non-dominated Sorting Genetic Algorithm II (NSGA-II) is used in [54] to design ATO CBTC speed profiles to minimise energy traction and in [55], its performance is compared with Multi-Objective Particle Swarm Optimisation (MOPSO) generating the Pareto curve for real cases of the Madrid Underground. Ant Colony Optimisation can be seen in [27], proposing a distance-based train trajectory optimising the train trajectory, and in [56] to achieve the optimal driving strategy. Finally, a MOPSO algorithm is used in [57] to design efficient speed profiles for the ATO of the trains of a metro line.
With the purpose to obtain accurate results, it is important not only to use detailed simulation models but also to take into account the possible uncertainties that emerge from the HSR operation. These uncertainties can be classified in three groups: limitations in the movement authority produced by signals or preceding trains [58,59], variability in the manual operation of the train [45], and the climate conditions. While the first two uncertainty sources have been widely studied, climate conditions have not been deeply considered in the eco-driving problem.
Most academic works focus on winter scenarios in Northern Europe to forecast delays or disruptions in traffic operation. Some studies have reviewed the issues during winter railway operations caused by snow and ice in Nordic countries [60]. In study [61], impact of humidity, temperature, ice, and snow depth on the occurrence of delays in high-speed operation in Sweden has been analysed through a Cox model and Markov chain model. A train delay prediction model with weather data has been designed in [62] and tested on two high-speed railway lines in China. The Nordland railway line has been assessed in [63], where extreme cold weather is a crucial factor related to delays and low punctuality. Other papers have focused on service delay time and its exposure time to bad weather [64]. A combination of fuzzy theory and rough sets under adverse weather conditions has been used to contribute to the forewarning method for train operation in [65]. Causal relationships between extreme weather patterns and derailments on railway turnouts have been investigated by means of Fuzzy Bayesian Network model [66]. The impact on freight railways in winter weather in Northern Europe has been discussed [67]. In study [68], an estimation of the effects of weather conditions on railway operator performance of passenger train services has been carried out, focusing on infrastructure disruptions. In study [69], a model to capture the cascade dynamics of delay propagation on railway networks under inclement weather has been presented.
As a general rule, climatological factors are considered random variables and instantly excluded from the models considered by railway operators. However, in HSR, the running resistance is the major factor determining the trains' energy consumption. Climatic conditions such as wind or air density can produce significant variations in running resistance value. Therefore, climatology plays an important role in the energy that high-speed trains demands. Previous reasons motivate this paper to generate a new model to quantify the climatological affection in railway operation, both in energy consumption and in running time. Climatology has a great variability due to the impossibility of having complete knowledge of all the variables that model its behaviour. Therefore, climatic conditions such as pressure, temperature, and wind have been modelled by means of fuzzy logic [70]. This way, this affection is included in the running resistance, given the vagueness in the knowledge of these variables. The main advantage of fuzzy logic modelling is the ability to handle imprecise or incomplete information. Additionally, a positive feature is its flexibility and simplicity to be developed and implemented to obtain efficient calculations. That is why many studies have applied fuzzy logic in train operation models. In study [71], the train service provider and the infrastructure provider have been modelled as software agents, and the negotiation for a train schedule and the associated track access charge as a prioritised fuzzy constrain satisfaction problem. In study [72], predictive fuzzy control is used for carrying out the train dwell-time control in an event-driven framework. A methodology for predictive fuzzy control in an event-driven environment for multi-objective decision making on DC railway systems has been presented in [73]. In study [74], an approach within the fuzzy framework for tackling the complexity of multidimensional service evaluations and the judgment on the goodness of the schedule has been proposed. In study [75], the authors investigated a passenger train timetable problem with fuzzy passenger demand. The paper presented by [76] exposed a framework to evaluate the logistics performance of intermodal freight transportation by applying fuzzy set techniques. A multilevelled distributed railway traffic fuzzy control system is proposed in [77]. An assistant support system for railway traffic control is described in [78] making use of fuzzy knowledge. A fuzzy optimisation model for rescheduling high-speed timetables is proposed in [11], focusing this work on punctuality and passengers' affection without considering energy minimisation. In studies [79,80], a fuzzy model is proposed for an Automatic Train Control system. In studies [45,59], uncertainty in the manual driving application of the speed regulation command is modelled using fuzzy numbers. In study [16], the optimal speed profile in urban DC railways is found with a GA, where the energy consumption and the running time are the fuzzy results. A fuzzy model of delays and punctuality constraints is proposed in [46] for the offline design of timetables when efficient driving is applied. The algorithm proposed in [54] is used for urban ATO operation including uncertainty in the mass, generating optimal energy-time Pareto curves of speed profiles. Fuzzy logic is applied in [81] to detect the railway wheelset conicity as a measure of the deterioration. In study [82], a fuzzy model is used to evaluate the potential risks of railway crossings. As can be seen, climatological parameters and their associated uncertainty have not been modelled previously; therefore, this possible affection to the railway operation has not been considered.
The main contribution of this paper is the introduction of climatological parameters in the eco-driving optimisation model. The proposed model allows an efficient driving design considering the uncertainty due to climate conditions by means of fuzzy modelling. Climatological parameters are not typically considered in the design of speed profiles, and they could affect the running time and energy consumption.
The optimisation model is based on a Genetic Algorithm combined with fuzzy parameters. This algorithm applies a search process that optimises a driving holding speed without a braking commands set that can be easily applied by human drivers. Applying these commands generates the efficient driving with punctuality requirements at the destination. Using a detailed simulation platform makes solving all possible scenarios possible when the climatological parameters vary according to the proposed model. This methodology can also be applied to quantify the error obtained in the expected energy saving if climate conditions are not initially considered in the optimisation of the railway operation. The model has been applied in a simulation test based on real data of a Spanish HSR, to compare scenarios where different wind and air pressure conditions are considered.
The paper is organised as follows: Section 2 describes the simulation model in which the weather parameters are included. The uncertainty associated with the climatological parameters is modelled in Section 3. Section 4 is a description of the optimisation model developed to generate efficient driving. Section 5 analyses the case study and presents the results obtained for a real high-speed railway line using the models proposed in previous sections. A discussion of the proposed procedure and the results are presented in Section 6. At the end of the document, the bibliographical references and previous academic studies from which information for the present analysis has been obtained are presented.

Simulation Model
The simulation model is described in this section. This model is detailed in [45], and the results provided by the simulation were compared with real data of a Spanish HSR to validate the simulator. Differences of 1.2% for the running time and differences of 0.4% when compared with the energy consumption were shown. The model was used from [45] to design eco-driving in a high-speed line, and measurements were recorded. As a result, energy savings between 20% and 34% were measured in commercial services. It was a detailed discrete event deterministic simulation model, that combines next-event and fixed-discrete time steps for advancing the simulation clock [83,84]. The simulation model takes several input data related to the rolling stock involved in the operation and related to the railway line:

O
Rolling stock.
i. Physical parameters: mass of the train, the mass of the load, length, maximum speed, rotatory inertia, and adhesion traction. ii.
Traction effort curve, braking effort curve, running resistance coefficients, energetic and power systems and auxiliary systems consumptions.
O Railway line. Includes Kilometric Points of relevant information as stations (K.P.), grades, track curvatures, speed limits, tunnels, and grade transitions considering the length of the train and the electric neutral zones.
The equation that models the motion of the train at each simulation step considering the forces acting on the train is: where • a is the acceleration of the train. • F m is the motor force. • R is the running resistance of the train defined by the Davis formula (Equation (7)). • F g is the force due to the railway grades. • F r is the force due to track curvature. • m eq is the equivalent mass.
By definition, the acceleration of the train in every step of the simulation is: The equivalent mass, effective mass or inertial mass can be calculated by: where • m 0 is the empty train mass. • m l is the load mass. • λ is the dimensionless rotating mass factor.
Force due to gradient (Equation (4)) shows acceleration due to gravity effect, considering uphill or downhill motion. Therefore, the force due to gradient profile can be written as: where Sustainability 2022, 14, 8645 6 of 26 • g is the gravity acceleration. • θ is the slope angle. Given the small values of slope angle in typical railway tracks, Equation (4) can be written as: where • p is the gradient of the track (mm/m or ‰).
Likewise, the force due to the track curvature can be expressed as follows: where • K r is an empirical constant defined by track gauge. • r is the radius of track curvature.
The traction effort or motor force (F m ) is bounded by a maximum value of the electrical traction effort curve and a maximum electrical braking effort curve both dependent on the train speed. The simulation model considers at every moment the speed limits established along the rail track by the railway administration. Speed limits are defined by stretches along the railway where the train must drive respecting the maximum speed limit considering the total length of the train. When a train passes through a neutral zone, auxiliary systems cannot be fed from a catenary, and motors apply at least a constant braking effort to maintain the charge of the batteries that feed auxiliary systems along neutral zones.
Running resistance is the model that defines and groups a series of resistances opposing the train's movement. These include mechanical resistances (rolling and internal friction), air intake resistances (for engine cooling and passenger air renewal), and aerodynamic resistance [85]. The running resistance is modelled with a second-degree equation dependent on the instantaneous velocity and defined by the Davis expression [86]: where • a is a constant related to the mechanical resistance to the motion. • b is a constant related to the resistance due to the air inlet in the train. • c is a constant related to the aerodynamic resistance. • k is the tunnel factor. • v train is the velocity of the train.
The Davis formulae and similar expressions of a polynomial nature are models of reality. Each term of this equation refers to different physical components of the total running resistance. However, it must be considered that these polynomials are fitted models of highly complex phenomena that interact themselves. Therefore, although for general purposes certain approximations can be estimated, in reality, they are not completely accurate and do not represent reality in its entirety. Likewise, within the running resistance model, there are parameterised coefficients (a, b and c) that represent different oppositions to the movement of the train. These values are experimental, complex to define, and depend on the physical characteristics of the rolling stock. In order to obtain them, manufacturers test their trains and their relationships with the resistances generated. These coefficients are the nominal values, so they have been evaluated under standard conditions of pressure and temperature given by the International Standard Atmosphere, ISA (ISO, 1975: T = 15 • C; ρ = 1.225 kg/m 3 ).
The term of the running resistance model (Equation (7)) that is not related to the effect of the air outside the train is called mechanical resistance. In this term, the coefficient a is involved. In the most general case, it is derived from the frictional resistance between mechanical components, rail-to-wheel contact, track irregularities, and energy losses in the traction and suspension equipment of the vehicles due to the oscillating or parasitic movements of the suspended mass.
The term of the train speed-dependent running resistance corresponds, mainly, to the resistance caused by the air intake of the train. In trains, there is a constant and noticeable air flow, necessary for cooling engines and other equipment as well as that required for air renewal for passengers. This term is related to the b coefficient of the running resistance model.
The so-called aerodynamic running resistance is the longitudinal force that opposes the train's movement as a consequence of the interaction between the train and the surrounding air with which it collides and envelops it. The coefficient c is responsible for modelling this resistance.
The tunnel factor k increases the quadratic term of the running resistance model when the train runs through a tunnel due to the increased air friction against the outer surface of the train. When a tunnel is not involved in the simulation, the factor takes the value 1.
Equations (1)-(7) describe the motion of the train in next-event and fixed-discrete simulation steps, limited by speed limits, calculating braking curves to stop at the arrival point or station point with comfort and safety requirements.
The model related to the energy consumed by the train measured at the pantograph E pantograph is expressed as shown in Equation (8) while the train energy consumption estimated at the electrical substation (E substation ) is defined by means of Equation (9): where • P pantograph is the electrical power consumed by the train measured at the pantograph. • P substation is the estimation of the power consumed at the substation.
Previous calculus requires some other equations to complete the model: where • P mec is the mechanical power. • P aux is the power consumed by the auxiliary systems. • η T is the electrical chain efficiency for the traction case.

Climatological Model
The following section focuses on the running resistance variation when the railway operation is developed under real climatological circumstances. Initially, the influence of the variation of density along the railway line due to altitude changes has been modelled. Then, the wind's influence on the running resistance model is incorporated. When motion equations are applied (Equations (1)- (7)) for simulating the train's driving in a railway line, it is supposed that pressure and temperature values are constant along the whole line under standard conditions, given by the International Standard Atmosphere [87]. Thus, Equation (7) is just modified by the velocity of the train and it incorporates coefficients, b and c given by the manufacturer. For that reason, previous coefficients are also constant along the railway operation in absence of wind and they are always the same for each line, no matter the geographical location and the season of the year considered for the railway operation. The purpose of this section is to include in the simulation model the influences on the running resistance with climatological considerations of pressure, temperature, and wind.
The railway operation can occur between points of different altitudes, and it is, therefore, to be expected that the model must deal with this variation along the railway line. In summary, throughout the journey, the air density will vary due to the change in altitude and this variation will modify the running resistance equation according to the model described below. When the train runs on a line without wind, the running resistance term can be defined by Equation (7). In the running resistance model, lineal and quadratic terms vary as a function of train speed. In the quadratic term, the coefficient c can be defined by the general aerodynamics equation as shown below, according to the quasi-steady theory [88]: where • c is the aerodynamic coefficient of the running resistance model (Equation (7)). From Equation (15), it is deduced that factor c of the Davis formula is dependent on the density of the air and, therefore, on its variation. Finally, with the previous information, the coefficient c is parameterised as follows, varying with density and considering standard weather conditions as a reference [89]: Finally, when all information related is considered, the running resistance model can be written as follows. However, this equation models the situation where the train is running without being influenced by the wind: In this manner, a main and essential task is to know the air density at every point or analyse the stretch along the railway line. The density is a function of the temperature and pressure of the air, and its relationship can be expressed by the equation of state defined in the ideal gas law as follows: where • ρ is the density. • p is the pressure.
T is the temperature.
Pressure and temperature values strongly depend on the meteorological conditions given a specific area due to its geographical characteristics and the year's season. Therefore, considering the climatological model, they can be significantly different from one area to another. This work analyses these differences, which are present in the variations along the railway line.
On the other side, the running resistance force is also significantly affected by the wind and its variation along the route. The influence of the wind is modelled by the combination of its intensity, or wind speed, and the angle of incidence in each stretch relative to the train (φ) [88][89][90]. The train speed relative to the wind is the train speed minus the wind speed vectors, as it can be seen in Figure 1. Figure 2 shows the reference frame considered in this work where angle φ (Equation (19)) is the principal parameter to model the wind influence in the railway operation.
where • β and φ are the angles defined in Figure 1 by the velocities. • v wind is the velocity of the wind. • v rel is the relative train speed of the train, graphically defined in Figure 1.
is the temperature.
Pressure and temperature values strongly depend on the meteorological conditions given a specific area due to its geographical characteristics and the year's season. Therefore, considering the climatological model, they can be significantly different from one area to another. This work analyses these differences, which are present in the variations along the railway line.
On the other side, the running resistance force is also significantly affected by the wind and its variation along the route. The influence of the wind is modelled by the combination of its intensity, or wind speed, and the angle of incidence in each stretch relative to the train ( ) [88][89][90]. The train speed relative to the wind is the train speed minus the wind speed vectors, as it can be seen in Figure 1. Figure 2 shows the reference frame considered in this work where angle (Equation (19)) is the principal parameter to model the wind influence in the railway operation.
where  and are the angles defined in Figure 1 by the velocities.  is the velocity of the wind.  is the relative train speed of the train, graphically defined in Figure 1.
Running resistance model, given by Equation (7), is valid just when the wind has not been considered, so the model has to be modified, now influenced by the speed of the where  is the density.  is the pressure.  is the ideal gas constant.  is the temperature.
Pressure and temperature values strongly depend on the meteorological conditions given a specific area due to its geographical characteristics and the year's season. Therefore, considering the climatological model, they can be significantly different from one area to another. This work analyses these differences, which are present in the variations along the railway line.
On the other side, the running resistance force is also significantly affected by the wind and its variation along the route. The influence of the wind is modelled by the combination of its intensity, or wind speed, and the angle of incidence in each stretch relative to the train ( ) [88][89][90]. The train speed relative to the wind is the train speed minus the wind speed vectors, as it can be seen in Figure 1. Figure 2 shows the reference frame considered in this work where angle (Equation (19)) is the principal parameter to model the wind influence in the railway operation.
where  and are the angles defined in Figure 1 by the velocities.  is the velocity of the wind.  is the relative train speed of the train, graphically defined in Figure 1.
Running resistance model, given by Equation (7), is valid just when the wind has not been considered, so the model has to be modified, now influenced by the speed of the Running resistance model, given by Equation (7), is valid just when the wind has not been considered, so the model has to be modified, now influenced by the speed of the train and wind vectors, as shown in Figure 1. The modification consists of replacing the train speed in Equation (7) by the train speed relative to the wind, which is equivalent to considering zero wind in the original equation: Besides and in addition to the previous model, it has to be considered that factor c also depends on the angle β as shown in Figure 1. Yaw angle can be calculated using Equation (21).
To get the final model for the running resistance of a train, the relative train speed must be calculated, determining the simulation model as follows: Equation (23) is proposed in [91]. This expression establishes the variation of the constant c according to the angle β. This is an empirical expression and becomes meaningless when β > 30 • . An angle with such magnitude for typical speed values in operation would imply extreme winds, which would affect the train operation or could lead to the total suspension of the service.
where • c(ρ) is the value of c when density variation is considered according to Equation (16).
Finally, the variable running resistance considering climatological factors is modelled as: When the train runs through a tunnel, the force exerted by the influence of the wind shall not be considered. Therefore, for the sections simulated in a tunnel, the expression for the running resistance must be modelled with an expression that considers the train speed as input data, including the tunnel factor in the quadratic term. As applicable, it is used a k parameter as the tunnel factor: Pressure, temperature and wind models have been discretised along the railway line. If the altitude suffers any variation along the journey, the model can deal with the change by measuring the impact of the air density variance. Two different reasons have justified the division or discretisation along the railway line. Firstly, the variability of the wind makes necessary a partition along the line where the climatological and geographical characteristics determine the main direction in which the wind blows. And secondly, the change of directions taken by the train along the railway line has been dealt implementing stretches where it can be supposed that the relative angle between the wind and the train movement is approximately constant along a straight line. Hence, the proposed wind model considers both the railway track angle and the incidence wind angle to be parametrised. In conclusion, previous considerations allow the simulation model of the train movement to deal with the wind behaviour's complexity without compromising the final results' accuracy.

Fuzzy Climatological Model
This paper quantifies the relevance of the climatological factors and their variation by modelling them with their associated uncertainty during railway operation by means of a climate fuzzy model.
The parameters modelled by fuzzy numbers [92] are defined in Section 2.2: the pressure and temperature defining, consequently, the density, and the intensity and the angle of incidence of the wind.
Among the many possibilities for the membership function of the fuzzy numbers: piece-wise linear [93], linear, hyperbolic, exponential and S-shaped [94,95], climatological factors have been defined as triangular fuzzy numbers. But the model proposed in this paper would be applicable if any other shape is utilised to model climatological parameters' uncertainty. Each of the fuzzy numbers, the temperature, barometric pressure, air density, wind intensity (or speed) and the angle of incidence of wind are defined by the shape shown in Figure 3. Therefore, a fuzzy number is defined by its core, its lower limit and its upper limit of the support.
by modelling them with their associated uncertainty during railway operation by means of a climate fuzzy model.
The parameters modelled by fuzzy numbers [92] are defined in Section 2.2: the pressure and temperature defining, consequently, the density, and the intensity and the angle of incidence of the wind.
Among the many possibilities for the membership function of the fuzzy numbers: piece-wise linear [93], linear, hyperbolic, exponential and S-shaped [94,95], climatological factors have been defined as triangular fuzzy numbers. But the model proposed in this paper would be applicable if any other shape is utilised to model climatological parameters' uncertainty. Each of the fuzzy numbers, the temperature, barometric pressure, air density, wind intensity (or speed) and the angle of incidence of wind are defined by the shape shown in Figure 3. Therefore, a fuzzy number is defined by its core, its lower limit and its upper limit of the support. Fuzzy pressure (), fuzzy temperature ( ), fuzzy wind speed ( ) and fuzzy angle of incidence ( ) have been defined for every division realised along the railway line in order to get a more detailed climatological model on railway track sections as explained in Section 2.2.
The running time is calculated from the indicated fuzzy parameters and the model input data (train and railway line, Section 2.1). This result will be a fuzzy parameter due to the uncertainty associated with some of the calculation parameters. In the proposed fuzzy model, the punctuality constraint is imposed on the fuzzy running time, which must be less or equal to the objective running time: , The punctuality constraint can be expressed either as a measure of the possibility of arriving at the scheduled time or as a measure of the necessity of fulfiling that objective. Thus, if it is expressed as a measure of the possibility , the expression is: , The corresponding necessity measure, in this case, is 0: 0.
If the punctuality constraint is expressed as a measure of the required necessity , (Figure 4) the expression is: Fuzzy pressure (p), fuzzy temperature (T), fuzzy wind speed (v wind ) and fuzzy angle of incidence (φ) have been defined for every division realised along the railway line in order to get a more detailed climatological model on railway track sections as explained in Section 2.2.
The running time is calculated from the indicated fuzzy parameters and the model input data (train and railway line, Section 2.1). This result will be a fuzzy parameter due to the uncertainty associated with some of the calculation parameters. In the proposed fuzzy model, the punctuality constraint is imposed on the fuzzy running time, which must be less or equal to the objective running time: The punctuality constraint can be expressed either as a measure of the possibility of arriving at the scheduled time or as a measure of the necessity of fulfiling that objective. Thus, if it is expressed as a measure of the possibility α ρ , the expression is: The corresponding necessity measure, in this case, is 0: If the punctuality constraint is expressed as a measure of the required necessity n p , ( Figure 4) the expression is: The corresponding possibility measure, in this case, is 1: Π (T r ≤ t obj ) = 1. The corresponding possibility measure, in this case, is 1: 1.

Eco-Driving Optimisation with Fuzzy Parameters
Having introduced the climatological fuzzy model and how it affects the train motion simulation model, the objective is now to present the driving optimisation model. This optimal driving will be constructed as a combination of efficient commands considering coasting points known as Command Matrix ( ). Then, the problem and proposed algorithm for the problem resolution will be presented: Genetic Algorithm with Fuzzy Parameters (GA-F).

Efficient Driving Commands
In order to design the eco-driving, it is necessary to make use of efficient commands. These commands are based on the coasting application and will be used in this paper with the speed regulation without braking strategy. This driving is the efficient version of the speed regulation strategy and it is easily applied by the drivers on long distance and highspeed railway lines. To execute this strategy, the constant regulation speed is applied as long as it is needed for the traction. It occurs in horizontal or non-steep track alignment. Moreover, braking is not applied if needed to maintain the cruising speed and the maximum speed limits allow it [23,24]. It occurs in steep downhill alignments. This way, the train uses the track grades to increase speed without consuming energy. It allows gaining time that can be used to drive the train at a lower speed at other track sections where energy is needed to maintain the cruising speed. Braking is applied just to observe maximum speed limitations and braking curves up to reductions of maximum speed profile or up to stopping points.
Commands are modelled in the form of Command Matrix ( ) [24]. The matrix is the result of the optimisation process and means the division of the railway line where the driver is required to apply certain driving commands. The is defined as a matrix formed with the values of the points where every stretch end and the value of the regulation speed without braking are applied in every section. The end (Kilometric Point) of the divisions forms the first column, and the regulation speed value without braking forms the second column. Thus, the initial point of a division is the end of the previous section except for the first division, where the initial point is the journey starting position. Furthermore, the end of the last division will be the coasting section's starting point until the braking curve is reached up to the station. The numbers of rows of the matrix plus 1 are the number of the sections defined. Figure 5 shows the design of the driving by the divisions along the railway line where the speed regulation without braking is applied.

Eco-Driving Optimisation with Fuzzy Parameters
Having introduced the climatological fuzzy model and how it affects the train motion simulation model, the objective is now to present the driving optimisation model. This optimal driving will be constructed as a combination of efficient commands considering coasting points known as Command Matrix (C m ). Then, the problem and proposed algorithm for the problem resolution will be presented: Genetic Algorithm with Fuzzy Parameters (GA-F).

Efficient Driving Commands
In order to design the eco-driving, it is necessary to make use of efficient commands. These commands are based on the coasting application and will be used in this paper with the speed regulation without braking strategy. This driving is the efficient version of the speed regulation strategy and it is easily applied by the drivers on long distance and high-speed railway lines. To execute this strategy, the constant regulation speed is applied as long as it is needed for the traction. It occurs in horizontal or non-steep track alignment. Moreover, braking is not applied if needed to maintain the cruising speed and the maximum speed limits allow it [23,24]. It occurs in steep downhill alignments. This way, the train uses the track grades to increase speed without consuming energy. It allows gaining time that can be used to drive the train at a lower speed at other track sections where energy is needed to maintain the cruising speed. Braking is applied just to observe maximum speed limitations and braking curves up to reductions of maximum speed profile or up to stopping points.
Commands are modelled in the form of Command Matrix (C m ) [24]. The matrix is the result of the optimisation process and means the division of the railway line where the driver is required to apply certain driving commands. The C m is defined as a matrix formed with the values of the points where every stretch end and the value of the regulation speed without braking are applied in every section. The end (Kilometric Point) of the divisions forms the first column, and the regulation speed value without braking forms the second column. Thus, the initial point of a division is the end of the previous section except for the first division, where the initial point is the journey starting position. Furthermore, the end of the last division will be the coasting section's starting point until the braking curve is reached up to the station. The numbers of rows of the matrix plus 1 are the number of the sections defined. Figure 5 shows the design of the driving by the divisions along the railway line where the speed regulation without braking is applied. The application of the Command Matrix calculated considers the operational requirements as the speed limits and the stopping points at the stations. Moreover, the optimisation model includes comfort demands for the passengers and punctuality constraints. Hence, the Command Matrix ( ) obtained as a solution to the problem will be efficient, safe, and comfortable driving.
The solutions given by the algorithm are feasible and easily understood, and executed by the driver because they represent realistic driving. The number of driving commands to be executed by the driver is a choice to be made by the railway operator according to the desired complexity during the railway operation. The definition of driving must consider that the more divisions are taken, the more complicated it will be to execute the eco-driving strategy.

Genetic Algorithm with Fuzzy Parameters: GA-F
Once the fuzzy numbers of the climatological model are given, the driving simulated will have an associated fuzzy energy consumption ( ) and a fuzzy running time ( ). When the efficient driving is designed, the objective is to find the set of commands that will produce the minimum value of energy consumption for the required running time ( ). The fitness function ( , ) is shown in Equation (29) and defines the evaluation of both the running time and energy consumption affected by the weighting factors considered ( , ) to minimise the energy consumption considering a design or target running time. Energy consumption associated with the fastest possible driving along the railway line with operational constraints is called . , where  and are weighting factors.  is the energy consumption when the flat-out driving is applied.  is the target running time.
The Genetic Algorithm (GA) used to solve the problem evolves a population of possible solutions (individuals) by means of iterations. An elite group, with the fittest solutions, survives from one generation to the next while the process is moving forward using as the genome of the individuals. Additionally, crossover and mutation operators are applied to generate and offspring of individuals from the elite group at each iteration. The fuzzy optimisation algorithm is the procedure for the calculation of the driving commands when the fuzzy parameters are included. The GA-F process is graphically described in Figure 6 and it is explained as follows: 1. First random generation of individuals ( ) equal to the population size . 2. Simulation of command matrix of the population. This step generates a running time and an energy consumption associated to each individual. (29)). The application of the Command Matrix calculated considers the operational requirements as the speed limits and the stopping points at the stations. Moreover, the optimisation model includes comfort demands for the passengers and punctuality constraints. Hence, the Command Matrix (C m ) obtained as a solution to the problem will be efficient, safe, and comfortable driving.

Evaluation of each individual by the fitness function (Equation
The solutions given by the algorithm are feasible and easily understood, and executed by the driver because they represent realistic driving. The number of driving commands to be executed by the driver is a choice to be made by the railway operator according to the desired complexity during the railway operation. The definition of driving must consider that the more divisions are taken, the more complicated it will be to execute the eco-driving strategy.

Genetic Algorithm with Fuzzy Parameters: GA-F
Once the fuzzy numbers of the climatological model are given, the driving simulated will have an associated fuzzy energy consumption (Ê) and a fuzzy running time (T r ). When the efficient driving is designed, the objective is to find the set of commands that will produce the minimum value of energy consumption for the required running time (t obj ).
The fitness functionF(T r ,Ê) is shown in Equation (29) and defines the evaluation of both the running time and energy consumption affected by the weighting factors considered (w e , w t ) to minimise the energy consumption considering a design or target running time. Energy consumption associated with the fastest possible driving along the railway line with operational constraints is called E f lat−out .
where • w e and w t are weighting factors. • E f lat−out is the energy consumption when the flat-out driving is applied. • t obj is the target running time.
The Genetic Algorithm (GA) used to solve the problem evolves a population of possible solutions (individuals) by means of iterations. An elite group, with the fittest solutions, survives from one generation to the next while the process is moving forward using C m as the genome of the individuals. Additionally, crossover and mutation operators are applied to generate and offspring of individuals from the elite group at each iteration. The fuzzy optimisation algorithm is the procedure for the calculation of the driving commands when the fuzzy parameters are included. The GA-F process is graphically described in Figure 6 and it is explained as follows: 1.
First random generation of individuals (C m ) equal to the population size N pop .

2.
Simulation of command matrix of the population. This step generates a running timê T r and an energy consumptionÊ associated to each individual.

3.
Evaluation of each individual by the fitness function (Equation (29)).

4.
Sorting of each individual in increasing order by means of their fitness value.

5.
Elite group selection. The first individuals N e survive in the next iteration population and the rest are eliminated. 6.
Offspring generation. Mutation and crossover operators are applied to the elite group to generate new individuals until the population reaches its size N pop . In this step, a number of mutations N m and a number of crossovers N c are applied. 7.
The process is repeated from step 2 until the iteration number (it) is equal to the maximum numbers of iterations defined (it max ), giving the best solution to the fittest in the last population.  4. Sorting of each individual in increasing order by means of their fitness value. 5. Elite group selection. The first individuals survive in the next iteration population and the rest are eliminated. 6. Offspring generation. Mutation and crossover operators are applied to the elite group to generate new individuals until the population reaches its size . In this step, a number of mutations and a number of crossovers are applied. 7. The process is repeated from step 2 until the iteration number ( ) is equal to the maximum numbers of iterations defined ( ), giving the best solution to the fittest in the last population. The proposed algorithm can be expressed in terms of the α-cuts of the fuzzy numbers [96]. As shown in Figure 4, if the punctuality requirement is expressed as a necessity measure np, the upper limit of the α-cuts of the fuzzy running time (equivalent to the upper limit of the α-cuts of the fuzzy running time) is necessary for the evaluation of the punctuality constraint. For this upper limit to be obtained, given the fuzzy numbers associated with the climatological parameters, it has to be analysed whether the running time is an increasing or a decreasing function in each of them.
As discussed in the model, the running time and energy consumption are dependent on the running resistance model. The running resistance force is a function of the fuzzy parameters introduced in this model: the air density and the wind during operation. In turn, the density is a function of the pressure and the temperature. The wind has been modelled by its intensity or speed and their angle of incidence related to the train's movement.
With equal , train and railway track data, the running time calculated is a function of the uncertainties considered: pressure (), temperature ( ), velocity of the wind ( ), and angle of incidence of the wind ( ). All climatological factors, running time, and energy consumption are defined as fuzzy parameters.
, , , , Function (Equation (30)) is monotonous with all climatological parameters. The greater the running resistance, the lower the acceleration and the greater the running time as a result. In that sense, function is increasing with the air density, and thus, according to Equation (18), it is monotonically increasing with the pressure () and monotonically decreasing with the temperature ( ). The proposed algorithm can be expressed in terms of the α-cuts of the fuzzy numbers [96]. As shown in Figure 4, if the punctuality requirement is expressed as a necessity measure n p , the upper limit of the α-cuts of the fuzzy running time (equivalent to the upper limit of the α-cuts of the fuzzy running time) is necessary for the evaluation of the punctuality constraint. For this upper limit to be obtained, given the fuzzy numbers associated with the climatological parameters, it has to be analysed whether the running time is an increasing or a decreasing function in each of them.
As discussed in the model, the running time and energy consumption are dependent on the running resistance model. The running resistance force is a function of the fuzzy parameters introduced in this model: the air density and the wind during operation. In turn, the density is a function of the pressure and the temperature. The wind has been modelled by its intensity or speed and their angle of incidence related to the train's movement.
With equal C m , train and railway track data, the running time calculated is a function of the uncertainties considered: pressure (p), temperature (T), velocity of the wind (v wind ), and angle of incidence of the wind (φ). All climatological factors, running time, and energy consumption are defined as fuzzy parameters.
Function F 1 (Equation (30)) is monotonous with all climatological parameters. The greater the running resistance, the lower the acceleration and the greater the running time as a result. In that sense, function F 1 is increasing with the air density, and thus, according to Equation (18), it is monotonically increasing with the pressure (p) and monotonically decreasing with the temperature (T).
Regarding the angle of incidence of the wind (φ), it is necessary to analyse every division of the railway line (explained at the end of Section 2.2). In this case, and taking Figure 1 as a reference, F 1 is defined differently depending on the situation. Moreover, the presence of the angle of incidence in Equation (22) has to be considered in order to analyse the behaviour of the function F 1 . In this way, the variation of the angle and, therefore, the change of cos φ in the term of Equation (22) determines the different situations in the high-speed railway scenario. If the angle of incidence lies between: (1) 0 • <φ < 180 • then F 1 is monotonically decreasing with the angle φ.
After considering the term of the angle of incidence, the speed is another factor related to the wind and its contribution to the railway operation. Considering Figure 1, the behaviour is as follows: (1) 0 • <φ < 90 • or 270 • <φ < 360 • then F 1 is monotonically increasing with the velocity of the wind. (2) 90 • <φ < 270 • then F 1 is monotonically decreasing with the velocity of the wind.
Likewise, the energy consumption associated to the railway operation is also a function of the climatological parameters:Ê Function F 2 (Equation (31)) is monotonous with all climatological parameters as F 1 which permits the application of α-cut calculations. Moreover, the dependence of F 2 with these parameters is the same as F 1 because the greater the running resistance, the greater the traction effort and energy consumption.
In conclusion, the upper and lower limits ofT r α-cuts, and the upper and lower limits ofÊ α-cuts can be calculated by means of the alfa-cut of the parameters. Equations (32)- (39) show the corresponding formulas for the 4 quadrants of the angle of incidence as follows, considering the situations presented previously: • For 0 • <φ < 90 • : • For 90 • <φ < 180 • : • For 180 • <φ < 270 • : • For 270 • <φ < 360 • : Figure 7 shows an example corresponding to the case 180 • <φ < 270 • of how the climatological fuzzy parameters are considered and the results associated with a specific necessity level of punctuality.  Figure 7 shows an example corresponding to the case 180° 270° of how the climatological fuzzy parameters are considered and the results associated with a specific necessity level of punctuality. Therefore, the proposed problem can be solved with the application of the GA-F by means of α-cuts with 1 . The parameter is the imposed punctuality constraint. The algorithm will search for the solution in the form of a Command Matrix with a fuzzy running time associated. Its upper α-cut gives the objective running time as a result, considering punctuality requirements and minimising the energy consumption.

Introduction
Real commercial service and infrastructure data of a Spanish high-speed railway line have been analysed in this paper. The railway line is operated by a train with an unladen mass of 322 t, a length of 200 m, a maximum traction effort of 200 kN, and a maximum speed of 300 km/h. This line has been divided as shown in the following figures to elaborate an appropriate data analysis as explained at the end of Section 2.2. The railway line has been divided into zones where the wind is considered constant as can be seen in Figure 8a, and into sections where the angle of the track (Table 1) is considered constant as shown in Figure 8b.  Therefore, the proposed problem can be solved with the application of the GA-F by means of α-cuts with α = 1 − n p . The parameter n p is the imposed punctuality constraint. The algorithm will search for the solution in the form of a Command Matrix with a fuzzy running time associated. Its upper α-cut gives the objective running time as a result, considering punctuality requirements and minimising the energy consumption.

Introduction
Real commercial service and infrastructure data of a Spanish high-speed railway line have been analysed in this paper. The railway line is operated by a train with an unladen mass of 322 t, a length of 200 m, a maximum traction effort of 200 kN, and a maximum speed of 300 km/h. This line has been divided as shown in the following figures to elaborate an appropriate data analysis as explained at the end of Section 2.2. The railway line has been divided into zones where the wind is considered constant as can be seen in Figure 8a, and into sections where the angle of the track ϕ rail (Table 1) is considered constant as shown in Figure 8b.  Figure 7 shows an example corresponding to the case 180° 270° of how the climatological fuzzy parameters are considered and the results associated with a specific necessity level of punctuality. Therefore, the proposed problem can be solved with the application of the GA-F by means of α-cuts with 1 . The parameter is the imposed punctuality constraint. The algorithm will search for the solution in the form of a Command Matrix with a fuzzy running time associated. Its upper α-cut gives the objective running time as a result, considering punctuality requirements and minimising the energy consumption.

Introduction
Real commercial service and infrastructure data of a Spanish high-speed railway line have been analysed in this paper. The railway line is operated by a train with an unladen mass of 322 t, a length of 200 m, a maximum traction effort of 200 kN, and a maximum speed of 300 km/h. This line has been divided as shown in the following figures to elaborate an appropriate data analysis as explained at the end of Section 2.2. The railway line has been divided into zones where the wind is considered constant as can be seen in Figure 8a, and into sections where the angle of the track (Table 1) is considered constant as shown in Figure 8b.

Fuzzy Climatological Parameters
The climatological parameters have been modelled taking into account the temperature, the wind intensity, and its angle of incidence collected by the State Meteorological Agency of the Spanish Government [97] in every one of the closest weather stations to the railway line for this case study, and the pressure is a function of the track altitude according to the International Standard Atmosphere [87].
The following tables collect the fuzzy climatological numbers modelled: temperature, pressure, intensity, and angle of incidence of the wind (Tables 2-5, respectively).

Flat-Out Driving
Initially, the flat-out driving simulations are presented for the Initial Case scenario, with no climatological parameters, and also for the scenarios when the climatological parameters are considered (both in winter and summer). The flat-out driving is the fastest possible in the journey fulfiling the speed limitations of the line. This type of driving has the fastest running time and the associated highest energy consumption.
The running time in the simulated Initial Case (no climatological parameters) is 2:18:49, and the associated energy consumption is 10.804 MWh. Table 6 presents the results.  Figure 9 shows the speed profile of the train for the Initial Case, when the flat-out driving is applied and climatological parameters are not impacting the operation.

Flat-Out Driving
Initially, the flat-out driving simulations are presented for the Initial Case scenario, with no climatological parameters, and also for the scenarios when the climatological parameters are considered (both in winter and summer). The flat-out driving is the fastest possible in the journey fulfiling the speed limitations of the line. This type of driving has the fastest running time and the associated highest energy consumption.
The running time in the simulated Initial Case (no climatological parameters) is 2:18:49, and the associated energy consumption is 10.804 MWh. Table 6 presents the results.  Figure 9 shows the speed profile of the train for the Initial Case, when the flat-out driving is applied and climatological parameters are not impacting the operation. The running time and the energy consumption shown in Table 6 for winter and summer climatological conditions, correspond to the value calculated for the core of the fuzzy sets of temperature, pressure, wind intensity, and angle of incidence. Figures 10 and 11 show the energy consumption of the flat-out driving and the associated running time when the uncertainty in the climatological parameters are included. The running time and the energy consumption shown in Table 6 for winter and summer climatological conditions, correspond to the value calculated for the core of the fuzzy sets of temperature, pressure, wind intensity, and angle of incidence. Figures 10 and 11 show the energy consumption of the flat-out driving and the associated running time when the uncertainty in the climatological parameters are included.
The impact of climatological parameters on energy consumption is more important than on the running time. In the summer scenario shown in Figure 10, the variation of energy consumption for the minimum and maximum values (support of the fuzzy numbers) is 14.8% in summer and 13.2% in winter. To analyse the impact of winter and summer conditions on the energy consumption, the core of the fuzzy energy consumption in summer and winter can be compared; this comparison shows that the energy consumption is 5.4% less in summer than in winter.
On the other hand, the impact on the running times is less significant. In summer, the variation is 0.18%, and 0.16% in winter. The variation between summer and winter conditions shows a difference of 0.1%.

Eco-Driving Design without Climatological Parameters
After analysing the impact of weather conditions on the flat-out driving, the driving optimisation model is run to design the eco-driving without taking into account the weather conditions, which is solved by means of a GA. This would be the typical approach The impact of climatological parameters on energy consumption is more important than on the running time. In the summer scenario shown in Figure 10, the variation of energy consumption for the minimum and maximum values (support of the fuzzy numbers) is 14.8% in summer and 13.2% in winter. To analyse the impact of winter and summer conditions on the energy consumption, the core of the fuzzy energy consumption in summer and winter can be compared; this comparison shows that the energy consumption is 5.4% less in summer than in winter.
On the other hand, the impact on the running times is less significant. In summer, the variation is 0.18%, and 0.16% in winter. The variation between summer and winter conditions shows a difference of 0.1%.

Eco-Driving Design without Climatological Parameters
After analysing the impact of weather conditions on the flat-out driving, the driving optimisation model is run to design the eco-driving without taking into account the weather conditions, which is solved by means of a GA. This would be the typical approach to solving the problem and will serve as a basis for comparison with the model proposed in this article. The objective travel time for the journey is 2 h and 45 min.
The result of the optimisation is the Command Matrix to apply to obtain the optimal driving, shown in Table 7.  Figure 12 shows the speed profile for the eco-driving obtained from the optimisation model. The driving commands produce the nominal driving and define the commercial eco-driving designed in a specific railway operation. This driving produces an energy consumption of 7.589 MWh.  Figure 12 shows the speed profile for the eco-driving obtained from the optimisation model. The driving commands produce the nominal driving and define the commercial eco-driving designed in a specific railway operation. This driving produces an energy consumption of 7.589 MWh. In Table 8, the energy consumption obtained by the eco-driving design and the associated energy saving compared to the flat-out driving can be seen. Thus, following this typical procedure, the energy saving that the railway operator would expect would be 29.76%. The substantial and significant energy saving is achieved due to the slack time included in the timetable (necessary to satisfy punctuality when a delay arises), and due to the eco-driving design to use the slack time optimally to minimise the energy consumption. Other studies show similar energy savings, 31.27% in [16] and 35.54% in [35]. The next step is the execution of the proposed GA-F model to optimise the driving considering the fuzzy climatological parameters (Tables 2-5) with an objective running time of 2:45:00. The punctuality constraint is expressed as a necessity measure 0.5. The fuzzy optimisation model has been applied for the summer and winter scenarios. Table 9 collects the results (energy consumption and running time) obtained by simulating for different α-cut with the eco-driving commands obtained with the optimisation In Table 8, the energy consumption obtained by the eco-driving design and the associated energy saving compared to the flat-out driving can be seen. Thus, following this typical procedure, the energy saving that the railway operator would expect would be 29.76%. The substantial and significant energy saving is achieved due to the slack time included in the timetable (necessary to satisfy punctuality when a delay arises), and due to the eco-driving design to use the slack time optimally to minimise the energy consumption. Other studies show similar energy savings, 31.27% in [16] and 35.54% in [35]. The next step is the execution of the proposed GA-F model to optimise the driving considering the fuzzy climatological parameters (Tables 2-5) with an objective running time of 2:45:00. The punctuality constraint is expressed as a necessity measure np = 0.5. The fuzzy optimisation model has been applied for the summer and winter scenarios. Table 9 collects the results (energy consumption and running time) obtained by simulating for different α-cut with the eco-driving commands obtained with the optimisation model. Specifically, the lower and upper limits of the α-cuts with α = 0 and with α = 0.5 are shown, as well as the core (α = 1). Table 9. Fuzzy results obtained with the eco-driving design applying GA-F with np = 0.5 in winter and in summer.

Season
Winter Summer For the fuzzy energy consumption results, the difference between the core values in summer and winter is 5.32%. The variation between the lower limit of α-cut = 0 and the upper limit of α-cut = 0 in winter is 12.16% and 13.85% in summer. On the other hand, the variation related to the running time is 00:01:08 hh:mm:ss for the winter scenario (0.69%) and 0:00:59 for the summer scenario (0.60%). Figures 13 and 14 are the fuzzy energy consumption and the fuzzy running times obtained by applying the proposed GA-F model with fuzzy climatological parameters. Figure 14 shows that the objective running time is achieved when α = 0.5 as the necessity measure for punctuality imposed was np = 0.5 = 1 − α. For the fuzzy energy consumption results, the difference between the core values in summer and winter is 5.32%. The variation between the lower limit of α-cut = 0 and the upper limit of α-cut = 0 in winter is 12.16% and 13.85% in summer. On the other hand, the variation related to the running time is 00:01:08 hh:mm:ss for the winter scenario (0.69%) and 0:00:59 for the summer scenario (0.60%). Figures 13 and 14 are the fuzzy energy consumption and the fuzzy running times obtained by applying the proposed GA-F model with fuzzy climatological parameters. Figure 14 shows that the objective running time is achieved when 0.5 as the necessity measure for punctuality imposed was np = 0.5 = 1 − .   To conclude, the typical eco-driving design procedure where weather conditions are not considered would provide expected energy savings of 29.76%, as previously indicated. With the application of the proposed GA-F model with fuzzy climatological parameters, the energy consumption of the new eco-driving designs is lower, and thus, the expected energy savings are higher (34.7% in summer and 31.03% in winter). Weather conditions have a significant impact on energy consumption and to a lesser extent, on journey time, and it is worth incorporating this information into the model in order to adapt driving to these pressure, temperature, and wind conditions on the journey under analysis.

Conclusions
This paper proposes a new model to design eco-driving in railway lines taking into account climatological conditions. The climatological parameters affecting the model are the temperature, pressure, the wind intensity, and its angle of incidence. The uncertainty associated with these climatological parameters is modelled by fuzzy numbers, and the optimisation problem is solved by means of a GA-F, Genetic Algorithm with fuzzy parameters. The analysis of the energy consumption and running time functions depending on these parameters permits an efficient calculation and resolution of the problem applying α-cut simulations.
The proposed model has been applied to a real Spanish high-speed railway line. First, To conclude, the typical eco-driving design procedure where weather conditions are not considered would provide expected energy savings of 29.76%, as previously indicated. With the application of the proposed GA-F model with fuzzy climatological parameters, the energy consumption of the new eco-driving designs is lower, and thus, the expected energy savings are higher (34.7% in summer and 31.03% in winter). Weather conditions have a significant impact on energy consumption and to a lesser extent, on journey time, and it is worth incorporating this information into the model in order to adapt driving to these pressure, temperature, and wind conditions on the journey under analysis.

Conclusions
This paper proposes a new model to design eco-driving in railway lines taking into account climatological conditions. The climatological parameters affecting the model are the temperature, pressure, the wind intensity, and its angle of incidence. The uncertainty associated with these climatological parameters is modelled by fuzzy numbers, and the optimisation problem is solved by means of a GA-F, Genetic Algorithm with fuzzy parameters. The analysis of the energy consumption and running time functions depending on these parameters permits an efficient calculation and resolution of the problem applying α-cut simulations.
The proposed model has been applied to a real Spanish high-speed railway line. First, the flat-out driving has been simulated and the important impact that the climatological parameters have on the energy consumption of the train and to a lesser extent on the running times have been shown.
Then, the efficient driving is designed without and with climatological parameters (applying the proposed GA-F). The eco-driving design is based on applying coasting commands and speed regulation without braking commands. In this case study, it has been shown that the expected energy savings obtained in this process are underestimated when eco-driving is designed without considering climatological parameters. Of course, depending on the specific characteristics of the railway line under study, the climatological conditions could be less or more favorable. However, it can be concluded that weather conditions have a significant impact on energy consumption and to a lesser extent, on journey time, and it is worth incorporating this information into the model in order to adapt driving to the pressure, temperature, and wind conditions on the journey under analysis. The final aim is to maximise energy savings considering as much as possible parameterised data involved in the rail traffic operations.
Future developments will be oriented to the application of the model in real-time to adapt the optimal driving to the current climatological conditions and include other sources of uncertainty related to the traffic conditions.