Land surface temperature (LST) plays an important role in Earth energy balance and water/carbon cycle processes and is recognized as an Essential Climate Variable (ECV) and an Essential Agricultural Variable (EAV). LST products that are issued from satellite observations mostly depict landscape-scale temperature due to their generally large footprint. This means that a pixel-based temperature integrates over various components, whereas temperature individual components are better suited for the purpose of evapotranspiration estimation, crop growth assessment, drought monitoring, etc. Thus, disentangling soil and vegetation temperatures is a real matter of concern. Moreover, most satellite-based LSTs are contaminated by directional effects due to the inherent anisotropy properties of most terrestrial targets. The characteristics of directional effects are closely linked to the properties of the target and controlled by the view and solar geometry. A singular angular signature is obtained in the hotspot geometry, i.e., when the sun, the satellite and the target are aligned. The hotspot phenomenon highlights the temperature differences between sunlit and shaded areas. However, due to the lack of adequate multi-angle observations and inaccurate portrayal or neglect of solar influence, the hotspot effect is often overlooked and has become a barrier for better inversion results at satellite scale. Therefore, hotspot effect needs to be better characterized, which here is achieved with a three-component model that distinguishes vegetation, sunlit and shaded soil temperature components and accounts for vegetation structure. Our work combines thermal infrared (TIR) observations from the Sea and Land Surface Temperature Radiometer (SLSTR) onboard the LEO (Low Earth Orbit) Sentinel-3, and two sensors onboard GEO (geostationary) satellites, i.e. the Advanced Himawari Imager (AHI) and Spinning Enhanced Visible and Infrared Imager (SEVIRI). Based on inversion with a Bayesian method and prior information associated with component temperature differences as constrained, the findings include: 1) Satellite observations throughout East Asia around noon indicate that for every 10 degrees change in angular distance from the sun, LST will on average vary by 0.6 K; 2) As a better constraint, the hotspot effect can benefit from multi-angle TIR observations to improve the retrieval of LST components, thereby reducing the root mean squared error (RMSE) from approximately 3.5 K, 5.8 K, and 4.1 K to 2.8 K, 3.5 K, and 3.1 K, at DM, EVO and KAL sites, respectively; 3) Based on a dataset simulated with a threedimensional radiative transfer model, a significant inversion error may result if the hotspot is ignored for an angular distance between the viewing and solar directions that is smaller than 30 degrees. Overall, considering the hotspot effect has the potential to reduce inversion noise and to separate the temperature difference between sunlit and shaded areas in a pixel, paving the way for producing stable temperature component products.

The Land Surface Temperature (LST) is well suited to monitor biosphere-atmosphere interactions in forests, as it depends on water availability and atmospheric/meteorological conditions above and below the canopy. Satellite-based LST has proven integral in observing evapotranspiration, estimating surface heat fluxes and characterising vegetation properties. Since the radiative regime of forests is complex, driven by canopy structure, components radiation properties and their arrangement, forest radiative temperatures are subject to strong angular effects. However, this depends on the scale of observation, where scattering mechanisms from canopy-to satellite-scales influence anisotropy with varying orders of magnitude. Given the heterogeneous and complex nature of forests, multi-angular data collection is particularly difficult, necessitating instrumentation distant enough from the canopy to obtain significant canopy brightness temperature and concurrent observations to exclude turbulence/atmospheric effects. Accordingly, current research and understanding on forest anisotropy at varying scales (from local validation level to satellite footprint) remain insufficient to provide practical solutions for addressing angular effects for upcoming thermal satellite sensors and associated validation schemes. This study presents a novel method founded in the optical remote sensing domain to explore the use of microcanopies that represent forests at different scales in the footprint of a multi-angular goniometer observing system. Both Geometric Optical (GO) and volumetric scattering dominated canopies are constructed to simulate impacts of anisotropy in heterogeneous and homogeneous canopies, and observed using a thermal infrared radiometer. Results show that heterogeneous canopies dominated by GO scattering are subject to much higher magnitudes of anisotropy, reaching maximum temperature differences of 3 degrees C off-nadir. Magnitudes of anisotropy are higher in sparse forests, where the gap fraction and crown arrangement (inducing sunlit/shaded portions of soil and vegetation) drive larger off-nadir differences. In dense forests, anisotropy is driven by viewing the maximum portion of sunlit vegetation (hotspot), where the soil is mostly obscured. Canopy structural metrics such as the fractional cover and gap fraction were found to have significant correlation with off-nadir differences. In more homogeneous canopies, anisotropy reaches a lower magnitude with temperature differences up to 1 degrees C, driven largely by volumetric scattering and components radiation properties. Optimal placement of instrumentation at the canopy-scale (more heterogeneous behaviour due to proximity to the canopy and small pixel size) used to validate satellite observations (more homogeneous behaviour due to larger pixel size) was found to be in cases of viewing maximum sunlit vegetation, for dense canopies. Given upcoming high spatial resolution sensors and associated validation schemes needed to benchmark LST and downstream products such as evapotranspiration, a better understanding of anisotropy over forests is critical to provide accurate, long-term and multi-sensor products.

