Protoplanetary disks around young Sun-like stars are the cradles of the vast majority of detected exoplanets. Probing these disks at multiple spatial scales is key to uncovering how planets form. The inner astronomical unit, the star-disk interaction region, is of utmost importance because most detected exoplanets occupy this zone.

We aim to spatially and spectrally resolve the inner disk and star-disk interaction region of the M0.3 T Tauri star DO Tau by combining two complementary techniques.

We used high-resolution near-infrared spectra from CFHT/SPIRou to constrain the magnetospheric star-disk interaction process and optical long-baseline interferometry with ESO VLTI/GRAVITY to determine the sizes of the K-band continuum and Br γ line emitting regions. From the SPIRou spectra, we measured the veiling in the YJHK bands along with the equivalent widths of the HeI λ 1083, Paβ, and Br γ emission lines, from which we estimated the mass accretion rate. We were able to monitor the time variability of these quantities thanks to our long-sequence of observations over about 40 days. We fit the GRAVITY visibilities in the continuum and the differential quantities in the line with geometrical models to obtain the orientation and the size of the inner disk as well as the size and the on-sky displacement of the Br γ emitting region. We derived a mass accretion rate of ∼10 ⁻⁸ −10 ⁻⁷ M _⊙ yr ⁻¹ , which confirms that this ∼0.5 M _⊙ star is a strong accretor. The HI and HeI lines exhibit strong variability on a daily timescale, consistent with the burster classification of DO Tau derived from its K2 light curve. We report a periodic modulation of the intensity of the redshifted high-velocity wings of the Br γ line profile. The modulation occurs at the rotational period of the star (5.128 d), which suggests the existence of corotating magnetospheric funnel flows. We derived an upper limit of 0.35 on the ratio between the magnetospheric truncation radius and the disk corotation radius, indicative of an ordered unstable accretion regime. The size of the Br γ line emitting region obtained from GRAVITY is quite small ( R _Brγ = 0.011 au ∼ 1.3 R _* ), and it is much smaller than the K-band continuum emitting region ( R _K = 0.09 au ∼ 11 R _* ). Such a compact Br γ emission region suggests that most of the line flux originates from the magnetospheric accretion region and/or from an inner wind close to the magnetosphere-disk interface. The on-sky displacements of the blue and red Br γ line velocity channels suggest a rotation pattern of the emitting gas, as they appear to be nearly aligned along the position angle of the disk. The inclination we derived for the inner disk (∼45-55°) differs from that of the outer disk inferred from the ALMA continuum (∼30°). This points toward a misalignment or warp of the outer disk that may originate from the suspected past encounter with the neighboring HV Tau system. Based on combining high-resolution spectroscopy and long baseline interferometry, we find that the T Tauri star DO Tau appears to be a strong accretor undergoing magnetospheric accretion in an ordered unstable regime, with a Br γ line emitting region as compac