Abstract

Surface albedo is a crucial parameter in land surface radiation budget. As bias exists between the model simulated and observed surface albedo, data assimilation is an important method to improve the simulation results. Moreover, surface albedo is associated with the wavelength of the sunlight. So, solar radiation partitioning is important to parameterize the surface albedo. In this paper, the moderate resolution imaging spectroradiometer- (MODIS-) retrieved direct visible, direct near-infrared, diffuse visible, and diffuse near-infrared surface albedos were assimilated into the integrated urban land model (IUM). The solar radiation partitioning method was introduced to parameterize the surface albedo. Based on the albedo data from MODIS and the solar radiation partitioning method, the surface albedo data set for the Beijing municipal area was generated. Based on the surface albedo data set and the IUM, the impacts of the surface albedo on the surface radiation budget were discussed quantitatively. Surface albedo is inversely proportional to the net radiation. For urban areas, after assimilation, the annual average net radiation decreases about 5.6%. For cropland, grassland, and forest areas, after assimilation, the annual average net radiations increase about 20.2%, 24.3%, and 18.7%, respectively.

1. Introduction

Surface albedo is a crucial parameter in land surface radiation budget [1]. The parameterization of surface albedo is important in land-atmosphere interaction and coupling [18]. Surface albedo is not only associated with the soil moisture, soil color, and solar zenith angle [915] but also associated with the wavelength of the solar radiation [16, 17]. Although progress has been made in surface albedo simulation by land surface models, biases [10, 18] still exist inevitably between the simulation and observation.

Compared with the model simulation, surface albedo data retrieved from remote sensing images have their advantages [1923] and are used and assimilated [21, 2427] into the land surface models. Among them, surface albedo from the moderate resolution imaging spectroradiometer (MODIS) [18, 2434] is one of the most widely used products. Unfortunately, for most of the assimilation methods [2527, 31], the effect of the solar radiation partition [17, 35] is not included. In fact, the surface albedos for visible and near-infrared solar radiations are quite different [17]; solar radiation partitioning is important to parameterize the surface albedo.

In this paper, the MODIS-retrieved direct visible, direct near-infrared, diffuse visible, and diffuse near-infrared surface albedos were assimilated into the integrated urban land model (IUM). The solar radiation partitioning method [17] was introduced to parameterize the surface albedo to consider the effect of the wavelength of the solar radiation. Based on the surface albedo data from MODIS and the solar radiation partitioning method, the surface albedo data set for the Beijing municipal area was generated. Based on the surface albedo data set and the IUM, the impacts of the surface albedo on surface radiation budget were discussed quantitatively.

2. Data and Method

2.1. Study Area

The study area of this paper is the Beijing municipal area. The research region is located at 39.3°–41.1°N, 115.2°–117.6°E. The terrain tilts from northwest to southeast (Figure 1(a)). The dominate land cover category (LUC) in each grid is based on MODIS 1 km resolution data in 2017 (Figure 1(b)). The LUCs were transferred from International Geosphere-Biosphere Programme (IGBP) classification to United States Geological Survey (USGS) classification using the relationship listed in Table 1. The Chinese Academy of Sciences 325 m-high Meteorology and Environmental Observation Tower (325 m tower) was located in downtown Beijing, and the longitude and latitude of the 325 m tower are 116.37 E and 39.97 N, respectively (Figure 1(a)).

2.2. Data

The 325 m tower data were used for the surface albedo data comparison. The tower is located in downtown Beijing, and the altitude of the foot of the tower is 49 m. The radiation fluxes including the upward and downward shortwave and longwave radiation are measured using the radiometer at the 47-meter height. The time ranges of the observed upward and downward shortwave radiation are from March to October 2018. The observed surface albedo is calculated as follows:where is the surface albedo; is the upward shortwave radiation (W m−2); and is the downward shortwave radiation (W m−2).

The albedos for direct visible, direct near-infrared, diffuse visible, and diffuse near-infrared solar radiations are from MODIS in 2015. The temporal and spatial resolution of the MODIS albedo data are 8d and 1 km, respectively. The MODIS albedo products were spliced and clipped firstly; then, they were interpolated linearly from 8d to 1d.

