Coal fires are a common and serious problem in most coal-bearing countries. Thus, it is very important to monitor changes in coal fires. Remote sensing provides a useful technique for investigating coal fields at a large scale and for detecting coal fires. In this study, the spreading direction of a coal fire in the Wuda Coal Field (WCF), northwest China, was analyzed using multi-temporal Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) thermal infrared (TIR) data. Using an automated method and based on the land surface temperatures (LST) that were retrieved from these thermal data, coal fires related to thermal anomalies were identified; the locations of these fires were validated using a coal fire map (CFM) that was developed via field surveys; and the cross-validation of the results was also carried out using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) thermal infrared images. Based on the results from longtime series of satellite TIR data set, the spreading directions of the coal fires were determined and the coal fire development on the scale of the entire coal field was predicted. The study delineated the spreading direction using the results of the coal fire dynamics analysis, and a coal fire spreading direction map was generated. The results showed that the coal fires primarily spread north or northeast in the central part of the WCF and south or southwest in the southern part of the WCF. In the northern part of the WCF, some coal fires were spreading north, perhaps coinciding with the orientation of the coal belt. Certain coal fires scattered in the northern and southern parts of the WCF were extending in bilateral directions. A quantitative analysis of the coal fires was also performed; the results indicate that the area of the coal fires increased an average of approximately 0.101 km² per year.

The reflectance anisotropy effect on albedo retrieval was evaluated using the Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution functions (BRDFs) product, and archetypal BRDFs. Shortwave-band archetypal BRDFs were established, and validated, based on the Anisotropy Flat indeX (AFX) and time series MODIS BRDF over tile h11v03. To generate surface albedo, archetypal BRDFs were used to fit simulated reflectance, based on the least squares method. Albedo was also retrieved based on the least root-mean-square-error (RMSE) method or normalized difference vegetation index (NDVI) based prior BRDF knowledge. The difference between those albedos and the MODIS albedo was used to quantify the reflectance anisotropy effect. The albedo over tile h11v03 for day 185 in 2009 was retrieved from single directional reflectance and the third archetypal BRDF. The results show that six archetypal BRDFs are sufficient to represent the reflectance anisotropy for albedo estimation. For the data used in this study, the relative uncertainty caused by reflectance anisotropy can reach up to 7.4%, 16.2%, and 20.2% for sufficient, insufficient multi-angular and single directional observations. The intermediate archetypal BRDFs may be used to improve the albedo retrieval accuracy from insufficient or single observations with a relative uncertainty range of 8–15%.

As sensor parameters and atmospheric conditions create uncertainties for at−sensor radiation detection, radiometric consistency among satellite images is difficult to achieve. Relative radiometric normalization is a method that can improve multi−image consistency with accurate pseudo−invariant features (PIFs), especially over large areas or long time series satellite images. Although there are algorithms that manually or automatically select PIFs, the spatial mismatch of satellite images can affect PIF extraction, particularly with artificial pixels. To alleviate this problem, we proposed to use Landsat−8 OLI as the reference image and Sentinel−2A as the subject image, to apply pseudo−invariant features−based algorithms with polygon features through the single−band and multiple−band regression. Compared to pseudo−invariant point features, hyperspectral library, and histogram matching approaches, the results demonstrate the superiority of the proposed algorithms with correlation coefficients of 0.9948 and 0.9945, and an RMSE of 0.0097 and 0.0095 with multiple− and single−band regression, respectively. We also found more accurate linear fitting and better shape matching through band scattering and reflectance frequency analysis. The proposed algorithms are a significant improvement in radiometric normalization, within artificial pixels, achieving spectral signature consistency.

Many previous studies have attempted to extract prior reflectance anisotropy knowledge from the historical MODIS Bidirectional Reflectance Distribution Function (BRDF) product based on land cover or Normalized Difference Vegetation Index (NDVI) data. In this study, the feasibility of the method is discussed based on MODIS data and archetypal BRDFs. The BRDF is simplified into six archetypal BRDFs that represent different reflectance anisotropies. Five-year time series of MODIS BRDF data over three tiles are classified into six BRDF archetype classes according to the Anisotropy Flat indeX (AFX). The percentage of each BRDF archetype class in different land cover classes or every 0.1-NDVI interval is determined. Nadir BRDF-Adjusted Reflectances (NBARs) and NDVIs simulated from different archetypal BRDFs and the same multi-angular observations are compared to MODIS results to study the effectiveness of the method. The results show that one land cover type, or every 0.1-NDVI interval, contains all the potential BRDF shapes and that one BRDF archetypal class makes up no more than 40% of all data. Moreover, the differences between the NBARs and NDVIs simulated from different archetypal BRDFs are insignificant. In terms of the archetypal BRDF method and MODIS BRDF product, this study indicates that the land cover or NDVI is not necessarily related to surface reflectance anisotropy.

Ciconia boyciana is a rare migrant bird in China and is listed as endangered species by IUCN. In order to identify their important stopover sites during migration, based on the Random Forest and the MaxEnt model, the occurrence records of Ciconia boyciana in spring migration and autumn migration in 2021 and 24 environmental variables were used to simulate the potentially suitable stopover sites. Compared hotspot stopover sites with current conservation area distribution, the strategies of the conservation establishment were proposed.

Considering that the complexity and dynamicity of CO2 emissions, the spatiotemporal distribution pattern of atmospheric CO2 and its drivers remain unclear. In this study, we used the Geographically Weighted Regression (GWR) method to analyze the comprehensive distribution of CO2 concentration in mainland China from 2015 to 2019. We considered the relationship between nine factors, including natural and human activities, and CO2 concentration. To clarify the correlation between CO2 concentration and drivers, we utilized Pearson’s correlation coefficient. Then, the GWR analysis revealed the spatial heterogeneity across provinces, which reflects the extent to which impact factors influence CO2 concentrations. Finally, we analysed CO2 concentration spatiotemporal variation characteristics and predicted future trends of CO2 concentration in 31 provinces in China. According to our research, GDP has a major impact on CO2 growth, while natural factors have a minor influence on CO2 concentration. Our study found significant regional differences in the effects of combined variables on CO2 concentrations, with monthly rotational patterns temporally and clustering of high growth rates spatially and CO2 concentration in mainland China will continue to steadily increase. The findings of this research are crucial for China’s future energy low-carbon transition and policy-making.