Radiative Transfer Model Simulations Research Articles

A machine learning method trained by radiative transfer model inversion for generating seven global land and atmospheric estimates from VIIRS top-of-atmosphere observations

The Visible Infrared Imaging Radiometer Suite (VIIRS) has observed the Earth since 2011 and will continue for several decades. Unfortunately, few high-level land surface products have been produced and publicly released. The numerical inversion of radiative transfer model (RTM) has long been used for the retrieval of land surface and atmospheric variables from satellite data; however, it is computationally intensive. Some methods (e.g., using look-up tables) are efficient but are most suitable for estimating individual variable. Various machine learning (ML) models have been developed, mostly trained by either ground measurements that may not be adequate to represent different conditions, or RTM simulations that may not be realistic due to model uncertainty and unnecessary combinations of illumination-viewing geometries and variable values. In this study, we replaced the numerical inversion scheme by the ML algorithm that is trained using RTM inversion for estimating multiple variables simultaneously. Our training dataset for the ML models was generated by inverting the land and atmospheric variables from VIIRS top-of-atmosphere (TOA) reflectance data using a coupled land-surface–atmosphere RT model. Two multi-output ML algorithms were explored: backpropagation neural networks and random forest (RF) regression. The best-performing RF model was then applied to estimate seven land and atmospheric variables globally from VIIRS TOA data: the leaf area index (LAI), incident photosynthetically active radiation (PAR), fraction of absorbed photosynthetically active radiation (FAPAR), incident shortwave radiation (ISR), land surface albedo, land surface reflectance, and TOA albedo. The inversion results were validated using ground measurements at 54 sites and showed comparable accuracy to the numerical inversion scheme that is suitable for only small regions. An experiment was also conducted to test this method on a much larger region for detailed evaluation. Finally, seven global estimates in 2013 were experimentally produced. Inter-comparison with the existing satellite products showed that the VIIRS estimates had high consistency with the LAI and FAPAR products of the Global LAnd Surface Satellite (GLASS) products suite. Moreover, the VIIRS officially released the VNP43MA3 albedo and VNP09A1 surface reflectance products, the retrieved estimates also matched them well. It was also shown that simultaneous estimation of all seven variables was more accurate and efficient than estimating individual variable sequentially.

Can Forest Fires Be an Important Factor in the Reduction in Solar Power Production in India?

The wildfires over the central Indian Himalayan region have attracted the significant attention of environmental scientists. Despite their major and disastrous effects on the environment and air quality, studies on the forest fires’ impacts from a renewable energy point of view are lacking for this region. Therefore, for the first time, we examine the impact of massive forest fires on the reduction in solar energy production over the Indian subcontinent via remote sensing techniques. For this purpose, we used data from the Moderate Resolution Imaging Spectroradiometer (MODIS), the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIPSO), the Satellite Application Facility on support to Nowcasting/Very Short-Range Forecasting Meteosat Second Generation (SAFNWC/MSG) in conjunction with radiative transfer model (RTM) simulation, in addition to 1-day aerosol forecasts from the Copernicus Atmosphere Monitoring Service (CAMS). The energy production during the first quarter of 2021 was found to reach 650 kWh/m2 and the revenue generated was about INR (Indian rupee) 79.5 million. During the study period, the total attenuation due to aerosols and clouds was estimated to be 116 and 63 kWh/m2 for global and beam horizontal irradiance (GHI and BHI), respectively. The financial loss due to the presence of aerosols was found to be INR 8 million, with the corresponding loss due to clouds reaching INR 14 million for the total Indian solar plant’s capacity potential (40 GW). This analysis of daily energy and financial losses can help the grid operators in planning and scheduling power generation and supply during the period of fires. The findings of the present study will drastically increase the awareness among the decision makers in India about the indirect effects of forest fires on renewable energy production, and help promote the reduction in carbon emissions and greenhouse gases in the air, along with the increase in mitigation processes and policies.

On the Use of COSMO-SkyMed X-Band SAR for Estimating Snow Water Equivalent in Alpine Areas: A Retrieval Approach Based on Machine Learning and Snow Models

