Applying functional near-infrared spectroscopy and eye-tracking in a naturalistic educational environment to investigate physiological aspects that underlie the cognitive effort of children during mental rotation tests

Spatial cognition is related to academic achievement in science, technology, engineering, and mathematics (STEM) domains. Neuroimaging studies suggest that brain regions’ activation might be related to the general cognitive effort while solving mental rotation tasks (MRT). In this study, we evaluate the mental effort of children performing MRT tasks by measuring brain activation and pupil dilation. We use functional near-infrared spectroscopy (fNIRS) concurrently to collect brain hemodynamic responses from children’s prefrontal cortex (PFC) and an Eye-tracking system to measure pupil dilation during MRT. Thirty-two healthy students aged 9–11 participated in this experiment. Behavioral measurements such as task performance on geometry problem-solving tests and MRT scores were also collected. The results were significant positive correlations between the children’s MRT and geometry problem-solving test scores. There are also significant positive correlations between dorsolateral PFC (dlPFC) hemodynamic signals and visuospatial task performances (MRT and geometry problem-solving scores). Moreover, we found significant activation in the amplitude of deoxy-Hb variation on the dlPFC and that pupil diameter increased during the MRT, suggesting that both physiological responses are related to mental effort processes during the visuospatial task. Our findings indicate that children with more mental effort under the task performed better. The multimodal approach to monitoring students’ mental effort can be of great interest in providing objective feedback on cognitive resource conditions and advancing our comprehension of the neural mechanisms that underlie cognitive effort. Hence, the ability to detect two distinct mental states of rest or activation of children during the MRT could eventually lead to an application for investigating the visuospatial skills of young students using naturalistic educational paradigms.

