Abstract
Introduction
More than 100 million containers are shipped across the globe on containerships per year. According to containerised trade data, the number reached approximately 160.5 million containers in 2019. 1 Based on this, container transportation has become even more important for global maritime trade. However, significant container shipping disasters where hundreds of containers were lost in a single event have occurred in recent years. 2 The disastrous fires and explosions on Maersk Honam,3,4 MSC Flaminia,5,6 Hyundai Fortune6,7 and Hanjin Pennsylvania,6,8,9 hull fracture on MSC Napoli5,10,11 and hull girder fracture on Mol Comfort,5,12 and the breaking of Rena in two,13,14 collapsed and fallen overboard containers on MSC Zoe15,16 have caused the worst maritime environmental disasters in the last decade. Besides the loss of containers severely damaging the marine environment, tragically, some crew members have died because of the accidents.
Each operational activity carried out onboard ships includes risks due to the nature of the work. Therefore, identifying the risk factors and minimising them to an acceptable level is paramount to enhancing the safety level. 17 Human error, technical, mechanical, structural failure, and environmental factors are common causes of marine accident risk. 18 As the regulatory body, International Maritime Organization (IMO) emphasises that the human factor plays a crucial role in accidents. 19 The statistics show that more than 80% of shipping casualties are directly related to human error.20–22 Thereby, human error contribution should be the core point of the quantitative risk analysis (QRA) in maritime operations. A variety of approaches that focus on human error probability (HEP) quantifications have also been implemented in different industries such as offshore,23–27 aviation, 28 railway,29–32 nuclear power plants33–35 and mining. 36
The maritime industry seeks to reduce losses in the future. However, risk assessments carried out apart from the crew safety performance shall be insufficient in analysing the potential threats. At this point, some impact factors related to the task, individuals or working environment should also be considered while evaluating the HEPs. These relative factors, 37 called performance shaping factors (PSFs), are of paramount influence on human performance negatively or positively. 32
The SLIM technique considering HEP assessments has been used to determine the human error contribution to operational risks22,37–40 in the maritime transportation industry. In this study, a quantitative risk analysis is performed by considering the possible human errors in the container loading operation process. In this context, this paper proposed a hybrid approach by incorporating Fault Tree Analysis (FTA) and Interval type-2 fuzzy-based SLIM to evaluate the human contribution to risks and the criticality of the loading operation activities in a container terminal. To achieve this goal, the paper is structured as follows: The first part presents the motivation behind the study and basic literature review on significant container shipping disasters. Because of the substantial role of each method in the study, a brief literature review and the theoretical background of the methods are provided in section 2. Section 3 offers the integration of the proposed approach, while Section 4 illustrates the exemplificative application of the proposed approach to risk of container loss in maritime transportation. Findings and extended discussion are presented in section 5. Finally, the conclusion and research contribution to maritime transport is included in the last section.
Methods
The hybrid approach is proposed to determine the contribution of human error to the risks related with the most critical vulnerabilities in the operational processes. In this context, the SLIM estimates the HEPs whilst the FTA perform a comprehensive risk assessment. Since there is an ambiguity with the crisp value of probability, the IT2FS deals with vagueness and subjectivity in using experts’ judgements.39,41,42
IT2FS
The concept of a type-2 fuzzy set was first introduced by Zadeh 43 as an extension of the idea of a conventional fuzzy set called a type-1 fuzzy set (T1FS).41,44 A fuzzy set states the degree to which an element belongs to a set. In case it is not possible to determine the membership of an element in a set as 0 or 1, the type 1 or type 2 fuzzy sets are utilised. The membership grade for each element of the type-2 fuzzy set (T2FS) is a fuzzy set in [0,1]. On the other hand, a type-1 is a fuzzy set where a membership grade is a crisp number in [0,1].45,46 The basic principle behind systems is the same for both Type-1 and Type-2. However, T2FS can better express a higher degree of fuzziness and provides more various parameters than T1FS.45,47
An interval type-2 fuzzy set (IT2FS) is a special case of the generalised T2FS 41 in which the membership grade of every domain point is a crisp set whose domain is some interval contained in [0,1]. 44 Mendel 48 proposed the interval type-2 fuzzy set to describe an imprecise linguistic term, linguistically and quantitatively. 49 The data collected from the experts’ linguistic expressions are subjective and have limitations. At this point, the IT2FS can cope with complex conditions and reflects uncertainties better.44,50,51 IT2FS is rather adequate for utilising in real-case applications compared to generalised T2FS 52 and is commonly used in decision-making problems.53,54 The IT2FS is applied almost all problems by reason of their reduced computational effort and feasibility.39,44 Following a description of the T2FS and the IT2FS, the below equations present the mathematical operations’ definitions and step-by-step developments, respectively.
In addition, the type-2 fuzzy set
Where
where