Satellite-derived land surface temperature (LST) is a directional variable and has significant angular anisotropy. This characteristic contributes to enhance the differences among different satellite-derived LST products, and therefore increases the challenge of using multi-sensor and multi-decadal data to provide a long-term and angleconsistent LST climate data record. The kernel-driven model can balance the interpretability and operability well, so that it is suitable for angular normalization of LST products. The calibration of the kernel-driven model depends on multi-angle data which is difficult to obtain due to the spatial-temporal heterogeneity of LST. In this study, a novel LST angular normalization method based on the kernel-driven model was proposed to correct the angular effect of satellite-derived LST product by constructing multi-angle LST dataset from one geostationary satellite (GOES-R/ABI) and four polar-orbiting satellites (Terra/MODIS, Aqua/MODIS, Metop/AVHRR, and SNPP/VIIRS). The dataset gathered more abundant angle information, i.e., LSTs from three different observation geometries for the same pixel. The kernel-driven model was calibrated using the multi-angle LST dataset in the Continental United States (CONUS) during the year 2020. The discrepancies of the root mean square difference between LST before and after angular normalization range from 0.14 K to 1.10 K over nine land cover types in the four seasons. Similar results are obtained when the calibrated kernel-driven model was further expanded to other years and areas (i.e., the CONUS in 2021 and East Asia in 2020). The LST angular normalization method was applied to correct the angular effect of MODIS LST product. The results indicate that there is a strong correlation between the spatial distribution of LST differences (LST before and after angular normalization) and view zenith angle (VZA). MODIS LSTs before and after angular normalization were compared with Landsat 8 LST and Sentinel-3 A LST in near-nadir viewing for January, April, July, and October 2020. The angular normalization reduced the root mean square error (RMSE) between MODIS LST and Landsat 8 LST by 0.94-2.06 K in different months and by 0.13-2.61 K over various land cover types. For Sentinel-3 A, the RMSE decreased by 0.30-0.64 K in different months. The accuracies of MODIS LST before and after angular normalization were further validated using in situ measurements at the six SURFRAD sites. There are large discrepancies between the RMSE of MODIS LST before and after angular normalization versus in situ LST for VZA >= 45 degrees. The largest discrepancy is up to approximately 1.3 K at the GWN site. The LST angular normalization method has the potential to provide an angle-consistent LST climate data record.