Spatial cognition is related to academic achievement in science, technology, engineering, and mathematics (STEM) domains. Neuroimaging studies suggest that brain regions' activation might be related to the general cognitive effort while solving mental rotation tasks (MRT). In this study, we evaluate the mental effort of children performing MRT tasks by measuring brain activation and pupil dilation. We use functional near-infrared spectroscopy (fNIRS) concurrently to collect brain hemodynamic responses from children's prefrontal cortex (PFC) and an Eye-tracking system to measure pupil dilation during MRT. Thirty-two healthy students aged 9-11 participated in this experiment. Behavioral measurements such as task performance on geometry problem-solving tests and MRT scores were also collected. The results were significant positive correlations between the children's MRT and geometry problem-solving test scores. There are also significant positive correlations between dorsolateral PFC (dlPFC) hemodynamic signals and visuospatial task performances (MRT and geometry problem-solving scores). Moreover, we found significant activation in the amplitude of deoxy-Hb variation on the dlPFC and that pupil diameter increased during the MRT, suggesting that both physiological responses are related to mental effort processes during the visuospatial task. Our findings indicate that children with more mental effort under the task performed better. The multimodal approach to monitoring students' mental effort can be of great interest in providing objective feedback on cognitive resource conditions and advancing our comprehension of the neural mechanisms that underlie cognitive effort. Hence, the ability to detect Introduction Spatial cognition is usually described as "the capacity to generate, retain, retrieve, and transform well-structured visual images in mind" (Lohman, 1994). This ability to understand and manipulate the objects' spatial characteristics has been associated with academic achievement in science, technology, engineering, and mathematics (STEM) domains (Wai et al., 2009;Khine, 2016). Psychology research has often measured spatial cognition via visuospatial tasks, including the mental rotation task (MRT) paradigm Metzler, 1971, 1988;Bruce and Hawes, 2015;Sladky et al., 2016). MRT requires the subject to determine whether pairs of figures exhibited are the same or mirror images. To accomplish this task, the subject creates a mental image and rotates it into a different orientation, which is related to spatial reasoning, visuospatial ability, and resembles everyday situations that require spatial orientation and the capacity to manipulate objects mentally (Guillot et al., 2012).
In cognitive neuroscience, the MRT has been applied to different ages and neural development stages. Studies with adolescents and adults indicate a relationship between MRT and mathematical achievement (Kyttälä and Lehto, 2008;Tolar et al., 2009;Wei et al., 2012). Visuospatial cognition studies on children highlight the relevance of spatial ability in many scientific domains and support intervention to develop spatial skills (Webb et al., 2007;Tzuriel and Egozi, 2010). However, there are very few studies that investigate the brain activity of children during MRT (Roberts and Bell, 2000;Wu et al., 2020). Indeed, many studies examine the neural substrate associated with MRT and its relationship to adults' learning process (Jordan et al., 2001;Harris and Miniussi, 2003;Zacks et al., 2003;Zacks, 2008). Neuroimaging studies have shown that some brain regions, including the superior parietal [Brodmann's area (BA) 40], premotor cortex (BA 6), and the frontal cortices [dorsolateral prefrontal cortex (dlPFC) BAs 9 and 46; ventrolateral prefrontal cortex (vlPFC) BA 44], are involved with bilateral activation regions during the MRT (Cohen et al., 1996;Hugdahl et al., 2006;Zacks, 2008;Wu et al., 2020). Other studies found that the modulation of activity in the frontal cortex reflects the participant's cognitive effort, i.e., mental resources required under cognitive demands (Prat et al., 2007;Ayaz et al., 2010Ayaz et al., , 2012aFishburn et al., 2014;Causse et al., 2017). For example, Causse et al. (2017) demonstrated that the functional near-infrared spectroscopy (fNIRS) was sensitive to task difficulty variations. Also, the prefrontal activation intensity sheds insight on the level of mental effort, indicating the amount of cognitive workload required by demanding tasks. Importantly, there are age-related differences in cognitively demands of mental processing and performance (McCabe et al., 2010). Older adults tend to express a higher effortful reaction and lower performance on cognitively demanding activity (Hess and Ennis, 2012;Westbrook et al., 2013), suggesting that poor performance could be a sign of depletion or reduced cognitive resources. It seems that greater levels of mental effort are related to lower performance (Hess and Ennis, 2012). Yet, little is known about how schoolchildren's cognitive effort relates to students' performance (Chevalier, 2018), and physiological measurements could be helpful to address this matter.
Eye-tracking (ET) is another reliable technological tool for research in the educational context. It provides a non-invasive and real-time measurement of pupil sizes from subjects during cognitive tasks (Lum et al., 2017;McGarrigle et al., 2017). Pupillometry studies indicate that pupil diameter changes in response to mental activity and increases with task difficulty in children (Boersma et al., 1970;Karatekin, 2004;Karatekin et al., 2007;Chatham et al., 2009;Johnson et al., 2014;Rojas-Líbano et al., 2019) and adults (Hess and Polt, 1964;Beatty and Kahneman, 1966;Granholm et al., 1996). Pupillometry measures during cognitive demanding mental arithmetic tasks suggest that pupil dilation may indicate mental activity, which includes mental arithmetic tasks (Hess and Polt, 1964;Ahern and Beatty, 1979;Jainta and Baccino, 2010;Klingner et al., 2011) and spatial ability tasks (Buckley et al., 2018;Campbell et al., 2018;Bauer et al., 2021a,b). For example, Toth and Campbell (2019) showed that the cognitive demand of the MRT elicited increases in participants' pupil diameters (Toth and Campbell, 2019). It is well known that the noradrenergic system has been implicated in mental effort and pupil dilation. However, the precise mechanism related to PFC activity and fluctuations in pupil size is still not fully understood (Bradley et al., 2008;Mathot, 2018).
There is growing recognition that fNIRS is moderately tolerant to motion artifacts and is a portable, cost-effective device (Ferrari and Quaresima, 2012;Curtin and Ayaz, 2018). All these together make fNIRS a suitable tool to use in conjunction with an eye-tracker device for monitoring the brain activity of children in a naturalistic setting such as realistic educational situations. Together, pupillometry and neuroimaging techniques are valuable for integrating mental effort signs and shedding light on the neural mechanisms that underlie cognitive load (Hosseini et al., 2017;Ísbilir et al., 2019). Although previous studies have investigated the cognitive effort's neural mechanism in adults, little is known about children's brain activity under visuospatial cognitive demands. Considering the spatial cognition relevance to school achievement (Giofrè et al., 2013), physiological measurements of students' mental effort during visuospatial tasks will make it possible to advance the understanding of cognitive aspects related to academic performance, specially geometry problemsolving development. We believe that traditional cognitive load measures could benefit from ET-fNIRS-based mental workload detection by combining additional information sources related to the mental effort of children during challenging tasks. Such a multimodal approach has the potential to contribute to a better understanding of young students' cognitive resources that underlie the learning process.
Multimodal applications can enhance physiological response detection from different modalities, which enable more accurate and robust measurements of cognitive states (Siegle et al., 2003;Shin et al., 2018), including mental effort (Laeng et al., 2012;Rozado and Andreas, 2015). Therefore, we propose the ET-fNIRS fusion as a promising framework for future investigations in educational research. Here, we sought to shed light on physiological responses related to children's cognitive resources during visuospatial tasks. More specifically, we wanted to improve our understanding of the relationship between pupil dilatation and hemodynamic PFC activity during a visuospatial task in children. Therefore, we use a multimodal approach by applying fNIRS and Eye-tracking simultaneously to investigate primary school children's mental effort expressed by brain hemodynamics changes and pupil size during the MRT. Based on the research cited above, we hypothesized that the higher bilateral PFC activity and pupil dilatation might be associated with higher mental effort during the spatial task and, consequently, lower children's performance (i.e., visuospatial scores).

