Reducing Uncertainties in Applying Remotely Sensed Land Use and Land Cover Maps in Land-Atmosphere Interaction: Identifying Change in Space and Time

Land use and land cover (LULC) data are a central component of most land-atmosphere interaction studies, but there are two common and highly problematic scale mismatches between LULC and climate data. First, in the spatial domain, researchers rarely consider the impact of scaling up fine-scale LULC data to match coarse-scale climate datasets. Second, in the temporal domain, climate data typically have sub-daily, daily, monthly, or annual resolution, but LULC datasets often have much coarser (e.g., decadal) resolution. We first explored the effect of three spatial scaling methods on correlations among LULC data and a land surface climatic variable, latent heat flux in China. Scaling by a fractional method preserved significant correlations among LULC data and latent heat flux at all three studied scales (0.5◦, 1.0◦, and 2.5◦), whereas nearest-neighbor and majority-aggregation methods caused these correlations to diminish and even become statistically non-significant at coarser spatial scales (i.e., 2.5◦). In the temporal domain, we identified fractional changes in croplands, forests, and grasslands in China using a recently developed and annually resolved time series of LULC maps from 1982 to 2012. Relative to common LULC change (LULCC) analyses conducted over two-time steps or several time periods, this annually resolved, 31-year time series of LULC maps enables robust interpretation of LULCC. Specifically, the annual resolution of these data enabled us to more precisely observe three key and statistically significant LULCC trends and transitions that could have consequential effects on land-atmosphere interaction: (1) decreasing grasslands to increasing croplands in the Northeast China plain and the Yellow river basin, (2) decreasing croplands to increasing forests in the Yangtze river basin, and (3) decreasing grasslands to increasing forests in Southwest China. Our study not only demonstrates the importance of using a fractional spatial rescaling method, but also illustrates the value of annually resolved LULC time series for detecting significant trends and transitions in LULCC, thus potentially facilitating a more robust use of remotely sensed data in land-atmosphere interaction studies.


