Abstract
Trait anxiety is characterized by heightened apprehension and worry about future events, with individuals experiencing it at varying levels. Identifying physiological markers of these differences could enhance our understanding of anxiety-related processes. We tested the hypothesis that resting-state rhythmic electroencephalography (EEG) activity (oscillations during alternating 1-minute eyes-open/eyes-closed rest conditions) predicts trait anxiety, as measured by the Behavioral Inhibition Scale (BIS). Multiple regression revealed that during eyes-open rest, higher frontal midline theta and lower frontal midline alpha activity significantly predicted BIS scores. A more stringent analysis confirmed that frontal midline theta and alpha activity met stricter criteria for rhythmicity. Interestingly, only parietal midline alpha power during eyes-open rest, but not oscillations, was predictive of BIS, suggesting that rhythmicity may differentiate the functional roles of frontal and posterior alpha in relation to anxiety. We also found that alpha asymmetry (power or oscillations) did not predict BIS scores. Our results suggest that heightened frontal midline theta during eyes-open rest may reflect increased baseline vigilance in individuals with high trait anxiety. Although speculative, lower frontal midline alpha may indicate greater cortical activation in regions involved in inhibition and cognitive control. Moreover, by differentiating power and oscillation effects, our findings help resolve ambiguities in related research.
1 Introduction
Trait anxiety is associated with a tendency toward increased apprehension and worry about future events (Barlow, 2004), which, in turn, can influence emotional and cognitive processes. There is accumulating evidence that some of these effects are driven by rhythmic brain activity. Here we focus on brain activity during rest (no task), when we expect anxiety-related apprehensive and worrying thoughts will be prominent. As we summarize below, evidence linking resting rhythmic neural activity to trait anxiety remains mixed, likely due to methodological differences in signal analysis, such as measures of rhythmic (oscillatory) activity also being strongly influenced by non-rhythmic activity, reducing the specificity of the measure, as well as variations in how the rest state is induced. In this study, we focus more directly on rhythmic activity (while also comparing it with the more conventional, less-selective power measure) during a rest state that is designed to plausibly allow anxiety-related thoughts and emotions to emerge.
Theories of anxiety. According to reinforcement sensitivity theory, two primary systems regulate emotion and behavior: the Behavioral Inhibition System (BIS) and the Behavioral Activation System (BAS; Gray, 1982). BIS manages motivational conflict by interrupting behavior in uncertain or risky situations, driving anxiety and prompting risk assessment. It heightens vigilance in response to perceived threats, especially when conflicting motivations, such as approach and avoidance, are present. Unlike the Fight-Flight-Freeze System (FFFS), which handles immediate fear responses to clear danger, the BIS governs anxiety-driven caution and deliberation (Gray & McNaughton, 2000). In contrast, the BAS drives approach behaviors, motivating individuals to pursue goals or rewards.
Individuals with higher levels of trait anxiety often exhibit an overactive BIS, leading to excessive behavioral inhibition and an increased likelihood of avoiding goal-directed behaviors to prevent possible adverse outcomes, such as declining promising opportunities due to fear of failure (Gray, 1982; Gray & McNaughton, 2000). Moreover, high BIS levels have been associated with greater trait anxiety and neuroticism (Campbell-Sills et al., 2004; Carver & White, 1994; Smits & Boeck, 2006). Although anxiety serves a necessary evolutionary purpose and behavioral inhibition is essential for appropriate functioning, an overactive BIS and elevated trait anxiety can become maladaptive, significantly impacting daily life. Therefore, identifying the neural correlates of BIS is crucial for advancing our understanding of anxiety-related processes and dysfunctions.
Despite the theoretical strength of BIS models, the connection between animal-based frameworks and human trait anxiety remains unclear. This study addresses these gaps by focusing on rhythmic neural activity, isolating its contributions from non-rhythmic sources, and refining the neural signals associated with trait anxiety. Our findings demonstrate that neural rhythms provide valuable insights into stable personality traits, advancing understanding of how intrinsic brain activity shapes individual differences in emotion, cognition, and behavior. Additionally, this work enhances methodologies for studying resting-state neural rhythms, with broader applications for exploring other personality dimensions and emotional traits influenced by similar neural mechanisms.
The value of investigating resting rhythmic brain activity related to anxiety. Electroencephalography (EEG) is commonly used to measure such neural correlates of anxiety, including rhythmic activity. Much of the existing research has focused on the connection between trait anxiety and neural activity during specific tasks or behaviors. In contrast, resting-state EEG activity, also viewed as a neural trait, reflects relatively stable neural patterns that are independent of task engagement (Näpflin et al., 2007; Nash et al., 2014). As such, linking trait anxiety to resting-state EEG may provide deeper insights into intrinsic anxiety-related neural processes, unaffected by the external demands of tasks or stimuli but potentially more influenced by internalized processes, such as intrusive thoughts, rumination, and worry.
The nature of the resting state remains debated, with some suggesting it represents mere “idling,” while others argue it involves complex cognitive processes such as daydreaming and rumination (Gusnard et al., 2001; Menon, 2023; Pfurtscheller, 1992). Evidence also suggests that sleep may be common during resting states, especially if they extend beyond 3 minutes (Tagliazucchi & Laufs, 2014). This debate highlights that the resting state is not homogenous; for some individuals, it may predominantly reflect stable, low-arousal neural activity, whereas for others, especially those prone to anxiety, it may involve more active cognitive processes, such as rumination. This variability suggests that mental traits can uniquely influence the nature of rest and how it is reflected in EEG recordings. Considering the high test–retest reliability of the resting-state EEG (Metzen et al., 2022; Nash et al., 2014), when recorded under 3 minutes, it becomes a valuable tool for assessing individual differences in the neural correlates of personality traits.
Methods of quantifying rhythmic activity, limitations, and solutions. Conventional spectral analysis measures total power which is known to combine both rhythmic (oscillatory) and non-rhythmic sources, collapsing these different processes together (He, 2014; Klimesch, 1999; Whitten et al., 2011). Distinguishing between rhythmic and non-rhythmic power sources is crucial as they reflect different physiological processes: rhythms indicate synchronized neural firing and network coordination, while non-rhythmic power can include transient events, or other non-oscillatory neural activity. To address these issues, we use the Better OSCillation detection (BOSC) method, which calibrates measures based on the signal’s background spectrum (“colored noise”) and applies additional thresholds for rhythmicity (Caplan et al., 2001). Importantly, this distinction may help clarify discrepancies in the literature, as combining rhythmic and non-rhythmic signals without distinguishing their nature could obscure the specific contributions of each to anxiety-related neural activity. Therefore, separating these signal types allows us to more accurately identify neural oscillations linked to behavioral inhibition and anxiety regulation, providing insights that may be overlooked when only total power is considered.
1.1 Resting-state alpha
Alpha (8–12 Hz) is a dominant rhythm in adult EEG and is most often and easily detected during rest (Berger, 1929). Traditionally, alpha power was believed to be inversely related to cortical activation in adults (Davidson, 1992). This idea was initially supported by the observation that, at posterior electrodes, alpha power is higher when individuals close their eyes (reduced visual cortex activation), and lower when their eyes are open (increased visual cortex activation). Further research confirmed that alpha activity is also suppressed by other, non-visual mental processes (Allen et al., 2004). This led to the hypothesis that different cortical areas exhibit their own alpha rhythms, representing a resting or “idling” state of those brain regions (Kropotov, 2009; Palva & Palva, 2007).
Beyond single-electrode analysis, researchers have also investigated global alpha (across-electrode mean activity) measurements. Consistent with the “idling” hypothesis, transitioning from a less aroused Eyes Closed (EC) condition to a more aroused Eyes Open (EO) condition decreases global alpha power (Barry et al., 2007). Evidently, this change is accompanied by increased skin conductance levels, a well-established physiological marker of arousal and stress. The negative correlation between EC resting global alpha power and skin conductance further supports the view of alpha activity as an inverse marker of arousal (Barry et al., 2020). Collectively, these studies suggest that higher alpha power is associated with idling or inactive states, whereas lower alpha power corresponds to more active, aroused states, which may also be linked to higher anxiety.
Surprisingly, a contrasting body of research indicates that higher global alpha power, both during rest and task presentation, is positively associated with trait anxiety (measured by Taylor Manifest Anxiety and State and Trait Anxiety Inventory, STAI) and behavioral inhibition (measured by the Gray–Wilson Personality Questionnaire; Knyazev et al., 2002, 2003, 2004). Overall, Knyazev and colleagues found that higher frequencies, such as alpha, were predominant in individuals with higher BIS scores, and lower frequencies (theta and below) were more prevalent in individuals with higher BAS scores. Notably, their analysis sometimes combined global mean alpha power across both eyes-open and eyes-closed resting states (Knyazev et al., 2003), or included only eyes-closed rest (Knyazev et al., 2004). This contrasts with studies that focus on specific resting states (eyes open vs. eyes closed), which are considered functionally independent, potentially introducing ambiguity in interpreting the results.
Resting-state frontal alpha asymmetry (right hemisphere alpha activity compared with left), and its association with affect and emotion, is another contentious subject. Evidence from decades of research suggests that frontal cortical asymmetry explains a large amount of individual variability in emotional behavior and affective style, or personality traits (Davidson, 1992; Sutton & Davidson, 1997). In line with this, earlier studies demonstrated that the right prefrontal cortex is a biological substrate of withdrawal behavior and negative affect, and higher BIS is related to greater left frontal alpha power compared with right frontal alpha power (Sutton & Davidson, 1997). In contrast, Harmon-Jones and Allen (1997) did not find BIS to be related to the greater left frontal alpha, and Coan and Allen (2003) found only a weak statistical trend at neighboring electrodes. Moreover, a recent review by Firth et al. (2024) concluded that findings about BIS and frontal alpha asymmetry are largely inconsistent, and mainly absent during rest. (See also Kołodziej et al. (2021), who found a similar null relationship between alpha asymmetry and depression.)
This inconsistency can be partially attributed to the complexity of the BIS system and its relationship with withdrawal behavior (Coan & Allen, 2003). Notably, Gray’s original Behavioral Inhibition System (Gray, 1982) was later divided into two distinct systems, the anxiety-related BIS and the fear-related Fight-Flight-Freeze System (FFFS) (Gray & McNaughton, 2000). Given that this distinction was not reflected in behavioral inhibition measuring scales, if withdrawal tendencies activating right frontal brain are triggered exclusively by threat cues and not anxiety, only certain aspects of the BIS scales might be linked to increased right frontal neural activity (Coan & Allen, 2003).
Although considerable attention has been given to posterior and global alpha, as well as frontal alpha asymmetry, frontal midline alpha power remains largely unexplored in relation to rest and its connection to anxiety and the behavioral inhibition system. Given the importance of anterior brain regions in negative emotions and anxiety, it is possible that beyond asymmetry, significant relationships may also emerge in midline electrodes (Gray & McNaughton, 2000).
The focus on the difference-based frontal asymmetry over single-electrode measures may stem not only from theoretical assumptions but also from methodological challenges. Consider, for example, that brain activity is usually referenced to a control condition, often a pre-stimulus rest state. If rest is the condition of interest, there is no obvious control to subtract out. This poses a difficulty interpreting resting-state EEG signal without comparison with a reference state. It also means there is no calibration of the measurement scale across participants, which limits removing large sources of subject variability simply due to measurement scale (which could be due to numerous effects of no interest, such as skull thickness). As a difference measure, asymmetry is thus appealing because the opposite hemisphere functions mathematically similarly to a control condition, providing some calibration. We speculate that because interpreting absolute power measurements is challenging without baseline values, this may have led researchers to favor asymmetry measures that essentially self-calibrate through comparisons between left and right hemispheric activity. However, we can now address the calibration problem with a method such as BOSC, which uses the signal’s own statistical properties to tune oscillatory detection thresholds (Caplan et al., 2001). This approach makes the oscillation measures comparable across participants and interpretable without additional mental-state reference. Importantly, this approach reveals how midline alpha activity contributes to anxiety and the behavioral inhibition system, complementing the insights gained from global and asymmetry measures.
1.2 Resting-state theta
Generally, minimal to no theta rhythm is detected during short periods of rest in human EEG recordings (Caplan et al., 2015). This poses a significant challenge for studying resting-state theta oscillations, especially when compared with the well-documented hippocampal theta in animals or task-related cortical theta in humans. Consequently, most studies have focused on task-related theta activity in relation to BIS and anxiety, whereas few have explored the connection between resting-state theta and anxiety. Yet, despite its apparent absence in most individuals, if the theta rhythm is not universally present but does legitimately occur in a subset of individuals, we may find that these sporadic theta oscillations reflect important individual differences. Our study aims to address this gap by investigating whether variability in resting-state theta could serve as a meaningful marker for anxiety-related processes.
On the other hand, the link between task-related theta activity and the behavioral inhibition system has been extensively studied, and the theoretical view of task-related theta activity might also extend to rest. The task-related theta rhythm was first investigated to understand the neuropsychology of anxiety in animal models (Gray, 1982). It has been suggested that theta power in the rodent hippocampus (over a more generous range, 4–12 Hz) is a neural hallmark of BIS and anxiety processes (Gray, 1982; Gray & McNaughton, 2000). The corresponding human model suggests that internal conflict between equally compelling but incompatible goals activates the septohippocampal system, which then resolves the goal conflict through recursive communication with the neocortex, facilitated by rhythmic activity in the theta range (Corr, 2008; Gray & McNaughton, 2000). Indeed, some studies indicate that goal conflict increases theta activity, aligning with the septohippocampal theory of anxiety (Cavanagh & Shackman, 2015; Moore et al., 2006; Neo et al., 2011; Shadli et al., 2021). However, these studies specifically relate to tasks activating goal conflict and behavioral inhibition, which may differ fundamentally from the concept of trait-like behavioral inhibition observed during rest without task engagement, as we investigate here. Additionally, variability in frequency ranges measured and the electrodes analyzed (elaborated in the discussion) introduces uncertainty about the precise role of theta rhythm in anxiety and BIS.
A few earlier studies related to theta activity and anxiety at rest offer results that seem to challenge BIS theory. For example, Mizuki et al. (1984, 1989) observed that individuals with higher trait anxiety (as measured by Maudsley Personality Inventory and STAI) exhibit lower FMT (4.5–7.5 Hz) during a mental arithmetic task and do not exhibit FMT during rest. Interestingly, while FMT was absent at rest, it appeared in individuals after administering anxiolytic drugs, allowing authors to suggest that the appearance of FMT reflects relief from anxiety. Later, Suetsugi et al. (2000) confirmed these findings, reporting that non-medicated patients with generalized anxiety disorder, who initially lacked FMT, began to exhibit it following treatment, coinciding with symptom reduction. It is important to note that the EEG recordings in some of the earlier studies were captured on paper and analog tape recorders. Given the limitations of these older recording methods, it may be helpful to approach these findings and interpretations with some caution. In contrast, a sensor-level connectivity analyses have shown increased midline theta coherence during rest in individuals with generalized social anxiety disorder compared with healthy individuals (Xing et al., 2016).
Another opposing line of research, as discussed in the section on resting-state alpha, indicates that both at rest and during task presentation, BIS was associated with power at higher frequencies, such as alpha, whereas lower frequencies, such as theta, were linked to BAS (Knyazev et al., 2002, 2003, 2004). As already mentioned, the use of global theta measures and the combination of eyes-open and eyes-closed rest conditions in their analysis adds further ambiguity to BIS and theta, as well as alpha research.
Despite the controversy, the theoretical framework and most of the task-related research agree that the theta rhythm could be a valuable correlate of human anxiety, although the precise relationship, cortical region, or the frequency range of the rhythm remain unresolved.
1.3 Goals
It should be evident from this summary that the findings related to alpha and theta activity in anxiety research are highly variable and apparently contradictory. This inconsistency may stem from methodological differences (such as different questionnaires used to measure anxiety), as well as the diverse experimental and resting conditions employed. For instance, methodological limitations in earlier studies might explain the absence of FMT at rest, potentially influencing other contentious findings. In addition, the electrodes recorded and analyzed vary substantially. Some of the studies linking resting state and anxiety combine eyes-open and eyes-closed conditions in their analysis, leading to ambiguous results and interpretations.
Moreover, most research on the relationship between theta activity and BIS focuses on task-based, rather than resting-state, EEG. However, behavioral inhibition during tasks that trigger internal conflict may differ fundamentally from the concept of behavioral inhibition observed during rest, where no specific task demands are present. Therefore, gaining a comprehensive understanding of resting-state theta and alpha in anxiety could offer valuable insights not only into rest itself but also into how task-related activity relates to anxiety, serving as a reference point and refining our interpretation of these connections. Lastly, as discussed above, we suggest that one of the most significant factors contributing to inconsistent findings regarding both alpha and theta rhythms during rest may be the failure to distinguish between rhythmic oscillations and non-rhythmic components within spectral power (Caplan et al., 2015; see also Klimesch, 1999).
We address these inconsistencies in several ways. By applying a more rigorous method that reliably differentiates rhythmic from non-rhythmic power changes, we can directly test whether the power and oscillation measures match or contradict. By using a sufficiently long time period (1 minute), we aim to capture true rest, avoiding responses to rapid state transitions, but by keeping these shorter than many studies, we also avoid the onset of sleep. Moreover, to effectively measure anxiety during more and less aroused resting states, we analyze eyes-open and eyes-closed conditions separately. We also focus on theoretically grounded electrode selection. Considering that BIS and anxiety are primarily associated with anterior brain regions, we aim to address the current gap in understanding how resting-state frontal midline alpha and theta power, as well as oscillations, relate to BIS and trait anxiety. We hypothesize that the frontal midline alpha oscillations, and even the rare frontal midline theta oscillations detected at rest, might account for individual differences in trait anxiety.
In sum, we aim to (1) test the hypothesis that frontal midline alpha and theta oscillations predict BIS, (2) test the hypothesis that frontal alpha asymmetry predicts BIS, (3) compare the conventional power measurement with the oscillation measures, (4) test the hypothesis that posterior alpha, potentially related to alertness or visual inattention, also relates to BIS.
2 Methods
2.1 Design of the current study
The resting-state EEG experiment was part of a larger unpublished study investigating the role of neural oscillations in associative memory and anxiety. Resting-state EEG recording was performed after a brief break following the completion of associative memory task (not described here), anxiety and demographic questionnaires. To measure anxiety, we used Behavioral Inhibition and Activation Scale (BIS/BAS; Carver & White, 1994), State and Trait Anxiety Inventory (STAI; Spielberger et al., 1983), and Ten Item Personality Inventory (TIPI; Gosling et al., 2003).
2.2 Participants
Sixty-nine participant (all self-reported right-handed, 36 female) undergraduate students enrolled in an introductory psychology course at the University of Alberta, aged 17–51 years (mean = 20.51, SD = 4.86), participated for course credit. Data from 17 participants were excluded from analyses due to excessive amounts of artifacts in the EEG; because trial rejection was not possible considering the nature of resting-state recording, we removed participant data with more than 20% of bad electrodes or more than 20% of bad ICA components. An additional five participants’ data were excluded from BIS analysis due to incomplete BIS questionnaire data. We asked participants to report whether they had been diagnosed with anxiety or any other mental health disorders (10 out of 51 responded yes, %). However, as we did not assess clinical anxiety ourselves, we did not use these self-reports in our analysis. Importantly, no participants were excluded from the analysis based on their responses.* Data from the remaining 51 participants for the STAI and TIPI and 46 participants for BIS were used for further analysis. All participants were required to either be fluent in English or have English as their first language (learned before the age of 6 years), have normal or corrected vision, and have normal or corrected hearing. Written informed consent was obtained prior to the experiment in accordance with a University of Alberta ethical review board.
2.3 Behavioral protocol
The session took place in an electrically shielded, sound-attenuated chamber. The task was presented in E-prime (Psychology Software Tools Inc., Pittsburgh, PA). After completing the associative recognition memory task, participants were given three anxiety questionnaires, followed by a short demographic questionnaire. After the completion of all the questionnaires, participants took a short break. After the break and receiving oral instructions, participants were shown the fixation cross and were given written instructions to keep their eyes open after the beep, until they heard the next beep (1-minute eyes-open condition), and then to close their eyes and open them only after they heard the next beep (1-minute eyes-closed condition). This was repeated twice, with four beep presentations in total, constituting 4 minutes of two 1-minute eyes-open and two 1-minute eyes-closed resting cycles (Fig. 1).
(Left) Experimental design: subjects received verbal and written instructions to begin with their eyes open and to close their eyes gently upon hearing a beep. The beep was presented every 1 minute (instructing participants to open or close their eyes) with four beep presentations in total, constituting 4 minutes (2 minutes eyes-open and 2 minutes eyes-closed) EEG recording. (Right) 21 electrodes of the International 10–20 system for EEG with highlighted electrodes were used in our analysis. Namely, red represents Fz and Pz electrodes used in alpha and theta power and oscillatory multiple regressions and correlations; gray represents F3 and F4 for mid-frontal alpha asymmetry, and orange represents F7 and F8 for lateral-frontal alpha asymmetry analyzes. Note: The EGI Geodesic-Sensor-Net used in the experiment consists of 256 electrodes; a 10–20 system is provided as an easier reference.
(Left) Experimental design: subjects received verbal and written instructions to begin with their eyes open and to close their eyes gently upon hearing a beep. The beep was presented every 1 minute (instructing participants to open or close their eyes) with four beep presentations in total, constituting 4 minutes (2 minutes eyes-open and 2 minutes eyes-closed) EEG recording. (Right) 21 electrodes of the International 10–20 system for EEG with highlighted electrodes were used in our analysis. Namely, red represents Fz and Pz electrodes used in alpha and theta power and oscillatory multiple regressions and correlations; gray represents F3 and F4 for mid-frontal alpha asymmetry, and orange represents F7 and F8 for lateral-frontal alpha asymmetry analyzes. Note: The EGI Geodesic-Sensor-Net used in the experiment consists of 256 electrodes; a 10–20 system is provided as an easier reference.
The BIS/BAS questionnaire consisted of 20 questions with a 7-point Likert scale. Seven questions assess BIS (anxiety), whereas the remaining items measure Behavioral Activation System subscales (BAS-Drive, BAS-Reward, and BAS-Fun Seeking). Therefore, we used only the seven anxiety-related items (items 2, 3, 6, 9, 12, 15, 18) representing a score range from 7 to 49 (mean: 37.5). Two items from BIS (12 and 18) were reverse scored. These reversed-scored items are also those that were scored as fear-related items when the BIS scale was divided into BIS-anxiety and BIS-FFFS measuring subscales. More precisely, according to Johnson et al. (2003), J-BIS (anxiety-measuring) subscale consists of items 2, 3, 6, 9, and 15, whereas J-FFFS (fear-measuring) subscale includes items 12 (even if something bad is about to happen to me, I rarely experience fear or nervousness) and 18 (I have very few fears compared with my friends).
The STAI consisted of 40 questions total, 20 for state and 20 for trait anxiety scales, with a 4-point Likert scale. We scored and used the trait anxiety scale, which represents a score range from 20 to 80 (mean: 51.7). The TIPI consisted of 10 questions total (2 for each personality type), with a 7-point Likert scale. We only scored and used two (items 4 and 9) questions measuring neuroticism. Item 9 was reverse scored. Respectively, the score ranges from 2 to 14 (mean: 8).
2.4 Electroencephalography (EEG) recording and analyses
EEG was recorded using a high-density 256-channel Geodesic Sensor Net (Electrical Geodesics Inc., a neurodiagnostic equipment company based in Eugene, OR), amplified at a gain of 1,000 and sampled at 500 Hz. Impedances were kept below 50 k, and EEG was initially referenced to the vertex electrode (Cz). Participants were shown the recording artifacts produced by eye blinks, jaw clenching, limb movement, and other common tics to demonstrate the importance of remaining still. Data were analyzed by custom MATLAB scripts in conjunction with the open-source EEGLAB toolbox (Delorme & Makeig, 2004). Artifacts were corrected via Independent Component Analysis, implemented in EEGLAB (Jung et al., 2000), in conjunction with the ICLabel toolbox (Pion-Tonachini et al., 2019). This procedure classifies the components that capture muscle and eye-blink-related artifacts. The selection of components was based on visual inspection of spatial topographies, time courses, and power spectral characteristics of all components. ICLabel was used as an additional tool to classify eye, muscle, or heart components with a probability , but the decision to exclude components was made after visually inspecting them and other suspicious components. The components accounting for stereotyped artifacts, including eye blinks, eye movements, and muscle movements, were removed from the data. Artifacts of line noise at 60 Hz were also removed. If more than 20% of ICs or electrodes needed to be removed, the subject data were not considered for further analysis. Event latencies were corrected with a time lag correction (18 ms) due to a known hardware calibration problem identified by EGI.
2.5 Oscillation analysis
The spectral analysis informs us about the different frequencies (cycle per second, Hz) and their power (measured in ) at a given time window of the recording (Cohen, 2014). Conventional spectral analysis methods are not selective for rhythmicity of the signal. To test whether the neural activity is truly rhythmic or a non-rhythmic (non-repeating) signal, oscillations were detected and quantified using BOSC (Caplan et al., 2001; Whitten et al., 2011). The first step in the BOSC process is to calculate power using Morlet wavelet transform (Grossmann & Morlet, 1985). The subsequent steps involve determining the additional thresholds (explained below) to control for true rhythmicity. To check whether our results might generalize to more conventional spectral analysis methods, we ran parallel statistical analyses using BOSC and a measure of wavelet power (before additional thresholds are applied), which is expected to be less selective and relatively more sensitive to non-rhythmic signals. For clarity in terminology, we use “power” when referring to wavelet power and “oscillations” when referring to activity identified by BOSC after the thresholds were applied, passing its more stringent criteria for rhythmicity.
Both BOSC and power analysis were run on the continuous resting-state EEG recording. Our main interest for the single electrode analysis lies in frontal midline activity, thus the Fz is a primary electrode of interest. However, due to the alpha signal being stronger at the parietal locations and associated with changes in arousal, we decided to test the Pz electrode as well. For the frontal alpha asymmetry analysis, we selected two mid-frontal (F4–F3) and two lateral-frontal (F8–F7) electrode pairs. Thus, the spectral analysis was performed at five frontal (Fz, F3, F4, F7, F8) and one parietal (Pz) electrodes† (Fig. 1), using a Morlet wavelet transform (Grossmann & Morlet, 1985) with a window of 6 cycles and frequency sampled in 20 logarithmic steps covering the bandwidth from 1 to 32 Hz.
With BOSC analysis, the signals are only classified as rhythmic if they exceed a particular power threshold, for a given frequency, , and a particular duration threshold for each frequency. To calculate the power threshold, the method takes into account and models the background “colored noise” spectrum of EEG. The colored noise spectrum is a basic property of EEG as well as other natural autocorrelated signals (Shlesinger & West, 1988). This spectrum is then fit with linear regression in log-log space. The power threshold () is set to the percentile of the theoretical probability distribution of power values at a given frequency, and the duration threshold () is set at each frequency to three consecutive cycles (). This ensures that detected oscillations are reliably different from the background noisy component of the EEG signal and does not include increases in spectral amplitude that are non-repeating. This method is thus more selective for truly rhythmic (repeating) activity than other methods (Caplan et al., 2001; Whitten et al., 2011).
The final measure derived from the BOSC method is termed , the proportion of time during which oscillations are detected at a given frequency. For example, means oscillations at the frequency of interest were detected during 80% of the recording (illustrated in Fig. 2A). The frequency bands of interest were defined as theta, 4–8 Hz, and alpha, 9–14 Hz. The power and oscillation measures within a band were derived by averaging power or within that particular bandwidth, respectively.
These figures demonstrate the detection of oscillatory activity (A) and its predictive relationship with BIS scores (B). (A) Example of the raw EEG traces recorded during eyes open rest at electrode Fz, with oscillations detected by BOSC in theta (in red) and alpha (in yellow) ranges in one high-BIS (top, score 48) and one low-BIS (bottom, score 23) participant. (B) The representation of the multiple regression model predicting BIS scores with theta and alpha oscillatory activity () measures recorded at frontal midline (Fz) electrode during eyes open and eyes closed resting states. Red bars represent beta coefficients for theta oscillations predicting BIS, and yellow bars represent beta coefficients for alpha oscillations predicting BIS.
These figures demonstrate the detection of oscillatory activity (A) and its predictive relationship with BIS scores (B). (A) Example of the raw EEG traces recorded during eyes open rest at electrode Fz, with oscillations detected by BOSC in theta (in red) and alpha (in yellow) ranges in one high-BIS (top, score 48) and one low-BIS (bottom, score 23) participant. (B) The representation of the multiple regression model predicting BIS scores with theta and alpha oscillatory activity () measures recorded at frontal midline (Fz) electrode during eyes open and eyes closed resting states. Red bars represent beta coefficients for theta oscillations predicting BIS, and yellow bars represent beta coefficients for alpha oscillations predicting BIS.
All trials were time locked to the onset of the beep presentation, and of all frequencies at each EO and EC trial was calculated using the time window of 1–60,000 ms latency post-stimulus, minus shoulders (12,000 ms) to avoid edge artifacts. This left 48,000 ms of each EO or EC phase. This shoulder subtraction also ensures that we are not measuring EEG activity right after the beep presentation to avoid the auditory signal processing and oculomotor processes (and artifact). Because the defining frequencies for a rhythm band (like theta and alpha) vary from study to study, it is also important to examine each frequency and its relationships with anxiety scores individually. To check the robustness of our frequency bands, we also examined our correlation analyses at all frequencies sampled over the 1–32 Hz range.
To measure frontal alpha asymmetry, the laterality quotient (LQ) was calculated by subtracting left frontal alpha power/oscillations from right frontal alpha power/oscillations—LQ = (Reznik & Allen, 2018). Thus, positive LQ values represent relatively higher right alpha power/oscillation (rightward asymmetry), whereas negative LQ values represent relatively higher left alpha power/oscillation (leftward asymmetry) (Metzen et al., 2022). If one assumes the inverse relationship between alpha power and cortical activity, positive LQ would reflect higher left cortical activity, whereas negative LQ would reflect higher right cortical activity. As already mentioned, we calculated alpha asymmetry at two pairs of electrodes, mid-frontal electrode pair (F4–F3) and lateral-frontal (F8–F7).
2.6 Statistical analysis
All statistical analyses were carried out using MATLAB and Statistic Toolbox Release 2020a (The MathWorks, Inc., Natick, Massachusetts, United States), and JASP (Jasp Team, University of Amsterdam, Netherlands) on power and for oscillations at the corresponding electrodes and time windows.
First, we conducted two multiple regression analyses, one at Fz and one at Pz, to explore whether BIS scores are predicted by alpha and theta oscillations during eyes-open and eyes-closed resting states. We then conducted two similar regressions (one at Fz and one at Pz) to test whether BIS scores are predicted by alpha and theta power during eyes-open and eyes-closed resting states. Next, we ran two additional regressions to test whether BIS scores are predicted by frontal alpha asymmetry (one with power and one with oscillations) at mid-frontal and lateral-frontal electrodes recorded during eyes-open and eyes-closed resting states. Finally, as an exploratory analysis, we conducted a post hoc regression combining Fz and Pz to test whether BIS scores are predicted by alpha and theta oscillations recorded during eyes-open and eyes-closed resting states, with a similar analysis for power measures. We combined Fz and Pz predictors in the regression to determine their unique contributions to predicting trait anxiety when considered simultaneously. We also conducted Pearson correlations between alpha and theta oscillations (and power) recorded at Fz and Pz electrodes during eyes open and eyes closed rest, and BIS, STAI and TIPI neuroticism scores. Correlations were computed both for the average values within each frequency band (theta and alpha) and for each individual sampled frequency (1 - 32 Hz). We also conducted a post hoc reliability analysis using bootstrapping with sub-sampled data segments ranging from 10% to 100% of the total duration.
Statistical tests are reported with both Classical and Bayesian approaches. Significance is assessed with but p values near-threshold are interpreted with caution. Bayes Factors are considered to provide support for the null hypothesis if or for the hypothesis if (Kass & Raftery, 1995). All the BFs are reported in the brackets in the Results section and on tables, but only a subset is highlighted in the text.
3 Results
3.1 BIS and resting-state neural oscillation and power model at Fz
Multiple regression for BIS and resting-state oscillations (the BOSC measure) model at Fz electrode is summarized in Table 1 and also visualized in Figure 2B. The model with four predictors (Fz Eyes Open, Eyes Closed, Eyes Open, Eyes Closed) was significant, with adjusted , , , . The (, , , ) and (, , , ) recorded during eyes-open rest are reliably different from zero. These findings suggest that higher mid-frontal theta oscillations and lower mid-frontal alpha oscillations recorded during eyes-open rest predict trait anxiety, as measured with behavioral inhibition scale.
Regression model for mid-frontal (Fz) and mid-parietal (Pz) resting-state neural oscillations (left) and power (right) predicting BIS scores.
Fz . | Oscillations . | Power . | ||||
---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . |
Eyes Open | 0.33 | 6.14 | 0.93 | 3.6 | ||
Eyes Closed | -0.05 | -0.32 | 1.13 | -0.52 | -1.21 | 1.36 |
Eyes Open | -0.46 | 16.65 | -0.94 | 7.07 | ||
Eyes Closed | 0.3 | 1.82 | 3.22 | 0.49 | 1.97 | 2.32 |
Fz . | Oscillations . | Power . | ||||
---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . |
Eyes Open | 0.33 | 6.14 | 0.93 | 3.6 | ||
Eyes Closed | -0.05 | -0.32 | 1.13 | -0.52 | -1.21 | 1.36 |
Eyes Open | -0.46 | 16.65 | -0.94 | 7.07 | ||
Eyes Closed | 0.3 | 1.82 | 3.22 | 0.49 | 1.97 | 2.32 |
Pz . | Oscillations . | Power . | ||||
---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . |
Eyes Open | 0.02 | 0.09 | 0.6 | 0.43 | 0.66 | 0.91 |
Eyes Closed | 0.14 | 0.79 | 0.81 | 0.25 | 0.39 | 1.08 |
Eyes Open | -0.4 | 2.44 | -1.04 | 2.75 | ||
Eyes Closed | 0.19 | 1.13 | 0.94 | 0.45 | 1.44 | 1.54 |
Pz . | Oscillations . | Power . | ||||
---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . |
Eyes Open | 0.02 | 0.09 | 0.6 | 0.43 | 0.66 | 0.91 |
Eyes Closed | 0.14 | 0.79 | 0.81 | 0.25 | 0.39 | 1.08 |
Eyes Open | -0.4 | 2.44 | -1.04 | 2.75 | ||
Eyes Closed | 0.19 | 1.13 | 0.94 | 0.45 | 1.44 | 1.54 |
Each predictor for oscillation analysis is the proportion of oscillations () of band-average neural activity, and each predictor in power analysis is the power of a band-average neural activity recorded at rest during eyes-open or eyes-closed conditions, respectively. ; . Overall, the Fz oscillation model explained 26% of the variance (), the Fz power model explained 21% of the variance (), the Pz oscillation model explained 12% of the variance (), and the Pz power model explained 14% of the variance ().
Multiple regression for BIS and resting-state power model is also summarized in Table 1. Similar to the multiple regression with oscillation measures, the model with all four power predictors was also significant, with adjusted , , , . The main predictors are the same in the power model as we found in the oscillation model. More precisely, the (, , , ) and (, , , ) recorded during eyes-open rest are reliably different from zero, also supported by Bayesian statistics. The fact that BOSC and power measures both predict the BIS scores with a similar pattern strengthens the claim that resting-state neural activity holds valuable information for anxiety processes. Moreover, BOSC measures suggest that there is legitimate rhythmicity underlying the effect.
3.2 BIS and resting-state neural oscillation and power models at Pz
Multiple regression for BIS and Pz resting-state model is summarized in Table 1. The model with four oscillatory predictors (Pz Eyes Open, Eyes Closed, Eyes Open, Eyes Closed) fell just short of significance, with adjusted , , , and . The only significant predictor from the Pz oscillation model was alpha recorded during eyes-open rest (, , , ). Contrary to the Pz oscillation model, the similar model with all four power predictors was significant with adjusted , , , and . As in the oscillation model, the only significant predictor in the power model was alpha recorded during eyes-open rest (, , , ). However, the Bayes Factors suggest there is not quite enough evidence for eyes-open alpha oscillation () or power () recorded at the Pz electrode to be included in the model, potentially suggesting that we may be recording more of classic “resting visual alpha” at Pz, leaving the frontal midline neural activity exclusively important in predicting anxiety.
3.3 BIS and resting-state neural oscillation and power models at Fz and Pz combined
Having observed that both Fz and Pz electrodes shared the significant predictors, we ran a post hoc, combined regression with all eight predictors (Fz: Eyes Open, Eyes Closed, Eyes Open, Eyes Closed, and Pz: Eyes Open, Eyes Closed, Eyes Open, Eyes Closed), summarized in Table 2. The model was significant, with adjusted , , , and . Interestingly, the only significant predictor in this combined analysis proved to be theta oscillations recorded at Fz electrode during eyes-open rest (, , , ).
The combined regression model for mid-frontal (Fz) and mid-parietal (Pz) resting-state neural oscillations (left) and power (right) predicting BIS scores.
. | Oscillation . | Power . | |||||
---|---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . | |
Fz | Eyes Open | 0.39 | 3.26 | 1.07 | 2.43 | ||
Eyes Closed | -0.2 | -0.88 | 0.69 | -1.41 | 1.86 | ||
Eyes Open | -0.45 | -1.31 | 2.03 | -0.54 | -1.0 | 2.16 | |
Eyes Closed | 0.42 | 1.34 | 1.17 | 0.68 | 1.41 | 1.5 | |
Pz | Eyes Open | -0.26 | -1.13 | 0.87 | -0.58 | -0.75 | 1.11 |
Eyes Closed | 0.3 | 1.42 | 0.92 | 1.7 | 2.08 | 1.92 | |
Eyes Open | 0.07 | 0.19 | 1.11 | -0.57 | -0.8 | 1.81 | |
Eyes Closed | -0.22 | -0.73 | 0.78 | -0.33 | -0.54 | 1.3 |
. | Oscillation . | Power . | |||||
---|---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . | |
Fz | Eyes Open | 0.39 | 3.26 | 1.07 | 2.43 | ||
Eyes Closed | -0.2 | -0.88 | 0.69 | -1.41 | 1.86 | ||
Eyes Open | -0.45 | -1.31 | 2.03 | -0.54 | -1.0 | 2.16 | |
Eyes Closed | 0.42 | 1.34 | 1.17 | 0.68 | 1.41 | 1.5 | |
Pz | Eyes Open | -0.26 | -1.13 | 0.87 | -0.58 | -0.75 | 1.11 |
Eyes Closed | 0.3 | 1.42 | 0.92 | 1.7 | 2.08 | 1.92 | |
Eyes Open | 0.07 | 0.19 | 1.11 | -0.57 | -0.8 | 1.81 | |
Eyes Closed | -0.22 | -0.73 | 0.78 | -0.33 | -0.54 | 1.3 |
Each predictor for oscillation analysis is the proportion of oscillations () for band-average neural activity, and each predictor in power analysis is the power of a band-average neural activity recorded at rest during eyes-open or eyes-closed conditions, respectively. ; . Overall, the Fz oscillation model explained 26% of the variance (), the Fz power model explained 21% of the variance (), the Pz oscillation model explained 12% of the variance (), and the Pz power model explained 14% of the variance ().
For power, similarly, the combined Fz and Pz regression model with eight predictors was significant, with adjusted , , , and . Interestingly, in the combined power model, the prediction was driven by two significant predictors (in contrast to the one significant predictor in the oscillation model), which are theta oscillations recorded at Fz during eyes-open rest (, , , ) and during eyes-closed rest (, , , ). Surprisingly, neither of these predictors were conclusive with Bayesian analysis.
This approach allowed us to see whether each region’s activity remained significant when controlling for the other. Our analysis showed that although parietal alpha power predicted BIS scores in isolation, it did not remain significant when frontal measures were included. This suggests that the relationship of the parietal alpha with BIS may be weaker or more context dependent relative to frontal midline activity. This distinction clarifies the dominant role of frontal midline rhythms in anxiety-related inhibition.
3.4 Reliability analysis
To address concerns regarding the reliability of the findings based on the 2-minute resting-state EEG epochs for each condition, we conducted a robustness check using bootstrapping with sub-sampled data segments ranging from 10% to 100% of the total recording duration. For each percentage of data used, we performed 100 bootstrapped linear regressions to predict BIS scores using mid-frontal (Fz) and mid-parietal (Pz) theta and alpha oscillations recorded during eyes-open (EO) and eyes-closed (EC) conditions.
As shown in Figure 3, the average beta values plateaued as the percentage of data increased. The results indicate that using between 30% and 50% of the total EEG recording duration provides reliable and stable predictions of trait anxiety. Beyond this range, increasing the data duration did not significantly improve the predictive power of the models.
Reliability analysis was performed using bootstrapping with sub-sampled EEG data segments ranging from 10% to 100% of the total recording duration (48 seconds for oscillation and 60 seconds for power analysis). For each percentage window of data used, we performed 100 bootstrapped linear regressions to predict BIS scores using mid-frontal (Fz) and mid-parietal (Pz) theta and alpha oscillations (A and C) and power (B and D) recorded during eyes-open (EO) and eyes-closed (EC) conditions. The error bars represent the 95% confidence intervals of the standardized beta coefficients.
Reliability analysis was performed using bootstrapping with sub-sampled EEG data segments ranging from 10% to 100% of the total recording duration (48 seconds for oscillation and 60 seconds for power analysis). For each percentage window of data used, we performed 100 bootstrapped linear regressions to predict BIS scores using mid-frontal (Fz) and mid-parietal (Pz) theta and alpha oscillations (A and C) and power (B and D) recorded during eyes-open (EO) and eyes-closed (EC) conditions. The error bars represent the 95% confidence intervals of the standardized beta coefficients.
These findings suggest that the original 1-minute (total of 2-minute per condition) EEG epochs used in the main analysis are sufficient to produce robust results. The bootstrapping approach further supports the stability of the findings, reducing the likelihood that the observed relationships between mid-frontal oscillations and BIS scores are due to chance or variability in the data.
3.5 Correlations between eyes-open resting-state oscillations and anxiety questionnaires
Although BIS scores are not reliably predicted by oscillations recorded at the Pz electrode, the negative correlations of all questionnaire scores with eyes-open resting states, both for frontal (Fz) (Fig. 4A) and parietal (Pz) electrodes, are significant in the range (Table 3). Interestingly, oscillations are not correlated with anxiety questionnaires; only a marginal but non-significant correlation is found between frontal midline recorded during eyes-open rest and BIS scores () (Fig. 4B).
Scatter plots for Pearson correlations () between eyes-open alpha (A) and theta (B) oscillations recorded at Fz electrode and BIS scores. Red dots represent individuals with self-reported anxiety diagnosis.
Scatter plots for Pearson correlations () between eyes-open alpha (A) and theta (B) oscillations recorded at Fz electrode and BIS scores. Red dots represent individuals with self-reported anxiety diagnosis.
Pearson correlations between alpha and theta oscillations recorded at Fz (left) and Pz (right) electrodes during eyes-open and eyes-closed rest, and BIS, STAI Trait, and TIPI neuroticism scores.
. | Fz . | Pz . | |||||||
---|---|---|---|---|---|---|---|---|---|
. | df . | EO . | EC . | EO . | EC . | EO . | EC . | EO . | EC . |
BIS | 45 | 0.25 | 0.25 | 0.12 | 0.22 | -0.17 | 0.17 | ||
STAI Trait | 50 | -0.002 | 0.03 | 0.05 | 0.15 | -0.15 | 0.06 | ||
TIPI | 50 | 0.02 | -0.02 | 0.05 | -0.25 | 0.16 | -0.16 | 0.15 |
. | Fz . | Pz . | |||||||
---|---|---|---|---|---|---|---|---|---|
. | df . | EO . | EC . | EO . | EC . | EO . | EC . | EO . | EC . |
BIS | 45 | 0.25 | 0.25 | 0.12 | 0.22 | -0.17 | 0.17 | ||
STAI Trait | 50 | -0.002 | 0.03 | 0.05 | 0.15 | -0.15 | 0.06 | ||
TIPI | 50 | 0.02 | -0.02 | 0.05 | -0.25 | 0.16 | -0.16 | 0.15 |
, .
We also plotted correlation values as functions of all frequencies at electrodes Fz and Pz (Fig. 5). Generally, the correlation results at all frequencies analyzed confirmed the results of the frequency band-specific (band-average) analyses. However, it also becomes evident that the effects are not limited to theta and alpha ranges. Frequencies in the range of 1–8 Hz (which includes theta band) are positively associated with BIS, whereas frequencies in the range of 8–32 Hz (which includes alpha band) are more negatively associated with the BIS scores. This type of analysis highlights a pitfall of arbitrarily defining frequency ranges for rhythms such as theta and alpha. For example, although theta is traditionally considered to span 4–8 Hz, the correlation at all sampled frequencies demonstrated in Figure 5A reveals that the 8 Hz frequency aligns more closely with the negative correlation trend of the “alpha peak” rather than a positive correlation pattern of the “theta peak.” Furthermore, the marginal but non-significant nature of the theta-band average correlation with BIS scores () suggests that removing 8 Hz from the theta range could potentially reveal an underlying significant correlation with more frequency specificity. Even though we refrain from making such adjustments post hoc, this is an important consideration for future research.
Pearson correlation () for BIS anxiety scores correlating with resting-state neural oscillations () recorded at Fz and Pz electrodes during eyes-open (A) and eyes-closed (B) conditions plotted for all frequencies. The dashed lines denote the significance thresholds ().
Pearson correlation () for BIS anxiety scores correlating with resting-state neural oscillations () recorded at Fz and Pz electrodes during eyes-open (A) and eyes-closed (B) conditions plotted for all frequencies. The dashed lines denote the significance thresholds ().
3.6 Alpha asymmetry and anxiety
The multiple regression for BIS and resting-state alpha oscillatory asymmetry model is summarized in Table 4 and also visualized in Figure 6. The model with all four predictors (F4–F3 asymmetry EO, F4–F3 asymmetry EC, F8–F7 asymmetry EO, and F8–F7 asymmetry EC) was not significant (), nor were any of the individual predictors. The same analysis on power instead of (Table 4) was also not significant (), nor were any of the particular predictors. These findings suggest that frontal alpha asymmetry is not a reliable predictor of BIS scores.
The representation of the multiple regression model predicting BIS scores with lateral-frontal (F8–F7) and mid-frontal (F4–F3) alpha oscillatory () asymmetry measured during eyes-open and eyes-closed resting states. Orange bars represent beta coefficients for lateral-frontal (F8–F7) alpha oscillatory asymmetry measures predicting BIS, and gray bars represent beta coefficients of mid-frontal (F4–F3) alpha oscillatory asymmetry measures predicting BIS scores.
The representation of the multiple regression model predicting BIS scores with lateral-frontal (F8–F7) and mid-frontal (F4–F3) alpha oscillatory () asymmetry measured during eyes-open and eyes-closed resting states. Orange bars represent beta coefficients for lateral-frontal (F8–F7) alpha oscillatory asymmetry measures predicting BIS, and gray bars represent beta coefficients of mid-frontal (F4–F3) alpha oscillatory asymmetry measures predicting BIS scores.
The regression model for mid-frontal (F4–F3) and lateral-frontal (F8–F7) resting-state neural oscillatory activity (left) and power asymmetry (right) predicting BIS scores.
. | Oscillations . | Power . | ||||
---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . |
F4–F3 asym EO | -0.2 | -1.14 | 0.59 | -0.61 | -1.28 | 0.3 |
F4–F3 asym EC | 0.18 | 1.0 | 0.55 | 0.56 | 1.16 | 0.29 |
F8–F7 asym EO | -0.3 | -1.77 | 1.66 | -0.02 | -0.04 | 0.24 |
F8–F7 asym EC | 0.06 | 0.35 | 0.45 | 0.08 | 0.15 | 0.26 |
. | Oscillations . | Power . | ||||
---|---|---|---|---|---|---|
Predictors . | beta . | t . | . | beta . | t . | . |
F4–F3 asym EO | -0.2 | -1.14 | 0.59 | -0.61 | -1.28 | 0.3 |
F4–F3 asym EC | 0.18 | 1.0 | 0.55 | 0.56 | 1.16 | 0.29 |
F8–F7 asym EO | -0.3 | -1.77 | 1.66 | -0.02 | -0.04 | 0.24 |
F8–F7 asym EC | 0.06 | 0.35 | 0.45 | 0.08 | 0.15 | 0.26 |
Each predictor represents the band-average neural activity ( or power) for the left electrode (F7 or F3) subtracted from the right electrode (F8 or F4), recorded during eyes-open or eyes-closed resting states. ; .
Following the Coan and Allen (2003) logic about BIS system complexity, as an exploratory step, we also looked into different BIS subscales split according to the model suggested by Johnson et al. (2003). To recap, it was suggested that if more right hemisphere alpha activity is triggered by the BIS threat-related subsystem, and not by anxiety-related processes, the relationship between full BIS scale and right frontal alpha activity would be sample- and situation-dependent Coan and Allen (2003).
To directly test this idea, we correlated scores from different BIS subscales with our mid-frontal (F4–F3) and lateral-frontal (F8–F7) asymmetry measures. We found greater lateral-frontal (F8–F7) oscillatory asymmetry during eyes-open rest to be negatively correlated with overall BIS scores (, ). However, after splitting the BIS scale, this correlation only stayed significant for the J-FFFS BIS subscales (, ), but not for the J-BIS one (, ). This dissociation suggests that individuals with more activated fear-related processes, but not exclusively anxiety-related processes, had lower rightward alpha asymmetry during eyes-open rest (Table 5). The relationship was not found between power asymmetry measures and BIS scales.
Pearson correlation between alpha oscillatory () asymmetry recorded at lateral-frontal (F8–F7) & mid-frontal (F4–F3) electrode pairs during eyes-open and eyes-closed rest with BIS scale and subscale scores.
. | Oscillations . | . | . |
---|---|---|---|
BIS | F7–F8 Asym EO | * | |
F7–F8 Asym EC | |||
F4–F3 Asym EO | |||
F4–F3 Asym EC | |||
J_BIS | F7–F8 Asym EO | ||
F7–F8 Asym EC | |||
F4–F3 Asym EO | |||
F4–F3 Asym EC | |||
J_FFFS | F7–F8 Asym EO | * | |
F7–F8 Asym EC | |||
F4–F3 Asym EO | |||
F4–F3 Asym EC |
. | Oscillations . | . | . |
---|---|---|---|
BIS | F7–F8 Asym EO | * | |
F7–F8 Asym EC | |||
F4–F3 Asym EO | |||
F4–F3 Asym EC | |||
J_BIS | F7–F8 Asym EO | ||
F7–F8 Asym EC | |||
F4–F3 Asym EO | |||
F4–F3 Asym EC | |||
J_FFFS | F7–F8 Asym EO | * | |
F7–F8 Asym EC | |||
F4–F3 Asym EO | |||
F4–F3 Asym EC |
Denotes .
4 Discussion
Taking an individual differences approach, we tested whether neural activity recorded at the single frontal midline (Fz) and parietal midline (Pz) electrodes would predict behavioral inhibition as a measure of trait anxiety. We found that higher BIS scores were associated with increased frontal midline theta and decreased frontal midline alpha oscillations during eyes-open rest. Additionally, resting-state parietal midline alpha power, but not oscillations, explained less variance in BIS scores (with inconclusive Bayes factors) than the frontal midline theta and alpha power and oscillations. Further, we discovered that frontal alpha asymmetry (power or oscillations) did not predict BIS scores. Notably, correlation analysis showed that frontal midline oscillatory alpha asymmetry during eyes-open rest was negatively correlated with the fear-related BIS subscale (J-FFFS), and not with the anxiety-related subscale (J-BIS).
4.1 Frontal midline theta and anxiety
Associating higher frontal midline theta with increased trait anxiety and the behavioral inhibition system aligns with Gray and McNaughton’s septohippocampal theory (Gray, 1982; Gray & McNaughton, 2000), extending its application to resting-state EEG activity. Earlier task-related studies found that human hippocampal “theta” frequency at right frontal electrodes (F8) in the range of 4–12 Hz (named goal-conflict-specific EEG rhythmicity (GCSR)) was positively correlated with neuroticism and trait anxiety and was reduced by all key anxiolytic drugs (Neo et al., 2011; Shadli et al., 2020, 2021). Shadli et al. (2021) also found that goal-conflict-related “theta” rhythm was higher in individuals with clinical anxiety levels. However, despite its theoretical significance, the terminology of referring to GCSR as “theta” rhythm warrants some consideration, given that its range (4–12 Hz) extends beyond the traditionally defined human theta band (4–8 Hz).
In the more classical view of human theta rhythms, meta-analysis by Cavanagh and Shackman (2015) concluded that frontal midline (Fz) theta (4–8 Hz) oscillations reflect increased cognitive control and relate to greater trait anxiety and avoidant behavior during goal-conflict arousing tasks. However, they could not find any clinical significance of this neural marker. In contrast, Cavanagh et al. (2017) found that error- and conflict-related neural oscillations in the theta band can reliably differentiate individuals with generalized anxiety disorder from healthy controls.
BIS theory also suggests that unresolved goal conflict subjectively manifests as anxious rumination (Corr, 2008; Gray & McNaughton, 2000). Andersen et al. (2009) experimentally linked this idea to theta activity, showing that personally meaningful goal conflict, potentially triggering anxious rumination, indeed enhanced the scalp-wide EEG theta (4–6 and 6–8 Hz) power and coherence, compared with less personally meaningful, nominal ruminations. Surprisingly, this increase was not associated with participants’ BIS levels. On the other hand, Moore et al. (2006) found that scalp-wide theta power and coherence effectively distinguished between high and low BIS participants. However, in this study, goal conflict only increased theta coherence across various regions, whereas theta power showed a non-significant trend, suggesting that the differentiation in BIS levels was primarily driven by theta coherence measures, and not power. More interestingly, despite the hypothesized frontal midline theta activation, neither study observed region-specific increases in theta power; significant results only emerged when scalp-wide averages were analyzed.
These studies collectively suggest that the septohippocampal system is activated during goal-conflict tasks in humans. Extending this perspective, our findings suggest that in individuals with high trait anxiety, this system may be more generally activated. It appears to remain in a heightened “ready” state to inhibit behaviors fraught with potential negative outcomes, even during rest. Additionally, it is arguably more plausible to detect periods of internal anxious rumination when individuals are resting rather than engaging in goal-conflict tasks. We speculate that individuals with higher BIS may be more prone to experiencing internalized anxious rumination during rest, and that increased FMT could reflect such intrusive thoughts and internalized scenarios of goal conflict. From this perspective, our finding of an association between resting-state theta rhythm and BIS scores might indicate a stronger link to these cognitive processes compared with the null correlation between task-induced anxious rumination and BIS reported by Andersen et al. (2009).
Our findings seem at odds with Caplan et al. (2015), who used the same BOSC method but found virtually no evidence of theta oscillations during rest. So, what accounts for the robust effects in the theta band at rest in our study? A key difference between the two studies lies in the duration of the resting-state recordings. In the Caplan et al. (2015) study, each eyes-open or eyes-closed segment lasted only 5 seconds. In contrast, we used much longer 1-minute segments for each condition. This extended duration may have given us a better opportunity to capture the kinds of internalized cognitive and emotional processes, such as anxiety-related rumination, that we associate with FMT activity. This aligns with the idea that when people ruminate or dwell on anxious thoughts, those thoughts typically unfold over a longer time frame. It is unlikely that such internalized processes would arise and stabilize within just 5 seconds. By using longer recording segments, we may have provided the necessary conditions for these ruminative thoughts and the associated theta activity to emerge. If frontal theta reflects prolonged internal processes such as rumination, then short recordings simply would not allow enough time for this signal to manifest. Taken together, our findings suggest that longer resting-state recordings may be essential to accurately detect theta activity related to internalized thought processes such as rumination.
Our results also differ from older studies that reported no frontal midline theta activity during rest and lower FMT in more anxious individuals performing mental arithmetic tasks (Mizuki et al., 1984; Suetsugi et al., 2000). As mentioned in the introduction, some of these earlier studies relied on EEG recordings captured on paper and analog tape recorders, which likely limited the detection of finer neural patterns, prompting us to interpret their findings carefully. Moreover, in these studies, anxiety was not measured with behavioral inhibition scales but with the Maudsley Personality Inventory and STAI. To our knowledge, this is the first paper analyzing resting frontal midline theta with behavioral inhibition scales and methods that are selective for not just spectral power but also rhythmicity.
Clinical Relevance. In a comprehensive review, McLoughlin et al. (2022) found that FMT is consistently elevated in individuals with certain psychiatric conditions, including generalized anxiety disorder (GAD), but primarily during tasks involving uncertainty or cognitive conflict rather than rest. By contrast, Kopanska et al. (2024), reviewing both resting and task-based EEG rhythms in neuropsychiatric disorders, found no consistent evidence linking theta oscillations with anxiety disorders.
Direct clinical research examining anxiety-related FMT during resting-state EEG remains limited. One notable exception is Xing et al. (2016) who found increased frontal midline theta coherence in individuals with generalized social anxiety disorder during rest, a pattern that aligns with our findings. Theta rhythms have also been linked to other mental health disorders, including depression, a common comorbidity of anxiety disorders. Gold et al. (2013) examined FMT and frontal alpha asymmetry in a clinically depressed sample and found that although reliability was adequate, FMT was only weakly correlated with anxiety measures. Zhang et al. (2022) also explored whether resting theta rhythm could serve as a useful biomarker for distinguishing between anxiety and depressive disorders, and found that theta power (4–7 Hz) was significantly stronger in individuals with depression than in those with anxiety, though this effect was limited to certain electrodes (e.g., F3, O2, T3, P3, P4, FP1, FP2, and F8). Frontal theta activity has additionally been explored as a therapeutic target. For example, a study by Li et al. (2022) demonstrated that continuous theta-burst stimulation (cTBS) over the right dorsolateral prefrontal cortex improved anxiety symptoms in patients with generalized anxiety disorder, accompanied by increases in alpha oscillations.
Collectively, these findings suggest that although theta activity may have transdiagnostic relevance, its relationship with clinical anxiety remains debated and may depend on factors such as experimental conditions, diagnostic category, and electrode site. Notably, most clinical studies rely on anxiety symptom scales, which differ substantially from the behavioral inhibition scale used here, as it targets a specific anxiety-related process (Gray & McNaughton, 2000; Shadli et al., 2021). Our findings complement this body of work by demonstrating that FMT is also elevated during eyes-open resting-state EEG in individuals with high trait anxiety, a known risk factor for developing psychiatric disorders. This supports the idea that FMT may serve not only as a clinical biomarker but also as a potential vulnerability marker, thereby extending the relevance of clinical observations to a resting, task-free context and highlighting its potential utility in early identification of individuals at risk for anxiety-related disorders.
4.2 Frontal midline alpha and anxiety
Measuring frontal midline alpha during rest is not common in EEG research. As far as we are aware, this is the first report of frontal midline alpha being investigated and associated with the behavioral inhibition system. Due to the higher popularity of the approach-withdrawal motivational model and its association with frontal alpha asymmetry, while there has been a large amount of research done in the alpha activity and BIS domain, researchers almost always use asymmetry coefficients in their analysis and disregard the single-electrode information. However, our findings clearly suggest that alpha at midline electrodes also convey valuable information and explain a substantial part of the variance in BIS scores. As mentioned in the introduction, there might be good reasons why earlier studies favored asymmetry measures over single-electrode ones. When subtracting right hemisphere neural activity from the left, the power measurement is calibrated, and some of the intrinsic variability of the EEG signal is corrected. However, analysis methods such as BOSC, which are calibrated using the signal’s own background power spectrum, make single-electrode analyses more feasible and comparable, and potentially more accurate than the right–left asymmetry measure.
Consider, now, the alpha “idling” hypothesis. Just as greater occipital alpha is robustly associated with the inhibition of task-irrelevant visual stimuli (Foxe & Snyder, 2011), it is plausible that greater frontal midline alpha reflects the inhibition of processes governed by the frontal brain regions, such as rumination, anxiety, cognitive control, and emotion regulation. Consequently, the observation that lower frontal midline alpha predicts behavioral inhibition or trait anxiety could be attributed to increased cortical activation in anterior brain regions involved in inhibition and cognitive control, such as prefrontal cortex and anterior cingulate cortex (Bush et al., 2000; Diamond, 2013). Further, following the assumption that alpha power is an inverse measure of cortical activation (Davidson, 1992), lower alpha power in these regions suggests their greater activation. However, these explanations should be interpreted with caution due to the limited spatial resolution of EEG and the aforementioned arguments against the idea that alpha activity is always an inverse of cortical activation.
Although posterior alpha activity has long been associated with arousal and alertness (Klimesch et al., 2007), our findings indicate that only parietal-recorded eyes-open alpha power predicted BIS scores, with no significant relationship found for oscillations. Furthermore, Bayesian analyses did not support the alpha power and BIS results, suggesting that the evidence for such a relationship is weak and insufficient to make a confident claim. Moreover, in a post hoc analysis combining Fz and Pz electrode predictors in one regression model, only frontal midline theta (both power and oscillation) emerged as the primary predictor of BIS, out of all eight predictors. Bayesian analysis further confirmed that frontal midline theta oscillations were the only well-supported predictor.
Given the statistically weak results for the parietal electrode power measures, it is possible that the Pz eyes-open alpha recorded in our study reflects classic visual alpha activity (Klimesch, 1999), which we speculate is less closely tied to anxiety-related neural mechanisms. Although speculative, it is also possible that rhythmicity differentiates the functional roles of frontal and posterior alpha in relation to anxiety. Anxiety-related frontal alpha and theta may be driven by rhythmic processes, whereas posterior alpha could reflect non-rhythmic power, contributing to the variability in results across studies. This distinction could help explain why global alpha measures, which collapse both electrodes and signal types, have produced inconsistent findings in the literature. Separating these components, as done with the BOSC method, offers a more nuanced understanding and could clarify the distinct contributions of oscillatory and non-oscillatory activity to anxiety-related neural dynamics. Future studies could further test this by specifically isolating rhythmic components to better understand their distinct role in anxiety-related neural activity.
4.3 Frontal alpha asymmetry and anxiety
Evidence also suggests that the right inferior frontal cortex is important in behavioral inhibition (Aron et al., 2014). The assumption that alpha power inversely reflects cortical activation would imply that reduced alpha oscillations at the right frontal electrodes would indicate heightened activity in this area. Yet, our frontal alpha asymmetry analysis did not support this line of reasoning. We examined frontal alpha asymmetry at two electrode pairs but found that asymmetry at neither location sufficiently explained individual variability in BIS scores. These null results align with the conclusions of a recent systematic review by Firth et al. (2024), who found that among the reviewed studies, most reported no consistent relationship between resting frontal alpha asymmetry and behavioral inhibition. This convergence of evidence suggests that frontal alpha asymmetry may not be a reliable neural correlate of BIS sensitivity in resting-state EEG.
Although the overall regression was not significant, there was a weak trend suggesting that lateral-frontal (F8–F7) alpha asymmetry during eyes-open rest might predict BIS scores. However, an exploratory correlational analysis, which involved dividing the BIS scale into two subscales, revealed that only the FFFS-related subscales was significantly associated with lateral-frontal (F8–F7) asymmetry. This finding may help explain the mixed results often reported in the literature. The updated septohippocampal theory of anxiety distinguishes between two sub-processes within the avoidance/withdrawal dimension (previously grouped under BIS). Since then, some researchers have suggested that the BIS scale does not align well with this updated theory and should be revised to clearly separate FFFS (fear-related) processes from BIS (anxiety-related) processes (Heym et al., 2008; Johnson et al., 2003). Given that earlier studies did not consider such a distinction, the inconsistency in findings could reflect the differential representation of these two constructs, with alpha asymmetry measures being differently associated with FFFS and BIS. This was evident in our exploratory analysis, where only the FFFS subscale, but not the BIS subscale, significantly correlated with alpha asymmetry recorded at lateral-frontal (F8–F7) electrodes.
Another possibility is that previous studies, using conventional power measures, may not have captured the rhythmic subset of activity showing frontal alpha asymmetry. If only non-rhythmic alpha asymmetry is associated with BIS or FFFS, but the rhythmic activity is not, these studies may have detected the true effect only some of the time, contributing to the mixed results. This may be addressed in the future by using both BOSC (or similar methods) and conventional power measures when studying frontal alpha asymmetry, as we did here. Although neither power- nor oscillatory-asymmetry reliably predicted BIS scores in our study, comparing these methods in future resting-state EEG analyses could provide additional clarity on the true nature of the effect.
4.4 Resting state: eyes-open versus eyes-closed rest
Our findings also reveal significant differences in resting-state power and oscillatory activity between eyes-open and eyes-closed conditions, reinforcing the idea that these two states are functionally distinct and should not be conflated in analyses. Although alpha-rhythm researchers since Berger (1929) have known that at least posterior alpha activity changes drastically when the participant opens or closes their eyes, many studies of resting-state activity still overlook these distinctions, sometimes even combining analysis of both states or using one as a baseline for the other. For instance, Knyazev et al. (2003) reported an opposite pattern, where relative alpha power was positively associated with behavioral inhibition (i.e., trait anxiety) as measured by the Gray–Wilson Personality Questionnaire and Eysenck Personality Inventory in mostly females. However, their analysis aggregated alpha power across both eyes-closed and eyes-open states and six electrode locations (global alpha). Similarly, Knyazev et al. (2004) found the same association in males, but only used eyes-closed rest to measure the “resting state.” Moreover, in further analysis, they used activity during eyes-closed rest as a baseline to then evaluate the activity of the eyes-open resting states (before beep, after beep, or after word presentation). Given the drastic differences in eyes-open and closed resting-state oscillatory patterns, such analysis might yield confounded results or overlook a nuanced relationship between resting states (eyes open and eyes closed), frontal and parietal neural activity, and trait anxiety.
The greater predictive power of eyes-open resting-state oscillations for anxiety is consistent with the idea that eyes-open rest may be a more arousing state for individuals with high trait anxiety, whereas eyes-closed rest may be more relaxing. If FMT serves as a biomarker of behavioral inhibition, individuals with high trait anxiety may exhibit increased FMT during eyes-open rest as they prepare to inhibit unwanted stimuli or behaviors. In eyes-closed rest, when the inhibitory system is more at ease, FMT activity may decrease. In line with this reasoning, future research on resting state and individual differences should examine both eyes-open and eyes-closed resting conditions to capture these variations.
Overall, our findings suggest that anxiety-related lower frequency power and oscillations are detectable during more vigilant and potentially meaningful resting states, such as eyes open, especially in individuals with high-trait anxiety. Although both frontal midline alpha and frontal midline theta oscillations were found to be reliable predictors of behavioral inhibition, our results show that frontal midline theta is statistically more robust. This strengthens the case for frontal midline theta oscillations, which are strongly supported by theoretical frameworks, as a reliable neural correlate of behavioral inhibition (Gray & McNaughton, 2000; Shadli et al., 2021).
5 Limitations and Broader Downstream Implications
A potential limitation of the current study is that self-reports of trait anxiety can be influenced by biases and measurement errors (Kreitchmann et al., 2019; Wall & Lee, 2022). Future research could benefit from integrating alternative measures of anxiety, such as heart rate variability (Chalmers et al., 2014) or cortisol (Zorn et al., 2017) to mitigate these limitations. However, it is important to note that these measures have their own limitations and are primarily biomarkers of state rather than trait anxiety (Bennett, 2022). This underscores the need for more robust and reliable measures of both trait and state anxiety, and we hope that our findings contribute to advancing the development of such measures in the field.
It is also important to clarify that when we discuss frontal midline alpha and theta oscillations predicting anxiety, we refer to anxiety as a psychological trait, not a clinical disorder. Our study measured trait anxiety in healthy individuals, and we do not claim that these oscillatory activities are linked to any specific anxiety psychopathologies.
Given that trait anxiety is a known risk factor for anxiety disorders, future studies could examine the link between specific anxiety diagnoses and frontal midline alpha and theta oscillations. This could help develop objective biomarkers to improve diagnostic accuracy, which remains a challenge in psychiatry, where diagnoses still rely heavily on subjective assessments. However, exploring this in clinical settings goes beyond the scope of the current study.
6 Conclusion
In conclusion, this study demonstrates that both frontal midline theta and frontal midline alpha oscillations play a significant role in predicting trait anxiety and behavioral inhibition, particularly during eyes-open rest, where individuals with high trait anxiety may experience heightened arousal. Supported by stronger statistical results and robust theoretical evidence, the findings linking frontal midline theta to the behavioral inhibition system provide meaningful insights into the neural mechanisms underlying anxiety. The results also highlight the importance of distinguishing between eyes-open and eyes-closed resting states in future studies. Additionally, the use of the BOSC method to separate rhythmic from non-rhythmic activity provided clearer insights into the neural dynamics linked to anxiety, emphasizing the need to account for both signal types when studying traits such as anxiety. Together, these findings provide valuable insights into the underlying neural mechanisms of anxiety and behavioral inhibition and highlight key considerations for future research in this area.
Data and Code Availability
The datasets and code used in this study are available from the authors upon reasonable request.
Author Contributions
Tamari Shalamberidze: experimental design, data collection, data analysis, discussion, methodology, conceptualization, writing draft, revision. Kyle Nash: conceptualization, discussion, supervision, reviewing and editing. Jeremy B. Caplan: conceptualization, discussion, methodology, software, supervision, reviewing and editing.
Funding
Partly supported by the Social Sciences and Humanities Research Council in Canada, and the Natural Sciences and Engineering Research Council of Canada.
Declaration of Competing Interest
The authors have no conflicts of interests to declare.
References
Our primary focus is on the general population and the variance within it, encompassing the entire spectrum of anxiety levels. We intentionally did not exclude participants with high or low anxiety scores to ensure that we captured the full range of variability, which is crucial for our study objectives. To provide reassurance that individuals with self-reported anxiety diagnoses did not disproportionately affect our results or exhibit a distinct association pattern with BIS, we highlighted these participants on the Figure 3 scatterplots.
It is common for researchers using EGI’s dense-array caps (256 electrodes with minimal setup time) to analyze averages of neighboring electrode clusters to improve the signal-to-noise ratio. However, in our laboratory’s experience, with sufficient data, cluster averages produce results very similar to single electrodes. To prioritize simplicity and reliability in identifying markers of trait anxiety, we follow our laboratory’s standard of analyzing single electrodes, as this approach is not only effective but also provides an easier pathway for potential future diagnostic applications. Notably, signal quality was indirectly enhanced by the full 256 electrodes contributing to the ICA-driven artifact correction procedures.