Retrieval Of Surface Reflectance Research Articles

Retrieval of pseudo-BRDF-adjusted surface reflectance at 440 nm from the Geostationary Environmental Monitoring Spectrometer (GEMS)

Abstract. In satellite remote sensing applications, enhancing the precision of level 2 (L2) algorithms relies heavily on the accurate estimation of the surface reflectance across the ultraviolet (UV) to visible (VIS) spectrum. However, the mutual dependence between the L2 algorithms and the surface reflectance retrieval poses challenges, necessitating an alternative approach. To address this issue, many satellite algorithms generate Lambertian-equivalent reflectivity (LER) products as a priori surface reflectance data; however, this often results in an underestimation of these data. This study is the first to assess the applicability of background surface reflectance (BSR), derived using a semi-empirical bidirectional reflectance distribution function (BRDF) model, in an operational environmental satellite algorithm. This study pioneered the application of the BRDF model to hyperspectral satellite data at 440 nm, aiming to provide more realistic preliminary surface reflectance data. In this study, the Geostationary Environment Monitoring Spectrometer (GEMS) data were used, and a comparative analysis of the GEMS BSR and GEMS LER retrieved in this study revealed an improvement in the relative root mean squared error (rRMSE) accuracy of 3 %. Additionally, a time series analysis across diverse land types indicated a greater stability exhibited by the BSR than by the LER. For further validation, the BSR was compared with other LER databases using ground-truth data, yielding superior simulation performance. These findings present a promising avenue for enhancing the accuracy of surface reflectance retrieval from hyperspectral satellite data, thereby advancing the practical applications of satellite remote sensing algorithms.

Time-Series Cross-Radiometric Calibration and Validation of GF-6/WFV Using Multi-Site

The GaoFen6 (GF-6) satellite, equipped with a wide full-swath (WFV) sensor, offers high spatial resolution and extensive coverage, making it widely utilized in agricultural and forestry classification, land resource monitoring, and other fields. Accurate on-orbit radiometric calibration of GF-6/WFV is crucial for these quantitative applications. Currently, the absolute radiometric calibration of GF-6/WFV relies primarily on vicarious calibration conducted by the China Center for Resources Satellite Data and Application (CRESDA). However, annual vicarious calibration may not adequately capture the radiometric performance of GF-6/WFV due to performance degradation. Therefore, increasing the frequency of on-orbit radiometric calibration throughout the lifetime of GF-6/WFV is essential. This study proposes a method for conducting long-term cross-radiometric calibrations of GF-6/WFV by taking the multispectral imager (MSI) onboard the Sentinel-2 satellite as a reliable reference sensor and the sites from RadCalNet as reference ground targets. Firstly, we conducted 62 on-orbit cross-radiometric calibrations of GF-6/WFV since its launch by tracking with the Sentinel-2/MSI sensor after correcting the discrepancy spectrum and solar zenith angle. Then, validation of cross-radiometric calibration results against RadCalNet products indicated an average absolute relative error between 3.55% and 4.64%. Cross-validation with additional reference sensors, including Landsat-8/OLI and MODIS, confirmed the reliability of calibration, demonstrating relative differences from GF-6/WFV of less than 5%. Furthermore, the overall uncertainty of the cross-radiometric calibration was estimated to be from 4.08% to 4.89%. Finally, trend analysis of the time-series radiometric performance was also conducted and revealed an annual degradation rate ranging from 0.57% to 2.31%. This degradation affects surface reflectance retrieval, introducing a bias of approximately 0.0073 to 0.0084. Our findings highlight the operational effectiveness of the proposed method in achieving long-time-series on-orbit radiometric calibration and degradation monitoring of GF-6/WFV. The study also demonstrates that the radiometric performance of GF-6/WFV is relatively stable and suitable for further quantitative applications, especially for long-term monitoring applications.

Quantitative assessment of the potential of optimal estimation for aerosol retrieval from geostationary weather satellites in the frame of the iAERUS‐GEO algorithm

