Research Article | | Peer-Reviewed

Quantitative Assessment of Ship Oil Spill Risk in Arctic Ice Area Based on MCDM Assessment Model

Received: 17 July 2025     Accepted: 4 August 2025     Published: 26 August 2025
Views:       Downloads:
Abstract

Given the intensifying effects of climate change and the escalating human activities in the Arctic, such as shipping and resource extraction, the likelihood of oil spills in the region has increased substantially. The assessment of oil spill risks in Arctic ice-covered regions requires careful consideration of multiple interacting factors, including environmental parameters, socioeconomic vulnerabilities, and emergency response capabilities. There is an urgent need for a systematic assessment approach to support informed risk management decisions. This research aims to develop a multi-criteria decision-making framework to quantify the key drivers of Arctic oil spill risk and classify risk levels, thereby providing a basis for targeted intervention strategies. The intricate nature of Arctic ecosystems, combined with extreme weather conditions and sparse infrastructure, renders conventional risk models insufficient. In response to this challenge, our study employs sophisticated, integrative methodologies to more effectively capture the interrelationships between various risk factors. By integrating the DEMATEL-ANP method, this paper analyzes the causal relationships and relative weights among 14 indicators across four interrelated network layers. These methods not only clarify the causal chains between factors but also account for the intricate interrelations and dependencies within the system, offering a more comprehensive and adaptable decision-making framework. This integration leverages the strengths of both methods, enhancing the depth and accuracy of the analysis. A fuzzy comprehensive evaluation (FCE) approach is applied to convert expert judgments and derived weights into 5 categorical risk levels. The fuzzy logic component handles uncertainties inherent in expert elicitation, ensuring robustness despite data scarcity. The results indicate that ecological and socioeconomic vulnerability contribute most significantly to overall risk, while environmental factors such as sea ice coverage severely constrain emergency response effectiveness. The weighted FCE score identifies the Arctic oil spill scenario as high risk. The integrated DEMATEL-ANP-FCE framework effectively addresses the challenges of limited data availability in complex risk environments and provides a scientific foundation for Arctic oil spill monitoring, emergency planning, and policy formulation.

Published in International Journal of Transportation Engineering and Technology (Volume 11, Issue 3)
DOI 10.11648/j.ijtet.20251103.11
Page(s) 98-110
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2025. Published by Science Publishing Group

Keywords

Arctic Oil Spill, Risk Assessment, Multi-criteria Decision-making, DEMATEL-ANP-FCE, Quantitative Risk Analysis, DANP Model, Fuzzy Comprehensive Evaluation