Introduction
Human activities have transformed a large proportion of the planet's land surface [1] through processes such as the deforestation of tropical forests, as well as intensified agricultural land use and urbanization in China and India. Human-induced land use and land cover change (LULCC) can alter surface roughness, surface wetness, the partitioning of surface energy between sensible and latent heat fluxes, and terrestrial carbon storage [2][3][4][5][6]. These changes are increasingly becoming a focus of concern because of their potential to influence the climate system [7,8] and, as a consequence, the Intergovernmental Panel on Climate Change (IPCC) has emphasized the importance of understanding the climate response to LULCC at local, regional, and global scales [9]. In particular, identifying and quantifying LULCC is crucial for a better understanding of land-atmosphere interaction and thereby climate change and variability [6].
Efforts to map land use and land cover (LULC) and its changing patterns using remotely sensed data have been a focus of much attention over the last three decades [10][11][12][13][14][15]. Remote sensing can potentially provide not only accurate and repeatable global LULC information, but also time series data that can be used to map change [16,17]. However, remotely sensed LULC maps usually have very different spatial resolutions from that of climate data, including atmospheric and oceanic variables. For both physical and observational reasons, climate data generally have much coarser spatial resolution than remotely sensed LULC data. For example, National Center for Environmental Prediction (NCEP)/ National Center for Atmospheric Research (NCAR) reanalysis, European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA- 40), and Climate Research Unit Time-series (CRU TS) high-resolution gridded datasets have a resolution of 2.5 • by 2.5 • , 1.125 • by 1.125 • , and 0.5 • by 0.5 • , respectively. While anthropogenic forces and other drivers (e.g., natural disturbance from fires or storms, sharp gradients of riparian vegetation) can dramatically alter LULC at fine scales, the atmosphere is often more well mixed, with properties changing more gradually across spatial scales. Simply resampling LULC maps into the same resolution as climate data is problematic, as LULCC is complex, with heterogeneous patterns that may not be evident in simple measures, such as dominant change type, in coarse resolution data.
Problems of spatial scale mismatches are well known to the remote sensing community. For example, Woodcock and Strahler [18] proposed local variance as a framework for examining the effects of scale and spatial resolution in remotely sensed images. They found that traditional spectral classification that produced a single class for each pixel was appropriate where the spatial resolution of the imagery was much finer than the objects in the scene. In contrast, for images where the spatial resolution was not sufficient to resolve the objects, a mixture model, which estimated the proportion of the classes within a pixel, was required [18]. In subsequent work, a more rigorous treatment was developed using variograms [19]. Cao and Lam [20] thoroughly reviewed scale and resolution effects that were relevant to geographic information system (GIS) and remote sensing, and usefully identified four conceptualizations of spatial scale: cartographic scale, geographic (observational) scale, operational scale, and measurement (resolution) scale. In a more recent study, Peng et al. [21] compared several downscaling methods for remotely sensed soil moisture products, including a satellite-based fusion method, a method using geoinformation data, and model-based methods. They concluded that each method has its own advantages and disadvantages, and none of the methods can be applied everywhere across the world without any calibration or improvements [21]. In land-atmosphere interaction studies, the potentially problematic scale effects of the spatial resolution of LULC data and issues of rescaling have been largely overlooked. For example, the Community Earth System Model (CESM) adopted LULC data from Moderate Resolution Imaging Spectroradiometer (MODIS) at 0.05 • (approximately 5 km) resolution [22], while the original spatial resolution of MODIS LULC data at that time was 1 km [23]. It is therefore important to identify the effect of different spatial scaling methods for LULC data on land-atmosphere interaction.
Another uncertainty in land-atmosphere interaction studies is due to insufficient LULC time series data to detect long-term LULCC patterns. Cao et al. [24] emphasized that longer time periods of LULC data are needed to quantify land-atmosphere interaction. Currently, there are few continuous annual LULC data, such as MODIS land cover data covering from 2001 onward [25] and European Space Agency (ESA) climate change initiative (CCI) land cover maps covering from 1992 to 2015 [26]. The primary method of extracting LULCC information from remotely sensed data has been through separate classifications of imagery of two different dates, which are overlaid to obtain change information [12]. For example, Rawat and Kumar [27] used 1990 and 2010 Landsat Thematic Mapper (TM) images to classify the land surface of Hawalbagh block, Uttarakhand, India, and thus mapped LULCC over a 20-year period. Similarly, Liu et al. [28] classified cropland based on 1990 and 2000 Landsat TM/ Enhanced Thematic Mapper Plus (ETM+) data to quantify the changing spatial patterns of cropland in China from 1990 to 2000. However, LULCC information from just two-time steps or even several time periods may obscure important overall trends in land cover change, especially in regions where land cover change is dynamic. For example, temporary fluctuations in a landscape where land cover change is common, such as occurs with the logging of short-rotation forests [29], may obscure important long-term trends in land cover change. Change that is non-significant over long time periods could bias analyses that use LULCC as an input or initial forcing for the biophysical and biogeochemical processes of land-atmosphere interaction [30][31][32][33][34][35]. Considering this limitation in recent LULCC detection studies, a spatiotemporal analysis using a long-term time series of LULC maps is an effective option. Such an approach can identify statistically significant spatiotemporal LULCC, as well as provide information regarding the types of LULCC and the dominant transitions among the LULC types.
In this study, we addressed these uncertainties related to the spatial and temporal resolutions of LULC data. Firstly, we used two regions in China ( Figure 1): a relatively homogenous region (North China plain), and a more heterogeneous region (Sichuan basin) to explore the effect of LULC spatial resolution on the relationship of LULC data and the land surface climatic variable of latent heat flux. We aggregated the original resolution of the LULC maps from 0.0833 • to 0.5 • , 1.0 • , and 2.5 • , which are the commonly used spatial resolutions of climatic variables. At each resolution, we investigated how the correlation of LULC data and latent heat flux is affected by three spatial aggregation methods. We used nearest-neighbor and majority-aggregation methods, which preserve the LULC maps as nominal maps, and a fractional method, which transforms the nominal LULC maps into fractional maps. Using the optimal scale and scaling method identified in these experiments within the spatial domain, we then explored the temporal changes from 1982 to 2012 for croplands, forests, and grasslands, which are the dominant LULC types in China. Specifically, we used a newly developed 31-year times series of annual LULC maps [36,37] to ask: where are the statistically significant LULCC regions in China based on the long-term time series LULC maps, and what are the dominant transitions among the different types of LULC?