AbstractSatellite remote sensing enables the study of atmospheric aerosols at large spatial scales, with geostationary platforms making this possible at sub‐daily frequencies. High‐temporal‐resolution aerosol observations can be made from geostationary data by using robust numerical inversion methods such as the widely‐used optimal estimation (OE) theory. This is the case of the instantaneous Aerosol and surfacE Retrieval Using Satellites in GEOstationary orbit (iAERUS‐GEO) algorithm, which successfully retrieves aerosol optical depth (AOD) maps from the Meteosat Second Generation weather satellite based on a simple implementation of the OE approach combined with the Levenberg–Marquardt method. However, the exact gain in inversion performances that can be obtained from the multiple and more advanced possibilities offered by OE is not well documented in the current literature. Against this background, this article presents the quantitative assessment of OE for the future improvement of the iAERUS‐GEO algorithm. To this end, we use a series of comprehensive experiments based on AOD maps retrieved by iAERUS‐GEO using different OE implementations, and ground‐based observations used as reference data. First, we assess the varying importance in the inversion process of satellite observations and a priori information according to the content of satellite aerosol information. Second, we quantify the gain of AOD estimation in log space versus linear space in terms of accuracy, AOD distribution and number of successful retrievals. Finally, we evaluate the accuracy improvement of simultaneous AOD and surface reflectance retrieval as a function of the regions covered by the Meteosat Earth's disk.

An Algorithm Developed for Smallsats Accurately Retrieves Landsat Surface Reflectance Using Scene Statistics

Closed-form Method for Atmospheric Correction (CMAC) is software that overcomes radiative transfer method problems for smallsat surface reflectance retrieval: unknown sensor radiance responses because onboard monitors are omitted to conserve size/weight, and ancillary data availability that delays processing by days. CMAC requires neither and retrieves surface reflectance in near real time, first mapping the atmospheric effect across the image as an index (Atm-I) from scene statistics, then reversing these effects with a closed-form linear model that has precedence in the literature. Five consistent-reflectance area-of-interest targets on thirty-one low-to-moderate Atm-I images were processed by CMAC and LaSRC. CMAC retrievals accurately matched LaSRC with nearly identical error profiles. CMAC and LaSRC output for paired images of low and high Atm-I were then compared for three additional consistent-reflectance area-of-interest targets. Three indices were calculated from the extracted reflectance: NDVI calculated with red (standard) and substitutions with blue and green. A null hypothesis for competent retrieval would show no difference. The pooled error for the three indices (n = 9) was 0–3% for CMAC, 6–20% for LaSRC, and 13–38% for uncorrected top-of-atmosphere results, thus demonstrating both the value of atmospheric correction and, especially, the stability of CMAC for machine analysis and AI application under increasing Atm-I from climate change-driven wildfires.

Spatially constrained atmosphere and surface retrieval for imaging spectroscopy

Imaging spectrometers provide important Earth surface data for terrestrial and aquatic ecology, hydrology, geology through retrieved surface reflectance spectra. Surface reflectance is generally recovered from measured radiance through model inversions of the surface and atmospheric system. Most retrieval approaches treat all pixels independently and ignore spatial correlations in the atmosphere for neighboring pixels, resulting in surface-related biases in the retrieved surface and atmospheric parameters. In this work, we present a mathematical framework to more accurately retrieve surface reflectance and atmospheric parameters by leveraging the expected spatial smoothness of the atmosphere. This framework retrieves the full spatial–spectral reflectance data cube in a computationally efficient manner and is to our knowledge the first approach capable of accurately finding this global solution for large imaging spectroscopy datasets. We implement this mathematical framework in an example proposed algorithm which we call spatially constrained optimal estimation (SCOE). We show that SCOE improves the overall surface reflectance retrieval error, reduces surface-related biases in the retrieved surface reflectance and atmospheric parameters, and is more robust to noise than a traditional pixel-by-pixel approach and an operational empirical line approach in simulated and experimental results. This reduction in surface reflectance error will impact science data products across application areas, making imaging spectroscopy algorithms more robust for global acquisitions.

Impact of Atmospheric Correction on Classification and Quantification of Seagrass Density from WorldView-2 Imagery.

