Abstract
Objective
To investigate speech features, human–machine alarms, and operator–system interaction for the estimation of cognitive workload in full-scale realistic simulated scenarios.
Background
Theories and models of cognitive workload are critical for the design and evaluation of human–machine systems. Unfortunately, there are very few nonintrusive cognitive workload measures available for realistic dynamic human–machine interaction.
Method
The study was conducted in a full-scope control room research simulator of an advanced nuclear reactor. Six crews, each consisting of three operators, participated in 12 scenarios. The operators rated their workload every second minute. Machine learning algorithms were trained to estimate operators’ workload based on crew communication, operator–system interaction, and system alarms.
Results
Random Forest (RF) utilizing speech and system features achieved an accuracy of 67% on test data. Utilizing speech features only, the accuracy achieved was 63%. The most important speech features were pitch, amplitude, and articulation rate. A 61% accuracy was achieved when alarms and operator–system interaction features were used. The most important features were the number of alarms and amount of operator–system interaction. Accuracy for algorithms trained for each operator ranged from 39% to 98%, with an average of 72%. For a majority of analyses performed, RF and extreme gradient boosting (XGB) outperformed other algorithms.
Conclusion
The results demonstrate that the features investigated and machine learning models developed provide a potential for the dynamic nonintrusive measurement of cognitive workload.
Application
The approach presented can be developed for nonintrusive workload measurement in real-world human–machine applications, simulator-based training, and research.
Introduction
The measurement of cognitive workload during human–machine interaction is a topic of considerable interest for design, training, and evaluation (Gopher & Donchin, 1986; Hart & Staveland, 1988; International Organization for Standardization [ISO], 2006; Reinerman-Jones et al., 2019; Young et al., 2015). Postsession assessment of workload has been popular due to the ease of application (Grier, 2015; Reid & Nygren, 1988). However, there is a need to understand the dynamic interaction between workload and performance, and to investigate how operators and teams manage and adapt to cognitive workload (Hancock, 2017; Moray, 1988; Rusnock & Borghetti, 2018; Vicente, 1997). For example, it is challenging to detect and understand tipping points where cognitive workload relates to performance impairment (Chen et al., 2012; Hancock, 2017). Research using relatively simple tasks and highly controlled settings has established the potential of estimating cognitive workload using psychophysiological and behavioral features (e.g., Humphrey & Kramer, 1994; Yap et al., 2010; Yin et al., 2007), for example, by training machine learning algorithms (Wilson & Russell, 2003). While these studies demonstrate a potential and provide important background for unobtrusive and continuous measures, the extent to which such approaches can be used to estimate operators’ cognitive workload during realistic human–machine work is unclear. Investigating the measurement of workload during realistic work is especially challenging because the team’s strategy, interactive behavior, and timing of critical actions might strongly influence the overall task load, the distribution of load within the team, and the data available for workload estimation (Vukovic et al., 2019).
Dynamic Estimation of Cognitive Workload
What approaches might be suitable for realistic human–machine studies and field applications? Studies on real-time workload measurement include the use of eye behavior (Ahlstrom & Friedman-Berg, 2006; Schulz et al., 2011), heart rate, and heart rate variability (Kennedy & Scholey, 2000), brain activity (Humphrey & Kramer, 1994; Wilson & Russell, 2003), galvanic skin response (Shi et al., 2007), electroencephalographic signals, and breathing rate (Ferreira, et al., 2014). Although recent research has investigated less intrusive methods for some of these metrics (e.g., Rahman et al., 2020), these methods tend to require equipment that may interfere with the operator’s tasks and controlled conditions not easily implemented for realistic human–machine work. Dynamic workload measurement approaches using speech features (Brenner, et al., 1985; Lively et al., 1993; Vukovic et al., 2019), system variables (Rahman et al., 2020), and user interaction (Chen et al., 2012) can be applied in a nonintrusive manner, and may be more suitable for workload measurement in realistic settings. We have therefore focused on the latter type of features in our study.
Several studies suggest that fundamental frequency, vocal intensity, and speaking rate are among the speech features most sensitive to cognitive load and are therefore appropriate candidates for workload estimation (Brenner et al., 1994; Huttunen, Keränen, Väyrynen, et al., 2011, Huttunen, Keränen, Pääkkönen, et al., 2011; Wittels et al., 2002). We have, therefore, included these features in our study. The fundamental frequency (f0), defined as the rate at which the vocal folds open and close (Lieberman, 2007), may be influenced by stress and level of arousal (Brenner et al., 1994; Lively et al., 1993; Luig & Sontacchi, 2010). Speaking rate has been associated with an acceleration of cognitive and motor processes, often emerging under stress (Brenner et al., 1994), and increased vocal intensity has been associated with stress (Brenner et al., 1994; Luig & Sontacchi, 2010). Also, structural aspects, such as time spent pausing, may reflect an increased level in cognitive processing (Esposito et al., 2007; Lin et al., 2013).
Research on user–system interaction has found correlations between workload and several behavioral features related to interaction with the user interface, for example, frequency of pauses and changes of direction in mouse movements (Arshad et al., 2013; Lin et al., 2013), number of air traffic controllers’ flight clearances (Tobaruela et al., 2014), and mouse and keyboard activities (Pimenta et al., 2016). Braarud and Kirwan (2011) reported information intensity in terms of number of alarms as a factor of human–machine task load. Alarms require operators’ attention and aid the team in initiating tasks in response (Wu & Li, 2013). To our knowledge, no research has been published on using alarms as an indicator of workload related to operators’ main tasks, although the related topic of alarm fatigue due to false or nuisance alarms has been researched within hospital settings (e.g., Cvach, 2012; Ruskin & Hueske-Kraus, 2015).
A question with no clear answer in the literature is the unique contribution from different types of workload indicators and the accuracy that can be achieved from combining these types. In general, combining different measurement methods may compensate for the limitations of each individual method (Cegarra & Chevalier, 2008). For example, Chen et al. (2012) combined speech features and user interaction to predict workload. Solovey, et al. (2014) found that including physiological data gave higher algorithm classification accuracy than vehicle sensor data alone. On the contrary, McDonald et al. (2020) reported that combining physiological data and driver behavior data did not improve classification accuracy above the accuracy achieved with driver behavior data alone.
Application of Machine Learning Algorithms
Training machine learning algorithms for dynamic workload measurement have been reported in the literature (Vukovic et al., 2019; Wilson & Russell, 2003). However, the literature provides limited support for selecting adequate machine learning algorithms for cognitive workload estimation; hence, we aimed to include a set of commonly applied algorithms to ensure that an adequate algorithm would be included in the study. Machine learning algorithms reported include k-nearest neighbor (KNN; Solovey, et al., 2014), Naïve Bayes (NB; Solovey, et al., 2014), neural networks (NNs; Wilson & Russell, 2003), Random Forest (RF; McDonald et al., 2020; Ragab et al., 2014), support vector machines (SVMs; McDonald et al., 2020; Vukovic et al., 2019), and extreme gradient boosting (XGB; Momeni et al., 2019). Results from comparing algorithms’ performance are scattered, but there is some evidence that RF is the highest performing algorithm for driver distraction classification (McDonald et al., 2020; Ragab et al., 2014). Other results include Momeni et al. (2019), who found XGB to outperform a set of algorithms in their classification of cognitive workload; Solovey, et al. (2014) found decision trees, logistic regression, multilayer perceptrons, and NB to be equally good algorithms for classifying driver workload, while KNN performed significantly worse.
An important question for algorithm performance is the relevance of the training data. Individual differences have been identified in the experience of cognitive workload (Braarud, 2020; Guastello et al., 2015; Matthews et al., 2015), and there are individual differences regarding most speech features. Algorithm training on data from the person to estimate workload for is referred to as person-dependent, while algorithm training on data from other individuals is referred to as person-independent (Chen et al., 2012). Earlier research suggests that person-dependent models should outperform person-independent models (van Segbroeck et al., 2014; Vukovic et al., 2019). In addition, in human–machine systems, the team is a relevant unit of analysis (Hutchins, 1995) and may also be a relevant unit for training algorithms. There might be patterns of speech within the team, and activities by one operator might impact the other team members’ workload.
Study Purpose
Existing literature provides limited insights into how to continuously measure cognitive workload in a way that is applicable to realistic human–machine work. The extent to which the techniques studied in relatively controlled settings can provide tenable results is unclear, as is the answer to the question “which are the most suitable machine learning techniques (algorithms) for this kind of estimation?” Furthermore, there is uncertainty about the combined effect of different types of features and their relative importance. To address this gap in the literature, this study investigated the potential of using nonintrusive features from speech, operator–system interaction, and system alarms for cognitive workload estimation in a full-scope research simulator. Teams of licensed control room operators, guided by authentic operating procedures, performed in realistic accident scenarios that were relatively demanding. The operators’ self-labeled cognitive workloads were used for supervised training of machine learning algorithms. We investigated the performance of six commonly applied algorithms in order to identify the best performing model. We applied the best performing model to provide insights into the approach’s efficiency for cognitive workload estimation, and we investigated the combined and relative importance of the different types of features. In addition, the study explored three different approaches to training data for machine learning estimation of cognitive workload.
Method
Participants
The Organisation for Economic Co-operation and Development (OECD) Halden Reactor Project (OECD NEA, 2020) established an agreement with a plant with a reactor design similar to the one simulated at the research site. Six control room crews from the plant volunteered for participation in the study, which lasted 3 working days for each crew. Each crew consisted of a supervisor, a reactor operator, and a turbine operator, all being male licensed operators. Their mean age was 39.2 years (SD = 11.9), ranging from 26 to 62 years, and their mean control room work experience was 10.6 years (SD = 10.8), ranging from 1 to 35 years. Five crews consisted of team members that work as a crew at their home plant, while one crew was assembled from different home plant crews. The study was reviewed and approved by the Halden Reactor Project Human Studies Review Committee, and was performed according to Halden Reactor Project’s human participant protection procedures.
Simulator and Scenarios
The study was performed at the Institute for Energy Technology’s (IFE’s) Halden Man-Machine Laboratory (HAMMLAB; Institute for Energy Technology [IFE], 2019), which is a full-scope research simulator based on an advanced nuclear power plant. The simulator has a fully computerized human–machine interface. Figure 1 shows the control room layout. The shift supervisor workstation is at the back, the reactor operator workstation is to the left, and the turbine operator workstation to the right. The large screen display at the front provides a plant overview.