Materials
The primary dataset used in this study is the continuous annual LULC maps of China, covering the period from 1982 to 2012. Figure 1 shows the LULC map of 2000 as an example. The data were generated by Reference [37] using a random forest classification of 19 phenological metrics derived from the Advanced Very High Resolution Radiometer (AVHRR) Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index third generation (NDVI3g) data. The phenological metrics include the start of growing season, end of growing season, maximum and minimum NDVI values, and so on. The classifier was trained using reference data derived from pixels within the MODIS land cover (MCD12Q1) maps that did not change over the period of 2001 to 2010. The LULCC maps have a 0.0833 • (approximately 8 km) spatial resolution and 13 classes, as listed in the legend of Figure 1. The overall accuracy of a simplified eight-class version of the 2012 LULC map was estimated to be 73.8% (refer to He, Lee, and Warner [37] for more information). In order to explore the variation in the relationship between land cover data and climatic variable, we used LULC maps in every six (or five for the last interval) years from the first year of the study period (1982,1988,1994,2000,2006, and 2011) for spatial rescaling experiments.
As a climatic variable for exploring the effects of rescaling, we chose latent heat flux, because it links the land surface condition to the atmosphere. Latent heat flux data was obtained from the FLUXNET-Multi-Tree Ensemble (MTE). FLUXNET-MTE datasets were estimated using a machine learning approach called model tree ensemble based on FLUXNET measurements, a long-term remotely sensed monthly fraction of absorbed photosynthetically active radiation (fAPAR) dataset, near surface air temperature from Climatic Research Unit (CRU), precipitation data from the Global Precipitation Climatology Center (GPCC), an estimate of the top of the atmosphere shortwave radiation, and information on land cover [38].
The land cover map used in the FLUXNET-MTE dataset did not vary over time [38]. In the MTE model, the land cover data were only used to stratify the data, and did not act as predictor variables for deriving data layers such as the latent heat flux in the regression equations [39]. The global estimates of latent heat flux were similar to other independent estimates [40][41][42] and the R-squared was 0.92 when correlated with catchment water balances [43]. The FLUXNET-MTE dataset has been widely used in hydrology and land-atmosphere interactions studies [43][44][45][46][47]. The spatial resolution of the latent heat flux datasets is 0.5 • by 0.5 • , over the period of 1982 to 2011. For consistency with the LULC data, latent heat flux during August in 1982August in , 1988August in , 1994August in , 2000August in , 2006, and 2011 was chosen for the analysis. We chose August because in that month there is a strong relationship between the proportion of croplands and latent heat flux. This is because for croplands August is a time of peak growth and high evapotranspiration, which is directly related to latent heat flux [48].
Although the latent heat flux derived from FLUXNET-MTE can be influenced by land cover types [49][50][51], our objective is to investigate how the relationships between latent heat flux and LULC data derived from different spatial scaling methods (i.e., fractional method, nearest-neighbor method, and majority-aggregation method) change over three different resolutions (i.e., 0.5 • , 1.0 • , and 2.5 • ). In doing so, we anticipate identifying the best spatial rescaling method, which can preserve the relationships between latent heat flux and LULC data among different scales.

Spatial Scaling Methods
To explore the suitability of different approaches for scaling LULC maps for land-atmosphere interaction studies, we compared three rescaling methods: nearest neighbor, majority aggregation, and a fractional method. The effectiveness of these methods was compared for the aggregation of LULC maps of China with an original 0.0833 • spatial resolution, which we rescaled to 0.5 • , 1.0 • , and 2.5 • spatial resolutions.
Nearest-neighbor and majority-aggregation methods, which are commonly used for categorical data such as LULC maps, preserve the nominal nature of LULC maps. In the nearest-neighbor method, the value of an output pixel is determined by the pixel nearest to it in the input data [52]. Majority-aggregation assigns the most frequently occurring value in the input data to the new pixel [37]. The fractional method is a proportional estimate of land cover for each class at the coarser resolution. Specifically, the spatial resolution of the LULC map is aggregated from 0.0833 • to 0.5 • , 1.0 • , and 2.5 • by calculating the proportion of the new coarse-scale pixel respectively covered by each class within a 6 × 6, 12 × 12, and 30 × 30 pixel region in the fine-scale data.
We chose bilinear interpolation to resample the latent heat flux data from the original 0.5 • to 1.0 • and 2.5 • spatial resolutions, because it is the method commonly used in the climate community to regrid interval data. In bilinear interpolation, the output pixel value is estimated through the linear interpolation of the four pixel values in two orthogonal directions within the input data [53].