Component temperature and emissivity are crucial for understanding plant physiology and urban thermal dynamics. However, existing thermal infrared unmixing methods face challenges in simultaneous retrieval and multicomponent analysis. We propose Thermal Remote sensing Unmixing for Subpixel Temperature and emissivity with the Discrete Anisotropic Radiative Transfer model (TRUST-DART), a gradient-based multi-pixel physical method that simultaneously separates component temperature and emissivity from non-isothermal mixed pixels over urban areas. TRUST-DART utilizes the DART model and requires inputs including at-surface radiance imagery, downwelling sky irradiance, a 3D mock-up with component classification, and standard DART parameters (e.g., spatial resolution and skylight ratio). This method produces maps of component emissivity and temperature. The accuracy of TRUST-DART is evaluated using both vegetation and urban scenes, employing Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images and DART-simulated pseudo-ASTER images. Results show a residual radiance error is approximately 0.05 W/(m2 & sdot;sr). In absence of the co-registration and sensor noise errors, the median residual error of emissivity is approximately 0.02, and the median residual error of temperature is within 1 K. This novel approach significantly advances our ability to analyze thermal properties of urban areas, offering potential breakthroughs in urban environmental monitoring and planning. The source code of TRUSTDART is distributed together with DART (https://dart.omp.eu).

Soil thermal conductivity (STC) plays a crucial role in regulating the energy distribution of both the surface and underground soil layers. It is widely applied in various fields, including engineering design, geothermal resource development and climate change research. A rapid and accurate estimation of STC remains a key focus in the study of soil thermodynamic parameters. However, the methods for estimating STC and their distinct characteristics have yet to be systematically reviewed. In this study, we used bibliometrics to comprehensively and systematically review the literature on STC, focusing on knowledge graph characteristics to analyze the development trend of calculation schemes. The main conclusions drawn from the study are as follows: (1) In recent years, most studies have been focused on soil thermal characteristics and their main contributing factors, the soil hydrothermal process in the Qinghai-Tibet Plateau, geothermal equipment and numerical simulations, and the exploration of geothermal resources. (2) A systematic review of various schemes indicates that no single scheme is universally applicable to all soil types. Moreover, a single parameterization scheme fails to meet the practical requirements of land surface process models. We evaluated the advantages and disadvantages of the traditional heat conduction schemes, parameterization schemes, and machine learning-based schemes and the findings suggest that a comprehensive scheme that integrates these three different schemes for STC simulations should be urgently developed.

The recent increase of the air temperature due to the global climate change is considered as one of the important reasons for the wildfires increase in the world, even in areas where the wildfires are not that common. In addition to the various physical damages adversely affecting the ecological balance, harmful gases and solid particles are released into the atmosphere due to wildfires, causing serious health problems. In this study, impacts of the most serious forest fire in modern history of the country lasting 16 days from 23rd of July 2022 in the National Park Bohemian Switzerland in the D & ecaron;& ccaron;& iacute;n district, Czech Republic, were investigated using remote sensing satellite datasets by cloud-based Google Earth Engine (GEE) platform. The normalized difference moisture index (NDMI), normalized burn ratio index (NBR), normalized difference vegetation index (NDVI), land surface temperature (LST) and soil moisture index (SMI) were calculated from Landsat-8 Operational Land Imager and Thermal Infrared Sensor (OLI and TIRS) dataset for the dates of 31st October 2021, 18th June 2022, and 31st October 2022. Relationship of the remote sensing indices were calculated to estimate the impacts of the wildfire. Furthermore, distribution of nitrogen dioxide (NO2) was extracted using Sentinel-5P TROPOMI (Tropospheric Monitoring Instrument) to observe changes before and after the forest fire in the study region. The burnt area approximately 13.20 km2 from the total area of 79.28 km2 was detected using different time series of the remote sensing indices in the national park.

Estimating the landscape and soil freeze-thaw (FT) dynamics in the Northern Hemisphere (NH) is crucial for understanding permafrost response to global warming and changes in regional and global carbon budgets. A new framework for surface FT-cycle retrievals using L-band microwave radiometry based on a deep convolutional autoencoder neural network is presented. This framework defines the landscape FT-cycle retrieval as a time-series anomaly detection problem, considering the frozen states as normal and the thawed states as anomalies. The autoencoder retrieves the FT-cycle probabilistically through supervised reconstruction of the brightness temperature (TB) time series using a contrastive loss function that minimizes (maximizes) the reconstruction error for the peak winter (summer). Using the data provided by the Soil Moisture Active Passive (SMAP) satellite, it is demonstrated that the framework learns to isolate the landscape FT states over different land surface types with varying complexities related to the radiometric characteristics of snow cover, lake-ice phenology, and vegetation canopy. The consistency of the retrievals is assessed over Alaska using in situ observations, demonstrating an 11% improvement in accuracy and reduced uncertainties compared to traditional methods that rely on thresholding the normalized polarization ratio (NPR).

Global warming accelerates permafrost degradation, compromising the reliability of critical infrastructure relied upon by over five million people daily. Additionally, permafrost thaw releases substantial methane emissions due to the thawing of swamps, further amplifying global warming and climate change and thus posing a significant threat to more than eight billion people worldwide. To mitigate this growing risk, policymakers and stakeholders need accurate predictions of permafrost thaw progression. Comprehensive physics-based permafrost models often require complex, location-specific fine-tuning, making them impractical for widespread use. Although simpler models with fewer input parameters offer convenience, they generally lack accuracy. Purely data-driven models also face limitations due to the spatial and temporal sparsity of observational data. This work develops a physics-informed machine learning framework to predict permafrost thaw rates. By integrating a physics-based model into machine learning, the framework significantly enhances the feature set, enabling models to train on higher-quality data. This approach improves permafrost thaw rate predictions, supporting more reliable decision-making for construction and infrastructure maintenance in permafrost-vulnerable regions, with a forecast horizon spanning several decades.

Permafrost, a major component of the cryosphere, is undergoing rapid degradation due to climate change, human activities, and other external disturbances, profoundly impacting ecosystems, hydroclimate, engineering geological stability, and infrastructure. In Northeast China, the thermal dynamics of Xing'an permafrost (XAP) are particularly complex, complicating the accurate assessment of its spatial extent. Many earlier mapping efforts, despite significant progress, fall short in accounting for some key local geo-environmental factors. Thus, this study introduces a new approach that incorporates four key driving factors-biotic, climatic, physiographic, and anthropogenic-by integrating multisource datasets and in situ observations. Four machine learning (ML) models [random forest (RF), support vector machine (SVM), logistic regression (LR), and extreme gradient boosting (XGB)] are applied to simulate permafrost distribution and probability, as well as to evaluate their performance. The results indicate that models' accuracy, ranked from highest to lowest, is as follows: RF (area under the curve (AUC) =0.88 and accuracy =0.81), XGB (0.86 and 0.77), LR (0.81 and 0.73), and SVM (0.76 and 0.66), with RF emerging as the most effective model for permafrost mapping in Northeast China. Analysis of the relationships between predictors and permafrost occurrence probability (POP) indicates that vegetation and snow cover exert nonlinear effects on permafrost, while human activities significantly reduce POP. Additionally, finer soil textures and higher soil organic matter content are positively correlated with increased POP. The modeling results, combined with field survey data, also show that permafrost is more prevalent in lowlands than in uplands, confirming the symbiotic relationship between permafrost and wetlands in Northeast China. This spatial variation is influenced by local microclimates, runoff patterns, and soil thermal properties. The primary sources of model error are uncertainties in the accuracy of multisource datasets at different scales and the reliability of observational data. Overall, ML models demonstrate great potential for mapping permafrost in Northeast China.

The European spruce bark beetle (Ips typographus) is an insect species that causes significant damage to Norway spruce (Picea abies) forests across Europe. Infestation by bark beetles can profoundly impact forest ecosystems, affecting their structure and composition and affecting the carbon cycle and biodiversity, including a decrease in net primary productivity (NPP), a key indicator of forest health. The primary objective of this study is to enhance our understanding of the interplay among NPP, bark beetle infestation, land surface temperature (LST), and soil moisture content as key components influencing the effects of climate change-related events (e.g., drought) during and after a drought event in the Bavarian Forest National Park in southeastern Germany. Earth observation data, specifically Landsat-8 TIR and Sentinel-2, were used to retrieve LST and leaf area index (LAI), respectively. Furthermore, for the first time, we incorporated a time series of high-resolution (20 m) LAI as a remote sensing biodiversity product into a process-based ecological model (LPJ-GUESS) to accurately generate high-resolution (20 m) NPP products. The study found a gradual decline in NPP values over time due to drought, increased LST, low precipitation, and a high rate of bark beetle infestation. We observed significantly lower LST in healthy Norway spruce stands compared to those infested by bark beetles. Likewise, low soil moisture content was associated with minimal NPP value. Our results suggest synergistic effects between bark beetle infestations and elevated LST, leading to amplified reductions in NPP value. This study highlights the critical role of integrating high-resolution remote sensing data with