Continuum damage formulations are commonly used for the simulation of diffuse fracture processes. Implicit gradient damage models are employed to avoid the spurious mesh dependencies associated with local continuum damage models. The C0-continuity of traditional finite elements has hindered the study of higher-order gradient damage approximations. In this contribution we use isogeometric finite elements, which allow for the construction of higher-order continuous basis functions on complex domains. We study the suitability of isogeometric finite elements for the discretization of higher-order gradient damage approximations