Alexander Ignatov

The first NOAA full-mission reanalysis (RAN1) of the sea surface temperature (SST) from the two Moderate Resolution Imaging Spectroradiometers (MODIS) onboard Terra (24 February 2000–present) and Aqua (4 July 2002–present) was performed. The dataset was produced using the NOAA Advanced Clear-Sky Processor for Ocean (ACSPO) enterprise SST system from Collection 6.1 brightness temperatures (BTs) in three MODIS thermal emissive bands centered at 3.7, 11, and 12 µm with a spatial resolution of 1 km at nadir. In the initial stages of reprocessing, several instabilities in the MODIS SST time series were observed. In particular, Terra SSTs and corresponding BTs showed three ‘steps’: two on 30 October 2000 and 2 July 2001 (due to changes in the MODIS operating mode) and one on 25 April 2020 (due to a change in its nominal blackbody temperature, BBT, from 290 to 285 K). Additionally, spikes up to several tenths of a kelvin were observed during the quarterly warm-up/cool-down (WUCD) exercises, when the Terra MODIS BBT was varied. Systematic gradual drifts of ~0.025 K/decade were also seen in both Aqua and Terra SSTs over their full missions due to drifting BTs. These calibration instabilities were mitigated by debiasing MODIS BTs using the time series of observed minus modeled (‘O-M’) BTs. The RAN1 dataset was evaluated via comparisons with various in situ SSTs. The data meet the NOAA specifications for accuracy (±0.2 K) and precision (0.6 K), often by a wide margin, in a clear-sky ocean domain of 19–21%. The long-term SST drift is typically less than 0.01 K/decade for all MODIS SSTs, except for the daytime ‘subskin’ SST, for which the drift is ~0.02 K/decade. The MODIS RAN1 dataset is archived at NOAA CoastWatch and updated monthly in a delayed mode with a latency of two months. Additional archival with NASA JPL PO.DAAC is being discussed.

The NOAA in situ Sea Surface Temperature (SST) Quality Monitor (iQuam) online system collects in situ SSTs from various sources, performs quality control (QC), and provides QC’ed data to users. Like many other in situ QCs, the iQuam QC employs comparisons with Level 4 SST analysis. However, the current daily L4 analyses do not capture the diurnal cycle, nor do they resolve the fine structure of SST in dynamic areas. As a result, high-quality in situ SSTs significantly deviating from the L4 SST may be rejected. This paper discusses the new Diurnal Reference Check (DRC), which addresses overscreening for buoys whose sampling frequency is sufficient for resolving the diurnal cycle. The DRC separates records from individual buoys into 24-h segments and characterizes each segment with the median nighttime (MNT) SST and the amplitude of the diurnal signal (ADS). The segment is rejected if the ADS is unrealistically large or if the difference between the MNT and L4 SST exceeds a geographically dependent threshold. The outliers are further screened out by comparison of individual in situ SSTs with the MNT. All thresholds are determined from the analysis of matchups with reprocessed NOAA SSTs from multiple low-orbiting satellites. The satellite matchups are also used to validate the QC results. The DRC minimizes the overscreening, increases the number of high-quality in situ data by ~5%, and reduces the QC reliance on the L4 analysis. In addition, a new retrospective satellite-based quality check is introduced to identify matchups, which are most useful for training SST algorithms and validation of reprocessed satellite data.

