In order to assess the linearity of the mechanisms subserving renal blood flow autoregulation, broad-band arterial pressure fluctuations at three different power levels were induced experimentally and the resulting renal blood flow responses were recorded. Linear system analysis methods were applied in both the time and frequency domain. In the frequency domain, spectral estimates employing FFT, autoregressive moving average (ARMA) and moving average (MA) methods were used; only the MA model showed two vascular control mechanisms active at 0.02-0.05 Hz and 0.1-0.18 Hz consistent with previous experimental findings [Holstein-Rathlou et al., Amer. J. Physiol., vol. 258, 1990.]. In the time domain, impulse response functions obtained from the MA model indicated likewise the presence of these two vascular control mechanisms, but the ARMA model failed to show any vascular control mechanism at 0.02-0.05 Hz. The residuals (i.e., model prediction errors) of the MA model were smaller than the ARMA model for all levels of arterial pressure forcings. The observed low coherence values and the significant model residuals in the 0.02-0.05 Hz frequency range suggest that the tubuloglomerular feedback (TGF) active in this frequency range is a nonlinear vascular control mechanism. In addition, experimental results suggest that the operation of the TGF mechanism is more evident at low/moderate pressure fluctuations and becomes overwhelmed when the arterial pressure forcing is too high.