The Halden Man–Machine Laboratory control room.
The simulation and test scenarios were controlled from an experimenters’ gallery adjacent to the control room. Two subject-matter experts implemented and controlled the scenarios according to a scenario script and role-played the plant staff external to the control room, including field operators, technical support, and plant management. The control room team’s orders to field operators were thus implemented in the simulator by the subject-matter experts during the scenario. The control room team communicated with external staff by telephone, as they would at their home plant. Both subject-matter experts were former shift supervisors and training instructors with extensive knowledge of the simulated plant and the study’s specific scenarios. The subject-matter experts were involved in the scenario design as members of the research team.
Each team participated in 12 scenarios (Supplemental Material provides an overview of each of the 12 scenarios). Scenarios 1 to 8 were designed to last for 15–25 min. Due to a second topic of the study regarding human factors validation methods (Braarud et al., 2020), Scenarios 9 to 12 were designed to last 30–40 min. Even though the duration of scenarios varied, the scenarios comprised similar types of operator behavior and thus provided a consistent set of data for the investigation of cognitive workload estimation. The general scenario structure consisted of initial work, comprising normal operation or periodic testing including minor plant system failures and aggravation deviations leading to actuation of plant safety functions (reactor protection systems). The scenarios included events commonly included in nuclear power plant safety analysis, such as loss of reactor coolant accidents, loss of all offsite power, loss of turbine condenser, and loss of main feedwater. The team applied event-based procedures to acquire an overview of the plant’s safety status and to develop a basis for selecting a strategy to mitigate the situation and control plant safety.
In order to make parts of scenarios complicated and cognitively demanding, malfunctions were simulated for a number of safety system components. In these scenarios, safety components did not take the correct position after automatic actuation. While individual component failures like these can occur several times during a plant’s lifetime, the combination of the main event and multiple safety component malfunctions makes the specific scenarios unlikely. Nevertheless, the general type of scenario was not unfamiliar to the operators because similar scenarios are included in their annual retraining. The subject-matter experts stated that the workload imposed by the scenarios were within the range experienced during simulator training. The teams were instructed to use work routines and procedures as they ordinarily would in simulator training at home. This included routines for communication, teamwork, and interacting with plant staff external to the control room.
Operator Workload Rating
Each operator was instructed to self-assess his cognitive load every second minute using a six-point rating scale. The scale applied in the study was developed specifically for the evaluation of control room work and was based on the concepts of spare capacity (Gopher & Donchin, 1986; Moray, 1967; Wickens & Tsang, 2015) and adaption to task demand (Hancock & Warm, 1989; Hockey, 1997; Moray, 1988; Rusnock & Borghetti, 2018; Vicente, 1997). A 2-min interval for periodic self-assessment was chosen based on the subject-matter experts’ evaluation of the structure and timing of the control room work. Also, Jordan (1992) suggested 2-min intervals for similar types of application. Table 1 shows a description of the scale values. This description was available at each control room workstation as a part of the operators’ instructions for workload rating.
Description of Workload Scale Levels
The rating scale was implemented on tablet devices positioned on the participants’ desks. The tablets provided an audio signal every second minute for the participants to rate workload.
In addition to the workload rating during the scenario, the NASA Task Load Index (Hart & Staveland, 1988) was administered postsession for each scenario. The study utilized an unweighted version, often referred to as the raw TLX (Byers et al., 1989; Nygren, 1991).
Session Recording
The simulator sessions were recorded using the laboratory’s data collection systems. Each operator wore a headset with a microphone. The recordings included separate audio tracks for each operator, the simulated plant’s scenario history (alarms, process parameters, and process events), operator process commands, navigation and interfaces accessed by the operator, video of each operator workstation, and an overview video of the control room.
Speech Feature Extraction
Speech features were extracted from the session audio recordings using the Praat tool (Boersma & van Heuven, 2001). The Praat library was run from Python using the Parselmouth application programming interface (Jadoul et al., 2018). Praat was chosen as the accuracy is on par with other popular acoustic analysis software packages (Burris et al., 2014) and the application programming interface enables convenient use of the extraction algorithms from Python. The parameter settings used were: silence threshold, −20 dB; minimum dip between peaks, 2 dB; minimum pause duration, 0.3 s; and time step, 0.01 s. A small sample of audio recordings was manually checked to verify that the algorithm’s count of syllables and pauses was correct, and that noise or other disturbances were not influencing the feature extraction. The fundamental frequency range was set to exclude periods without speech when calculating mean values. Table 2 presents the features extracted. An aggregated, or mean, value for f0 was calculated for each operator for every 2-min segment of the scenario corresponding to the workload rating times.
Speech Features Extracted
Plant Alarms and Operator Activity Extraction
The identification of plant alarm and operator activity features was discussed with the two subject-matter experts on control room work and based on the available variables in the simulator logs as well as the design of the human–system interface (HSI). The initial discussions considered using the behavior of plant process parameters and the status of safety functions as task demand indicators. However, it was noted that plant alarms already cover status and encapsulate easily accessible variables covering the reason for requiring the operator’s attention. Alarms were therefore chosen as indicators of incoming task load. We identified three categories of features that we assumed to be related to the operators’ cognitive workload: incoming alarms and the acknowledgment of these, interface navigation, and commands directly manipulating plant systems. Most of the simulator’s alarms were specified as either belonging to the scope of the Reactor Operator or the scope of the Turbine Operator, while the remaining alarms were not specified for a specific operator. For our purpose, the former was defined as “operator alarms” and the latter as “crew alarms.” Regarding HSI activity, the calculation included all types of operator navigation recorded by the simulator’s logging system. Therefore, all plant alarms and operator activity features were extracted from the simulator logs and calculated as counts of occurrence during each scenario segment. Table 3 presents the features extracted from the simulator session recordings. A value was calculated for each operator for every 2-min scenario segment.
System Features Extracted
Data Pre-Processing, Algorithm Training, and Evaluation
For the classification algorithms, the six-point workload scale was aggregated to three levels, representing low, medium, and high cognitive workload. Currently, much of the reported research on machine learning classification of cognitive workload seems to be limited to two (Vukovic et al., 2019) or three (Charfuelan & Kruijff, 2013; McDonald et al., 2020; Yap et al., 2010; Yin et al., 2007) levels. A six-level classification is more challenging than a three-level classification. To be able to compare the study’s results with results reported in the literature, we chose to aggregate the scale to three levels for the machine learning classification. To produce labeled speech and plant system data, each 2-min period of speech and plant system development prior to the workload rating was labeled with the corresponding self-assessed workload.
All algorithms were trained in R (R Core Team, 2019) using the “caret” package (Kuhn, 2008; Kuhn et al., 2019) as the framework. The R caret package was used to split data into training and test sets, normalizing all features and controlling the algorithm training. Table 4 presents the six machine learning algorithms that were selected based on approaches commonly reported in earlier literature and our awareness of commonly applied classification algorithms within data sciences. The NN algorithm applied was limited to one hidden layer. The table also shows the parameters tuned for each algorithm. The parameters are described in the references provided in Table 4.
Machine Learning Algorithms Applied
The data were partitioned into algorithm training and algorithm test sets by the caret “createDataPartitionfunction.” This function performed random partitioning within the levels of the outcome variable, balancing the class distributions of the operators’ workload rating within the training and test sets. The split was set to 80% training data and 20% test data, except for the training and evaluation of person-dependent models. Person-dependent models utilized the data for a given operator, so in order to obtain a reasonable number of test cases, the split was set to 70% training and 30% test data. The data were centered and scaled using caret functionality. To identify the best model, models were trained using 10-fold cross-validation (James et al., 2013) repeated five times with accuracy as the training metric. However, that approach would inflate the degrees of freedom and increase type I error if used to compare models statistically. Therefore, for model comparison based on model variability, models were trained using twofold cross-validation repeated five times (Dietterich, 1998). Models were tested by estimating cognitive workload from the withheld test data. A model’s accuracy was determined by comparing the estimated classification with the operators’ actual workload classification. Accuracy ranged from 0% to 100% correct classification.
Results
As described in the previous sections, we used different algorithm training approaches for identifying the best model and for statistically comparing model performance. Please note, therefore, that the results from these two approaches are not directly comparable.
Overview of Operators’ Cognitive Workload Classification
The operator’s self-rating of cognitive workload during the 12 scenarios resulted in a total of 2780 labeled 2-min scenario segments: 1202 segments were labeled low, 1406 were medium, and 172 were high cognitive workload. The occurrences of workload categories were unbalanced with relatively few high workload segments.
Table 5 shows that the operators rated scenarios segments quite differently. For example, the Crew 1 Supervisor rated 20% of the segments as low workload and 70% of segments as medium, while the Crew 6 Supervisor rated 97% of segments as low and 3% as medium workload. The observed differences might result both from individual differences in the experience of cognitive workload and from the unique progression of the scenarios due to the crew’s strategy and timing of actions.
Counts of Self-Rated Cognitive Workload Segments
Note. Op. = Operator.
To compare the operators’ classification of workload during the scenario with the postsession NASA TLX rating, we averaged the 2-min workload ratings to produce a scenario average score. For a simple comparison, we rescaled workload scores by dividing the score by the respective scale length, providing a proportional score ranging between 1/scale length to 1. For the TLX rating, we calculated the average of all six TLX items. Figure 2 shows the TLX average and the online average per scenario.