Exploring the Relationship between LULC Data and Latent Heat Flux at Different Spatial Resolutions
For this section, we focused on croplands, because croplands are the dominant LULC type in the North China plain and the Sichuan basin ( Figure 1). For fractional LULC maps, we explored relationships between LULC data and latent heat flux by exploring the correlation between the fraction of croplands and latent heat flux. The correlation coefficient is [54]: where X is the fraction of croplands at different spatial scales in each pixel, and Y is the corresponding August latent heat flux value for each pixel. The significance of the correlation coefficient was tested by the Student's t-test.
For nominal LULC maps, the relationship between LULC data and latent heat flux was tested by the Wilcoxon rank sum test. The Wilcoxon rank sum test was used to identify the significance of the difference in the distributions of the latent heat flux values in croplands compared to non-croplands. The Wilcoxon rank sum test is a non-parametric test and is suitable for small samples [55]. We chose the Wilcoxon rank sum test because our sample size was proportionally reduced as the spatial scale coarsened.

Linear Regression Trend Analysis
We tested the usefulness of the fractional rescaling method for evaluating trends in land cover change. We aggregated our 31-year LULC maps using the fractional method, and in doing so, generated fractional maps for croplands, forests, and grasslands, respectively. These 31-year fractional maps were summarized through linear regression trend analysis and spatial pattern correlation analysis to determine the long-term LULCC patterns.
Linear regression analysis is a statistical method for analyzing the relationship between two or more variables by evaluating the degree to which one variable can be predicted or explained by the others [56]. The simple linear regression model is: Parameters a and b are the regression coefficients, which are estimated by a least squares method [57]; e is the regression residual. In this study, we were interested in b, the slope of the regression line, which characterized how the Y variables (fraction of croplands, forests, and grasslands, respectively) change over time (X). The significance of the trends was estimated using a Student's t-test. A trend analysis was performed to determine the spatiotemporal changes in croplands, forests, and grasslands for each grid cell independently, and regions where LULCC is significant were thereby identified. Statistically non-significant regions were masked out, and consequently only pixels significant at the 5% level were shown in the maps.

Spatial Pattern Correlation Analysis
Within four regions where statistically significant changes were identified, we performed spatial pattern correlation analysis to investigate the transitions among the LULC types (e.g., grasslands transitioning to croplands, or croplands to forests). The correlation coefficient here is the same as Equation (1). However, the X and Y have a different meaning. Here, X and Y are the trend values (b in the Equation (2)) of croplands, forests, and grasslands at each grid cell, respectively. The significance of correlation coefficient was quantified using an adjusted Student's t-test, accounting for spatial autocorrelation through a correction of the degrees of freedom of the samples [58].

Categorical and Factional LULC Maps at Different Spatial Resolutions in the North China Plain and the Sichuan Basin
The North China plain is a relatively homogenous region, with croplands being dominant. In contrast, the Sichuan basin is a relatively heterogeneous region, with croplands, mixed forests, and grasslands being prevalent (Figure 1). The overall pattern of croplands in the North China plain and the Sichuan basin is evident in both categorical LULC maps and fractional maps at all three spatial scales ( Figure 2). However, due to coarsening of the spatial scale, the heterogeneous patterns of croplands are smoothed, especially at the 2.5 • scale. Generally speaking, nearest-neighbor, majority-aggregation, and fractional methods capture the croplands patterns better at the finest scale, 0.5 • , with notable smoothing at the scales of 1.0 • and 2.5 • .