The 3rd full-mission reanalysis (RAN3) of global sea surface temperature (SST) with a 750 m resolution at nadir is available from VIIRS instruments flown onboard two JPSS satellites: NPP (February 2012–present) and N20 (January 2018–present). Two SSTs, ‘subskin’ (sensitive to skin SST) and ‘depth’ (proxy for in situ SST at depth of 20 cm), were produced from brightness temperatures (BTs) in the VIIRS bands centered at 8.6, 11 and 12 µm during the daytime and an additional 3.7 µm band at night, using the NOAA Advanced Clear Sky Processor for Ocean (ACSPO) enterprise SST system. The RAN3 dataset is fully archived at NASA JPL PO.DAAC and NOAA CoastWatch, and routinely supplemented in near real time (NRT) with a latency of a few hours. Delayed mode (DM) processing with a 2 months latency follows NRT, resulting in a more uniform science quality SST record. This paper documents and evaluates the performance of the VIIRS RAN3 dataset. Comparisons with in situ SSTs from drifters and tropical moorings (D+TM) as well as Argo floats (AFs) (both available from the NOAA iQuam system) show good agreement, generally within the NOAA specifications for accuracy (±0.2 K) and precision (0.6 K), in a clear-sky domain covering 18–20% of the global ocean. The nighttime SSTs compare with in situ data more closely, as expected due to the reduced diurnal thermocline. The daytime SSTs are also generally within NOAA specs but show some differences between the (D+TM) and AF validations as well as residual drift on the order of −0.1 K/decade. BT comparisons between two VIIRSs and MODIS-Aqua show good consistency in the 3.7 and 12 µm bands. The 11 µm band, while consistent between NPP and N20, shows residual drift with respect to MODIS-Aqua. Similar analyses of the 8.6 µm band are inconclusive, as the performance of the MODIS band 29 centered at 8.6 µm is degraded and unstable in time and cannot be used for comparisons.

The 40+ years-long sea surface temperature (SST) dataset from 4 km Global Area Coverage (GAC) data of the Advanced Very High-Resolution Radiometers (AVHRR/2s and/3s) flown onboard ten NOAA satellites (N07/09/11/12/14/15/16/17/18/19) has been created under the NOAA AVHRR GAC SST Reanalysis 2 (RAN2) Project. The data were reprocessed with the NOAA Advanced Clear Sky Processor for Ocean (ACSPO) enterprise SST system. Two SST products are reported in the full ~3000 km AVHRR swath: ‘subskin’ (highly sensitive to true skin SST, but debiased with respect to in situ SST) and ‘depth’ (a closer proxy for in situ data, but with reduced sensitivity). The reprocessing methodology aims at close consistency of satellite SSTs with in situ SSTs, in an optimal retrieval domain. Long-term orbital and calibration trends were compensated by daily recalculation of regression coefficients using matchups with drifters and tropical moored buoys (supplemented by ships for N07/09), collected within limited time windows centered at the processed day. The nighttime Sun impingements on the sensor black body were mitigated by correcting the L1b calibration coefficients. The Earth view pixels contaminated with a stray light were excluded. Massive cold SST outliers caused by volcanic aerosols following three major eruptions were filtered out by a modified, more conservative ACSPO clear-sky mask. The RAN2 SSTs are available in three formats: swath L2P (144 10-min granules per 24 h interval) and two 0.02° gridded (uncollated L3U, also 144 granules/24 h; and collated L3C, two global maps per 24 h, one for day and one for the night). This paper evaluates the RAN2 SST dataset, with a focus on the L3C product and compares it with two other available AVHRR GAC L3C SST datasets, NOAA Pathfinder v5.3 and ESA Climate Change Initiative v2.1. Among the three datasets, the RAN2 covers the global ocean more completely and shows reduced regional and temporal biases, improved stability and consistency between different satellites, and in situ SSTs.

Marine boundary layer (MBL) clouds, ubiquitous to the world’s oceans, help govern the radiative balance of Earth’s climate system. Satellite remote sensing provides the most practical means to monitor cloud worldwide. Whereas visible‐based detection of MBL clouds from environmental satellites is relatively straightforward during the daytime, the night presents challenges. In certain conditions, the conventional infrared (IR) methods used for nocturnal cloud detection, such as the commonly used 11–3.9 μm brightness temperature difference, offer poor thermal and spectral contrast between clouds and the clear‐sky background—resulting in missed clouds. A less explored question is to what extent these IR techniques overstate the cloud field. The Day/Night Band (DNB), a novel low‐light sensor carried on the Joint Polar Satellite System constellation, is helping to address this question. By way of its daytime‐analog moonlight reflectance cloud detection, the DNB reveals situations where IR‐based techniques yield false‐alarm low clouds. These problems occur in conditions of cool surface and a warm/moist lower atmosphere, prevalent in areas of coastal upwelling, tidal mixing, river estuaries, and along oceanic fronts and cyclonic eddies. We show how differential sensitivity to MBL moisture can trigger algorithmic thresholds for cloud detection. Presented here are examples illustrating the problem and an evaluation of our hypothesis against idealized radiative transfer simulations. We consider the implications of such artifacts, including climate data records of sea surface temperature which rely upon IR‐based nocturnal cloud masks. The results provide a framework for designing algorithmic improvements to nighttime MBL cloud detection.

