Authors: Heitor Martinez-Grueira (Centro Militar de Farmacia de la Defensa, Ministry of Defense, Colmenar Viejo, Madrid, Spain), Rafael Asorey-Cacheda (Department of Information and Communications Technologies, Universidad Politécnica de Cartagena, Cartagena, Región de Murcia, Spain), Antonio-Javier Garcia-Sanchez (Department of Information and Communications Technologies, Universidad Politécnica de Cartagena, Cartagena, Región de Murcia, Spain), Joan Garcia-Haro (Department of Information and Communications Technologies, Universidad Politécnica de Cartagena, Cartagena, Región de Murcia, Spain)
Categories: Research Article
Source: PLOS ONE
Authors: Heitor Martinez-Grueira, Rafael Asorey-Cacheda, Antonio-Javier Garcia-Sanchez, Joan Garcia-Haro
The Naval Ship Code (NSC) was enacted in 2009 to standardize regulations for NATO member naval forces, and a study commissioned by the Spanish Navy General Staff (EMA) aimed to identify the factors that influence onboard personnel’s ability to move during an evacuation process. This study validated the soundness of the safety protocols implemented on navy vessels and highlighted the impact of certain characteristics of the embarked military contingent, such as body mass index, age, and seniority. It also found that such characteristics could act as distinctive factors among the embarked contingents in the evacuation of a military vessel. The study quantified the effect of these intervening characteristics, confirming the need for different displacement models for each of the study contingents to improve ship evacuation maneuvers. The findings of this study provide insights into the behavior of different embarked contingents during the evacuation process and can inform the development of more effective safety protocols for military naval operations. The starting hypothesis is that certain characteristics of the embarked military contingent have a decisive influence on their displacement capacity during the evacuation process. This hypothesis has been expanded in the sense that these same characteristics can act as differentiating elements among the embarked contingents evacuating a military vessel. It is possible to quantify the influence of these characteristics and implement a displacement model applicable in escape, evacuation, and rescue processes. Thus, the specific characteristics of a study contingent will be reflected in its displacement model. In this article we find that while members of the landing force (LF) show greater displacement capacity through a longitudinal corridor (around 10%), their ability to overcome other passage elements present on the study vessel is reduced (around 30%) compared to members of the vessel’s own crew.
Passenger and crew safety on civilian ships lacked proper technical or legislative consideration until well into the 20^th^ century [1–3]. On the other hand, this topic has always been present in the military environment due to its highly-qualified personnel, whose duties can involve risks to the ship and themselves. This singularity has put the Navy at the forefront of technical and technological innovation in the field of maritime safety.
The change in the Spanish Armed Forces (FAS) model at the beginning of this century caused naval safety to take a back seat until the NSC and its chapter VII dealing with Escape, Evacuation, and Rescue (EER) was issued. This subject is strengthened in the current version of the NSC [4], although it is still general and, in many cases, incompatible with the military environment as it collects data and references from civil regulation standards in the framework of the International Maritime Organization (IMO) [5, 6].
This is one of the reasons justifying the study carried out by the Spanish Military Operational Research Cabinet (GIMO), given the scarce or nonexistent application of personnel distribution or the standards in the NSC for military vessels based again on the IMO regulations for civilian ships. Among its particularities, it is worth mentioning the unrealistic ratios of gender and age for a military vessel, as well as its reference to passengers with reduced mobility or the distinction between passengers and crew [7].
Thus, it was necessary to undertake a rigorous study that would include a preliminary phase, during which a range of tests would be executed to gather essential data (affiliation, anthropometric, professional, behavioral, etc.) about the Navy’s personnel. The data collected and the conclusions of the study would consolidate a knowledge base from which to deploy new, more specific, and ambitious research lines [8].
This paper follows the logical assumption that different contingents present different displacement models, which implies acquiring new data once the obsolescence of the existing data had been confirmed. Thus, the same range of tests designed in the original study was executed by GIMO, focusing on the most common obstacles or passageways of a military vessel. This approach enriched the research with data obtained aboard an amphibious vessel using a Statistical Design of Experiments (SDE) [9], allowing for the sizing of sample contingents composed of the vessel’s crew (VC) and a landing force of Marine Infantry personnel.
By using the SDE methodology, it was possible to identify the population’s distinctive features (body mass index –BMI–, age, and seniority) that truly play a role in the study process and, at the same time, to determine the extent of that role. Then, a displacement model for each of the reference contingents was developed.
An ANOVA analysis was performed on the data obtained from the study contingents, where the null hypothesis of equal means was rejected, leading to the conclusion that there was not enough data to assume that both contingents belong to the same population. Thus, an independent displacement model for each contingent is reasonably justified. The model consists of functions that quantify the displacement capacity (time) of an individual belonging to a specific contingent based on their particular characteristics and the different obstacles on the military vessel studied during the evacuation process [10].
The rest of this paper is organized as follows. Section “Material and methods” describes the material and methods. Section “Data collection” presents the data collection process. Section “Dataset integration” shows how the dataset was integrated. Section “Comparing tests between populations” compares and analyzes the results achieved for each of the populations. In Section “Summary of the results”, a summary of the results is presented and discussed. Finally, Section “Conclusions” provides the closing remarks.
To begin the new research, an SDE was planned at an early stage [11] to determine the size of the necessary population sample and identify the characteristics that affect their ability to move during evacuation from the naval platform. We concluded that BMI, age, and seniority [12] decisively influence the movement of a contingent through the passageway elements tested [13, 14].
The current Military Personnel Law, enacted in 2007, foresaw the gradual aging of members of the Spanish Armed Forces who maintain a temporary professional relationship. Since age has a significant impact on the movement model of a contingent during an evacuation process, the data from previous investigations could be outdated.
Due to this risk, collecting updated data is necessary to determine whether the previously identified characteristics continue to influence the deployment of military personnel on board and to what extent (see example of Fig 1). To achieve this, the same range of tests planned in the original study would be executed to gain the necessary insight and implement a new displacement model for the reference personnel if necessary.
![Fig 1: Staff age by category and gender.Spanish Navy staff age by category and gender in 2017 [15].](pone.0316266.g001.jpg)
To conduct the range of tests, the Spanish Navy was asked to designate the amphibious assault ship “Galicia” (L-51) as the study vessel during the course of regular national maneuvers in the Gulf of Cadiz and all participating military personnel provided informed consent. This type of ship was used since the platform’s and its personnel’s characteristics needed to match the planned research. This would be challenging due to the vessel’s dimensions, the size of the reference population, and the heterogeneity of the personnel on board.
This type of ship performs very diverse functions requiring specific personnel to provide the necessary capabilities (landing force, flight units, medical personnel, staff, etc.). Furthermore, the deployment of this personnel serves a specific and temporary function that depends on the duration of their operations. This temporality limits the knowledge of the general layout of the ship and training and instruction in the vessel’s emergency procedures.
Among the personnel on board this type of ship, the landing force (LF), composed of members of the Marine Corps, stands out for its size, which sometimes doubles the size of the VC. Moreover, its temporality is limited to the duration of the planned amphibious operation, after which this force will leave the vessel. This temporality greatly restricts the knowledge this force has of the general layout of the ship and its emergency response protocols, compared to the VC’s level of training and instruction.
However, considering the LF as mere passengers would be an error, as it is a military contingent and, therefore, has basic instruction and training in the vessel’s emergency plans and procedures. In accordance with the EMA’s regulations [16] for this type of ship, the embarked contingents are mandatory participants in periodically scheduled training programs, with the objective of instructing them in knowledge of the vessel and its procedures, including the EER plan.
Drawing from extensive experience in data collection and analysis, a comprehensive set of tests was designed, taking into account prior investigations [17, 18]. Using the Statistical Design of Experiments (SDE), we determined the necessary sample size, thus eliminating the necessity to involve the entire population on board, particularly during active maneuvers. The SDE methodology allows for the formation of control groups (population samples) based on members’ specific characteristics. As a result, the influence of these characteristics on the study process is quantified, ensuring a representative study without the need for more resource-intensive methods, such as stratified random sampling [9].
To carry out an SDE, it is necessary to know the study population and its characteristics. With this information, the control groups that will participate in the tests to collect data (time measurements) are formed to make the test representative. This technique also allows us to determine the degree of involvement of each of the study factors and, thus, conclude which factors should be discarded if they have little or no influence.
During the tests, the time taken by each member of the control groups to pass through the following common passageways or overcome obstacles along the evacuation routes of the reference platform was



