The noradrenergic locus coeruleus (LC) is a controller of brain and behavioral states. Activating LC neurons en masse by electrical or optogenetic stimulation promotes a stereotypical "activated" cortical state of high-frequency oscillations. However, it has been recently reported that spontaneous activity of LC cell pairs has sparse yet structured time-averaged cross-correlations, which is unlike the highly synchronous neuronal activity evoked by stimulation. Therefore, LC population activity could consist of distinct multicell ensembles each with unique temporal evolution of activity. We used nonnegative matrix factorization (NMF) to analyze large populations of simultaneously recorded LC single units in the rat LC. NMF identified ensembles of spontaneously coactive LC neurons and their activation time courses. Since LC neurons selectively project to specific forebrain regions, we hypothesized that distinct ensembles activate during different cortical states. To test this hypothesis, we calculated band-limited power and spectrograms of local field potentials in cortical area 24a aligned to spontaneous activations of distinct LC ensembles. A diversity of state modulations occurred around activation of different LC ensembles, including a typical activated state with increased highfrequency power as well as other states including decreased high-frequency power. Thus-in contrast to the stereotypical activated brain state evoked by en masse LC stimulation-spontaneous activation of distinct LC ensembles is associated with a multitude of cortical states.