Earlier studies of spatial variability in sea surface temperature (SST) using ship-based radiometric data suggested that variability at scales smaller than 1 km is significant and affects the perceived uncertainty of satellite-derived SSTs. Here, we compare data from the Ball Experimental Sea Surface Temperature (BESST) thermal infrared radiometer flown over the Arctic Ocean against coincident Moderate Resolution Imaging Spectroradiometer (MODIS) measurements to assess the spatial variability of skin SSTs within 1-km pixels. By taking the standard deviation, σ, of the BESST measurements within individual MODIS pixels, we show that significant spatial variability of the skin temperature exists. The distribution of the surface variability measured by BESST shows a peak value of O(0.1) K, with 95% of the pixels showing σ < 0.45 K. Significantly, high-variability pixels are located at density fronts in the marginal ice zone, which are a primary source of submesoscale intermittency near the surface. SST wavenumber spectra indicate a spectral slope of −2, which is consistent with the presence of submesoscale processes at the ocean surface. Furthermore, the BESST wavenumber spectra not only match the energy distribution of MODIS SST spectra at the satellite-resolved wavelengths, they also span the spectral slope of −2 by ~3 decades, from wavelengths of 8 km to <0.08 km.

We evaluate the reliability and basic characteristics of observations of extreme DW events from current operational geostationary satellite sea surface temperature (SST) products through examination of three months of diurnal warming (DW) estimates derived by different methodologies from the Spinning Enhanced Visible and Infrared Imager on Meteosat-11, Advanced Himawari Imager on Himawari-8, and Advanced Baseline Imager on the Geostationary Operational Environmental Satellite (GOES)-16. This work primarily focuses on the following research questions: (1) Can these operational SST products accurately characterize extreme DW events? (2) What are the amplitudes and frequencies of these events? To answer these, we compute distributions of DW and DW exceedance and compare them amongst the different methods and geostationary sensors. Examination of the DW estimates demonstrates several challenges in accurately deriving distributions of the DW amplitude, particularly associated with estimating the “foundation” temperature and uncertainties in cloud screening. Overall, the results suggest that current geostationary sensors can reliably assess extreme DW, but the estimates are sensitive to the computational methods applied. We thus suggest careful processing/screening of the SST retrievals. We find a value of 3 K, corresponding to the 99th percentile, provides a potential practical threshold for extreme warming, but events of at least 6 K were reliably observed. Warming in excess of 6 K occurred somewhere an average of 47% of the time, and its probability at a given location was of O(10⁻⁶).

Available observations and a theoretical simulation are used to explore the consistency and relationship between sea surface temperature (SST) and sea ice concentration (SIC) within open-ocean-sea ice mixed satellite pixels as a function of grid resolution. The maximum limiting SST value for a specified SIC and spatial resolution is first examined within collocated satellite-derived products contained within existing Level 4 SST analyses distributed using the data specification from the Group for High Resolution Sea Surface Temperature. The shape of the interdependence is further validated with manually quality-controlled buoy SST and SIC collocations. A parametric equation for the limiting SST value is derived from simulations of a mixed ocean/ice pixel with specified ice fraction and a linear SST gradient extending away from the ice edge. The exponential curve matching the observed interdependence suggests a maximum 5 km pixel-averaged SST at SIC values approaching zero between 6 and 8 °C. This maximum value is significantly greater than the previously assumed limiting values of ~3 °C and the corresponding SST gradient is larger than those typically observed with satellite SST products, but agrees well with recent Saildrone SST observations near ice. The curve provides a conservative limit with which inconsistent SST/SIC pairings can be identified, not only near the ice edge but at intermediate ice concentrations. Application of the filter improves the agreement between the SST/SIC relationship in satellite products and available Saildrone observations as well as the internal consistency of the different satellite products.