The atmospheric forcing data used to drive the land surface model are originated from the Global Land Data Assimilation System (GLDAS) [36]. The initial field data are also originated from GLDAS. The temporal resolution of GLDAS is 3 h. The spatial resolution of GLDAS data is 0.25°, which was interpolated to 1 km by using a bilinear interpolation method.

2.3. Land Model

We employed the integrated urban land model (IUM) [37] to study the impact of the surface albedo on radiation budget. IUM was developed based on the common land model (CoLM) [38]. IUM integrates the land surface models for urban and natural land surfaces.

2.4. Surface Albedo Parameterization

The surface albedo for bare soil [15] in the IUM could be described as follows:where is the solar elevation angle (o); is the volumetric soil moisture; and are the albedos for visible and near-infrared solar radiations, respectively; and is the albedo for bare soil. The surface albedo for urban in the IUM could be described as follows:where and are the albedos for visible and near-infrared solar radiations, respectively, and is the albedo for urban. The surface albedo for vegetated surfaces is the same as that in the CoLM. A two-stream approximation method is used to parameterize it. The grid point surface albedo is calculated as follows:where is the albedo for vegetation and , , and are the fractional cover for soil, urban, and vegetation, respectively.

In this paper, the moderate resolution imaging spectroradiometer- (MODIS-) retrieved direct visible, direct near-infrared, diffuse visible, and diffuse near-infrared surface albedos were assimilated into the integrated urban land model (IUM). The solar radiation partitioning method was introduced to parameterize the surface albedo. The assimilation algorithm is as follows:where , , , and are direct visible, diffuse visible, direct near-infrared, and diffuse near-infrared solar radiations (W m−2), respectively; they are calculated based on Weiss and Norman [35] and corrected to consider the effect of cloud cover [17]. , , , and are albedos for direct visible, direct near-infrared, diffuse visible, and diffuse near-infrared solar radiations, respectively; they are retrieved from MODIS. The surface albedo data set for the Beijing municipal area was generated based on the assimilation algorithm and MODIS data. The spatial and temporal resolution of the surface albedo data are 1 km and 1d, respectively.

3. Results and Discussion

First, in order to validate the accuracy of the surface albedo simulations before and after assimilation, we compared them with the observed surface albedos in the 325 m tower site. Figures 2(a) and 2(b) are the monthly and seasonal average surface albedo simulations before and after assimilation and the observations. The biases between the monthly average surface albedo simulations before and after assimilation and the observations are listed in Table 2. In the whole assimilation period, after assimilation, the albedo simulations are larger than those before assimilation. After assimilation, the biases are decreased apparently for all the months and seasons compared with those before assimilation.

Then, we compared the spatial distribution of the annual and seasonal average surface albedo simulations before (Figure 3) and after (Figure 4) assimilation in the Beijing municipal area. The differences (after assimilation minus before assimilation) between the simulation results before and after assimilation are plotted in Figure 5. For urban areas, the surface albedo simulations after assimilation are usually larger than those before assimilation. However, for other LUCs, the surface albedo simulations after assimilation are usually smaller than those before assimilation.

Then, we compared the annual and seasonal average surface albedo simulations before and after assimilation for different LUCs. Here, to simplify the comparison, we merged the LUCs of the research region into four main LUCs, that is, urban, cropland, forest, and grassland (Table 3). For most of the urban areas, the annual average surface albedo simulations after assimilation are larger than those before assimilation (Figure 6). However, for cropland, grassland, and forest areas, the annual average surface albedo simulations after assimilation are smaller than those before assimilation (Figure 6). Positive correlations exist for all the four LUCs, and the correlation coefficients between the simulations before and after assimilation for urban, cropland, grassland, and forest LUCs are 0.447, 0.182, 0.557, and 0.570, respectively.

