Satellite-derived surface upward longwave radiation (SULR) is essential for monitoring the global surface radiation budget, ecological processes and climate change. However, the widely-used SULR products derived from thermal infrared (TIR) remote sensing exhibit spatial discontinuities because TIR signals cannot penetrate cloud cover. Conventional cloud-sky SULR estimation approaches often utilize post-processed reanalysis data as inputs, which could not meet the real-time requirement of operational system. This study proposes a lightweight cloud-sky SULR real-time estimation method for the Fengyun-4A (FY-4A) geostationary satellite using a Light Gradient Boosting Machine (LightGBM) model. The daytime cloud-sky SULR is estimated by applying the established relationship between auxiliary variables and clear-sky SULR to cloudy conditions, while the nighttime cloud-sky SULR values are estimated by applying the determined relationship between input variables and a publicly accessible, gap-filled SULR product. The model inputs include (1) spatial-temporal location record data, (2) multiple surface characteristic parameters generated from previous-year data, and (3) two categories of operational FY-4A radiation products, with both components being available in real-time. Validation against six Heihe Watershed Allied Telemetry Experimental Research (HiWATER) sites demonstrates that the reconstructed cloud-sky SULR achieves acceptable RMSE (MBE) values of 33.4 W/m2(1.5 W/m2) for daytime and 25.2 W/m2(4.7 W/m2) for nighttime conditions. Therefore, the proposed lightweight method could improve the spatial coverage of current FY-4A SULR product and further promote real-time SULR-related applications.

Thermal radiation directionality (TRD) describes the anisotropic signature in the thermal infrared domain, leading to significant uncertainties in current land surface temperature (LST) and surface upward longwave radiation (SULR) products. The kernel-driven model (KDM) is considered as the most potential tool to correct TRD effects due to its good tradeoff between physical accuracy and computational efficiency. However, the application of existing 4-parameter KDMs is limited due to the requirement of simultaneously ≥4 multi-angle observations. By combining a diurnal temperature cycle model, the time-evolving kernel driven model (TEKDM) had achieved significant TRD elimination effect over the overlapping region of two Geostationary Operational Environmental Satellite (i.e., GOES-16 and GOES-17) LST products. However, the performance of TEKDM in SULR TRD elimination is still not clear. In this letter, the TEKDM was extended to the overlapping region of FengYun-4A (FY-4A) and Himawari-8 satellite observations, in order to correct the directional SULR (SULRD) to hemispherical integrated SULR (SULRH). Then, a step-by-step SULR downscaling method based on multiple linear regression was conducted forSULRDandSULRH(from 4 km to 40 m). The downscaled SULRD and SULRH values were validated by the pyrgeometer-observed hemispherical SULR of 14 in-situ sites within a heterogeneous region. The validation result shows that the TEKDM could eliminate the TRD effect of FY-4A and Himawari-8 SULR with an RMSE decrease of 4.47 W/m2(16.5%) and 4.03 W/m2(15.1%), respectively. Therefore, the TEKDM also performs well for the TRD correction of SULR products of two geostationary satellites.

Radiative transfer models (RTMs) designed to reproduce the anisotropy of surface brightness temperature (BT) are particularly useful for applications on Earth’s energy budget when using remote sensing (RS) datasets. Despite the fact that several thermal infrared (TIR) RTMs have been developed, a quantitative analysis comparing the benefits and limits of these models remains necessary. Herein, three modeling frameworks (physical hybrid, analytical parameterization, and kernel driven) have been evaluated comparatively for homogeneous vegetation, a row-planted crop, and a sparse forest. Airborne measurements and the discrete anisotropy radiative transfer (DART) model simulations were retained as the benchmark. Forward modeling and inverse fitting schemes were proposed for the sake of comparison. Results reveal that: 1) in the forward modeling scheme, from airborne measurements, the hybrid model performs better with root-mean-squared errors (RMSEs) of$0.17~^{\circ }$C,$1.57~^{\circ }$C, and$0.38~^{\circ }$C for homogenous, row-planted vineyard, and sparse forest scenes, respectively; the analytical model appears similar performant ($0.17~^{\circ }$C,$0.40~^{\circ }$C) for the homogeneous and sparse forest scenes, but less performant ($2.39~^{\circ }$C) for the row-planted scene and 2) in the inverse fitting scheme, the uncertainties (95% of probability) of model coefficients and predicted directional anisotropies were considered. The kernel-driven model has fewer modeling constraints and statistically performs better for the homogeneous and sparse forest scenes with RMSEs of$0.07~^{\circ }$C and$0.19~^{\circ }$C, respectively, whereas it is less efficient for the row-planted scene with RMSE of$0.80~^{\circ }$C. This study highlights the differences in accuracy between models of different complexity and provides reference information for researchers to improve existing models and for users to choose their best modeling solution.

In order to correct the land surface temperature (LST) errors caused by angular effects to normalize the multisource data, various models have been proposed, especially the kernel-driven model shows good prospects for development. However, there are fewer studies on the validation of thermal radiation directionality (TRD) characteristics of features based on real measurements. In this article, a method based on unmanned aerial vehicle (UAV) observation of TRD characteristics of different surfaces is proposed to validate the TRD, which uses the circling flight mode to obtain the thermal radiation characteristics of features in multiple directions. Observations show that there is a significant hotspot effect in the three vegetation scenes, and the TRD and dispersion degree are different due to the differences in vegetation structure. For the realization of further applications of UAV-based methods for measuring TRD, the UAV-based angular observations are evaluated against the simulation results of the current kernel-driven model with high accuracy, and the results show that the correlation coefficients are all greater than 0.78, the coefficients of determination are greater than 0.60, and the root-mean-square error (RMSE) is less than 1. However, the ability of the model to extrapolate forward varies considerably in different vegetation scenes. The present study provides an effective and potential way to observe the TRD, which can lay a practical foundation for the theoretical study and model simulation of the TRD.

The accurate acquisition of land surface temperature (LST) on complex mountainous surfaces has always been a difficult problem and hot topic in thermal infrared remote sensing inversion, and the uncertainty caused by the radiation angle effect is one of the important factors hindering the accurate inversion of LST. Researchers have proposed a variety of models to simulate and eliminate the influence of the angle effect, among which, the kernel-driven model has a bright prospect of development. However, fewer studies have been conducted to observe the effects of terrain and land cover on thermal radiation directionality (TRD) properties based on measured data. This paper intends to carry out observations on a small spatial scale of a complex mountainous area using an unmanned aerial vehicle (UAV) to investigate the specific effects of different slopes, aspects, and land cover on the TRD characteristics. The measured results show that the intensity of thermal radiation anisotropy is actively correlated with the complexity of surface structure, and the influence of slope on TRD bias is in the form of a staged “S”, where its intensity increases slowly and then rapidly before slowing down again; the dispersion of thermal radiation in TRDs of different aspects is affected by the duration and intensity of solar radiation, and there is a time lag effect. Meanwhile, in order to evaluate the accuracy of different radiation directionality models, this paper further evaluates the currently more recognized kernel-drive model named LSF-Chen and the thermal equivalent slope kernel-driven (TESKD) model based on the TRD measurements from UAVs. The results show that the correlation coefficients of simulation results and measurements are all greater than 0.6, and the RMSEs are all less than 2K, and that the two methods both have a better simulation of the TRD effect, but the TESKD is better overall in terms of accuracy and methodological details. Through this study, a new method of applying UAVs to capture the thermal direction of the complex surface in mountainous areas is proposed, which provides methodological support for the extraction and accurate simulation of the TRD characteristics on the complex mountainous areas.

Land surface temperature (LST) is a vital parameter for energy budget and land surface process models. With the development of high spatial resolution remote sensing technologies, especially the wide application of unmanned aerial vehicle (UAV) in land surface observations, the acquisition of high-resolution thermal infrared (TIR) data has made the accurate extraction of varied LST products possible. However, the complex heterogeneous structure of land surface may cause adjacency effect, i.e., multiple scattering of radiation and geometric occlusion, which is inevitable in high-resolution UAV-based observations. The existing TIR temperature retrieval methods mainly focus on the one-dimensional (1D) scenarios, lacking consideration for the three-dimensional (3D) structure of actual ground objects, resulting in significant errors in extracting LST from high-resolution data. To address these limitations, we propose a novel retrieval method to eliminate the influence of multiple scattering caused by finer spatial resolution on the accuracy of LST. This method introduces a 3D radiative transfer model to account for the radiative transfer in a heterogeneous scenario, based on remotely sensed imagery, corresponding 3D structure of the ground surface and component attributes. An optimization strategy is adopted to iteratively converge toward physically consistent LST results. The proposed retrieval method is validated using the UAV-based TIR images and in-situ surface temperature measurement data obtained from the Huailai remote sensing test site in Hebei province, China. Simulated TIR image datasets of typical vegetation and urban scenes were additionally employed to validate the applicability at different values of key parameters. The factors influencing the adjacency effect were extensively analyzed. Results show that 1) the effects of adjacent objects on LST can result in an overestimate exceeding 0.9 K in cases of typical vegetation scene for TIR observations when the spatial resolution was finer than 1 m; 2) the proposed LST retrieval method based on a 3D radiative transfer model can significantly reduce the influence of adjacency effect on the accuracy of LST; 3) in addition to the sky view factor (SVF), the irradiance from the adjacent objects can also have a significant impact on the accuracy of the temperature retrieval, especially when the emissivity is relatively low. This retrieval method provides a solution for high-resolution near-surface UAV/airborne TIR data and a promising framework for enhancing the LST accuracy using multi-source geographic information assistance.

