In this paper, we propose a Multilayer Markovian model for change detection in registered aerial image pairs with large time differences. A Three Layer Markov Random Field takes into account information from two different sets of features namely the Modified HOG (Histogram of Oriented Gradients) difference and the Gray-Level (GL) Difference. The third layer is the resultant combination of the two layers. Thus we integrate both the texture level as well as the pixel level information to generate the final result. The proposed model uses pairwise interaction retaining the sub-modularity condition for energy. Hence a global energy optimization can be achieved using a standard min-cut/ max flow algorithm ensuring homogeneity in the connected regions.