Relationships of Categorical and Fractional LULC Data with Latent Heat Flux at Different Spatial Resolutions in the North China Plain and the Sichuan Basin
The correlation coefficients for the association of fraction of croplands and latent heat flux in both the North China plain and the Sichuan basin change minimally as the spatial scale coarsened for all six years, and in each case are significant at the 1% level ( Table 1). The difference in mean correlation coefficient of the six years among the three spatial resolutions within each of these two regions is quite small (≤0.08). The significant positive correlations are associated with the increased evapotranspiration during the peak growing time of crops. In general, the correlations in the North China plain are larger than those in the Sichuan basin, which may be due to the existence of more croplands in the North China plain (Figure 1). For the relationships between categorical LULC data and latent heat flux, the distribution of latent heat flux in croplands and non-croplands is different at the 1% significance level for all six years at both the 0.5 • and 1.0 • scales in the North China plain and the Sichuan basin for both nearest-neighbor and majority-aggregation methods (Tables 2 and 3). At the 2.5 • scale, the significance level varies among years in the North China plain. The difference of the distribution is significant in 2006 and 2011 at the 1% level, in 1994 at the 5% level, in 1988 and 2000 at the 10% level, and not significant in 1982 for both nearest-neighbor and majority-aggregation methods ( Table 2). In the Sichuan basin, the difference of the distribution is not significant at the 2.5 • resolution for all six years, excepting 1994 for the nearest-neighbor method and 2011 for the majority-aggregation method, at the 5% level (Table 3). In addition, the significance level decreases more in the more heterogeneous Sichuan basin than in the relatively homogenous North China plain.  In summary, the fractional method preserves the significant relationship between LULC data and latent heat flux for all six years and all three resolutions in both the North China plain and the Sichuan basin, while the nearest-neighbor and majority-aggregation methods cause these correlations to diminish and even become statistically insignificant at coarser scales in some years. Therefore, we recommend rescaling using fractional maps in land-atmosphere studies. Although the different aggregation methods for LULC data lead to different relationships between LULC data and latent heat flux, the bilinear interpolation of latent heat flux may also influence the relationships. It would be interesting to know if the different relationships are due to only the LULC data aggregation or the aggregation methods for both the LULC data and climate data; this needs to be explored further.

