Next Article in Journal
Preliminary Evaluation of FY-3E Microwave Temperature Sounder Performance Based on Observation Minus Simulation
Previous Article in Journal
Reconstruction of Vegetation Index Time Series Based on Self-Weighting Function Fitting from Curve Features
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relative Merits of Optimal Estimation and Non-Linear Retrievals of Sea-Surface Temperature from MODIS

by
Malgorzata D. Szczodrak
* and
Peter J. Minnett
Department of Ocean Sciences, Rosenstiel School of Marine and Atmospheric Science, University of Miami, 4600 Rickenbacker Causeway, Miami, FL 33149, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(9), 2249; https://doi.org/10.3390/rs14092249
Submission received: 21 March 2022 / Revised: 28 April 2022 / Accepted: 5 May 2022 / Published: 7 May 2022
(This article belongs to the Section Ocean Remote Sensing)

Abstract

:
We compared the results of an Optimal Estimation (OE) based approach for the retrieval of the skin sea surface temperature (SSTskin) with those of the traditional non-linear sea surface temperature (NLSST) algorithm. The retrievals were from radiance measurements in two infrared channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) on the NASA satellite Aqua. The OE used a reduced state vector of SST and total column water vapor (TCWV). The SST and atmospheric profiles of temperature and humidity from ERA5 provided prior knowledge, and we made reasonable assumptions about the variance of these fields. An atmospheric radiative transfer model was used as the forward model to simulate the MODIS measurements. The performances of the retrieval approaches were assessed by comparison with in situ measurements. We found that the OESST reduces the satellite–in situ bias, but mostly for retrievals with an already small bias between in situ and the prior SST. The OE approach generally fails to improve the SST retrieval when that difference is large. In such cases, the NLSST often provides a better estimate of the SST than the OE. The OESST also underperforms NLSST in areas that include large horizontal SST gradients.

1. Introduction

Sea surface temperature (SST) is a critically important geophysical parameter in the Earth System. The Committee on Earth Observation Satellites (CEOS) has declared SST to be an Essential Climate Variable (ECV) that has a high impact on the requirements of the UNFCCC (United Nation Framework Convention on Climate Change) and the IPCC (Intergovernmental Panel on Climate Change) [1]. SST derived from satellite measurements offers the best source of global, repeated fields, along with accompanying uncertainty characteristics. For several decades now, SST has been routinely derived from radiances measured by a range of infrared (IR) or microwave (MW) satellite sensors [2]. The satellite-derived SST is an important element of the global observing system for numerical weather prediction (NWP) and for climate research. The requirements for the accuracy of SST are very demanding, as is the need for error analysis to fulfill the requirements for the generation of climate data records (CDRs) [3]. The concept of a CDR was formally introduced in a National Academy of Sciences report [4] as being “a data set designed to enable study and assessment of long-term climate change, with ‘long-term’ meaning year-to-year and decade-to-decade change. Climate research often involves the detection of small changes against a background of intense, short-term variations…. The production of CDRs requires repeated analysis and refinement of long-term data sets, usually from multiple data sources”. The generation of a CDR of global SST poses stringent requirements on the accuracy of derived SST: An absolute temperature uncertainty of 0.1 K and a decadal stability of ~0.04 K [5].
Over the years since satellite SST retrievals began [2], the most common approach to the retrieval of SST from satellite IR sensors employs algorithms based on the “split-window” technique [6,7,8] where satellite measurements are made in two IR channels typically centered near 11 and 12 μm wavelengths where the cloud-free atmospheric transmissivity is relatively high, but showing a wavelength dependence. The SST is derived by a slightly nonlinear combination of the brightness temperatures (BTs) in these two channels and some baseline sea surface temperature estimate (bsst). This is known as the non-linear SST algorithm (NLSST; ibid). For the Moderate Resolution Imaging Spectroradiometer (MODIS) which is flown on the NASA Aqua and Terra satellites [9] the NLSST has the form:
NLSST = c 1 + c 2   T 11 + c 3     ( T 11 T 12 )     bsst + c 4 sec ( satz ) ( T 11 T 12 ) + c 5 m i r r o r s i d e + c 6   satz + c 7     satz 2  
where T1 and T12 are brightness temperatures (BTs) in the 11 and 12 μm channels, bsst is the baseline SST, satz is the satellite zenith angle, mirrorside depends on which side of the two-sided, paddle-wheel scan mirror directs radiation to the radiometer, and c1 to c7 are coefficients [10]. Similar formulations exist for other multi-channel IR instruments including the Advanced Very High Resolution Radiometer (AVHRR) [11], the (Advanced) Along-Track Scanning Radiometer ((A)ATSR) [12,13], the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) [14], the Advanced Baseline Imager (ABI) [15], the Advanced Himawari imager (AHI) [16], and the Visible Infrared Imaging Radiometer Suite (VIIRS) [17]. The channel BT difference, T11-T12, is mainly sensitive to the water vapor content of the atmosphere. The coefficients c1 to c7 are derived by regression of the BTs to in situ measurements of SST (e.g., [8,11,18,19]), that are included in the Match-Up Data Bases (MUDBs) [11]. The coefficients may alternatively be derived by regression of channel BTs simulated using radiative transfer modeling and the value of the SST specified as the boundary condition of the model (e.g., [12,20,21]).
Since the NLSST algorithms rely on coefficients that are generally obtained by least-squares regression to in situ data (moored buoy or drifter measurements) they produce good estimates of SST in conditions that are represented by the conditions implicit in the MUDB records used to generate the coefficients. However, when the sea or atmospheric conditions depart from the average state, the NLSST estimates can be quite erroneous. The existence of regional and seasonal biases in the NLSST is well documented [11,22,23]. A good example of this bias is the retrieval over the Northern Atlantic Ocean during Saharan Air Layer (SAL) outbreaks [24]. Typically, the atmosphere in this region is moist and such conditions are dominantly represented in the MUDB used to derive the coefficients of the NLSST formulation. During SAL events, moist air is replaced by a layer of very dry air of continental origin, often including dust aerosols from the Saharan region, and the established coefficients are not representative of these conditions, resulting in larger retrieval errors [25].
One way of reducing the NLSST bias is to select from multiple sets of coefficients those that are specific to different seasons and geographic regions [11,17,19]. Another way is to model explicitly the effects of anomalous atmospheric conditions, including SAL and aerosol effects, using radiative transfer calculations and subtract the modeled bias from the retrieval; this is typically used for retrievals under SAL conditions [25,26,27,28].
In general, in the ocean–atmosphere system, there are fundamentally two leading modes of variability represented in the measurements of the two IR channels used in the NLSST formulation. These are the SST itself and the total atmospheric moisture content (total column water vapor, TCWV) and their variability is captured by the retrieval coefficients derived by linear regression, meaning the NLSST estimates cannot be improved without introducing additional information into the retrieval system. This is accomplished when regional or seasonal NLSST coefficients are used in the retrieval. Prior information on the mean and variability of the SST and atmosphere of a region or a season is then effectively embedded in the coefficients.
More recently methods based on Bayesian approach to inverse problems [29] have been proposed and used to retrieve SST from multichannel IR measurements [30]. The term often used to refer to this approach is Optimal Estimation (OE). In the OE approach, prior knowledge of the state of the ocean and the atmosphere at the time and place of each satellite measurement is used to derive the SST estimates (OESST) and as such could be viewed as an additional level of the coefficient tuning of the NLSST approach.
The OE retrieval of SST has gained popularity in recent years with the promise of improving on the biases inherent in the NLSST approach. A very good agreement between the OE estimates and in situ buoy measurement of average bias of −0.06 K and −0.01 K have been achieved by Merchant et al. (2008) for AVHRR on the Metop-A satellite and by Merchant et al. [31,32] for SEVIRI on Meteosat 9. The statistics of these biases are indeed good, being several times lower than those for corresponding SST derived by NLSST algorithms using global coefficients. Additionally, these retrievals specifically exclude conditions where the OE approach might not work well: Only nighttime measurements were considered, and measurements affected by SAL were excluded.
With an expectation of OE generally outperforming the NLSST retrievals, we focus on comparisons between the two. The OE approach comes with its own set of uncertainties and sources of error which, if not addressed, undermine the expected gain in accuracy or even usefulness of the retrieval. Koner et al. ([33]) proposed an alternative SST retrieval method: The Modified Total Least Squares (MTLS) and found that OE retrieval was worse than the MTLS solution, and sometimes worse than the a priori. They argued that OE is unable to produce the best retrieval for a fairly linear and moderately ill-conditioned problem of SST retrieval, relating it to the estimation of errors (which the MTLS avoids) being difficult in practice. Indeed, the error characteristics are typically assessed by “expert judgment” based on the best available knowledge [32]. Additionally, since the OE requires unbiased priors, care must be taken to remove bias from the priors before the retrieval. In order to remove the degree of arbitrariness from the problem, Merchant et al. (ibid) proposed an OE-based approach to the derivation of the biases and covariance parameters and in application to SST retrieval from SEVIRI measurements achieved improved results over an OE application with non-optimized parameters.
In this work we are not seeking optimal values of the covariance parameters but focus on examining how possible values of these parameters influence the usefulness of the OE method to the SST retrieval from MODIS. We show that in the case of MODIS SST the main factor contributing to the poor OE retrieval is the bias of the prior SST, while all other combined biases play a much smaller role. To assess the relative consequences of factors other than the bias of prior SST on retrieved OESST, we first conducted the OE retrieval with prior SSTs from the temperature measured in situ corrected by the mean skin effect. This eliminates the prior SST bias from the retrieval problem. Since we consider the matching in situ measurements of SST as the reference temperature, any difference between the OESST in this case and the in situ SST is attributable to all the other biases. We found the unresolved biases account on average for about 0.02 K of the OE and in situ temperature difference which could be viewed as a limit on the accuracy of OESST. We then evaluate the accuracy of the OESST with prior SST from ERA5 reanalysis.
We will also address the question: By how much is our knowledge of the SST improved in the OE retrieval in comparison with the a priori knowledge, and how does it compare to the NLSST retrieval?
We briefly review the basic concepts of the Bayesian retrieval in Section 2 and its application to a linear model in Section 3. In Section 4 we describe the application of the OE approach to the retrieval of SST from the MODIS data set used in this study. Section 5 introduces measures used to assess the information gain of the retrieval. In Section 6 we discuss the results of the OE retrieval and identify situations that impair the performance of OE. We compare OE and the NLSST retrievals and suggest when one of these methods is more likely to provide a better estimate of SST than the other in Section 7. Summary and conclusions are given in Section 8.