This study aims at estimating the dry snow water equivalent (SWE) by using X-band SAR data from the COSMO-SkyMed (CSK) satellite constellation. Time series of CSK acquisitions have been collected during the dry snow period in the Alto Adige test site, in the Italian Alps, during the winter seasons from 2013 to 2015 and from 2019 to 2021. The SAR data have been analyzed and compared with the in-situ measurements to understand the X-band SAR sensitivity to SWE, which has been further assessed by Dense Media Radiative Transfer (DMRT) model simulations. The sensitivity analysis provided the basis for addressing the SWE retrieval from the CSK data, by exploiting two different machine learning (ML) techniques, namely Artificial Neural Networks (ANN) and Support Vector Regression (SVR). To ensure a statistical independence of training and validation processes, the algorithms are trained and tested using SWE predictions of the fully distributed snow model AMUNDSEN as reference data and are subsequently validated on the experimental dataset. Due to its influence on the CSK estimates, the effect of forest canopy was accounted for in the analysis. Depending on the algorithm, the validation resulted in a correlation coefficient 0.78 ≤ R ≤ 0.91, and a Root Mean Square Error 55.5 mm ≤ RMSE ≤ 87.4 mm between estimated and in-situ SWE. Further analysis and validation are needed; however, the obtained results seem suggesting the Cosmo-SkyMed constellation as effective tool for the retrieval of the dry snow water equivalent in alpine areas.

Assessing Interactions Between Crop Biophysical Parameters and X-Band Backscattering Using Empirical Data and Model Sensitivity Analysis

Active microwave remote sensing data at different frequencies can provide crucial information on crop morphology and conditions, thus effectively supporting agronomic management at different scales. Despite the ever-increasing availability of spaceborne platforms and the extensive research developed throughout more than two decades, some knowledge gaps still await to be filled toward operational use, dealing with SAR backscatter response to crop-specific features and seasonal dynamics, including the effects of agronomic practices. In this work, we used variance-based global sensitivity analysis (GSA) as a quantitative framework for investigating the sensitivity of X-band backscattering to agronomic and morphological features typical of two different crops maize and rice. To this end, we jointly exploited empirical data on crop status and growth, high-resolution TerraSAR-X (TSX) data, and microwave radiative transfer model (RTM) simulations. Phenology-informed simulations allowed us to quantify the contributions of different scattering mechanisms for the two crops under varying observation setups, to assess the sensitivity of X-band backscattering to morphostructural crop biophysical parameters (BPs) (and their interactions), and to evaluate the effects of crop biomass on backscatter across growth stages. In particular, multidimensional GSA outputs accounting for model input correlations through Shapley effects provided a comprehensive suite of information on the relative proportion of total backscatter variance explained by a range of parameters and a quantitative description of the different behavior in vertical and horizontal polarization (changing throughout plant growth) in paddy rice, and the mixed contribution of canopy density and leaf angle distribution (depending on the incident angle) in maize.

L-Band Microwave Satellite Data and Model Simulations Over the Dry Chaco to Estimate Soil Moisture, Soil Temperature, Vegetation, and Soil Salinity

The Dry Chaco in South America is a semi-arid ecoregion prone to dryland salinization. In this region, we investigated coarse-scale surface soil moisture (<i>SM</i>), soil temperature, soil salinity and vegetation, using L-band microwave brightness temperature (<i>T</i>_B) observations and retrievals from the Soil Moisture Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite missions, Catchment Land Surface Model (CLSM) simulations, and in situ measurements within 26 sampled satellite pixels. Across these 26 sampled pixels, the satellite-based <i>SM</i> outperformed CLSM <i>SM</i> when evaluated against field data, and forward L-band <i>T</i>_B simulations derived from in situ <i>SM</i> and temperature performed better than those derived from CLSM estimates when evaluated against SMOS <i>T</i>_B observations. The surface salinity for the sampled pixels was on average only 4 mg&#x002F;g and only locally influenced the <i>T</i>_B simulations, when including salinity in the dielectric mixing model of the forward radiative transfer model (RTM) simulations. To explore the potential of retrieving salinity together with other RTM parameters to optimize <i>T</i>_B simulations over the entire Dry Chaco, the RTM was inverted using 10 years of multi-angular SMOS <i>T</i>_B data and constraints of CLSM <i>SM</i> and temperature. However, the latter modeled <i>SM</i> was not sufficiently accurate and factors such as open surface water were missing in the background constraints, so that the salinity retrievals effectively represented a bulk correction of the dielectric constant, rather than salinity per se. However, the retrieval of vegetation, scattering albedo and surface roughness resulted in realistic values.

Radiative Transfer Model Simulations for Ground-Based Microwave Radiometers in North China

