Submitted:
18 April 2024
Posted:
19 April 2024
You are already at the latest version
Abstract
Keywords:
1. Introduction
- For the first time, a high accuracy results implementation of Sequence to Sequence (Seq2seq) Ensemble Deep Transfer Learning for day-ahead (1-24 hours) forecasting is conducted on three distinct datasets from islands of the Greek power system. Although the training dataset of Rhodes exhibits somehow different behavior, compared to the other two datasets, the proposed algorithms perform very satisfactory results. This fact further enhances the performance of the proposed strategies and models. The characteristics of the Rhodes dataset lead to a more robust and comprehensive evaluation of the models, as it introduces variability and challenges that may not be present in the other datasets. This diversity in behavior across datasets provides a more realistic and thorough assessment of the models’ capabilities.
- The results obtained indicate that Deep Transfer Learning (DTL) could provide particular value to both Transmission System Operators (TSO) and Distribution System Operators (DSO), within various regions of the Greek system.
- The application of the models is performed on actual load data with minimal data preprocessing, a fact that creates optimistic conclusions regarding their applicability under real-time conditions.
2. Materials and Methods
2.1. Dataset Analysis
2.2. Data Preprocessing
- Anomaly Detection: Anomaly detection in timeseries involves establishing a baseline of normal behavior through statistical methods or machine learning algorithms, extracting relevant features, and training a model on labeled data to distinguish normal patterns from anomalies. Due to instances of zero consumption during specific hours, likely caused by network faults, these particular values were set equal to the corresponding values from one week prior. This adjustment was made to address the challenge of unexpected situations in the data, as algorithms may struggle to account for such anomalies. The goal is to ensure the optimal training of each model by handling these irregularities in the dataset.
- Filling missing values: Filling missing values in timeseries data is a crucial preprocessing step for anomaly detection. Since anomalies are often identified based on patterns and trends in the data, it’s essential to address gaps caused by missing values.
- Min-Max Scaling: This preprocessing method applied to all datasets in this paper involves Min-Max Scaling, which normalizes data points to a range between 0 and 1. To achieve this, two distinct scalers were employed—one for input and another for output datasets. The primary rationale behind utilizing Min-Max Scaling is its ability to enhance the efficiency of training deep learning models during the training phase, facilitating faster convergence to the optimal solution of the loss function.
- One-Hot Encoding: With this process, numerical data are transformed to cyclical, through trigonometric equations. In this study, day of the week, hour of the day and month of the year were converted to sin and cosine formulation. Figure 4 represents the day of the week transformed in sin and cosine format.
2.3. Feature Creation
- Power in hourly resolution: The sequence of 168 hours of load values for 7 days/one week.
- Cos of Day of Week: The sequence of 168 values of Day of the Week (0-6) converted by One Hot Encoding to cosine type.
- Sin of Day of Week: The sequence of 168 values of Day of the Week (0-6) converted by One Hot Encoding to sin type.
- Cos of Hour of Day: The sequence of 168 values of Hour of Day (0-23), converted by One Hot Encoding to cosine type.
- Sin of Hour of Day: The sequence of 168 values of Hour of Day (0-23), converted by One Hot Encoding to sin type.
- Cos of Month of Year: The sequence of 168 values of Month of the Year (1-12) converted by One- Hot Encoding to cosine type.
- Sin of Month of Year: The sequence of 168 values of Month of the Year (1-12) converted by One- Hot Encoding to sin type.
- IsWeekend: The sequence of 168 values of a dummy variable named “Is Weekend”, with value equal to 0 for working days and 1 for weekends and holidays.
3. Methodology
3.1. Deep Learning Models
3.1.1. Multilayer Perceptron
3.1.2. Convolutional Neural Network
3.1.3. Ensemble Learning Model
3.2. Deep Transfer Learning Forecasting Strategies
- For DTL Case 1, only the dataset of Lesvos was used for the first training phase of the models, specifically for the time period from 2019-01-01 01:00:00 to 2021-12-31 23:00:00. Then, for the second phase, i.e., fine-tuning, the time period from 2019-01-01 01:00:00 to 2021-12-31 23:00:00 of the Chios dataset is utilized.
- For DTL Case 2, the first training phase of the models was based on the datasets from Lesvos and Rhodes, and more specifically, for the time period from 2019-01-01 01:00:00 to 2021-12-31 23:00:00 for each dataset. Then, for the second phase, i.e., fine-tuning, the same time period as in DTL Case 1 is used, from 2019-01-01 01:00:00 to 2021-12-31 23:00:00 of the Chios dataset.
- For DTL Case 3, the training dataset that was used covers the time period from 2019-01-01 01:00:00 to 2022-12-31 23:00:00 from the Lesvos dataset without fine tuning. The choice of Lesvos alone was made because it exhibits a higher correlation and similarity with the corresponding dataset of the Chios timeseries, compared to the Rhodes island.
- For MTDL, the datasets from all three islands were used simultaneously for training. Specifically, the datasets of Rhodes and Lesvos for the time period 2019-01-01 01:00:00 to 2022-12-31 23:00:00, and for Chios dataset, for the time period 2019-01-01 01:00:00 to 2021-12-31 23:00:00 were used.
3.2.1. Deep Transfer Learning Case 1
3.2.2. Deep Transfer Learning Case 2
3.2.3. Deep Transfer Learning Case 3
3.2.4. Multi-Task Deep Learning
- N is the total number of tasks.
- X are the input data.
- is the output for task i.
- are the parameters of the neural network model.
3.3. Optimization Framework
3.4. Evaluation Metrics
3.5. Software environment
4. Results Analysis
4.1. Multilayer Perceptron Results
4.2. Convolutional Neural Networks Results
4.3. Ensemble Learning Model Results
4.4. Results Comparison
4.4.1. Variation of MAE for Each Strategy
- For the DTL strategy with fine-tuning and training data from the Lesvos timeseries, DTL Case 1, it is observed that the MLP and Ensemble models outperform the CNN model.
- In the forecasting strategy with fine-tuning and training data from both the Rhodes and Lesvos timeseries, DTL Case 2, it seems that the MLP and Ensemble models exhibit comparable performance, except for October, where MLP outperforms.
- Regarding for the use of pre-trained models without fine-tuning, DTL Case 3, the CNN significantly lags behind the other two models, with MLP consistently exhibiting the highest prediction accuracy.
- Finally, regarding the Multi-task Deep Learning application strategy, MTDL, it is observed that the MLP model consistently shows inferior results compared to the other two models for all months. Here, the CNN and Ensemble achieve similar accuracy.
4.4.2. Variation of MAPE for Each Strategy
4.4.3. Variation of RMSE for Each Strategy
4.5. Aggregated Results
- In DTL Case 1, it is observed that the Ensemble model achieves the best prediction for the month of June, presenting an MAPE of 5.29%.
- In DTL Case 2, again, the Ensemble model achieves the best prediction, which pertains to the month of February and exhibits an MAPE of 5.31%.
- Regarding the DTL Case 3, the Ensemble model achieves the best prediction in January with an MAPE of 7.85%.
- In MTDL, the CNN model manages the best prediction in January, corresponding to an MAPE of 5.62%.
4.6. Results Discussion
5. Conclusion and Future Study Proposals
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
Abbreviations
| BO | Bayesian Optimization |
| CNN | Convolutional Neural Network |
| DNN | Deep Neural Network |
| DTL | Deep Transfer Learning |
| EDA | Exploratory Data Analysis |
| EDL | Ensemble Deep Learning |
| MAE | Mean Absolute Error |
| MAPE | Mean Absolute Percentage Error |
| MDTL | Multi-task Deep Transfer Learning |
| MLP | Multilayer Perceptron |
| R2 | R-Squared |
| RMSE | Root Mean Squared Error |
| Seq2Seq | Sequence-to-Sequence |
References
- Meng, S.; Li, C.; Tian, C.; Peng, W.; Tian, C. Transfer learning based graph convolutional network with self-attention mechanism for abnormal electricity consumption detection. Energy Reports 2023, 9, 5647–5658. [Google Scholar] [CrossRef]
- Antoniadis, A.; Gaucher, S.; Goude, Y. Hierarchical transfer learning with applications to electricity load forecasting. International Journal of Forecasting 2023. [Google Scholar] [CrossRef]
- Yang, C.; Wang, H.; Bai, J.; He, T.; Cheng, H.; Guang, T.; Yao, H.; Qu, L. Transfer learning enhanced water-enabled electricity generation in highly oriented graphene oxide nanochannels. Nature Communications 2022, 13, 6819. [Google Scholar] [CrossRef] [PubMed]
- Dong, Y.; Xiao, L. A Transfer Learning Based Deep Model for Electrical Load Prediction. 2022 IEEE 8th International Conference on Computer and Communications (ICCC), 2022, pp. 2251–2255. [CrossRef]
- Li, D.; Li, J.; Zeng, X.; Stankovic, V.; Stankovic, L.; Xiao, C.; Shi, Q. Transfer learning for multi-objective non-intrusive load monitoring in smart building. Applied Energy 2023, 329, 120223. [Google Scholar] [CrossRef]
- Peirelinck, T.; Kazmi, H.; Mbuwir, B.V.; Hermans, C.; Spiessens, F.; Suykens, J.; Deconinck, G. Transfer learning in demand response: A review of algorithms for data-efficient modelling and control. Energy and AI 2022, 7, 100126. [Google Scholar] [CrossRef]
- Laitsos, V.; Vontzos, G.; Bargiotas, D. Investigation of Transfer Learning for Electricity Load Forecasting. 2023 14th International Conference on Information, Intelligence, Systems & Applications (IISA). IEEE, 2023, pp. 1–7. [CrossRef]
- Wu, D.; Lin, W. Efficient Residential Electric Load Forecasting via Transfer Learning and Graph Neural Networks. IEEE Transactions on Smart Grid 2023, 14, 2423–2431. [Google Scholar] [CrossRef]
- Syed, D.; Zainab, A.; Refaat, S.S.; Abu-Rub, H.; Bouhali, O.; Ghrayeb, A.; Houchati, M.; Bañales, S. Inductive Transfer and Deep Neural Network Learning-Based Cross-Model Method for Short-Term Load Forecasting in Smarts Grids. IEEE Canadian Journal of Electrical and Computer Engineering 2023, 46, 157–169. [Google Scholar] [CrossRef]
- Laitsos, V.; Vontzos, G.; Bargiotas, D.; Daskalopulu, A.; Tsoukalas, L.H. Enhanced Automated Deep Learning Application for Short-Term Load Forecasting. Mathematics 2023, 11, 2912. [Google Scholar] [CrossRef]
- Santos, M.L.; García, S.D.; García-Santiago, X.; Ogando-Martínez, A.; Camarero, F.E.; Gil, G.B.; Ortega, P.C. Deep learning and transfer learning techniques applied to short-term load forecasting of data-poor buildings in local energy communities. Energy and Buildings 2023, 292, 113164. [Google Scholar] [CrossRef]
- Arvanitidis, A.I.; Bargiotas, D.; Daskalopulu, A.; Kontogiannis, D.; Panapakidis, I.P.; Tsoukalas, L.H. Clustering informed MLP models for fast and accurate short-term load forecasting. Energies 2022, 15, 1295. [Google Scholar] [CrossRef]
- Luo, T.; Tang, Z.; Liu, J.; Zhou, B. A Review of Transfer Learning Approaches for Load, Solar and Wind Power Predictions. 2023 Panda Forum on Power and Energy (PandaFPE), 2023, pp. 1580–1584. [CrossRef]
- Li, S.; Wu, H.; Wang, X.; Xu, B.; Yang, L.; Bi, R. Short-term load forecasting based on AM-CIF-LSTM method adopting transfer learning. Frontiers in Energy Research 2023, 11, 1162040. [Google Scholar] [CrossRef]
- Chan, S.; Oktavianti, I.; Puspita, V. A deep learning cnn and ai-tuned svm for electricity consumption forecasting: Multivariate time series data. 2019 IEEE 10th Annual Information Technology, Electronics and Mobile Communication Conference (IEMCON). IEEE, 2019, pp. 0488–0494. [CrossRef]
- Kontogiannis, D.; Bargiotas, D.; Daskalopulu, A.; Arvanitidis, A.I.; Tsoukalas, L.H. Structural ensemble regression for cluster-based aggregate electricity demand forecasting. Electricity 2022, 3, 480–504. [Google Scholar] [CrossRef]
- Jung, S.M.; Park, S.; Jung, S.W.; Hwang, E. Monthly electric load forecasting using transfer learning for smart cities. Sustainability 2020, 12, 6364. [Google Scholar] [CrossRef]
- Al-Hajj, R.; Assi, A.; Neji, B.; Ghandour, R.; Al Barakeh, Z. Transfer Learning for Renewable Energy Systems: A Survey. Sustainability 2023, 15, 9131. [Google Scholar] [CrossRef]
- Nivarthi, C.P. Transfer Learning as an Essential Tool for Digital Twins in Renewable Energy Systems. arXiv preprint arXiv:2203.05026, arXiv:2203.05026 2022. [CrossRef]
- Miraftabzadeh, S.M.; Colombo, C.G.; Longo, M.; Foiadelli, F. A Day-Ahead Photovoltaic Power Prediction via Transfer Learning and Deep Neural Networks. Forecasting 2023, 5, 213–228. [Google Scholar] [CrossRef]
- Vontzos, G.; Laitsos, V.; Bargiotas, D. Data-Driven Airport Multi-Step Very Short-Term Load Forecasting. 2023 14th International Conference on Information, Intelligence, Systems & Applications (IISA). IEEE, 2023, pp. 1–6. [CrossRef]
- Yang, M.; Liu, Y.; Liu, Q. Nonintrusive residential electricity load decomposition based on transfer learning. Sustainability 2021, 13, 6546. [Google Scholar] [CrossRef]
- Li, K.; Wei, B.; Tang, Q.; Liu, Y. A Data-Efficient Building Electricity Load Forecasting Method Based on Maximum Mean Discrepancy and Improved TrAdaBoost Algorithm. Energies 2022, 15, 8780. [Google Scholar] [CrossRef]
- Publication of NII Daily Energy Planning Data | HEDNO. https://deddie.gr/en/themata-tou-diaxeiristi-mi-diasundedemenwn-nisiwn/leitourgia-mdn/dimosieusi-imerisiou-energeiakou-programmatismou/. Accessed: (2023-01-10).



