Mapping the seagrass distribution and density in the underwater landscape can improve global Blue Carbon estimates. However, atmospheric absorption and scattering introduce errors in space-based sensors' retrieval of sea surface reflectance, affecting seagrass presence, density, and above-ground carbon () estimates. This study assessed atmospheric correction's impact on mapping seagrass using WorldView-2 satellite imagery from Saint Joseph Bay, Saint George Sound, and Keaton Beach in Florida, USA. Coincident in situ measurements of water-leaving radiance (), optical properties, and seagrass leaf area index (LAI) were collected. Seagrass classification and the retrieval of LAI were compared after empirical line height (ELH) and dark-object subtraction (DOS) methods were used for atmospheric correction. DOS left residual brightness in the blue and green bands but had minimal impact on the seagrass classification accuracy. However, the brighter reflectance values reduced LAI retrievals by up to 50% compared to ELH-corrected images and ground-based observations. This study offers a potential correction for LAI underestimation due to incomplete atmospheric correction, enhancing the retrieval of seagrass density and above-ground Blue Carbon from WorldView-2 imagery without in situ observations for accurate atmospheric interference correction.

Evaluation of surface reflectance retrieval over diverse surface types using SREM algorithm in varied aerosol conditions for coarse to medium resolution data from multiple spaceborne sensors

ABSTRACT Surface reflectance (SR) is one of the most important components for deriving different biophysical parameters from remotely sensed data. For accurate retrieval of SR, it is essential to accurately estimate the degree of atmospheric alteration to the target signals reaching at-sensor. Simplified and robust SR estimation method (SREM) is one such approximate estimation method that could demonstrate the retrieval of SR from top-of-atmosphere (ToA) signals, without using any precalculated LUTs or aerosol information. The present study attempts to assess the accuracy of the SREM algorithm over varied Indian land surfaces using data from Landsat 8, Sentinel 2, AWiFS and LISS III. The results indicate that while SREM performs quite well for Sentinel 2 over most land surface types, the performance is good for Landsat 8 as well except over water and in coastal aerosol band. SREM retrievals also show high correlation for AWiFS and LISS III but the estimated values are much lower than level 2 SR due to smaller ToA reflectance values. Furthermore, while aerosol load does impact the performance adversely, it is least impactful for Sentinel 2 and LISS III datasets and most over Landsat 8, particularly over water. Additionally, while AWiFS and LISS III SR retrievals are underestimated at all spectral bands, retrieval is slightly underestimated from green to SWIR band in case of Sentinel 2 and SWIR band in case of Landsat 8. Overall, the findings demonstrate that while SREM will need to be modified to adapt to LISS III and AWiFS retrievals, it can, nevertheless, be readily used as a simplified method for quick and approximate estimation of SR over the Indian region using data from Landsat 8 and Sentinel 2.

Evaluation of Low-Cost Radiometer for Surface Reflectance Retrieval and Orbital Sensor’s Validation

This paper evaluates the Arable Mark 2 sensor, an automated and low-cost radiometer, for its potential to retrieve surface reflectance data and validate orbital sensors such as the Landsat-8 (L8) Operational Land Imager (OLI) Level 2 product. While orbital sensors are widely used for monitoring solar radiation changes, managing natural resources, and understanding climatic trends, atmospheric effects can make it challenging to obtain accurate measurements. Equipped with multiple sensors, including long-wave and short-wave radiometers, the Arable Mark 2 sensor can measure upwelling and downwelling irradiance to calculate surface reflectance. To assess the accuracy and consistency of the Arable Mark 2 sensor, the study performed a cross-calibration using a ground truth measurement collected with the Analytical Spectral Device (ASD) as the reference point. Additionally, a spectral band adjustment factor (SBAF) was applied across the calibrated Arable surface reflectance to compare it against the orbital sensor. An automated library aided in calculating SBAF for the days with unavailable hyperspectral data. The study found that the Arable Mark 2 sensor can provide accurate surface reflectance data that can be used for orbital sensor validation. The Arable sensor was successfully calibrated against the ASD FieldSpec with an average difference of less than 1/10 reflectance unit (reflectance unit = 0.01) for the blue, green, yellow, and red bands. The red-edge and NIR-1 bands showed an average difference of less than 1/2 reflectance units, while the NIR-2 band had an average difference of less than 1/10 reflectance unit of calibration accuracy. The calibrated Arable surface reflectance data was then compared against orbital sensor surface reflectance data, and the results showed good agreement between the two datasets. The study concludes that the low-cost and automated nature of the Arable Mark 2 sensor makes it a promising tool for surface reflectance retrieval and orbital sensor validation.