2. Bayesian Inversion Concepts

Imagine that we have a system that can be described by state vector x, and in order to gain some knowledge about this system, we perform a measurement which results in a vector y. The act of measurement maps the state space into the measurement space. In the presence of a measurement error, the mapping is not exact, i.e., a point in the state space maps into a region in the measurement space that is determined by the probability density function (pdf) of the measurement error ϵ. This mapping can be represented by a forward function F(x) or, in practice, often by a forward model F, and written as:
      y = F ( x ) + ϵ      
Inversely, a measurement of y does not represent a single point but rather a region in the state space that can be described by some pdf. If we have some prior (imperfect) knowledge of the state before the measurement, then that information can be used to constrain the region representing the solution. This prior knowledge can often also be described by a pdf. Bayes’ Theorem shows how the pdf of an imperfect measurement maps back to the state space and combines with the pdf of the imperfect prior knowledge to improve our understanding of the system,
P ( x | y ) = P ( y | x ) P ( x ) / P ( y )  
where P(x|y) is the pdf of the state being described by vector x provided the measurement is y, P(x) is the pdf of the prior state, P(y|x) is the conditional pdf of obtaining measurement y provided the state vector is x, and P(y) is pdf of obtaining a measurement y. P(y) can be obtained by integrating P(y|x)P(x) over all states x, but is only a normalizing factor and is typically not necessary.
Conceptually, the Bayesian approach can be summarized as follows:
i.
we have some knowledge of the state before the measurement is even made and this prior knowledge is expressed by a prior pdf of the state variables
ii.
we have a forward model that maps the state variable into the measurement space
iii.
we know the pdf of the measurement errors
iv.
we use Equation (3) to calculate posterior pdf by augmenting the prior pdf of the state vector with the measurement.
The result of the Bayesian inversion is thus a pdf of all possible states. This ‘solution’ informs us how the measurement improves our knowledge of the state but to obtain a retrieval that can be more readily compared with other methods, one particular state needs to be selected from all the possible posterior states. This would be an optimal state according to some metric, and often is selected as the most probable state or the expected value state of the posterior pdf.

3. Linear Forward Model and Gaussian PDFs

The distribution of error in many real-world measurements adheres fairly closely to the Gaussian distribution and since the Gaussian formulation is mathematically expedient it is often convenient to use it as an approximation for both the measurement error pdf and the prior pdf.
The forward model F can typically be linearized around some reference state and then Equation (2) becomes
y = K x + ϵ  
where K is a matrix of partial derivatives Kij = ∂Fi/∂xj (a Jacobian), sometimes referred to as the kernel.
In the case of a linear F and Gaussian measurement error and prior pdfs, the posterior pdf is also Gaussian with expected value and covariance given by
  x ^ = x a + ( K T S ϵ 1 K + S a 1 ) 1 K T S ϵ 1 ( y F ( x a ) )
and
      S ^ 1 = ( K T S ϵ 1 K + S a 1 )  
where xa and Sa are the expected value and the error covariance matrix of the prior pdf and the Sϵ is the covariance matrix of the measurement error.
With the above assumptions, Equation (5) indicates that the retrieved state is a sum of the prior state and a correction which is directly related to the difference between the measurement and the model simulated values for the a priori state. In the context of SST retrieval from satellite IR instruments, the measurement is the radiance measured in the spectral channels of the satellite instrument while the model measurement is the simulated value of radiance in those channels calculated by radiative transfer for a selected a priori state.

4. Gain and Averaging Representation

Equation (5) represents a regularized least square problem with Sa−1 acting as the regularization matrix. The regularization is necessary when the inverse problem is underdetermined or ill-posed, i.e., the matrix KTSϵ−1 K has a high condition number, as it often is the case in remote sensing inversion problems. In OE the strength of the regularization and thus the ability to obtain a stable solution is determined by the statistics of the a priori state in relation to the matrix KTSϵ−1 K.
Equation (5) is often written as
  x ^ = x a + G ( y F ( x a ) )    
where G, the gain matrix, gives the sensitivity of the retrieval to the observations: G =   x ^ / y . Any errors in the measurement are multiplied by the gain and propagated to the solution. Due to the regularization the retrieval obtained by Equation (5) varies from the state that is most consistent with the pure measurement. The retrieved OE state results from the combination of the measurement and the priori information and the ratio at which they are combined depends on the strength of the regularization. In another transformation Equation (5) can be written as
        x ^ = A x a + ( I n A ) x a    
where matrix A = GK is the averaging kernel matrix and represents the sensitivity of the retrieved state to the true state ( A = x ^ / x ) and In is an identity matrix of size n, where n is the dimension of the state. The rows of the averaging kernel represent the weighting functions, which determine to what degree the result depends on the components of the true state.
It should be instructive to examine the regularization and the matrices G and A for any retrieval problem as they will inform one about the stability of the OE solution and reveal how much the measurement contributes to the improvement of our knowledge of the state vector as given by the a priori, or if a given measurement contributes to the solution at all. What is more, this can be evaluated even before the measurements or retrievals are made. For example, since the gain matrix is determined by the relative magnitude of KTSϵ−1 K and Sa1 if the measurement error is large so that KTSϵ−1 K << Sa1, then the gain will be close to zero and x ^     x a   . Inversely, if the measurement is so precise that KTSϵ−1 K >> Sa1, the a priori information becomes irrelevant. Correspondingly, in terms of the averaging matrix A, if A is close to the identity matrix the retrieval will be mostly determined by the measurement.
What do we know about K, Sϵ−1 and Sa1 for the problem of SST retrieval from MODIS 11 and 12 μm channel radiances, and under what condition can we expect the OE retrieval to be worthwhile? Answering these questions is the subject of the next section.

5. Application to SST Retrievals

5.1. Satellite Dataset

We apply the OE retrieval and the above concepts to a subset of MODIS measurements from the Aqua satellite with the intention of establishing the relative merits of the OE approach compared to the NLSST. The MODIS data used in this study are a subset of the MUDB MODIS R2019 data set, where R2019 is the identifier for the set of cloud-screening and atmospheric correction algorithms, sometimes referred to as Version 6 or Collection 6 [10]. Our subset spans two years, 2015 and 2016, from the area of the northern Atlantic Ocean and the Mediterranean Sea selected because of the relatively large number of match-ups; a wide range of SSTs; extreme values of horizontal SST gradients, ranging from the edges of the Gulf Stream to central gyres; and the high variability of atmospheric conditions including those influenced by different continental airmasses advected over the ocean. In addition to in situ SST measurements matched to the SSTskin derived with the NLSST algorithm, the MUDB contains the 11 and 12 μm radiances measured by the instrument in channels 31 and 32, and a number of other fields relevant to the SST retrieval. Each entry in the MUDB includes assessment of the confidence in the retrieved value expressed by the quality flag (Q). The best match-ups have Q = 0 and the slightly worse Q = 1 (these are match-ups that would otherwise be assigned Q = 0 but for the most part correspond to large satellite zenith angles meaning long atmospheric path lengths). Quality flags > 1 typically indicate pixels likely contaminated by radiance from clouds or dense aerosols and are not considered reliable. Here we only consider match-ups with Q = 0 or 1. For the purpose of this study we use a random sample of 30,595 Q = 0 data points (out of a random sample of 100,000 Q = 0 data points, 30,596 fall into the geographical area under consideration) and a random sample of 15,749 Q = 1 data points (out of 50,000).

5.2. Prior Knowledge

Our experience with the NLSST algorithm tells us that the radiances measured in MODIS channels 31 and 32 are most sensitive to the SST and the channel radiance difference is most sensitive to the TCWV [11,19]. With an ideal (error free) measurement of just two quantities one could hope to retrieve at most two components of the state vector, so we limit the representation of the state to a reduced state vector x = [SST, TCWV].
Our prior knowledge is the SST and TCWV fields taken from the European Center for Medium-Range Weather Forecasts (ECMWF) ERA 5 dataset, the fifth generation of ECMWF atmospheric reanalyses of the global climate [34]. Thus, we assume that the prior pdf of x is Gaussian with the mean [SST0, TCWV0] given by ERA5 and some prior variance that we will have to specify.
There are two relevant sea surface temperature fields in the ERA5 data set. First, there is the foundation SST (“sst” in ERA5) defined as temperature at a sufficient depth that there is no diurnal heating and which is taken from the Operational Sea-Surface Temperature and Sea Ice Analysis (OSTIA) analysis [35,36]. Under moderate to high wind speeds, and at night, this SST should best represent the in situ buoy measurements that are also made at depth (typically 20 cm) but it differs from the SSTskin obtained by satellite measurements that are sensitive to the temperature of the top few tens of micrometers of the ocean surface. The average offset has been estimated to be −0.17 K [37]. Thus, we use this value with the ERA5 SST to obtain the prior SSTskin. During daytime, the foundation temperature might not represent the subskin ocean temperature since diurnal heating which can reach several degrees or more [38,39,40] is not included. Using the foundation temperature as a prior in OE, even if corrected for the skin effect, will then necessarily lead to a biased OESST, especially for the daytime retrievals. A solution to this problem would be to correct the foundation temperature for the diurnal signal. To evaluate the diurnal correction, we drew random samples of 50,000 data points for each Q = 0 and Q = 1 subsets of the MUDB and examined the differences between the ERA5 SST and the in situ buoy measurements SSTb (both adjusted to skin temperature using the −0.17 K offset) separating nighttime and daytime measurements. In situ measurements from drifting or moored buoys and from ship mounted radiometers serve to validate the satellite retrieved SST and ultimately the results of the OESST retrieval will be assessed by comparison to these in situ measurements. If there is a bias between the prior SST and the in situ measurements, the OE results also will be biased so it is important to remove the bias before the retrieval.
One would expect the diurnal heating applied to the foundation SST to depend on time of the day or correspondingly on solar zenith angle, but this relationship turns out to be quite weak for the MODIS AQUA dataset. This could be due to the sun-synchronous orbit of the AQUA satellite leading to the MODIS data being acquired over a narrow range of local times (around 01:30 and 13:30 at the Equator). As a result, we will use a simple nighttime and daytime offset to account for the missing diurnal signal of the ERA5 SST foundation temperature. The offsets are calculated as an average of the robust mean differences between ERA5 SST and the in situ measurements of the two random samples drawn from the MUDB set. The nighttime correction to the ERA5 SST is −0.015 K and the daytime correction is −0.220 K.
The ERA5 data also contain a skin temperature variable (“skt”), which we denote here as SKT to differentiate it from the foundation SST. This product is based on a model that describes the difference between skin and foundation temperature as a combination of the cool skin effects and diurnal warming effects [40,41]. Most previous OESST retrieval studies and implementations are based on using the foundation temperature reduced by skin effect as the prior SST. Merchant et al. [30] use a correction for the skin effect of −0.2 K; here, as given above, we follow Donlon et al. [37] and use −0.17 K. In principle, the SKT should be a better prior for the OE retrieval so here we examine the utility of both of these fields.
For the purpose of the forward radiative transfer modeling, we also obtain from the ERA5 the atmospheric profiles of the temperature, water vapor, and ozone. The spatial resolution of the data is 0.4 × 0.4 degree with a temporal resolution of 6 h for the atmospheric profiles and SKT, but daily for the foundation temperature. ERA5 data can be downloaded through the Copernicus Climate Data Store (CDS) infrastructure. More details about the ERA5 SST product can be found in Hirahara et al. (2016) [42].

