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).

Northeastern China (NEC) is the largest grain base in China. Improving understanding of the effect of climate change on grain production over NEC is conducive to providing immediate response strategies for grain production. In this study, the relationships of the maize production with the dry state during the different maize growth stage have been investigated using the year-to-year increment method. Results showed that the severe drought that occurred from the jointing to maturity period have exerted severe effects on the maize growth. Further analysis indicated that the sea surface temperature (SST) anomalies over North Atlantic and Maritime Continent in later spring are the important factors affecting the summer droughts over NEC. The late spring SST anomaly over North Atlantic can excite the Rossby waves from the western North Atlantic and propagate eastward to NEC. The snow anomaly over western Siberia in late spring and the soil moisture anomaly over NEC in summer are key factors linking the SST anomaly to drought over the NEC. On the other hand, the Maritime Continent SST anomaly in late spring can modulate the activity of the East Asian jet stream via the East AsiaPacific (EAP) teleconnection, which can provide the favorable conditions for the soil moisture reduction over NEC. Eventually, a predictive model for maize yield over NEC is successfully developed by using the predictive indices of the North Atlantic and the Maritime Continental SST during late spring. Both the cross-validation and independent sample tests show that the calibrated prediction model is robust and exhibits high skill in predicting maize yield over NEC.

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.

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.

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

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 degrees C, 1.57 degrees C, and 0.38 degrees C for homogenous, row-planted vineyard, and sparse forest scenes, respectively; the analytical model appears similar performant (0.17 degrees C, 0.40 degrees C) for the homogeneous and sparse forest scenes, but less performant (2.39 degrees 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 degrees C and 0.19 degrees C, respectively, whereas it is less efficient for the row-planted scene with RMSE of 0.80 degrees 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.

Land surface temperature (LST) is an essential climate variable (ECV) which can be estimated from appropriate measurements of the surface thermal infrared (TIR) radiance. LST varies on a very short time scale and closely depends on the illumination and scan angles considered. To fully exploit LST products, a method for reconstructing the temporal profile and the angular dependence at the same time is proposed here. A combined visible- thermal envelope method (VT-KDTC) is built using kernel-driven (KD) and diurnal temperature cycle (DTC) models, referring to the surface structure and thermal factors, respectively. To demonstrate the reliability of the approach, TIR data from the geostationary satellite Himawari 8 are combined with visible and near-infrared (VNIR) data from the polar orbit satellite Sentinel-3A/3B. In addition to satellite observations, a synthetic dataset from the Soil Canopy Observation, Photochemistry and Energy Fluxes (SCOPE) model is also generated. Considering an anisotropy model in addition to the DTC model leads to a method displaying a better ability to simulate LSTs with a root mean squared error (RMSE) of 0.48 K against the original satellite results, compared to only the DTC model up to 1.44 K. By utilizing the field measurements as a reference, the reconstructed results are improved with a total bias of 0.72K and an RMSE of 2.58 K. Compared to the original results without correction, approximately 41% and 10% decreases are obtained in bias and RMSE, respectively. Our proposed method can also achieve LST downscaling supported by the higher spatial resolution of VNIR data when the temperature difference is assumed to be homogeneous within the coarse pixels. Thus, a simple achievable solution can be used for temperature reconstruction to enhance the quality of the LST product.