## 1 Introduction  
Land surface albedo describes the fraction of incoming solar radiation reflected by the surface of the land. It influences the surface energy budget [1], and it is essential for global and regional estimation of energy and mass exchanges between the Earth's surface and the atmosphere [2; 3; 4; 5; 6]. An accuracy of 0.02-0.05 is required by climate models for global surface albedo [7; 8]. To monitor the spatio-temporal changes in land surface albedo, albedo products are routinely generated from satellite data, such as the Polarization and Directionality of the Earth's Reflectance (POLDER) [9; 10; 11; 12], the Medium Resolution Imaging Spectrometer (MERIS) [13], the Clouds and the Earth's Radiant Energy System (CERES) [11; 12], the Visible Infrared Imaging Radiometer Suite (VIIRS) [14; 15], and the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) [16]. The Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Earth Observation System (EOS) Terra and Aqua satellites routinely provided data to derive the land surface shortwave and visible albedos [17; 18; 19] used to calibrate and improve albedo parameterizations for land, weather, and climate models [20; 21; 22; 23; 24; 25; 26]. Assessment of the accuracy of these products is important because it is critical to the scientific community for various applications. Feedback from this activity will help improve the generation of these products [27]. Direct comparison with the ground-based observations of albedo values is the commonly used method to assess the accuracy of remote sensing albedo products [28; 29; 30]. The tower observation is compared directly with the satellite product [31; 32; 33] based on the assumption that the satellite and tower observation have the same footprint or the landscape is homogeneous. Researchers subsequently recognized that the land surface albedo varies strongly in space and across seasons, because of which land surface homogeneity is now examined prior to the \"point-to-pixel\" comparison. Only sites whose representativeness is adequate for satellite pixel scales are selected in the direct validation and \"heterogeneous\" sites are excluded [34]. The method most suited to a comparison for all kinds of landscapes is using higher spatial resolution data as intermediate data [27; 35; 36]. [PERSON], et al. [37] used airborne hyperspectral imagery to validate MODIS albedo product in snow-covered areas; [PERSON], et al. [10] used the convolved albedo onboard the Formosat-2 Taiwanese satellite as a reference to evaluate the newly released MCD43D product. However, since no global high-resolution albedo product (at a level of tens of meters) is available, validation using intermediate data has been conducted at a limited number of locations. Prior to evaluating the remote sensing product, identifying the spatial representativeness of the products is essential. It can help ground sampling point settings and match remote sensing data from multiple sources. MODIS albedo products are retrieved using observations covering a large area that depends on the view zenith angles. Although observations are weighted by angular coverage before albedo retrieval, the actual coverage of the pixel is always larger than the nominal spatial resolution [38; 39]. Efforts have been made to characterize the effective resolution of the MODIS gridded product [40]. [PERSON], et al. [10] used Formost-2 data at a resolution of 8 m to characterize the equivalent point spread function of MODIS albedo at a 1 km pixel. The Full Width at Half Maximum (FWHM), recognized as the effective resolution, has been confirmed to represent the footprint of MODIS data for accurate validation. [PERSON] et al. [41] used extensive time series data (2003-2014) at a large size of the linear natural target in the Netherlands to analyze the effective spatial resolution of the MODIS albedo product (the spatial scale the pixel represents).