In this paper, we propose a novel deep domain adaptation method based on graph neural network (GNN) for multitemporal hyperspectral remote sensing images. In GNN, graphs are constructed for source and target data, respectively. Then the graphs are utilized in each hidden layer to obtain features. GNN operates on graph structure and the relations between data samples can be exploited. It aggregates features and propagate information through graph nodes. Thus, the extracted features have an improved smoothness in each spectral neighborhood which is beneficial to classification. Furthermore, the domain-wise correlation alignment and classwise correlation alignment are jointly embedded in GNN network to achieve a joint distribution adaptation performance. By introducing the joint correlation alignment strategy in GNN, the extracted features can not only be aligned between domains but also have a superior discriminability in each domain. This domain adaptation network is named as joint correlation alignment based graph neural network (JCGNN). Experiments using multitemporal Hyperion and NCALM datasets demonstrate the effectiveness of the proposed method. Index Terms-Hyperspectral remote sensing, graph neural network, domain adaptation, classification Recently, deep learning has been applied to domain adaptation because of its excellent feature representation ability. Compared with shallow domain adaptation approaches, end-toend deep domain adaptation approaches transfer more effective knowledge by embedding adaptation modules in the network architecture. Maximum Mean Discrepancy (MMD) is a popular distribution alignment strategy, which is utilized in many deep domain adaptation methods. Tzeng et al. [17] proposed the deep domain confusion method, which firstly regularizes the single adaptation layer of deep neural network using linear-kernel MMD. Similar to MMD, correlation alignment could also measure the distribution distance between domains, and it attempts to align the second-order statistics of the source and target features. Sun et al. [18] extended the correlation alignment to deep architectures and proposed the correlation alignment for deep domain adaptation (D-CORAL) approach, so that a nonlinear transformation that aligns the correlations of features between the two domains is obtained. Besides,