1. Introduction
With the deepening of economic globalization and the accelerating retreat of Arctic sea ice due to global climate change, the Arctic region is increasingly emerging as a strategic maritime corridor of global significance. The opening of Arctic shipping routes-particularly the Northern Sea Route (NSR)-has the potential to reduce transit distances between Asia and Europe by up to 40%, thereby significantly lowering transportation costs, fuel consumption, and carbon emissions . These advantages have attracted growing interest from the international shipping industry, prompting a notable increase in vessel traffic through Arctic waters.
However, this surge in navigational activity has concurrently heightened the risk of marine accidents such as collisions, groundings, and mechanical failures, many of which can lead to oil spills. In the ecologically fragile Arctic environment, oil spills pose severe and often irreversible threats to marine ecosystems, particularly the habitats of marine mammals and the livelihoods of Indigenous communities that depend on them. The region’s limited self-purification capacity and the complex pathways of pollutant dispersion under ice-covered conditions further exacerbate the environmental consequences of such incidents .
Furthermore, oil discharged from vessels due to operational errors or maritime accidents may enter the Arctic food chain if not promptly contained, ultimately contributing to regional and even global ecological disruptions . As a result, the development of a scientifically robust, multi-dimensional risk assessment framework-one that accounts for spill characteristics, environmental conditions, and socio-ecological vulnerability-is urgently needed. Such a framework is essential for quantifying the risks of oil spills in Arctic ice-covered waters, identifying priority risk areas, and supporting proactive environmental governance and emergency response planning in the face of increasing Arctic maritime activity.
2. Literature Review
Due to the unique geographical and climatic conditions of the Arctic region-such as sea ice drift, extreme weather, and limited accessibility-oil spill emergency response and ecological restoration remain particularly challenging. Scholars at home and abroad have conducted relevant studies, primarily focusing on constructing multi-dimensional quantitative models to assess the environmental, economic, and social impacts of oil spill incidents . Methods such as Analytic Network Process (ANP) and Bayesian Networks have been employed to conduct multi-factor integrated evaluations . However, these models often lack a detailed characterization of ice conditions, especially in terms of the coupling mechanisms between sea ice coverage, ice concentration, ice thickness, and emergency response deployment capacity.
In the field of risk assessment modeling, Bagheri et al applied the Delphi method to evaluate the impact of oil spills in the Arctic , while Lin et al employed fuzzy mathematics to analyze oil spill risks in the region . Tabesh et al utilized Monte Carlo simulations , and Cai et al adopted the Analytic Hierarchy Process (AHP) to perform risk assessments for Arctic oil spill scenarios. The Battelle Memorial Institute, through its Geneva Research Center, adopted the Decision-Making Trial and Evaluation Laboratory (DEMATEL) method in project research. This efficient approach integrates expert consensus (via brainstorming) to precisely identify causal relationships among multiple indicators . Wang et al proposed an innovative model combining satisfaction-importance analysis and DEMATEL to identify factors potentially hindering project performance . Ou Yang et al developed an information security risk control model that integrates VIKOR, DEMATEL, and ANP to resolve feedback conflicts and dependencies among evaluation indicators in organizational environments . In efforts to improve marine insurance regulation concerning Arctic oil spill incidents, Afenyo et al developed a socioeconomic impact model to assess the implications of oil spills on equitable compensation from a global governance perspective .
Through reviewing the current body of domestic research, it is observed that most studies have focused on ecological impact evaluations and post-incident management recommendations. However, limited attention has been given to the proactive regulatory role of emergency response capacity during the initial phase of oil spill events. Furthermore, factors such as human settlement density, which is crucial to understanding the long-term social and economic consequences, have not been systematically analyzed. The sensitivity of ecological systems and the degree of regional economic development also play a critical role in amplifying or mitigating spill-related damage over time. In addition, the physicochemical properties of different oil types directly influence their diffusion pathways and removal efficiency in ice-covered marine environments , yet these key interdependencies are often underrepresented in mainstream risk assessment frameworks.
In response, this study proposes an integrated risk assessment framework that combines the DEMATEL-ANP method with the Fuzzy Comprehensive Evaluation (FCE) model. The framework is structured around four interrelated network layers: oil spill characteristics, environmental conditions, emergency response capacity, and ecological-socioeconomic vulnerability, encompassing 14 key influencing factors. The DEMATEL approach is first employed to identify the causal relationships among factors and construct a directed influence graph. Subsequently, ANP is used to compute the global weights of each factor, revealing systemic interaction mechanisms-such as the suppressive effect of emergency response capacity on spill characteristics and the delaying impact of ice conditions on response time.
Furthermore, the FCE method incorporates triangular fuzzy numbers to account for the subjectivity and uncertainty of expert judgment, transforming qualitative assessments into five-tiered quantitative risk indicators. This integrated framework offers theoretical and practical support for Arctic maritime environmental governance, early warning system development, and the enhancement of liability and compensation mechanisms for oil spill incidents.
3. Materials and Methods
3.1. Index System for Arctic Oil Spill Risk Assessment
Figure 1. Arctic study area.
The Scott Inlet cold seep, located off the coast of Baffin Island in the Nunavut region, is recognized as one of the most prominent subsea petroleum seepage sites, continuously releasing hydrocarbons into the Arctic marine environment. With the increase in maritime traffic, the surrounding area has emerged as a critical shipping hub. Owing to its unique geographical location, this region has been identified in multiple studies as a representative case area for oil spill response under extreme environmental conditions , as shown in the region highlighted in Figure 1.
In risk assessment modeling, the development of a scientifically sound and well-structured index system is the first and most critical step, as its quality directly affects the effectiveness and applicability of the entire evaluation process. In the case of Arctic oil spills, an event characterized by high complexity and uncertainty, it is essential to identify key impact factors systematically across multiple dimensions, including environmental conditions, operational constraints, and socioeconomic sensitivity. Based on the typical structure of the Analytic Network Process, this research constructs a multi-layered network model to represent the risk assessment mechanism of Arctic oil spills. Through a comprehensive review of relevant literature, Arctic oil spill risk was defined as the primary control layer, under which four interrelated network levels were established and further refined into 14 specific evaluation indicators, as shown in Figure 2. The network layer representing oil spill characteristics (S1) includes oil type (S11), spill volume (S12), and sea ice coverage rate (S13), which collectively reflect the scale and potential hazard of the spill. The environmental conditions of ice-covered regions (S2) capture the specific Arctic context that contributes to oil spill events; indicators such as ice concentration (S21), ice thickness (S22), and wind speed (S23) describe the physical and meteorological characteristics of the region. The emergency response capacity (S3) layer encompasses response time (S31), availability of equipment (S32), personnel deployment (S33), and oil recovery efficiency (S34), providing a measure of the operational effectiveness in mitigating spill impacts. Additionally, the ecological and socio-economic vulnerability (S4) layer incorporates indicators such as ecological sensitivity (S41), population density (S42), degree of economic development (S43), and navigational priority (S44), reflecting the region’s potential exposure and resilience to oil spill incidents.
Figure 2. Structure of the indicator system.
3.2. Weight Calculation Using Integrated DEMATEL-ANP Method
The Decision-Making Trial and Evaluation Laboratory (DEMATEL) is a graph-based multi-criteria decision-making (MCDM) method that explores interdependencies and effective influence values among evaluation factors by constructing causal relationship diagrams . Given the non-independence of evaluation indicators in Arctic ice-region oil spill risk assessment, the DEMATEL method is employed to identify and visualize the causal interactions among indicators. The resulting influence structure is then used as the foundation for calculating indicator weights using the Analytic Network Process (ANP), thereby forming a combined DEMATEL-ANP approach. Unlike the Analytic Hierarchy Process (AHP), which assumes mutual independence among indicators when calculating weights, the ANP method allows for interdependencies and feedback among elements .
In this study, a 0-4 scale was adopted to capture the degree of influence a single indicator exerts on another. A panel of four polar safety experts was invited to complete the pairwise influence assessments. The internal consistency of the expert questionnaire was tested using Cronbach’s alpha, yielding a reliability coefficient of 0.97, indicating excellent consistency and reliability of the evaluation scale.
3.2.1. Procedure for Applying the DEMATEL Approach
This method effectively integrates graph theory and matrix tools, providing robust support for analyzing complex relationships among factors within a system. The specific computational steps are as follows:
(1) Construct the Direct Influence Matrix
A pairwise scoring system (0-4) is used to assign values to the indicators, where 0 to 4 represent No Influence, Weak Influence, Moderate Influence, Strong Influence, and Very Strong Influence, respectively. The direct influence matrix X is obtained:
Xn×n=x11x12x1nx21x22x2nxn1xn2xnn(1)
(2) Compute the Normalized Matrix
The direct influence matrix is normalized by dividing each element by the maximum row sum, yielding the normalized matrix N:
N=Xmax1≤i≤nj=1nxij(2)
(3) Compute the Total Influence Matrix
The total influence matrix T is derived from the normalized matrix N, where I is the identity matrix:
T=N(I-N)-1(3)
(4) Calculate Influence Measures
Influence Degree (R): Sum of the row elements in T, indicating the total influence exerted by a factor on others.
Ri=j=1ntij, (i=1,2,…,n)(4)
Dependence Degree (D): Sum of the column elements in T, indicating the total influence received by a factor from others.
Dj=i=1ntji, (j=1,2,…,n)(5)
Cause Degree (R - D): If positive, the factor is a cause (net influencer); if negative, it is an effect (net receiver).
Centrality (R + D): Reflects the factor’s overall importance in the system.
(5) Plot the Cause-Effect Diagram
A causal diagram is drawn with Centrality (R + D) as the x-axis and Cause Degree (R - D) as the y-axis. Centrality quantifies a factor’s relative prominence within the network structure, with elevated values signifying greater systemic importance. Cause Degree evaluates a factor’s function as a causal driver within the system, where higher values reflect stronger influence over other interconnected elements.
3.2.2. Weight Derivation via the Integrated DEMATEL-ANP Method
(1) Computation of the Unweighted Supermatrix
Following pairwise comparisons of element importance, normalized eigenvectors wi are derived for each element group. These eigenvectors represent local priorities within their respective clusters. The unweighted supermatrix Wij is constructed by assembling these normalized vectors, where each column sums to unity, reflecting the relative importance of element ei in cluster j as a secondary criterion. The mathematical formulation is:
Wij=wi1j1wi1j2wi1jnwi2j1wi2j2wi2jnwinj1winj2winjn(6)
(2) Derivation of the Weighted Supermatrix
The unweighted supermatrix Wij captures intra-cluster priorities but lacks inter-cluster comparability. To address this, a cluster-weighted matrix Hij is constructed from normalized eigenvectors of inter-cluster pairwise comparisons. The weighted supermatrix Wa is obtained by:
Wa=Hij×Wij (i,j=1,2,…,n)(7)
(3) Calculation of the Limit Supermatrix
To derive stable weights, the weighted supermatrix Wa is raised to a sufficiently large power k until convergence (Markov chain principle), yielding the limit supermatrix W:
W=limk→∞ (Wa)k(8)
(4) Computation of Integrated Weights
The final weights ω are obtained by multiplying the limit supermatrix-derived element weights w with the criterion-layer weights γ:
ω=wγ(9)
3.3. Development of the Fuzzy Risk Evaluation Model
In view of the inherent fuzziness and uncertainty in oil spill risk assessment within the Arctic region, this research adopts the Fuzzy Comprehensive Evaluation (FCE) method to enhance the scientific rigor and rationality of risk identification and classification .
(1) Establishment of Fuzzy Evaluation Factor Sets and Weight Sets
Based on the hierarchical structure of indicators and their respective weights, the evaluation factor sets and weight sets are constructed as shown in Table 1.
Table 1. Factor sets and weight sets.

