Ultra-high spatial resolution fMRI often uses 3D-EPI for data acquisition. However, 3D-EPI can be susceptible to artefacts from inter-shot variability, such as motion and physiological noise. Here, we present an improved hybrid radial-Cartesian 3D EPI sampling strategy, TURBINE, to generate ultra-high spatial resolution fMRI data at conventional scan times at 7T. TURBINE enables intrinsic correction of global and z-dependent shot-to-shot phase variations, and self-navigated motion correction. Using these features, along with a temporally regularized reconstruction, we demonstrate robust BOLD activation in 0.67mm isotropic resolution (16mm slab, TRvol = 2.32s) and 0.8x0.8x2.0mm (whole-brain, TRvol = 2.4s) acquisitions.