On-orbit rapid detection of land surface anomalies is important for ensuring ecological security and human safety. Land surface anomalies (e.g., fire, industrial heat source, and deforestation, etc.) are often accompanied by different degrees of thermal anomalies. Existing methods for detecting thermal anomalies have focused primarily on high-temperature anomalies, without available approach for detecting widespread low-temperature anomalies. Here, a Novel Method based on Constructed Reference land surface temperatures (LST) for on-orbit remote sensing detection of various Thermal Anomalies (NMCRTA) is proposed and further evaluated using Landsat 8 LST product. In this method, we first construct an fitted reference temperature based on LST spatiotemporal trend surface modeling and a real reference temperature based on contextual averaging. Then, the difference (including the step of removing atmospheric effects) between on-orbit observed LST and fitted reference LST, and the difference between on-orbit observed LST and real reference LST are calculated, respectively. Finally, these two temperature differences are utilized to detect thermal anomalies using corresponding thresholds. The results indicate that the NMCRTA can effectively detect deforestation, newly constructed buildings, and river drying, with an overall F1-score of 0.867, in a 100 × 100 km region scale. Meanwhile, the NMCRTA exhibited excellent accuracy in detecting fires, deforestation, and landslides at the 15 × 15 km scene scale, achieving F1-scores of 0.943, 0.857, and 0.791, respectively. Furthermore, the NMCRTA can continuously capture different thermal anomaly events associated with a newly constructed industrial heat source and perform well in nighttime. The NMCRTA is promising for future on-orbit remote sensing detection of various land surface anomalies, as a valuable supplement to optical on-orbit detection method.

As urbanization accelerates, more and taller buildings and less greenery are closely related to changes in the urban thermal environment (UTE). Knowledge of spatial and temporal variations of UTE is becoming increasingly concerning, and this can be measured with the land surface temperature (LST). Satellite observation of LST is an important tool for monitoring; the strong thermal anisotropy limits the use of satellite thermal infrared (TIR) data. Hitherto, the poor investigation was focused on the modeling and analysis of urban thermal anisotropy (UTA), especially in mountainous urban areas with multi-slope environments. These areas exhibit a distinctive “roof effect”, which is defined as the radiative transfer effect between the roof and the adjacent wall due to the slope that results in different heights between the roofs; multiple scattering has also been changed. Although an analytical thermal anisotropy model for the urban canopy over sloping terrain (AU3SM) has been proposed, its inability to effectively account for roof effects and multi-scattering mechanisms limits its daytime TIR observation applicability. To address these limitations, we developed an enhanced AU3SM that considers the roof effect and multiple scattering, which is labeled AU3SM-RS. The model was evaluated using measurements based on unmanned aerial vehicles (UAVs) in the mountainous city of Chongqing, China, with values of the root mean square error (RMSE) and coefficient of determination (R2) of 0.83 K and 0.93 in UTA. Comparison with a graphic processing unit-based solution for the faster 3-D radiative transfer model (GRay) further validates the model’s reliability with RMSE and R2values of 0.12 K and 0.96, respectively. Simulations in a certain scenario reveal that as the slope increases, the roof effect increases and the multiple scattering effect decreases in UTA and brightness temperature (BT), ignoring the roof effect and multiple scattering can result in maximum UTA biases of approximately 0.54, 0.48, and 0.72 K, BT biases approximately 1.02, 1.62, and 2.4 K at 5°, 15°, and 30° slopes, the biases due to neglecting the second scattering are very slight compared to the first scattering. Under certain conditions with a slope of 10°, the wider roof and narrower roadway, a related more dramatic roof effect; the narrower and deeper street canyon, a related more dramatic scattering effect. The proposed model is an efficient computational tool to assess UTA in mountainous areas quickly.

Land Surface Temperature (LST) estimation relies on precise measurements of surface thermal infrared (TIR) radiance. LST, classified as an essential climate variable (ECV), exhibits rapid temporal variations within specific spatial scales, closely tied to illumination and scan angles. To enhance the utility of LST products, this study proposes a novel approach for concurrently reconstructing the temporal profile and angular dependence. The suggested method, known as the visible-thermal envelope method, integrates kernel-driven (KD) and diurnal temperature cycle (DTC) models, addressing surface structure and thermal factors, respectively. Additionally, this method facilitates LST downscaling by leveraging the higher spatial resolution of visible and near-infrared (VNIR) data, assuming temperature differences are homogeneous within coarse pixels. To validate the reliability of the proposed approach, TIR data from the geostationary satellite Himawari 8 are amalgamated with VNIR data from the polar-orbit satellite Sentinel-3A/3B. When compared to field measurements, the reconstructed results exhibit improvements with a total bias of 0.55 K and Root Mean Square Error (RMSE) of 2.14 K. Notably, in contrast to the original uncorrected results, this correction results in an approximate 50% reduction in bias and a 10% decrease in RMSE.

As one of the most important spectral invariants, the recollision probability (p) is a crucial physical parameter that describes the interaction process between the photon and the canopy. This study discusses the applicability of the recollision probability in the thermal infrared (TIR) domain by Monte Carlo simulation. By extending the recollision probability from the VNIR to the TIR, an Improved Canopy Brightness Temperature model based on the p theory (ICBT-P) is introduced to simulate the canopy thermal radiation directionality. A cross-validation between the analytical ICBT-P model and 3D DART model was conducted, yielding R2values over 0.99 in two homogeneous scenarios where the leaf and soil components exhibit a significant temperature difference. In other two homogeneous scenarios where the leaf and soil components have the same temperature, the R2values are over 0.98. The average brightness temperature difference across the four typical scenes is 0.08K. This new model (i.e., ICBT-P) has the potential to be used for the canopy temperature estimation and the component temperature separation of leaves and soil.

Daily mean, maximum, and minimum land surface temperatures (Tmean, Tmax, and Tmin) are fundamental parameters for geophysical studies. Machine learning methods provide powerful tools to estimate these three parameters from instantaneous land surface temperature (LST) observations. This paper evaluated eight machine learning methods for two Moderate Resolution Imaging Spectroradiometer (MODIS) LST products (MxD11A1 and MxD21A1) and their joint application. The preliminary results indicate the support vector regression and linear regression methods exhibit satisfactory performance and fusion of two MODIS LST products is beneficial for estimating Tmeanand Tmax. Notably, the MxD11A1 nighttime LST product performs best for estimating Tmin.

Land surface temperature (LST) is an essential parameter for geoscience studies, and the 1-km polar-orbiting satellite LST products are widely used due to their global coverage on a daily basis. Aimed to jointly utilize all 1-km LST products to thoroughly quantify the LST temporal tendency within one day, comprehensive evaluation of existing 1-km LST products over the same sites is meaningful. In this study, taking sparsely vegetated sites as an example, ten 1-km polar-orbiting LST products were evaluated against the in situ measurements of ten USCRN sites with high spatial representativeness (from December 2019 to November 2020 covering four seasons). The evaluated LST products included three Moderate Resolution Imaging Spectroradiometer (MODIS) LST products (i.e., MOD11, MYD11, and MYD21), two Sea and LST Radiometer (SLSTR) LST products (i.e., Sentinel-3A and Sentinel-3B), one LST product from Advanced Very High Resolution Radiometer (AVHRR) onboard Metop-B, three Visible Infrared Imaging Radiometer Suite (VIIRS) LST products (i.e., VNP21, S-NPP EDR, and NOAA-20 EDR), and one LST product from Visible and Infra-Red Radiometer (VIRR) onboard FY-3B. The results indicate that: 1) the RMSE is much higher in the daytime (2.5–4.0 K) than in the nighttime (1.7–2.5 K); 2) the daytime MBE is negative (from −0.3 to −2.8 K) and this underestimate trend is significantly slighter for nighttime MBE (from −1.7 to 0.2 K); and 3) the daytime RMSE varies from 2.7 to 4.6 K in summer, 2.7–4.5 K in spring, 2.5–3.7 K in winter, and 2.0–3.4 K in autumn. Furthermore, we found that the severe daytime RMSE is related to the thermal anisotropy amplitude, which is around 4 K in summer, 2 K in spring and autumn, and 1 K in winter. Correcting the thermal anisotropy of daytime LST is an essential step that needs to be taken seriously before conducting high-quality joint applications, such as diurnal temperature cycle (DTC) modeling.

