Bayesian cloud detection for 37 Years of Advanced Very High Resolution Radiometer (AVHRR) Global Area Coverage (GAC) data

Cloud detection is a source of significant errors in retrieval of sea surface temperature (SST). We apply a Bayesian cloud detection scheme to 37 years of Advanced Very High Resolution Radiometer (AVHRR) Global Area Coverage (GAC) data, which is an important source of multi-decadal global SST information. The Bayesian scheme calculates a probability of clear-sky for each image pixel, conditional on the satellite observations and prior probablity. We compare the cloud detection performance to the operational Clouds from AVHRR Extended algorithm (CLAVR-x), as a measure of improvement from reduced cloud-related errors. To do this we use sea surface temperature differences between satellite retrievals and in-situ observations from drifting buoys and the Global Tropical Moored Buoy Array (GTMBA). The Bayesian scheme reduces the absolute difference between the mean and median SST biases and reduces the standard deviation of the SST differences by ~10 % for both daytime and nighttime retrievals. These reductions are indicative of removing cloud contaminated outliers in the distribution, as these fall only on one side of the distribution forming a cold tail. At a probability threshold of 0.9 typically used to determine a binary cloud mask for SST retrieval, the Bayesian mask also reduces the robust standard deviation by ~5-10 % during the day, in comparison with the operational cloud mask. This shows an improvement in the central distribution of SST differences for daytime retrievals

Satellite-based time-series of sea-surface temperature since 1981 for climate applications

A climate data record of global sea surface temperature (SST) spanning 1981–2016 has been developed from 4 × 10^12 satellite measurements of thermal infra-red radiance. The spatial area represented by pixel SST estimates is between 1 km^2 and 45 km^2. The mean density of good-quality observations is 13 km^−2 yr^−1. SST uncertainty is evaluated per datum, the median uncertainty for pixel SSTs being 0.18 K. Multi-annual observational stability relative to drifting buoy measurements is within 0.003 K yr^−1 of zero with high confidence, despite maximal independence from in situ SSTs over the latter two decades of the record. Data are provided at native resolution, gridded at 0.05° latitude-longitude resolution (individual sensors), and aggregated and gap-filled on a daily 0.05° grid. Skin SSTs, depth-adjusted SSTs de-aliased with respect to the diurnal cycle, and SST anomalies are provided. Target applications of the dataset include: climate and ocean model evaluation; quantification of marine change and variability (including marine heatwaves); climate and ocean-atmosphere processes; and specific applications in ocean ecology, oceanography and geophysics

Satellite-based time-series of sea-surface temperature since 1980 for climate applications

A 42-year climate data record of global sea surface temperature (SST) covering 1980 to 2021 has been produced from satellite observations, with a high degree of independence from in situ measurements. Observations from twenty infrared and two microwave radiometers are used, and are adjusted for their differing times of day of measurement to avoid aliasing and ensure observational stability. A total of 1.5 × 1013 locations are processed, yielding 1.4 × 1012 SST observations deemed to be suitable for climate applications. The corresponding observation density varies from less than 1 km−2 yr−1 in 1980 to over 100 km−2 yr−1 after 2007. Data are provided at their native resolution, averaged on a global 0.05° latitude-longitude grid (single-sensor with gaps), and as a daily, merged, gap-free, SST analysis at 0.05°. The data include the satellite-based SSTs, the corresponding time-and-depth standardised estimates, their standard uncertainty and quality flags. Accuracy, spatial coverage and length of record are all improved relative to a previous version, and the timeseries is routinely extended in time using consistent methods

An Analysis of Proximal Volcanic Ash Emissions

Volcanic ash is a product of explosive volcanic eruptions and refers to those particles < 2 mm in size that are ejected from a volcano. Once erupted into the atmosphere, ash can be transported vast distances. Satellite remote sensing has provided us with the tools to map and monitor these plumes. These methods are restricted in viewing only that portion that is optically transparent, due to the requirement to observe the interaction of ground upwelling radiation with the plume. For this study, methods have been developed to map not only the opaque portion of the plume, but in higher spatial resolution than has been attempted previously. This has been done by using the unique emissivity spectra that are produced by materials in the thermal infrared portion of the electromagnetic spectrum. By using an end-member linear unmixing model with a spectral emissivity library of different volcanic ash types, opaque plume bearing pixels of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data have been processed. These end-member results have shown the applicability of this technique, with sampled volcanoes mapped with very low root mean square (RMS) errors. Volcanoes with unknown composition were also mapped with the library with varying degrees of success. Further expansion of the laboratory spectral library will allow more accurate assessment of these volcanoes. The collection of per pixel emissivity spectra of opaque plumes was also attempted using an experimental multispectral FLIR thermal camera apparatus, in order to collect data on the rising volcanic ash column. These data showed that the camera can produce emissivity spectra, however further calibration and correction is required in order to make a volcanological assessment of the eruptive products. Finally, a method of better tracking disconnected ash plumes is assessed. In cases where a plume is detected disconnected from or distally from the source volcano, the application of geostatistical methods to backward trajectory model data can yield the source location. This could then be used to helped point the ASTER sensor off-axis towards erupting targets, thus providing a greater quantity of data to be analyzed

