<figure><img src="image_1.png" /></figure>

HAL
open science

# Effects of emotional valence and arousal on the voice perception network

Patricia E G Bestelmeyer, Sonja Kotz, Pascal C Belin

► To cite this version:

Patricia E G Bestelmeyer, Sonja Kotz, Pascal C Belin. Effects of emotional valence and arousal on the voice perception network. Social Cognitive and Affective Neuroscience, 2017, 12 (8), pp.1351-1358. 10.1093/scan/nsx059. hal-02008909

HAL Id: hal-02008909

https://amu.hal.science/hal-02008909

Submitted on 6 Feb 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Effects of emotional valence and arousal on the voice perception network

Patricia E. G. Bestelmeyer,¹ Sonja A. Kotz,²,³ and Pascal Belin⁴,⁵,⁶

¹School of Psychology, Bangor University, Bangor, Gwynedd, UK, ²Department of Neuropsychology and Psychopharmacology, Maastricht University, Maastricht, The Netherlands, ³Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany, ⁴Institut des Neurosciences de La Timone, UMR 7289, CNRS & Université Aix-Marseille, Marseille, France, ⁵Institute of Neuroscience and Psychology, University of Glasgow, Glasgow, UK, and ⁶International Laboratory for Brain, Music and Sound Research, University of Montréal & McGill University, Montréal, Canada

Correspondence should be addressed to Patricia E. G. Bestelmeyer, School of Psychology, Brigantia Building, Penrallt Road, Bangor University, Gwynedd LL57 2AS, UK. E-mail: p.bestelmeyer@bangor.ac.uk.

## Abstract

Several theories conceptualise emotions along two main dimensions: valence (a continuum from negative to positive) and arousal (a continuum that varies from low to high). These dimensions are typically treated as independent in many neuroimaging experiments, yet recent behavioural findings suggest that they are actually interdependent. This result has impact on neuroimaging design, analysis and theoretical development. We were interested in determining the extent of this interdependence both behaviourally and neuroanatomically, as well as teasing apart any activation that is specific to each dimension. While we found extensive overlap in activation for each dimension in traditional emotion areas (bilateral insulae, orbitofrontal cortex, amygdalae), we also found activation specific to each dimension with characteristic relationships between modulations of these dimensions and BOLD signal change. Increases in arousal ratings were related to increased activations predominantly in voice-sensitive cortices after variance explained by valence had been removed. In contrast, emotions of extreme valence were related to increased activations in bilateral voice-sensitive cortices, hippo-campi, anterior and midcingulum and medial orbito- and superior frontal regions after variance explained by arousal had been accounted for. Our results therefore do not support a complete segregation of brain structures underpinning the processing of affective dimensions.

**Key words:** arousal; emotion; fMRI; valence; voice

## Introduction

Emotions pervade every aspect of our life by influencing our physiology, perceptions and behaviours. Despite the prominence and effects of emotions on even basic sensory perception we know relatively little about how emotions are represented in the brain. The two most widely cited theories suggest that emotions are either characterised as discrete entities (i.e. basic emotion categories; Panksepp, 1982; Ekman, 1992) or as independent dimensions, arousal and valence (Russell, 1980). Arousal (or

intensity) is the level of autonomic activation that an event creates, and ranges from calm (or low) to excited (or high). Valence, on the other hand, is the level of pleasantness that an event generates and is defined along a continuum from negative to positive.

Basic emotion theories propose that humans are evolutionarily endowed with a limited set of emotions. Thus the behavioural and physiological expression of each emotion should be separate and arise from unique neural substrates (Ekman,
1992). While popular, this theory has had limited scientific support from animal, developmental, cross-cultural and physiological studies (Carroll et al., 1999). Neuroimaging studies using traditional univariate analyses have generally not supported the existence of consistent and specifically localised, distinct brain regions for the perception of different emotion categories (e.g. see Phan et al., 2002; Lindquist et al., 2012 for meta-analyses, and Hamann, 2012 for review). More recently, any studies confirming unique and consistent neural signatures using multivariate pattern analysis approaches have been fiercely criticised based on a misunderstanding of the statistics and of basic emotion theory (Clark-Polner et al., 2017, but see e.g., Kotz et al., 2013 for a cautious approach).

Dimensional theories, on the other hand, posit that any affective experience can be defined along at least two dimensions (Wundt, 1924). The essence of these theories is that emotions are multicomponent phenomena. Early support for a multi-dimensional organisation of the emotional space comes from psychophysiological studies. For example, startle reflex amplitudes are greatest for negatively valenced photos but decrease with positive emotional content in photos (Anders et al., 2004; Adank et al., 2012). Similarly, skin conductance amplitudes co-vary with increasing arousal ratings and electromyographic activity correlates with valence ratings (Lang et al., 1993). Additional support for a multicomponent view of emotion comes from pupillometry whereby pupil size is larger during emotionally negative and positive sounds than during the presentation of neutral sounds (Partala and Surakka, 2003).

While there is support for dimensional theories of emotion, the exact number and labelling of the dimensions as well as the relationship between these dimensions is less clear. Russell (1980) proposed that the dimensions of valence and arousal are *independent and linear*, i.e. how pleasant one is feeling gives no information about how calm or activated one feels. Intuitively, however, valence and arousal often go together (e.g. very negative events are also typically high in arousal), which is perhaps the reason why arousal has frequently been interpreted as reflecting the intensity of an emotion (Schlosberg, 1954; Schachter and Singer, 1962; Otten and Moskowitz, 2000). A relationship of this kind would imply valence and arousal being quadratically related and this weak but consistent relationship between the two dimensions has been confirmed in behavioural studies (for words: e.g. Bradley and Lang, 1994; Redondo et al., 2007; Vo et al., 2009; Kanske and Kotz, 2010 and pictures: e.g. Cuthbert et al., 1996; Bradley and Lang, 1999; see also Kuppens et al., 2013). Due to this relationship it is difficult to manipulate the two dimensions independently and find the neural underpinnings of each dimension separately.