Land surface temperature (LST) is an important parameter in many research fields. As a temperature characteristic, thermal radiation has an obvious directionality. In urban areas, the directional anisotropy (DA) in surface temperatures can reach 10 K during airborne measurements, limiting the applicability of urban surface temperature products. Notably, the intricate 3-D structures and heterogeneous temperature distributions in urban environments significantly influence this anisotropy. In recent decades, numerous models have been developed for urban thermal anisotropy (UTA) analysis; however, most of these models predominantly focus on horizontal surfaces, with little consideration of mountainous architectures. To balance the calculation efficiency and model complexity, this study introduces an analytical thermal anisotropy model applicable to urban areas. This model considers the effects of 3-D structures and slopes and is labeled as the AU3SM. The multiangle data used for the analysis are recorded by an unmanned aerial vehicle (UAV) with a multicircle observation scheme in Chongqing, China. Comparisons of the AU3SM simulated data with these measurement data yield root mean square error (RMSE) and coefficient of determination ($R^{2}$) values of 0.57 K and 0.89, respectively. Furthermore, similar comparisons with the discrete anisotropic radiative transfer (DART) model yield$R^{2}$and RMSE values of 0.98 and 0.12 K, respectively. Simulations reveal that slope influences the UTA, and the hotspot lies on the opposite side of the solar direction; ignoring slope values of 5°, 10°, 15°, and 25° results in UTA maximum biases of approximately 1.2, 3.5, 7, and 7.5 K, respectively, under certain conditions. These findings demonstrate that the proposed model can accomplish rapid UTA assessments in mountainous areas.

Early warning systems and emergency management for disasters, environmental pollution, and illegal development require timely and accurate Earth surface anomaly detection (ESAD). Remote sensing, which uses satellites to observe the Earth’s surface, is an emerging approach to address this need. However, current remote sensing methods for ESAD are limited by their focus on specific anomalies, reliance on high-level satellite data, and the demand for significant computational and storage resources. In this article, we present a novel framework for general and on-orbit ESAD, which combines large vision models and lightweight priors. Our method characterizes images with a large vision model that is highly generalizable, reducing the dependency on high-level data, and compressing the prior base with sampling techniques, facilitating transmission and on-orbit storage. For on-orbit detection, we use a dictionary look-up style method for efficient anomaly prediction, enabling detection on satellites with limited computation resources. We evaluate our framework on typical scenarios and compare it with popular change detection (CD)-based and anomaly detection (AD)-based methods. Our results show that our framework achieves good performance while reducing the prior size by at least 95.56 times. Moreover, our framework can handle unpaired data, providing a chance to detect anomalies in the absence of near-term and paired images. Our framework has the potential to support the development and applications of general, on-orbit ESAD. The code and dataset are available at the following site:https://github.com/YummyWaffle/ESAD.

Land surface temperature (LST) is an important variable in Earth science research and can be measured using various thermal infrared (TIR) observations. Due to variations in data sources and inversion algorithms, LST products yield inconsistent results, further affecting subsequent applications (e.g., drought and vegetation monitoring). Although many evaluation studies have been conducted, most of them have focused on product validation or differences between inversion algorithms. There is a lack of analysis of the data sources, which is important for the data fusion. Therefore, a consistency analysis was conducted herein. Mainstream polar-orbiting satellite data were selected, including data from the Moderate Resolution Imaging Spectroradiometer (MODIS), Sea and Land Surface Temperature Radiometer (SLSTR), and Visible Infrared Imaging Radiometer Suite (VIIRS). The same inversion algorithm (split-window) for LST was employed across all datasets, thereby ensuring that differences in satellite data were the primary factor. Following validation based on in situ measurements, the polar-orbiting LST results were intercompared and analyzed with the LST results derived from the Himawari-8 Advanced Himawari Imager (AHI) observations. The results indicated that 1) similar conclusions were obtained from the intercomparison results and the ground-based validation results, with root mean square errors (RMSEs) for intercomparison results ranging from 2.903 K to 3.353 K; 2) based on the intercomparison results, regression analysis revealed that surface temperature status, land cover and vegetation information, and angular factors had a significant impact on the evaluation results, with t tests yielding p values less than 0.05 for all of these factors; and 3) based on a decision tree analysis, the contributions of angular factor, surface temperature status, and land surface structure were 50.9%, 31.3%, and 17.8%, respectively. These findings enhance knowledge regarding the impact of different satellite data on LST inversion results and emphasize the necessity of preprocessing before the joint application of satellite data, such as angle normalization and radiometric calibration.

There is a fundamental tension between metadata scalability and POSIX semantics within distributed file systems. The bottleneck lies in the coordination, mainly locking, used for ensuring strong metadata consistency, namely, atomicity and isolation. CFS is a scalable, fully POSIX-compliant distributed file system that eliminates the metadata management bottleneck via pruning the scope of critical sections for reduced locking overhead. First, CFS adopts a tiered metadata organization to scale file attributes and the remaining namespace hierarchies independently with appropriate partitioning and indexing methods, eliminating cross-shard distributed coordination. Second, it further scales up the single metadata shard performance by single-shard atomic primitives, shortening the metadata requests' lifespan and removing spurious conflicts. Third, CFS drops the metadata proxy layer but employs the light-weight, scalable client-side metadata resolving. CFS has been running in the production environment of Baidu AI Cloud for three years. Our evaluation with a 50-node cluster and microbenchmarks shows that CFS simultaneously improves the throughput of baselines like HopsFS and InfiniFS by 1.76--75.82× and 1.22--4.10×, and reduces their average latency by up to 91.71% and 54.54%, respectively. Under cases with higher contention and larger directories, CFS' throughput benefits expand by one order of magnitude. For three real-world workloads with data accesses, CFS introduces 1.62--2.55× end-to-end throughput speedups and 35.06--62.47% tail latency reductions over InfiniFS.

Land surface temperature (LST) is a vital parameter for studying global ecological, climatic, and environmental changes. A variety of regional and global scale LST products have been produced based on satellite remote sensing. Therefore, reliable retrieval accuracy is crucial for the application of LST products. A ground measurement processing method for medium and high spatial resolution LST products is proposed in this paper, to solve the spatial scale issues between the ground radiometer's field-of-view and the satellite pixel in T-based validation. The method is divided into two steps: LST quality control and spatial scale transformation. The Thermal Airborne Spectrographic Imager (TASI) data was used to evaluate the accuracy of the method. The results showed that this method can improve the reliability of the ground LST measurements.

Optical canopy models that connect land surface properties and satellite-observed radiance must be validated before being used. These models include the bidirectional reflectance distribution function (BRDF) models in the visible and near-infrared domains, and directional brightness temperature (DBT) models in the thermal infrared domain. Near-surface experiments have been extensively conducted to evaluate the modeling accuracy, including ground-, tower-, and aircraft-based measurements. Indeed, it should be noted that in situ measured BRDF/DBT results are sensitive to the experiment protocol, such as sensor moving orientation, flight height, and sampling frequency. A practical tool for optimizing the in situ measurement protocols is needed in the community of remote sensing modeling. For that, we devised a virtual experiment framework based on the discrete anisotropic radiative transfer (DART) 3-D radiative transfer model that is capable of simultaneously simulating both the BRDF/DBT pattern and the images acquired by in situ cameras. Here, as an optimization case, we use it to determine the optimal sensor flight orientation over heterogeneous vegetated canopies (a row-planted scene with three solar angles and a discrete scene with three solar angles) for measuring their DBT distribution. Results showed considerable errors (i.e., image-extracted DBT minus DART-simulated DBT) exist for sensor flight orientation along the canopy rows ($R^{2}$= 0.24 and root mean square error (RMSE) = 4.32 K), and they become much smaller ($R^{2}$= 0.94 ~ 0.98 and RMSE = 0.82 ~ 1.03 K) in other typical orientations (e.g., cross row plane, solar principal plane, and cross solar principal plane). The critical azimuth offset relative to the row direction that can ensure an acceptable RMSE < 1 K is quantified as atan(3*Unitwidth/Scenesize) based on a series of intensive simulations by this new tool. However, the RMSE of the discrete scene is not sensitive to the flight orientation. Such accuracy differences in various protocols were experimentally verified over row-planted maize using a 4-D tower in Huailai, Hebei, China. The result highlights the great potential of this newly designed DART-based virtual experiment to optimize near-surface experiment protocols.

All-weather Land Surface Temperature (LST) derived from passive microwave (PMW) sensors has significant implications for characterization on the physical processes of surface energy and water balance at local through global scales. However, the PMW sensors (e.g., the AMSR-E) suffer from swath gaps, cannot provide completely spatial-gapless observations. The existing Multi-temporal Feature Connection-CNN (MTFC-CNN) method caused obvious traces of ‘gaps’ when the sample number is small or features are not rich. This paper proposes a Sample Optimized-MTFC (SO-MTFC) seamless reconstruction method based on analyzing the periodicity and complementarity of AMSR-E swath gaps. Sample optimization includes two aspects: sample enhancement and single-cycle mask strategy. Taking China’s landmass as the study area, experimental results show that the original AMSR-E LSTs and the reconstructed AMSR-E LSTs are basically connected seamlessly. Validation against with the MODIS LSTs show that the daytime (nighttime) RMSEs of the original and the reconstructed AMSR-E LSTs are 3.87 K (2.57 K) and 4.76 K (2.96 K), respectively; while the corresponding daytime (nighttime) R2are 0.88 (0.94) and 0.73 (0.90), respectively. Validation against with the six in-situ LSTs show that the RMSEs and R2of reconstructed AMSR-E LSTs against in-situ LSTs are almost consistent with those of the original AMSR-E LST. The ablation study proves the effectiveness of the sample optimization. These findings indicated the SO-MTFC achieved a good reconstruction effect. Compared with the MTFC-CNN, the SO-MTFC got higher scores in visually and quantitatively. The SO-MTFC can potentially be implemented with other satellite PMW sensors to produce completely spatial-seamless PMW LST records on a global scale.