Atmospheric Correction Inter-comparison eXercise, ACIX-II Land: An assessment of atmospheric correction processors for Landsat 8 and Sentinel-2 over land

The correction of the atmospheric effects on optical satellite images is essential for quantitative and multi-temporal remote sensing applications. In order to study the performance of the state-of-the-art methods in an integrated way, a voluntary and open-access benchmark Atmospheric Correction Inter-comparison eXercise (ACIX) was initiated in 2016 in the frame of Committee on Earth Observation Satellites (CEOS) Working Group on Calibration & Validation (WGCV). The first exercise was extended in a second edition wherein twelve atmospheric correction (AC) processors, a substantially larger testing dataset and additional validation metrics were involved. The sites for the inter-comparison analysis were defined by investigating the full catalogue of the Aerosol Robotic Network (AERONET) sites for coincident measurements with satellites' overpass. Although there were more than one hundred sites for Copernicus Sentinel-2 and Landsat 8 acquisitions, the analysis presented in this paper concerns only the common matchups amongst all processors, reducing the number to 79 and 62 sites respectively. Aerosol Optical Depth (AOD) and Water Vapour (WV) retrievals were consequently validated based on the available AERONET observations. The processors mostly succeeded in retrieving AOD for relatively light to medium aerosol loading (AOD < 0.2) with uncertainties <0.08, while the overall uncertainty values were typically 0.23 ± 0.15. Better performances were observed for WV retrievals with >90% of the results falling within the suggested empirical specifications and with the Root Mean Square Error (RMSE) being mostly <0.25 g/cm2. Regarding Surface Reflectance (SR) validation two main approaches were followed. For the first one, a simulated SR reference dataset was computed over all of the test sites by using the 6SV (Second Simulation of the Satellite Signal in the Solar Spectrum vector code) full radiative transfer modelling (RTM) and AERONET measurements for the required aerosol variables and water vapour content. The performance assessment demonstrated that the retrievals were not biased for most of the bands. The uncertainties ranged from approximately 0.003 to 0.01 (excluding B01) for the best performing processors in both sensors' analyses. For the second one, measurements from the radiometric calibration network RadCalNet over La Crau (France) and Gobabeb (Namibia) were involved in the validation. The performance of the processors was in general consistent across all bands for both sensors and with low standard deviations (<0.04) between on-site and estimated surface reflectance. Overall, our study provides a good insight of AC algorithms' performance to developers and users, pointing out similarities and differences for AOD, WV and SR retrievals. Such validation though still lacks of ground-based measurements of known uncertainty to better assess and characterize the uncertainties in SR retrievals.

Assessing the Cloud Adjacency Effect on Retrieval of the Ground Surface Reflectance from MODIS Satellite Data for the Baikal Region

The cloud adjacency effect on surface reflectance retrievals for the region of the Russian Federation with coordinates 51–54∘ N, 103–109∘ E including the southern part of Lake Baikal for the period of 1–23 July 2021 is assessed in this paper. The method is based on the computer program for statistical simulation of radiative transfer in the atmosphere with the stochastic cloud field including a deterministic gap of a given radius. The results of this program are then used in the interpolation formula. Masks of cloudless pixels, for which the cloud adjacency effect (CAE) changes the ground surface reflectance by more than 0.005, are constructed. The analysis of the resulting CAE radii shows that the average radius is 13.7 km for MODIS band 8, 11.2 km for band 3, 8.4 km for band 4, 7.2 km for band 1, and 7 km for band 2. For the considered MODIS images and bands, the pixels with strong CAE make up from 2.8 to 100% of the total number of cloudless pixels. The correlation coefficients between the initial data and the CAE radius suggest that the cloud optical depth, cloud cover index, and ground surface reflectance exert the major influence on the considered images. A simplified approximation equation for the CAE radius as a function of the cloud optical depth, cloud cover index, and surface reflectance is derived. The analysis of the approximation shows that for the considered images, the CAE radius decreases nearly linearly with wavelength for low reflective surfaces. However, for high reflective surfaces, its wavelength dependence is nonlinear.

Evaluation of MODIS Dark Target AOD Product with 3 and 10 km Resolution in Amazonia