Participants
The local Ethics Committee of UFABC approved all aspects of our experiment. The experiment was performed following all relevant guidelines and federal regulations. All subjects participated voluntarily and without any financial compensation. Thirty-two healthy students (fifteen boys and sixteen girls, 9-11 years old) from the 5th grade of a Brazilian elementary public school participated in this study. Parents/legal guardians of the young children provided written informed consent and were informed verbally of the purpose of the study and the safety of the experiments. All participants had normal or corrected to normal vision, and did not report present or previous neurological or psychiatric conditions.

Eye-tracking data acquisition
The pupil diameter was recorded using an ASL Mobile Eye-5 (Applied Sciences Laboratory, Bedford, MA, United States) with a capture rate of 60 Hz. The glasses were mounted with one small camera over the right eye. The ASL Mobil Eye-5 auto-calibrates to lighting conditions in the environment and compensates for head movements.

Functional near-infrared spectroscopy data acquisition
The hemodynamic signals were recorded using the continuous-wave NIRSport System (NIRx Medical Technologies, Glen Head, NY, United States). The cap montage was an array of eight pairs of optodes (sources and detectors), resulting in 28 channels (twenty long-range channels and eight short-distance channels) placed bilaterally over the prefrontal cortex (Figure 1), including the brain regions of the premotor cortex (BA 6), and dlPFC (BAs 9, 10, and 46; Koessler et al., 2009) (see Supplementary material 1 and Supplementary Figure 1 for more details). Source and detector distance were 30 mm with optodes positioned on the measuring cap regarding the 10-10 international system (Jurcak et al., 2007). The illumination sources emitted two wavelengths of near-infrared light (760 and 850 nm), and the signals were recorded at a sampling rate of 7.91 Hz. A dark cloth was placed over the measuring cap to reduce noise from the external light. We used the software NIRStar 15.2 (NIRx Medical Technologies, Glen Head, NY, United States) to record the brain activity data and to evaluate the signal quality of the channels.

Geometry problem-solving test
To assess the geometry problem-solving performance, we selected the questions from the Brazilian Basic Education Assessment System (SAEB). It is a national test that assesses the mathematical and language skill levels of Brazilian students from the first and last year of elementary school, Schematic representation of the fNIRS optodes positioned over the scalp. The position of the optodes follows the universal configuration of the 10-10 EEG international system. Channels are composed of neighboring optodes represented by the black edges, sources in red (S1-S8) corresponding to F3, AF7, AF3, Fz, Fpz, AF4, AF8, F4, and detectors in blue (D1-D7) matching F5, F1, Fp1, AFz, F2, Fp2, F6. # channel.
5th and 9th grades, respectively (Carnoy et al., 2015). We selected and printed eleven geometry multiple-choice problems from the SAEB test public database (SAEB, 2011). The test was designed to measure the general geometry problem-solving skills of 5th graders, and we adapted the questions in such a way that the students did not know the results beforehand.