Operators’ TLX average and online workload rating average. Mean and standard deviation. TLX = Task Load Index.
Figure 2 shows that the mean ratings for the TLX and online ratings follow a similar pattern from Scenarios 1 to 8. The average online ratings were slightly above or at a similar level as the mean TLX rating. The standard deviation of both the TLX ratings and the online rating reveals substantial individual differences. This corresponds to the literature regarding individual difference in the experience of cognitive workload (Braarud, 2020; Guastello et al., 2015; Matthews et al., 2015). The correlation between the TLX average and the average online workload rating was r = .55 (p < .001). The NASA TLX is a postsession evaluation based on the operator’s perception and memory of the whole scenario while the online rating focuses on the immediate perception of 2-min scenario segments. Some level of disagreement between these two measures must therefore be expected. Based on the similar pattern of means and individual variation, along with the moderate correlation between the two measures, we find that the online rating measures captured an element of cognitive workload in the control room.
Overview of Extracted Features
This section provides the features’ distribution across the operator-rated cognitive workload levels. For brevity, we selected the five most important speech features and three most important system features identified in the “Variable Importance” section later. The uneven selection, five versus three, was based on the study including 14 speech features but only six system features. Figure 3 shows the feature distribution for the reactor operator, supervisor, and turbine operator positions, respectively. The features were scaled with a minimum-maximum range suitable for non-normal distribution.