Ground-based microwave radiometer profilers (MWRPs) are widely used to provide high-temporal resolution atmospheric temperature and humidity profiles. The quality of the observed brightness temperature (TB) from MWRPs is key for retrieving accurate atmospheric profiles. In this study, TB simulations derived from a radiative transfer model (RTM) were used to assess the quality of TB observations. Two types of atmospheric profile data (conventional radiosonde and ERA5 reanalysis) were combined with the RTM to obtain TB simulations, then compared with corresponding observations from three MWRPs located in different places in North China to investigate the influence of input atmospheric profiles on TB simulations and evaluate the quality of TB observations from the three MWRPs. The comparisons of the matching samples under clear-sky conditions showed that TB simulations derived from both radiosonde and ERA5 profiles were very close to the TB observations from most of the MWRP channels; however, the correlation was lower and the bias was obvious at 51.26 GHz and 52.28 GHz, which indicates that the oxygen absorption component in the RTM needs to be improved for lower-frequency temperature channels. The difference in location of the radiosonde and MWRP sites affected the TB simulations for the water vapor channels, but had little impact on temperature channels that are insensitive to humidity. Comparisons of both simulations (ERA5 and Radiosonde) and the corresponding TB observations from the three sites indicated that the water vapor channels observation quality for the MWRP located in southern Beijing needs improvement. For the two types of profile data, ERA5 profiles have a more positive effect on TB simulations in the water vapor channels, such as enhanced consistence, reduced bias and standard deviation between simulations and observations for those MWRPs located away from the radiosonde station. Therefore, hourly ERA5 data are an optimal option in terms of compensating for limited radiosonde measurements and enhancing the monitoring quality of MWRP observations within 24 h.

The ultraviolet radiation environment and shielding in pit craters and cave skylights on Mars

Caves and their entrances have been proposed as habitable environments and regions that could have preserved evidence of life, mostly due to their natural shielding from the damaging ionizing and non-ionizing radiation present on the surface. However, no studies to date have quantitatively determined the shielding offered by these voids on Mars. This paper describes the ultraviolet (UV) radiation environment in such environments by means of extensive radiative transfer (RT) model simulations applied to representative void geometries on Mars. The spatial and temporal characterization of the radiation field inside each void geometry is also performed. The RT model was previously tested against in situ Mars Science Laboratory UV observations on the Martian surface; the model showed excellent performance in terms of predicting the UV irradiances. The results presented here strongly suggest that pit craters and cave skylights are effectively shielded from the damaging UV radiation found on the Martian surface. In particular, locations on the floors of pit craters and skylights without cavernous spaces strongly attenuate UV irradiance at latitudes higher than ~45°; observing in several cases irradiances of the order of ~2% of those values found on the surface. Numerical simulations of cave entrances show a reduction even more than two orders of magnitude in UV radiation, both in the maximum instantaneous and cumulative doses, throughout the year and at any location of the planet. The attenuated UV irradiance in cavernous spaces remains stable in a 10 Myr timescale, while at the same time the photosynthetically active radiation is higher than the minimum required for Earth-like phototrophs. The intermediate radiation environment between the damaging radiation on the surface and the permanent darkness of a hypothetical cave offered by voids on Mars may represent favorable environments for habitability without constraining the type of energy source for potential as-yet unknown Martian organisms.

Impact of Aerosol and Cloud on the Solar Energy Potential over the Central Gangetic Himalayan Region

We examine the impact of atmospheric aerosols and clouds on the surface solar radiation and solar energy at Nainital, a high-altitude remote location in the central Gangetic Himalayan region (CGHR). For this purpose, we exploited the synergy of remote-sensed data in terms of ground-based AERONET Sun Photometer and satellite observations from the MODerate Resolution Imaging Spectroradiometer (MODIS) and the Meteosat Second Generation (MSG), with radiative transfer model (RTM) simulations and 1 day forecasts from the Copernicus Atmosphere Monitoring Service (CAMS). Clouds and aerosols are one of the most common sources of solar irradiance attenuation and hence causing performance issues in the photovoltaic (PV) and concentrated solar power (CSP) plant installations. The outputs of RTM results presented with high accuracy under clear, cloudy sky and dust conditions for global horizontal (GHI) and beam horizontal irradiance (BHI). On an annual basis the total aerosol attenuation was found to be up to 105 kWh m−2 for the GHI and 266 kWh m−2 for BHI, respectively, while the cloud effect is much stronger with an attenuation of 245 and 271 kWh m−2 on GHI and BHI. The results of this study will support the Indian solar energy producers and electricity handling entities in order to quantify the energy and financial losses due to cloud and aerosol presence.