| Month | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Model | Jan | Feb | Mar | Apr | May | June | July | Aug | Sep | Oct | Nov | Dec | |
| MLP | MAE | 2.16 | 1.56 | 2.01 | 1.38 | 1.12 | 1.12 | 1.57 | 2.12 | 1.40 | 1.39 | 1.64 | 1.42 |
| MAPE | 7.38 | 6.24 | 7.68 | 7.98 | 6.97 | 6.05 | 7.01 | 9.23 | 7.66 | 9.04 | 9.39 | 7.00 | |
| RMSE | 2.79 | 1.98 | 2.59 | 1.78 | 1.41 | 1.48 | 2.15 | 2.76 | 1.78 | 1.70 | 2.06 | 1.79 | |
| R2 | 0.83 | 0.87 | 0.82 | 0.67 | 0.73 | 0.74 | 0.63 | 0.08 | 0.75 | 0.53 | 0.70 | 0.82 | |
| CNN | MAE | 1.69 | 1.52 | 1.73 | 1.31 | 1.06 | 1.24 | 1.63 | 2.52 | 1.43 | 1.75 | 1.54 | 1.38 |
| MAPE | 5.57 | 6.07 | 6.60 | 7.57 | 6.55 | 6.69 | 7.46 | 10.95 | 7.79 | 11.39 | 8.85 | 6.79 | |
| RMSE | 2.37 | 2.03 | 2.28 | 1.78 | 1.38 | 1.74 | 2.25 | 3.19 | 1.90 | 2.20 | 1.98 | 1.79 | |
| R2 | 0.88 | 0.87 | 0.86 | 0.68 | 0.74 | 0.64 | 0.59 | 0.11 | 0.70 | 0.21 | 0.73 | 0.82 | |
| ELM | MAE | 1.76 | 1.42 | 1.77 | 1.16 | 0.86 | 0.98 | 1.38 | 2.20 | 1.22 | 1.23 | 1.37 | 1.14 |
| MAPE | 5.87 | 5.72 | 6.77 | 6.73 | 5.36 | 5.29 | 6.24 | 9.57 | 6.67 | 8.04 | 7.82 | 5.62 | |
| RMSE | 2.42 | 1.84 | 2.31 | 1.57 | 1.12 | 1.36 | 1.91 | 2.80 | 1.58 | 1.57 | 1.73 | 1.51 | |
| R2 | 0.88 | 0.89 | 0.86 | 0.75 | 0.83 | 0.78 | 0.71 | 6.25 | 0.80 | 0.61 | 0.80 | 0.87 | |
| Month | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Model | Jan | Feb | Mar | Apr | May | June | July | Aug | Sep | Oct | Nov | Dec | |
| MLP | MAE | 1.91 | 1.51 | 1.91 | 1.34 | 0.98 | 1.11 | 1.44 | 2.04 | 1.28 | 1.28 | 1.51 | 1.32 |
| MAPE | 6.50 | 6.07 | 7.29 | 7.79 | 6.24 | 6.01 | 6.61 | 8.87 | 6.95 | 8.38 | 8.66 | 6.51 | |
| RMSE | 2.44 | 1.94 | 2.47 | 1.80 | 1.29 | 1.50 | 1.98 | 2.67 | 1.64 | 1.63 | 1.89 | 1.71 | |
| R2 | 0.87 | 0.88 | 0.84 | 0.67 | 0.77 | 0.73 | 0.69 | 0.15 | 0.78 | 0.57 | 0.75 | 0.84 | |
| CNN | MAE | 1.85 | 1.40 | 1.91 | 1.43 | 1.17 | 1.37 | 2.14 | 2.40 | 1.66 | 2.39 | 1.90 | 1.46 |
| MAPE | 6.26 | 5.60 | 7.29 | 8.30 | 7.27 | 7.38 | 9.85 | 10.44 | 9.02 | 15.51 | 10.90 | 7.23 | |
| RMSE | 2.56 | 1.83 | 2.43 | 1.80 | 1.46 | 1.79 | 2.76 | 2.96 | 2.06 | 2.96 | 2.36 | 1.90 | |
| R2 | 0.86 | 0.89 | 0.84 | 0.67 | 0.71 | 0.62 | 0.39 | 0.02 | 0.66 | 0.08 | 0.61 | 0.80 | |
| ELM | MAE | 1.76 | 1.32 | 1.82 | 1.20 | 0.95 | 1.14 | 1.57 | 2.10 | 1.30 | 1.71 | 1.51 | 1.27 |
| MAPE | 6.03 | 5.31 | 6.94 | 6.95 | 5.92 | 6.17 | 7.23 | 9.10 | 7.09 | 11.17 | 8.68 | 6.26 | |
| RMSE | 2.38 | 1.72 | 2.33 | 1.57 | 1.23 | 1.52 | 2.12 | 2.66 | 1.65 | 2.12 | 1.89 | 1.65 | |
| R2 | 0.88 | 0.91 | 0.85 | 0.75 | 0.80 | 0.72 | 0.64 | 0.16 | 0.78 | 0.27 | 0.75 | 0.85 | |
| Month | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Model | Jan | Feb | Mar | Apr | May | June | July | Aug | Sep | Oct | Nov | Dec | |
| MLP | MAE | 2.56 | 2.65 | 2.58 | 2.51 | 2.12 | 2.01 | 2.12 | 2.42 | 2.03 | 2.10 | 2.34 | 2.37 |
| MAPE | 8.75 | 10.60 | 9.87 | 14.54 | 13.12 | 10.86 | 9.74 | 10.51 | 11.05 | 13.71 | 13.43 | 11.69 | |
| RMSE | 3.14 | 3.32 | 3.38 | 3.06 | 2.56 | 2.48 | 2.71 | 3.03 | 2.44 | 2.62 | 2.91 | 2.94 | |
| R2 | 0.79 | 0.64 | 0.69 | 0.05 | 0.10 | 0.27 | 0.41 | -0.10 | 0.52 | -0.11 | 0.41 | 0.51 | |
| CNN | MAE | 3.02 | 3.69 | 3.53 | 4.81 | 3.62 | 2.71 | 2.49 | 2.95 | 2.96 | 4.55 | 4.45 | 4.59 |
| MAPE | 10.34 | 14.79 | 13.48 | 27.82 | 22.40 | 14.63 | 11.45 | 12.82 | 16.09 | 29.63 | 25.47 | 22.62 | |
| RMSE | 4.39 | 5.37 | 5.10 | 6.25 | 4.75 | 3.70 | 3.27 | 3.67 | 3.58 | 5.66 | 5.93 | 6.59 | |
| R2 | 0.58 | 0.05 | 0.29 | -2.96 | -2.08 | -0.64 | 0.14 | -0.61 | -0.05 | -4.18 | -1.45 | -1.45 | |
| ELM | MAE | 2.29 | 2.80 | 2.68 | 2.93 | 2.41 | 1.98 | 1.92 | 2.44 | 2.23 | 3.09 | 3.18 | 3.21 |
| MAPE | 7.85 | 11.22 | 10.22 | 16.96 | 14.91 | 10.68 | 8.84 | 10.58 | 12.16 | 20.10 | 18.18 | 15.84 | |
| RMSE | 3.07 | 3.76 | 3.65 | 3.82 | 3.06 | 2.61 | 2.55 | 3.09 | 2.67 | 3.78 | 4.05 | 4.27 | |
| R2 | 0.80 | 0.54 | 0.64 | -0.48 | -0.28 | 0.19 | 0.48 | -0.14 | 0.42 | -1.30 | -0.14 | -0.03 | |
| Month | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Model | Jan | Feb | Mar | Apr | May | June | July | Aug | Sep | Oct | Nov | Dec | |
| MLP | MAE | 3.21 | 2.74 | 2.81 | 1.98 | 1.66 | 1.69 | 1.98 | 2.67 | 1.88 | 2.22 | 2.75 | 2.80 |
| MAPE | 10.97 | 10.99 | 10.73 | 11.45 | 10.31 | 9.10 | 9.06 | 11.50 | 10.23 | 14.45 | 15.75 | 13.80 | |
| RMSE | 3.75 | 3.24 | 3.36 | 2.47 | 2.00 | 2.16 | 2.61 | 3.43 | 2.32 | 2.64 | 3.31 | 3.27 | |
| R2 | 0.70 | 0.66 | 0.69 | 0.38 | 0.46 | 0.44 | 0.45 | -0.40 | 0.56 | -0.12 | 0.24 | 0.40 | |
| CNN | MAE | 1.67 | 1.69 | 0.80 | 1.49 | 1.28 | 1.20 | 1.46 | 2.14 | 1.22 | 1.60 | 1.75 | 1.57 |
| MAPE | 5.62 | 6.78 | 6.86 | 8.63 | 7.97 | 6.47 | 6.58 | 9.30 | 6.66 | 10.46 | 10.03 | 7.75 | |
| RMSE | 2.30 | 2.22 | 2.45 | 1.99 | 1.62 | 1.67 | 1.98 | 2.87 | 1.51 | 1.88 | 2.04 | 2.02 | |
| R2 | 0.89 | 0.84 | 0.84 | 0.60 | 0.64 | 0.67 | 0.69 | 0.01 | 0.82 | 0.43 | 0.70 | 0.77 | |
| ELM | MAE | 2.03 | 1.85 | 2.03 | 1.35 | 0.96 | 1.17 | 1.49 | 2.28 | 1.14 | 1.00 | 1.38 | 1.56 |
| MAPE | 6.86 | 7.42 | 7.77 | 7.80 | 5.96 | 6.33 | 6.86 | 9.89 | 6.22 | 6.56 | 7.92 | 7.69 | |
| RMSE | 2.61 | 2.28 | 2.56 | 1.80 | 1.21 | 1.57 | 2.05 | 2.96 | 1.46 | 1.30 | 1.84 | 1.95 | |
| R2 | 0.85 | 0.83 | 0.82 | 0.67 | 0.80 | 0.71 | 0.66 | -0.05 | 0.83 | 0.73 | 0.77 | 0.79 | |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).