Risk Assessment Model and Sensitivity Analysis of Ordinary Arterial Highways Based on RSR–CRITIC–LVSSM–EFAST
Abstract
:1. Introduction
- Due to the limited data available, most of the current research focuses on the study of expressways, and there are few studies on the prediction of accident risk levels on ordinary arterial highways. However, there are many ordinary arterial highways in China, and their risk classification is an urgent problem to be solved. At the same time, the research on the risk prediction of ordinary arterial highways mostly focuses on the qualitative point of view.
- Some models have a single risk source and lack practicality. In the analysis of the traffic accident severity, the influencing factors of traffic accidents in the current research are relatively small, and there are many factors that increase accident severity. The application scope and applicability of the model have not been fully explored.
- Road risk assessment and operation management methods used alone have obvious shortcomings and special conditions.
2. Risk Assessment Model Based on RCLE
2.1. Safety Risk Assessment Index System
2.2. Risk Judgment Matrix
2.3. Risk Assessment Model
3. Case Study
3.1. Profile
3.2. Evaluating Indicators
3.3. Risk Assessment
4. Analysis of Evaluation Results
4.1. LSSVM Fitting Results
4.2. EFAST Global Sensitivity Analysis
4.3. Analysis and Discussion
5. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
References
- Yang, Y.; Jin, L. Visualizing Temporal and Spatial Distribution Characteristic of Traffic Accidents in China. Sustainability 2022, 14, 13706. [Google Scholar] [CrossRef]
- Raghvendra, V.C.; Joseph, H.S. Coordinability and Consistency: Application of Systems Theory to Accident Causation and Prevention. J. Loss Prevent. Proc. 2015, 33, 200–212. [Google Scholar]
- Abraha, H.H.; Liyanage, J.P. Review of Theories and Accident Causation Models: Understanding of Human-Context Dyad Toward the Use in Modern Complex Systems. In Proceedings of the 7th World Congress on Engineering Asset Management, Deajeon City, Republic of Korea, 8–9 October 2012; Springer: Berlin, Germany, 2015; pp. 17–32. [Google Scholar]
- Kageyama, I.; Kuriyagawa, Y.; Haraguchi, T.; Kaneko, T.; Nishio, M.; Watanabe, A. Study on Road Friction Database for Traffic Safety: Construction of a Road Friction-Measuring Device. Inventions 2022, 7, 90. [Google Scholar] [CrossRef]
- Khattak, A.; Almujibah, H.; Elamary, A.; Matara, C.M. Interpretable Dynamic Ensemble Selection Approach for the Prediction of Road Traffic Injury Severity: A Case Study of Pakistan’s National Highway N-5. Sustainability 2022, 14, 12340. [Google Scholar] [CrossRef]
- Abdullah, P.; Sipos, T. Drivers’ Behavior and Traffic Accident Analysis Using Decision Tree Method. Sustainability 2022, 14, 11339. [Google Scholar] [CrossRef]
- Nasibeh, A.F.; Anna, S.; Ehsan, R.; Jaime, A.C. Risk Assessment of Occupational Injuries Using Accident Severity Grade. Saf. Sci. 2015, 76, 160–167. [Google Scholar]
- Ye, X.; Pendyala, R.M.; Shankar, V.; Konduri, K.C. A Simultaneous Equations Model of Crash Frequency by Severity Level for Freeway Sections. Accid. Anal. Prev. 2013, 57, 140–149. [Google Scholar] [CrossRef]
- Mujalli, R.O.; De, O.J. A Method for Simplifying the Analysis of Traffic Accidents Injury Severity on Two-lane Highways Using Bayesian Networks. J. Saf. Res. 2011, 42, 317. [Google Scholar] [CrossRef]
- Msallam, M.; Asi, I.; Abudayyeh, D. Safety Evaluation (Skid Resistance) of Jordan’s National Highway Network. Jordan J. Civ. Eng. 2019, 11, 2017–2165. [Google Scholar]
- Hussein, N.; Hassan, R.; Fahey, T. Effect of Pavement Condition and Geometrics at Signalised Intersections on Casualty Crashes. J. Saf. Res. 2021, 76, 276–288. [Google Scholar] [CrossRef]
- Khattk, A.J.; Shamayleh, H. Highway Safety Assessment Through Geographic Information System-based Data Visualization. Comput. Civil Eng. 2005, 19, 407–411. [Google Scholar] [CrossRef]
- Transportation Association of Canada. Drag Model Approach: Developments and Applications; Transportation Research Board: Washington, DC, USA, 2002. [Google Scholar]
- Charly; Mathew. Evaluation of Driving Performance in Relation to Safety on an Expressway Using Field Driving Data. Transp. Lett. 2020, 12, 340–348. [Google Scholar] [CrossRef]
- Orsini, F.; Gecchele, G.; Gastaldi, M.; Rossi, R. Large-scale Road Safety Evaluation Using Extreme Value Theory. IET Intell. Transp. Sy. 2020, 14, 1004–1012. [Google Scholar] [CrossRef]
- Alfonso; Montella. Safety Reviews of Existing Roads: Quantitative Safety Assessment Methodology. Transport. Res. Rec. 2018, 1922, 62–72. [Google Scholar]
- Ait-Mlouk, A.; Agouti, T. DM-MCDA: A Web-based Platform for Data Mining and Multiple Criteria Decision Analysis: A Case Study on Road Accident. Soft. X. 2019, 10, 100323. [Google Scholar] [CrossRef]
- Ekmekcioğlu, Ö.; Koc, K.; Özger, M. Stakeholder Perceptions in Flood Risk Assessment: A Hybrid Fuzzy AHP-TOPSIS Approach for Istanbul, Turkey. Int. J. Disast. Risk Re. 2021, 60, 102327. [Google Scholar] [CrossRef]
- Dai, B.; Li, D.; Zhang, L.; Liu, Y.; Zhang, Z.; Chen, S. Rock Mass Classification Method Based on Entropy Weight–TOPSIS–Grey Correlation Analysis. Sustainability 2022, 14, 10500. [Google Scholar] [CrossRef]
- Raza, A.; Ali, M.U.; Ullah, U.; Fayaz, M.; Alvi, M.J.; Kallu, K.D.; Zafar, A.; Nengroo, S.H. Evaluation of a Sustainable Urban Transportation System in Terms of Traffic Congestion—A Case Study in Taxila, Pakistan. Sustainability 2022, 14, 12325. [Google Scholar] [CrossRef]
- Zhao, Y.; Li, J.; Ying, X. Study on Risk of Long-Steep Downgrade Sections of Expressways Based on a Fuzzy Hierarchy Comprehensive Evaluation. Appl. Sci. 2022, 12, 5924. [Google Scholar] [CrossRef]
- Ntafloukas, K.; McCrum, D.P.; Pasquale, L. A Cyber-Physical Risk Assessment Approach for Internet of Things Enabled Transportation Infrastructure. Appl. Sci. 2022, 12, 9241. [Google Scholar] [CrossRef]
- Yari, M.; Monjezi, M.; Bagherpour, R. A Novel Investigation in Blasting Operation Management Using Decision Making Methods. Rud. Geol. Naft. Zb. 2014, 29, 69–79. [Google Scholar]
- Yu, X.; Suntrayuth, S.; Su, J. A Comprehensive Evaluation Method for Industrial Sewage Treatment Projects Based on the Improved Entropy-TOPSIS. Sustainability 2020, 12, 6734. [Google Scholar] [CrossRef]
- Yari, M.; Bagherpour, R.; Khoshouei, M.; Pedram, H. Investigating a Comprehensive Model for Evaluating Occupational and Environmental Risks of Dimensional Stone Mining. Rud. Geol. Naft. Zb. 2020, 35, 101–109. [Google Scholar] [CrossRef]
- Hayaty, M.; Tavakoli Mohammadi, M.R.; Rezaei, A.; Shayestehfar, M.R. Risk Assessment and Ranking of Metals Using FDAHP and TOPSIS. Mine Water Environ. 2014, 33, 157–164. [Google Scholar] [CrossRef]
- Yari, M.; Bagherpour, R.; Almasi, N. An Approach to the Evaluation and Classification of Dimensional Stone Quarries with an Emphasis on Safety Parameters. Rud. Geol. Naft. Zb. 2016, 31, 15–26. [Google Scholar] [CrossRef] [Green Version]
- Lu, X.; Wang, Y.P.; Ziehn, T.; Dai, Y. An Efficient Method for Global Parameter Sensitivity Analysis and its Applications to the Australian Community Land Surface Model (CABLE). Agr. Forest Meteorol. 2013, 182–183, 292–303. [Google Scholar] [CrossRef]
- Maxd; Morris. Factorial Sampling Plans for Preliminary Computational Experiments. Technometrics 1991, 33, 161–174. [Google Scholar] [CrossRef]
- Cukier, R.I.; Levine, H.B.; Shuler, K.E. Nonlinear Sensitivity Analysis of Multiparameter Model Systems. J Comput. Phys. 1978, 26, 1–42. [Google Scholar] [CrossRef]
- Sobol, I.M. Sensitivity Estimates for Nonlinear Mathematical Models. Math. Mod Comput. Exp. 1993, 1, 407–414. [Google Scholar]
- Li, Y.; Huang, C.; Lu, L. Global Sensitivity Analysis of SEBS Model Parameters Based on EFAST Method. Remote Sen. Tech. Appl. 2014, 29, 719–726. [Google Scholar]
- Saltelli, A.; Tarantola, S.; Chan, K. A Quantitative Model-Independent Method for Global Sensitivity Analysis of Model Output. Technometrics 1999, 41, 39–56. [Google Scholar] [CrossRef]
- Broughton, J.; Knowles, J. Providing the Numerical Context for British Casualty Reduction Targets. Safety Sci. 2010, 48, 1134–1141. [Google Scholar] [CrossRef]
- Hagenzieker, M.P.; Commandeur J., J.F.; Bijleveld, F.D. The history of Road Safety Research: A Quantitative Approach. Transport. Res. F-Traffic. 2014, 25, 150–162. [Google Scholar] [CrossRef]
- Zhu, Y.; Liu, B.; Wang, S.; Yang, Z. Development of a FBG Stress Sensor for Geostress Measurement Using RSR Method in Deep Soft Fractured Rock Mass. Appl. Sci. 2022, 12, 1781. [Google Scholar] [CrossRef]
- Hu, K.; Ge, H.; Li, H.; Xie, S.; Xu, S. Rate-Dependent Hysteresis Modeling and Displacement Tracking Control Based on Least-Squares SVM for Axially Pre-Compressed Macro-Fiber Composite Bimorph. Materials 2022, 15, 6480. [Google Scholar] [CrossRef]
- Hu, L.; Zhang, C.; Zhao, X.; Liu, F.; Lyu, Y.; Xue, Y. A Method for Assessing the Risks of Freeway Segments with Combined Horizontal and Vertical Curves. J. Transport Inform. Saf. 2022, 238, 30–41. [Google Scholar]
- Zhang, J.; Wang, W.; Li, T.; Luo, H. Global Sensitivity Analysis of Highway Beam Bridge Based on EFAST Method. J. Civ. Eng. Manag. 2020, 37, 51–56. [Google Scholar]
- Jiang, L.; Cai, Q.; Zhang, Y.; Shi, H. Global parameter sensitivity analysis of passive system in nuclear power plant based on FAST. J. Naval Univ. Eng. 2018, 30, 65–70. [Google Scholar]
- Tajnik, S.; Luin, B. Impact of Driver, Vehicle, and Environment on Rural Road Crash Rate. Sustainability 2022, 14, 15744. [Google Scholar] [CrossRef]
- Ivajnšič, D.; Horvat, N.; Žiberna, I.; Kotnik, E.K.; Davidović, D. Revealing the Spatial Pattern of Weather-Related Road Traffic Crashes in Slovenia. Appl. Sci. 2021, 11, 6506. [Google Scholar] [CrossRef]
Content | Name | Applicable Conditions |
---|---|---|
MCDM of road risk assessment [18,27] | AHP | The analytic hierarchy process is based on the evaluator’s understanding of the nature and elements of the evaluation problem. It is a more qualitative analysis and judgment method than the general quantitative method. When we want to solve more common problems, the number of indicators selected is likely to increase. |
Fuzzy evaluation method | Fuzzy evaluation deals with fuzzy evaluation objects using precise digital means, and it can make a more scientific, reasonable, and practical quantitative evaluation of fuzzy information. However, the method is complex, and the determination of the index weight vector is subjective. | |
TOPSIS | TOPSIS is a ranking method that approximates the ideal solution. By obtaining the proximity of each evaluation scheme to the optimal scheme, it is used as the basis for the advantages and disadvantages of the scheme. | |
Road risk object identification [19,20,21,22,23,24,25,26] | SVD | SVD is a widely used algorithm in the field of machine learning. It is used for feature decomposition in dimensionality reduction algorithms. It can also be used in recommendation systems, natural language processing, and other fields. In traffic, it can be used to calculate the weight of the judgment matrix in MCDM. |
Cluster analysis | Cluster analysis can be applied to a variety of research scopes, such as regional planning and risk object identification. The determination of the category level of the same type of variable is subjective and objective. | |
RSR method | The RSR method is comprehensive, can show small changes, and is not sensitive to outliers. It is an effective means of comparing and finding relationships by sorting and grading each evaluation object and finding out the advantages and disadvantages. | |
Expert method | The expert method is highly subjective in identifying road risk objects, as it is based on experience. | |
Neural network method | The neural network method is widely used in traffic volume predictions, regression analysis, clustering analysis, etc. It has strong work randomness. | |
Global sensitivity analysis [29,30,31,32,33] | Morris | Morris is based on statistical theory, including the scatter diagram method, the correlation coefficient method, the regression analysis method, and so on. It has strong applicability to linear monotonic model analysis. |
EFAST | EFAST method is a global sensitivity analysis method based on the FAST method combined with the Sobol method. The integral required to calculate the sensitivity index becomes a single variable, saving calculation time. Each order’s sensitivity index can be obtained to evaluate the coupling between several indexes. | |
Sobol | Sobol is a sensitivity analysis method based on variance decomposition. As a typical global sensitivity analysis method, the Sobol method can only evaluate the coupling effect of each index with all other indexes. |
Accident Severity | Risk Possibility | ||
---|---|---|---|
Acceptable Risk (Low Risk) | Tolerable Risk (Medium Risk) | Intolerant Risk (High Risk) | |
Light loss | Low risk of accidents, minor consequences ( * | Accidents not frequent, minor consequences ( ** | Accident-prone, minor consequences ( *** |
Acceptable loss | Low risk of accidents, acceptable consequences ( ** | Accidents not frequent, acceptable consequences ( *** | Accident-prone, acceptable consequences ( **** |
Heavy loss | Low risk of accidents, serious consequences ( *** | Accidents not frequent, serious consequences ( **** | Accident-prone, serious consequences ( ***** |
Index | Meaning | Evaluation Methods |
---|---|---|
X1 | Overload vehicle ratio | . The number of overloaded vehicles, , and the total number of vehicles, , in the sampling survey. |
X2 | Heavy-duty vehicle ratio | Statistics of vehicles in the traffic flow by vehicle type, of which heavy vehicles include large trucks, very large trucks, and container trucks. |
X3 | Ratio of drivers with driving experience less than 5 years | . Drivers with more than 5 years of driving experience in sample survey, ; total drivers’ sample, |
X4 | AADT | Provided by traffic flow observation stations in highway networks. |
X5 | Basic performance of pavement | The basic performance of the pavement includes flatness, ; skid resistance, ; lane width, ; and shoulder width, . |
X6 | Proportion of horizontal and vertical bad linear sections | Poor linear sections of horizontal and vertical sections include sharp bends, steep slopes, continuous downhills, poor sight distance sections, and their combinations. |
X7 | Average road test hazard level | The roadside danger degree is divided into four grades according to the width of the roadside clear area, slope grade, and roadside dangerous goods. |
X8 | Village roads ratio | . Length of road through village, ; total length of road, |
X9 | The ratio of bad weather affecting road sections | Effect of cloudy, rainy, snowy, fog, high-temperature, freezing, dust, and other adverse weather conditions such that the driver‘s line of sight is blocked; the road adhesion coefficient decreased and the road length accounted for the proportion of the total road length. |
X10 | Number of special hydrological and geological disaster points | According to data from monitoring stations, the affected road sections of disasters such as earthquakes, landslides, collapses, debris flows, and roadbed subsidence are counted. |
X11 | Traffic signs coverage | Number of traffic signs per kilometer of road included. |
X12 | Fence coverage | Proportion of road length covered by roadside guardrails. |
X13 | Number of deceleration facilities | The deceleration setting includes deceleration markings, vibration deceleration belts, three-dimensional deceleration markings, etc. |
X14–X19 | Severity dimension indicators represent the number of traffic accidents, car accident rates, percentage of serious accidents, death tolls, car accident death rates, and direct economic loss. | Severity dimension indicators can be provided by the Highway Authority. |
Index | Ordinary Arterial Highway Number | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
G310 | G108 | G327 | G210 | S209 | S107 | S108 | S201 | S205 | S305 | |
X1 | 0.0145 | 0.0278 | 0.0136 | 0.0249 | 0.1652 | 0.2108 | 0.1725 | 0.2519 | 0.1568 | 0.1818 |
X2 | 0.1071 | 0.1009 | 0.1803 | 0.1201 | 0.1317 | 0.2245 | 0.1146 | 0.1869 | 0.1245 | 0.1221 |
X3 | 0.5861 | 0.5365 | 0.5219 | 0.5488 | 0.6017 | 0.5932 | 0.5346 | 0.6105 | 0.6033 | 0.6045 |
X4 | 1341 | 1928 | 4485 | 5624 | 6007 | 1668 | 1874 | 5546 | 2815 | 1457 |
X5 | 1.78 | 0.77 | 1.66 | 1.52 | 1.99 | 1.66 | 2.39 | 2.74 | 1.21 | 1.32 |
X6 | 0.1047 | 0.0431 | 0.0229 | 0.0781 | 0.171 | 0.1481 | 0.0125 | 0.0199 | 0.1321 | 0.0257 |
X7 | 0.1925 | 0.1604 | 0.1568 | 0.2465 | 0.6126 | 0. 1176 | 0.381 | 0.0901 | 0.3925 | 0.2489 |
X8 | 0.037 | 0.2819 | 0.1547 | 0.1324 | 0.3025 | 0.4025 | 0.2198 | 0.264 | 0.4554 | 0.3215 |
X9 | 0.5056 | 0.0707 | 0.1357 | 0.3458 | 0.312 | 0.4025 | 0.2198 | 0.264 | 0.4554 | 0.3217 |
X10 | 42 | 31 | 25 | 16 | 8 | 66 | 21 | 4 | 52 | 31 |
X11 | 4.5945 | 3.0433 | 4.1571 | 3.1245 | 3.3978 | 1.775 | 2.3432 | 2.1869 | 2.456 | 2.647 |
X12 | 0.2885 | 0.3652 | 0.4863 | 0.4332 | 0.5997 | 0.2449 | 0.487 | 0.1131 | 0.4882 | 0.3214 |
X13 | 18 | 16 | 12 | 13 | 6 | 5 | 8 | 7 | 2 | 1 |
X14 | 45 | 82 | 12 | 20 | 121 | 245 | 33 | 26 | 21 | 84 |
X15 | 0.0568 | 0.625 | 0.0128 | 0.0267 | 0.5519 | 4.0242 | 0.04824 | 0.3214 | 0.2524 | 0.6932 |
X16 | 0.5 | 0.2837 | 0.2347 | 0.3645 | 0.0331 | 0.0286 | 0.0606 | 0.3461 | 0 | 0.1247 |
X17 | 2 | 3 | 1 | 4 | 11 | 18 | 9 | 5 | 0 | 6 |
X18 | 0.0284 | 0.226 | 0.1234 | 0.3415 | 0.0182 | 0.1807 | 0.0292 | 0.0445 | 0 | 0.2315 |
X19 | 15,150 | 668.47 | 6549 | 2142 | 235 | 114.08 | 92.22 | 120.15 | 2.3 | 203 |
Index | Weight | Index | Weight |
---|---|---|---|
X1 | 0.0548 | X11 | 0.0452 |
X2 | 0.0487 | X12 | 0.0476 |
X3 | 0.0530 | X13 | 0.0490 |
X4 | 0.0797 | X14 | 0.0392 |
X5 | 0.0525 | X15 | 0.0408 |
X6 | 0.0529 | X16 | 0.0744 |
X7 | 0.0573 | X17 | 0.0416 |
X8 | 0.0418 | X18 | 0.0617 |
X9 | 0.0465 | X19 | 0.0672 |
X10 | 0.0461 |
Ordinary Arterial Highway Number | Possibility Dimension | Severity Dimension | Individual Dimension |
---|---|---|---|
G310 | 0.705 | 0.582 | 0.666 |
G108 | 0.738 | 0.674 | 0.719 |
G327 | 0.732 | 0.768 | 0.745 |
G210 | 0.629 | 0.643 | 0.632 |
S209 | 0.462 | 0.824 | 0.575 |
S107 | 0.397 | 0.534 | 0.448 |
S108 | 0.714 | 0.884 | 0.764 |
S201 | 0.506 | 0.792 | 0.591 |
S205 | 0.425 | 0.988 | 0.601 |
S305 | 0.557 | 0.742 | 0.608 |
Variable | Regression Equation | Correlation Coefficient r | Significance Level p |
---|---|---|---|
0.997 | 0.00001519 | ||
0.9679 | 0.00003481 | ||
0.9815 | 0.000008897 |
Light Risk (>0.791) | Low Loss (0.638~0.791) | Medium Risk (0.485~0.638) | Higher Level of Risk (0.332~0.485) | High Risk (<0.332) | |
---|---|---|---|---|---|
Light loss (>0.967) | S205 | ||||
Low level of loss (0.800~0.967) | S108 | S209 | |||
Medium loss (0.632~0.800) | G108 | G327 | G210, S201, S305 | ||
Higher level of loss (0.464~0.632) | G310 | S107 | |||
Heavy loss (<0.464) |
Index | Mean Value * | Standard Deviation | Diversity | Dispersion Pattern |
---|---|---|---|---|
X1 | 1.388 | 0.191 | 0.137 | Uniform |
X2 | 1.082 | 0.131 | 0.121 | Uniform |
X3 | 1.085 | 0.078 | 0.072 | Uniform |
X4 | 1 | 0.03 | 0.03 | Uniform |
X5 | 1 | 0.012 | 0.012 | normal |
X6 | 1.015 | 0.05 | 0.05 | normal |
X7 | 0.999 | 0.044 | 0.043 | normal |
X8 | 0.788 | 0.044 | 0.044 | Uniform |
X9 | 1.09 | 0.085 | 0.108 | Uniform |
X10 | 1.131 | 0.031 | 0.028 | Uniform |
X11 | 1.034 | 0.032 | 0.022 | normal |
X12 | 1.321 | 0.045 | 0.036 | normal |
X13 | 1 | 0.065 | 0.065 | normal |
Availability | Small Covariance | Medium Covariance | Large Covariance |
---|---|---|---|
YS107 | 0.99 | 0.099 | 0.98 |
YS108 | 0.96 | 0.97 | 0.96 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 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 (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Wang, J.; Ma, C.; Wang, S.; Lu, X.; Li, D. Risk Assessment Model and Sensitivity Analysis of Ordinary Arterial Highways Based on RSR–CRITIC–LVSSM–EFAST. Sustainability 2022, 14, 16096. https://doi.org/10.3390/su142316096
Wang J, Ma C, Wang S, Lu X, Li D. Risk Assessment Model and Sensitivity Analysis of Ordinary Arterial Highways Based on RSR–CRITIC–LVSSM–EFAST. Sustainability. 2022; 14(23):16096. https://doi.org/10.3390/su142316096
Chicago/Turabian StyleWang, Jianjun, Chicheng Ma, Sai Wang, Xiaojuan Lu, and Dongyi Li. 2022. "Risk Assessment Model and Sensitivity Analysis of Ordinary Arterial Highways Based on RSR–CRITIC–LVSSM–EFAST" Sustainability 14, no. 23: 16096. https://doi.org/10.3390/su142316096
APA StyleWang, J., Ma, C., Wang, S., Lu, X., & Li, D. (2022). Risk Assessment Model and Sensitivity Analysis of Ordinary Arterial Highways Based on RSR–CRITIC–LVSSM–EFAST. Sustainability, 14(23), 16096. https://doi.org/10.3390/su142316096