Factor sets

Weight sets

S = [S1, S2, S3, S4]

WS=[wS1,wS2,wS3,wS4]

S1 = [S11, S12, S13]

WS1=[wS11,wS12,wS13]

S2 = [S21, S22, S23]

WS2=[wS21,wS22,wS23]

S3 = [S31, S32, S33, S34]

WS3=[wS31,wS32,wS33,wS34]

S4 = [S41, S42, S43, S44]

WS4=[wS41,wS42,wS43,wS44]

(2) Construction of the oil Spill Risk Evaluation Set
Table 2. Risk scale distributions.

Risk level

Extremely high

High

Medium

Low

Ultra-low

Score

5

4

3

2

1

Score range

[4.5, 5]

[4, 4.5]

[3, 4]

[2, 3]

[0, 2]

The risk level refers to the set of all potential risk outcomes associated with the indicators of Arctic oil spill impacts. Depending on the characteristics of the indicators and practical circumstances, the risk criteria can be classified into different levels. In this research, the risk levels are defined as V = {Extremely High, High, Medium, Low, Ultra-low}, with each grade corresponding to a specific score, as detailed in Table 2.
(3) Formation of the Fuzzy Evaluation Matrix
Multiple evaluators are invited to assess the indicator layer based on the evaluation grade set. After quantifying the indicators, the membership degree Fij for the i-th factor with respect to the j-th evaluation grade is calculated as the proportion of evaluators who assign the j-th grade to the i-th factor relative to the total number of evaluators. The fuzzy relation matrix is then constructed as follows:
Fij=F11F12F13F1jF21F22F23F2jFi1Fi2Fi3Fij(10)
(4) Calculation of Membership Degrees
The weighted average operator model M(∙,) is employed to multiply the weight set W with the fuzzy matrix F, yielding the membership degree matrix for the higher-level indicators. This process is repeated iteratively by multiplying the higher-level weight sets with their corresponding membership degree matrices until the membership degree matrix Z for the target layer is obtained.
(5) Computation of Final Scores
Figure 3. Flowchart of the Arctic oil spill risk assessment model.
The comprehensive evaluation score E is derived by multiplying the target-layer membership degree matrix Z with the predefined score vector V:
Score=k=1nZkVk(11)
The specific methodological steps are illustrated in Figure 3.
4. Results
4.1. Interdependencies Among Risk Factors
The integrated influence matrix for the selected indicator layer is presented in Table 3. Each cell in the matrix reflects the degree of influence exerted by the row indicator on the corresponding column indicator.
Based on the matrix in Table 3, both the centrality-causality chart and the causal relationship diagram of the indicators are derived, as shown in Figure 4. In Figure 4(b), arrows indicate the direction of influence, pointing from the influencing factor to the influenced one. In Figure 4(a), The indicators within the environmental conditions group exhibit positive causality values, indicating that they act as causal factors affecting other groups. In contrast, the indicators in the emergency response capacity group have negative causality values, suggesting they belong to the result group and are predominantly influenced by other factor groups.
Figure 4. DEMATEL results. (a) Centrality-causality diagram. (b) Causal relationship network.
Indicators with higher centrality values, such as oil recovery efficiency, demonstrate greater overall importance within the evaluation system. Centrality analysis allows for the identification of relatively significant indicators in Arctic oil spill risk assessment, while causality analysis helps determine those indicators with broader influence across the system.
Table 3. Total influence matrix of the indicator layer.

