Learners' use of video controls in educational videos provides implicit signals of cognitive processing and instructional design quality, yet the lack of scalable and explainable predictive models limits instructors' ability to anticipate such behavior before deployment. We propose a scalable, interpretable pipeline for predicting population-level watching, pausing, skipping, and rewinding behavior as proxies for cognitive load from video content alone. Our approach leverages multimodal large language models (MLLMs) to compute embeddings of short video segments and trains a neural classifier to identify temporally fine-grained interaction peaks. Drawing from multimedia learning theory on instructional design for optimal cognitive load, we code features of the video segments using GPT-5 and employ them as a basis for interpreting model predictions via concept activation vectors. We evaluate our pipeline on 77 million video control events from 66 online courses. Our findings demonstrate that classifiers based on MLLM embeddings reliably predict interaction peaks, generalize to unseen academic fields, and encode interpretable, theory-relevant instructional concepts. Overall, our results show the feasibility of cost-efficient, interpretable pre-screening of educational video design and open new opportunities to empirically examine multimedia learning theory at scale.
Figures & tables
Figure 1: Pipeline for predicting learners’ interactions with online learning videos and explaining predictions from multimedia learning theory.
Field
#Courses
#Course runs
#Learners
Avg. #Videos
Avg. Length [min]
Mathematics
18
42
49,311
62.7
12.1
Neuroscience
15
28
78,067
43.7
13.5
Electrical Engineering
10
27
43,841
41.9
17.9
Physics
10
18
60,924
37.8
17.0
Management
9
15
19,451
44.3
10.9
Others
23
44
105,671
41.6
12.0
Table 1: Statistics of online learning data in our study. "Others" include Architecture, Biology, Computer Science, Education, Environmental Science, and Microengineering.
Figure 2: Three modalities of video segments around t are encoded by pre-trained transformers. A neural classifier predicts if Signalv(t) is among the top K% at timepoint t .
Modality
Feature
Description
Scale
Lit.
Ft
Formula
Math notation beyond single symbols is visible.
{0,1}
[ 11 ]
Ft
Instructor
Instructor’s head is visible in the video.
{0,1}
[ 9 ]
Ft
Screen
The user interface of a computer screen is embedded inside the image (e.g., a code editor).
{0,1}
[ 11 ]
Ft
Structured Information Visualization
A diagram, graph, schematic drawing, or table is on the slide (no formulas, no GUI).
{0,1}
[ 11 ]
Ft
Text Object
A printed or handwritten sentence, derivation, or bullet point is visible (no labels, titles, footers, text in screenshots or code editors, references).
{0,1}
[ 11 ]
Ft
Visual Complexity
The amount, complexity, and diversity of textual and graphical content objects visible. Ignore instructors and typical slide elements like logos or titles as elements.
1-5
[ 11 ]
Table 2: Rubric for coding CTML features of the video lecture content.
Signal →
Watchedv(t)
PausedAtv(t)
RewoundTov(t)
SkippedFromv(t)
Embeddings ↓
Top 10%
Top 5%
Top 10%
Top 5%
Top 10%
Top 5%
Top 10%
Top 5%
eT
0.58±0.02
0.58±0.03
0.62±0.02
0.62±0.02
0.59±0.02
0.61±0.02
0.58±0.02
0.57±0.02
eS
0.50±0.01
0.52±0.02
0.53±0.01
0.54±0.01
0.56±0.03
0.56±0.02
0.57±0.03
0.57±0.03
eF,t
0.59±0.02
0.59±0.03
0.64±0.03
0.68±0.03
0.60±0.03
0.62±0.04
0.57±0.03
0.57±0.02
eF,(t−10,t,t+10)
0.63±0.02
0.63±0.04
0.67±0.04
0.70±0.04
0.67±0.03
0.69±0.06
0.62±0.03
0.62±0.05
eV,{1,32,64}
0.65±0.03
0.67±0.04
0.71±0.03
0.71±0.03
0.71±0.02
0.76±0.03
0.62±0.04
0.62±0.03
Table 3: Predictive performance in terms of AUC ( ± std across 10 seeds) of moment-to-moment video interaction prediction from multimodal video content embeddings; Lift@K% ( ± std) is additionally reported for the final multimodal model ex .
Figure 3: AUC ( ± std, 5 seeds) for prediction of top 5% learner-video interaction moments on fields unseen during training as a measure of generalization of our classifier.
Signal →
Watchedv(t)
PausedAtv(t)
RewoundTov(t)
SkippedFromv(t)
Setting
Top 10%
Top 5%
Top 10%
Top 5%
Top 10%
Top 5%
Top 10%
Top 5%
Reference eˉx + Shared Weights in h1
0.67±0.02
0.69±0.03
0.74±0.02
0.74±0.02
0.72±0.02
0.76±0.03
0.64±0.02
0.66±0.03
Reference eˉx + No sharing in h1
0.66±0.02
0.68±0.03
0.72±0.02
0.74±0.02
0.71±0.02
0.77±0.03
0.61±0.03
0.65±0.03
No Reference eˉx
0.66±0.02
0.68±0.03
0.70±0.03
0.71±0.03
0.71±0.02
0.73±0.03
0.64±0.03
0.63±0.04
1. Qwen-VL Layer eTV,1
0.63±0.04
0.61±0.03
0.64±0.02
0.65±0.02
0.65±0.03
0.65±0.03
0.58±0.02
0.60±0.03
32. Qwen-VL Layer eTV,32
0.67±0.02
0.69±0.02
0.72±0.02
0.73±0.02
0.72±0.03
0.76±0.03
0.63±0.03
0.64±0.03
Table 4: AUC ( ± std, 10 seeds) for different neural architectures (top), three different layers of Qwen-VL (middle), and a smaller vision language model (bottom).
Figure 4: Learner-video interactions at 6,000 video moments (50% from top 5% ranks) by CTML features coded by GPT-5 (inter-rater agreement in the right panel). For example, moments without a formula have an average PausedAtv(t) rank of 67% within their video, whereas moments with a formula have a higher average rank of 71%.
CTML Features
Multimodal Embeddings ex
Sample Size →
500
1,000
2,000
5,000
500
1,000
2,000
5,000
19,000
PausedAtv(t)
0.66 ± 0.02
0.69 ± 0.01
0.68 ± 0.01
0.68 ± 0.0
0.62 ± 0.02
0.67 ± 0.01
0.70 ± 0.0
0.72 ± 0.0
0.77 ± 0.01
RewoundTov(t)
0.73 ± 0.01
0.72 ± 0.0
0.73 ± 0.0
0.73 ± 0.0
0.68 ± 0.02
0.71 ± 0.01
0.73 ± 0.0
0.74 ± 0.0
0.78 ± 0.01
Table 5: Classifier performance (AUC ± 5-fold CV std) when using CTML features versus embeddings to predict learners’ interactions across training dataset sizes.
Figure 5: TCAV values of CTML concepts and activations in our classifier H . Significant (*) values above 0.5 mean that the classifier was positively sensitive to the concept present in the activation space.