Land surface temperature (LST) is a vital parameter for studying global ecological, climatic, and environmental changes. Although various LST retrieval algorithms have been proposed, including split-window (SW), dual-window (DW), three-channel (TC), and dual-angle (DA) algorithms, few studies have compared these algorithms using the same satellite observations. The Sea and Land Surface Temperature Radiometer (SLSTR) onboard Sentinel-3A provides a unique opportunity to conduct this comparison owing to its dual-angle viewing capability and multiple thermal infrared (TIR) and mid-infrared (MIR) channels. Here, we implemented two SW algorithms, one DW algorithm, two TC algorithms and one DA algorithm for the SLSTR data. The LST retrievals from these six algorithms were validated, along with the SLSTR operational LST product based on an emissivity-implicit SW algorithm. Temperature-based and radiance-based validation methods were used to evaluate different LST retrievals across different land cover types. The results indicated that the proposed SW algorithm had the highest accuracy, followed by the Pérez-Planells SW and the official algorithms. The overall root-mean-square errors (RMSEs) of these three SW algorithms were 1.42 K, 1.79 K and 2.05 K, respectively. The three algorithms involving the MIR channel (one DW and two TC algorithms) were more suitable for nighttime LST retrieval and had similar performances to the three SW algorithms, with a nighttime RMSE of approximately 1.36 K. The LST retrieval accuracy of the DA algorithm had the highest uncertainty and was closely related to the angular variation in surface emissivity and brightness temperature. The findings of this study contribute to a better understanding of the different LST retrieval algorithms and facilitate potential improvements in the official LST retrieval algorithm for SLSTR.

Land surface temperature (LST) is an important parameter that characterizes the energy balance of the land surface, and it is widely used in various research fields. This paper proposes an operational split-window (SW) algorithm for use with the Chinese Fengyun-3 (FY-3) series satellite data, with the purpose of generating long-term global LST products. The algorithm primarily involves three steps. First, the brightness temperatures of the FY-3 Visible and Infra-Red Radiometer (VIRR) were recalibrated using historical recalibration coefficients to improve the accuracy of the absolute radiometric calibration. Second, daily dynamic emissivity maps were estimated using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global emissivity dataset (GED) and vegetation/snow cover products based on the vegetation cover method. Finally, the coefficients of the SW algorithm were simulated using MODTRAN 5 combined with the SeeBor V5.0 atmospheric profile library and ASTER spectral library, and then the coefficients were stratified by the view zenith angle and atmospheric water vapor content to improve the fitting accuracy. The proposed SW algorithm was integrated into the MUlti-source data SYnergized Quantitative (MUSYQ) remote sensing production system to then generate FY-3 VIRR LST products. Ten land surface sites from the HiWATER and SURFRAD networks and nine water surface sites from the National Data Buoy Center (NDBC) were used to evaluate the accuracy of the FY-3 VIRR LST products. The results demonstrated that the accuracy of the historical recalibration coefficients of the FY-3A/B VIRR is higher than that of the operational calibration coefficients for LST retrieval. The evaluation results revealed that the FY-3A VIRR LST products (2009-2013) had a bias of 0.13 K and an RMSE of 2.77 K, and the FY-3B VIRR LST products (2011-2020) had a bias of -0.07 K and an RMSE of 2.83 K. These results demonstrate that the proposed operational SW algorithm has reasonable accuracy and can be used to produce global LST products from the FY-3 VIRR data.

Three-dimensional (3D) radiative transfer (RT) models are frequently recognized as a prerequisite when using high spatial resolution remote sensing data in heterogeneous surfaces. However, most studies of 3D RT models have been restricted to limited applications due to the low computational efficiency. Therefore, this study proposed a graphic processing unit (GPU)-based solution for the ray tracing 3D RT model. A state-of-the-art graphics and compute application programming interface, Vulkan, was introduced to implement the RT process. A bounding box method was adopted for the computation acceleration. By comparison with a central processing unit (CPU)-based solution, the performance efficiency of the proposed solution is significantly better: the simulation time of a GPU model is significantly reduced by more than 99% when facing a large-scale simulation mission. The simulation accuracy of the two solutions is similar, with root mean squared errors (RMSEs) lower than 0.005, 0.032 and 0.31 K for the red, near-infrared (NIR) and brightness temperature images, respectively. An evaluation based on airborne multiangle measurements also indicated that the accuracy of the proposed solution was satisfactory for simulating the red and NIR bidirectional reflectance factor and brightness temperature directional anisotropies, with RMSEs lower than 0.003, 0.020 and 0.20 K, respectively, when treating the whole scene as a pixel. Considering the simulation accuracy and efficiency, a GPU-based model will be an important supplement to the CPU model.

There are five widely used kernel-driven models in the thermal infrared domain designed for the angular correction of land surface temperature (LST), including three-parameter Roujean Lagouarde (RL), Vinnikov, RossThick-LiSparseR (Ross–Li), LiStrahlerFriedl-LiDenseR (LSF-Li), and four-parameter Vinnikov-RoujeanLagouarde (Vinnikov-RL). Their fitting accuracies with hundreds of observation angles (i.e., sufficient angle) were studied; however, the fitting ability of these five models with limited observation angles is unknown, which makes it difficult to choose the appropriate one in applications. To solve this problem, 30 600 groups of multiangle directional brightness temperature (DBT) datasets were simulated by the unified optical-thermal 4-stream model considering scattering by arbitrary inclined leaves (4SAIL) model considering ten different leaf area index values, three leaf inclination distribution functions, two hotspot factors, 17 different component temperatures, five solar zenith angles, and six solar azimuth angles. Each group contains DBT values in 21 960 viewing directions [i.e., 61 viewing zenith angle (VZA)$\times $360 viewing azimuth angle (VAA)]. We assume that all limited observations are in the plane with VAA = 180°/0° and VZA changing from −60° to 60° with a step of 10°. There are 13 candidate angles to be selected. Five, seven, nine, and 11 angle sampling schemes include 225, 400, 225, and 36 limited multiangle combinations, respectively. Each combination was used to drive these five kernel-driven models to fit 21 960 DBTs for 30 600 groups of 4SAIL simulations. The root-mean-square error (RMSE) of each combination and mean RMSE of all 886 combinations were used to assess the overall fitting ability of five kernel-driven models. In addition, 1 k errors were added to the driven DBTs to evaluate the models’ robustness. Four groups of airborne measured DBTs were adopted to validate the assessment conclusions. Results show that the recommended order of these five models driven by 5–11 multiangle DBTs is Vinnikov-RL, LSF-Li, Vinnikov, Ross–Li, and RL when the driven DBTs do not contain errors; Vinnikov-RL, Vinnikov, LSF-Li, Ross–Li, and RL when the driven DBTs contain 1k errors; and Vinnikov-RL, LSF-Li, Ross–Li, RL, and Vinnikov for four groups of airborne measured datasets.

In the thermal infrared domain, there are four widely used 3-parameter kernel-driven models (RL, Vinnikov, Ross-Li and LSF-Li), and four newly proposed 4-parameter kernel-driven models (Vinnikov-RL, Vinnikov-Chen, LSF-RL and LSF-Chen). The fitting abilities of these eight models with limited observations are unknown which makes it difficult to choose the best one in practical applications. In this study, aimed to comprehensively evaluate all eight kernel-driven models, 1700 groups of multi-angle directional brightness temperature (DBT) datasets were simulated by the 4SAIL model considering 10 different leaf area indexes, 1 7 different component temperatures, 2 hotspot factors and 5 solar zenith angles. There are 13 angles from -60° to 60° with a step of 10° in the solar principal plane. Five, seven, nine and eleven angle groups produce 225, 400, 225, and 36 limited-angle combinations, respectively. The total 886 combinations were used to drive these eight kernel-driven models. Then, the fitted DBTs were compared with the 4SAIL simulated DBTs. Results show that LSF - Li always has the least fitting mean RMSE (0.24), followed by Ross-Li (0.33), Vinnikov (0.35), and RL (0.42) for 3-parameter models. The 4-parameter kernel-driven models have the same level accuracy (mean RMSE ≈ 0.10) and are much better than these 3-parameter models.