Half a century of satellite remote sensing of sea-surface temperature

Sea-surface temperature (SST) was one of the first ocean variables to be studied from earth observation satellites. Pioneering images from infrared scanning radiometers revealed the complexity of the surface temperature fields, but these were derived from radiance measurements at orbital heights and included the effects of the intervening atmosphere. Corrections for the effects of the atmosphere to make quantitative estimates of the SST became possible when radiometers with multiple infrared channels were deployed in 1979. At the same time, imaging microwave radiometers with SST capabilities were also flown. Since then, SST has been derived from infrared and microwave radiometers on polar orbiting satellites and from infrared radiometers on geostationary spacecraft. As the performances of satellite radiometers and SST retrieval algorithms improved, accurate, global, high resolution, frequently sampled SST fields became fundamental to many research and operational activities. Here we provide an overview of the physics of the derivation of SST and the history of the development of satellite instruments over half a century. As demonstrated accuracies increased, they stimulated scientific research into the oceans, the coupled ocean-atmosphere system and the climate. We provide brief overviews of the development of some applications, including the feasibility of generating Climate Data Records. We summarize the important role of the Group for High Resolution SST (GHRSST) in providing a forum for scientists and operational practitioners to discuss problems and results, and to help coordinate activities world-wide, including alignment of data formatting and protocols and research. The challenges of burgeoning data volumes, data distribution and analysis have benefited from simultaneous progress in computing power, high capacity storage, and communications over the Internet, so we summarize the development and current capabilities of data archives. We conclude with an outlook of developments anticipated in the next decade or so

Land surface temperature and evapotranspiration estimation in the Amazon evergreen forests using remote sensing data