Several neuroimaging studies have attempted to describe the independent neural networks underlying valence and arousal using very different designs, analyses, tasks and sensory stimuli and often without awareness of the findings in the behavioural literature. Thus many neuroimaging researchers have assumed a linear and orthogonal relationship between the two dimensions and designed their studies accordingly. While most studies report that arousal and valence are coded by distinct subsystems (Anderson et al., 2003; Anders et al., 2004; Dolcos et al., 2004; Lewis et al., 2007; Nielen et al., 2009; Colibazzi et al., 2010; Sieger et al., 2015) there is little consensus in the literature, which neural systems serve these functions. Nevertheless, lesion (Öngür and Price, 2000) and neuroimaging studies (Zuckerman et al., 1990; Viinikainen et al., 2010) tend to agree on a central role of the medial or orbital prefrontal cortex in the processing of emotional valence with both positive and

negative valence coded in overlapping neural regions (Lindquist et al., 2012). Reports frequently, but not consistently, cite the amygdala as the region that underpins the processing of arousal (Anderson et al., 2003; Small et al., 2003; Lewis et al., 2007). In addition to difficulties in replicating certain regions, another recent complication is that BOLD signal may not actually vary linearly with linear modulations of the valence spectrum (Viinikainen et al., 2010, 2012).

We were interested in localising the overlap between valence and arousal as well as examining the effects specific to each affective dimension on neural activity when processing non-linguistic, vocal expressions. To this end we created vocal morphs between neutral and various emotional expressions to broadly sample the valence and arousal space. We then computed two parametric models whereby each of the regressors of interest (i.e. ratings of valence and arousal) was orthogonalised with respect to the other in the design matrix. These two models allowed us to examine any overlap between valence and arousal activations as well as any areas specific to each dimension. We predicted quadratic relationships between valence ratings and fMRI signal in bilateral orbitofrontal cortex after variance accounted for by arousal ratings had been removed. In contrast, we expected increased subcortical activation with increased arousal ratings. Given previous behavioural studies, we predicted some overlap in the neural networks dealing with the processing of arousal and valence.

# Materials and methods

## Participants

Twenty-six healthy volunteers from the under- and postgraduate community of the University of Glasgow took part (12 females, mean age = 23.73, range = 18–36 years, s.d. = 5.61). All participants reported normal hearing and were reimbursed £12 for their time. Informed consent was obtained from all individuals and the study protocol was approved by the local ethics committee.

## Experimental stimuli

Recordings were taken from the Montreal Affective Voices (Belin et al., 2008) in which actors were instructed to produce emotional interjections using the vowel/a/. We used neutral expressions as well as expressions of fear, anger and pleasure from two identities (one female) to avoid identity-specific effects of emotion expression. These two identities had the highest recognition rates for all four expressions. Anger-fear, pleasure-fear, pleasure-anger continua as well as continua between each emotional expression and a neutral/a/were created separately for each identity, in nine steps that corresponded to –10/110%, 5/95%, 20/80%, 35/65%, 50/50%, 65/35%, 80/20%, 95/5% and 110/–10%, e.g. pleasure/fear. –10/110% and 110/–10% were caricatured versions of the morphed expressions. We therefore had six continua and a total of 54 morphs per male and female identity. The duration of the vocalisations within each continuum was kept constant and ranged between 0.6 and 1.0s across different continua (the varied stimulus durations are incorporated in our models). We used STRAIGHT (Kawahara and Matsui, 2003); see also http://www.wakayama-u.ac.jp/~kawahara/STRAIGHTadv/index_e.html for more information) for stimulus generation and Psychtoolbox-3 (Brainard, 1997; Kleiner et al., 2007) for stimulus presentation and response recording for the fMRI and validation study. Toolboxes and experimental scripts were run on MatlabR2007b (Mathworks, Inc).
STRAIGHT performs an instantaneous pitch-adaptive spectral smoothing of each stimulus for separation of contributions to the voice signal arising from the glottal source (including f0) vs supralaryngeal filtering (distribution of spectral peaks, including the first formant frequency, F1). Voice stimuli were decomposed by STRAIGHT into five parameters: f0 (the perceived pitch of the voice), frequency, duration, spectrotemporal density and aperiodicity. Each parameter can be manipulated independently. For each voice we manually identified one time landmark with three frequency landmarks (corresponding to the first three formants) at the onset of phonation and the same number of landmarks at the offset of phonation. Morphed stimuli were then generated by re-synthesis based on the interpolation (linear for time; logarithmic for F0, frequency, and amplitude) of these time-frequency landmark templates (see also Schweinberger et al. (2014) for a recent discussion of the voice morphing technique). Stimuli were normalised in energy (root mean square) before and after morphing. Acoustic analyses to illustrate the validity of the morphing technique with emotional sounds are published elsewhere (Bestelmeyer et al., 2010).

Sixty-five undergraduates (45 female) at Bangor University rated all morphs on valence and arousal in return for course credit on an analogue scale from 0 to 1 (from 'extremely negative'/'extremely low in arousal' to 'extremely positive'/'extremely high in arousal'). Stimuli were presented randomly (but blocked for identity). Half of our participants were asked to rate the morphs on valence first and then arousal, while the other half conducted the blocks in reverse order. For illustration purposes we 'min-max' normalised and averaged the data across the two identities for each morph type and rating task (Figure 1). For the parametric modulation analyses we computed 'min-max' normalised valence and arousal regressors separately for each morph identity.