5.3. Forward Model

As the forward model we use the Radiative Transfer for TOVS (RTTOV) v12.1 [43,44], where TOVS means the TIROS Operational Vertical Sounder, and TIROS means Television InfraRed Observation Satellite. While initially developed to simulate the measurements of TOVS, it is now used to simulate the measurements of many satellite radiometers, including MODIS on Aqua. For every match-up in the MODIS MUDB we find corresponding ECMWF data by interpolating the model fields to the time and position of the matched-up measurement. Using RTTOV we compute MODIS channel 31 and 32 radiance measurements as well as the Jacobian matrix K. The forward model is of course only an approximation and thus is an additional source of error.

5.3.1. Measurement Error Covariance Sϵ

We followed Merchant et al. [30] in assuming that the measurements in channels 31 and 32 are independent and include the forward modeling errors (also independent) in the combined covariance matrix. If ξ = ξ31 = ξ32 is the uncertainty of radiance measurement in MODIS channels 31 and 32 and ξm is the forward model error in these two channels, then
    S ϵ = [ e 31 2 0 0 e 32 2 ]      
where
          e 31 2 = e 32 2 = ξ 2 + ξ m 2 sec ( θ ) 2    
and θ is the satellite zenith angle. The MODIS instrument specifications list the NeDT in IR channels at about 0.05 K which, for wavelengths of 11–12 μm, corresponds to 0.007 W/(m2 sr μm) in radiance units which we take as the estimate of ξ. Merchant et al. [30] estimate ξm to be ~0.15 K for the Advanced Very High Resolution Radiometer (AVHRR) on the first Meteorological Operational satellite (Metop-A) of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). Our comparison of RTTOV simulations of MODIS channel radiances with ERA5 atmospheric profiles and SST taken from in situ buoy measurements with MODIS measurements for a sample of 29,515 quality 0 matchups yields a root mean square of the difference of 0.02 W/m2-μm-sr which we take to represent the ξm.

5.3.2. Prior Covariance Sa

While the ERA5 SST and TCWV provide the mean values of the prior SST and TCWV distributions, the prior error covariance must still be specified. While not necessarily true, the SST and TCWV are typically considered independent and thus the error covariance matrix Sa can be represented by a diagonal matrix with variance of prior SST and TWC along the diagonal. Variance, or standard error of either SST or TCWV are not part of the ERA5 data so we need to estimate them separately. As stated above, the ERA5 SST is based on the Operational Sea Surface Temperature and Ice Analysis (OSTIA). OSTIA uses satellite SST data provided by international agencies via the Group for High Resolution SST (GHRSST) to produce daily global SST fields. The performance of OSTIA is assessed by comparison with in situ SST measurements. The average standard deviation of the differences between the OSTIA and in situ SST from Argo floats is ~0.5 K for the area of the Northern Atlantic Ocean although it varies regionally and can be much larger in areas of high SST variability [36]. For the subset of the MUBD dataset considered here we calculated the robust standard deviation of the difference between the ERA5 SST and the in situ temperature measured by buoys and obtained values ~0.5 K, in agreement with the evaluation of Good et al., (2020) [36]. We use robust statistics to avoid contamination by large outliers as these are almost certainly spurious. We take this standard deviation to represent the covariance of the prior SST, eSST2 = 0.52 K2. This value is in line with values used by others, e.g., [30].
The error in the prior TCWV is harder to estimate. Chen and Liu [45] compared TCWV derived from ECMWF with radiosonde measurements over the period of 1979 to 2014 and with GPS and Microwave Satellite measurements over the years 2000 to 2014, and 1992 to 2014, respectively. They found RMS differences between the ECMWF and the three types of measurements in range between 7% to 20% in global average with the smallest values for the satellite microwave measurements. Merchant et al. [30] used numerical weather prediction (NWP) forecast fields from the ARPEGE (Action de Recherche Petite Echelle Grande Echelle) operational model of Météo-France and assumed eTCWV = 0.25 × TCWV.
We settle for eTCWV = 0.20 × TCWV after Chen and Liu [45] but we will explore how different values of eTCWV and eSST affect the retrieval. As stated before, we assume the SST and TCWV to be independent and thus the prior covariance matrix is:
  S a = [ e S S T 2 0 0 e T C W V 2 ]    

6. Information Content Aspects of the OE Retrieval

6.1. Degrees of Freedom

In the case of two-channel measurements in the SST retrieval, if there were no measurement error, we should be able to determine exactly two independent variables: SST and TCWV. This is a problem with two degrees of freedom. However, given measurement error, the retrieved values will have some uncertainty associated with them and thus the effective number of degrees of freedom will be less than two. This can be seen as a component of the state whose variability is smaller than the measurement error and cannot therefore be distinguished from other states. Often the variability of the state components cannot be directly compared with the measurement error as they are in different spaces, so some transformation is necessary. The aim is to compare the measurement covariance with the prior covariance of the forward model simulated measurement. The effective number of degrees of freedom is the number of singular values of matrix
    K ˜ = S ϵ 1 / 2 K S a 1 / 2  
which are greater than about 1.
The number of degrees of freedom can further be partitioned into those related to the signal (ds) and those related to measurement noise (dn). The ds and dn can be expressed in terms of the singular values of the matrix K ˜ as
  d s = i = 1 m λ i 2 1 + λ i 2  
and
d n = i = 1 m 1 1 + λ i 2    
where λi are the singular values of matri   K ˜ x, and m is the number of independent measurements.

6.2. Shannon Information Content

The Shannon information content is a concept in information theory that relates to the entropy of the probability density function. Qualitatively, it can be thought of as a factor by which knowledge of a quantity is improved by making a measurement. In our case, the question is whether the satellite IR measurement improves our knowledge of SST over our prior knowledge. Formally, if P1(x) is a pdf that describes knowledge before the measurement and P2(x) describes knowledge after the measurement, and if S is entropy, then the information content of the measurement is the reduction in entropy, H:
H = S ( P 1 ) S ( P 2 )    
For Gaussian pdfs
S ( P ) = ln ( | S | ) 2 + c o n s t    
where |S| is the determinant of the covariance matrix. And so
  H = ln ( | S ^ 1 S a | ) 2   .
The Shannon information content can also be expressed in terms of the singular values of matrix K ˜ ,
H = 1 / 2 i = 1 m ln ( 1 + λ i 2 )
A much more detailed explanation and discussion of the concepts summarized above is given by Rodgers [29].

6.3. Retrieval Potential