The trapezoidal membership function of IT2FS.
A trapezoidal interval type-2 fuzzy set:
where
Mathematical operations using between two IT2FSs for further calculations are also as given below39,42,55:
For the addition operation:
For the subtraction operation:
For the multiplication operation:
For the arithmetic operations:
SLIM
The SLIM 56 was first introduced to estimate the probability of success of specific human actions in nuclear power plants. 57 The fundamental rationale of the SLIM is that the success likelihood of a task is based on the combined effects of a set of performance shaping factors (PSFs) which has a considerable influence on human performance. 58 The SLIM is a simple and flexible approach24,37,59 that makes use of domain expert judgement to select and weigh the PSFs according to their perceived contribution in a given task for estimating HEPs. 60 Accordingly, the core and crucial step is the formation of a committee of experts to generate the relevant data reliably. Following the quantification of PSFs, a Success Likelihood Index (SLI) is obtained utilising experts’ judgements for each action of the specific task.22,61 Subsequently, the SLI value is calibrated with the human error data to predict the HEP value. The main steps of the method are expressed as follows: (i) PSF derivation, (ii) PSF rating, (iii) PSF weighting, (iv) SLI determination and (v) HEP calculation
The below equation is utilised in the SLI determination process.
In the equation above,
Accordingly, the conversion of the SLIs to HEP values is achieved by a logarithmic relationship represented in equation (8).
In equation (8),
FTA
Fault Tree Analysis (FTA) is one of the most crucial logic and probabilistic techniques extensively utilised for reliability evaluation and probabilistic risk assessment of complex systems.62–64 The technique generates a mechanism for efficient system-level risk assessments. As a top-down and deductive failure analysis, 59 the technique identifies the sub-systems essential for the operation of a complex system. 65
Visualising a conventional fault tree comprises three major graphic symbols: events, logical gates and transfer symbols.66–68 Several sequential fault combinations that cause the undesired event called the ‘top event’ (TE) are depicted at different system levels. The TE is of enormous significance for the complex system due to cause catastrophic consequences for humans, commodity, and the environment.
69
Therefore, a fault tree is directly focused on the top event of the tree. In line with this purpose, the fault tree represents the logical interrelationships of basic events (BEs), which trigger the main event when they co-occur, and employs Boolean algebra rules. These rules are utilised to acquire one form of the fault tree, called the minimal cut set (MCS), that allows qualitative and quantitative assessments to be performed simply. The MCS specifies the system’s structural vulnerability.
69
The logical gates utilised to represent the relationships of events express the relationship type of the input events needed for the output event. The quantification of probabilities occurs according to the MCSs describing the relationships between BEs using ‘AND’ and ‘OR’ gates. Accordingly, the equation (9) is utilised to obtain the occurrence probability of the top event associated with the ‘AND’ gate, where
Associated with the ‘OR’ gate event, the equation (10) is utilised to acquire the top event’s occurrence probability:
The MCSs and overall failure probability of the top event are needed to calculate once the occurrence probabilities of BEs and IEs are gathered. The following equations are used for MCSs.70,71
The below equations are utilised to calculate the occurrence probability of TE.71,72
In the FTA technique, the FV-I (Fussell Vesely Importance Measure) method is utilised to ascertain the importance value of BEs and MCs constructing the TE.3,73 The following equation is used for the FV-I.
where
Integration of methodologies
The integration of methodologies for comprehensive risk analysis is provided in this section. The FTA is combined with the IT2FS-SLIM approach. In this context, Figure 2 illustrates the conceptual framework of the integrated method.

