Full-spectrum Sun-induced chlorophyll fluorescence (SIF) offers profound physiological insights into plant functional status compared to single-band SIF. We propose a Fourier series-based method (FSM) for retrieving full-spectrum SIF, aiming to address the limitations of existing methods, such as reliance on reflectance training datasets and the limited spectral range of retrieved SIF spectrum. The core principle of the FSM involves modeling reflectance as a wavelength-dependent function, which can be approximated by successive summations using high-order expansions of the Fourier series. The performance of the FSM was thoroughly evaluated through a combination of simulations and field measurements. The findings illustrate FSM’s capability to achieve high-precision full-spectrum SIF retrieval, with an average relative root-mean-square error (RRMSE) of 2.468% based on synthetic data. Moreover, the corresponding RRMSE values in the O2-A and O2-B bands, at 1.1% and 3.724%, respectively, indicate accuracy comparable to the spectral fitting method (SFM) and advanced FSR (aFSR) methods and superior to the SpecFit method. In the field full-spectrum SIF retrieval, FSM exhibited improved reflectance reconstruction and produced more reasonable results for the diurnal variation of full-spectrum SIF. The diurnal comparison of single-band SIF at both Italian and German sites further highlights the close alignment between FSM-retrieved SIF and the SFM SIF, with$R^{2}$values exceeding 0.96 and a maximum RMSE of 0.118 mW/m2/sr/nm. Conversely, the aFSR method encountered challenges stemming from an under-representation of the training dataset, resulting in the maximum RMSE at the Italian site reaching 0.506 mW/m2/sr/nm, along with a minimum$R^{2}$of 0.809. The FSM demonstrates the promising potential for full-spectrum SIF retrieval, accompanied by fewer limitations.

The monthly mean land surface temperature (MMLST) reflects more stable intra- and interannual temperature variations, and therefore, it has a wider range of applications than instantaneous land surface temperature (LST). This study aimed to generate a high-resolution global MMLST product by temporally upscaling the Moderate Resolution Imaging Spectroradiometer (MODIS) 1-km instantaneous LST. First, six current methods were comprehensively evaluated using cross-validation technology. These six methods are the cross combinations of two temporal aggregation schemes: the average by observations (ABO) and average by days (ABD), and three conversion models: the diurnal temperature cycle model (DTC), the simple average of two instantaneous LSTs (TSA), and a weighted average model for multiple instantaneous LSTs (MWA). The analysis with measurements from 235 flux stations worldwide revealed that the choice of conversion model considerably affected the overall retrieval accuracy, whereas the influence of the aggregation scheme was minor. From the conversion model standpoint, MWA performed best, followed by DTC, and finally TSA; this order remained the same even if DTC and TSA were improved with mean bias correction. Notably, the errors of ABDMWA decreased as the number of daily mean LST (NOD) increased, whereas the errors of ABOMWA were not related to NOD. Accordingly, we deduced that the optimal strategy for estimating MMLST is using ABOMWA when NOD is$\ge 20$. Subsequently, we adopted this combination method to process MODIS instantaneous LSTs and produced a global 1-km MMLST dataset for the years 2003–2020. The validation showed a satisfactory accuracy with a root mean square error (RMSE) of 1.6 K. The intercomparison with MMLSTs from geostationary (GEO) satellites (containing complete LST daily cycle) presented a good agreement (biases < 0.3 K and STDs < 2 K). Compared with atmospheric infrared sounder (AIRS) L3 monthly standard physical retrieval (AIRS3STM) product which had the same temporal span, the newly generated product exhibited a high consistency in reflecting temporal variations of global temperature. Most importantly, it had a prominently better ability to retrieve spatial details of temperature variations due to its higher resolution. Our new method and product show promising prospects for applications in global change studies, where accurate spatially resolved MMLST data are one of the fundamental geophysical variables required.

Accurate retrieval of land surface temperature (LST) over urban areas is of great significance for urban thermal environment monitoring. In previous studies, most of the urban LST retrieval methods were developed based on the assumption of a flat surface without considering the influence of urban 3-D geometry structure, which has a significant impact on the retrieval accuracy of LST over urban areas. In this study, a radiative transfer equation (RTE)-based single-channel method was developed to retrieve LST with urban geometry effect correction from the Landsat 8 thermal infrared (TIR) data in band 10. The increase in adjacent radiance from the surrounding pixels and the decrease in atmospheric downwelling radiance caused by urban geometry structure were taken into account in this method. Because it is difficult to directly validate the retrieval accuracy of LST over urban areas usingin situLST measurements, the performance of the RTE-based LST retrieval method was evaluated via comparing brightness temperature (BT) at the top of the atmosphere (TOA) simulated by the discrete anisotropic radiative transfer (DART) model and the urban RTE over three subregions. There is a good agreement between BT at the TOA simulated by the DART model and the urban RTE, with a root-mean-squared error (RMSE) of less than 0.25 K. The variations in LST retrieved with urban geometry effect correction over different local climate zones (LCZs) were analyzed. In general, built-up LCZs have relatively higher LST than land cover LCZs. The differences between LST retrieved without/with urban geometry effect correction over different LCZs are greater than 0.2 K. The largest average LST difference over built-up LCZs is approximately 0.9 K, whereas that over land cover LCZs is approximately 0.65 K. LST retrieved without/with urban geometry effect correction was used to calculate urban heat island intensity (UHII) in terms of the LCZ-based method. The results indicate that UHII calculated from LST with urban geometry effect correction is lower than that calculated from LST without urban geometry effect correction, with an average difference of approximately 0.5 K.

