Few lossy compression methods consider both the global and local multidimensional coupling correlations, which could lead to information loss in data compression. Here we develop an adaptive lossy compression method, Adaptive-HGFDR, to capture both the global and local variation of multidimensional coupling correlations and improve approximation accuracy. The method can achieve good compression performances for most flux variables with significant spatiotemporal heterogeneity.