Revisiting the use of red and near-infrared reflectances in vegetation studies and numerical climate models

Surface reflectance data acquired in red and near-infrared spectra by remote sensing sensors are traditionally applied to construct various vegetation indices (VIs), which are related to vegetation biophysical parameters. Most VIs use pre-defined weights (usually equal to 1) for the red and NIR reflectance values, therefore constraining particular weights for red and NIR during the VI design phase, and potentially limiting capabilities of the VI to explain an independent variable. In this paper, we propose an approach to estimate biophysical variables, such as Leaf Area Index (LAI), Canopy Chlorophyll Content (CCC) and Fraction of Photosynthetically Active Radiation (FPAR) absorbed by green vegetation, represented as linear combinations of the red and NIR reflectances with weights determined empirically from observations and radiative transfer model (PROSAIL) simulations. The proof of concept is first tested on available close-range observations over maize and soybean crops in Nebraska, USA. The empirical results compare well with those from PROSAIL model simulations. The proposed LAI model is then used with data from Landsat 8, Sentinel-2 and Planet/Dove, and the results are validated with in situ LAI measurements in Ukraine. We show that the weights on red and NIR reflectances are vegetation-specific and stable in time. The approach is further tested on crops and forests in the conterminous USA and on a global scale using MODIS LAI and FPAR products as proxies for “ground observations”. These LAI and FPAR, however, are not independently measured but derived from the corresponding remotely sensed reflectances, which precludes recommending a final set of the weights/coefficients for the users, and, thus, should be considered mostly for demonstrating the concept. The results for crop types, other than maize and soybean, and for all forests are conceptual and need to be tested with real ground data. It was, however, encouraging to see that the derived maps of coefficients/weights exhibit regular patterns over the globe compatible with those of vegetation classes and crop types. Tedious and thorough work on compiling available in situ measurements on various crops and forests needs to be accomplished prior to large-scale applications, and the method needs to be further tested and proven that it works at a large scale.The proposed parameterization may be attractive for global studies of various sub-classes of vegetation, once the parameter coefficients are established, validated, tabulated and their stability verified. Ultimately, this approach may provide quantification of vegetation traits for the past decades and be a useful asset for climate models that include satellite-derived land cover classifications and vegetation variables for simulating surface fluxes.This is a conceptual paper, with a proof-of-concept supported by observations over two crops, for which we had close-range observations. It is not a technical note, which would provide users with a recommended set of coefficients for global applications. Our intent was to develop a paradigm, which could ultimately be useful in global models.

Mapping of thermal power plant emitted atmospheric carbon dioxide concentration using AVIRIS-NG data and atmospheric radiative transfer model simulations

Our study demonstrates the potential of Airborne Visible/Infrared Imaging Spectrometer Next Generation (AVIRIS-NG) hyperspectral data and radiative transfer modeling to retrieve the atmospheric carbon dioxide (CO2) concentration from point source emission of one of India’s major coal-fed super thermal power plants at Kota, Rajasthan, India. AVIRIS-NG data with synchronous in situ measurements were collected as a part of the first ISRO-NASA joint airborne hyperspectral science campaign in India. The method utilized in our study includes theoretical simulations of at-sensor radiance in AVIRIS-NG spectral bands through the atmospheric radiative transfer model MODTRAN. Simulations were performed for the specific scene-sensor-atmospheric conditions pertaining to AVIRIS-NG overpass over Kota site. To eliminate the interfering effect of atmospheric water vapor with CO2 concentrations, simulations pertaining to variable water vapor values varying from 0.25 to 4.5 g / cm2 were also carried out, which in turn were utilized to normalize the water vapor effects. Based upon the simulation results pertaining to the specific absorption bands of CO2 in the shortwave infrared spectral range of AVIRIS-NG radiance data, a theoretical model was established between continuum interpolated band ratio (CIBR) and CO2 concentration. The CIBR − CO2 model coefficients for each of the water vapor subranges were then applied to AVIRIS-NG L1 radiance data to obtain the CO2 concentration map of the study area. A distinct CO2 plume could be detected from the coal-fed Kota Super Thermal Power Station with the CO2 concentration (XCO2) of the order of more than 500 ppmv near the power plant stacks. A range of XCO2 from 400 to 550 ppmv was observed within the scene. Our study has provided promising results in mapping the atmospheric CO2 content from the point source using the high-resolution airborne hyperspectral AVIRIS-NG data.