The techniques and analyses employed by remote sensing provide key information about atmospheric particle properties at regional and global scales. However, limitations in optical spectral models used to represent the different types of aerosols in the atmosphere and their effects (direct and indirect) are still one of the major causes of sources of uncertainties and substantial impacts in climate prediction. There are no studies yet in South America, especially in the Amazon Basin, that have evaluated the advantages, disadvantages, inconsistencies, applicability, and suitability of the MODIS sensor (Moderate Resolution Imaging Spectroradiometer) destined for monitoring the ambient aerosol optical thickness over rivers and continents. In this study, the results of the DT (Dark Target) algorithm for products with 3 km and 10 km resolutions were systematically evaluated for six sites in the Amazon rainforest. The comparisons between the products were carried out with the AERONET (Aerosol Robotic Network) measurements, which were used as reference. Statistical parameters between AERONET vs. MODIS were also evaluated based on biomass burning records in the site regions. Here, the DT 10 km product showed satisfactory performance for the Amazon region, with observations between the expected error (EE) limits above 66%, in addition to R &gt; 0.8 and RMSE &lt; 0.3. However, the regional analysis for the two sites in the central and southern regions of the Amazon basin did not have the same performance, where the results showed an EE of 24 and 47%, respectively. The DT 3 km product did not perform well in any site, with an EE below 50%. Both products overestimated the AOD, but the 3 km product overestimated it approximately four times more due to its algorithm setup. Thus, we recommend the 10 km product for general analysis in Amazonia. Regional biomass burning records showed a direct relationship with the AERONET vs. MODIS DT with overestimation of both products. All variations between products and sites were justified based on the difficulty of retrieving surface reflectance and the model selected for local aerosols. Improvements in the optical spectral model currently implemented in the algorithms, with more realistic representations of the main types of the aerosol present in the Amazon Basin, may contribute to better performance among the evaluated products.

Unified Topographic and Atmospheric Correction for Remote Imaging Spectroscopy

Models linking surface characteristics within incident solar radiation are inexorably dependent on the topography of the given region. To date, however, most operational surface reflectance retrievals treat this dependence by assuming a flat terrain, leading to significant deviations in the estimated reflectance. Here, we demonstrate that incorporating dynamic topography directly into the joint surface and atmospheric model during retrievals has several advantages. First, it allows for a more complete physical accounting of downwelling illumination, providing more accurate estimates of the absolute magnitude of reflectance. Second, it facilitates a superior resolution of the atmospheric state, most notably due to the confounding influence of atmospheric aerosols and unresolved topographic effects. Our methodology utilizes a practical, high-fidelity, model-driven approach to separate out diffuse and direct irradiation and account for topographic effects during the joint inversion of atmosphere and surface properties. We achieve this by enhancing the atmosphere/surface inversion to account for the radiative transfer effects of surface slope. We further demonstrate how uncertainties in topographic features can be quantified and leveraged within our formulation for a more realistic posterior uncertainty estimates. Our results demonstrate that the inclusion of topographic effects into the retrieval model reduces errors in the reflectance of an only moderately rugged terrain by more than 15%, and that a post hoc accounting of topography cannot achieve these same results.

Spatial Surface Reflectance Retrievals for Visible/Shortwave Infrared Remote Sensing via Gaussian Process Priors

Remote Visible/Shortwave Infrared (VSWIR) imaging spectroscopy is a powerful tool for measuring the composition of Earth’s surface over wide areas. This compositional information is captured by the spectral surface reflectance, where distinct shapes and absorption features indicate the chemical, bio- and geophysical properties of the materials in the scene. Estimating this surface reflectance requires removing the influence of atmospheric distortions caused by water vapor and particles. Traditionally reflectance is estimated by considering one location at a time, disentangling atmospheric and surface effects independently at all locations in a scene. However, this approach does not take advantage of spatial correlations between contiguous pixels. We propose an extension to a common Bayesian approach, Optimal Estimation, by introducing atmospheric correlations into the multivariate Gaussian prior. We show how this approach can be implemented as a small change to the traditional estimation procedure, thus limiting the additional computational burden. We demonstrate a simple version of the technique using simulations and multiple airborne radiance data sets. Our results show that the predicted atmospheric fields are smoother and more realistic than independent inversions given the assumption of spatial correlation and may reduce bias in the surface reflectance retrievals compared to post-process smoothing.

Aerosol models from the AERONET database: application to surface reflectance validation

Abstract. Aerosols play a critical role in radiative transfer within the atmosphere, and they have a significant impact on climate change. In this paper, we propose and implement a framework for developing an aerosol model using their microphysical properties. Such microphysical properties as the size distribution, the complex refractive index, and the percentage of sphericity are derived from the global AERosol RObotic NETwork (AERONET). These measurements, however, are typically retrieved when almucantar measurement procedures are performed (i.e., early mornings and late afternoons with clear sky) and might not have a temporal correspondence to a satellite overpass time, so a valid validation of satellite-derived products cannot be carried out. To address this problem of temporal inconsistency of satellite and ground-based measurements, we developed an approach to retrieve these microphysical properties (and the corresponding aerosol model) using the optical thickness at 440 nm, τ440, and the Ångström coefficient between 440 and 870 nm, α440–870. Such aerosol models were developed for 851 AERONET sites within the last 28 years. Obtained results suggest that empirically microphysical properties can be retrieved with uncertainties of up to 23 %. An exception is the imaginary part of the refractive index ni, for which the derived uncertainties reach up to 38 %. These specific parametric models of aerosol can be used for the studies when retrieval of microphysical properties is required as well as validation of satellite-derived products over land. Specifically, we demonstrate the usefulness of the aerosol models to validate surface reflectance records over land derived from optical remote sensing sensors. We then quantify the propagation of uncertainties in the surface reflectance due to uncertainties with the aerosol model retrieval that is used as a reference from radiative transfer simulations. Results indicate that individual aerosol microphysical properties can impact uncertainties in surface reflectance retrievals between 3.5 × 10−5 to 1 × 10−3 (in reflectance units). The overall impact of microphysical properties combined yields an overall uncertainty in surface reflectance &lt; 0.004 (in reflectance units). That corresponds, for example, to 1 to 3 % of the retrieved surface reflectance in the red spectral band (620–670 nm) by the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. These uncertainty values are well below the specification (0.005 + 0.05ρ; ρ is the retrieved surface reflectance) used for the MODIS atmospheric correction.

Deriving Surface Reflectance From Visible/Near Infrared and Ultraviolet Satellite Observations Through the Community Radiative Transfer Model

Earth's surface reflectance is an important parameter affecting ultraviolet (UV) and visible (VIS) radiance calculations at the top of the atmosphere because many UV and VIS channels can acquire information about the surface and atmosphere. This article provides the theoretical basis for deriving the surface reflectance from satellite-measured UV and VIS observations at window and lower sounding channels with the help of the community radiative transfer model (CRTM) and collocated atmospheric profiles such as ozone, water vapor, and aerosols. Cirrus cloud may be included in the calculation as long as the observations contain enough reflected radiation from the surface. An explicit equation with three scalar parameters <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\alpha $</tex-math></inline-formula> , <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\beta $</tex-math></inline-formula> , and <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\delta $</tex-math></inline-formula> is obtained for users to calculate Lambertian surface reflectance from the observation. The expressions for the three parameters are somewhat complicated and computationally expansive. We found a simple and smart way that can exactly calculate the three parameters with quasi-linear functions. Numerical experiments using the CRTM simulations have demonstrated the algorithm accuracy for the surface reflectance retrieval better than 2.0E-14. As a case study, measured surface reflectance and the derived surface reflectance over desert from satellite UV measurements are compared. The derived surface reflectance from Suomi National Polar-orbiting Partnership. Visible Infrared Imaging Radiometer Suite (VIIRS) observations and the VIIRS reflectance product are compared as well. In addition, this methodology can also be used to calculate microwave and infrared surface emissivity with scatterings and solar radiation by adding the surface Planck radiance at the surface temperature.

Impact of Tree Crown Transmittance on Surface Reflectance Retrieval in the Shade for High Spatial Resolution Imaging Spectroscopy: A Simulation Analysis Based on Tree Modeling Scenarios

