Functional MRI modeling is challenged by long-range coupling, non-linearity, and lack of detailed physiological information. I will review our progress in modeling the spatio-temporal dynamics of fMRI including hemodynamic deconvolution, blind deconvolution, and brain state decoding based on spatio-temporal kernel methods.