The first full-mission global AVHRR FRAC sea surface temperature (SST) dataset with a nominal 1.1 km resolution at nadir was produced from three Metop First Generation (FG) satellites: Metop-A (2006-on), -B (2012-on) and -C (2018-on), using the NOAA Advanced Clear Sky Processor for Ocean (ACSPO) SST enterprise system. Historical reprocessing (‘Reanalysis-1’, RAN1) starts at the beginning of each mission and continues into near-real time (NRT). ACSPO generates two SST products, one with global regression (GR; highly sensitive to skin SST), and another one with piecewise regression (PWR; proxy for depth SST) algorithms. Small residual effects of orbital and sensor instabilities on SST retrievals are mitigated by retraining the regression coefficients daily, using matchups with drifting and tropical moored buoys within moving time windows. In RAN, the training windows are centered at the processed day. In NRT, the same size windows are employed but delayed in time, ending four to ten days prior to the processed day. Delayed-mode RAN reprocessing follows the NRT with a two-month lag, resulting in a higher quality and a more consistent SST record. In addition to its completeness, the newly created Metop-FG RAN1 SST dataset shows very close agreement with in situ data (including the fully independent Argo floats), well within the NOAA specifications for accuracy (global mean bias; ±0.2 K) and precision (global standard deviation; 0.6 K) in a ~20% clear-sky domain (percent of clear-sky SST pixels to the total of ice-free ocean). All performance statistics are stable in time, and consistent across the three platforms. The Metop-FG RAN1 data set is archived at the NASA JPL PO.DAAC and NOAA NCEI. This paper documents the newly created dataset and evaluates its performance.

In situ sea surface temperatures (SST) are the key component of the calibration and validation (Cal/Val) of satellite SST retrievals and data assimilation (DA). The NOAA in situ SST Quality Monitor (iQuam) aims to collect, from various sources, all available in situ SST data, and integrate them into a maximally complete, uniform, and accurate dataset to support these applications. For each in situ data type, iQuam strives to ingest data from several independent sources, to ensure most complete coverage, at the cost of some redundancy in data feeds. The relative completeness of various inputs and their consistency and mutual complementarity are often unknown and are the focus of this study. For four platform types customarily employed in satellite Cal/Val and DA (drifting buoys, tropical moorings, ships, and Argo floats), five widely known data sets are analyzed: (1) International Comprehensive Ocean-Atmosphere Data Set (ICOADS), (2) Fleet Numerical Meteorology and Oceanography Center (FNMOC), (3) Atlantic Oceanographic and Meteorological Laboratory (AOML), (4) Copernicus Marine Environment Monitoring Service (CMEMS), and (5) Argo Global Data Assembly Centers (GDACs). Each data set reports SSTs from one or more platform types. It is found that drifting buoys are more fully represented in FNMOC and CMEMS. Ships are reported in FNMOC and ICOADS, which are best used in conjunction with each other, but not in CMEMS. Tropical moorings are well represented in ICOADS, FNMOC, and CMEMS. Some CMEMS mooring reports are sampled every 10 min (compared to the standard 1 h sampling in all other datasets). The CMEMS Argo profiling data set is, as expected, nearly identical with those from the two Argo GDACs.