S11

S12

S13

S21

S22

S23

S31

S32

S33

S34

S41

S42

S43

S44

S11

0.0706

0.1301

0.0593

0.0562

0.0391

0.0553

0.2398

0.1981

0.2426

0.268

0.2032

0.1738

0.1969

0.1378

S12

0.0664

0.1279

0.0818

0.0543

0.0382

0.0517

0.2376

0.1982

0.2411

0.2181

0.177

0.1742

0.2215

0.1635

S13

0.0969

0.1817

0.0972

0.1659

0.114

0.1354

0.3507

0.2967

0.351

0.357

0.2418

0.2072

0.2697

0.2505

S21

0.0934

0.1761

0.1669

0.0894

0.1119

0.1322

0.3428

0.2895

0.3204

0.3485

0.2354

0.2016

0.2627

0.2451

S22

0.0926

0.1736

0.1398

0.1356

0.0599

0.1054

0.3352

0.2837

0.3361

0.3413

0.2309

0.1972

0.2575

0.239

S23

0.0789

0.1504

0.126

0.1223

0.0518

0.0676

0.2764

0.2303

0.2781

0.2813

0.2038

0.1972

0.2281

0.215

S31

0.1416

0.2518

0.0804

0.0742

0.0527

0.072

0.2295

0.2781

0.3074

0.3295

0.2729

0.2464

0.3008

0.228

S32

0.0808

0.1778

0.0677

0.0634

0.0446

0.0616

0.2873

0.1722

0.2684

0.2914

0.2213

0.2178

0.2455

0.1806

S33

0.1439

0.2288

0.1401

0.1345

0.0844

0.1299

0.3233

0.2929

0.2525

0.3512

0.2626

0.2306

0.2895

0.2215

S34

0.1816

0.3028

0.1814

0.174

0.1186

0.167

0.3943

0.3345

0.3964

0.3056

0.2766

0.2635

0.3314

0.2568

S41

0.1261

0.2158

0.0979

0.0927

0.0708

0.1132

0.273

0.2248

0.2749

0.2775

0.1531

0.1939

0.2013

0.163

S42

0.1074

0.2336

0.0713

0.0656

0.0459

0.0632

0.2805

0.2824

0.3072

0.2849

0.2106

0.1576

0.2832

0.1914

S43

0.1119

0.2674

0.0767

0.0701

0.0494

0.0669

0.2942

0.2945

0.3211

0.2982

0.2207

0.2634

0.201

0.2235

S44

0.1222

0.2879

0.1432

0.1362

0.1103

0.1065

0.3365

0.3069

0.3401

0.3407

0.2486

0.2649

0.3275

0.1823

4.2. Weighting Results of the Index System Based on DANP
The weight results derived from expert scoring and subsequent ANP calculations are presented in Table 4. To ensure the reliability and domain relevance of the results, four experts were invited to participate in the evaluation process. These experts possess extensive experience in the fields of Arctic environmental science, maritime safety, oil spill emergency response, and ecological risk assessment. Their judgments provide a comprehensive and interdisciplinary foundation for the weighting process.
As shown in Table 4, the weighting results derived from the DANP methodology reveal a clear hierarchy of importance among the Arctic oil spill risk assessment indicators. Ecological and Socioeconomic Vulnerability emerges as the most critical dimension, accounting for 44.83% of the total weight, with Economic Exposure (S43) representing the single most influential factor at 16.21% of the comprehensive weight. This dominance reflects the profound long-term consequences oil spills can have on Arctic communities and ecosystems, where economic activities are particularly vulnerable to marine environmental disruptions. The Emergency Response Capability category follows in significance with 29.17% of the total weight, where Response Time (S31) and Personnel Capacity (S33) stand out as crucial operational factors. Environmental Conditions, particularly Sea Ice Thickness (S22) and Concentration (S21), collectively contribute 19.99% to the overall risk assessment, underscoring the unique challenges posed by Arctic marine conditions. Interestingly, while Oil Spill Characteristics represent the smallest category at just 6.00%, the Ice Area Coverage indicator (S13) within this group carries substantial relative importance (76.07%), highlighting how Arctic-specific conditions fundamentally alter the risk profile of spills compared to temperate regions. The comprehensive weighting structure demonstrates how traditional oil spill risk factors are recontextualized in the Arctic environment, with greater emphasis placed on ecological fragility, response limitations, and socioeconomic consequences rather than just the physical characteristics of the spill itself. This weighting scheme provides a nuanced framework for prioritizing risk mitigation strategies that account for both the unique environmental conditions and the heightened vulnerability of Arctic social-ecological systems.
Figure 5. Weight histograms for the indicator impact factors.
Table 4. Weights of oil spill risk assessment indicators.

Control layer

Network layer

Relative weight

Indicator element

Relative weight

Comprehensive weight

Arctic oil spill hierarchical model

Oil spill characteristics (S₁)