The covariance matrices Sa and Sϵ together with the Jacobian matrix K determine the retrieval process (i.e., how much weight is assigned to prior state vs. the measurement through Equation (8), so once these covariance matrices are known and the Jacobians are calculated, we can glimpse what we can expect from the OE in the particular retrieval problem we are considering. While the interpretation of scalar weight is straightforward it is less so for matrix weights appearing in Equation (8). However, since we are dealing with just two variables it is easy to write out the matrix equation explicitly, thus
  ( S S T O E T C W V O E ) = ( A 11 S S T + ( 1 A 11 ) S S T a + A 12 ( T C W V T C W V a ) A 22 T C W V + ( 1 A 22 ) T C W V a + A 21 ( S S T S S T a )     )  
and since we are interested primarily in SST we need to focus on just two elements of the A matrix A11 and A12.
The values of Sa and Sϵ laid out in Section 5 represent what we consider the most likely values of these parameters for our case of SST retrieval from MODIS. We will assess the relative contributions of the true and prior state to the OE retrieval for this case by examining the elements of matrix A, but we will also consider how these contributions change if the conditions of the retrieval change. In particular, we will consider two additional scenarios: (a) Poorer/better model or measurement scenario when either the radiative transfer model used or the instrument is less/more accurate which will affect the matrix Sϵ; and (b) when the a priori state covariance is smaller/larger than in our case (different a priori fields could be used in the retrieval).
Figure 1 shows how the two elements of the averaging matrix respond to changes in Sϵ and Sa.
Sϵ will change if either ξm or ξ31,32 change. We will only consider changes in ξm as (1) ξ31,32 come from instrument specification and for given instrument should not be a subject to change; and (2) the effect of changes in either ξ31,32 or ξm on matrix A will follow a similar pattern. On the other hand, ξm could change if, for example, a different RTM were to be used to compute the simulated radiances. Figure 1a,b show how A11 and A12 depend on ξm, Figure 1c,d show dependence on eSST, and Figure 1e,f on eTCWV. The red dots in the plots indicate values for the parameters we will use in this study as described in Section 5. For these values, as can be seen in Figure 1, the results of the OE retrieval will be dominated by the true state with only a small contribution coming from the prior state. This corresponds to only a weak regularization in Equation (5) and the OE solution being quite sensitive to measurement errors.
In Figure 2 we consider a range of different retrieval scenarios. Figure 2a,b show A11 and A12 for a range of ξm and eSST values, Figure 1c shows logarithm of the condition number of the matrix that is to be inverted in Equation (5), and Figure 1d shows the standard deviation of the retrieved SST distribution. The dashed lines are isolines of eSST; the red line and the red dot are those from Figure 1 representing the values accepted as most likely in our case. As expected, for any value of eSST, an increase in ξm (i.e., a less accurate RTM model) will lead to an increased contribution to the OESST from the prior SST and as the prior SST becomes more accurate, the measurement becomes less relevant. The contribution from the cross term initially increases as the model or measurement becomes less accurate but then drops and the more accurate prior SST lowers the impact of this term on the final retrieval. This contribution will also be small if the prior TCWV is close to the true value independent of the ξm or eSST. Larger values of ξm or eSST also lead to a stronger regularization of the retrieval problem as expressed by the condition number in Figure 2c. The standard variation of the posterior SST distribution is for small eSST independent of the ξm as expected, since the solution is almost entirely determined by the prior state, it then increases with ξm as eSST increases.
The red dots in Figure 2 indicate that for our particular retrieval problem we expect the solution for SST to be driven by measurements with only a weak influence of the prior state. Where the prior state contributes is to the retrieval of the TCWV as can be seen in Figure 3 which shows the weights A22 and A21.

7. The Retrieval Process and Results

7.1. Retrieval with ERA5 “sst” as Prior SST

Figure 4a,b shows histograms of (a) nighttime, and (b) daytime OESST–SSTb for ERA5 SST0 prior with the nighttime or daytime correction. The average statistics of the OESST–SSTb comparison given in the plots (figures in the brackets are robust values) are very good but they are not very different from the SST0–SSTb (see Figure 6a,b in Section 7.2). Figure 4c,d shows corresponding night and daytime histograms for the NLSST–SSTb. The main motivation behind developing new retrieval methods is to improve the accuracy of the retrievals especially for cases where the established methods perform poorly and there is an improvement in the robust statistics for our sample when going from NL to OESST although this is very small for daytime matchups. However, in comparison with NLSST most of the improvement appears to take place near the center of the histograms, i.e, the accuracy of OE is improved mostly for cases when the NL atmospheric correction is already effective. But the tails of the histograms appear to extend further for the OE than the NLSST.
Figure 5 illustrates the fraction of the dataset with absolute value of the SST–SSTb above a given threshold for OESST, NLSST, and also for ERA5 and ERA5 with the night/day correction. For both the (a) nighttime and (b) daytime retrieval, OE does the best job for absolute differences between the retrieval and the in situ temperature below about 0.8 K; however, it does only slightly better than the corrected ERA5 prior. The fraction of retrievals with differences larger than 0.8 K is smallest for the NLSST. In total, about 76% of the OESST data fall outside of the 0.1 K range of CDR target and about 55% fall outside the 0.2 K accuracy of the drifter buoy measurements. The corresponding percentages for the NLSST retrievals are 82% and 65% so with respect to this metric there is an advantage in using OESST rather than the NLSST.

7.2. Using ERA5 “skt” as Prior SST

The ERA5 Skin Temperature variable might seem a more appropriate prior SST than the foundation temperature and it should be able to be used directly without the need for corrections for the skin effect or for diurnal signals (although with the SKT only available four times a day interpolation in time might still be necessary). The SKT, however, is a modeled product that comes with its additional uncertainties [46,47]. Previous comparisons of SKT with skin temperature measured by the shipborne Marine-Atmospheric Emitted Radiance Interferometer (M-AERI, [48]) showed average differences of −0.237 K for nighttime and −0.172 K for daytime measurements [47]. These values are not only larger than the CDR requirement for the SST accuracy, but taking into account that the estimated accuracy of the M-AERI SSTskin is better than 0.05 K indicates a real departure of the ERA5 skin temperature model from the measurements. In Figure 6 we compare the nighttime and daytime differences between ERA5 SST (corrected for the skin effect and with added diurnal correction) and ERA5 SKT and the SSTb in our match-up dataset. Both nighttime and daytime SKT–SSTb differences are larger than SST–SSTb differences which makes the ERA5 SKT a poorer candidate for the SST0 in the OESST retrievals than the ERA5 SST, perhaps a surprising result.

7.3. Geographical Distribution of ERA5 SST0–SSTb

We focus on the large differences between the ERA5 SST0 and the buoy measurements since these contribute most to the OESST–SSTb offset, and so are the retrievals in most need of improvement. Figure 7 shows a map with the locations of match-up points colored by the value of ERA5 SST0–SSTb over contours of SST0 for 1 January 2015, the first day of our dataset. (While the contours of SST0 are for one time instance, the SST0–SSTb are for the entire subset of the matchups.) Panel (a) of Figure 7 shows locations of SST0–SSTb > 1 K differences, and panel (b) the locations of match-ups with negative SST0–SSTb < −1 K. These data points correspond to the tails of the distributions in Figure 4. For comparison, match-ups with |SST0–SSTb| < 0.2 K are shown in panel (c). Different symbols correspond to different types of in situ measurements: Moored buoys, drifting buoys, and ship radiometers.
Figure 7 reveals that the majority of match-ups with large |SST0–SSTb| are found in ocean areas with strong horizontal gradients in SST. Since ERA5 SST represents the temperature of an area of 0.4 × 0.4 degrees in latitude and longitude, it might not represent well any individual buoy, especially those in areas with strong horizontal gradients. In such conditions the OESST might not provide good estimates of the SST. In Table 1 we compare the statistics of OESST–SSTb with NLSST–SSTb in bins of SST0–SSTb for night (a) and daytime (b) Q = 0 measurements to see if and under which conditions the NLSST can provide more accurate estimates of the SST than the OE. The StDev for OE retrievals is calculated as the square root of the OESST variance given by Equation (6) with eSST = 0.5 K. The ΔNLSST is an estimate of NLSST measurement error calculated through the propagation of the uncertainties in the MODIS channel measurements, e31 and e32, through the NLSST equation given in Equation (1). The graphical representation of the mean, RMS, and standard deviation from Table 1 daytime data is shown in Figure 8 (nighttime data follow a very similar pattern and are not shown).
As can be seen in Table 1 and Figure 8 in terms of mean value, OE retrievals are more accurate than those of the NLSST provided the difference between the prior SST and buoy measurements is in the range ± 0.1 K for daytime, with the exception of nighttime negative differences where OE slightly outperforms NLSST up to a difference of −0.25 K and daytime positive differences where the NLSST is better starting bin 0.05 to 0.1 K. The RMS values are consistently smaller for the OE up to ± 1 K difference between SST0 and SSTb. Additionally, the mean value NLSST–SSTb appears to be largely independent of the SST0–SSTb and although RMS does increase with increasing SST0–SSTb, it does so at a slower rate than OESST–SSTb. These results indicate that in areas of the ocean with large SST variability or where strong horizontal SST gradients are prevalent, the OE approach will not necessarily improve on the NLSST if the a priori SST field does not correctly capture the variability. Unlike the diurnal bias, this spatial sampling bias is more difficult to correct, other than improving the spatial resolution of the prior SST0.

7.4. Information Measures

As outlined in Section 5 we calculated the number of degrees of freedom ds and the Shannon information content H for the retrieval with ERA5 SST0 with the diurnal correction and with eSST = 0.5 K. The average ds is 1.28 with standard deviation of 0.25 indicating that in practice only one element of the state vector can be retrieved.
The average Shannon information content of the retrieval is 4.69 with standard deviation of 0.47 for SST but just 0.30 and 0.30 for TCWV. So, while the measurement does improve our knowledge of the SST with respect to the a priori, it does not contribute significantly to our knowledge of the TCWV.

8. Comparison with NLSST

Table 1 demonstrates that OESST can be quite erroneous compared to buoy measurements if the prior SST0 is similarly biased, and the NLSST can provide more accurate SST estimates in such situations. One of the main reasons for considering the benefit of switching from NLSST to OESST is the prospect of improving on the biases in the NLSST that are due to the coefficients not adequately representing the atmospheric conditions during the measurement. Table 2 presents the improvement of OESST over NLSST in bins of NLSST–SSTb similar to Table 1 for (a) night and (b) daytime Q = 0 measurements. Mean, RMS, and standard deviation values from Table 2b are shown in Figure 9.
In most cases the mean value of |OESST-SSTb| is indeed smaller than |NLSST-SSTb| which could be considered an improvement as, on average, the OESST retrievals appear to agree better with the in situ measurements than those of the NLSST algorithm. However, in the majority of cases, SST0 agrees with the in situ SST even better than do the OESST retrievals which puts in question the need of performing the OE retrieval where SST0 values are available and understood to be accurate. The smaller mean values of |OESST–SSTb| (and |SST0–SSTb|) are accompanied by large values of the RMS error indicating large scatter about the mean. The RMS errors are typically significantly smaller for the NLSST–SSTb for all bins up to [−0.5, −0.25] K on the negative side and up to the last bin of [1.0, 5.0] K for the positive differences daytime data, and two last bins for the positive differences nighttime data.
Figure 10 shows the geographical distribution of (a) positive, and (b) negative NLSST–SSTb greater or more negative than ±1 K, respectively. It appears that the retrievals with large positive differences are distributed similarly to the large SST0–SSTb differences, while the large negative differences are distributed quite uniformly throughout the whole domain. This might indicate different causes for the positive and negative biases in the NLSST retrievals: Possibly cloud contamination causing the negative bias and land influence on the coastal atmospheres leading to a positive bias.
As shown in Section 6, the OE cannot work well for situations with large SST0–SSTb. However, since SSTb is not available for all satellite measurements, one cannot rely on this metric to determine if the OE approach is beneficial in any particular case. Figure 11 shows the dependence of OESST–SSTb and NLSST–SSTb on NLSST–SST0. It reveals a strong, almost linear relationship between OESST–SSTb and NLSST–SST0 with a correlation coefficient of −0.64 indicating that NLSST–SST0 might be used to determine whether the OE is likely to improve on the NLSST retrieval. There is also a correlation between NLSST–SSTb and NLSST–SST0 but it is slightly weaker with a correlation coefficient of 0.60.

9. Discussion and Conclusions

New approaches to the retrieval of SST from satellite measurements are introduced to improve on the existing algorithms in terms of increasing the accuracy of the retrievals, eliminating known biases and/or expanding coverage. The main promise of the OE approach is that it would eliminate, or at least reduce, the bias of the NLSST retrievals due to the inability of the NLSST atmospheric correction algorithm to compensate for the effects of atypical atmospheric conditions. However, the OE method comes with its own set of problems; a thorough review is given by von Clarmann et al. [49].
Firstly, the OE approach requires an unbiased prior SST field. In practice, the prior SST comes from either numerical reanalysis, as in our case, or from numerical weather prediction models such as used, for example, by Merchant et al. [30,50] or from operation applications. These prior SST fields are limited in both temporal and spatial resolution, typically being issued 4 times daily and with spatial resolution of tenths of a degree which represents a rather large surface compared to ~1 × 1 km typical pixel size of satellite infrared radiometers or a point buoy measurement. Thus, the prior fields are inherently likely to be biased and so, before OE can be attempted, the biases should be removed. Merchant [32] applies the OE to estimate the prior biases and possible measurement bias before applying the OE to retrieve SST from SEVIRI measurements. One could ask if, after that step, is the OE retrieval still necessary? Here, we estimated the biases directly to apply corrections to the a priori SST field. The results are summarized:
(a)
Through retrievals using situ buoy measurements used as the prior SST in the OE, we estimated that in our application the biases not related to those in the prior SST are on order of 0.02 K, which is small in comparison to the accuracy of buoy measurements of ~0.2 K [51,52,53,54] or the requirement of a CDR of SST of 0.1 K [5]. Newer surface drifters provide more accurate subsurface temperature measurements, with a target of 0.01 K [55,56], but thus far relatively few have been deployed.
(b)
It is relatively simple to estimate the consequences of neglecting diurnal signals that are missing from the foundation SST of the ERA5 data. A significant improvement in daytime OE SST retrievals is achieved by applying a simple diurnal correction, reducing the mean difference between OESST and buoy measurements from −0.25 K to −0.05 K.
(c)
Perhaps, somewhat surprisingly, this simple correction to ERA5 foundation SST provides better prior for the OE retrieval than the ERA5 skin temperature (SKT). The SKT–SSTb bias is larger than the ERA5 SST–SSTb bias even before the diurnal correction, and for both nighttime and daytime matchups.
(d)
Despite the improvement in the mean after application of the diurnal correction, there remain large values of OESST–SSTb. The large disagreements between OESST and SSTb occur mostly in the areas of the ocean where there are large horizontal gradients in the SST, i.e., coastal areas or at the surface expressions of oceanic thermal fronts. The limited spatial resolution of the numerical models providing prior SST fields often fails to adequately represent the spatial variability of SST in those areas. The spatial sampling bias is difficult to correct without an improvement in the resolution of the prior SST. We use the nearest ERA5 SST value for the a priori and more advanced interpolation could improve the estimates but might also introduce undesirable artifacts in areas of strong spatial variability.
(e)
In areas with strong horizontal thermal gradients, the NLSST often provides better estimates of SST than the OE approach.
(f)
In general, OESST improves on NLSST in terms of mean offset from in situ measurements, but not for the RMS, and improvement in the mean is driven mostly by the accuracy of the underlying prior SST with the OE retrieval not departing much from the prior.
(g)
The OESST departure from the prior is on average quite small so if the SST0–SSTb is large then the agreement between the retrieval and the in situ is not much improved.
(h)
On average, for very large differences between NLSST and in situ measurements, the OE provides a better estimate of SST but the OESST retrieval is not as accurate as the corrected prior.
In conclusion, we find that OE can improve estimates of SST from MODIS top-of-atmosphere channel radiance measurements if the biases are effectively removed from the prior SST field. This is not an easy task when it comes to biases resulting from spatial sampling, and often in areas of large SST variability the NLSST provides better retrievals than the OE. We proposed a diagnostic test to determine if the OE approach would be useful by evaluating the difference between the NLSST and the SST0: Cases with large NLSST–SST0 are not likely to lead to improved SST estimated by OE.
One of the main aspirations for the OE was to improve SST retrievals in situations where the atmospheric conditions depart from those implicit in the derivation of the NLSST retrieval coefficients, such as in areas of the oceans affected by large amounts of terrestrial aerosols that frequently move offshore typically from desert areas, and may span ocean basins having marked temporal and spatial variability [57,58,59,60]. The presence of aerosols in the marine atmosphere influences the IR radiation measured by satellite sensors. However, in the OE approach with the reduced state vector [SST, TCWV], this aerosol effect is not accounted for in the differences between the measured and modeled radiances which leads to a bias in the retrieved SST, as does using non-optimised coefficients lead to a bias in the NLSST algorithm in such conditions. Thus, to improve the OESST retrievals under aerosol layers, the OE would need to account explicitly for the presence of aerosols. That would of course also bring new uncertainties into the problem.
In view of these OE limitations, OESST might not always provide the best estimates of the SST. In particular, research related to smaller scale spatial variability of thermal currents, or SST in coastal zones may benefit from utilizing NLSST-generated satellite SST fields. Moreover, given the importance of oceanic fronts on atmospheric circulation especially in relation to extratropical cyclones and atmospheric rivers that develop mostly along the sharp fronts on SST [61,62,63,64] and the role of thermal fronts in marine ecosystems [65], the NLSST seems better suited to the detection of such fronts. This is of even greater importance as the positions of oceanic fronts change on several different time scales [66] and possibly shift in response to climate change [67,68].

Author Contributions

M.D.S.—Conceptualization, Formal analysis, Methodology, Visualization, Writing—original draft. P.J.M.—Conceptualization, Funding acquisition, Methodology, Project Administration, Writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by NASA, grant 80NSSC18K1515.

Data Availability Statement

The MODIS matchup data used here are a subset of those available from https://seabass.gsfc.nasa.gov/search/sst, accessed on 5 October 2021. The ERA5 data are available at https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview, accessed on 21 December 2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IPCC. Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2021; p. 3949. Available online: https://www.ipcc.ch/report/ar6/wg1/ (accessed on 21 March 2022).
  2. Minnett, P.J.; Alvera-Azcárate, A.; Chin, T.M.; Corlett, G.K.; Gentemann, C.L.; Karagali, I.; Li, X.; Marsouin, A.; Marullo, S.; Maturi, E.; et al. Half a century of satellite remote sensing of sea-surface temperature. Remote Sens. Environ. 2019, 233, 111366. [Google Scholar] [CrossRef]
  3. NRC. Climate Data Records from Environmental Satellites; National Academy of Sciences: Washington, DC, USA, 2004; p. 150. [Google Scholar] [CrossRef]
  4. NRC. Issues in the Integration of Research and Operational Satellite Systems for Climate Research: II. Implementation; National Academy of Sciences: Washington, DC, USA, 2000; Available online: http://www.nap.edu/catalog/9966.html (accessed on 21 March 2022).
  5. Ohring, G.; Wielicki, B.; Spencer, R.; Emery, B.; Datla, R. Satellite Instrument Calibration for Measuring Global Climate Change: Report of a Workshop. Bull. Am. Meteorol. Soc. 2005, 86, 1303–1313. [Google Scholar] [CrossRef]
  6. Prabhakara, C.; Dalu, G.; Kunde, V.G. Estimation of sea-surface temperature from remote sensing in 11 to 13.5mm window region. J. Geophys. Res. 1974, 80, 5039–5044. [Google Scholar] [CrossRef]
  7. McMillin, L. Estimation of sea-surface temperatures from two infrared window measurements with different absorption. J. Geophys. Res. 1975, 80, 5113–5117. [Google Scholar] [CrossRef]
  8. Walton, C.C.; Pichel, W.G.; Sapper, J.F.; May, D.A. The development and operational application of nonlinear algorithms for the measurement of sea surface temperatures with the NOAA polar-orbiting environmental satellites. J. Geophys. Res. 1998, 103, 27999–28012. [Google Scholar] [CrossRef]
  9. Barnes, W.L.; Xiong, X.; Salomonson, V.V. Status of Terra MODIS and Aqua MODIS. Adv. Space Res. 2003, 32, 2099–2106. [Google Scholar] [CrossRef]
  10. Kilpatrick, K.A.; Podesta, G.; Walsh, S.; Evans, R.; Minnett, P.J. Implementation of Version 6 AQUA and TERRA SST Processing; University of Miami: Miami, FL, USA, 2014; p. 14. Available online: https://oceancolor.gsfc.nasa.gov/reprocessing/r2014/sst/MODIS_V6_SST_changes.pdf (accessed on 21 March 2022).
  11. Kilpatrick, K.A.; Podestá, G.P.; Evans, R.H. Overview of the NOAA/NASA Pathfinder algorithm for Sea Surface Temperature and associated Matchup Database. J. Geophys. Res. 2001, 106, 9179–9198. [Google Scholar] [CrossRef]
  12. Závody, A.M.; Mutlow, C.T.; Llewellyn-Jones, D.T. A radiative transfer model for sea surface temperature retrieval for the along-track scanning radiometer. J. Geophys. Res. Ocean. 1995, 100, 937–952. [Google Scholar] [CrossRef]
  13. Blackmore, T.A.; O’Carroll, A.G.; Saunders, R.W.; Aumann, H. A Comparison of Sea Surface Temperatures from the Advanced Along-track Scanning Radiometer and the Atmospheric Infrared Sounder Instruments; UK Met Office: Exeter, UK, 2007; p. 21. [Google Scholar]
  14. Romaguera, M.; Sobrino, J.A.; Olesen, F.S. Estimation of sea surface temperature from SEVIRI data: Algorithm testing and comparison with AVHRR products. Int. J. Remote Sens. 2006, 27, 5081–5086. [Google Scholar] [CrossRef]
  15. Petrenko, B.; Ignatov, A.; Kihai, Y.; Pennybacker, M. Optimization of Sensitivity of GOES-16 ABI Sea Surface Temperature by Matching Satellite Observations with L4 Analysis. Remote Sens. 2019, 11, 206. [Google Scholar] [CrossRef] [Green Version]
  16. Park, K.-A.; Woo, H.-J.; Chung, S.-R.; Cheong, S.-H. Development of Sea Surface Temperature Retrieval Algorithms for Geostationary Satellite Data (Himawari-8/AHI). Asia Pac. J. Atmos. Sci. 2019, 56, 187–206. [Google Scholar] [CrossRef]
  17. Minnett, P.J.; Kilpatrick, K.A.; Podestá, G.P.; Evans, R.H.; Szczodrak, M.D.; Izaguirre, M.A.; Williams, E.J.; Walsh, S.; Reynolds, R.M.; Bailey, S.W.; et al. Skin Sea-Surface Temperature from VIIRS on Suomi-NPP—NASA Continuity Retrievals. Remote Sens. 2020, 12, 3369. [Google Scholar] [CrossRef]
  18. McClain, E.P.; Pichel, W.G.; Walton, C.C. Comparative performance of AVHRR-based multichannel sea surface temperatures. J. Geophys. Res. 1985, 90, 11587–11601. [Google Scholar] [CrossRef]
  19. Kilpatrick, K.A.; Podestá, G.; Walsh, S.; Williams, E.; Halliwell, V.; Szczodrak, M.; Brown, O.B.; Minnett, P.J.; Evans, R. A decade of sea surface temperature from MODIS. Remote Sens. Environ. 2015, 165, 27–41. [Google Scholar] [CrossRef]
  20. Llewellyn-Jones, D.T.; Minnett, P.J.; Saunders, R.W.; Zavody, A.M. Satellite multichannel infrared measurements of sea surface temperature of the N.E. Atlantic Ocean using AVHRR/2. Q. J. R. Meteorol. Soc. 1984, 110, 613–631. [Google Scholar] [CrossRef]
  21. Merchant, C.J.; Harris, A.R.; Murray, J.; Zavody, A.M. Toward the elimination of bias in satellite retrievals of skin sea surface temperature. 1: Theory. modelling and inter-algorithm comparison. J. Geophys. Res. 1999, 104, 23565–23578. [Google Scholar] [CrossRef]
  22. Minnett, P.J. The regional optimization of infrared measurements of sea-surface temperature from space. J. Geophys. Res. 1990, 95, 13497–13510. [Google Scholar] [CrossRef]
  23. Kumar, A.; Minnett, P.J.; Podesta, G.; Evans, R.H.; Kilpatrick, K. Analysis of Pathfinder SST algorithm for global and regional conditions. Proc. Indian Acad. Sci. Earth Planet. Sci. 2000, 109, 395–405. [Google Scholar] [CrossRef] [Green Version]
  24. Luo, B.; Minnett, P.J.; Nalli, N.R. Infrared satellite-derived sea surface skin temperature sensitivity to aerosol vertical distribution—Field data analysis and model simulations. Remote Sens. Environ. 2021, 252, 112151. [Google Scholar] [CrossRef]
  25. Luo, B.; Minnett, P.J.; Gentemann, C.; Szczodrak, G. Improving satellite retrieved night-time infrared sea surface temperatures in aerosol contaminated regions. Remote Sens. Environ. 2019, 223, 8–20. [Google Scholar] [CrossRef]
  26. Merchant, C.J.; Embury, O.; Le Borgne, P.; Bellec, B. Saharan dust in nighttime thermal imagery: Detection and reduction of related biases in retrieved sea surface temperature. Remote Sens. Environ. 2006, 104, 15–30. [Google Scholar] [CrossRef]
  27. Le Borgne, P.; Péré, S.; Roquet, H. Night time detection of Saharan dust using infrared window channels: Application to NPP/VIIRS. Remote Sens. Environ. 2013, 137, 264–273. [Google Scholar] [CrossRef]
  28. Bogdanoff, A.S.; Westphal, D.L.; Campbell, J.R.; Cummings, J.A.; Hyer, E.J.; Reid, J.S.; Clayson, C.A. Sensitivity of infrared sea surface temperature retrievals to the vertical distribution of airborne dust aerosol. Remote Sens. Environ. 2015, 159, 1–13. [Google Scholar] [CrossRef]
  29. Rodgers, C.D. Inverse Methods for Atmospheric Sounding: Theory and Practice; World Scientific: Singapore, 2000; p. 238. Available online: https://books.google.com/books?id=p3b3ngEACAAJ (accessed on 21 March 2022).
  30. Merchant, C.J.; LeBorgne, P.; Marsouin, A.; Roquet, H. Optimal estimation of sea surface temperature from split-window observations. Remote Sens. Environ. 2008, 112, 2469–2484. [Google Scholar] [CrossRef]
  31. Merchant, C.J.; Le Borgne, P.; Roquet, H.; Marsouin, A. Sea surface temperature from a geostationary satellite by optimal estimation. Remote Sens. Environ. 2009, 113, 445–457. [Google Scholar] [CrossRef]
  32. Merchant, C.J.; Saux-Picart, S.; Waller, J. Bias correction and covariance parameters for optimal estimation by exploiting matched in-situ references. Remote Sens. Environ. 2020, 237, 111590. [Google Scholar] [CrossRef]
  33. Koner, P.K.; Harris, A.; Maturi, E. A Physical Deterministic Inverse Method for Operational Satellite Remote Sensing: An Application for Sea SurfaceTemperature Retrievals. IEEE Trans. Geosci. Remote Sens. 2015, 53, 5872–5888. [Google Scholar] [CrossRef]
  34. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  35. Donlon, C.J.; Martin, M.; Stark, J.; Roberts-Jones, J.; Fiedler, E.; Wimmer, W. The Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) system. Remote Sens. Environ. 2012, 116, 140–158. [Google Scholar] [CrossRef]
  36. Good, S.; Fiedler, E.; Mao, C.; Martin, M.J.; Maycock, A.; Reid, R.; Roberts-Jones, J.; Searle, T.; Waters, J.; While, J.; et al. The Current Configuration of the Ostia System for Operational Production of Foundation Sea Surface Temperature and Ice Concentration Analyses. Remote Sens. 2020, 12, 720. [Google Scholar] [CrossRef] [Green Version]
  37. Donlon, C.J.; Minnett, P.J.; Gentemann, C.; Nightingale, T.J.; Barton, I.J.; Ward, B.; Murray, J. Toward improved validation of satellite sea surface skin temperature measurements for climate research. J. Clim. 2002, 15, 353–369. [Google Scholar] [CrossRef] [Green Version]
  38. Gentemann, C.L.; Minnett, P.J.; LeBorgne, P.; Merchant, C.J. Multi-satellite measurements of large diurnal warming events. Geophys. Res. Lett. 2008, 35, L22602. [Google Scholar] [CrossRef] [Green Version]
  39. Castro, S.L.; Wick, G.A.; Buck, J.J.H. Comparison of diurnal warming estimates from unpumped Argo data and SEVIRI satellite observations. Remote Sens. Environ. 2014, 140, 789–799. [Google Scholar] [CrossRef] [Green Version]
  40. Minnett, P.J. Radiometric measurements of the sea-surface skin temperature—The competing roles of the diurnal thermocline and the cool skin. Int. J. Remote Sens. 2003, 24, 5033–5047. [Google Scholar] [CrossRef]
  41. Fairall, C.; Bradley, E.; Godfrey, J.; Wick, G.; Edson, J.; Young, G. Cool-skin and warm-layer effects on sea surface temperature. J. Geophys. Res. 1996, 101, 1295–1308. [Google Scholar] [CrossRef]
  42. Hirahara, S.; Balmaseda, M.A.; de Boisseson, E.; Hersbach, H. Sea Surface Temperature and Sea Ice Concentration for ERA5; ECMWF: Reading, UK, 2016; p. 27. Available online: https://www.ecmwf.int/node/16555 (accessed on 21 March 2022).
  43. Saunders, R.; Hocking, J.; Turner, E.; Rayer, P.; Rundle, D.; Brunel, P.; Vidot, J.; Roquet, P.; Matricardi, M.; Geer, A.; et al. An update on the RTTOV fast radiative transfer model (currently at version 12). Geosci. Model Dev. 2018, 11, 2717–2737. [Google Scholar] [CrossRef] [Green Version]
  44. Hocking, J.; Rayer, P.; Rundle, D.; Saunders, R.; Matricardi, M.; Geer, A.; Brunel, P.; Vidot, J. RTTOV v12 Users Guide; EUMETSAT: Darmstadt, Germany, 2018; p. 142. Available online: https://www.nwpsaf.eu/site/download/documentation/rtm/docs_rttov12/NWPSAF-MO-DS-029_Prod_Spec_RTTOV12.pdf (accessed on 21 March 2022).
  45. Chen, B.; Liu, Z. Global water vapor variability and trend from the latest 36 year (1979 to 2014) data of ECMWF and NCEP reanalyses, radiosonde, GPS, and microwave satellite. J. Geophys. Res. Atmos. 2016, 121, 11,442-11,462. [Google Scholar] [CrossRef]
  46. ECMWF. IFS Documentation—Cy43r1. Part IV: Physical Processes; ECMWF: Reading, UK, 2016; p. 223. Available online: https://www.ecmwf.int/sites/default/files/elibrary/2016/17117-part-iv-physical-processes.pdf#section.H.10 (accessed on 21 March 2022).
  47. Luo, B.; Minnett, P. Evaluation of the ERA5 Sea Surface Skin Temperature with Remotely-Sensed Shipborne Marine-Atmospheric Emitted Radiance Interferometer Data. Remote Sens. 2020, 12, 1873. [Google Scholar] [CrossRef]
  48. Minnett, P.J.; Knuteson, R.O.; Best, F.A.; Osborne, B.J.; Hanafin, J.A.; Brown, O.B. The Marine-Atmospheric Emitted Radiance Interferometer (M-AERI), a high-accuracy, sea-going infrared spectroradiometer. J. Atmos. Ocean. Technol. 2001, 18, 994–1013. [Google Scholar] [CrossRef]
  49. Von Clarmann, T.; Degenstein, D.A.; Livesey, N.J.; Bender, S.; Braverman, A.; Butz, A.; Compernolle, S.; Damadeo, R.; Dueck, S.; Eriksson, P.; et al. Overview: Estimating and reporting uncertainties in remotely sensed atmospheric composition and temperature. Atmos. Meas. Techn. 2020, 13, 4393–4436. [Google Scholar] [CrossRef]
  50. Merchant, C.J.; Harris, A.R.; Roquet, H.; Le Borgne, P. Retrieval characteristics of non-linear sea surface temperature from the Advanced Very High Resolution Radiometer. Geophys. Res. Lett. 2009, 36, L17604. [Google Scholar] [CrossRef] [Green Version]
  51. O’Carroll, A.G.; Eyre, J.R.; Saunders, R.W. Three-Way Error Analysis between AATSR, AMSR-E, and In Situ Sea Surface Temperature Observations. J. Atmos. Ocean. Technol. 2008, 25, 1197–1207. [Google Scholar] [CrossRef] [Green Version]
  52. Gentemann, C.L. Three way validation of MODIS and AMSR-E sea surface temperatures. J. Geophys. Res. Ocean. 2014, 119, 2583–2598. [Google Scholar] [CrossRef]
  53. Kennedy, J.J.; Smith, R.O.; Rayner, N.A. Using AATSR data to assess the quality of in situ sea-surface temperature observations for climate studies. Remote Sens. Environ. 2012, 116, 79–92. [Google Scholar] [CrossRef]
  54. Kennedy, J.J. A review of uncertainty in in situ measurements and data sets of sea surface temperature. Rev. Geophys. 2014, 51, 2013RG000434. [Google Scholar] [CrossRef]
  55. Poli, P.; Lucas, M.; O’Carroll, A.; Le Menn, M.; David, A.; Corlett, G.K.; Blouch, P.; Meldrum, D.; Merchant, C.J.; Belbeoch, M. The Copernicus Surface Velocity Platform drifter with Barometer and Reference Sensor for Temperature (SVP-BRST): Genesis, design, and initial results. Ocean. Sci. 2019, 15, 199–214. [Google Scholar] [CrossRef] [Green Version]
  56. Le Menn, M.; Poli, P.; David, A.; Sagot, J.; Lucas, M.; O’Carroll, A.; Belbeoch, M.; Herklotz, K. Development of Surface Drifting Buoys for Fiducial Reference Measurements of Sea-Surface Temperature. Front. Mar. Sci. 2019, 6, 578. [Google Scholar] [CrossRef] [Green Version]
  57. Prospero, J.M.; Carlson, T.N. Vertical and areal distribution of Saharan dust over the western Equatorial North Atlantic Ocean. J. Geophys. Res. 1972, 77, 5255–5265. [Google Scholar] [CrossRef] [Green Version]
  58. Adams, A.M.; Prospero, J.M.; Zhang, C. CALIPSO-Derived Three-Dimensional Structure of Aerosol over the Atlantic Basin and Adjacent Continents. J. Clim. 2012, 25, 6862–6879. [Google Scholar] [CrossRef]
  59. Zuidema, P.; Alvarez, C.; Kramer, S.J.; Custals, L.; Izaguirre, M.; Sealy, P.; Prospero, J.M.; Blades, E. Is Summer African Dust Arriving Earlier to Barbados? The Updated Long-Term In Situ Dust Mass Concentration Time Series from Ragged Point, Barbados, and Miami, Florida. Bull. Am. Meteorol. Soc. 2019, 100, 1981–1986. [Google Scholar] [CrossRef]
  60. Tegen, I.; Schepanski, K. The global distribution of mineral dust. IOP Conf. Ser. Earth Environ. Sci. 2009, 7, 012001. [Google Scholar] [CrossRef]
  61. Minobe, S.; Kuwano-Yoshida, A.; Komori, N.; Xie, S.-P.; Small, R.J. Influence of the Gulf Stream on the troposphere. Nature 2008, 452, 206. [Google Scholar] [CrossRef]
  62. Minobe, S.; Miyashita, M.; Kuwano-Yoshida, A.; Tokinaga, H.; Xie, S.-P. Atmospheric Response to the Gulf Stream: Seasonal Variations. J. Clim. 2010, 23, 3699–3719. [Google Scholar] [CrossRef] [Green Version]
  63. Parfitt, R.; Czaja, A.; Minobe, S.; Kuwano-Yoshida, A. The atmospheric frontal response to SST perturbations in the Gulf Stream region. Geophys. Res. Lett. 2016, 43, 2299–2306. [Google Scholar] [CrossRef] [Green Version]
  64. Liu, X.; Ma, X.; Chang, P.; Jia, Y.; Fu, D.; Xu, G.; Wu, L.; Saravanan, R.; Patricola, C.M. Ocean fronts and eddies force atmospheric rivers and heavy precipitation in western North America. Nat. Commun. 2021, 12, 1268. [Google Scholar] [CrossRef]
  65. Belkin, I.M.; Cornillon, P.C.; Sherman, K. Fronts in Large Marine Ecosystems. Prog. Oceanogr. 2009, 81, 223–236. [Google Scholar] [CrossRef]
  66. Joyce, T.M.; Kwon, Y.-O.; Yu, L. On the Relationship between Synoptic Wintertime Atmospheric Variability and Path Shifts in the Gulf Stream and the Kuroshio Extension. J. Clim. 2009, 22, 3177–3192. [Google Scholar] [CrossRef]
  67. Hand, R.; Keenlyside, N.S.; Omrani, N.-E.; Bader, J.; Greatbatch, R.J. The role of local sea surface temperature pattern changes in shaping climate change in the North Atlantic sector. Clim. Dyn. 2019, 52, 417–438. [Google Scholar] [CrossRef] [Green Version]
  68. Chapman, C.C.; Lea, M.-A.; Meyer, A.; Sallée, J.-B.; Hindell, M. Defining Southern Ocean fronts and their influence on biological and physical processes in a changing climate. Nat. Clim. Change 2020, 10, 209–219. [Google Scholar] [CrossRef]