Mental rotations test
The MRT was adapted for the 9-11 years old children. The task consisted of 70 slides with images in two dimensions following the geometrical patterns (abstract figures). There were two kinds of abstract figures (Figure 2) in the stimuli session: (1) one gray abstract figure placed on the center of the slide with a gray arrow indicating the left or right direction of rotation to the next slide, and (2) a pair of gray abstract figures identified as "A" and "B" with a similar form, but different orientations. Only one of the options was a 90 • rotation of the figure presented before. The task requires children to imagine the movement of rotation in the first slide to decide, in the next slide, which of the choices, "A" or "B, " match the stimuli presented before. There were ten trials with three slides of stimuli, three slides of alternatives, and one slide of the control condition in which subjects were instructed to fixate on the cross until the onset of the stimuli. The child had 30 s to rest, 7 s to imagine the rotated geometric figure, and 3 s to choose the alternative, three times per block, yielding 30 s of stimuli. All the sessions had the same level of difficulty to counterbalance the different phases of the task. We used the software NIRStim (NIRx Medical Technologies, Glen Head, NY, United States) to present the visual stimuli and collect the fNIRS data.

Experimental design
The MRT was performed in the students' usual place at school. Each child participated in the task individually in a session that lasted approximately 30 min. First, students were randomly called and asked to answer eleven geometry questions from the Brazilian basic education assessment system (SAEB). After solving all the eleven mathematical problems with multiple choices printed on paper, participants were assigned to another room fitted with the fNIRS and Eye-tracker. Two experimenters controlled the equipment and observed the child during the task. Participants were instructed about the experiment and asked to answer a practice trial of MRT to make sure that they understood the task. If the child made a mistake, the experimenter would repeat the trial. After the instruction trials, experimenters set up the fNIRS and eyetracker. Students were then told that they would next play on the computer screen the same kind of game as in the Schematic diagram of the experimental procedure. Fixation cross (the 30 s) was followed by the first cue (7 s) as one gray abstract figure placed on the center of the slide with a gray arrow indicating the left or right direction of rotation. A stimulus "response" was then presented (3 s) as a pair of gray abstract figures identified as "A" and "B" with a similar form, but different orientations at 90 • according to the cue (to the right or left visual field) so the participants could answer verbally. The child had 30 s to rest, 7 s to imagine the rotated geometric figure, and 3 s to choose the alternative, three times per block, yielding 30 s of stimuli.
instruction trials of MRT. The stimuli were presented, followed by the choice "A" and "B, " in which students answered verbally. Participants' answers were collected during the experiment for further analysis.

Eye-tracking data analysis
Five children were excluded from the sample due to excessive missing pupil data on the MRT (bad frames), yielding data from twenty-six children. The data were inspected for artifacts such as eye blinks or undetected signals during the task. The high-frequency noise also was also considered an eye tracker artifact. Therefore, we removed the noise from the data by smoothing it using a moving median (Gao et al., 2007). Then we implemented a computer algorithm to automate the artifact removal process for all subjects. We considered recordings insufficient when less than 40% of data remained after cleaning the data to remove artifacts. Children's median pupil dilation variation (%) was calculated over each second of the task compared to the previous control condition. Then, we determined the median across the experimental trials. The 10 s interval preceding each stimuli block was treated as the baseline. Finally, a t-test was applied to test significant differences in the pupil dilation mean between control and task conditions.