Amazonian tropical forests play a significant role in global water, carbon and energy cycles. Considering the relevance of this biome and the climate change projections which predict a hotter and drier climate for the region, the monitoring of the vegetation status of these forests becomes of significant importance. In this context, vegetation temperature and evapotranspiration (ET) can be considered as key variables. Vegetation temperature is directly linked with plant physiology. In addition, some studies have shown the existing relationship between this variable and the CO2 absorption capacity and biomass loss of these forests. Evapotranspiration resulting from the combined processes of transpiration and evaporation links the terrestrial water, carbon and surface energy exchanges of these forests. How this variable will response to the changing climate is critical to understand the stability of these forests. Satellite remote sensing is presented as a feasible means in order to provide accurate spatially-distributed estimates of these variables. Nevertheless, the use of satellite passive imagery for analysing this region still has some limitations being of special importance the proper cloud masking of the satellite data which becomes a difficult task due to the continuous cloud cover of the region. Under the light of the aforementioned issues, the present doctoral thesis is aimed at estimating the land surface temperature and evapotranspiration of the Amazonian tropical forests using remote sensing data. In addition, as cloud screening of satellite imagery is a critical step in the processing chain of the previous magnitudes and becomes of special importance for the study region this topic has also been included in this thesis. We have mainly focused on the use of data from the Moderate Resolution Imaging Spectroradiometer (MODIS) which is amongst major tools for studying this region. Regarding the cloud detection topic, the potential of supervised learning algorithms for cloud masking was studied in order to overcome the cloud contamination issue of the current satellite products. Models considered were: Gaussian Naïve Bayes (GNB), Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Random Forests (RF), Support Vector Machine (SVM) and Multilayer Perceptron (MLP). These algorithms are able to provide a continuous measure of cloud masking uncertainty (i.e. a probability estimate of each pixel belonging to clear and cloudy class) and therefore can be used under the light of a probabilistic approach. Reference dataset (a priori knowledge) requirement was satisfied by considering the collocation of Cloud Profiling Radar (CPR) and Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) observations with MODIS sensor. Model performance was tested using three independent datasets: 1) collocated CPR/CALIOP and MODIS data, 2) MODIS manually classified images and 3) in-situ ground data. For the case of satellite image and in-situ testing, results were additionally compared to current operative MYD35 (version 6.1) and Multi-Angle Implementation of the Atmospheric Correction (MAIAC) cloud masking algorithms. These results showed that machine learning algorithms were able to improve MODIS operative cloud masking performance over the region. MYD35 and MAIAC tended to underestimate and overestimate the cloud cover, respectively. Amongst the models considered, LDA stood out as the best candidate because of its maximum accuracy (difference in Kappa coefficient of 0.293/0.155 (MYD35 /MAIAC respectively)) and minimum computational associated. Regarding the estimation of land surface temperature (LST), the aim of this study was to generate specific LST products for the Amazonian tropical forests. This goal was accomplished by using a tuned split-window (SW) equation. Validation of the LST products was obtained by direct comparison between LST estimates as derived from the algorithms and two types of different LST observations: in-situ LST (T-based validation) and LST derived from the R-based method. In addition, LST algorithms were validated using independent simulated data. In-situ LST was retrieved from two infrared radiometers (SI-100 and IR-120) and a CNR4 net radiometer, situated at Tambopata test site (12.832 S, 62.282 W) in the Peruvian Amazon. Apart from this, current satellite LST products were also validated and compared to the tuned split-window. Although we have mainly focus on MODIS LST products which derive from three different LST algorithms: split-window, day and night (DN) and Temperature Emissivity Separation (TES), we have also considered the inclusion of the Visible Infrared Imaging Radiometer Suite (VIIRS) sensor. In addition, a first assessment of the Sea and Land Surface Temperature Radiometer (SLSTR) is presented. Validation was performed separately for daytime and nighttime conditions. For MODIS sensor, current LST products showed Root Mean Square Errors (RMSE) in LST estimations between 2 K and 3K for daytime and 1 K and 2 K for nighttime. In the best case (with a restrictive cloud screening) RMSE errors decrease to values below 2K and around 1 K, respectively. The proposed LST showed RMSE values of 1K to approximately 2 K and 0.7-1.5 K (below 1.5 K and below 1 K in the best case) for daytime and nighttime conditions, thus improving current LST MODIS products. This is also in agreement with the R-based validation results, which show a RMSE reduction of 0.7 K to 1.7 K in comparison to MODIS LST products. For the case of VIIRS sensor daytime conditions, VIIRS-TES algorithm provides the best performance with a difference of 0.2 K to around 0.3 K in RMSE regarding the split window algorithm (in the best case it reduces to 0.2 K). All VIIRS LST products considered have RMSE values between 2 K and 3K. At nighttime, however VIIRS-TES is not able to outperform the SW algorithm. A difference of 0.7 K to 0.8 K in RMSE is obtained. Contrary to MODIS and the SW LST products, VIIRS-TES tends to overestimate in-situ LST values. Regarding SLSTR sensor, the L2 product provides a better agreement with in-situ observations than the proposed algorithm (daytime difference in RMSE around 0.6 K and up 0.07 K at nighttime). In the estimation of the ET, we focused on the evaluation of four commonly used remote-sensing based ET models. These were: i) Priestley-Taylor Jet Propulsion Laboratory (PT-JPL), ii) Penman-Monteith MODIS operative parametrization (PM-Mu), iii) Surface Energy Balance System (SEBS), and iv) Satellite Application Facility on Land Surface Analysis (LSASAF). These models were forced using remote-sensing data from MODIS and two ancillary meteorological data sources: i) in-situ data extracted from Large-Scale Biosphere-Atmosphere Experiment (LBA) stations (scenario I), and ii) three reanalysis datasets (scenario II), including Modern-Era Retrospective analysis for Research and Application (MERRA-2), European Centre for Medium-range Weather Forecasts (ECMWF) Re-Analysis-Interim (ERA-Interim), and Global Land Assimilation System (GLDAS-2.1). Performance of algorithms under the two scenarios was validated using in-situ eddy-covariance measurements. For scenario I, PT-JPL provided the best agreement with in-situ ET observations (RMSE = 0.55 mm/day, R = 0.88). Neglecting water canopy evaporation resulted in an underestimation of ET measurements for LSASAF. SEBS performance was similar to that of PT-JPL, nevertheless SEBS estimates were limited by the continuous cloud cover of the region. A physically-based ET gap-filling method was used in order to alleviate this issue. PM-Mu also with a similar performance to PT-JPL tended to overestimate in-situ ET observations. For scenario II, quality assessment of reanalysis input data demonstrated that MERRA-2, ERA-Interim and GLDAS-2.1 contain biases that impact model performance. In particular, biases in radiation inputs were found the main responsible of the observed biases in ET estimates. For the region, MERRA-2 tends to overestimate daily net radiation and incoming solar radiation. ERA-Interim tends to underestimate both variables, and GLDAS-2.1 tends to overestimate daily radiation while underestimating incoming solar radiation. Discrepancies amongst these inputs resulted in large absolute deviations in spatial patterns (deviations greater than 500 mm/year) and temporal patterns