Relative humidity (RH) is one of the key variables for understanding the water, energy, and carbon exchange between the Earth and the atmosphere. Traditional methods for deriving RH from remotely sensed data usually require ground meteorological observations or are limited to clear-sky conditions, thereby making it a significant challenge to obtain spatially complete RH under all-sky conditions, especially over the regions with sparse meteorological instruments for observation. To this end, a new approach for deriving all-sky RH entirely based on Moderate Resolution Imaging Spectroradiometer (MODIS) data was proposed in the present study. Two key assumptions in the approach under cloudy conditions are that the actual water vapor is linearly related to the total precipitable water vapor (PWV) and that air temperature is linearly related to land surface temperature (LST). Results from a total of 30 AmeriFlux stations proved the aforementioned assumptions based on MODIS data collected over a study period of three years from 2009 to 2011. For different aridity conditions, RH retrieval revealed reasonable accuracy with a root-mean-square error (RMSE) of approximately 15.3% over an arid and semiarid region, whereas a comparable RMSE of 17.0% was obtained over a humid area. Further results also indicated that the aforementioned linear relationships were generally temporally stable, thereby indicating that the proposed method can be used to obtain all-sky RH at a regional or global scale entirely based on MOD06_L2-derived LST and MOD05_L2-derived PWV data given that the assumed linear relationships can be easily determined by historical MOD07_L2-derived atmospheric profiles.

Accurate estimation of land surface temperature (LST) is crucial for ecological environment monitoring and climate change studies in mountainous areas. The current LST retrieval algorithms were developed without accounting for the topographic effect, which can only be used to retrieve LST over relatively flat surfaces. Due to the impact of 3-D structure of mountainous surfaces, rugged terrain makes the processes of thermal radiation more complex. In this study, a radiative transfer equation (RTE)-based single-channel algorithm was proposed to retrieve LST with topographic effect correction from the Landsat 8 thermal infrared (TIR) data in mountainous areas. This algorithm accounts for the changes in the thermal radiation components in the TIR RTE caused by the topographic effect. According to the analysis of simulation data, sky-view factor (SVF), atmospheric water vapor content, surface emissivity of target pixel, and average LST of the surrounding terrain have significant influence on the magnitude of the topographic effect. The differences between the LST retrieved without/with topographic effect correction from the Landsat 8 TIR data are related to SVF. The topographic effect should be taken into account in the LST retrieval algorithm when SVF is smaller than 0.7. The largest LST difference of approximately 1 K occurs in the deep valley. The results indicate that LST without topographic effect correction could be overestimated to be as high as 1 K. Due to a lack ofin situLST measurements, the performance of the LST retrieval algorithm in mountainous areas was only evaluated by comparing the brightness temperature (BT) at the top of the atmosphere (TOA) simulated by the DART+MODTRAN model and the TIR RTE over mountainous surfaces at three subregions. There is a good consistency between BT at the TOA simulated by the DART+MODTRAN model and the TIR RTE over mountainous surfaces at the three subregions, with a root-mean-squared error (RMSE) of less than 0.23 K.

Evapotranspiration (ET) is one of the most critical parameters in water- and energy-related domains. Two basic assumptions with respect to soil-moisture variation have been widely investigated for the retrieval of ET based on the trapezoid methods. Specifically, soil moisture within the surface and root-zone layers was assumed to vary synchronously in most of the earlier analyses. However, several recent investigations assumed that soil moisture within the upper soil layer should be dried up before the root-zone layer is stressed. To this end, the retrieval of ET under different aridity conditions over North American grasslands was investigated with the two assumptions, and the estimated ET was assessed using the flux data collected from eight AmeriFlux sites. Based on the available data from 2002 to 2018, results showed that the “asynchronous-assumed” method can obtain better ET estimates than the “synchronous-assumed” method over semiarid and subhumid areas, whereas the “synchronous-assumed” method can obtain better ET estimates in humid areas. Moreover, because of the different closure techniques used for the ET correction, no consistent conclusions could be found for the arid conditions to determine which trapezoid was better. Specifically, it was found that the cases of surface soil with zero water availability that were defined by the asynchronous-assumed trapezoid method rarely occur, even in arid areas, which indicated that the critical boundary that determines whether the root-zone layer begins to be water-stressed may need to be redefined.