The conceptual framework of the integration.
Construction of a FT diagram
The first step of the hybrid approach is to construct a fault tree addressing the events’ interaction resulting in container loss. In the process, the FT is developed with references from containership accidents (which occurred last two decades) databases and investigation reports, as well as previous literature, and with the assistance of a group of marine experts. The experts familiar with containership cargo operations on board are involved as consultants due to the lack of failure probability data in the maritime industry. 69 Failures related to crew performance, environmental factors, technical and mechanical failures, and equipment functions are considered altogether for an effective FTA.
Data derivation under the IT2FS-SLIM approach
This section presents the data derivation process to evaluate human error contribution to the operational risks. The evaluation of HEPs in the maritime industry is regarded as onerous due to the scarcity of numerical data.69,75 The IT2FS-based SLIM approach can generate HEPs, particularly in cases where a lack of numerical data exists. In the SLIM, the marine experts provide professional judgement to bridge the gap. Under the hybrid approach, the probabilities for each human error-related basic event are acquired. Accordingly, the main steps of the process and their brief explanations are as follows.
Computing TE and MCSs failure probabilities
The IT2F-based SLIM approach to performing HEP assessments provides probabilistic outcomes for risk assessment in maritime transportation. The HEPs obtained by utilising the IT2F-SLIM steps are incorporated into the FT of container loss. Based on these outcomes, the failure probability of all BEs is calculated. Thereby, the overall likelihood of the top event (TE) and MCSs are computed for detailed risk analysis.
Model application: The case of container loss risk
This paper evaluates the container loss probability in containership cargo operations based on an FTA structure under IT2F-SLIM approach is developed to conduct a comprehensive risk analysis.
Problem statement
Several factors ranging from rough seas and heavy weather conditions to more catastrophic events such as collision, explosion, grounding, and hull damage can result in containers being lost at sea. 76 Apart from mentioned events, the likelihood of having other major hazard events such as listing, capsizing, structural fracture, and stack collapse leading to container loss is also significant during the cargo operations at the port period. In this study, containership loading operation is selected to illustrate the applicability of the proposed hybrid approach since it has potential risks for the safety of a container ship, its crew and cargo, shore-based workers, port facilities and the marine environment.
In accordance with non-mandatory and mandatory regulations issued by authorities, to avoid unwanted events, significant items must be checked by the watchkeeping team regularly. Ship stability values (GM, bending moment, torsion moment, drafts, trim and shearing force), stowage plan, visibility line, specific containers such as IMDG, reefers and, OH/OOG, lashing gear, lashings of containers and hatch covers demands great attention 77 throughout the containership cargo operation. In this context, crew performance plays a considerable part in risk analysis in identifying what errors lead to or contribute to the top event. However, whilst determining the human error contributions in the shipboard operations the human error should be treated as a combined outcome of some factors onboard the ship. Besides, failure can sometimes be beyond the crew’s control, although rare. Shipper-related issues (i.e. mis declared cargo and incorrectly/poor container packing), port-related issues (issues with hoisting cranes and port storage, poorly stacking containers and poor arrangement of weight distribution) and environmental conditions are also relevant factors in losing containers.
Analysis of respondents
Accident data sets, investigation reports, and empirical studies are the ideal, and key sources for human error prediction. 58 However, the data on maritime transportation is scarce or incomplete due to commercial reasons. 69 To meet this challenge, the SLIM utilises qualified experts’ judgements in the decision-making process to predict human errors. In this study, the appraisal of human error contribution to ship operations is evaluated with the participation of 10 qualified experts with substantial seagoing and working experience in containership transportation. Two out of these marine experts also have working experience as operation manager in container terminals. The following criteria were determined to form an expert group in this research; (i) minimum oceangoing Master licence, (ii) minimum 10 years of experience onboard container ship and (iii) physically participated in cargo handling operation on board container ship. At this point, Table 1 contains the profile details of marine experts. The marine experts make professional judgements expressing the PSFs impacts on each human error-related basic event utilising the linguistic statements of defined type-2 fuzzy sets.
Marine experts’ profile details.
Data derivation under the IT2FS-SLIM approach
This section summarises how the HEP data is derived to perform quantitative risk analysis. Since the loss of container operational risk is a concern, Table 2 illustrates the fundamental container handling tasks throughout the operation at a container terminal.
Task analysis for container handling operation.
In the study, seven PSFs used are captured from the recent study associated with containership handling operations. 38 Since it has paramount importance to derive appropriate PSFs rather than all PSFs, experience, stress, fatigue, training, time limitation, complexity and safety culture were specified by the Elicitation Review Team (ERT) as effective PSFs on crew performance during the loading operation. A brief description of each PSF included in the HEP assessment is given below, respectively.
Stress: Negative effect upon seafarer performance to complete the task correctly due to increased anxiety and pressure.
Experience: Familiarity with the task and knowledge.
Training: Expansion of knowledge, performance, and capability of seafarers by activities or actions organised by ship management.
Fatigue: Extreme tiredness caused by mental/physical workload or illness.
Time Limitation: Amount of time required for the seafarer to complete the relevant task.
Complexity: The measure of task difficulty identifies interrelated and interdependent task components.
Safety Culture: Both individual or group perceptions, attitudes and values that reflect ship management’s commitment to safety.
The further step is to determine the PSF rating for each task. The PSFs are rated by marine experts due to the lack of failure data in the shipping industry. The marine experts nominated a rate for each determined task according to the 1–9 linear scale, which reflects their relative judgements. The geometric means of ratings of 10 experts participating in the survey were obtained to simplify the calculation. Accordingly, Table 3 illustrates PSF rates for each task.
Geometric means of PSF ratings based on the marine experts’ evaluations.
After having determined PSFs, the weighting process is performed. The IT2Fs are used for the weighting process of PSFs since it is capable of handling inaccurate information in a logically correct manner. In this context, Table 4 demonstrates the IT2FSs number, and their membership functions related to the linguistic terms for determining the PSFs’ importance weight. 42 The next step is to calculate the defuzzified values of PSFs weights. In this context, linguistic variables are converted to the IT2FSs to quantitatively transform the judgements of marine experts. Once the average IT2Fs values are calculated, the defuzzification is conducted using equation (1). Table 5 shows IT2FS, crisp and normalised values of PSFs. 38
Lingusitic terms and their corresponding IT2FSs.
Calculated average IT2F values.
The HEP values are calculated using equations (7) and (8) where
Calculated HEP values for cargo handling operation.
Quantitative risk assessment for container loss
This section performs quantitative risk analysis for container loss by systematically predicting human error contributions to the operational risks. To achieve this purpose, the FT is constructed by reviewing accident investigation reports, literature, and marine experts’ judgement. In the constructed FTA, 30 basic events that will be effective in the realisation of the top event have been determined. At this point, the environmental conditions have been ignored since no environmental obstacle hinders the present real-time containership cargo operation, and the human error contribution was the focal point. Table 7 illustrates the TE, BE and IE for container loss risk in this context.
Fault tree events for the loss of containers.
Three main events cause the top event identified as container loss in the fault tree. These are the failures associated with cargo (IE01), failures associated with lashing (IE02) and failures associated with cargo handling (IE03). Having just one of these three main intermediate events is sufficient to cause container damage. Therefore, IE01, IE02, and IE03 are linked to the TE with the ‘OR’ gate. Accordingly, Figure 3 depicts the FT diagram for container loss during cargo handling operations in maritime transportation.

