Brain parcellation is one of the most important issues in functional MRI (fMRI) data analysis. This parcellation allows establishing homogeneous territories that share the same functional properties. This paper presents a model-based approach to perform a subject-level parcellation into hemodynamic territories with similar hemodynamic features which are known to vary between brain regions. We specifically investigate the use of the Joint Parcellation-Detection-Estimation (JPDE) model initially proposed in [1] to separate brain regions that match different hemodynamic response function (HRF) profiles. A hierarchical Bayesian model is built and a variational expectation maximiza-tion (VEM) algorithm is deployed to perform inference. A more complete version of the JPDE model is detailed. Validation on synthetic data shows the robustness of this model to varying signal-to-noise ratio (SNR) as well as to different initializations. Our results also demonstrate that good parcellation performance is achieved even though the parcels do not involve the same amount of activation. On real fMRI data acquired in children during a language paradigm, we retrieved a parcellation along the superior temporal sulcus of the left hemisphere that matches the gradient of activation dynamics already reported in the literature.Read less <