As a surface component, the tree trunk affects the top-of-canopy (TOC) emissivity and thermal infrared (TIR) radiance over a forest with fewer leaves, which is important for the inversion of land surface temperatures (LSTs) and further applications such as predicting forest fires and monitoring drought conditions. Therefore, the tree trunk effect was analyzed in this article using a thermal radiation directionality model, in which the forest structure was considered by the geometric optical (GO) theory and the spectral invariance theory was introduced into the GO framework for the single-scattering effect between components. The model used was evaluated using unmanned aerial vehicle (UAV)-based measurements with root-mean-square errors (RMSEs) lower than 0.25 °C for directional anisotropies (DAs) of brightness temperatures (BTs). Comparison with a 3-D radiative transfer model, discrete anisotropic radiative transfer (DART), also indicated an acceptable tool of the proposed model for the trunk effect with RMSEs lower than 0.003 °C and 1.2 °C for DAs of emissivity and BTs, respectively. In this study, the root-mean-squared difference (RMSD) levels between the vegetation-soil and vegetation-trunk-soil canopies, which were viewed as an equivalent indicator of the trunk effect, were provided for the TOC emissivity and BTs as well as their DAs, by combination with the changes in the leaf area index (LAI), stand density, trunk shape, and component temperatures, which can help identify the cases in which the trunk effect should be considered. According to a comprehensive analysis, for cases with sparse stand density ( α0.04), the tree trunk should be considered for a BT RMSD level lower than 0.5 °C when the LAI value was lower than 0.6. The corresponding LAI value was 0.8 for an RMSD level of BT DA lower than 0.3 °C. Moreover, for the cases with low soil emissivity, the difference in the TOC emissivity with and without trunk can reach up to 0.035, and the RMSD was still larger than 0.01 when the stand density and LAI were 0.05 and 0.6, respectively.

In this study, two collection 6 (C6) Moderate Resolution Imaging Spectroradiometer (MODIS) level-2 land surface temperature (LST) products (MYD11_L2 and MYD21_L2) from the Aqua satellite were evaluated using temperature-based (T-based) and radiance-based (R-based) validation methods over barren surfaces in Northwestern China. The ground measurements collected at four barren surface sites from June 2012 to September 2018 during the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment were used to perform the T-based evaluation. Ten sand dune sites were selected in six large deserts in Northwestern China to carry out an R-based validation from 2012 to 2018. The T-based validation results indicate that the C6 MYD21 LST product has a better accuracy than the C6 MYD11 product during both daytime and nighttime. The LST is underestimated by the C6 MYD11 products at the four T-based sites during the daytime, with a mean bias of -2.82 K and a mean RMSE of 3.82 K, whereas the MYD21 LST product has a mean bias and RMSE of -0.51 and 2.53 K, respectively. The LST is also underestimated at night by the C6 MYD11 products at the four T-based sites, with a mean bias of -1.40 K and a mean RMSE of 1.72 K, whereas the MYD21 LST product has a mean bias and RMSE of 0.23 and 1.01 K, respectively. For the R-based validation, the MYD11 results are associated with large negative biases during both daytime and nighttime at three sand dune sites and biases within 1 K at the other seven sites, whereas the MYD21 results are more consistent at all ten sand dune sites, with a mean bias of 0.45 and 0.70 K for daytime and nighttime, respectively. The emissivities for these two products in MODIS bands 31 and 32 were compared with each other and then compared with the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) emissivity and laboratory emissivity. The results indicate that the emissivities in MODIS bands 31 and 32 of MYD11 at the four T-based and three of the R-based validation sites are overestimated and result in LST underestimation, whereas the emissivities of MYD21 are more consistent with the laboratory emissivity. Besides, an experiment was carried out to demonstrate that the physically retrieved dynamic emissivity of the MYD21 product can be utilized to improve the accuracy of the split-window (SW) algorithm for barren surfaces, making it a valuable data source for retrieving LST from different remote sensing data.