Structure of fault tree for the loss of container at port.
From the FT diagram and logic gates, TE (container loss) occurrence probability was calculated by applying equations (9) and (10), respectively. Based on the results, the occurrence probability of TE is found to be 5.54E-01. Accordingly, the MCSs, their occurrence probabilities, and the V–FIM list of MCSs are depicted in Table 8 (equations (11)–(13)).
Ranking of basic events according to Fussel-Vessely importance.
Findings and extended discussion
In light of the comprehensive risk assessment for container loss during the loading operation, the top event occurrence probability was calculated as 5.54E-01 which is a rather high. The obtained results show that 55 out of 100 cases may result in container loss due to the paramount contribution of human error during the loading operation. Since the fault tree structure is a graphic model representing the logical interrelationships of basic events, the possibility of each BE that includes human errors resulting in container loss was calculated to achieve TE occurrence probability. At this point, BE6 (1.38E-01), BE7 (1.20E-01) and BE21 (1.14E-01) with the highest HEP values were found to be the most contributory basic events increasing the risk of TE, respectively.
Further, the occurrence probabilities of the MCSs, the smallest combination of the BEs, were also calculated to identify the structural vulnerability of the system. Based on the results, BE4 (Misdeclaration/under declaration of the actual type/materials of Cargo) and BE5 (Misdeclaration/under declaration of the actual weight of the container) were the basic events that derive the most MCSs (four MCSs for each) among the others.
Lashing gear is a crucial item that needs to be checked by the watchkeeping team properly. Unlocked hatch cover cleats and loose lashings can cause a container stack to move and force on the adjacent stacks while the vessel is underway. Even worse, the forces on the adjacent stacks shall gradually increase and put the lashing equipment under additional load when the vessel rolls. Accordingly, any failure on lashing gear results in container loss due to stack collapse. However, the increasing effect of factors such as fatigue and limited time, makes the crew more vulnerable to errors, unavoidably.
One of the most significant goals of safe container handling is to minimise the occurrence probability of leaks, spills, or damage. Leakage is a crucial problem in the storage and transport of containers because it may corrode other stacked containers or produce toxic or inflammable fumes if they especially contain dangerous goods. Further, one of the essential parts of the planning is the confirmation that the permissible sequences of masses in stacks are not exceeded. Nevertheless, the weight of the leakage container becomes lighter as time goes by, resulting in container loss due to stack collapse. The primary cause of leakage is rough and inattentive container handling that causes structural damage during cargo operation, in general. Hence, each stowed container should be kept under strict control against any leakage throughout the handling process. At this point, safety culture, fatigue and training were determined as influential factors on human performance in the event of failure.
As for the misdeclared/undeclared cargo, the consequences can be catastrophic in some cases, an example being the disaster that resulted in the loss of the containership ‘Sea Elegance’ in 2003. 9 The report of the preliminary enquiry revealed that the fire and then explosion onboard originated in a container containing Calcium Hypochlorite that had not been declared. 78 Tragically, the disaster resulted in the death of one crew member and extensive cargo and vessel damage.
The disastrous explosion occurred in a cargo hold of the containership Hanjin Pennsylvania in 20026,8,9 is another unfortunate example of the significance of the subject. The containers filled with fireworks have been mis declared on the manifest. Thereby, the containers listed as having non-hazardous content were incorrectly stacked at the bottom of the hold and did not segregate as appropriate. The ship stayed afloat, but the disaster resulted in the death of two crew members and a substantial loss of cargo.
The consequences of underdeclared weights of containers led to a profound contribution to the catastrophic hull failure of MSC Napoli in 2007.5,10,11 Essentially, the vessel encountered rough seas that caused her to pitch heavily when on the passage in the English Channel. Following that, a catastrophic failure was suffered from her hull in the way of her engine room and then broke in two. The report by the MAIB (2008) stated a number of factors that contributed to the hull structure failure including the underdeclared weight of containers. All MSC Napoli’s containers were weighed again for investigation when beached in the UK, and the total weight of the 137 containers was 312 tonnes heavier than on the manifest. The load on the hull had increased by whipping effect and her hull already did not have sufficient buckling strength in way of the engine room. Although the detected non-compliance level was not evaluated as high, the report by the MAIB 79 identified it as concerning in the occurrence of this catastrophic event.
Conclusion
As a result of container losses from container ships, the maritime industry has taken the issue of safe stowage and securing of containers rather seriously because of the growing global concern over marine disasters. Since the tragic events caused the worst environmental disasters last two decades, the issue of container losses at ships is closely associated with environmental and economic aspects of the maritime transportation industry. At this point, identifying the causes of container losses can provide actionable solutions to reduce losses in future.
Despite the technological improvements, maritime operations remain dangerous for port facilities, vessels, the environment, and human life. Based on this, analysing the operational risk factors, and minimising the threats to an acceptable level is vital to enhance safety. Even though technical and mechanical failures are common causes increasing the risks, human error is found to be the most frequent and significant cause of marine accidents according to the conclusions drawn by the investigation reports.
This paper proposes a hybrid approach incorporating FTA and IT2FS-based SLIM to highlight the overriding importance of human-oriented failures in containership operations. In light of the extended risk analysis on real-time containership loading operation, the occurrence probability of the container loss was found to be 5.54E-01 which is considerably high. In the study, the importance of various factors was also identified as triggering human errors that should be addressed including ineffective safety culture, inadequate experience, fatigue, and limited time. Further, that the proposed approach can effectively be applied to identifying the operational vulnerabilities and critical human errors is concluded.
The fundamental limitation of the research is the scarcity of data. In the framework of the HEP assessment process that should contain both relevant data and real case studies, it is rather difficult to obtain empirical data in the maritime industry. Nevertheless, real data should be captured to validate the acquired results. A set of numerical simulations may also be carried out via risk analysis software in potential future research. This study is expected to provide qualitative and quantitative data on container transportation safety and insight into what measures may be necessary to decrease future losses by quantifying the potential failures in loading operations.