We also compared the monthly average surface albedos for these four LUCs (Figure 7). For urban areas, the surface albedo simulations after assimilation are larger than those before assimilation for all the 12 months. However, for cropland, grassland, and forest areas, the surface albedo simulations after assimilation are smaller than those before assimilation for all the 12 months. For cropland, grassland, and forest areas, the surface albedo simulations before assimilation are relatively small in summer and large in winter; on the contrary, the surface albedo simulations after assimilation are relatively large in summer and small in winter. For cropland, grassland, and forest areas, the surface albedo simulations before and after assimilation are relatively large in November because of the snowfall and December because of the snow cover (Figure 8).

Surface albedo is associated with the upward solar radiation, so it is a crucial parameter in surface radiation budget. Figure 9 is the monthly average net radiation simulations before and after assimilation for the four main LUCs. Figure 10 is the scatter plot of the annual average albedo and net radiation simulation differences before and after assimilation for the four main LUCs (after assimilation minus before assimilation). Albedo is inversely proportional to the net radiation. For urban areas, where the surface albedo simulations after assimilation are larger than those before assimilation, the net radiation simulations after assimilation are smaller than those before assimilation. However, for other LUCs, the net radiation simulations after assimilation are larger than those before assimilation.

Net radiation is not only associated with the upward solar radiation but also with the upward longwave radiation. As the upward longwave radiation is dependent on the surface radiative temperature, we compared the monthly average surface radiative temperature for these four LUCs (Figure 11). Figure 12 is the scatter plot of the annual average albedo and upward longwave radiation simulation differences before and after assimilation for the four main LUCs (after assimilation minus before assimilation). Albedo is inversely proportional to the upward longwave radiation. As the solar albedo increases, the solar radiation received by the land decreases; as a result, the surface radiative temperature and the upward longwave radiation decrease too.

The monthly average surface radiation budget simulations before and after assimilation for the four LUCs are shown in Figure 13. A quantitative analysis was performed for the impacts of the surface albedo on the surface radiation budget (Table 4). For urban areas, compared with that before assimilation, after assimilation, the annual average net radiation decreases about 5.6%. For cropland, grassland, and forest areas, compared with those before assimilation, after assimilation, the annual average net radiations increase about 20.2%, 24.3%, and 18.7%, respectively.

4. Summary

In this paper, the MODIS-retrieved direct visible, direct near-infrared, diffuse visible, and diffuse near-infrared surface albedos were assimilated into the IUM. The solar radiation partitioning method was introduced to parameterize the surface albedo. Based on the albedo data from MODIS and the solar radiation partitioning method, the surface albedo data set for the Beijing municipal area was generated. The surface albedo data were validated in the 325 m tower which was located in downtown Beijing. The result indicates that, after assimilation, the simulation results of the surface albedo are improved apparently.

Surface albedo is a crucial parameter in the surface radiation budget. The results indicate that the surface albedo is inversely proportional to the net radiation. For urban areas, compared with that before assimilation, after assimilation, the annual average net radiation decreases about 5.6%. For cropland, grassland, and forest areas, compared with those before assimilation, after assimilation, the annual average net radiations increase about 20.2%, 24.3%, and 18.7%, respectively.

In the near future, the observed four components of the solar radiation should be used to validate the solar radiation partitioning methods. More observational sites should be used to validate the surface albedo data set. The generation algorithm of the surface albedo data set will be extended to a larger area such as the whole of China or the world.

Data Availability

The observational data from the 325 m tower station are from the Institute of Atmospheric Physics, Chinese Academy of Sciences, and are available in the supplementary materials.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant nos. 41875125 and 41705086. The author thanks the Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences (CAS), for 325 m Beijing meteorology tower observation data; MODIS for the LUC (http://glcf.umd.edu/data/lc/) and the surface albedo data (https://ladsweb.modaps.eosdis.nasa.gov/search); and NASA for the GLDAS data (http://disc.sci.gsfc.nasa.gov/hydrology/data-holdings).

Supplementary Materials

The observational data from the Chinese Academy of Sciences 325 m-high Meteorology and Environmental Observation Tower (325 m tower) are provided as the Supplementary Materials. (Supplementary Materials)