Functional near-infrared spectroscopy data analysis
The hemodynamic signals were processed and analyzed by the software nirsLAB v2019.04 (NIRx Medical Technologies, Glen Head, NY, United States). The signals were truncated before the first block and after the end of the last block. Then, the modified Beer-Lambert law was applied to convert the optical signals of each wavelength (760 and 850 nm) in concentration changes of oxy-Hb and deoxy-Hb. The general linear model (GLM) method (Huppert, 2016) was applied to analyze changes in both oxy-Hb and deoxy-Hb amplitudes during the MRT. In the individual-level analysis, we obtained the activation beta coefficients for each channel of the participants. Then, we used a group GLM to combine all participants' beta coefficients yielding the group statistical activation map (t-test). The rate of accepted false-positive results was set at 5%. Bonferroni correction was applied for multiple comparisons for 20 longdistance channels (resulting p-value threshold: <0.0025). The hemodynamic signals were extracted from the channels that presented significant activation in the group analysis. Finally, we calculated the signal block average among the participants to evaluate the channel's hemodynamic response. The spatial representation of the brain regions related to each channel was rendered from the BrainNet Viewer toolbox (Xia et al., 2013) 1 .

Behavior data analysis
We analyzed the children's performance on the geometry problem-solving test and the MRT by counting the correct answers on tasks. The scores on both visuospatial tasks were correlated by Spearman's correlation coefficient. We also performed Spearman's correlation analysis to quantify the potential association between behavioral data with pupillary responses and cortical activation. In all cases, p < 0.05 (uncorrected) was considered statistically significant. When the short-distance channels were statistically significant to behavioral data, we considered the long-distance correspondent channel insufficient quality data for correlation analysis.

Pupillary response
The timeline analysis of the pupil size data indicated three distinct periods of pupillary variation (%) from the beginning (0 s) until the end of the MRT (p = 0.0006), as shown in Figure 3. First, there was a period without any significant stimuli (10 s of resting state as the baseline), followed by an increase in the children's pupil size with the stimulus display. The results show a peak of pupil dilation around 10 s after the children imagined the figure. It suggests a pupil dilatation follows a visuospatial task. The responses to the mental rotation stimuli appeared to elicit more pupil dilation than the neutral stimuli. It indicates that it is possible to identify two distinct mental states of rest or activation during the MRT. There is no significant correlation between the mean amplitude of pupil dilation and the number of correct answers, neither in geometric problem-solving nor in the MRT (p > 0.05 in all cases).

Brain-behavior correlation
The Spearman's correlation analysis revealed correlation between MRT scores and oxy-Hb variation at channel 25 (r = 0.43; p = 0.01; Figure 5B) and channel 5 (r = 0.41; p = 0.02), but we considered the channel 5 as insufficient quality data because the corresponding short-distance was also significantly correlated to the task performance (Channel 7; r = 0.38; p = 0.03; see Supplementary material 3). There is also significant correlation between geometry performance and deoxy-Hb variation at channel 8 (r = 0.40; p = 0.02; Figure 5C), at short-distance channels 4 (r = 0.33; p = 0.007), 7 (r = 0.41; p = 0.02), 24 (r = 0.39; p = 0.02) and oxy-Hb variation at short-channel 7 (r = 0.39; p = 0.03; see Supplementary material 3). All other comparisons were not significant. The spatial representation of the brain regions related to each channel is shown in Figure 5A, for the mental rotation tasks scores and the geometry test scores (channels respectively, highlighted in purple and orange). Channel 25 is on the left prefrontal cortex (PFC), and channel 8 is on the left dorsolateral prefrontal cortex (dlPFC).

Cortical activation
The results revealed a significant deactivation (Deoxy-Hb) of the left dlPFC (Ch.6; p-value corrected = 0.002), and the right dlPFC (Ch.22; p-corrected = 0.001; and Ch.23; p-value corrected <0.001, see all results in Supplementary material 2). Brain regions represented in Figure 6A show that the channels associated with the correlation between Oxy-Hb variation and MRT are placed on dlPFC. Mean changes in oxy-Hb levels in the other prefrontal regions that we evaluated were not significant. Figure 6B shows the mean group oxy-Hb and deoxy-Hb signal change during the MRT for the channels 6,22,and 23 (Ch.6,Ch.22,and Ch.23,respectively) indicating that the oxy-Hb concentration increase while the changes in deoxy-Hb were in the opposite direction, as expected.