Validation of baseline and modified Sentinel-2 Level 2 Prototype Processor leaf area index retrievals over the United States

The Sentinel-2 Level 2 Prototype Processor (SL2P) is made available to users for the retrieval of vegetation biophysical variables including leaf area index (LAI) from Multispectral Instrument (MSI) data within the Sentinel Application Platform (SNAP). A limited number of validation exercises have indicated SL2P LAI retrievals frequently meet user requirements over agricultural environments, but perform comparatively poorly over heterogeneous canopies such as forests. Recently, a modified version of SL2P was developed, using the directional area scattering factor (DASF) to constrain retrievals as an alternative to regularisation (SL2P-D). Whilst SL2P makes use of prior information on expected canopy conditions, SL2P-D is trained using uniform distributions of input parameters to define radiative transfer model (RTM) simulations. Using in situ measurements available through the Copernicus Ground Based Observations for Validation (GBOV) service, we performed an extensive validation of SL2P and SL2P-D LAI retrievals over 19 sites throughout the United States. For effective LAI (LAIe), SL2P demonstrated good overall performance (RMSD = 0.50, NRMSD = 31%, bias = −0.10), with all LAI retrievals meeting the Sentinels for Science (SEN4SCI) uncertainty requirements over homogeneous canopies (cultivated crops, grasslands, pasture/hay and shrub/scrub), whilst underestimation occurred over heterogeneous canopies (deciduous forest, evergreen forest, mixed forest, and woody wetlands). SL2P-D retrievals demonstrated reduced bias, slightly improving overall performance when compared with SL2P (RMSD = 0.48, NRMSD = 30%, bias = −0.05), indicating its retrieval approach appears to offer some advantages over regularisation using prior information, especially at LAIe > 3. Additionally, SL2P-D resulted in 32% more valid retrievals than SL2P, with the largest differences observed at LAIe < 1. Validation against in situ measurements of LAI as opposed to LAIe yielded similar patterns but poorer performance (RMSD = 1.08 to 1.13, NRMSD = 49% to 52%, bias = −0.64 to −0.68) because the RTM used by SL2P and SL2P-D does not account for foliage clumping. In addition to the retrievals themselves, we examined the relationship between predicted uncertainties and observed differences in retrieved and in situ LAI. With respect to LAIe, SL2P’s predicted uncertainties were conservative, underestimating observed differences in only 35% of cases, whilst those for LAI were unbiased.

Comparison of Absorbed and Intercepted Fractions of PAR for Individual Trees Based on Radiative Transfer Model Simulations

The fraction of absorbed photosynthetically active radiation (fAPAR) is a key parameter for estimating the gross primary production (GPP) of trees. For continuous, dense forest canopies, fAPAR, is often equated with the intercepted fraction, fIPAR. This assumption is not valid for individual trees in urban environments or parkland settings where the canopy is sparse and there are well-defined tree crown boundaries. Here, the distinction between fAPAR and fIPAR can be strongly influenced by the background and large illumination variations due to multi-scattering and shadows of buildings. This study investigates the radiative budget of PAR bands using a coupled leaf-canopy radiative transfer model (PROSPECT-DART), considering a suite of tropical tree species over a wide range of assumed leaf chlorophyll contents. The analyses simulate hyperspectral images (5 nm bandwidth) of individual tree crowns for the selected background (concrete vs. grass) and illumination conditions. We then use an artificial neural network-based method to partition sunlit vs. shaded leaves within each crown, as the latter have lower fAPAR and fIPAR values. Our results show fAPAR of sunlit leaves decreases with the ratio of diffuse to direct scene irradiance (SKYL), while SKYL has minimal influence for shaded leaves. Both fAPAR and fIPAR decrease at more oblique solar zenith angles (SZA). Higher values of fAPAR and fIPAR occur with concrete backgrounds and the influence of the background is larger at higher diffuse ratio and solar zenith angles. The results show that fIPAR is typically 6–9% higher than fAPAR, and up to 14% higher for sunlit leaves with a concrete background at SKYL = 0. The differences between the fIPAR and fAPAR also depend on the health condition of the leaves, such as chlorophyll content. This study can improve the understanding of urban individual trees fAPAR/fIPAR and facilitate the development of protocols for fAPAR field measurements.

