A RASSCF approach to simulate the O K-edge XANES spectra of uranyl is employed, utilizing three models that progressively improve the representation of the local crystal environment. Simulations successfully reproduce the observed three-peak profile of the experimental spectrum and confirm peak assignments made by Denning. The [UO2Cl4]2– model offers the best agreement with experiment, with peak positions (to within 1 eV) and relative peak separations accurately reproduced. Establishing a direct link between a specific electronic transition and peak intensity is complicated, as a large number of possible transitions can contribute to the overall peak profile. Furthermore, a relationship between oxygen character in the antibonding orbital and the strength of the transition breaks down when using a variety of orbital composition approaches at larger excitation energy. Covalency analysis of the U–O bond in both the ground- and excited-state reveals a dependence on the crystal environment. Orbital composition analysis reveals an underestimation of the uranium contribution to ground-state bonding orbitals when probing O K-edge core-excited states, regardless of the uranyl model employed. However, improving the environmental model provides core-excited state electronic structures that are better representative of that of the ground-state, validating their use in the determination of covalency and bonding.