Physical models simulating the radiative budget (RB) and remote sensing (RS) observation of three-dimensional (3D) landscapes are critical to better understand human and natural components of the Earth system and further develop RS technology. DART is one of the most comprehensive 3D models of Earth-atmosphere optical radiative transfer (RT), from ultraviolet (UV) to thermal infrared (TIR). It simulates the optical signal of proximal, aerial and satellite imaging spectrometers and laser scanners, the 3D RB and solar induced chlorophyll fluorescence (SIF) signal, for any urban or natural landscape and any experimental or instrument configuration. It is freely available for research and teaching activities (https://dart.omp.eu). Here, five recent advances are presented. 1) Atmosphere RT. 2) RT in non repetitive topography. 3) Monte Carlo modelling for fast RS image simulation of large landscapes. 4) SIF modelling for vegetation simulated as facets and turbid cells. 5) RS image inversion for mapping the optical properties of urban material and the urban radiative budget.

The in situ measurement of the leaf area index (LAI) from gap fraction is often affected by terrain slope. Path length correction (PLC) is commonly used to mitigate the topographic effect on the LAI measurements. However, the terrain-induced uncertainty and the accuracy improvement of the PLC for LAI measurements have not been systematically analyzed, hindering the establishment of an appropriate protocol for LAI measurements over mountainous regions. In this article, the above knowledge gap was filled using a computer simulation framework, which enables the estimated LAI before and after PLC to be benchmarked against the known and precise model truth. The simulation was achieved by using CANOPIX software and a dedicatedly designed ray-tracing method for continuous and discrete canopies, respectively. Simulations show that the slope distorts the angular pattern of the gap fraction, i.e., increasing the gap fraction in the down-slope direction and reducing it in the up-slope direction. The horizontally equivalent hemispheric gap fraction from the PLC can reconstruct the azimuthally symmetric angular pattern of the real horizontal surface. The azimuthally averaged gap fraction for sloping terrain can both be underestimated or overestimated depending on the LAI and can be successfully corrected through PLC. The topography-induced uncertainty in LAI measurements is found to be ~14.3% and >20% for continuous and discrete canopies, respectively. This uncertainty can be, respectively, reduced to ~1.8% and <; 7.3% after PLC, meeting the up-to-date uncertainty threshold of 15% established by the Global Climate Observing System (GCOS). Closer analysis shows that the topographic effect is influenced by fractional crown cover, and the largest uncertainty which corresponds to extensively clumping canopy can reach nearly up to 50%. The accuracy of the estimated LAI after PLC safely meets the GCOS uncertainty threshold even for this extreme case. This study demonstrates the necessity of a topographic correction for LAI measurements and the applicability of PLC for reconstructing the horizontally equivalent gap fraction and improving the LAI measurements over sloping terrains. The results of this article throw light on the design of a protocol for LAI measurements over mountainous regions.

Because of noises from land surface temperature (LST) inversion and different spatial resolution between nadir and oblique views, we proposed a combined algorithm for soil and vegetation temperatures with Sea and Land Surface Temperature Radiometer (SLSTR) data, in which a Bayesian strategy was adopted and retrieval results from a multi-pixel algorithm were selected as a priori information in a multi-angle algorithm. The SLSTR LST that inverted by a general split-window algorithm was evaluated first to understand inversion noise. Then, the combined algorithm was evaluated using a synthetic dataset, which displayed a robust performance than the multiangle and the multi-pixel algorithm.

There are four widely-used kernel driven models in the thermal infrared domain, including LSF-Li, RL, Ross-Li and Vinnikov. They have great potentials to achieve the angular effect correction for land surface temperature products. RL and Ross-Li are direct extensions of kernel models in the visible and near infrared region and other two models were directly designed for the thermal infrared region. Their fitting abilities should be evaluated comprehensively at first. Some evaluation works have been reported based on simulated datasets or measured dataset with limited view angles. In this paper, two measured multi-angle datasets with over 104directions were selected to do evaluation. Both of them were extracted from the airborne measurement using spatial-temporal average method. We found that LSF-Li is always the best one, followed by Ross-Li, RL and Vinnikov.

In this paper, in order to produce long term fully global land surface temperature (LST) product, the generalized split-window (GSW) algorithm and dual-window (DW) algorithm was used to retrieve LST from different geostationary (GEO) satellite data, including the FY-2E/4A, MTSAT-2/Himawari-8, MSG2, and GOES13/15. First, the coefficients of the GSW and DW algorithm were obtained from a simulation database constructed using the MODTRAN 5.2 and the SeeBor V5.0 atmospheric profile database. Second, the emissivity was estimated using the vegetation cover method, with the bare soil component emissivity derived from the ASTER global emissivity dataset (GED). Finally, the LST results of FY-4A AGRI and Himawari-8 AHI were retrieved and cross-validated. The results show that the LST algorithms developed in this work are capable of generating accurate high temporal resolution LST retrieval from global GEO satellite data.

In this study, we present an operational algorithm to retrieve the land surface temperature (LST) from MODIS thermal infrared data using physically retrieved emissivity product. This algorithm was implemented in the MUlti-source data SYnergized Quantitative (MuSyQ) remote sensing product system. First, the emissivity in the MODIS two split-window channels was estimated using the vegetation cover method, with the bare soil component emissivity derived from the ASTER global emissivity dataset. Then, the LST was retrieved using a modified generalized split-window algorithm. The MuSyQ MODIS LST product and the C6 MxD11 LST product were evaluated using ground measurements collected from four barren surface sites in Northwest China during the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment. The evaluation results indicate that the MuSyQ LST products provide better accuracy than the C6 MxD11 product during both daytime and nighttime at all four sites. This study demonstrates that physically retrieved emissivity products are a useful source for LST retrieval over barren surfaces.

The effects of Thermal Radiation Directionality (TRD) were originally evidenced through experiments in 1962, showing that two sensors simultaneously measuring temperature of the same scene may get significantly different values when the viewing geometry is different. Many models were developed to simulate the TRD effect aimed at mimicking the observed thermal infrared radiance with consideration of canopy structure, component emissivities and temperatures, and Earth surface energy exchange processes. In this paper, the models of vegetation canopies are classified into five groups. The basic assumption, and the disadvantages/advantages of them are summarized. Then, the recent progresses on the TRD modeling for vegetation canopy are given.

This letter presents an experiment to explore the nonisothermal effects on temperature and emissivity separation (TES). The innovation of this experiment lies in its design, which highlights the contrast between isothermal and nonisothermal conditions in emissivity measurements. We artificially created a sharply contrasting nonisothermal soil surface using liquid nitrogen cooling and solar heating. The iterative spectrally smooth TES (ISSTES) algorithm was used to process the experimental data. The analyzed results of the experimental data show that the nonisothermal conditions have a significant effect on TES. The bias of retrieved emissivity increases with the component temperature difference as well as with wavelength. The bias around the split window band can reach up to 1% when the difference of the component temperature is 40K. Considering that 1% error in emissivity can cause approximately 1K error of retrieved land surface temperature (LST), the nonisothermal effects on emissivity cannot be ignored. We hope that this experiment will arouse attention of the nonisothermal effects on TES and call for more efforts to be devoted to this issue in the future.

The directional anisotropy of canopy emissivity and brightness temperature in the thermal infrared band has widely been studied. However, the contribution of different scattering orders has been an open scientific question for many years. The recently proposed CE-P model enables us to analytically evaluate the different scattering orders. Herein, we derive expressions for the first double collisions (DCE12) and first triple collisions (DCE123). Our result shows that DCE123can simulate the observed emissivity with an error less than 0.001 and that DCE12is reasonably accurate when leaf emissivity is greater than 0.96. Numerical analysis shows that the contribution of quadruple or greater collisions can be ignored when the leaf (soil) emissivity is no less than 0.90. Furthermore, we develop the CBT-P model to simulate the directional brightness temperature (DBT) based on the new optimized CE-P model (DCE123) and validate it by 4SAIL (4-Stream Radiative Transfer Theory of Scattering by Arbitrary Inclined Leaves) and DART (Discrete Anisotropic Radiative Transfer) models. Both of isothermal (soil temperature is equal to leaf temperature) and nonisothermal (soil temperature is higher than leaf temperature) cases are considered. The maximum differences between the CBT-P model and 4SAIL (DART) are less than 0.35 K (0.42 K), the average differences between CBT-P and 4SAIL (DART) are less than 0.10 K (0.13 K), and the R2is over 0.99 (0.95) with component emissivities larger than 0.90 and the difference between soil and leaf temperatures less than 20 K. The directional anisotropy of DBT is dominated by the zero-scattering and the single scattering terms according to the new developed CBT-P model.

Many physical models have been proposed to simulate the directional anisotropy in the thermal infrared (TIR) region over vegetation canopies to produce angular corrected directional brightness temperature or land surface temperature. However, too many input parameters obstruct their operational use. Semiempirical kernel-driven models are designed to be a tradeoff between physical accuracy and operationality. Recently, four kernel-driven models have been proposed: the first two are direct extensions of kernel models in the visible- and near-infrared region and the last two were directly designed for the TIR region. In this paper, 153 continuous and 153 discrete canopies with varying structures and temperature distributions were considered in order to evaluate their accuracies against two physical models (4SAIL and DART). Their error distribution, scatterplots, and directional anisotropy patterns are compared. LSF-Li model, followed by Ross-Li, Vinnikov, and RL model, gave the best fitting results for all the scenes. The R2of all four kernel models can reach up to 0.82 for discrete scenes; however, the kernel-driven models underestimate the hotspot effect from continuous scenes; therefore, further improvements are necessary for operational use with future TIR satellite missions.

Land surface temperature and emissivity (LST&E) are key variables for a wide variety of surface-atmosphere studies, and it is very important to validate the accuracy of different LST&E retrieval algorithms. In this paper, the accuracy of the Advanced Spaceborne Thermal Emission and Reflection (ASTER) temperature emissivity separation (TES) algorithm with the water vapor scaling (WVS) method and the ASTER standard LSE&E products (without WVS) was validated using 12 ASTER scenes from May 2012 to September 2012 with concurrent ground LST and emissivity measurements collected in an arid area in northwest China during the Heihe Watershed Allied Telemetry Experimental Research experiment. Both the National Center for Environmental Prediction (NCEP) and MOD07 atmospheric profile products were employed to perform the atmospheric correction. The results showed that the WVSTES LSTs retrieved from the two profiles both demonstrate good accuracies, with average biases of 0.34 and 0.24 K and average root-mean-square errors (RMSEs) of 1.52 and 1.46 K for the NCEP and MOD07 profiles, respectively. When the WVS was not applied, we obtained an average bias of 1.26 K and an average RMSE of 2.05 K for the AST08 product. The emissivities from the WVSTES algorithm with the two profiles both showed good agreements with the ground CE312 measurements, with mean differences of the five ASTER bands that were less than 0.01. AST05 showed anomalous emissivity spectra and underestimated the emissivity values of graybody surfaces when compared with the ground data, with mean differences of greater than 0.015, which were more obvious for cases with high water vapor. This paper demonstrated that the WVS method is critical for retrieving ASTER TES LST with accuracies within 1.5 K and emissivity within 0.015 for a wide range of atmospheric conditions and land surface types.

In this study, to improve the accuracy of land surface temperature (LST) products over barren surfaces, we present an operational algorithm to retrieve the LST from Moderate-Resolution Imaging Spectroradiometer (MODIS) thermal infrared data using physically retrieved emissivity products. The LST algorithm involved two steps. First, the emissivity in the two MODIS split-window (SW) channels was estimated using the vegetation cover method, with the bare soil component emissivity derived from the ASTER global emissivity data set. Then, the LST was retrieved using a modified generalized SW algorithm. This algorithm was implemented in the MUlti-source data SYnergized Quantitative (MuSyQ) remote sensing product system. The MuSyQ MODIS LST product and the Collection 6 MODIS LST product (MxD11_L2) were compared and validated using ground measurements collected from four barren surface sites in Northwest China during the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment from June 2012 to December 2015. In total, 2268 and 2715 clear-sky samples were used in the validation for Terra and Aqua, respectively. The evaluation results indicate that the MuSyQ LST products provide better accuracy than the C6 MxD11 product during both daytime and nighttime at all four sites. For the daytime results, the LST is underestimated by the C6 MxD11 products at all four sites, with a mean bias of -1.78 and -2.86 K and a mean root-mean-square error (RMSE) of 3.16 and 3.94 K for Terra and Aqua, respectively, whereas the mean biases of the MuSyQ LST products are within 1 K, with a mean bias of -0.26 and -1.03 K and a mean RMSE of 2.45 and 2.71 K for Terra and Aqua, respectively. For the nighttime results, the LST is also underestimated by the C6 MxD11 products at all four sites, with a mean bias of -1.60 and -1.26 K and a mean RMSE of 1.93 and 1.60 K for Terra and Aqua, respectively, whereas the mean biases of the MuSyQ LST products are 0.16 and 0.58 K and the mean RMSEs are 1.12 and 1.25 K for Terra and Aqua, respectively. The results indicate that the underestimation of the C6 MxD11 LST product at all four sites mainly results from the overestimation of the emissivities in MODIS bands 31 and 32. This study demonstrates that physically retrieved emissivity products are a useful source for LST retrieval over barren surfaces and can be used to improve the accuracy of global LST products.

In this study, two latest Collection 6 (C6) MODIS level-2 LST products (MxD11_L2 and MxD21_L2) from both the Terra and Aqua satellites were evaluated against ground measurements collected in an arid area of northwest China during the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment. Ground measurements of four and a half years at four barren surface sites were used to carry out the evaluation, which took place from June 2012 to December 2016. The results show that the C6 MxD21 LST products demonstrate a better accuracy than C6 MxD11 LSTproducts both at daytime and nighttime for the four sites. For the daytime result, C6 MxD11 products underestimate the LST at the four barren surface sites, with an average bias of -1.77 K for Terra and -2.63 K for Aqua, while the average biases of MxD21 LST products are much smaller, with an average bias of 0.34 K for Terra and -0.67 K for Aqua, respectively. For the nighttime result, C6 MxD11 products also underestimate the LST at the four sites, with an average bias of -1.41 K for Terra and -1.01 K for Aqua, while the average biases of K LST products are 0.14 K for Terra and 0.54 K for Aqua, respectively.

In the paper, a geometric optical model is proposed for a vegetation-trunk-soil scene. The effect of tree trunks was analyzed by comparing directional brightness temperatures (BTs) between vegetation-soil and vegetation-trunk-soil scenes. The comparison result reveals the tree trunk can cause directional BTs as a whole lower because of its shadow and shaded area. Therefore, the tree trunk should be considered when retrieving temperatures from thermal infrared observations over a leafless forest. Efforts using measured TIR data requires to be done in the future.

Multi-angle information extraction is very important for the validation of land surface models, such as bi-directional reflectance distribution function (BRDF) model and directional brightness temperature (DBT) model. The experiment can be done in the near surface, on the plane, and on the satellite. The spatial-temporal average method is widely-used on the scale of plane. In this paper, we try to extend it to extract the multi-angle information from near surface observation sensors. We found this method can obtain good result over homogeneous land surface, but the observation protocol is critical over heterogeneous land surface. For instance, the observation plane should be perpendicular to the row direction. In addition, we found that the hot spot area will be shaded by the observation platform which leads to the BRDF result to be affected seriously.

A new directional canopy emissivity model (CE-P) based on spectral invariants can separate the multiple scattering effect and single scattering in vegetation canopy. So we can further evaluate the contribution of scattering effect to the canopy emissivity and brightness temperature based on CE-P model. Numerical analysis shows that the contribution of more than three times scattering can be ignored when the leaf (soil) emissivity is no less than 0.90. Then, we optimize CE-P model and obtain the expressions containing the first twice collisions (ε2) and first three times collisions (ε3). The result shows that ε3 can simulate the emissivity in any case with an error less than 0.001. Furthermore, we simulate the brightness temperature distribution using the optimized model and compare it with DART model. The difference between them is less than 0.3K and the R2of them is over 0.96 in all of the selected samples.

High resolution thermal infrared (TIR) remote sensing data are widely used in many fields, such as environmental protection, land resources, urban construction and other scientific applications. This article presents a high resolution thermal infrared satellite concept to meet above requirements. The satellite can operate in a sun synchronous orbit, and provide high spatial resolution data (~5m) with a 2-day revisit in 19 spectral bands. The payload of the satellite can obtain radiance data which are subsequently processed into measurements of atmospheric and surface properties, such as land surface temperature (LST), water surface temperature, evapotranspiration, column water vapor, material composition and others. For users with urgent needs, remote sensing data of the payload can be processed by utilizing the Intelligent Onboard Data Processing (IODP) technology, so that event information is obtained onboard. Through the direct data downlink, event information is transmitted to ground stations immediately, so application users can quickly respond. In addition, because of the small swath of the high resolution satellite, the payload must be installed on an agile platform, which can effectively enlarge swath by multi-stripes splicing imaging. In general, this article provides an overview of the technical approach and application direction of high resolution thermal infrared satellites, which can be used as a reference for satellite applications and engineering development.

In this paper, we proposed a temperature and emissivity separation (TES) algorithm for the simultaneous retrieval of land surface temperature and emissivity (LST&E) from the thermal infrared data of Chinese GaoFen-5 (GF-5) satellite's Multiple Spectral-Imager (MSI) payload. In order to improve the accuracy of the TES algorithm, a water vapor scaling (WVS) method for atmospheric correction was adopted. The Seebor V5.0 global atmospheric profile database and MODTRAN 5 were used to simulate the WVS coefficients. A total of 11 ASTER scenes were used to simulate the MSI images and concurrent ground measurements acquired in the HiW ATER experiment were used to validate the algorithm. The results showed that the bias and root mean square error (RMSE) in the retrieved LST were 0.47 K and 1.70 K, respectively, and the absolute emissivity differences between MSI and the ground measurements were smaller than 0.01 for the four MSI TIR bands, which demonstrated that the proposed algorithm can be used to retrieve high accurate and high spatial resolution LST&E from GF-5 MSI data.

A new directional canopy emissivity model (CE-P) based on spectral invariants is proposed in this paper. First, we prove the existence of the spectral invariant properties in the thermal infrared (TIR) band using a Monte Carlo model. Based on it, the equation of the new model is derived from the perspective of absorption. In this expression, single-scattering and multiscattering effects are separated analytically in the TIR band. We find that the overall contribution of multiple scatterings is less than 0.005 when the component emissivities are over 0.90, and the overall contribution decreases with increasing leaf or soil emissivity. Furthermore, the new model can avoid the logical difficulty encountered when using the traditional cavity effect factor to simulate the emissivity of a sparse vegetation canopy. The results of 4SAIL and Discrete Anisotropic Radiative Transfer (DART) are selected to do cross validation. The CE-P can achieve a high accuracy compared with 4SAIL and DART, with an absolute bias less than 0.002 when the leaf (soil) emissivity is equal to 0.98 (0.94). Four widely used analytical models are selected for comparison. The resulting accuracies of these models are ordered from CE-P to REN15, FR97, FR02, and VALOR96 with the most serious error up to 0.002, 0.002, 0.007, 0.013, and 0.014, respectively. Three main conclusions are obtained through the sensitivity analysis: the multiscattering between vegetation and the background can be ignored when the leaf (soil) emissivity is no less than 0.94 (0.90), the second and higher order scattering within the vegetation can also be ignored when the leaf (soil) emissivity is no less than 0.94 (0.90), and the single-scattering effect within the canopy should be considered which can be calculated using three view factors.

Soil semiempirical dielectric models (SEMs) are powerful, and they are generally considered a useful hybrid of both empirical and physical models. In this paper, the Wang-Schmugge dielectric model is improved to more accurately estimate the relative complex dielectric constants (CDCs) of moist soils. Instead of the Debye relaxation spectrum of liquid water located outside of the soil (i.e., free out-of-soil water) adopted in the Wang-Schmugge model, the Debye relaxation formula related to the free-water component inside the soil [i.e., free soil water (FSW)], which is correlated with the soil texture, is employed in the improved SEM. In addition, the effective conductivity loss term related to both soil texture and soil moisture is introduced to explain the ionic conductivity losses of FSW. Since the soil moisture influence is reduced at high frequencies, the effective conductivity loss term related to only the soil texture is also analyzed for 14-18 GHz. As in the Wang-Schmugge model, the relative CDC of bound soil water varies with the soil volumetric moisture content when the soil moisture is lower than the maximum bound water fraction in the new model, which takes a different approach than the Mironov mineralogy-based SEM. The proposed model obtains better fitting results than the three most widely employed SEMs. The improved model exhibits a significantly improved accuracy with a higher correlation coefficient (R2), a closer 1:1 relationship, and a lower root-mean-square error, including in the L-band, and especially in the imaginary part of the L-band.

Land surface temperature (LST) is often needed for using remotely sensed data to study the surface energy budget and hydrological cycle. However, LST is challenging to measure and simulate because of its high sensitivity to atmospheric instability and solar angle, particularly over large-scale heterogeneous scenes. We propose a model that combines radiosity theory and an energy budget method for surface temperatures; we also explore the anisotropic behavior of row-planted crop emissions. The surface thermodynamic equilibrium state is fulfilled via the interaction between the 3-D radiative transfer calculations of the thermal-region radiosity-graphics combined model and the energy balance equation. Despite its shortcomings, such as the time-consuming calculations, the proposed model is feasible according to the results of an intercomparison and validation analysis. The intercomparison shows that the model exhibits similar performance, in terms of surface temperature calculations, to that of the soil-canopy observation, photochemistry and energy balance model (root-mean-square differences) of 0.59 °C and 1.77 °C for the leaf and soil components, respectively. Excellent agreement with the observed directional variation over summer maize canopies is also obtained, with R2values exceeding 0.6 and a mean RMSE of 0.32 °C. Thus, we recommend the new combined model as an option for explaining directional anisotropy due to its potential application to 3-D scenes.

Surface upward longwave radiation (SULR) is a significant component of the surface radiation budget and is closely linked with evapotranspiration, soil moisture, and surface cooling on clear nights. Therefore, accurately estimating SULR is essential to better understand its spatiotemporal dynamics or to characterize the thermal environment of a given land surface. Currently, most methods for estimating SULR (including the physical and hybrid methods) fail to account for the thermal anisotropy, which can introduce significant errors into the calculation. We previously proposed the combined algorithm that considers the thermal anisotropy to more accurately estimate the SULR. However, this proposed method has several shortcomings. For example, it considers the directionality of the emissivity and the effective temperature separately under the support of a parametric directional emissivity model. However, the directional emissivity model is not maturely developed for different land surface types, especially on non-vegetated surfaces. And the separation of land surface temperature and emissivity may undermine the estimation accuracy. Furthermore, this proposed method requires a series of input parameters that is not always available, limiting its applicability. In this paper, we present a refined algorithm that uses a kernel-driven model and the technique of band conversion to calculate the SULR directly based on surface-leaving radiances. This direct physical algorithm is then applied to the Wide-angle infrared Dual-mode line/area Array Scanner data set and validated using longwave radiation data collected by automatic meteorological stations from the Heihe Watershed Allied Telemetry Experimental Research experiment. The results of these tests suggest that the direct algorithm works effectively. The root-mean-square error (RMSE) and mean bias error (MBE) of the direct algorithm on maize surfaces are 4.417 and 0.474 W · m-2, respectively. When the thermal anisotropy is incorporated, the RMSE and absolute MBE decrease by a maximum of 4.734 and 7.414 W·m-2, respectively. Different land types yield different results: for vegetable surfaces, the estimation biases of the direct model are approximately -2 W · m-2, whereas orchard surfaces yield biases are between -2 and -3.5 W · m-2, and village surfaces yield biases exceeding -10 W · m-2. These differences can be attributed to the varying effects of the kernel-driven model across different types of land surfaces. The RMSE and absolute MBE obtained using the direct algorithm are slightly smaller (0.587 and 1.685 W·m-2, respectively) than those obtained using the combined algorithm; they are also smaller than the results of the traditional temperature-emissivity algorithm (by 8.7 and 11.7 W · m-2, respectively).

In this study, a single-channel parametric model (SC-PM) algorithm were used to produce 300m LST product from HJ-1B IRS data. The NCEP atmospheric profiles and a parametric model were used for atmospheric correction. In order to improve the accuracy of the land surface emissivity (LSE), the 1km ASTER Global Emissivity Dataset (GED) and self-developed 5-day 1km vegetation cover product were used for estimating the LSE based on the Vegetation Cover Method. Two years of HJ-1B IRS LST product in Heihe River basin (Gansu province, China) from June 2012 to June 2014 were generated. The LST products were evaluated against ground observations collected during the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment. Four barren surface sites and ten vegetated sites were chosen for the evaluation. The results show that the produced HJ-1B IRS LST products demonstrate a good accuracy, with an average bias of 0.10 K and an average root mean square error (RMSE) of 2.43 K for all the sites during daytime. In addition, the biases are within 1K for the four barren surface sites. This indicate that using ASTER GED can produce reliable LST products from HJ-1B IRS data, especially for the barren surfaces.

Land surface component temperatures are important inputs in longwave radiation and evapotranspiration estimation models. Most component temperature inversion approaches focus only on two components, namely, soil and leaves, because space-based multiangle observations are lacking. This approach is inconsistent with ground-based measurements, which suggest that the temperatures of sunlit and shaded soil may significantly differ. This paper explores a three-component temperature inversion scheme that uses airborne multiangle thermal infrared observations to decrease the difference between the retrieved data and the actual subpixel temperature distribution. The FR97 model, which is an analytical directional brightness temperature model that was modified by dividing the soil component into sunlit and shaded portions, is adopted to calculate the matrix of component effective emissivity, which links multiangular observations and component temperatures. The new forward model and the inversion scheme are assessed using simulated data sets from the Scattering by Arbitrarily Inclined Leaves (4SAIL) model. The results indicate that the modified FR97 model provides good precision and that the inversion scheme based on the modified FR97 model is appropriate because of the model's simplicity and accuracy and the inversion's low sensitivity to noise. The inversion scheme is validated using airborne data collected by the wide-angle infrared dual-mode line/area array scanner over an area planted with maize and ground measurements collected during the Heihe Watershed Allied Telemetry Experimental Research campaign. The results indicate that the root mean square errors of the component temperatures of the leaves, sunlit soil, and shaded soil were 0.72 °C, 1.55 °C, and 2.73 °C, respectively. Because of the modified FR97's straightforward form and acceptable precision, we recommend this new retrieval scheme as an option for retrieving the component temperatures of leaves, sunlit soil, and shaded soil.

Surface upward longwave radiation (SULR) is an important component of the surface energy balance and is closely related to land surface temperature and emissivity. The estimation of SULR plays an important role in the study of surface energy circulation and climate change. State-of-the-art methods to estimate SULR, including the physical method and the hybrid method, are conducted without considering directional thermal radiation (DTR), which may induce a large error in the estimation, particularly over sparsely vegetated surfaces. In this paper, we modified the physical temperature–emissivity algorithm by combining a directional emissivity model (FRA97) and a kernel-driven DTR model to estimate the SULR of vegetated surfaces while considering the thermal directionality of the land surface. The most suitable kernel-driven model and an angle combination of the DTR were selected from six kernel-driven models and five angular combinations. The sensitivity of the proposed algorithm to the input parameters was also analyzed. The proposed algorithm was then validated with the Wide-angle infrared Dual-mode line/area Array Scanner (WiDAS) data set and longwave radiation data of automatic meteorological stations from the Heihe Watershed Allied Telemetry Experimental Research experiment. The results showed that the five-angle combination with large-angle intervals performs the best. When the leaf area index (LAI) is less than 1.2, the RossThick-LiSparseR model performs the best; when LAI is larger than 1.2, the RossThick-LiDenseR model is the most accurate. The SULR is not sensitive to surface downward longwave radiation and LAI, is slightly sensitive to leaf and soil emissivity at certain LAIs, and is highly sensitive to DTR, which may greatly affect the accuracy of the estimated SULR. The root-mean-square error (RMSE) and the mean bias error (MBE) of the SULR estimated using the WiDAS data and the proposed algorithm are 5.618 and −1.642 W/m2, respectively, thereby improving the estimation accuracy by as much as 7.479 and 10.511 W/m2at most in terms of RMSE and MBE, respectively, compared with the results calculated without considering the DTR.

The dielectric properties of vegetation reflect the coupling between the electromagnetic and physical properties of the material. Understanding dielectric behavior is essential for remote sensing applications. The traditional approach of measuring broadleaf permittivity involves exerting pressure on stacked leaves, which not only requires a complicated system to measure the pressure but also introduces a degree of uncertainty by changing the internal leaf microstructure. We propose a measurement technique that avoids the aforementioned shortcomings by using a vacuum packaging machine to remove the air from a sample bag, which is the first contribution of this paper. The proposed technique is validated against corn, cotton, and soybean leaf permittivity measurements. The second contribution of this paper involves how temperature is treated when modeling permittivity. The Debye-Cole dual-dispersion model, which is influential in vegetation remote sensing applications, is almost exclusively applied at 22 °C. However, it is assumed that this model can be extrapolated to other temperatures. Our analysis indicates that this model may overestimate permittivity at 28 °C, with mean relative errors (MREs) of up to 22% and 37% for relative permittivity and dielectric loss. Thus, we adopt the structure of the model due to its soundness and modify the model parameters. Overall, we used four different methods, among which the optimal method reduced the MRE by no more than -2.5%.

We compare five slope correction methods developed by Walter et al., Montes et al., Schleppi et al., España et al., and Gonsamo et al. (referred to as WAL, MON, SCH, ESP, and GON, respectively) using artificial fisheye pictures simulated by graphics software and a lookup table (LUT) retrieval method. The LUT is built by simulating the directional gap fraction as a function of leaf area index (LAI) and average leaf inclination angle (ALIA) using the Poisson law. LAI and ALIA estimates correspond to the case of the LUT that provides the lowest root-mean-square error between the observed gap fractions after slope correction and the simulated ones. Three LAI values (1.5, 3.5, and 5.5), four ALIA values (26.8°, 45°, 57.5°, and 63.2°), and three slope angles (0°, 20°, and 50°) constituted 36 samples of random scenes. ESP is recommended because its results are accurate and independent on the leaf angle distribution (LAD), while GON only performs well for spherical LAD. The three other methods present less good performances with underestimation or overestimation of LAI and/or ALIA depending on the LAD, and the recommended order for them is MON, SCH, and WAL.

A new geometric optical model is proposed in this letter to simulate the directional brightness temperature (DBT) distribution over mixed scenes of continuous crop and road. The DBT distributions of the crop and road zones are separately calculated, and the road zone consists of a road and adjacent crop sides. A road distribution polar map is designed to show all of the roads of different lengths, widths, and orientations in the scene. The airborne multiangle data set of the thermal infrared band that was acquired during the Heihe Watershed Allied Telemetry Experimental Research experiment is used for validation. The results demonstrate that the proposed model can simulate the DBT of a heterogeneous scene$(90\times 90\ \hbox{m}^{2})$with a root-mean-square error equal to 1.1 K and good trend similarity.

Leaf area index is usually estimated by gap fraction measurement. The mean projection of the unit leaf area on a plane normal to view direction is dependent with the view zenith angle except the view zenith angle is equal to 57.5°. Researchers tried to use this particular angle to overcome the influence of the leaf inclination distribution function. However, the ring of 57.5° is relative to the normal of slope on the slope ground while the measure platform is usually designed based on the flat coordinate. A general angle conversion strategy between slope coordinate and flat coordinate is required by the measurement on the slope. The conversion method was derived in this paper and the 57.5° rings of different slope angles were drawn and compared.

Vegetation clumping index is required for canopy reflectance modeling and leaf area index (LAI) retrieval. Frequently used is the average clumping index. The angle effect is normalized in the definition of average clumping index. Clumping index angular distribution with different canopy spatial pattern and the same architectural features of single plant is modeled and analyzed in this paper. The real structure of single plant and the canopy spatial pattern are simulated by modified L-systems. The gap fraction is calculated by Radiosity-Graphics combined Model (RGM) and the clumping index is calculated through a formula which is derived from Nilson's classic formula to solving gap fraction. It can be seen from the comparison of gap fraction of six scenes that the maximum is appeared in the direction of zenith and the shape between forward and backward is roughly symmetrical. It can be seen from the comparison of clumping index of six scenes that the maximum is not always appeared in the direction of zenith and the shape between forward and backward is not always roughly symmetrical. The results indicate that the asymmetry of canopy spatial pattern will bring the asymmetry of the clumping index angular distribution and the clumping index angular distribution depends on both the architectural features of single plant and the canopy spatial pattern.

We propose a new model to calculate the Directional Brightness Temperature (DBT) over the canopy of the row-planted corn. The most advantage of this new model is to simulate the DBT over the canopy cross the whole growth circle, both row structured and homogeneous structured. As the geometric structure of canopy changes greatly while the biomass grows from sparse to continuous, this new model is a better option as a tool to analysis remote sensing data.