Effect of solar zenith angle on satellite cloud retrievals based on O2–O2 absorption band

ABSTRACT Recent studies have attempted to monitor global air quality with geostationary satellites. However, the monitoring accuracy is affected by a high solar zenith angle (SZA) during the day-to-night transition period, since the retrieval algorithms commonly use shortwave radiation. This study investigates the inaccuracy of background cloud retrievals for air quality at high SZAs that commonly use the O2–O2 absorption band at 477 nm. To this end, the differential optical absorption spectroscopy (DOAS) method with a radiative transfer model (RTM) simulations is used to quantify the amount of O2–O2 absorption by the cloud properties. The results show that the deviation of cloud retrievals from that retrieved at zero SZA increases exponentially with increasing SZA. The cloud deviation with increasing SZA can be larger for a smaller effective cloud fraction (ECF) and higher cloud centroid pressure (CCP). At an SZA of 80°, the ECF is deviated up to 0.06 (28%) and the CCP up to – 90 hPa (–13%). This cloud deviation originates from the reduction in the 477 nm reflectance in the DOAS calculation, which is primarily owing to the in-cloud absorption effect in addition to the out-of-cloud O2–O2 absorption effect. Therefore, minimizing these absorption effects in the DOAS method during day-to-night transition is crucial for improving cloud retrievals and, eventually, the overall reliability of air quality monitoring.

A Survey of Active Learning for Quantifying Vegetation Traits from Terrestrial Earth Observation Data.

The current exponential increase of spatiotemporally explicit data streams from satellitebased Earth observation missions offers promising opportunities for global vegetation monitoring. Intelligent sampling through active learning (AL) heuristics provides a pathway for fast inference of essential vegetation variables by means of hybrid retrieval approaches, i.e., machine learning regression algorithms trained by radiative transfer model (RTM) simulations. In this study we summarize AL theory and perform a brief systematic literature survey about AL heuristics used in the context of Earth observation regression problems over terrestrial targets. Across all relevant studies it appeared that: (i) retrieval accuracy of AL-optimized training data sets outperformed models trained over large randomly sampled data sets, and (ii) Euclidean distance-based (EBD) diversity method tends to be the most efficient AL technique in terms of accuracy and computational demand. Additionally, a case study is presented based on experimental data employing both uncertainty and diversity AL criteria. Hereby, a a simulated training data base by the PROSAIL-PRO canopy RTM is used to demonstrate the benefit of AL techniques for the estimation of total leaf carotenoid content (Cxc) and leaf water content (Cw). Gaussian process regression (GPR) was incorporated to minimize and optimize the training data set with AL. Training the GPR algorithm on optimally AL-based sampled data sets led to improved variable retrievals compared to training on full data pools, which is further demonstrated on a mapping example. From these findings we can recommend the use of AL-based sub-sampling procedures to select the most informative samples out of large training data pools. This will not only optimize regression accuracy due to exclusion of redundant information, but also speed up processing time and reduce final model size of kernel-based machine learning regression algorithms, such as GPR. With this study we want to encourage further testing and implementation of AL sampling methods for hybrid retrieval workflows. AL can contribute to the solution of regression problems within the framework of operational vegetation monitoring using satellite imaging spectroscopy data, and may strongly facilitate data processing for cloud-computing platforms.

A Novel Detection Algorithm for Impulse Spike Interferences in Cross-Track Infrared Sounder Short-Wave Spectral Radiance