Figure 1. Elements A11 and A12 of the averaging matrix A as functions of (a,b) ξm, (c,d) eSST, and (e,f) eTCWV.
Figure 1. Elements A11 and A12 of the averaging matrix A as functions of (a,b) ξm, (c,d) eSST, and (e,f) eTCWV.
Remotesensing 14 02249 g001
Figure 2. Elements (a) A11 and (b) A12 of the averaging matrix A, and (c) logarithm of the condition number of the matrix that is to be inverted in Equation (5), and (d) the standard deviation of the retrieved SST distribution for a range of values of ξm and eSST. The dashed lines are isolines of eSST; the red line and the red dot are those from Figure 1 representing the values of ξm and eSST accepted as most likely in our case.
Figure 2. Elements (a) A11 and (b) A12 of the averaging matrix A, and (c) logarithm of the condition number of the matrix that is to be inverted in Equation (5), and (d) the standard deviation of the retrieved SST distribution for a range of values of ξm and eSST. The dashed lines are isolines of eSST; the red line and the red dot are those from Figure 1 representing the values of ξm and eSST accepted as most likely in our case.
Remotesensing 14 02249 g002
Figure 3. Elements (a) A22 and (b) A21 of the averaging matrix A for a range of values of ξm and eSST.
Figure 3. Elements (a) A22 and (b) A21 of the averaging matrix A for a range of values of ξm and eSST.
Remotesensing 14 02249 g003
Figure 4. Histograms of (a) nighttime, and (b) daytime OESST–SSTb for ERA5 SST0 prior with the nighttime or daytime correction, and (c) and (d) nighttime and daytime histograms of NLSST-SSTb.
Figure 4. Histograms of (a) nighttime, and (b) daytime OESST–SSTb for ERA5 SST0 prior with the nighttime or daytime correction, and (c) and (d) nighttime and daytime histograms of NLSST-SSTb.
Remotesensing 14 02249 g004
Figure 5. The fraction of the dataset with absolute value of the SST–SSTb above a given threshold for OESST, NLSST, and for ERA5 and ERA5 with the night/day correction: (a) Nighttime and (b) daytime data.
Figure 5. The fraction of the dataset with absolute value of the SST–SSTb above a given threshold for OESST, NLSST, and for ERA5 and ERA5 with the night/day correction: (a) Nighttime and (b) daytime data.
Remotesensing 14 02249 g005
Figure 6. Histograms of difference between ERA5 SST with the diurnal correction and SKT and buoy measurements: (a,b) SST–SSTb for nighttime and daytime matchups; and (c,d) SKT-SSTb for nighttime and daytime, respectively.
Figure 6. Histograms of difference between ERA5 SST with the diurnal correction and SKT and buoy measurements: (a,b) SST–SSTb for nighttime and daytime matchups; and (c,d) SKT-SSTb for nighttime and daytime, respectively.
Remotesensing 14 02249 g006
Figure 7. Locations of match-ups where ERA5 SST0 (with diurnal correction)–SSTb is (a) >1.0 K, (b) <−1.0 K, and (c) between −0.2 and 0.2 K, over contours of SST0 on 1 January 2015. Symbols indicate the type of in situ measurement.
Figure 7. Locations of match-ups where ERA5 SST0 (with diurnal correction)–SSTb is (a) >1.0 K, (b) <−1.0 K, and (c) between −0.2 and 0.2 K, over contours of SST0 on 1 January 2015. Symbols indicate the type of in situ measurement.
Remotesensing 14 02249 g007
Figure 8. Graphical representation of the statistics of SST0-SSTb, OESST-SSTb, and NLSST-SSTb from Table 1b.
Figure 8. Graphical representation of the statistics of SST0-SSTb, OESST-SSTb, and NLSST-SSTb from Table 1b.
Remotesensing 14 02249 g008
Figure 9. Graphical representation of the statistics of SST0-SSTb, OESST-SSTb, and NLSST-SSTb from Table 2b.
Figure 9. Graphical representation of the statistics of SST0-SSTb, OESST-SSTb, and NLSST-SSTb from Table 2b.
Remotesensing 14 02249 g009
Figure 10. As Figure 7, but for NLSST-SSTb.
Figure 10. As Figure 7, but for NLSST-SSTb.
Remotesensing 14 02249 g010
Figure 11. Correlation between OESST-SSTb (blue) and NLSST-SSTb (red) with NLSST-SST0.
Figure 11. Correlation between OESST-SSTb (blue) and NLSST-SSTb (red) with NLSST-SST0.
Remotesensing 14 02249 g011
Table 1. (a) Statistics of retrieved SST–SSTb differences in bins of SST0–SSTb for Q = 0 nighttime match-ups. Temperature differences are in K. OE StDev is the retrieved standard deviation of OESST distribution and ΔNLSST is the estimate of the error of NLSST obtained by propagation of error through the NLSST equation. (b) As (a), but for daytime.
Table 1. (a) Statistics of retrieved SST–SSTb differences in bins of SST0–SSTb for Q = 0 nighttime match-ups. Temperature differences are in K. OE StDev is the retrieved standard deviation of OESST distribution and ΔNLSST is the estimate of the error of NLSST obtained by propagation of error through the NLSST equation. (b) As (a), but for daytime.
a
SST0–SSTb
Bins
SST0–SSTbOESST–SSTbNLSST–SSTbN
Cases
MeanRMSStdMeanRMSStdOE
StDev
MeanRMSStdΔNL
SST
−5.0 > −1.0−1.5971.6660.472−1.5161.5830.4560.055−0.0670.7430.7400.284929
−1.0 > −0.5−0.6970.7110.141−0.6780.6970.1610.081−0.2250.7120.6750.3041057
−0.5 > −0.25−0.3500.3570.070−0.3480.3640.1080.107−0.285 0.6160.5460.3292147
−0.25 > −0.1−0.1700.1750.042−0.1700.1920.0880.117−0.2180.5340.4870.3402720
−0.1 > −0.05−0.0750.0760.014−0.0810.1130.0790.117−0.2010.5260.4860.3391126
−0.05 > 0−0.0250.0290.014−0.0340.0910.0850.122−0.1710.5220.4930.3421180
0 > 0.050.0250.0290.0150.0110.0800.0790.115−0.1440.4900.4900.3381128
0.05 > 0.10.0750.0760.0140.0570.1020.0840.118−0.1280.4690.4690.3401058
0.1 > 0.250.1690.1740.0430.1470.1700.0850.116−0.1110.4720.4720.3402453
0.25 > 0.50.3530.3600.0700.3120.3280.0990.105−0.0420.4940.4940.3322102
0.5 > 1.00.6800.6940.1370.6090.6260.1570.0830.0560.5730.5730.3111038
1.0 > 5.01.8241.9870.7891.6871.8470.7450.0630.3860.9160.9160.284569
b
SST0–SSTb
Bins
SST0–SSTbOESST–SSTbNLSST–SSTbN
Cases
MeanRMSStdMeanRMSStdOE
StDev
MeanRMSStdΔNL
SST
−5.0 > −1.0−1.6271.7100.525−1.4901.5690.4920.0670.1680.7150.6950.281932
−1.0 > −0.5−0.6890.7030.137−0.6490.6770.1930.102−0.0560.7800.7780.313996
−0.5 > −0.25−0.3600.3670.070−0.3450.3670.1300.122−0.1780.6270.6020.3351537
−0.25 > −0.1−0.1700.1760.043−0.1640.2000.1130.124−0.1470.5560.5360.3381757
−0.1 > −0.05−0.0750.0770.015−0.0710.1260.1040.123−0.1240.5140.4990.335666
−0.05 > 0−0.0250.0290.014−0.0190.0990.0970.121−0.0760.5810.5760.334771
0 > 0.050.0250.0290.0150.0170.1080.1070.122−0.0920.4990.4990.335785
0.05 > 0.10.0750.0770.0140.0700.1230.1010.119−0.0550.5080.5080.335750
0.1 > 0.250.1710.1760.0420.1550.1810.1110.110−0.0160.4920.4920.3291874
0.25 > 0.50.3550.3620.0700.3300.3510.1180.0980.0850.5450.5450.3171711
0.5 > 1.00.6730.6860.1330.6250.6460.1670.0840.2390.6640.6640.293840
1.0 > 5.01.7771.9430.7861.6511.8160.7540.0620.5491.0751.0750.264453
Table 2. (a) Statistics of retrieved SST–SSTb differences in bins of NLSST–SSTb for Q = 0 nighttime match-ups. OE StDev is the retrieved standard deviation of OESST distribution and ΔNLSST is the estimate of the error of NLSST obtained by propagation of error through the NLSST equation. (b) As (a), but for daytime.
Table 2. (a) Statistics of retrieved SST–SSTb differences in bins of NLSST–SSTb for Q = 0 nighttime match-ups. OE StDev is the retrieved standard deviation of OESST distribution and ΔNLSST is the estimate of the error of NLSST obtained by propagation of error through the NLSST equation. (b) As (a), but for daytime.
a
NLSST–SSTb
Bins
SST0–SSTbOESST–SSTbNLSST–SSTbN
Cases
MeanRMSStdMeanRMSStdOE
StDev
MeanRMSStdΔNL
SST
−5.0 > −1.0−0.2400.6790.636−0.3700.7090.6040.114−1.4151.4830.4430.3291009
−1.0 > −0.5−0.1240.4850.469−0.1670.4690.4380.124−0.6930.7060.1380.3382416
−0.5 > −0.25−0.0590.5070.504−0.0800.4820.4750.110−0.3620.3690.0710.3313028
−0.25 > −0.1−0.0450.5130.511−0.0530.4820.4790.104−0.1740.1790.0430.3292416
−0.1 > −0.05−0.0140.5140.514−0.0150.4780.4780.105−0.0740.0760.0150.331850
−0.05 > 0−0.0330.6000.599−0.0330.5570.5570.102−0.0240.0280.0150.329808
0 > 0.05−0.0150.5770.577−0.0180.5360.5360.0970.0250.0280.0280.328739
0.05 > 0.1−0.0140.6410.640−0.0100.5990.5990.0950.0730.0750.0740.324797
0.1 > 0.250.0250.5690.5680.0310.5310.5300.0980.1710.1770.1770.3261971
0.25 > 0.50.0690.7680.7650.0810.7140.7100.0980.3590.3660.3660.3252000
0.5 > 1.00.0460.8970.8960.0680.8380.8350.0940.6790.6930.6930.3211190
1.0 > 5.00.5081.5041.4160.5391.4371.3320.0811.4961.6331.6330.391292
b
NLSST–SSTb
Bins
SST0–SSTbOESST–SSTbNLSST–SSTbN
Cases
MeanRMSStdMeanRMSStdOE
StDev
MeanRMSStdΔNL
SST
−5.0 > −1.0−0.2720.6210.559−0.3860.6510.5230.120−1.3191.380.3630.321585
−1.0 > −0.5−0.1580.5380.515−0.2040.5240.4820.127−0.7070.7210.1420.3321563
−0.5 > −0.25−0.0860.5230.516−0.1070.4930.4810.113−0.3640.3710.0720.3261862
−0.25 > −0.1−0.0530.5230.521−0.0530.4930.4870.107−0.1730.1780.0420.3261508
−0.1 > −0.05−0.0290.5320.531−0.0360.4990.4980.104−0.0740.0760.0150.323585
−0.05 > 0−0.0490.5520.550−0.0380.5160.5140.106−0.0230.0280.0150.326636
0 > 0.05−0.0530.5350.532−0.0510.5070.5050.1020.0250.0280.0280.328514
0.05 > 0.1−0.0940.7090.702−0.0680.6500.6460.1050.0730.0740.0740.322590
0.1 > 0.25−0.0380.6070.605−0.0150.5620.5610.1000.1730.1780.1780.3181622
0.25 > 0.5−0.0040.7720.7720.0220.7150.7150.1010.3630.3670.3700.3191735
0.5 > 1.00.0230.9550.9550.0690.8860.8830.1000.6920.7050.7050.3141330
1.0 > 5.00.1461.2061.2000.2321.1341.1100.0891.5701.6881.6880.293542
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Szczodrak, M.D.; Minnett, P.J. Relative Merits of Optimal Estimation and Non-Linear Retrievals of Sea-Surface Temperature from MODIS. Remote Sens. 2022, 14, 2249. https://doi.org/10.3390/rs14092249

AMA Style

Szczodrak MD, Minnett PJ. Relative Merits of Optimal Estimation and Non-Linear Retrievals of Sea-Surface Temperature from MODIS. Remote Sensing. 2022; 14(9):2249. https://doi.org/10.3390/rs14092249

Chicago/Turabian Style

Szczodrak, Malgorzata D., and Peter J. Minnett. 2022. "Relative Merits of Optimal Estimation and Non-Linear Retrievals of Sea-Surface Temperature from MODIS" Remote Sensing 14, no. 9: 2249. https://doi.org/10.3390/rs14092249

APA Style

Szczodrak, M. D., & Minnett, P. J. (2022). Relative Merits of Optimal Estimation and Non-Linear Retrievals of Sea-Surface Temperature from MODIS. Remote Sensing, 14(9), 2249. https://doi.org/10.3390/rs14092249

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop