The main objective of this research is to apply unmanned aerial system (UAS) data in synergy with field spectroradiometry for the accurate radiometric correction of Landsat-8 (L8) and Sentinel-2 (S2) imagery. The central hypothesis is that imagery acquired with multispectral UAS sensors that are well calibrated with highly accurate field measurements can fill in the scale gap between satellite imagery and conventional in situ measurements; this can be possible by sampling a larger area, including difficult-to-access land covers, in less time while simultaneously providing good radiometric quality. With this aim and by using near-coincident L8 and S2 imagery, we applied an upscaling workflow, whereby: (a) UAS-acquired multispectral data was empirically fitted to the reflectance of field measurements, with an extensive set of radiometric references distributed across the spectral domain; (b) drone data was resampled to satellite grids for comparison with the radiometrically corrected L8 and S2 official products (6S-LaSRC and Sen2Cor-SNAP, respectively) and the CorRad-MiraMon algorithm using pseudo-invariant areas, such as reflectance references (PIA-MiraMon), to examine their overall accuracy; (c) then, a subset of UAS data was used as reflectance references, in combination with the CorRad-MiraMon algorithm (UAS-MiraMon), to radiometrically correct the matching bands of UAS, L8, and S2; and (d) radiometrically corrected L8 and S2 scenes obtained with UAS-MiraMon were intercompared (intersensor coherence). In the first upscaling step, the results showed a good correlation between the field spectroradiometric measurements and the drone data in all evaluated bands (R² > 0.946). In the second upscaling step, drone data indicated good agreement (estimated from root mean square error, RMSE) with the satellite official products in visible (VIS) bands (RMSE_VIS < 2.484%), but yielded poor results in the near-infrared (NIR) band (RMSE_NIR > 6.688% was not very good due to spectral sensor response differences). In the third step, UAS-MiraMon indicated better agreement (RMSE_VIS < 2.018%) than the other satellite radiometric correction methods in visible bands (6S-LaSRC (RMSE < 2.680%), Sen2Cor-SNAP (RMSE < 2.192%), and PIA-MiraMon (RMSE < 3.130%), but did not achieve sufficient results in the NIR band (RMSE_NIR < 7.530%); this also occurred with all other methods. In the intercomparison step, the UAS-MiraMon method achieved an excellent intersensor (L8-S2) coherence (RMSE_VIS < 1%). The UAS-sampled area involved 51 L8 (30 m) pixels, 143 S2 (20 m) pixels, and 517 S2 (10 m) pixels. The drone time needed to cover this area was only 10 min, including areas that were difficult to access. The systematic sampling of the study area was achieved with a pixel size of 6 cm, and the raster nature of the sampling allowed for an easy but rigorous resampling of UAS data to the different satellite grids. These advances improve human capacities for conventional field spectroradiometry samplings. However, our study also shows that field spectroradiometry is the backbone that supports the full upscaling workflow. In conclusion, the synergy between field spectroradiometry, UAS sensors, and Landsat-like satellite data can be a useful tool for accurate radiometric corrections used in local environmental studies or the monitoring of protected areas around the world.

Radiometric correction is a prerequisite for generating high-quality scientific data, making it possible to discriminate between product artefacts and real changes in Earth processes as well as accurately produce land cover maps and detect changes. This work contributes to the automatic generation of surface reflectance products for Landsat satellite series. Surface reflectances are generated by a new approach developed from a previous simplified radiometric (atmospheric+topographic) correction model. The proposed model keeps the core of the old model (incidence angles and cast-shadows through a digital elevation model [DEM], Earth–Sun distance, etc.) and adds new characteristics to enhance and automatize ground reflectance retrieval. The new model includes the following new features: (1) A fitting model based on reference values from pseudoinvariant areas that have been automatically extracted from existing reflectance products (Terra MODIS MOD09GA) that were selected also automatically by applying quality criteria that include a geostatistical pattern model. This guarantees the consistency of the internal and external series, making it unnecessary to provide extra atmospheric data for the acquisition date and time, dark objects or dense vegetation. (2) A spatial model for atmospheric optical depth that uses detailed DEM and MODTRAN simulations. (3) It is designed so that large time-series of images can be processed automatically to produce consistent Landsat surface reflectance time-series. (4) The approach can handle most images, acquired now or in the past, regardless of the processing system, with the exception of those with extremely high cloud coverage. The new methodology has been successfully applied to a series of near 300 images of the same area including MSS, TM and ETM+ imagery as well as to different formats and processing systems (LPGS and NLAPS from the USGS; CEOS from ESA) for different degrees of cloud coverage (up to 60%) and SLC-off. Reflectance products have been validated with some example applications: time series robustness (for a pixel in a pseudoinvariant area, deviations are only 1.04% on average along the series), spectral signatures generation (visually coherent with the MODIS ones, but more similar between dates), and classification (up to 4 percent points better than those obtained with the original manual method or the CDR products). In conclusion, this new approach, that could also be applied to other sensors with similar band configurations, offers a fully automatic and reasonably good procedure for the new era of long time-series of spatially detailed global remote sensing data

The use of Pseudoinvariant Areas (PIA) makes it possible to carry out a reasonably robust and automatic radiometric correction for long time series of remote sensing imagery, as shown in previous studies for large data sets of Landsat MSS, TM, and ETM+ imagery. In addition, they can be employed to obtain more coherence among remote sensing data from different sensors. The present work validates the use of PIA for the radiometric correction of pairs of images acquired almost simultaneously (Landsat-7 (ETM+) or Landsat-8 (OLI) and Sentinel-2A (MSI)). Four pairs of images from a region in SW Spain, corresponding to four different dates, together with field spectroradiometry measurements collected at the time of satellite overpass were used to evaluate a PIA-based radiometric correction. The results show a high coherence between sensors (r² = 0.964) and excellent correlations to in-situ data for the MiraMon implementation (r² > 0.9). Other methodological alternatives, ATCOR3 (ETM+, OLI, MSI), SAC-QGIS (ETM+, OLI, MSI), 6S-LEDAPS (ETM+), 6S-LaSRC (OLI), and Sen2Cor-SNAP (MSI), were also evaluated. Almost all of them, except for SAC-QGIS, provided similar results to the proposed PIA-based approach. Moreover, as the PIA-based approach can be applied to almost any image (even to images lacking of extra atmospheric information), it can also be used to solve the robust integration of data from new platforms, such as Landsat-8 or Sentinel-2, to enrich global data acquired since 1972 in the Landsat program. It thus contributes to the program’s continuity, a goal of great interest for the environmental, scientific, and technical community.

Month-to-month air temperature (T_air) surfaces are increasingly demanded to feed quantitative models related to a wide range of fields, such as hydrology, ecology or climate change studies. Geostatistical interpolation techniques provide such continuous and objective surfaces of climate variables, while the use of remote sensing data may improve the estimates, especially when temporal resolution is detailed enough. The main goal of this study is to propose an empirical methodology for improving the month-to-month T_air mapping (minimum and maximum) using satellite land surface temperatures (LST) besides of meteorological data and geographic information. The methodology consists on multiple regression analysis combined with the spatial interpolation of residual errors using the inverse distance weighting. A leave-one-out cross-validation procedure has been included in order to compare predicted with observed values. Different operational daytime and nighttime LST products corresponding to the four months more characteristic of the seasonal dynamics of a Mediterranean climate have been considered for a thirteen-year period. The results can be considered operational given the feasibility of the models employed (linear dependence on predictors that are nowadays easily available), the robustness of the leave-one-out cross-validation procedure and the improvement in accuracy achieved when compared to classical T_air modeling results. Unlike what is considered by most studies, it is shown that nighttime LST provides a good proxy not only for minimum T_air, but also for maximum T_air. The improvement achieved by the inclusion of remote sensing LST products was higher for minimum T_air (up to 0.35 K on December), especially over forests and rugged lands. Results are really encouraging, as there are generally few meteorological stations in zones with these characteristics, clearly showing the usefulness of remote sensing to improve information about areas that are difficult to access or simply with a poor availability of conventional meteorological data.