0.0600

Petroleum type (S11)

0.0631

0.0038

Oil quantity (S12)

0.1763

0.0106

Ice area coverage (S13)

0.7607

0.0457

Environmental conditions (S₂)

0.1999

Sea ice concentration (S21)

0.3255

0.0651

Sea ice thickness (S22)

0.3451

0.0690

Wind speed (S23)

0.3294

0.0659

Emergency response capability (S₃)

0.2917

Response time (S31)

0.2902

0.0847

Reaction facility (S32)

0.1578

0.0460

Personnel capacity (S33)

0.2786

0.0813

Efficiency of oil spill recovery (S34)

0.2734

0.0798

Ecological and socio-economic vulnerability (S₄)

0.4483

Sensitive habitat (S41)

0.1233

0.0553

Human settlement density (S42)

0.2758

0.1236

Economic exposure (S43)

0.3615

0.1621

Shipping lane priority (S44)

0.2394

0.1073

For enhanced visualization of the weighting distribution among different evaluation indicators, the quantified weight values are graphically illustrated in Figure 5.
4.3. Fuzzy Evaluation Matrix and Risk Analysis
To further enhance the reliability and comprehensiveness of the risk evaluation, 30 evaluators were invited to assess each indicator based on predefined evaluation criteria. The panel of evaluators was composed of a diverse and interdisciplinary group, including academic researchers in Arctic environmental science, government officials from maritime regulatory agencies, emergency response practitioners, and industry experts in oil spill prevention and control. This ensured that the evaluations incorporated both theoretical insights and practical field experience. Risk factors related to an oil spill were evaluated, and their score distributions are detailed in Table 5. As a multi-level index system, the hierarchical structure of this paper requires multi-level fuzzy comprehensive evaluation to obtain the final evaluation results. According to the statistical results, the fuzzy evaluation of Network layer S is carried out, and its weight set and fuzzy matrix are WS and FS respectively.
Table 5. Questionnaire data statistics table.

Risk level

Extremely high risk

High risk

Medium risk

Low risk

Ultra-low risk

S11

5

15

6

2

2

S12

21

5

3

1

0

S13

15

6

6

2

1

S21

3

12

6

4

5

S22

19

7

4

0

0

S23

0

6

17

3

4

S31

13

11

3

2

1

S32

4

21

2

2

1

S33

11

15

1

1

2

S34

17

4

4

3

2

S41

13

8

7

2

0

S42

17

8

4

1

0

S43

20

4

3

2

1

S44

12

10

4

2

2

WS=0.060.19990.29170.4483(12)
FS=0.51430.21310.18240.06080.02950.25110.27660.29780.07630.09820.40390.39260.08530.06650.05170.54650.23440.13360.05750.028(13)
ZS=WS×FS =0.44390.28770.15520.06410.049(14)
score=0.44390.28770.15520.06410.049×54321=4.01(15)
Upon synthesizing the four sub-layers, the system-level assessment (S layer) indicates a high-risk classification. This suggests that even though certain components exhibit only moderate risk, oil spill incidents in the Arctic’s ice-covered context still present substantial potential hazards. The result underscores the cumulative and interactive effects of multiple factors in amplifying the overall risk level.
5. Discussion
This study systematically evaluated the contributions of 14 indicators to Arctic oil spill risk using a multi-criteria decision-making framework. Among the four major categories, spill characteristics, environmental conditions, emergency response capacity, and ecological and socioeconomic vulnerability, the indicators under ecological and socioeconomic vulnerability exhibited the highest weights. In particular, population density (S42), economic development level (S43), and navigation priority (S44) were identified as dominant contributors to overall risk. Similar to findings from tropical coastal regions, where oil spill vulnerability increases with higher population density and economic reliance on marine resources , Arctic coastal communities may experience similar patterns of socioeconomic sensitivity, particularly in areas with clustered settlements and limited adaptive capacity. In protected coastal regions, areas with intensive human activity, such as fishing and tourism, also exhibited heightened vulnerability . The proposed evaluation framework successfully distinguishes indicators with high systemic importance, thereby facilitating evidence-based prioritization of risk management efforts.
The DEMATEL analysis revealed distinct causal structures among indicators. Ice coverage rate and oil type were identified as high-causality, upstream variables, while indicators related to emergency response (e.g., equipment availability, response time) were more often classified as downstream or result-oriented. These findings underscore how environmental constraints can diminish the operational effectiveness of human interventions. Unlike traditional AHP approaches that assume indicator independence, the integrated DEMATEL-ANP method captures feedback loops and interdependencies among variables. Similar modeling strategies have been applied in maritime logistics and energy planning, validating the utility of ANP for structurally complex, risk-driven systems . In this study, the global weights derived from ANP reflect both the direct influence and network-mediated effects of each indicator, leading to a more robust and context-sensitive risk quantification.
A key strength of extending the DEMATEL-ANP framework with Fuzzy Comprehensive Evaluation (FCE) lies in its ability to translate expert-elicited indicators into graded risk levels, as validated by our findings. Without FCE, the model would only produce weights, not actionable risk categories. Applying FCE enabled us to categorize regions into discrete risk levels (e.g., medium, high), revealing how small variations in key indicators can shift overall risk tiers. This outcome aligns with existing applications in high-uncertainty environments. Similarly, flood risk assessment studies combining DEMATEL-ANP with fuzzy evaluation showed that traditional threshold-based categorizations often mask intermediate risk gradations vital for dynamic response planning . By enabling us to model risk as a continuum rather than a binary state, FCE enhances the interpretability and flexibility of our risk assessment. This proves especially useful in Arctic oil spill scenarios, where operational complexity and environmental variability demand nuanced, threshold-sensitive decision tools.
The integrated DEMATEL-ANP-FCE framework offers a transparent and adaptable approach to assessing oil spill risks in the Arctic, especially under conditions of limited data and logistical constraints. While the current model operates as a static evaluation mechanism, it effectively identifies high-priority risk factors. These insights support strategic investments in early monitoring, emergency preparedness, and targeted policy interventions. Previous research has validated the feasibility of combining ANP and fuzzy evaluation for Arctic oil spill scenarios, yielding structured and interpretable indicator hierarchies.
To improve practical applicability, future research could integrate spatial and temporal elements. Coupling the model with GIS-based visualization of risk zones would help identify spatial clusters of high-risk Arctic areas and inform targeted interventions . Implementing dynamic oil-spill trajectory models that incorporate ice movement, ocean currents, and wind-for example, using discrete-element simulations used in Arctic ice-condition scenarios-can significantly enhance temporal accuracy . Furthermore, incorporating Indigenous knowledge and stakeholder input through community maps and traditional-use data-such as those developed by the Alaska Ocean Observing System-enriches the risk framework with local experience and boosts its credibility and inclusiveness.
6. Conclusions
This research proposes a comprehensive assessment method for oil spill risk in the Arctic ice-covered regions. By establishing an evaluation framework encompassing four dimensions, including spill characteristics, environmental conditions, emergency response capacity, and ecological-socioeconomic vulnerability, an integrated model combining DANP (Decision-Making Trial and Evaluation Laboratory with Analytic Network Process) and FCE (Fuzzy Comprehensive Evaluation) is employed to achieve quantitative risk assessment. Expert consultation with Arctic maritime safety specialists and an extensive literature review were used to inform the model development. Findings reveal that ecological and socioeconomic vulnerability contributes most significantly to the overall risk, with indicators such as population density, economic development level, and shipping priority playing decisive roles. The research also reveals that the distinctive environmental factors in the Arctic, especially ice coverage, substantially restrict the efficacy of emergency responses. The FCE method computes the comprehensive membership degrees across different risk levels by multiplying the fuzzy evaluation matrix with the weight vector, and then derives the final risk score through weighted summation based on predefined level values. The resulting Arctic oil spill risk score is 4.01, corresponding to a high-risk classification. The proposed model offers a scientific basis for Arctic maritime safety management and policymaking. However, limitations remain, such as the lack of consideration for spatiotemporal dynamics. Future research will aim to address these gaps by integrating GIS-based spatial analysis with dynamic oil spill diffusion modeling, incorporating Indigenous ecological knowledge, and applying machine learning algorithms to enhance the model’s accuracy and applicability. These advancements are expected to promote more intelligent and precise risk management for Arctic oil spill scenarios.
Abbreviations