## fMRI paradigm

The experimental run consisted of a fast event-related design and was always acquired first. It involved five phases during which we presented 54 morphs of the female identity and five phases during which we presented 54 morphs of the male identity. During each of these phases a vocal morph was presented every 2 s at the beginning of each TR. The order of the phases and the order of morphs within the phases was randomised using Matlab's *rand* function. Each phase of 54 morphs was separated by 20 TRs of silences (40 s) before beginning the next phase (see supplementary online figure for an illustration of the paradigm). Participants were asked to perform a three-

alternative forced choice task, in which each affective morph had to be categorised as either anger, fear or pleasure by means of an MRI compatible response box (FORP; Current Designs, Inc.). Participants were asked to react as fast as possible and to keep their eyes closed for the duration of all scans.

The experimental run was followed by a voice localiser scan. This block design uses vocal and environmental sounds to reliably identify the temporal voice areas (TVAs) using the vocal vs non-vocal contrast (for details see Belin et al., 2000). All sounds were presented binaurally using the electrostatic NNL headphone system (NordicNeuroLab, Inc.) at an intensity of 80 dB SPL(C). Lastly, we collected a detailed anatomical scan.

## Image acquisition

All scans were acquired in a 3.0 Tesla Siemens Tim Trio scanner using a 32-channel head coil. A whole brain T1-weighted anatomical scan was performed using fast gradient echo sequence known as T1 'Magnetization Prepared Rapid Gradient Echo' (MPRAGE) consisting of 192 axial slices of 1 mm thickness with an inplane resolution of 1×1×1 (Field of view (FOV) = 256) and a matrix of 256×256 performed at the end of the experimental session. Two T2*-weighted functional scans were acquired using an interleaved ascending sequence consisting of 32 slices of 3 mm thickness (0.3 mm gap) with an inplane resolution of 3×3×3 (FOV = 210), an acquisition matrix of 70×70 and a flip angle of 77°. The experimental scan (TR = 2 s, TE = 30 ms) consisted of 740 volumes while the voice localiser scan consisted of 310 volumes (otherwise same scanning parameters as the experimental scan). All MRI data were analysed using SPM12 (Wellcome Department of Cognitive Neurology, University College London, http://www.fil.ion.ucl.ac.uk/spm/). Pre-processing of the data consisted of AC-PC alignment of the anatomical images (and application of the re-orientation matrix to all functional images acquired in the same session). Functional scans were corrected for head motion (trilinear interpolation) by aligning all scans to the first scan of the last run and a mean image was created. The anatomical scan was co-registered to the mean image. Functional (3 mm isotropic voxels) and anatomical (1 mm isotropic voxels) data were transformed to Montreal Neurological Institute (MNI) space after segmentation of the anatomical scan. Normalised data were spatially smoothed by applying a Gaussian kernel of 8 mm full width at half maximum (FWHM).

## Image analysis

We used parametric modulations to analyse the data (Büchel et al., 1998) with two parametric modulators loading on the

<figure><img src="image_3.png" /><figcaption>Fig. 1 Behavioural results. Behavioural ratings of (A) arousal and (B) valence of the stimulus set. (C) Scatterplot illustrating the quadratic relationship between valence and arousal ratings.</figcaption></figure>
sound onset regressor. The reported results are based on whole-brain analyses rather than a region of interest approach. Based on previous work we computed linear and quadratic expansions of each regressor (Lewis et al., 2007; Viinikainen et al., 2010, 2012). Note that in SPM each parametric modulator is orthogonalised with respect to the previous ones to only account for variance not explained by previous regressors. We computed two models. For the first model we entered arousal as first parametric modulator and valence as the second. This allowed us to examine the linear and quadratic relationships between BOLD signal and valence after variance explained by arousal had been removed. The second model consisted of the reverse order of parametric modulators and allowed us to examine linear and quadratic relationships between BOLD and arousal after variance explained by valence had been removed. Since our two regressors of interest are not naturally orthogonal, both regressors will share variance. We therefore also examined the activations of the first parametric modulators of each model and illustrate the overlap in Figure 2.

In addition to the two parametric modulators, both models contained the six movement parameters as the last six regressors. At the second level of analysis we computed an F-test across linear and quadratic contrasts in a full factorial design for arousal (in the valence then arousal model; model 1) and for valence (in the arousal then valence model; model 2).

To illustrate the parameter estimates in Figure 3, we used SPM12's built-in functions to extract the beta estimates at the peak maxima within a sphere of 6 mm. Behavioural responses did not average to discrete steps but rather, resulted in a range of numbers between 0 and 1. We therefore had to group the averaged behavioural responses into five bins of 20%. Results are illustrated on an average anatomical scan using MRIcron (Rorden et al., 2007). Anatomical cluster location was assessed with xjview (8.1; http://www.alivelearn.net/xjview) and cross checked with Duvernoy's brain atlas (Duvernoy, 1999) to ensure accuracy. Activations are illustrated and reported in text and table at a threshold of P < 0.001; extent threshold (k) > 50 with FWE-correction at the cluster level.

# Results

## Behavioural

Results of the behavioural data from independent raters is summarised in Figure 1. Valence and arousal ratings were averaged separately as a function of morph step and continuum type. Figure 1A illustrates the averaged arousal ratings for each morph step and continuum type. Figure 1B shows the averaged valence ratings. Using non-linear regression we found a significant quadratic relationship between valence and arousal ratings ($R^2 = 0.601$; $F(2, 107) = 79.193$; P < 0.0001; Figure 1C).

<figure><img src="image_4.png" /><figcaption>Fig. 2 Overlapping activation maps for significant correlations between BOLD signal and arousal (blue) or valence (red) ratings (but without first removing variance explained by valence and arousal, respectively).</figcaption></figure>

## Neuroimaging

We used parametric modulation analyses to investigate areas that respond to a systematic change in valence or arousal with and without removing variance explained by the other dimension. We found significant positive correlations between arousal ratings and BOLD signal and between valence ratings and BOLD signal without first removing variance explained by the other dimension (see Table 1 for significant clusters for each dimension and BOLD signal). We then examined the shared variance or overlap (pink) between arousal (blue) and valence (red; Figure 2). Substantial overlap was apparent for bilateral insulae, superior temporal gyri (STG), bilateral amygdalae, caudate, precuneus, medial orbito- and superior frontal cortex as well as inferior and mid frontal cortex.

We were ultimately interested in which regions are specific to each emotional dimension. We therefore examined correlations between arousal after variance explained by valence had been accounted for (Figure 3A and B; model 1) and vice versa (Figure 3C and D; model 2). The first model, Figure 3A and B, revealed a correlation between BOLD and arousal ratings in bilateral superior temporal cortex with the cluster extending to bilateral temporal pole (right peak cluster: $F(2, 50) = 50.07$, $k = 1019$, MNI coordinates: $63 -7 -5$; left peak cluster: $F(2, 50) = 38.11$, $k = 803$, MNI coordinates: $-54 -10 -2$), right mid-temporal cortex with the cluster extending into mid and superior occipital cortex (peak cluster: $F(2, 50) = 27.27$, $k = 423$, MNI coordinates: $48 -731$), left superior and midoccipital cortex (peak cluster in left: $F(2, 50) = 13.29$, $k = 146$, MNI coordinates: $-9 -977$) and bilateral caudate (peak cluster in left: $F(2, 50) = 16.10$, $k = 142$, MNI coordinates: $-12 2610$). The line graphs indicate that BOLD signal increased with sounds that were rated as higher in arousal. Figure 3B reveals the overlap between this model and the TVA.

The second model, Figure 3C and D, revealed a correlation between BOLD and valence ratings in right superior temporal cortex (peak cluster: $F(2, 50) = 55.00$, $k = 2242$, MNI coordinates: $66 -10 -5$), left mid and superior temporal cortex (peak cluster in left: $F(2, 50) = 43.30$, $k = 2246$, MNI coordinates: $-66 -284$), bilateral hippocampi and parahippocampi (left peak cluster: $F(2, 50) = 28.06$, $k = 192$, MNI coordinates: $-33 -34 -11$; right peak cluster: $F(2, 50) = 18.72$, $k = 120$, MNI coordinates: $33 -34 -14$), bilateral midcingulum (peak cluster in right: $F(2, 50) = 25.79$, $k = 778$, MNI coordinates: $12 -2240$), bilateral superior midfrontal cortex (peak cluster in left: $F(2, 50) = 22.01$, $k = 309$, MNI coordinates: $-6 2346$), bilateral anterior cingulum (peak cluster in left: $F(2, 50) = 21.12$, $k = 453$, MNI coordinates: $029 -5$), bilateral cerebellum (left: $F(2, 50) = 15.69$, $k = 49$, MNI coordinates: $-9 -79 -29$; right: $F(2, 50) = 17.37$, $k = 71$, MNI coordinates: $9 -76 -23$), left midoccipital ($F(2, 50) = 15.70$, $k = 103$, MNI coordinates: $-42 -7634$), right midtemporal and midoccipital cortex (peak cluster: $F(2, 50) = 14.88$, $k = 209$, MNI coordinates: $39 -6713$) and left inferior parietal cortex ($F(2, 50) = 17.60$, $k = 104$, MNI coordinates: $-33 -4940$). The line graphs indicate a positive quadratic relationship between BOLD and sounds from the extreme ends of the valence continuum. Figure 3D reveals the overlap between this model and the TVA.

## Discussion

The aim of the current study was to localise and quantify the overlap between the dimensions of affect as well as examine any effects specific to each affective dimension on neural activity during the processing of vocal expressions. Whole-brain
<figure><img src="image_5.png" /><figcaption>Fig. 3 Activation maps of significant correlations between BOLD signal and (A) arousal after variance explained by valence has been removed. Significant relationships are evident in bilateral superior temporal gyri and right mid-occipital cortex. Overlap between this activation (blue) and that of the voice localiser (yellow) is shown in (B). Activation maps of significant correlations between BOLD signal and (C) valence after variance explained by arousal has been removed. Significant relationships are evident in bilateral superior temporal gyri, inferior frontal gyri, hippocampi, medial orbitofrontal, superior frontal gyri and midcingulum. Overlap between this activation (red) and that of the voice localiser (yellow) is shown in (D). Activations in (A) and (C) are illustrated on a T1-weighted average structural template. The line graphs illustrate BOLD signal change in response to the behavioural ratings (from 0 to 1; low to high arousal and negative to positive valence, respectively, in five bins of 0.2).</figcaption></figure>

analyses revealed significant overlap between activations for arousal and valence in traditional emotion areas (bilateral insulae, orbitofrontal cortex and amygdalae) as well as modality-specific, auditory areas. To investigate the specific relationship between BOLD signal and each affective dimension, we computed two parametric models, in which we orthogonalised each dimension with respect to the other. Increased arousal ratings correlated with increased BOLD signal in bilateral temporal voice areas (superior temporal sulcus/gyrus; STS/STG) and mid-occipital gyri after variance explained by valence had been accounted for (Figure 3A). In contrast, BOLD signal was increased with extremely positive or negative valence ratings compared to more neutral valence ratings in bilateral prefrontal cortex, mid- and anterior cingulum, hippocampi and parahippocampi but also bilateral temporal voice areas (STS/STG; Figure 3C). Our results suggest that the representations of valence and arousal are in part underpinned by sensory specific brain areas, the temporal voice areas, but also by brain areas (e.g. orbitofrontal cortex), which have been reported by various other studies of different modalities (Anderson et al., 2003; Chikazoe et al., 2014; Skerry and Saxe, 2015). In contrast to previous literature, our findings are not in line with the notion that valence and arousal are independent or orthogonal dimensions of emotion as previously thought and are not subserved by entirely distinct brain regions (Anderson et al., 2003; Small et al., 2003; Anders et al., 2004; Dolcos et al., 2004; Lewis et al., 2007; Nielen et al., 2009; Colibazzi et al., 2010).

Pernet et al. (2015) recently reported in a sample of over 200 participants that the temporal voice areas (TVA) consist of three bilateral voice-sensitive clusters. These clusters are located along bilateral posterior, mid and anterior STG, respectively.

While the activations to orthogonalised acoustic contributions to valence and arousal overlapped with much of the TVA (see also Ethofer et al., 2012), it is noteworthy that acoustic contributions to arousal were processed more anteriorly than the acoustic contributions to valence, which were processed more posteriorly (as indicated by each respective peak cluster). Our study is in line with previous research highlighting that the TVA is in charge of processing the underlying acoustic cues to emotion (Ethofer et al., 2006; Wiethoff et al., 2008; Frühholz et al., 2012; Bestelmeyer et al., 2014) and additionally alludes to an interesting dissociation in functionality between anterior and more posterior portions of the TVA.

Traditionally, fMRI studies have compared BOLD signal in response to high vs low valenced stimuli and frequently reported separate brain regions for each contrast (Dolcos et al., 2004; Nielen et al., 2009) in predefined regions of interest. These findings have led to speculations that valence may consist of two independent continua, one for neutral to negative and one for neutral to positive valence. With newer, parametric statistical methods it has become apparent that the valence dimension is unlikely to consist of two independent continua. There is now more and more support for the notion that brain regions respond to modulations in valence not linearly but quadratically (Lewis et al., 2009; Viinikainen et al., 2012), i.e. brain activity increases from neutral to most extreme emotion irrespective of whether the emotion is positive or negative. Our results are in line with this notion. Arousal, on the other hand, is still thought of as a linear continuum with activations increasing with increased arousal ratings. Our data show that activations in bilateral STS/STG increase largely linearly up to a plateau.
<table><caption>Table 1. Areas revealing significant modulations with arousal and valence (before removing variance by the other dimension)</caption><thead><tr><th rowspan="2">Modulator</th><th rowspan="2">Anatomical definition</th><th colspan="3">Peak voxel coordinates</th><th rowspan="2">F-value</th><th rowspan="2">Cluster size</th></tr><tr><th>x</th><th>y</th><th>Z</th></tr></thead><tbody><tr><td rowspan="15">AROUSAL</td><td>Left inferior frontal operculum</td><td>-45</td><td>20</td><td>31</td><td>30.41</td><td rowspan="3">940</td></tr><tr><td>Left precentral</td><td>-42</td><td>8</td><td>34</td><td>27.91</td></tr><tr><td>Left inferior frontal (pars triangularis)</td><td>-39</td><td>32</td><td>19</td><td>27.40</td></tr><tr><td>Right inferior frontal (pars triangularis)</td><td>45</td><td>23</td><td>28</td><td>28.65</td><td rowspan="3">556</td></tr><tr><td>Right inferior frontal operculum</td><td>45</td><td>8</td><td>28</td><td>17.21</td></tr><tr><td>Right midfrontal</td><td>42</td><td>38</td><td>16</td><td>13.51</td></tr><tr><td>Left precuneus</td><td>-6</td><td>-64</td><td>46</td><td>23.88</td><td rowspan="3">339</td></tr><tr><td>Left inferior parietal</td><td>-30</td><td>-49</td><td>40</td><td>22.07</td></tr><tr><td>Left inferior parietal</td><td>-36</td><td>-55</td><td>46</td><td>21.31</td></tr><tr><td>Right amygdala</td><td>24</td><td>-4</td><td>-17</td><td>22.55</td><td rowspan="3">623</td></tr><tr><td>Right superior temporal</td><td>48</td><td>-1</td><td>-11</td><td>20.71</td></tr><tr><td>Right superior temporal</td><td>51</td><td>-10</td><td>1</td><td>15.46</td></tr><tr><td>Right insula</td><td>30</td><td>26</td><td>-2</td><td>20.56</td><td>120</td></tr><tr><td>Left insula</td><td>-30</td><td>20</td><td>1</td><td>20.44</td><td>161</td></tr><tr><td>Left inferior frontal (pars triangularis)</td><td>-48</td><td>17</td><td>7</td><td>11.36</td><td rowspan="3">395</td></tr><tr><td>Left superior temporal</td><td>-48</td><td>-16</td><td>-2</td><td>19.76</td></tr><tr><td>Left superior temporal</td><td>-42</td><td>-13</td><td>-8</td><td>17.21</td></tr><tr><td>Left superior temporal</td><td>-54</td><td>-7</td><td>4</td><td>15.34</td><td rowspan="3">1082</td></tr><tr><td>Right caudate</td><td>6</td><td>11</td><td>-8</td><td>19.39</td></tr><tr><td>Left medial orbitofrontal</td><td>0</td><td>38</td><td>-17</td><td>17.75</td></tr><tr><td>Left medial orbitofrontal</td><td>-9</td><td>53</td><td>-11</td><td>17.44</td><td rowspan="3">76</td></tr><tr><td>Right inferior parietal</td><td>36</td><td>-49</td><td>40</td><td>13.62</td></tr><tr><td>Right angular</td><td>33</td><td>-58</td><td>43</td><td>12.37</td></tr><tr><td rowspan="15">VALENCE</td><td>Left supplementary motor area/medial frontal</td><td>-6</td><td>20</td><td>46</td><td>51.43</td><td rowspan="3">1775</td></tr><tr><td>Right supplementary motor area</td><td>6</td><td>20</td><td>46</td><td>36.45</td></tr><tr><td>Right insula</td><td>30</td><td>23</td><td>1</td><td>28.46</td></tr><tr><td>Left inferior frontal (pars opercularis)</td><td>-42</td><td>20</td><td>31</td><td>35.37</td><td rowspan="3">1533</td></tr><tr><td>Left precentral</td><td>-39</td><td>8</td><td>31</td><td>30.34</td></tr><tr><td>Left insula</td><td>-30</td><td>20</td><td>1</td><td>30.32</td></tr><tr><td>Left medial frontal</td><td>-6</td><td>50</td><td>16</td><td>29.50</td><td rowspan="3">2322</td></tr><tr><td>Left anterior cingulum</td><td>0</td><td>29</td><td>-5</td><td>28.41</td></tr><tr><td>Left medial orbitofrontal</td><td>-6</td><td>35</td><td>-14</td><td>25.33</td></tr><tr><td>Left midtemporal</td><td>-45</td><td>-64</td><td>19</td><td>26.50</td><td rowspan="3">237</td></tr><tr><td>Left midtemporal</td><td>-57</td><td>-61</td><td>16</td><td>14.68</td></tr><tr><td>Left midtemporal</td><td>-30</td><td>-52</td><td>19</td><td>10.07</td></tr><tr><td>Left inferior parietal</td><td>-33</td><td>-49</td><td>40</td><td>26.19</td><td rowspan="3">272</td></tr><tr><td>Left superior temporal</td><td>-51</td><td>-34</td><td>16</td><td>24.21</td></tr><tr><td>Left insula</td><td>-39</td><td>-13</td><td>-2</td><td>18.91</td></tr><tr><td>Left midtemporal</td><td>-60</td><td>-7</td><td>-14</td><td>18.77</td><td rowspan="3">68</td></tr><tr><td>Left precuneus</td><td>-9</td><td>-64</td><td>46</td><td>22.00</td></tr><tr><td>Right precuneus</td><td>9</td><td>-64</td><td>46</td><td>14.95</td></tr><tr><td>Right inferior parietal</td><td>36</td><td>-46</td><td>40</td><td>17.41</td><td rowspan="3">111</td></tr><tr><td>Right angular</td><td>33</td><td>-58</td><td>43</td><td>11.64</td></tr><tr><td>Left posterior cingulate</td><td>-6</td><td>-49</td><td>28</td><td>16.70</td></tr><tr><td>Left midcingulum</td><td>-12</td><td>-31</td><td>40</td><td>15.22</td><td rowspan="3">302</td></tr><tr><td>Left midcingulum</td><td>-15</td><td>-25</td><td>46</td><td>14.59</td></tr><tr><td>Right midcingulum</td><td>15</td><td>-22</td><td>34</td><td>12.25</td></tr><tr><td>Right midcingulum</td><td>15</td><td>-34</td><td>49</td><td>9.99</td><td>57</td></tr></tbody></table>

Using single-cell recordings Sieger et al. (2015) have recently shown that neurons within the same nucleus can independently code valence or arousal. In other words, overlap of activations to valence and arousal in functional neuroimaging studies does not exclude the possibility that neurons within this region are dedicated for the processing of one specific affective dimension. Thus, while the current study does not support the notion of distinct brain regions underpinning each affective dimension, it is still possible that our overlapping activations in these frontal areas contain specialised neurons for affective dimensions.

Bestelmeyer et al. (2014) have recently shown empirical support for Schirmer and Kotz’s (2006) notion of a multi-step process in the perception of auditory emotion. This fMRI study demonstrated that the temporal voice areas and amygdala likely deal with the acoustic features of the vocal emotion while the cognitive representation of emotion is underpinned by bilateral insulae, precuneus, superior frontal gyri and midcingulum. It is likely that auditory cues to valence and arousal are processed in bilateral temporal voice areas while cognitive representations of the positivity and intensity of the emotion are processed in dimension-specific, possibly supra-modal areas.
Due to differences in experimental designs, stimulus material and analyses in the past it remains unclear, which of these regions, if any, deal with affective dimensions at a supra-modal level.

Our study is the first to quantify the overlap in activation for each affective dimension as well as demonstrate activation specific to each dimension. Yet many open questions remain. Naturally, and as in our data set, arousal for negative expressions is typically higher than for positive ones (Ito et al., 1998; Ito and Cacioppo, 2005). It is unclear whether the response of the TVA would be modulated further with stimuli perfectly matched for arousal. Furthermore, while our study used a large number of carefully controlled stimuli the addition of null events and jittering of the stimulus onsets could have further optimised our design.

## Conclusion

Our results are not in line with the theory that valence and arousal are strictly independent dimensions of emotion (Russell, 1980). We found substantial correlations between behavioural ratings of valence and arousal as well as anatomical overlap of regions, which code modulations in valence and arousal, particularly in bilateral insulae, orbitofrontal cortex and amygdalae. We also noted significant correlations between brain activity and the behavioural ratings that were specific to each dimension. We found that increased arousal ratings correlated positively with brain activity predominantly in regions, which code acoustic properties of the sound, i.e. the temporal voice areas. In addition, we found that extreme valence ratings activated a more widespread network again involving the bilateral temporal voice areas but also bilateral hippocampi, anterior and midcingulum, medial orbito- and superior frontal regions.

## Supplementary data

Supplementary data are available at SCAN online.

## Funding

This work was supported by the Economic and Social Research Council/Medical Research Council grant (RES-060-25-0010). The authors declare no competing financial interests.

*Conflict of interest*. None declared.

## References

Adank, P., Davis, M.H., Hagoort, P. (2012). Neural dissociation in processing noise and accent in spoken language comprehension. *Neuropsychologia*, **50**, 77–84.

Anders, S., Lotze, M., Erb, M., Grodd, W., Birbaumer, N. (2004). Brain activity underlying emotional valence and arousal: a response-related fMRI study. *Human Brain Mapping*, **23**, 200–9.

Anderson, A.K., Christoff, K., Stappen, I., et al. (2003). Dissociated neural representations of intensity and valence in human olfaction. *Nature Neuroscience*, **6**, 196–202.

Belin, P., Fillion-Bilodeau, S., Gosselin, F. (2008). The Montreal affective voices: a validated set of nonverbal affect bursts for research on auditory affective processing. *Behavior Research Methods*, **40**, 531–9.

Belin, P., Zatorre, R.J., Lafaille, P., Ahad, P., Pike, B. (2000). Voice-selective areas in human auditory cortex. *Nature*, **403**, 309–12.

Bestelmeyer, P.E.G., Maurage, P., Rouger, J., Latinus, M., Belin, P. (2014). Adaptation to vocal expressions reveals multistep perception of auditory emotion. *Journal of Neuroscience*, **34**, 8098–105.

Bestelmeyer, P.E.G., Rouger, J., DeBruine, L.M., Belin, P. (2010). Auditory adaptation in vocal affect perception. *Cognition*, **117**, 217–23.

Bradley, M.M., Lang, P.J. (1994). *Affective norms for english words (ANEW): stimuli, instruction manual and affective ratings*. Technical report C-1. Gainsville, FL: The Center for Research in Psychophysiology, University of Florida.

Bradley, M.M., Lang, P.J. (1999). Fearfulness and affective evaluations of pictures. *Motivation and Emotion*, **23**, 1–13.

Brainard, D.H. (1997). The psychophysics toolbox. *Spatial Vision*, **10**, 433–6.

Büchel, C., Holmes, A.P., Rees, G., Friston, K.J. (1998). Characterizing stimulus-response functions using nonlinear regressors in parametric fMRI experiments. *Neuroimage*, **8**, 140–8.

Carroll, J.M., Yik, M.S.M., Russell, J.A., Feldman Barrett, L. (1999). On the psychometric principles of affect. *Review of General Psychology*, **3**, 14–22.

Chikazoe, J., Lee, D.H., Kriegeskorte, N., Anderson, A.K. (2014). Population coding of affect across stimuli, modalities and individuals. *Nature Neuroscience*, **17**, 1114–22.

Clark-Polner, E., Johnson, T.D., Feldman Barrett, L. (2017). Multivoxel pattern analysis does not provide evidence to support the existence of basic emotions. *Cerebral Cortex*, **27**, 1944–8.

Colibazzi, T., Posner, J., Wang, Z., et al. (2010). Neural systems subserving valence and arousal during the experience of induced emotions. *Emotion*, **10**, 377–89.

Cuthbert, B.N., Bradley, M.M., Lang, P.J. (1996). Probing picture perception: activation and emotion. *Psychophysiology*, **33**, 103–11.

Dolcos, F., LaBar, K.S., Cabeza, R. (2004). Dissociable effects of arousal and valence on prefrontal activity indexing emotional evaluation and subsequent memory: an event-related fMRI study. *Neuroimage*, **23**, 64–74.

Duvernoy, H.M. (1999). *The Human Brain: Surface, Blood Supply, and Three-Dimensional Sectional Anatomy*. Wien: Springer Verlag.

Ekman, P. (1992). An argument for basic emotions. *Cognition and Emotion*, **6**, 169–200.

Ethofer, T., Anders, S., Wiethoff, S., et al. (2006). Effects of prosodic emotional intensity on activation of associative auditory cortex. *Neuroreport*, **17**, 249–53.

Ethofer, T., Bretscher, J., Gschwind, M., Kreifelts, B., Wildgruber, D., Vuilleumier, P. (2012). Emotional voice areas: anatomic location, functional properties, and structural connections revealed by combined fMRI/DTI. *Cerebral Cortex*, **22**, 191–200.

Frühholz, S., Ceravolo, L., Grandjean, D. (2012). Specific brain networks during explicit and implicit decoding of emotional prosody. *Cerebral Cortex*, **22**, 1107–17.

Hamann, S. (2012). Mapping discrete and dimensional emotions onto the brain: controversies and consensus. *Trends in Cognitive Science*, **16**, 458–66.

Ito, T.A., Cacioppo, J.T. (2005). Variations on a human universal: individual differences in positivity offset and negativity bias. *Cognition & Emotion*, **19**, 1–26.

Ito, T.A., Cacioppo, J.T., Lang, P.J. (1998). Eliciting affect using the International Affective Picture System: trajectories through evaluative space. *Personality and Social Psychology Bulletin*, **24**, 855–79.

Kanske, P., Kotz, S.A. (2010). Leipzig affective norms for German: a reliability study. *Behavior Research Methods*, **42**, 987–91.
Kawahara, H., Matsui, H. (2003). Auditory morphing based on an elastic perceptual distance metric in an interference-free time-frequency representation. In: 2003 Ieee International Conference on Acoustics, Speech, and Signal Processing, Vol I, Proceedings, 256-9.

Kleiner, M., Brainard, D., Pelli, D. (2007). What's new in Psychtoolbox-3? 36 ECVP Abstract Supplement 14.

Kotz, S.A., Kalberlah, C., Bahlmann, J., Friederici, A.D., Haynes, J.D. (2013). Predicting vocal emotion expressions from the human brain. *Human Brain Mapping*, **34**, 1971–81.

Kuppens, P., Tuerlinckx, F., Russell, J.A., Barrett, L.F. (2013). The relation between valence and arousal in subjective experience. *Psychological Bulletin*, **139**, 917–40.

Lang, P.J., Greenwald, M.K., Bradley, M.M., Hamm, A.O. (1993). Looking at pictures – affective, facial, visceral and behavioral reactions. *Psychophysiology*, **30**, 261–73.

Lewis, J.W., Talkington, W.J., Walker, N.A., et al. (2009). Human cortical organization for processing vocalizations indicates representation of harmonic structure as a signal attribute. *Journal of Neuroscience*, **29**, 2283–96.

Lewis, P.A., Gritchley, H.D., Rotshtein, P., Dolan, R.J. (2007). Neural correlates of processing valence and arousal in affective words. *Cerebral Cortex*, **17**, 742–8.

Lindquist, K.A., Wager, T.D., Kober, H., Bliss-Moreau, E., Barrett, L.F. (2012). The brain basis of emotion: a meta-analytic review. *Behavioral and Brain Sciences*, **35**, 121–43.

Nielen, M.M.A., Heslenfeld, D.J., Heinen, K., et al. (2009). Distinct brain systems underlie the processing of valence and arousal of affective pictures. *Brain and Cognition*, **71**, 387–96.

Öngür, D., Price, J.L. (2000). The organization of networks within the orbital and medial prefrontal cortex of rats, monkeys and humans. *Cerebral Cortex*, **10**, 206–19.

Otten, S., Moskowitz, G.B. (2000). Evidence for implicit evaluative in-group bias: affect-biased spontaneous trait inference in a minimal group paradigm. *Journal of Experimental Social Psychology*, **36**, 77–89.

Panksepp, J. (1982). Toward a general psychobiological theory of emotions. *Behavioral and Brain Sciences*, **5**, 407–22.

Partala, T., Surakka, V. (2003). Pupil size variation as an indication of affective processing. *International Journal of Human-Computer Studies*, **59**, 185–98.

Pernet, C.R., McAleer, P., Latinus, M., et al. (2015). The human voice areas: spatial organization and inter-individual variability in temporal and extra-temporal cortices. *Neuroimage*, **119**, 164–74.

Phan, K.L., Wager, T., Taylor, S.F., Liberzon, I. (2002). Functional neuroanatomy of emotion: a meta-analysis of emotion activation studies in PET and fMRI. *Neuroimage*, **16**, 331–48.

Redondo, J., Fraga, I., Padron, I., Comesna, M. (2007). The Spanish adaptation of ANEW (Affective Norms for English words). *Behavior Research Methods*, **39**, 600–5.

Rorden, C., Karnath, H.-O., Bonilha, L. (2007). Improving lesion-symptom mapping. *Journal of Cognitive Neuroscience*, **19**, 1081–8.

Russell, J.A. (1980). A circumplex model of affect. *Journal of Personality and Social Psychology*, **39**, 1161–78.

Schachter, S., Singer, J.E. (1962). Cognitive, social, and physiological determinants of emotional state. *Psychological Review*, **69**, 379–99.

Schirmer, A., Kotz, S.A. (2006). Beyond the right hemisphere: brain mechanisms mediating vocal emotional processing. *Trends in Cognitive Science*, **10**, 24–30.

Schlosberg, H. (1954). Three dimensions of emotion. *Psychological Review*, **61**, 81–8.

Schweinberger, S.R., Kawahara, H., Simpson, A.P., Skuk, V.G., Zaeske, R. (2014). Speaker perception. *Wiley Interdisciplinary Reviews-Cognitive Science*, **5**, 15–25.

Sieger, T., Serranova, T., Ruzicka, F., et al. (2015). Distinct populations of neurons respond to emotional valence and arousal in the human subthalamic nucleus. *Proceedings of the National Academy of Sciences USA*, **112**, 3116–21.

Skerry, A.E., Saxe, R. (2015). Neural representations of emotion are organized around abstract event features. *Current Biology*, **25**, 1945–54.

Small, D.M., Gregory, M.D., Mak, Y.E., Gitelman, D., Mesulam, M.M., Parrish, T. (2003). Dissociation of neural representation of intensity and affective valuation in human gustation. *Neuron*, **39**, 701–11.

Viinikainen, M., Jaaskelainen, I.R., Alexandrov, Y., Balk, M.H., Autti, T., Sams, M. (2010). Nonlinear relationship between emotional valence and brain activity: evidence of separate negative and positive valence dimensions. *Human Brain Mapping*, **31**, 1030–40.

Viinikainen, M., Katsyri, J., Sams, M. (2012). Representation of perceived sound valence in the human brain. *Human Brain Mapping*, **33**, 2295–305.

Vo, M.L.H., Conrad, M., Kuchinke, L., Urton, K., Hofmann, M.J., Jacobs, A.M. (2009). The Berlin affective word list reloaded (BAWL-R). *Behavior Research Methods*, **41**, 534–8.

Wiethoff, S., Wildgruber, D., Kreifelts, B., et al. (2008). Cerebral processing of emotional prosody – influence of acoustic parameters and arousal. *Neuroimage*, **39**, 885–93.

Wundt, W. (1924). *An Introduction to Psychology*. London: Allen & Unwin.

Zuckerman, M., Hodgins, H., Miyake, K. (1990). The vocal attractiveness stereotype – replication and elaboration. *Journal of Nonverbal Behavior*, **14**, 97–112.