Fractional Maps of Croplands, Forests, and Grasslands
Based on the aforementioned analyses, the fractional method preserves the significant relationship between LULC data and latent heat flux, and the finer scale (i.e., 0.5 • resolution in this study) best captures the spatial patterns. Therefore, we used the fractional method to resample LULC maps to the spatial resolution of 0.5 • to identify spatiotemporal changes of croplands, forests, and grasslands from 1982 to 2012 in China. Figure 3 shows the mean fractional maps of croplands, forests, and grasslands averaged over the period of 1982 to 2012. The maps are color coded by the fraction of croplands, forests, and grasslands in each 0.5 • grid cell. Figure 3a indicates that croplands dominate in the Northeast China plain, the Hai river basin, the Yellow river basin, the Yangtze river basin, and the Sichuan basin (refer to Figure 1 for the locations of these regions). All of these regions have been recognized as important Chinese agricultural regions since the 11th century [28,[59][60][61].
Forests dominate along the borders of the Northeast China, including the mountain ranges of Daxing'anling, Xiaoxing'anling, and Changbaishan, as well as southeast and southwest mountainous regions (Figure 3b). This overall spatial pattern of the distribution of forests in China is consistent with the findings of Yin et al. [62] and Fenning [63]. Grasslands are distributed in a belt, stretching from the Northeast China plain, adjacent to Inner Mongolia, to the Tibetan Plateau. This region covers the major grassland regions of China, including the Hulun Buir and Xilin Gol grasslands in Inner Mongolia and the Naqu grassland in Tibet (Figure 3c). The additional grasslands identified in this study, such as the Narat and Bayinbuluke grasslands in northern Xinjiang province, in far northwestern China, were also identified by Kang et al. [64].

Spatiotemporal Changes of Croplands, Forests, and Grasslands during the Last Three Decades
The linear regression trend analysis indicates how the fraction of croplands, forests, and grasslands changed from 1982 to 2012 at each grid cell (Figure 4). The color bar represents the slope of the linear regression model over time (i.e., years) for locations where the change is statistically significant at the 5% level. Croplands significantly increased in the Northeast China plain and the Yellow river basin and decreased in the Yangtze river basin (Figure 4a). The highest rate of increase, greater than 1.5%/year, occurred in the Northeast China plain. The increase in croplands in the Northeast China plain may be due to the increase in population, and the associated increase in demand for food [65]. Furthermore, policy regulations such as scrapping the agriculture tax in 2006 and increasing agricultural investment in Northeast China to improve cropland yields in 2008, as well as scientific and technological progress in agriculture including increased access to fertilizer and irrigation and the mechanization of agriculture, may also have all contributed to the increase in croplands [66]. There are also some small isolated areas of significantly increased croplands in the northern part of Xinjiang province, in northwestern China, which may be attributed to the successful promotion of modern agronomic technology [61,67]. The decreased croplands in the Yangtze river basin may be due to increasing forests (Figure 4b), as well as the expansion of industry in that region, which has resulted in arable lands being converted to built-up areas [61]. The overall changing patterns of croplands, especially the increase in North China and the decrease in South China, are consistent with those reported by Liu et al. [62]. In contrast to the trend for croplands, forests significantly decreased along the borders of Northeast China, and increased in the Hai river basin, the Yangtze river basin, and Southwest China (Figure 4b).
The increasing trend of forest cover in the Yangtze river basin may be due to the national forest protection project after the catastrophic flood of 1998 [61] and the construction of the forest shelter-belt system [65], which converted other types of LULC, such as croplands, into forest. The increasing trend of forests in southern mountainous areas may be attributed to the "Grain for Green" program and related afforestation activities [65,68]. The reclamation of forests to control sandstorms in the Hai river basin, including the cities of Beijing and Tianjin, is a likely cause of the increasing forests in areas neighboring these two cities.
Grasslands significantly decreased in the Northeast China plain, the Yellow river basin, and Southwest China, and increased in some parts of Western China (Figure 4c). Grassland degradation, which is closely associated with desertification, is a long-term environmental concern in China. The main factors causing grassland degradation are over grazing and cultivation, processes that are common in, for example, Inner Mongolia [69]. However, a revision of the Grassland Law in 2002, which controls all aspects of livestock grazing, may yet result in the recovery of grasslands in Western China (Figure 4c) [69]. The pattern of decreasing grasslands also seems to correspond with an increase in croplands and forests. Figure 5 summarizes the relative frequency of fraction of croplands, forests, and grasslands for grid cells where the LULCC trend is significant. The histograms show that for all three LULC types, regions with the lowest fractions of LULC (i.e., 0-20%) have the largest number of grid cells with significantly changed LULC, whether the trend is decreasing or increasing, indicating that change tends to occur in cover classes that are a minority proportion of the landscape, rather than the dominant fraction. For croplands and forests, increasing trends are the most common, even where forests and croplands already dominate. For grasslands, decreasing trends dominate, especially for areas where grasslands constitute a low percentage of the land cover even before the additional loss of grassland. In contrast, in areas of moderate percentages of grassland cover (40-80%), increasing and decreasing trends somewhat balance, resulting in only a small net loss of grasslands in such areas.  Figure 4. Blue color represents a significantly decreasing trend and red denotes a significantly increasing trend.

Transitions between Croplands, Forests, and Grasslands
To examine the underlying relationships among the changes in croplands, forests, and grasslands during the last three decades, we performed spatial pattern correlation analysis, drawing only on grid cells with trends significant at the 5% level. We studied the three LULC types in the selected regions (shown as rectangles in Figure 4): the Northeast China plain (Rectangle 1 in Figure 4), the Yellow river basin (Rectangle 2), the Yangtze river basin (Rectangle 3), and Southwest China (Rectangle 4). Scatterplots and spatial correlation coefficients for the linear regression trends (%/year) of croplands, forests, and grasslands are shown in Figure 6.
The spatial correlation coefficients between croplands, forests, and grasslands in the four regions show very different associations. In the Northeast China plain and the Yellow river basin, grid cells with increasing croplands generally correspond to decreasing grasslands, with a significant negative correlation of r = −0.934 in the Northeast China plain (Figure 6b) and r = −0.894 in the Yellow river basin (Figure 6e). Trends of increased forests are significantly correlated with trends of decreased grasslands in the Yellow river basin (Figure 6f) and the Northeast China plain (Figure 6c). From Figure 6, together with the spatial trends shown in Figure 4, we can draw the general inference that in the Northeast China plain and the Yellow river basin the dominant transition was the replacement of grasslands by croplands. In addition, in the Yellow river basin, grasslands also appear to have been converted to forest cover, though the pattern is less strong than that observed with croplands. In the Yangtze river basin, the only significant association is between decreasing croplands and increasing forests (r = −0.545) (Figures 4 and 6g). In Southwest China, the transition between forests and grasslands has the strongest association (r = −0.668) (Figure 6l), followed by a transition between croplands and forests (r = −0.433) (Figures 4 and 6j).  Figure 4. Each data point represents a trend significant at the 5% level over the 31-year for a particular 0.5 • grid cell.

Conclusions
Remotely sensed LULC maps have been widely used in land-atmosphere interaction studies, and this research identified key uncertainties related to scale mismatches in the spatial and temporal domains. Both categorical LULC maps and fractional maps capture the overall pattern of croplands in the relatively homogenous North China plain and the comparatively heterogeneous Sichuan basin at all of the three spatial resolutions studied (0.5 • , 1.0 • , and 2.5 • ). The significant positive correlation between the fraction of croplands and latent heat flux in August of 1982August of , 1988August of , 1994August of , 2000August of , 2006, and 2011 at the three spatial resolutions is quite similar in the two study areas, although the difference in r values between 0.5 • and 2.5 • is greater in the Sichuan basin than in the North China plain. For the categorical LULC maps, the distribution of latent heat flux values over croplands and non-croplands is significantly different at the 1% level at the 0.5 • and 1.0 • spatial resolutions for both nearest-neighbor and majority-aggregation methods in North China plain and Sichuan basin for all of the six years. However, at the 2.5 • spatial resolution, the significance level varies among years. Nevertheless, both the nearest-neighbor and majority-aggregation methods are particularly problematic in regions with more heterogeneous and fine-grained spatial patterns of LULC. In contrast, for both the relatively homogeneous North China plain and heterogeneous Sichuan basin, the fractional maps preserved the relationships between LULC data and latent heat flux at all of the spatial scales. Therefore, we recommend rescaling using fractional maps in observational land-atmosphere studies. In summary, the fractional scaling method preserves significant correlations among LULC data and latent heat flux at all three studied scales (0.5 • , 1.0 • , and 2.5 • ), while nearest-neighbor and majority-aggregation methods cause these correlations to diminish and even become statistically non-significant at coarser scales.
Based on the optimal spatial scaling method, this study explored spatiotemporal changes in croplands, forests, and grasslands from 1982 to 2012 through 31-year LULC maps, which were aggregated to 0.5 • using the fractional method. The annual LULC maps enable a time-series analysis of changes in LULC, which is the key attribute compared to the previous LULCC studies of China. Croplands significantly increased in the Northeast China plain and the Yellow river basin, and decreased in the Yangtze river basin. Forests significantly decreased along the border of Northeast China, and increased in the Yangtze river basin and Southwest China. Grasslands significantly decreased in the Northeast China plain, the Yellow river basin, and Southwest China, and increased in some parts of Western China. The analysis of transitions among significant changes in croplands, forests, and grasslands shows that decreasing grasslands in the Northeast China plain and the Yellow river basin are significantly associated with increasing croplands. Similarly, increasing forests in the Yangtze river basin are significantly associated with decreasing croplands. In Southwest China, decreasing grasslands and, to a lesser extent, decreasing croplands, are significantly correlated with the increasing forests.
The time-series analysis of annual LULC maps in China suggests a complex pattern of spatiotemporal changes in croplands, forests, and grasslands from 1982 to 2012. The analysis facilitates the identification of regions where change is extensive and statistically significant, such as North China, where croplands have increased significantly, and South China, where forest have increased significantly. Compared to previous studies, which relied on overlaying LULC maps covering a short-time period [70], the statistically significant patterns of spatiotemporal changes from our study provide a robust approach for capturing LULCC in China.
The identified LULCC patterns in this study have the potential to influence the regional climate in China. The increasing forest in the Yangtze river basin and Southwest China, associated respectively with decreasing of croplands and grasslands, is likely decreasing the region's albedo, because of the higher albedo of croplands and grasslands compared to forest. A decreased albedo would induce higher net radiation, and consequently promote higher temperatures [7]. On the other hand, the conversions from croplands and grasslands to forests may reduce surface air temperature due to the cooling effects from the higher evapotranspiration of forests [71,72]. The transitions of LULC types within the regions of extensive significant LULCC revealed by this study can potentially alter surface heat and moisture conditions and thereby induce changes in the regional climate system. Previous studies have explored the climatic impacts of LULCC in China. For example, Cao et al. [24] used the weather research and forecasting (WRF) model to investigate how different LULCC types affect regional climate in the agro-pastoral transitional zone of North China. Ma et al. [73] simulated afforestation impacts on the regional climate in Jiangxi province, China using WRF. Fu [32] applied a regional integrated environmental model system (RIEMS) to investigate the effects of human-induced land cover change on the East Asia monsoon. However, most previous studies have not been able to incorporate long-term temporal LULCC patterns due to lack of long-term time series LULC data. Incorporating long-term LULCC pattern, such as that produced in this study, could therefore be useful for future empirical and modeling studies investigating how LULCC affects regional climate.