All these tests were conducted with normal lighting to determine the time taken by each individual under specific environmental conditions during their execution.
After collecting the measurements from the tests carried out by members of the sample populations (LF and VC) on board the reference vessel, sufficient data was provided to build the SDE project, generating a displacement model for each test. These data also went through a previous phase of data filtering and analysis.
Once the collected data were validated, the SDE project involved constructing a matrix of the multifactorial experimental plan to reference the experiments that make up each study test (corridor displacement, opening and closing of watertight doors, going through a hatch between decks) by combining the significance levels of their study factors (BMI, age, seniority).
As a preliminary step to the construction of the experimental plan matrix, it was necessary to build a table (see Table 1), where each experiment is identified by an inverted binary system that combines the assigned significance levels (high = 1 or low = −1) for each study factor (A for BMI, B for age, C for seniority). This combination determined the characteristics of the experiments and, therefore, the members who made up the control group (sample population) participating in each test.
The levels of significance (high, low) assigned to each study factor (BMI, age, and seniority) were defined by the mean value of each attribute reported by the sample population.
Based on the combination of factors and degree of significance, each experiment was provided with the data collected from the control group whose members verified the defined characteristics (example experiment A = 1, B = −1, C = 1; individuals with higher BMI, lower age, and more seniority than their respective means). Given the size of the participating sample population, it was determined that each control group would be composed of a total of 10 individuals.
Therefore, since each test had eight experiments and given the control groups were made up of 10 people, each response matrix (data matrix) consisted of eighty measurements (population sample).
On the basis of the data provided by the members of this sample population, the means attained for their study factors (BMI, age, and seniority) determined the levels of significance (see Table 2). These characteristics indicated an individual being assigned to a particular control group depending on the value of their attributes with respect to the factors (BMI, age, and seniority).
Thus, the measurements obtained for the individuals of each control group of the different experiments in each trial constituted the SDE’s response matrix. This technique allowed us to quantify the basic contributions and effects of the study factors to obtain the regression function that models the displacement capacity [19] for an LF subject based on their individual attributes (see Tables 3–5).
Based on the data provided by the members of this sample population, the averages obtained for their study factors (BMI, age, and seniority) determined the significance levels (see Table 6). These characteristics also established whether an individual was assigned to one control group or another, depending on the value of their attributes with respect to the factors (BMI, age, and seniority).
To achieve this, the measurements collected from the individuals in each control group of the various experiments conducted in each trial constituted the values in the response matrix of the SDE. Using this method, the basic contributions and effects of the study factors were quantified to obtain a regression function that models the displacement capacity [14] of a member of the crew based on their individual characteristics (see Tables 7–9).
The application of SDE allowed us to form control groups (population samples) that were representative of the original population and its inherent characteristics. These groups executed the planned tests to gather sufficient data to develop a displacement model. The captured data was then subjected to various tests [20, 21], refinement processes (Chauvenet’s criterion [22]), and a regression analysis based on the Shieffle equation to create the displacement function for a reference contingent. In this way, the displacement function estimates the time taken to overcome an obstacle or go through a passageway on the study platform, based on the intrinsic characteristics of an individual (BMI, age, and seniority) from the reference contingent during the EER process.
Given the disparate data collected from the different tests conducted on various passageways or obstacles on the platform, the same steps were taken for each contingent and each tested element or obstacle. This ensures a range of displacement models that quantifies the time it takes for a member of any contingent to move along the established evacuation routes on the platform. Therefore, an SDE was carried out for each test of the reference contingents, resulting in two separate lines of data analysis (LF and VC), to which refinement and regression techniques were applied.
After the data collected in the executed tests (longitudinal corridor and going through a hatch between decks and a watertight door) were cleaned and analyzed, regardless of the reference contingent, a displacement model was created for each test. Therefore, since the control groups are representative groups of their respective contingents (VC and LF), each contingent has a particular displacement model for each passage element or obstacle, modeled by regression functions that estimate the transit time of one of its members based on his/her particular characteristics (BMI, age, and seniority).
Even assuming that the tests conducted by each reference contingent were independent due to their distinct characteristics, a rigorous investigation required a new phase of data refinement and analysis supported by the following tests to accept or reject this
Having listed the tests used for the data cleansing and analysis of the data collected in the trials with the control groups consisting of members of the crew of the study vessel (see Table 10), the following observations can be
Based on these assumptions, the first two trials are immediately ruled out, and as a result of the ANOVA, the third trial is also discarded, thus rejecting the null hypothesis of equal means for the different trials.
According to these results and since the data from the different trials do not pass the ANOVA analysis, they do not validate their null hypothesis and, therefore, the available data do not allow us to assume that they have the same population mean, thus justifying the application of the independent displacement model for each trial constructed from the regression analysis using the SDE.
The same tests listed for the data cleansing and analysis phase of the trials with the control groups consisting of LF members were applied and the results (see Table 11) reveal the
As a result of these premises and the application of the ANOVA test, the null hypothesis of equal means for the different tests is rejected. Therefore, according to the results and since the tests do not have the same population mean, the development of a separate displacement model for each of them is justified, using the regression analysis process defined by the SDE.
Likewise, based on the results and given that the data from the different trials do not pass the ANOVA analysis, the null hypothesis is not validated, and, therefore, the available data do not allow us to assume that they have the same population mean. This supports the use of the independent displacement model for each trial, prepared from the regression analysis following the SDE.
In view of these results, the preparation of independent displacement models for reference passage or obstacle elements is justified. Furthermore, it is possible to verify whether the characteristics of the study contingents (VC vs. LF) or their respective levels of instruction and training in the general layout of the vessel and emergency procedures are substantially different among their members.
Thus, the data captured in the different trials were compared aggregately and independently of the contingent to which their control groups belong (VC and LF), obtaining the results that substantiate the following statements (see Table 12):
Therefore, the data from the different trials do not pass the ANOVA test and do not validate the null hypothesis, which means that the available data do not allow us to assume that they have the same population mean. This explains the application of an independent displacement model for each group, which is composed of the regression functions constructed for each trial in accordance with the SDE.
Based on the characteristics considered by the SDE for an individual (BMI, age, and seniority) and regardless of his/her originating group, to quantify displacement when going through a passage or overcoming an obstacle on the study platform, a displacement model is required that considers the characteristics of the individual in question for each passage or obstacle.
In accordance with the differences observed in instruction and training of the members of the study contingents (VC vs. LF), their displacement models will also be different and, therefore, the characteristics considered for an individual will have different weights as well, acting as differentiating elements in the study contingents.
The differences observed by the ANOVA between the reference contingents, although not very significant in the corridor displacement test (around 10%), are indeed quite remarkable in the tests involving opening and closing watertight doors and hatches between decks, where the average times employed can vary up to thirty percent (30%) of the average time used.
In summary, given the characteristics of the study, each reference contingent has a different displacement model composed of functions applicable to each passage element or obstacle (longitudinal corridor, watertight door, and hatch), as can be expected due to their uniqueness.
The following figures show our conclusions regarding the differences among the members of the study contingents in going through the passage or overcoming obstacles on the reference platform. While the LF members are faster going through the corridor (Fig 5), the VC members are quicker going down the hatch and opening and closing the watertight door (Figs 6 and 7).