Gliders are low‐power autonomous underwater vehicles used to obtain oceanic measurements in vertical sections. Assimilation of glider temperature and salinity into coastal ocean circulation models holds the potential to improve the ocean subsurface structure estimate. In this study, the impact of assimilation of glider observations is studied using a four‐dimensional variational (4DVAR) data assimilation and forecast system set offshore of Oregon and Washington on the U.S. West Coast. Four test cases are compared: (1) no assimilation, (2) assimilation of glider temperature and salinity data alone, (3) assimilation of the glider data in combination with the surface observations including satellite sea surface temperature, sea surface height, and high‐frequency radar surface velocities, and (4) assimilation of the surface data alone. It is found that the assimilation of glider observations alone creates unphysical eddies in the vicinity of the glider transect. As a consequence, the forecast errors in the surface velocity and temperature increase compared to the case without data assimilation. Assimilation of surface and subsurface observations in combination prevents these features from forming and reduces the errors in the forecasts for the subsurface fields compared to the other three experiments. These improvements persisted in 21‐day forecasts run after the last data assimilation cycle.

Monitoring of the diurnal warming cycle in sea surface temperature (SST) is one of the key tasks of the new generation geostationary sensors, the Geostationary Operational Environmental Satellite (GOES)-16/17 Advanced Baseline Imager (ABI), and the Himawari-8/9 Advanced Himawari Imager (AHI). However, such monitoring requires modifications of the conventional SST retrieval algorithms. In order to closely reproduce temporal and spatial variations in SST, the sensitivity of retrieved SST to SSTskin should be as close to 1 as possible. Regression algorithms trained by matching satellite observations with in situ SST from drifting and moored buoys do not meet this requirement. Since the geostationary sensors observe tropical regions over larger domains and under more favorable conditions than mid-to-high latitudes, the matchups are predominantly concentrated within a narrow range of in situ SSTs >2 85 K. As a result, the algorithms trained against in situ SST provide the sensitivity to SSTskin as low as ~0.7 on average. An alternative training method, employed in the National Oceanic and Atmospheric Administration (NOAA) Advanced Clear-Sky Processor for Oceans, matches nighttime satellite clear-sky observations with the analysis L4 SST, interpolated to the sensor’s pixels. The method takes advantage of the total number of clear-sky pixels being large even at high latitudes. The operational use of this training method for ABI and AHI has increased the mean sensitivity of the global regression SST to ~0.9 without increasing regional biases. As a further development towards improved SSTskin retrieval, the piecewise regression SST algorithm was developed, which provides optimal sensitivity in every SST pixel. The paper describes the global and the piecewise regression algorithms trained against analysis SST and illustrates their performance with SST retrievals from the GOES-16 ABI.

The Suomi National Polar-orbiting Partnership Program Visible Infrared Imaging Radiometer Suite (VIIRS) thermal emissive bands (TEB) have been performing well during nominal operations since launch. However, small but persistent calibration anomalies are observed in all TEBs during the quarterly blackbody (BB) warm-up/cool-down (WUCD) events. As a result, the time series of daytime sea surface temperature (SST) (derived from bands M15-M16) show warm spikes on the order of 0.25 K. This paper suggests that VIIRS TEB WUCD biases are band dependent, with daily-averaged biases about -0.04 and 0.05 K for I4 and I5, and -0.05, -0.05, 0.11, 0.09, and 0.05 K for M12-M16, respectively. Two correction methods--Ltrace and WUCD-C--have been implemented and evaluated using colocated observations from the Cross-track Infrared Sounder (CrIS), radiative transfer simulations, and SST retrievals. Also an error in the National Oceanic and Atmospheric Administration operational processing was identified and fixed. Both correction methods effectively minimize WUCD-induced SST anomalies. The Ltrace method works well for I5, M12, and M14-M16, with residual biases about 0.01 K. The WUCD-C method, on the other hand, performs well to correct WUCD biases in all TEBs, with residual biases also about 0.01 K. However, it introduces warm biases relative to CrIS at cold scene temperatures, which requires further study. Applying nonequal BB thermistor weights improves calibration at BB temperature set points, but its impact on daily-averaged WUCD biases is small. The proposed methodologies may also be applied to the VIIRS onboard the follow-on Joint Polar Satellite System satellites.