Due to deficiencies in the operational spike detection, many moderate and strong impulse spike interferences constantly remain nondetected in the cross-track infrared sounder (CrIS) sensor data record (SDR) data at the short-wave infrared (SWIR) band over the South Atlantic Anomaly (SAA) region. This study characterizes imaginary radiance spectrum features in the presence of four intensities of spikes at three different hitting locations through radiative transfer model simulations. It is found that periodic ringing patterns appear in the imaginary spectral radiance. Following it, we establish a novel spike detection algorithm based on the spiky feature of power spectral density (PSD) in detrended and normalized imaginary SWIR CrIS spectral radiances. The performance of the algorithm is assessed by applying it to the multiple years of the SNPP and NOAA-20 CrIS full spectral resolution (FSR) radiance data. Approximate 300 and 500 spike events per day are detected for SNPP and NOAA-20 CrIS individually. The spikes can obviously increase the uncertainty in brightness temperature measurements at SWIR window channels that is determined by hitting location of the interferogram arms and the impulse intensities. Strong impulse spikes can trigger spectral radiance uncertainty by several times in SWIR window channels. Even the weakest impulse event can increase spectral radiance uncertainty by 4% in SWIR window channels. In addition, our simulations and one case study of CrIS data show that the spikes can degrade the quality of the CrIS data at the SWIR band with the largest impact at absorbing channels and the smallest at the window channels. The resultant brightness temperature change are 3 K for absorbing channels and up to 1 K for window channels in the SWIR band, demonstrating the significance flagging the CrIS data that are contaminated by the moderate and strong spikes that are missed in the operational data processing.

Modeling the directional anisotropy of fine-scale TIR emissions over tree and crop canopies based on UAV measurements

Land surface temperature (LST) is a vital parameter for the achievement of the surface energy budget and in thorough investigations of water cycle processes. Lightweight thermal infrared (TIR) sensors onboard unmanned aerial vehicles (UAVs) are rapidly becoming key instruments for extracting high-resolution LSTs given the flexibility they offer in capturing different scales. With this expansion, there has been increasing concern regarding the growing demand to obtain a mapping of normalized LST given the directional anisotropy (DA) of surface fine-scale emissions. To date, this topic suffers from a lack of deep analysis and practical solutions for characterizing the DA of fine-scale TIR data from UAV measurements over tree and crop canopies. In this paper, the first objective was to understand the pattern of brightness temperatures (BTs) DAs at a high spatial resolution by using UAV-based multiangle observations and three-dimensional (3D) radiative transfer model simulations. This study highlighted the need for first performing an angular normalization of the BTs of fine-scale pixels prior to any application, as these were easily affected by adjacent pixels and displayed broad spatial variability from 0.5 °C to 5.0 °C due to 3D occlusion. The second objective of the present study was to appraise the reliability of a modified kernel-driven model, in comparison to the model from which it was derived, with an additional kernel designed to mimic the adjacency effect, plus, a quadratic function used to simplify the estimate of the directional emissivity kernel. The root mean square error of the best fit between the measured UAV dataset and the modified kernel-driven model was approximately 0.65 °C, which proves its efficiency since the DA indexes of the BTs were about 1.40 °C. This outlined the role of the model to normalize from directional effects the camera image pixels and thereby deliver fine-scale BTs. In addition, results from LESS simulations also demonstrated the good performance of the modified kernel-driven model for simulating the DAs of thermal emissions for both tree and row-planted scenes.Index Terms—Land surface temperature, UAV, directional anisotropy, high spatial resolution.

Evaluation of the MODIS LAI/FPAR Algorithm Based on 3D-RTM Simulations: A Case Study of Grassland

Uncertainty assessment of the moderate resolution imaging spectroradiometer (MODIS) leaf area index (LAI) and the fraction of photosynthetically active radiation absorbed by vegetation (FPAR) retrieval algorithm can provide a scientific basis for the usage and improvement of this widely-used product. Previous evaluations generally depended on the intercomparison with other datasets as well as direct validation using ground measurements, which mix the uncertainties from the model, inputs, and assessment method. In this study, we adopted the evaluation method based on three-dimensional radiative transfer model (3D RTM) simulations, which helps to separate model uncertainty and other factors. We used the well-validated 3D RTM LESS (large-scale remote sensing data and image simulation framework) for a grassland scene simulation and calculated bidirectional reflectance factors (BRFs) as inputs for the LAI/FPAR retrieval. The dependency between LAI/FPAR truth and model estimation serves as the algorithm uncertainty indicator. This paper analyzed the LAI/FPAR uncertainty caused by inherent model uncertainty, input uncertainty (BRF and biome classification), clumping effect, and scale dependency. We found that the uncertainties of different algorithm paths vary greatly (−6.61% and +84.85% bias for main and backup algorithm, respectively) and the “hotspot” geometry results in greatest retrieval uncertainty. For the input uncertainty, the BRF of the near-infrared (NIR) band has greater impacts than that of the red band, and the biome misclassification also leads to nonnegligible LAI/FPAR bias. Moreover, the clumping effect leads to a significant LAI underestimation (−0.846 and −0.525 LAI difference for two clumping types), but the scale dependency (pixel size ranges from 100 m to 1000 m) has little impact on LAI/FPAR uncertainty. Overall, this study provides a new perspective on the evaluation of LAI/FPAR retrieval algorithms.

