## 5 Conclusions  
Arctic ecosystems have undergone great changes in the context of climate change. GPP is one of the most crucial indicators of the response of ecosystems to climate change. However, few efforts have been devoted to exploring the spatial variation and phenological characteristics of GPP in the Arctic. In response to this challenge, this study investigated  
Figure 16: Interannual variation trend of GPP (g C m\({}^{-2}\) year\({}^{-1}\)) distribution with latitude (**a**), elevation (m) (**b**) over the Arctic. The gray belt refers to the confidence interval of 95%, the blue line refers to the fit line, and the boxplots denote the distribution of the annual-averaged GPP with the latitude. (**c**) Land cover-dependent interannual variation trend of GPP (g C m\({}^{-2}\) year\({}^{-1}\)) over the Arctic. FOR, SHR, SAV, GRA, WET, SNO, and BAR denote the forests, shrublands, savannas, grasslands, permanent wetlands, permanent snow and ice, barren, and water bodies, respectively. the spatial distribution as well as the seasonal (phenological) and interannual variations of GPP in the Arctic using MOD17A2H. Furthermore, the GPP variation trends with land cover types, latitude, and elevation were also explored. In order to ensure that the results were reliable, the accuracy of MOD17A2H was first evaluated using in situ measurements from FLUXNET. This study found that MOD17A2H generally underestimates GPP over the land cover types investigated in this study, and its accuracy tends to be higher over low latitude but lower over high latitude. However, the overall accuracy suggests that MOD17A2H is consistent with the FLUXNET 2015 dataset (RMSE = 1.26 g C m\({}^{-2}\) d\({}^{-1}\), R = 0.8, Bias = \(-\)0.32), and MOD17A2H can represent the phenological characteristics of GPP (RMSE = 8.9 days, R = 0.62). Based on MOD17A2H, it was demonstrated that the maximum GPP occurred in July. In addition, the spatial distribution of GPP is related to land cover types; for example, forests and savannas have relatively high carbon storage capacity from April to August. By comparing the GPP variation with latitude and elevation, it was shown that GPP generally decreases as the latitude and elevation increase. However, the phenomenon is not evident for latitudes in the range (62\({}^{\circ}\) N, 66\({}^{\circ}\) N) and elevation lower than 700 m. The overall trend of GPP in the Arctic is greater than zero and is dominated by the variation of vegetation productivity in June. Furthermore, the response to climate change is different across these land cover types; for example, forests are most sensitive to climate warming. The distribution of the interannual trend in GPP across latitudes and elevations is consistent with the changes in GPP as a function of latitude and elevation. This study is helpful for understanding the spatiotemporal distribution characteristics of GPP over the Arctic as well as the response of ecosystems to climate change. Nevertheless, the results need to be validated with different satellite products. The number of sites used for validation is limited. Therefore, the presented conclusion about the accuracy of MOD17A2H may not be transferable to other regions. Another limitation is that only individual factors such as land cover type, latitude, and elevation were considered in this paper. Other factors, such as air temperature, precipitation, and snow, that are related to vegetation growth status need to be explored further. This will also be our focus in the future. Conceptualization, [PERSON] and J.W. ; methodology, D.M. and X.W. ; validation, X.L., C.M. and X.W. ; formal analysis, D.M. and X.W; investigation, J.W. ; resources, X.W. ; data curation, D.M. ; writing--original draft preparation, D.M. and X.W. ; writing--review and editing, D.M., X.W. and X.M. ; visualization, J.W. ; supervision, C.M. All authors have read and agreed to the published version of the manuscript.