Dispersion of eight features aggregated for each operator position.
Figure 3 suggests several hypotheses regarding how the features relate to levels of cognitive workload; however, the development of such hypotheses is outside the scope of this paper. Two caveats apply regarding the interpretation of the dispersions presented in Figure 3. The operators classified few segments as high cognitive workload and consequently the dispersion of features for high cognitive workload is relatively more uncertain than for low and medium cognitive workload. Due to substantial individual differences, the development of hypotheses about features should include an investigation of person-level data beyond the aggregated visualizations per position presented here.
Classification of Workload: Identification of Best Performing Algorithm
The six machine learning algorithms were trained using each of the three input data types (speech, system data, or combined), resulting in 18 model types. For each model type, the best model was selected. The accuracy of a given model was determined by using it to classify the 20% data withheld. Figure 4 shows the models’ accuracy grouped by data type along with 95% confidence intervals. The 95% interval is a binomial proportion confidence interval, which is based on the accuracy resulting from the binomial proportion of hits (correct prediction; Brown et al., 2001).

Machine learning models’ accuracy on test data. Black lines represent 95% confidence intervals.
The differences in accuracy on test data were small between the models trained on speech data and models trained on system data. For the best performing algorithm (RF), accuracy achieved with speech data were 63% while accuracy achieved with system data were 61%. Training the algorithms on both speech and system data (combined) resulted in increased accuracy; however ,the increase was small. The accuracy of the RF algorithm, combined data, was 67%. The models’ accuracies can be compared to a random classifier resulting in an accuracy of 44%.
Comparing Model Performance Within Each Data Type
For a closer investigation of algorithm differences, we analyzed the variability in accuracy during training of the models (Dietterich, 1998). Based on Dietterich (1998), the models were developed by splitting the training data in two folds and repeating this procedure five times—resulting in 10 models with an accompanying calculated training accuracy. A t-test was performed for pairwise comparison of models with a Bonferroni correction of p values (Kuhn et al., 2019). The Bonferroni correction of a .05 criterion within each data type resulted in a .0033 criterion. Table 6 shows the results.
Pairwise Difference Between Models’ Accuracy and Statistical Significance
Note. *p < .0033. KNN = K-nearest neighbors; NB = Naïve Bayes; NN = neural networks; RF = Random Forest; SVML= SVM with Linear kernel; XGB = extreme gradient boosting.
Overall, the RF and XGB algorithms performed the best. For speech data, RF outperformed all models except for XGB, while XGB outperformed all models except for RF and NN. For system data, the best performing algorithms were RF, XGB, and NN. The were no significant differences for pairwise differences involving RF, XGB, and NN. For combined data, both RF and XGB outperformed all algorithms.
Confusion Matrix Across Data Sources
To provide further insight into model accuracy, we developed a confusion matrix with algorithm estimation and the operators’ classification—incorporating the accuracy according to the data sources (McDonald et al., 2020). For brevity, we focus here on the RF algorithm only.
Figure 5 shows poor performance on the classification of high cognitive workload. This challenge persisted for both data sources as well as for the combined data. Nevertheless, although the overall accuracy of the speech and system variables RF algorithms was almost identical, the speech data algorithm classified a higher proportion of high workload correctly, while the system variables algorithm more frequently misclassified high cognitive workload as low. Figure 5 also shows that the algorithm trained on all data improved the classification of both low and medium workload in comparison with the algorithms trained only on speech or system data; however, the combined data algorithm performed worse than the speech data algorithm at classifying high cognitive workload. Although high workload classification was challenging, the RF algorithm did not misclassify the operators’ low or medium workload as high workload. An investigation of the result of each of the six algorithms showed that they were all weaker at classifying high cognitive workload than lower workload.