With the advancement of high spatial resolution imaging spectroscopy, an accurate surface reflectance retrieval is needed to derive relevant physical variables for land cover mapping, soil, and vegetation monitoring. One challenge is to deal with tree shadows using atmospheric correction models if the tree crown transmittance Tc is not properly taken into account. This requires knowledge of the complex radiation mechanisms that occur in tree crowns, which can be provided by coupling the physical modeling of canopy radiative transfer codes (here DART) and the 3D representations of trees. First in this study, a sensitivity analysis carried out on DART simulations with an empirical 3D tree model led to a statistical regression predicting Tc from the tree leaf area index (LAI) and the solar zenith angle with good performances (RMSE ≤ 4.3% and R2 ≥ 0.91 for LAI ≤ 4 m2.m−2). Secondly, more realistic 3D voxel-grid tree models derived from terrestrial LiDAR measurements over two trees were considered. The comparison of DART-simulated Tc from these models with the previous predicted Tc over 0.4–2.5 µm showed three main sources of inaccuracy quoted in order of importance: (1) the global tree geometry shape (mean bias up to 21.5%), (2) the transmittance fraction associated to multiple scattering, Tscat (maximum bias up to 13%), and (3) the degree of realism of the tree representation (mean bias up to 7.5%). Results showed that neglecting Tc leads to very inaccurate reflectance retrieval (mean bias &gt; 0.04), particularly if the background reflectance is high, and in the near and shortwave infrared – NIR and SWIR – due to Tscat. The transmittance fraction associated to the non-intercepted transmitted light, Tdir, can reach up to 95% in the SWIR, and Tscat up to 20% in the NIR. Their spatial contributions computed in the tree shadow have a maximum dispersion of 27% and 8% respectively. Investigating how to approximate Tdir and Tscat spectral and spatial variability along with the most appropriate tree 3D modeling is crucial to improve reflectance retrieval in tree shadows when using atmospheric correction models.

Satellite-observed changes of surface spectral reflectances due to solar farming and the implication for radiation budget

Increased deployment of solar farms in the last decade is coincident with continued monitoring of surface spectral reflectances by the moderate resolution imaging spectroradiometers (MODIS) aboard NASA Aqua satellite. By analyzing the surface spectral reflectances retrieved at seven MODIS shortwave bands, this study quantifies the changes of surface spectral reflectances caused by solar farms commissioned in the southwestern US. Our case study using one solar farm indicates a 20%–25% reduction of surface reflectance over the seven MODIS bands due to the solar panel installation. Radiative transfer calculation shows that such a reduction in surface spectral reflectance leads to a ∼23% decrease in the upward shortwave broadband flux at the surface and a ∼14%–18% decrease in the clear-sky reflected shortwave flux at the top of atmosphere (TOA). Though the MODIS surface reflectance retrievals can be affected by other factors such as aerosol or thin-cirrus contaminations, five out of six solar farm sites analyzed here show unambiguously detectable changes of surface reflectances due to solar panel installations. The brighter the surface is at a given MODIS band before the solar-farm construction, the larger the spectral reflectance reduction tends to be. If all the bright deserts, which occupies about 4% of the global surface, were covered with the solar panels similar to what has been studied here, the estimated instantaneous TOA shortwave radiative forcing would be no more than 1.1 W m−2.

Evaluation of atmospheric correction methods for low to high resolutions satellite remote sensing data

