Recent advances in technologies have lead to a vast influx of data on movements, based on discrete recorded position of animals or fishing boats, opening new horizons for future analyses. However, most of the potential interest of tracking data depends on the ability to develop suitable modelling strategies to analyze trajectories from discrete recorded positions. A serious modelling challenge is to infer the evolution of the true position and the associated spatio-temporal distribution of behavioural states using discrete, error-prone and incomplete observations. In this paper, a Bayesian Hierarchical Model (HBM) using Hidden Markov Process (HMP) is proposed as a template for analyzing fishing boats trajectories based on data available from satellite-based vessel monitoring systems (VMS). The analysis seeks to enhance the definition of the fishing pressure exerted on fish stocks, by discriminating between the different behavioural states of a fishing trip, and also by quantifying the relative importance of each of these states during a fishing trip. The HBM approach is tested to analyse the behaviour of pelagic trawlers in the Bay of Biscay. A hidden Markov chain with a regular discrete time step is used to model transitions between successive behavioural states (e.g., fishing, steaming, stopping (at Port or at sea)) of each vessel. The parameters of the movement process (speed and turning angles) are defined conditionally upon the behavioural states. Bayesian methods are used to integrate the available data (typically VMS position recorded at discrete time) and to draw inferences on any unknown parameters of the model. The model is first tested on simulated data with different parameters structures. Results provide insights on the potential of HBM with HMP to analyze VMS data. They show that if VMS positions are recorded synchronously with the instants at which the process switch from one behavioural state to another, the estimation method provides unbiased and precise inferences on behavioural states and on associated movement parameters. However, if the observations are not gathered with a sufficiently high frequency, the performance of the estimation method could be drastically impacted when the discrete observations are not synchronous with the switching instants. The model is then applied to real pathways to estimate variables of interest such as the number of operations per trip, time and distance spent fishing or travelling. (C) 2010 Elsevier B.V. All rights reserved.