We study the evolution of quantum correlations, quantified by the geometric discord, of two excitonic quantum dot qubits under the influence of the phonon environment. We show that the decay of these correlations differs substantially form the decay of entanglement. Instead of displaying sudden death type behavior, the geometric discord shows a tendency to undergo transitions between different types of decay, is sensitive to non-local phase factors, and may already be enhanced by weak environment-mediated interactions. Hence, two-qubit quantum correlations are more robust under decoherence processes, while showing a richer and more complex spectrum of behavior under unitary and non-unitary evolution.PACS numbers: 03.65. Ta, 63.20.Kd, 78.67.Hc The study of quantum correlations in realistic systems has, for a long time, been limited to the study of entanglement, due to the fact that straightforward methods of calculating the amount of correlations in a two-qubit system have only been available for some entanglement measures, such as the concurrence [1,2] or negativity [3,4]. Although entanglement itself is a very important resource for a number of applications [5], including quantum computation, quantum cryptography, or teleportation, separability (the lack of entanglement) does not automatically exclude the presence of quantum correlations [6]. This is in particular the reason why quantum computation models relying on mixed, separable (not entangled) states [7][8][9] are possible.The quantum discord [10, 11] is a measure of quantum correlations (see, however Refs. [11,12] for Holevo-type and thermodynamic based measures) which captures correlations beyond entanglement; it is defined as the difference of two classically equivalent formulas for mutual information and is non-negative. Due to the null volume of the set of zero-discord states [13], discord measures are not expected to undergo sudden death which is characteristic for entanglement evolutions [14][15][16]. The geometric measure of the discord describes the amount of correlations in a quantum system by finding the minimal Hilbert-Schmidt distance to the set of zero-discord states [17]. Recently, a lower [17] and an upper [18] bound on the geometric discord which can be calculated from a two-qubit density matrix have been found, which substantially simplifies the problem of studying the evolution of the quantum discord and opens the path for a qualitative and quantitative description of the decay of quantum correlations in realistic open quantum systems.In this paper we study the evolution of the lower and upper bounds of the geometric discord of two exciton quantum dot (QD) qubits interacting with an open phonon environment in order to capture the physical aspects of decoherence effects on quantum correlations. The interactions present in the system and the resulting dynamics are well understood. The experimentally observed evolution on picosecond timescales [19,20] can be described by pure dephasing within the independent boson model [20,21]. Super-Oh...