The objective of atmospheric correction is to retrieve surface reflectance from the top of atmosphere (TOA) reflectance. However, estimating surface reflectance from the TOA reflectance satellite data requires knowledge about the state of the atmosphere (e.g., water vapor and ozone) and the contribution of aerosols to the atmospheric path radiance. Thus, obtaining precise measurements of these parameters, which is very difficult, is crucial for accurate estimation of surface reflectance. The SREM (Simplified and Robust Surface Reflectance Estimation Method) is a physical-based atmospheric correction method based on the Radiative transfer (RT) equations of the second simulation of the Satellite Signal in Solar Spectrum (6SV). Essentially the SREM is a simplified version of 6SV which does not require Aerosol Optical Depth (AOD), aerosol type, water vapor, and ozone. An initial study showed accuracy comparable to the Landsat operational Surface Reflectance Products (SRProd) which is generated through different RT models using AOD, water vapor, and ozone data. To further validate the SREM under varying atmospheric conditions and at different spatial resolutions, an independent Reference Surface Reflectance (SRRef) dataset was generated using the AERONET (Aerosol Robotic Network) measurements as input to the 6SV RT model. The surface reflectances estimated by SREM (SRSREM) and SRProd from Planet Scope (PS, at 3 m spatial resolution), Sentinel-2 AB (S2AB) Multi-spectral Instrument (MSI, at 10 to 60 m spatial resolution), and Landsat-8 (L8) operational Land Imager (OLI, at 30 m spatial resolution) were validated against SRRef. Results showed that SRSREM performed similar to the SRProd of PS, S2AB MSI, and L8 OLI against SRRef. An inferior performance (R of 0.35 and 0.57) of L8 OLI's SRProd in the coastal blue (SB1) and blue (SB2) bands was observed, compared to SREM. The comparison of SRSREM with SRProd reveals the robustness of SREM, without using AOD, water vapor, and ozone data, for estimation of surface reflectance for all RT models tested. For some dates, SRRef and the SRProd under-corrected and produced higher values than the TOA reflectance, even when the atmosphere was clear but this was not the case for SREM. Analysis of surface reflectance estimation in shadowed areas revealed that the SRRef and SRProd had mainly negative values in coastal blue and blue bands for L8 OLI, while no negative SR value was observed for SREM in any band. These results recommend the utilization of SREM for the provision of surface reflectance products across a range of sensors

Uncertainty quantification for a global imaging spectroscopy surface composition investigation

Airborne and orbital imaging spectroscopy can facilitate the quantification of chemical and physical attributes of surface materials through analysis of spectral signatures. Prior to analysis, estimates of surface reflectance must be inferred from radiance measurements in a process known as atmospheric correction, which compensates for the distortion of the electromagnetic signal by the atmosphere. Inaccuracies in the correction process can alter characteristic spectral signatures, leading to subsequent mischaracterization of surface properties. Global observations pose new challenges for mapping surface composition, as varied atmospheric conditions and surface biomes challenge traditional atmospheric correction methods. Recent work adopted an optimal estimation (OE) approach for retrieving surface reflectance from observed radiance measurements, providing the reflectance estimates with a posterior probability. This work incorporates these input probabilities to improve the accuracy of surface feature measurements. We demonstrate this using a generic feature-fitting method that is applicable to a wide range of Earth surface studies including geology, ecosystem studies, hydrology and urban studies. Specifically, we use a probabilistic framework based on generalized Tikhonov-regularized least squares, a rigorous formulation for appropriate weighting of features by their observation uncertainty and leveraging of prior knowledge of material abundance for improving estimation accuracy. We demonstrate the validity of this procedure and quantify the increase in model performance by simulating expected accuracies in the reflectance estimation. To evaluate global uncertainties in mineral estimation, we simulate observations representative of the expected global range of atmospheric water vapor and aerosol levels, and characterize the sensitivity of our procedure to those quantities.

DRONE DATA ATMOSPHERIC CORRECTION CONCEPT FOR MULTI- AND HYPERSPECTRAL IMAGERY – THE DROACOR MODEL

Abstract. Remote sensing with unmanned aerial vehicles (UAVs) is a fast and cost-efficient tool for mapping and environmental monitoring. The sensors are operated at low flight altitudes, usually below 500 m above ground, leading to spatial resolutions up to the centimeter range. This type of data causes new challenges in atmospheric compensation and surface reflectance retrieval. Based on these specific boundary conditions, a new drone based atmospheric correction concept (DROACOR) is proposed, which is designed for currently available UAV based sensors. It is suited for multispectral visible/near infrared sensors as well as hyperspectral instruments covering the 400–1000 nm spectral region or the 400–2500 nm spectrum. The goal of the development is a fully automatic processor which dynamically adjusts to the given instrument and the atmospheric conditions. Optionally, irradiance measurements from simultaneously measured cosine receptors or from in-field reference panels can be taken into account to improve the processing quality by adjusting the irradiance parameter or by performing an in-flight vicarious calibration. Examples of DROACOR processing results are presented for a multispectral image data set and a hyperspectral data set, both acquired at variable flight altitudes. The resulting spectra show the applicability of the methods for both sensor types and an accuracy level below 2.5% reflectance units.