All these differences are verified by the ANOVA test, although different speeds are not noticeable in Fig 5 (displacement through corridor) with a difference of around ten percent, but they are significant in Figs 6 and 7 (overcoming watertight hatch and deck hatch), where the differences are around thirty percent of the average times taken.
Thus, the displacement models are not compatible and each contingent member will have to use his/her corresponding displacement model according to their contingent and the passage or obstacle in question.
Data collected in the tests carried out (opening and closing a longitudinal corridor, overcoming a watertight door, and going down a hatch between decks) by the control groups of both reference contingents do not show equality in their population means. Thus, the application of the corresponding displacement function for the reference obstacle is justified to estimate the time taken by an individual with specific characteristics (BMI, age, and seniority) to overcome it.
Thus, the initial assumptions are confirmed, as the uniqueness of the study contingents results in significant differences depending on their level of knowledge about the general layout of the vessel and their training and instruction in the safety plans and procedures of the reference ship. All of this is reflected in their displacement model, which acts as a differentiating element in their ability to move during the execution of an EER process on board the study platform.
Considering the study’s findings, it is recommended that relevant jurisdictions adopt policies that promote data-driven approaches to decision making. Thus, ship evacuation plans should take into consideration and adjust to the demographic characteristics of the embarked population. In particular, these authorities should mandate regular data collection and validation to ensure the accuracy and relevance of data sets used in policy planning and plan development.