Neural oscillations are the intrinsic characteristics of brain activities. Traditional electrophysiological techniques (e.g., the steady-state evoked potential, SSEP) have provided important insights into the mechanisms of neural oscillations in the high frequency ranges (>1 Hz). However, the neural oscillations within the low frequency ranges (<1 Hz) and deep brain areas are rarely examined. Based on the advantages of the low frequency blood oxygen level dependent (BOLD) fluctuations, we expected that the steady-state BOLD responses (SSBRs) would be elicited and modulate low frequency neural oscillations. Twenty six participants completed a simple reaction time task with the constant stimuli frequencies of 0.0625 Hz and 0.125 Hz. Power analysis and hemodynamic response function deconvolution method were used to extract SSBRs and recover neural level signals. The SSEP-like waveforms were observed at the whole brain level and at several task-related brain regions. Specifically, the harmonic phenomenon of SSBR was task-related and independent of the neurovascular coupling. These findings suggested that the SSBRs represent non-linear neural oscillations but not brain activations. In comparison with the conventional general linear model, the SSBRs provide us novel insights into the non-linear brain activities, low frequency neural oscillations, and neuroplasticity of brain training and cognitive activities.