Aerosol Direct Radiative Effects over China Based on Long-Term Observations within the Sun–Sky Radiometer Observation Network (SONET)

To investigate aerosol radiative effects, the Sun–Sky Radiometer Observation Network (SONET) has performed long-term observations of columnar atmospheric aerosol properties at 20 distributed stations around China. The aerosol direct radiative forcing (RF) and efficiency (RFE, the rate at which the atmosphere is forced per unit of aerosol optical depth) were estimated using radiative transfer model simulations based on the ground-based observations dating back to 2009. Results of multi-year monthly average RF illustrated that: the dust-dominant aerosol population at arid and semi-arid sites exerted moderate cooling effects (−8.0~−31.2 W/m2) at the top and bottom of atmosphere (TOA and BOA); RF at continental background site was very weak (−0.8~−2.9 W/m2); fine-mode dominant aerosols at urban and suburban sites exerted moderate cooling effects (−9.3~−24.1 W/m2) at TOA but more significant cooling effects (−19.4~−50.6 W/m2) at BOA; RF at coastal sites was comparable with values of urban sites (−5.5~−19.5 W/m2 at TOA, and −15.6~−44.6 W/m2 at BOA), owing to combined influences by marine and urban–industrial aerosols. Differences between RFE at TOA and BOA indicated that coarse-mode dominant aerosols at arid, semi-arid, and continental background sites were less efficient to warm the atmosphere; but fine-mode dominant aerosols at urban, suburban, and coastal sites were shown to be more efficient to heat the atmosphere.

Comparison of SLSTR Thermal Emissive Bands Clear-Sky Measurements with Those of Geostationary Imagers

The Sentinel-3 series satellites belong to the European Earth Observation satellite missions for supporting oceanography, land, and atmospheric studies. The Sea and Land Surface Temperature Radiometer (SLSTR) onboard the Sentinel-3 satellites was designed to provide a significant improvement in remote sensing of skin sea surface temperature (SSTskin). The successful application of SLSTR-derived SSTskin fields depends on their accuracies. Based on sensor-dependent radiative transfer model simulations, geostationary Geostationary Operational Environmental Satellite (GOES-16) Advanced Baseline Imagers (ABI) and Meteosat Second Generation (MSG-4) Spinning Enhanced Visible and Infrared Imager (SEVIRI) brightness temperatures (BT) have been transformed to SLSTR equivalents to permit comparisons at the pixel level in three ocean regions. The results show the averaged BT differences are on the order of 0.1 K and the existence of small biases between them are likely due to the uncertainties in cloud masking, satellite view angle, solar azimuth angle, and reflected solar light. This study demonstrates the feasibility of combining SSTskin retrievals from SLSTR with those of ABI and SEVIRI.

Development of a Land Surface Temperature Retrieval Algorithm from GK2A/AMI

Land surface temperature (LST) is an important geophysical element for understanding Earth systems and land–atmosphere interactions. In this study, we developed a nonlinear split-window LST retrieval algorithm for the observation area of GEO-KOMPSAT-2A (GK2A), the next-generation geostationary satellite in Korea. To develop the GK2A LST retrieval algorithm, radiative transfer model simulation data, considering various impacting factors, were constructed. The LST retrieval algorithm was developed with a total of six equations as per day/night and atmospheric conditions (dry/normal/wet), considering the effects of diurnal variation of LST and atmospheric conditions on LST retrieval. The emissivity of each channel required for LST retrieval was calculated in real-time with the vegetation cover method using the consecutive 8-day cycle vegetation index provided by GK2A. The indirect validation of the results of GK2A LST with Moderate Resolution Imaging Spectroradiometer (MODIS) LST Collection 6 showed a high correlation coefficient (0.969), slightly warm bias (+1.227 K), and root mean square error (RMSE) (2.281 K). Compared to the MODIS LST, the GK2A LST showed a warm bias greater than +1.8 K during the day, but a relatively small bias (&lt;+0.7 K) at night. Based on the results of the validation with in situ measurements from the Tateno station of the Baseline Surface Radiation Network, the correlation coefficient was 0.95, bias was +0.523 K, and RMSE was 2.021 K.