Confusion matrix: data sources.
Variable Importance
Estimation of feature importance provides information on the relative importance of the variables for the models. Feature importance was estimated using the mean decrease in accuracy when removing a variable from the model (Kuhn et al., 2019). Figures 6 and 7 show the RF model features’ relative importance for the 10 most important speech features and the system features respectively in estimating cognitive workload. The most important feature is set to 100% and the other features importance are scaled accordingly.

Relative importance of speech features for estimation of cognitive workload measured by mean decrease in accuracy.

Relative importance of operator–system interaction features for estimation of cognitive workload measured by mean decrease in accuracy.
The results show that the operator’s pitch (f0), amplitude, and articulation rate are the most important features for estimating cognitive workload. Both mean level and variability were important for pitch and amplitude. The top three features—pitch, articulation rate, and amplitude—suggest that raising the pitch, speaking faster, and speaking louder are relevant indicators of cognitive workload.
The system data included only six features. The results show that plant system alarms and human–system interaction are the most important features. Also, acknowledging operator alarms contributed substantially. The low contribution of commands probably relates to the relatively few plant system and plant object commands actuated by the operators during the scenarios. The advanced and highly automated plant simulated implies a high proportion of cognitive work and that relatively few, but important, commands are executed.
When combining speech and system data for algorithm training, the results show that the five most important variables were all speech features. The top five features were identical to those identified as most important for the speech only data. The most important system variable was operator alarms, which was the sixth most important variable.
Person-Dependent and Person-Independent Algorithm Training
We trained three model types applicable to the estimation of workload of each individual operator. One “generic” model based on training data from all other operators expect the one in question, a “team” model for each of the six teams based on the training data from all three operators within a team, and a “person” model for each of the 18 operators based on training data from the given operator only. From “generic,” “team,” to “person” models, the relevance of the training data increases while the number of training cases reduces—gradually more relevant but less data available. All three model types were tested on each operator’s test data—the 30% withheld data for each operator. Thus, within each model type, the analysis resulted in 18 model accuracies.
Figure 8 shows that the person and team models performed better overall than the generic person-independent models despite having less data available for model training. The mean accuracy of the generic model was .58; the team model, .71; and the person model, .72. A nonparametric Friedman test showed that accuracy was significantly affected by the training approach, χ2(2) = 9.24, p < .01. Post-hoc comparisons using a Wilcoxon test with a Bonferroni-adjusted α level (.017) showed that both the person model and the team model were significantly different from the generic model.

