Gadolinium-doped ceria or gadolinium-stabilized ceria (GDC) is an important technical material due to its ability to conduct O 2− ions, e.g., used in solid oxide fuel cells operated at intermediate temperature as an electrolyte, diffusion barrier, and electrode component. We have synthesized Ce 1-x Gd x O 2-y :Eu 3+ (0 ≤ x ≤ 0.4) nanoparticles (11−15 nm) using a scalable spray pyrolysis method, which allows the continuous large-scale technical production of such materials. Introducing Eu 3+ ions in small amounts into ceria and GDC as spectroscopic probes can provide detailed information about the atomic structure and local environments and allows us to monitor small structural changes. This study presents a novel approach to structurally elucidate europium-doped Ce 1-x Gd x O 2-y :Eu 3+ nanoparticles by way of Eu 3+ spectroscopy, processing the spectroscopic data with the multiway decomposition method parallel factor (PARAFAC) analysis. In order to perform the deconvolution of spectra, data sets of excitation wavelength, emission wavelength, and time are required. Room temperature, time-resolved emission spectra recorded at λ ex = 464 nm show that Gd 3+ doping results in significantly altered emission spectra compared to pure ceria. The PARAFAC analysis for the pure ceria samples reveals a high-symmetry species (which can also be probed directly via the CeO 2 charge transfer band) and a low-symmetry species. The GDC samples yield two low-symmetry spectra in the same experiment. High-resolution emission spectra recorded under cryogenic conditions after probing the 5 D 0 − 7 F 0 transition at λ ex = 575−583 nm revealed additional variation in the low-symmetry Eu 3+ sites in pure ceria and GDC. The total luminescence spectra of CeO 2-y :Eu 3+ showed Eu 3+ ions located in at least three slightly different coordination environments with the same fundamental symmetry, whereas the overall hypsochromic shift and increased broadening of the 5 D 0 − 7 F 0 excitation in the GDC samples, as well as the broadened spectra after deconvolution point to less homogeneous environments. The data of the Gd 3+containing samples indicates that the average charge density around the Eu 3+ ions in the lattice is decreased with increasing Gd 3+ and oxygen vacancy concentration. For reference, the Judd−Ofelt parameters of all spectra were calculated. PARAFAC proves to be a powerful tool to analyze lanthanide spectra in crystalline solid materials, which are characterized by numerous Stark transitions and where measurements usually yield a superposition of different contributions to any given spectrum.