Discussion
This study aimed to investigate the 9-to 11-year-old schoolchildren's cognitive effort expressed by changes in pupil diameter and PFC neural activity under spatial cognitive demand. Our findings suggest that the visuospatial task was able to demand students' mental effort, which was detectable by the multimodal approach using pupillometry and optical neuroimage. To our knowledge, this is the first study to apply fNIRS and Eye-tracking to assess the cognitive resources required under mental demands through cortical hemodynamic activity and the pupil responses of children during the MRT. 10.3389/fnhum.2022.889806 FIGURE 3 Pupil diameter variation timeline during the MRT. The red line indicates the beginning of the task (stimuli), which had three different figures to rotate per block, yielding 30 s of stimuli (p = 0.0006) and 10 s of baseline. The solid line indicates the mean (µ) and the dashed line indicates ± 1 standard error (σ). Scatter plot between the number of correct responses on MRT and geometry problem-solving test (scores).
The results indicate a significant positive correlation between the children's MRT and geometry test scores. There are also significant positive correlations between dlPFC hemodynamic signals and visuospatial task performances (MRT and geometry problem-solving scores). Our results are consistent with dlPFC activity related to higher-level cognitive processing (Miller and Cohen, 2001;Vassena et al., 2014Vassena et al., , 2019. Moreover, we found significant activation in the amplitude of deoxy-Hb variation on the dlPFC and that pupil diameter increased during the MRT, suggesting Correlation graphs and fNIRS channels positions. (A) 3D-rendered brain template (ICBM152) with left (L) and right (R) hemispheres from BrainNet Viewer toolbox (Xia et al., 2013). The fNIRS channels are positioned across the cortex. The red points represent sources, the blue points are detectors, and the green lines are channels. The purple highlight indicates the channel associated with the correlation between oxy-Hb variation and mental rotation task (MRT). The orange highlight indicates the channel associated with the correlation between Deoxy-Hb variation and geometry problem-solving (GPS). (B) Scatter plots between the number of correct responses on MRT and oxy-Hb variation at channel 25 (r = 0.43; p = 0.01). (C) Scatter plots between the number of correct responses on the geometry test and deoxy-hemoglobin (deoxy-Hb) variation at channel 8 (r = 0.4; p = 0.02). Group level HbO cortical activation map. (A) 3D-rendered brain template (ICBM152) with left (L) and right (R) hemispheres from BrainNet Viewer toolbox (Xia et al., 2013). The fNIRS channels are positioned across the cortex. The red points represent sources, the blue-points are detectors, and the green lines are channels. The purple highlight indicates the channel associated with the correlation between deoxy-Hb variation and mental rotation task (MRT). Thresholded SPM t image, p-value = 0.0025 for deoxy-Hb. (B) Block-averaged hemodynamic response (mean ± SD) computed by averaging the oxy-Hb (red) and deoxy-Hb (blue) signals during the MRT for channels 6, 22, and 23. that both physiological responses are related to mental effort processes during the visuospatial task. Our findings indicate that children with more mental effort under the task performed better, which was unexpected. According to the general cognitive load theory (Sweller and Chandler, 1991;Sweller et al., 1998), better problem solvers tend to express less mental effort than novice learners' performances, as demonstrated in preview studies (Çakır et al., 2011;Shaw et al., 2018). A possible explanation is an age-related difference in our study. Unlike adults, children do not have well-developed visuospatial skills, so they need more activation to solve geometric problems and perform well. Despite the potential to investigate learning and development with eyetracking, few studies used task-evoked pupil changes analysis on children (Eckstein et al., 2017). So, it is possible to apply pupillometry to elucidate students' cognitive processes and development.
Spatial cognition improves with age in childhood (Orde, 1996) and declines with age in adulthood (Pak, 2001). Furthermore, studies investigating cognitive efficiency in young and older adults show age-related differences in mental effort and performance (Baltes et al., 1999;Westbrook et al., 2013;Westbrook and Braver, 2015). Older adults tend to express a higher effortful reaction and lower performance on work memory tasks (McCabe et al., 2010;Westbrook et al., 2013) and math problemsolving (Hess and Ennis, 2012). On the other hand, studies with children demonstrated that young children primarily engage the frontal cortex when solving numerical tasks, decreasing by age and differences in brain development (Houdé et al., 2010;Artemenko et al., 2018;Soltanlou et al., 2018a,b). In this sense, we believe that adults would have different findings from children because of neural development. Such information has important implications for understanding students' spatial cognition and developing an intuitive understanding of fundamental geometric concepts.
A significant positive correlation was observed between MRT and geometry problem-solving performance. It indicates that children with good performance in figure rotation also present better performances in our geometry test based on the Brazilian mathematics assessment. It corroborates that our visuospatial task requires mental rotation skills associated with geometry problem-solving abilities. Geometry is a fundamental approach to interpreting the physical environment and everyday situations (Hwang et al., 2009). It contains basic mathematical concepts and demands spatial cognition during the learning process of shapes, size, quantity, and length (Biber et al., 2013;Hertanti et al., 2019). Previous studies have shown that spatial cognition is a good predictor of mathematics achievement, specifically geometry (Delgado and Prieto, 2004). Similarly, our results are in line with the idea that spatial cognition is relevant to geometry performance. Our analysis correlation between MRT and geometry problemsolving underscores the critical role that spatial cognition plays in mathematical achievement and has implications for educational practices.
There seems to be a lack of evidence for neural correlates of children's cognitive effort and MRT, maybe for neuroimaging technique limitations on movement artifacts. However, the last few years have shown an increasing technology development that allows studies in schools with children (Soltanlou et al., 2018a), which includes naturalistic experimental settings with fNIRS (Ferrari and Quaresima, 2012;Balardin et al., 2017;Barreto et al., 2021) and Eye-Tracking (Epelboim and Suppes, 2001;Bolden et al., 2015) outside the laboratory environment. Our study is in line with the idea that fNIRS and eyetracking are suitable for measuring the cognitive processes of schoolchildren in naturalistic settings (Mücke et al., 2018), and eye-tracking is a valuable device to use to measure physiological activations during spatial cognition stimuli in a familiar environment.
In addition to measuring cortical oxygenation, pupil diameter is also an important physiological measure of cognitive effort. Previous studies have suggested that changes in pupil size could be an index of how intensely the processing system is operating (Just and Carpenter, 1993;Laeng et al., 2012). The locus coeruleus (LC) is a brainstem nucleus that represents the primary source of norepinephrine (NE) in the brain (Kardon, 2005). The LC-NE system modulates the autonomic nervous system that influences pupil diameter through parasympathetic constriction and sympathetic dilation pathway (McDougal and Gamlin, 2008). The seminal study of Aston-Jones and Cohen (2005) demonstrated that pupil size is straightly related to fluctuations of the LC-NE activity. The same noradrenergic system has been implicated in arousal (Waterhouse, 2003;Aston-Jones, 2005; and mental effort (Bradley et al., 2008;Mathot, 2018). In fact, the LC widely projects throughout the PFC regions (Aston-Jones and Cohen, 2005) and seems to module neuronal network activity in cognitively demanding tasks (Robbins and Arnsten, 2009;Suttkus et al., 2021).
Although the LC-NE system presumably promotes the relationship between pupil dilation and mental effort, the mechanism in which the PFC activity areas influence the changes in pupil size is still not elucidated. Uncertainty about the mechanisms through the relationship between pupil size and PFC activity during cognitively demanding tasks emphasizes the need for further investigation into this matter. Multimodal studies have the potential to investigate how different areas of the nervous system are associated, revealing new levels of organization and functioning. For example, a recent multimodal ET-fNIRS study demonstrated a functional relationship between the lateral PFC and the LC that may vary at different difficulty levels during the work memory task (Yeung et al., 2021). It is already known that pupil dilates with cognitive load in adults (Van Der Meer et al., 2010;Van Der Wel and van Steenbergen, 2018) and children (McGarrigle et al., 2017). Our data relates schoolchildren's pupil dilation during the MRT, suggesting a task-evoked increasing effort by spatial cognition demand. Further research should investigate the neural network specialization in visuospatial demand tasks across different age stages and neural development.
Interestingly, there are pieces of evidence that the brain activation of children can change after training. Soltanlou et al. (2018a) demonstrated that 2-week multiplication training with children reduced brain activation on the frontoparietal network. Furthermore, neuroimaging studies have suggested that children show decreased brain activation of the frontal cortex with increasing age (Houdé et al., 2010;Artemenko et al., 2018), which indicates reduced reliance on cognitive effort and attentional resources. Indeed, the cognitive effort is related to the degree of engagement with demanding tasks, and high engagement tends to enhance performance and attention (Kaplan and Berman, 2010). Considering the challenges to measuring the cognitive load (Paas et al., 2003), there is an increasing interest in the educational research field to apply a combination of measurements to obtain reliable and objective physiological data of cognitive effort in real-time (Kruger and Doherty, 2016).
The multimodal methodology has been used recently as a framework for the physiological measurement of cognitive load (Hosseini et al., 2017;Ísbilir et al., 2019;Larmuseau et al., 2020;Vanneste et al., 2021;Yeung et al., 2021). Studies applying pupillometry-fMRI concurrently measure have shown that pupil dilation was temporally correlated with activation of adults' brain regions implicated in workload and cognitive control during working memory and decision-making tasks (Siegle et al., 2003;Satterthwaite et al., 2007). Similarly, we measured cortical activity and pupillary response concurrently using fNIRS and eye-tracking glasses to examine the neural systems linked to pupil dilation under mental rotation tasks. We also explored which PFC regions activate more during the visuospatial task. We believe that pupil dilations could be an index of mental effort differences indicating the mental effort linked to the allocation of cognitive resources on cognitively demanding tasks. We hypothesize that the change in pupil diameter could suggest how hard the student is trying to solve the problem, indicating the performance slightly. However, we did not apply different level difficulties in the MRT trials, which is a relevant limitation of our study. We found some shortchannels signals significantly correlated to the task performance, which indicate extracerebral perfusion (e.g., forehead skin vessels) may possibly be an interference in our hemodynamic signals. Although such short channels do not correspond to long-distance channels that correlated significantly with the task performance, we considered this limitation worthy of mention.
One point of remark is that two fNIRS statistical analyses were conducted. The first focused on the significance of the contrast between task and rest. The second had as concern whether the amplitude of activation was correlated with behavioral data. We could use the results of the first analysis to select the channels for the second one. However, we opted to consider all channels in the latter and to be conservative by applying Bonferroni correction for multiple comparisons. We justify this choice since the two analyses are not only complementary, but they are testing different null hypotheses. The second point of remark is that although both oxy-and deoxy-Hb changes are consequences of the same neurophysiological process (the neurohemodynamic coupling process), the signal-to-noise ratio and susceptibility to systemic artifacts are not the same and may also depend on brain location. Thus, the statistical findings are not always reciprocal since the effect sizes are different. Moreover, for the same reasons, the lack of statistical power could explain why the correlation with task performance was only found at the right dlPFC.

Conclusion
Our findings suggest that children with more mental effort during MRT performed better in the visuospatial and geometry problem-solving tests. The multimodal approach to monitoring students' mental effort can be of great interest in providing objective feedback on cognitive resource conditions and advancing our comprehension of the neural mechanisms that underlie cognitive effort. Hence, the ability to detect two distinct mental states of rest or activation of children in real-time during the MRT could eventually lead to an application for investigating visuospatial skills of young students using naturalistic educational paradigms with a multimodal approach by applying the fNIRS and Eye-Tracker devices.

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement
The studies involving human participants were reviewed and approved by the UFABC. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author contributions
RS and JR designed the study, collected and analyzed the data, and wrote the manuscript. AA and CB wrote and reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.