Accuracy of cognitive workload estimation from generic, team, and person models.
Discussion
An operator’s cognitive workload during full-scope simulation scenarios could be classified from combining speech features and system features with an accuracy of 67%. Speech features and system features achieved similar accuracy, 63% versus 61%. The most important speech variables were pitch (f0), articulation rate, and amplitude, and the most important system variables were alarms and operator–system activities. Training algorithms specifically for a given person or for a given team of three operators improved mean accuracy to 72% and 71%, respectively.
Machine Learning Algorithm Performance
The differences in accuracy between algorithms were moderate and consistent across all three data types. RF and XGB achieved the highest classification accuracy. However, KNN and NN achieved similar levels of accuracy for the system data. The relatively good performance of RF corresponds with results reported by McDonald et al. (2020), while the relatively good performance of XGB corresponds with results reported by Momeni et al. (2019).
The 67% accuracy of RF trained on the combined data, and the 71% and 72% accuracy of the team and person models, respectively, are a promising result for the classification of cognitive workload in realistic accident scenarios. The accuracy achieved is below the accuracy of workload classification reported for simple individual tasks but similar to several results reported from realistic settings. For example, Yin et al. (2007) reported an accuracy of 78.9% for the relatively simple Stroop task, and Yap et al. (2010) reported an accuracy 84.4%. Charfuelan and Kruijff (2013) achieved a classification accuracy of 74% for realistic urban search and rescue training. McDonald et al. (2020) reported accuracy ranging from 24% to 73% for various machine learning algorithms used to classify driver distraction. However, Chen et al. (2012) demonstrated that relatively high accuracy can be achieved in realistic settings. They reported an accuracy of 82.2% over three levels of workload for operators of an emergency communication center. Compared to our study, Chen et al. (2012) identified and annotated shorter workload segments—10-s segments compared to our 2-min segments. This might help to determine a segment’s cognitive workload level. Longer segments pose a higher chance of including varying workload levels, so they are potentially more challenging to label. Further, Chen et al. (2012) utilized postsession segment coding by subject-matter experts. This might result in more consistent labeling compared to participant self-labeling.
The challenge of correctly classifying high cognitive workload was seen for all algorithms. We interpret this challenge as representing questions regarding measurement rather than about the type of algorithm or type of data. The proportion of high workload observed was low and concentrated on a few individual operators. Professional operators are trained to cope with high workload, so more extreme cases of cognitive workload might be necessary to classify workload as high. Related to this issue is the procedure for subjective labeling of workload. Even though the procedure was extensively explained to the operators before the data collection, future studies could verify the comprehension of the procedure during the data collection. The latter conflicts with the conventional experimental procedure of avoiding influence on participants’ responses during data collection but might be a measure needed to increase the quality of training data for machine learning algorithms.
Data Types, Important Speech, and System Variables in Cognitive Workload Estimation
The accuracies achieved with speech data and system data individually provide convergent evidence for speech features and system features being relevant for cognitive workload estimation. However, the lack of improved accuracy when combining these two relatively different sources of data is a surprising result. From a theoretical viewpoint, we expected speech features to indicate the effort of the operator in managing workload and to be highly influenced by individual differences between operators. On the other hand, alarms represent the task load from the plant process perspective and most alarms were common across crews.
That surprising result aside, we propose that the study demonstrates a potential for both speech data and system data to be used for estimating cognitive workload during realistic human–machine work. The top three speech features indicate that raising pitch, speaking faster, and speaking louder are relevant indicators of cognitive workload. This finding aligns with prior research (Brenner et al., 1994; Luig & Sontacchi, 2010). Alarms and operator–system interaction were found to be important system features, indicating a potential for using this type of feature to classify and measure cognitive workload.
Training Approaches and Accuracy Achieved
The study found that person-dependent and team-dependent models resulted in higher classification accuracy, with 72% and 71%, respectively, compared to a generic person-independent model (58% accuracy). The improved accuracy of person models corresponds with earlier findings (van Segbroeck et al., 2014; Vukovic et al., 2019) and may be related to both operator differences regarding features extracted and the operators’ self-rating of cognitive workload. Interestingly, the team models resulting from training the algorithm on a given team’s three operators performed on average as well as the person models. The results suggest that developing person models, whenever feasible, is better than using generic models for use in dynamic, realistic, human–machine research. The accuracy of the team models is interesting since this option provides training data from all team members as opposed to the person model involving a single operator.
Study Limitations and Future Research
An inherent challenge with simulator studies of realistically demanding human–machine work is limited controllable scenario progression. Increased testing and fine-tuning of scenarios with the participation of actual operators could provide scenario segments that professional operators would classify more consistently as high workload. In addition, scales for cognitive workload classification could be calibrated with each participant, both before and during data collection.
The relevance of system variables for estimating cognitive workload might depend on the context. For example, the utilization of alarms assumes that operators attend to activated alarms and utilize alarm information in their control of the system. This would be expected for most nuclear control room work where alarms are frequently relevant to the crew’s main tasks. However, one can imagine a context where operators perceive alarms as less relevant, either due to the alarm content or due to poorly designed alarm systems. Reviewing work procedures and/or discussing with users and subject-matter experts would be important to identify to what extent system variables potentially can serve as workload indicators.
The choice of ground truth measure guiding the classification is an important assumption for modeling. We chose a subjective scale because it is currently the most widely used (Matthews et al., 2015) and was feasible to obtain during the simulation. However, subjective scales have been criticized for lack of correspondence with performance-based and physiological workload measures (Matthews et al., 2020). Primary task measures (Hancock & Matthews, 2019; Young et al., 2015) remain challenging to define and implement for complex systems due to limited sensitivity to levels of workload (Hancock & Matthews, 2019; Wu & Li, 2013), while secondary task measures have been successfully applied to realistic process control (Wu & Li, 2013). Secondary task measures, such as the Detection-Response Task (International Organization for Standardization [ISO], 2016), have been extensively applied for cognitive workload measurement in driving environments (Strayer et al., 2019), military (Hollands et al., 2019), and laboratory settings (Castro et al., 2019). Physiological measures, for example, heart rate and heart rate variability (Mulder, 1992), have been found to be sensitive to cognitive workload in human–machine control room settings (Gao et al., 2013; Wulvik et al., 2020). In addition to ground truth, several of the above measures represent nonintrusive approaches for cognitive workload estimation. Future research could compare the approach presented here with these approaches to investigate the merits of this study’s approach relative to other approaches.
The study included male operators only, which might question the generalizability across gender. This limitation could be more relevant for the estimation of workload from speech features than for the estimation from system variables.
This study demonstrates the potential of utilizing speech features, user interaction, and alarms to estimate cognitive workload in realistic human–machine settings. However, if considering a fully developed approach, the accuracy achieved can be classified as moderate only. For practical reasons, this study covered a subset of the speech and user interaction features presented in the literature. Examples of relatively comprehensive feature extraction include Charfuelan and Kruijff’s (2013) extraction of speech features and Arshad et al. (2013) extraction of user interaction features. We believe that further development of the specification and extraction of features specific to complex human–machine interaction could improve the features’ sensitivity to cognitive workload, and thereby increase estimation accuracy.
Application Potential
The approach presented is applicable to many domains. System states, user interaction, and data from human sensors are becoming easily available, and high-quality verbal communication can be recorded even in team settings. Applications include cognitive load indication in complex training situations (van Merriënboer & Sweller, 2005) and workload assessment in human–machine design and evaluation (Reinerman-Jones et al., 2019). The approach investigated is nonintrusive, so it does not affect team tasks and can be applied to complex research and applied settings without requiring cumbersome, dedicated, data registration equipment (Chen et al., 2012). The approach can potentially be applied for real-time estimation of workload and can be suitable for user-adapted design.
Conclusion
This study investigated the classification of operators’ cognitive workload through machine learning by combining speech and system features. The achieved classification accuracies demonstrate the approach’s potential for dynamic nonintrusive cognitive workload estimation in realistic settings. However, we believe future studies should advance the specification and extraction of speech and system features to potentially increase classification accuracy. Future studies should investigate improved ground truth measures for classification, for example, by utilizing secondary tasks or physiological measures.
The findings suggest that the most promising cognitive workload indicators are the speech features, pitch (fundamental frequency), articulation rate, and amplitude, and the system feature of number of activated alarms. RF and XGB outperformed other algorithms in terms of classification accuracy, and training using person-dependent classification models improved accuracy over generic models. The nonintrusive approach presented is widely applicable to human–machine research and applied human factors.
Key Points
There are few nonintrusive measures available for evaluating cognitive workload in realistic dynamic human–machine interaction.
Earlier research on simple tasks has established a relationship between speech features, user–system interaction, and system variables with cognitive workload.
This study found, through the use of machine learning algorithms, that speech features and system variables occurring as part of a control room team’s natural work in a realistic full-scope simulator could classify the cognitive workload experienced by operators.
The types of variables used to classify cognitive workload are available and easily accessible in many research and applied settings. The approach is, therefore, feasible for many domains for the measurement and evaluation of cognitive workload.
Supplemental Material
Supplementary Material 1 - Supplemental material for An Investigation of Speech Features, Plant System Alarms, and Operator–System Interaction for the Classification of Operator Cognitive Workload During Dynamic Work
Supplemental material, Supplementary Material 1, for An Investigation of Speech Features, Plant System Alarms, and Operator–System Interaction for the Classification of Operator Cognitive Workload During Dynamic Work by Per Ø. Braarud, Terje Bodal, John E. Hulsund, Michael N. Louka, Christer Nihlwing, Espen Nystad, Håkan Svengren and Emil Wingstedt in Human Factors: The Journal of Human Factors and Ergonomics Society
Footnotes
Acknowledgment
The simulator data utilized in this study were originally collected for work funded by the OECD Halden Reactor Project, Norway.
Supplemental Material
The online supplemental material is available with the manuscript on the HF website.
Author Biographies
Per Ø. Braarud is a senior scientist in the Human-Centered Digitalization department at Institute for Energy Technology (IFE) in Halden, Norway. He received his MSc in cognitive psychology from the Norwegian University of Science and Technology in 1994.
Terje Bodal is a senior engineer in the Applied Data Science department at Institute for Energy Technology (IFE) in Halden, Norway. He received his BSc in informatics from Østfold University College, Norway, in 1994.
John E. Hulsund is a senior scientist in the Applied Data Science department at Institute for Energy Technology (IFE) in Halden, Norway. He received his MSc in experimental particle physics from the University of Bergen, Norway, in 1993.
Michael N. Louka is the department head of the Virtual and Augmented Reality department at Institute for Energy Technology (IFE) in Halden, Norway. He received his BSc(hons) and Associateship in computing and IT from University of Surrey, UK, in 1994.
Christer Nihlwing is a principal engineer in the Control Room and Interaction Design department at Institute for Energy Technology (IFE) in Halden, Norway.
Espen Nystad is a senior scientist in the Human-Centered Digitalization department at Institute for Energy Technology (IFE) in Halden, Norway. He received his MSc in psychology from the Norwegian University of Science and Technology in 2002.
Håkan Svengren is a principal engineer in the Control Room and Interaction Design department at Institute for Energy Technology (IFE) in Halden, Norway.
Emil Wingstedt is a research scientist in the Applied Data Science department at Institute for Energy Technology (IFE) in Halden, Norway. He received his MSc in theoretical physics from the Norwegian University of Science and Technology in 2002.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