Land surface temperature (LST) is a key parameter for many fields of study. Currently, LST retrieved from satellite thermal infrared (TIR) measurements is attainable with an accuracy of about 1 K for most natural flat surfaces. However, over urban areas, TIR measurements are influenced by 3-D structures and their radiation that could degrade the performance of existing LST retrieval algorithms. Therefore, quantitative models are needed to investigate such impact. Current 3-D radiative transfer models are generally based on time-consuming numerical integrations whose solutions are not analytical, and are therefore difficult to exploit in the methods of physical retrieval of LST in urban areas. This article proposes an analytical TIR radiative transfer model over urban (ATIMOU) areas that considers the impact of 3-D structures and their radiation. The magnitude of this impact on TIR measurements is investigated in detail, using ATIMOU, under various conditions. Simulations show that failure to acknowledge this impact can potentially introduce a 1.87-K bias to the ground brightness temperature for street canyon whose ratio “wall height/road width” is 2, wall and road temperature is 300 K, wall emissivity is 0.906, and road emissivity is 0.950. This bias reaches 4.60 K if road emissivity decreases to 0.921, and road temperature decreases to 260 K. ATIMOU is also compared to the discrete anisotropic radiative transfer (DART) model. Small mean absolute error of 0.10 K was found between the models regarding the simulated ground brightness temperatures, indicating that ATIMOU is in good agreement with DART.

The change of global thermal environment plays an important role in land surface processes. In this study, global thermal environment was analyzed using the vertically polarized brightness temperature at 36.5 GHz. The daily brightness temperature from 2003 to 2010 were decomposed using the annual temperature cycle (ATC) model, and the annual cycle parameters (ACPs) were obtained. The results show that the brightness temperature decreases with the increasing latitudes respectively for the northern hemisphere and the southern hemisphere. The land covered by vegetation is colder than the desert and barren. Some plateaus lead to lower brightness temperature than surrounding areas. In addition, the atmospheric and ocean circulation also affect global brightness temperature. The ACPs from brightness temperature can generally characterize the global thermal environment.

Passive microwave data are much less affected by clouds than TIR data for the retrieval of land surface temperature (LST), providing its unique advantages in global mapping of LST. In this study, a physically-based algorithm for LST retrieval was applied to AMSR2 global brightness temperature data. The performances of this algorithm applied on different land cover types were further evaluated against nighttime MYD11A1 thermal infrared LST products. The results showed that (i) the overall accuracy of the algorithm is about 5.42 K by root mean square error (RMSE) and 2.99 K by bias against MODIS LST during nighttime; (ii) the algorithm overestimates the LST over all land types. The overestimation is most evident over barren/sparsely vegetated surfaces. The algorithm shows that the algorithm has a robust performance comparing with MODIS LST and could be applied to estimate LST effectively.

Evapotranspiration (ET) plays a key role for energy transfer and water circulation in the biosphere, lithosphere, hydrosphere, cryosphere and atmosphere. In present study, spatially complete ET over the Heihe river basin, Northwest of China, was estimated from the synergistic use of MODIS (MODerate-resolution Imaging Spectroradiometer) data and CLDAS (China Meteorological Administration Land Data Assimilation) gridded meteorological data from June 1 to September 15 in 2012. For the estimation of ET over clear-sky pixels, a pixel-to-pixel pattern of land surface temperature (LST)-vegetation index (VI) feature space was developed where meteorological data were used to determine the dry and wet edges for each pixel; whereas the traditional Penman-Monteith equation was implemented to obtain ET over clouds pixels. Finally, ground ET measurements collected at two sites (corn and orchard) were used to evaluate the estimated results, root mean square error (RMSE) of 77.2W/m2and 74.9W/m2can be obtained for the two sites, respectively, indicating that spatially complete ET can be derived from currently available satellite images and meteorological data.

Dust storms greatly affect the environment and the resources in arid and semi-arid areas. In this paper, we study the process of an Asian dust storm in the north of China in 2005 by means of Terra and Aqua MODIS data, which could monitor dust storms both in daytime and nighttime. We use three thermal infrared (TIR), 8.5 mum, 11 mum, and 12 mum, and a tri-spectral technique presented by Ackerman (1997). We find that combining brightness temperature differences BT8.5- BT11 with negative BT11- BT12 can detect dust storms in China. Furthermore, the pre-processing of images, such as cloud screening, is important for dust storm monitoring. Compared to daytime visible and near infrared RGB images, the TIR method predicted dust storm very well and could also obtain useful information in nighttime. We also analyzed dust storm motion by four MODIS images within 24 hours. Our results showed that the dust storm motions from inter land to ocean and relies on the main wind direction. This information is very important for studying the processes of dust storm.