MCDM

Multi-criteria Decision-making

ANP

Analytic Network Process

AHP

Analytic Hierarchy Process

FCE

Fuzzy Comprehensive Evaluation

DEMATEL

Decision Making Trial and Evaluation Laboratory

DANP

DEMATEL-based Analytic Network Process

Acknowledgments
I would like to express my sincere gratitude to Professor Yaoming Wei for his academic guidance.
Author Contributions
Zihan Zeng is the sole author. The author read and approved the final manuscript.
Funding
This work is not supported by any external funding.
Data Availability Statement
The data is available from the corresponding author upon reasonable request.
Conflicts of Interest
The author declares no conflicts of interest.
References
[1] Makarova, I., Makarov, D., Buyvol, P., Barinov, A., Gubacheva, L., Mukhametdinov, E., Mavrin, V. Arctic Development in Connection with the Northern Sea Route: A Review of Ecological Risks and Ways to Avoid Them. Journal of Marine Science and Engineering. 2022, 10(10), 1415.
[2] Liu, M., Kronbak, J. The Potential Economic Viability of Using the Northern Sea Route (NSR) as an Alternative Route between Asia and Europe. Journal of Transport Geography. 2010, 18(3), 434-444.
[3] Helle, I., Mäkinen, J., Nevalainen, M., Afenyo, M., Vanhatalo, J. Impacts of Oil Spills on Arctic Marine Ecosystems: A Quantitative and Probabilistic Risk Assessment Perspective. Environmental Science & Technology. 2020, 54(4), 2112-2121.
[4] Jonsson, H., Sundt, R. C., Aas, E., Sanni, S. The Arctic is No Longer Put on Ice: Evaluation of Polar cod (Boreogadus saida) as a Monitoring Species of Oil Pollution in Cold Water. Marine Pollution Bulletin. 2010, 60(3), 390-395.
[5] Krapivin, V. F., Phillips, G. W. Application of a Global Model to the Study of Arctic Basin Pollution: Radionuclides, Heavy Metals and Oil Hydrocarbons. Environmental Modelling & Software. 2001, 16(1), 1-17.
[6] Afenyo, M., Ng, A. K. Y., Jiang, C. A Multiperiod Model for Assessing the Socioeconomic Impacts of Oil Spills During Arctic Shipping. Risk Analysis. 2022, 42(3), 614-633.
[7] Chen, X., Liu, S., Liu, R. W., Wu, H., Han, B., Zhao, J. Quantifying Arctic Oil Spilling Event Risk by Integrating an Analytic Network Process and a Fuzzy Comprehensive Evaluation Model. Ocean & Coastal Management. 2022, 228, 106326.
[8] Arzaghi, E., Abbassi, R., Garaniya, V., Binns, J., Khan, F. An Ecological Risk Assessment Model for Arctic Oil Spills from a Subsea Pipeline. Marine Pollution Bulletin. 2018, 135, 1117-1127.
[9] Bagheri, M., Zaiton Ibrahim, Z., Mansor, S., Manaf, L. A., Akhir, M. F., Talaat, W. I. A. W., Beiranvand Pour, A. Land-Use Suitability Assessment Using Delphi and Analytical Hierarchy Process (D-AHP) Hybrid Model for Coastal City Management: Kuala Terengganu, Peninsular Malaysia. ISPRS International Journal of Geo-Information. 2021, 10(9), 621.
[10] Lin, Z., Ayed, H., Bouallegue, B., Tomaskova, H., Jafarzadeh Ghoushchi, S., Haseli, G. An Integrated Mathematical Attitude Utilizing Fully Fuzzy BWM and Fuzzy WASPAS for Risk Evaluation in a SOFC. Mathematics. 2021, 9(18), 2328.
[11] Tabesh, M., Roozbahani, A., Hadigol, F., Ghaemi, E. Risk Assessment of Water Treatment Plants Using Fuzzy Fault Tree Analysis and Monte Carlo Simulation. Iranian Journal of Science and Technology, Transactions of Civil Engineering. 2022, 46(1), 643-658.
[12] Cai, S., Fan, J., Yang, W. Flooding Risk Assessment and Analysis Based on GIS and the TFN-AHP Method: A Case Study of Chongqing, China. Atmosphere. 2021, 12(5), 623.
[13] Kahrama, C., Gulbay, M., Kabak, O. Applications of fuzzy sets in industrial engineering: A topical classification. In Fuzzy Applications in Industrial Engineering, Kahraman, C., Ed., Springer: Berlin, Heidelberg, Germany; 2006, pp. 1-55.
[14] Wang, W. C., Lin, Y. H., Lin, C. L., Chung, C. H., Lee, M. T. DEMATEL-Based Model to Improve the Performance in a Matrix Organization. Expert Systems with Applications. 2012, 39(5), 4978-4986.
[15] Ou Yang, Y. P., Shieh, H. M., Tzeng, G. H. A VIKOR Technique Based on DEMATEL and ANP for Information Security Risk Control Assessment. Information Sciences. 2013, 232, 482-500.
[16] Afenyo, M., Jiang, C., Ng, A. K. Y. Climate Change and Arctic Shipping: A Method for Assessing the Impacts of Oil Spills in the Arctic. Transportation Research Part D: Transport and Environment. 2019, 77, 476-490.
[17] French-McCay, D. P., Tajalli-Bakhsh, T., Jayko, K., Spaulding, M, L., Li, Z. Validation of Oil Spill Transport and Fate Modeling in Arctic Ice. Arctic Science. 2018, 4(1), 71-97.
[18] Heshka, N. E., Ridenour, C., Saborimanesh, N., Xin, Q., Farooqi, H., Brydie, J. A Review of Oil Spill Research in Canadian Arctic Marine Environments. Marine Pollution Bulletin. 2024, 209(Pt B), 117275.
[19] Shieh, J. I., Wu, H. H., Huang, K. K. A DEMATEL Method in Identifying Key Success Factors of Hospital Service Quality. Knowledge-Based Systems. 2010, 23(3), 277-282.
[20] Görener, A. Comparing AHP and ANP: An Application of Strategic Decisions Making in a Manufacturing Company. International Journal of Business and Social Science. 2012, 3(11), 194-208.
[21] Chen, J. F., Hsieh, H. N., Do, Q. H. Evaluating Teaching Performance Based on Fuzzy AHP and Comprehensive Evaluation Approach. Applied Soft Computing. 2015, 28, 100-108.
[22] Câmara, S. F., Pinto, F. R., da Silva, F. R., de Oliveira Soares, M., De Paula, T. M. Socioeconomic Vulnerability of Communities on the Brazilian Coast to the Largest Oil Spill (2019-2020) in Tropical Oceans. Ocean & Coastal Management. 2021, 202, 105506.
[23] Da Silva, F. R., Schiavetti, A., Malhado, A. C. M., Ferreira, B., de Paula Sousa, C. V., Vieira, F. P., Pinto, F. R., de Souza, G. B. G., Olavo, G., dos Santos, J. B. Q., Campos-Silva, J. V., de Oliveira Júnior, J. G. C., Messias, L. T., Barbosa Filho, M. L. V., Accioly, M., Fabré, N. N., Abdallah, P. R., Lopes, P. F. M., de Kikuchi, R. K. P., Câmara, S. F., Batista, V., Soares, M. O. Oil Spill and Socioeconomic Vulnerability in Marine Protected Areas. Frontiers in Marine Science. 2022, 9, 859697.
[24] Lin, W. C. Maritime Environment Assessment and Management Using through Balanced Scorecard by Using DEMATEL and ANP Technique. International Journal of Environmental Research and Public Health. 2022, 19(5), 2873.
[25] Fazli, S., Kiani Mavi, R., Vosooghidizaji, M. Crude Oil Supply Chain Risk Management with DEMATEL-ANP. Operational Research. 2015, 15(3), 453-480.
[26] Ma, X., Wang, Y., Teng, Z., Li, S. Urban Flood Risk Assessment Based on DEMATEL-ANP Hybrid Fuzzy Evaluation and Hydrodynamic Model. Water. 2025, 17(10), 1494.
[27] Li, W., Liang, X., Lin, J., Guo, P., Ma, Q., Dong, Z., Liu, J., Song, Z., Wang, H. Numerical Simulation of Ship Oil Spill in arctic Icy Waters. Applied Sciences. 2020, 10(4), 1394.
Cite This Article
  • APA Style

    Zihan, Z. (2025). Quantitative Assessment of Ship Oil Spill Risk in Arctic Ice Area Based on MCDM Assessment Model. International Journal of Transportation Engineering and Technology, 11(3), 98-110. https://doi.org/10.11648/j.ijtet.20251103.11

    Copy | Download

    ACS Style

    Zihan, Z. Quantitative Assessment of Ship Oil Spill Risk in Arctic Ice Area Based on MCDM Assessment Model. Int. J. Transp. Eng. Technol. 2025, 11(3), 98-110. doi: 10.11648/j.ijtet.20251103.11

    Copy | Download

    AMA Style

    Zihan Z. Quantitative Assessment of Ship Oil Spill Risk in Arctic Ice Area Based on MCDM Assessment Model. Int J Transp Eng Technol. 2025;11(3):98-110. doi: 10.11648/j.ijtet.20251103.11

    Copy | Download

  • @article{10.11648/j.ijtet.20251103.11,
      author = {Zeng Zihan},
      title = {Quantitative Assessment of Ship Oil Spill Risk in Arctic Ice Area Based on MCDM Assessment Model
    },
      journal = {International Journal of Transportation Engineering and Technology},
      volume = {11},
      number = {3},
      pages = {98-110},
      doi = {10.11648/j.ijtet.20251103.11},
      url = {https://doi.org/10.11648/j.ijtet.20251103.11},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.ijtet.20251103.11},
      abstract = {Given the intensifying effects of climate change and the escalating human activities in the Arctic, such as shipping and resource extraction, the likelihood of oil spills in the region has increased substantially. The assessment of oil spill risks in Arctic ice-covered regions requires careful consideration of multiple interacting factors, including environmental parameters, socioeconomic vulnerabilities, and emergency response capabilities. There is an urgent need for a systematic assessment approach to support informed risk management decisions. This research aims to develop a multi-criteria decision-making framework to quantify the key drivers of Arctic oil spill risk and classify risk levels, thereby providing a basis for targeted intervention strategies. The intricate nature of Arctic ecosystems, combined with extreme weather conditions and sparse infrastructure, renders conventional risk models insufficient. In response to this challenge, our study employs sophisticated, integrative methodologies to more effectively capture the interrelationships between various risk factors. By integrating the DEMATEL-ANP method, this paper analyzes the causal relationships and relative weights among 14 indicators across four interrelated network layers. These methods not only clarify the causal chains between factors but also account for the intricate interrelations and dependencies within the system, offering a more comprehensive and adaptable decision-making framework. This integration leverages the strengths of both methods, enhancing the depth and accuracy of the analysis. A fuzzy comprehensive evaluation (FCE) approach is applied to convert expert judgments and derived weights into 5 categorical risk levels. The fuzzy logic component handles uncertainties inherent in expert elicitation, ensuring robustness despite data scarcity. The results indicate that ecological and socioeconomic vulnerability contribute most significantly to overall risk, while environmental factors such as sea ice coverage severely constrain emergency response effectiveness. The weighted FCE score identifies the Arctic oil spill scenario as high risk. The integrated DEMATEL-ANP-FCE framework effectively addresses the challenges of limited data availability in complex risk environments and provides a scientific foundation for Arctic oil spill monitoring, emergency planning, and policy formulation.},
     year = {2025}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Quantitative Assessment of Ship Oil Spill Risk in Arctic Ice Area Based on MCDM Assessment Model
    
    AU  - Zeng Zihan
    Y1  - 2025/08/26
    PY  - 2025
    N1  - https://doi.org/10.11648/j.ijtet.20251103.11
    DO  - 10.11648/j.ijtet.20251103.11
    T2  - International Journal of Transportation Engineering and Technology
    JF  - International Journal of Transportation Engineering and Technology
    JO  - International Journal of Transportation Engineering and Technology
    SP  - 98
    EP  - 110
    PB  - Science Publishing Group
    SN  - 2575-1751
    UR  - https://doi.org/10.11648/j.ijtet.20251103.11
    AB  - Given the intensifying effects of climate change and the escalating human activities in the Arctic, such as shipping and resource extraction, the likelihood of oil spills in the region has increased substantially. The assessment of oil spill risks in Arctic ice-covered regions requires careful consideration of multiple interacting factors, including environmental parameters, socioeconomic vulnerabilities, and emergency response capabilities. There is an urgent need for a systematic assessment approach to support informed risk management decisions. This research aims to develop a multi-criteria decision-making framework to quantify the key drivers of Arctic oil spill risk and classify risk levels, thereby providing a basis for targeted intervention strategies. The intricate nature of Arctic ecosystems, combined with extreme weather conditions and sparse infrastructure, renders conventional risk models insufficient. In response to this challenge, our study employs sophisticated, integrative methodologies to more effectively capture the interrelationships between various risk factors. By integrating the DEMATEL-ANP method, this paper analyzes the causal relationships and relative weights among 14 indicators across four interrelated network layers. These methods not only clarify the causal chains between factors but also account for the intricate interrelations and dependencies within the system, offering a more comprehensive and adaptable decision-making framework. This integration leverages the strengths of both methods, enhancing the depth and accuracy of the analysis. A fuzzy comprehensive evaluation (FCE) approach is applied to convert expert judgments and derived weights into 5 categorical risk levels. The fuzzy logic component handles uncertainties inherent in expert elicitation, ensuring robustness despite data scarcity. The results indicate that ecological and socioeconomic vulnerability contribute most significantly to overall risk, while environmental factors such as sea ice coverage severely constrain emergency response effectiveness. The weighted FCE score identifies the Arctic oil spill scenario as high risk. The integrated DEMATEL-ANP-FCE framework effectively addresses the challenges of limited data availability in complex risk environments and provides a scientific foundation for Arctic oil spill monitoring, emergency planning, and policy formulation.
    VL  - 11
    IS  - 3
    ER  - 

    Copy | Download

Author Information