Validation of the Suomi NPP VIIRS Ice Surface Temperature Environmental Data Record

Continuous monitoring of the surface temperature is critical to understanding and forecasting Arctic climate change; as surface temperature integrates changes in the surface energy budget. The sea-ice surface temperature (IST) has been measured with optical and thermal infrared sensors for many years. With the IST Environmental Data Record (EDR) available from the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (NPP) and future Joint Polar Satellite System (JPSS) satellites; we can continue to monitor and investigate Arctic climate change. This work examines the quality of the VIIRS IST EDR. Validation is performed through comparisons with multiple datasets; including NASA IceBridge measurements; air temperature from Arctic drifting ice buoys; Moderate Resolution Imaging Spectroradiometer (MODIS) IST; MODIS IST simultaneous nadir overpass (SNO); and surface air temperature from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis. Results show biases of  ́0.34;  ́0.12; 0.16; ́3.20; and  ́3.41 K compared to an aircraft-mounted downward-looking pyrometer; MODIS; MODIS SNO; drifting buoy; and NCEP/NCAR reanalysis; respectively; root-mean-square errors of 0.98; 1.02; 0.95; 4.89; and 6.94 K; and root-mean-square errors with the bias removed of 0.92; 1.01; 0.94; 3.70; and 6.04 K. Based on the IceBridge and MODIS results; the VIIRS IST uncertainty (RMSE) meets or exceeds the JPSS system requirement of 1.0 K. The product can therefore be considered useful for meteorological and climatological applications.


Introduction
The Arctic has been warming at a greater rate than anywhere else on Earth, a trend that is projected to continue over the next century [1].Continuous measurements of the surface temperature are therefore important to the understanding, monitoring, and forecasting Arctic climate change.The surface temperature of sea ice, hereafter referred to as the "ice surface temperature (IST)", is the controlling factor for sea ice growth and melt.Furthermore, the strength of low-level atmospheric temperature inversions depends on the IST.IST and precipitable water are interdependent [2], as both the lower tropospheric temperature structure and water vapor content determine the downwelling longwave radiation.Consequently, changes in sea ice characteristics, temperature inversion strength, precipitable water, and associated cloud cover all feed back to the surface temperature.Therefore, accurate and spatially robust surface temperature measurements are required to better understand feedbacks within the Arctic climate system and to better quantify the degree to which the Arctic is warming at a greater rate than the global average, known as "Arctic amplification" [3,4].
Satellite imagers with longwave infrared bands have been used to measure the ice surface temperature under clear conditions for many years.Key and Haefliger [5] presented an algorithm for estimating IST using split-window thermal channels at 11 and 12 µm of the Advanced Very High Resolution Radiometer (AVHRR) sensors.Key et al. [6] presented some improvements to the Key and Haefliger [5] procedure.The method was applied to a multi-decadal time series of AVHRR to study recent climate change in the Polar Regions [7][8][9][10][11].Hall et al. [12] adopted the Key et al. [6] method for use with the Moderate-Resolution Imaging Spectroradiometer (MODIS) onboard the Terra and Aqua satellites.Their validation yielded a root mean square error (RMSE) of 1.2 K when compared to near-surface air temperature.Methods have also been developed for the estimation of land surface temperature from remotely sensed satellite-based thermal infrared data [13].The accuracy of those retrievals has been assessed in a variety of validation studies [13][14][15][16].
IST is now produced with the Visible Infrared Imager Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (NPP) satellite.It is one of the official VIIRS products, or Environmental Data Records (EDR).VIIRS has some advantages over its heritage sensors [17] in that the split window bands at 11 and 12 µm in "moderate" resolution have higher spatial resolution than both AVHRR and MODIS, with approximately 1.6 km for VIIRS at the edge of the scan (750 m at nadir) compared to 4.8 km for MODIS (1 km at nadir).VIIRS also has a wider swath than MODIS, approximately 3000 km versus 2330 km. Figure 1 gives an example of a daily composite of the VIIRS IST EDR on 1 March 2015.It shows the lowest IST north of the Canadian Archipelago and Greenland, higher temperatures north of Alaska, and the warmest areas around the Kara Sea and Bering Strait.This study presents the validation of the VIIRS IST EDR with multiple datasets, including aircraft measurements, 2 m air temperature from Arctic drifting buoys, MODIS IST, and surface air temperatures from an atmospheric reanalysis product.This work significantly extends the preliminary validation study of Key et al. [18] by employing additional IST product inter-comparisons, thereby providing a robust evaluation of the operational VIIRS IST product.

Data and Methods
The VIIRS IST EDR provides clear-sky surface skin temperature retrievals at VIIRS moderate resolution (750 m) for both day and night over snow and ice covered water surfaces.The VIIRS IST EDR algorithm is based on the split-window technique of Yu et al. [19].This statistical regression method is implemented as: IST " a 0 `a1 T 11 `a2 pT 11 ´T12 q `a3 psec pzq ´1q where T 11 and T 12 are the VIIRS brightness temperature at 10.8 µm and 12.0 µm in moderate resolution bands M15 and M16 respectively, z is the satellite zenith angle, and a 0 , a 1 , a 2 , and a 3 are regression coefficients.Different sets of coefficients are used for daytime and nighttime retrievals.Day is defined by solar zenith angles (SZA) less than 85 degrees.Retrievals are done only under clear sky conditions, which are identified with the VIIRS cloud mask [20].The brightness temperature difference term in Equation ( 1) is for the correction of atmospheric water vapor absorption; the term that includes the satellite zenith angle accounts for the atmosphere path length.In this study, VIIRS IST EDR data produced by the Interface Data Processing Segment (IDPS) over the period August 2012 to July 2015 are used.
The algorithm described above is similar to the ice surface temperature algorithm first developed by Key and Haefliger [5] and later revised by Key et al. [6].This approach is also similar to those used for sea surface temperature (SST) estimation.Other related algorithms are described by Comiso [21] and Lindsay and Rothrock [22].
The IDPS operational IST product [23] is based on prelaunch regression coefficients [24] derived from Global Synthetic Data [25,26] that includes the effects of VIIRS sensor relative spectral response and sensor noise, based on VIIRS F1 prelaunch, sensor characterization test measurements.At-sensor, clear sky, top-of-atmosphere (TOA) synthetic radiances for the VIIRS 10.8 µm and 12.0 µm IR bands were computed using the MODerate resolution atmospheric TRANsmission (MODTRAN) [27] radiative transfer model with inputs derived from physical properties based on NCEP Global Forecast System (GFS) snow/ice mask, surface skin-temperature and atmospheric data.VIIRS sensor effects were incorporated into the simulated radiances based on a VIIRS sensor model [28].

Validation Datasets
The datasets used to validate the VIIRS IST are summarized in Table 1.IST data sets from MODIS, specifically the Terra MODIS (MOD29) and Aqua MODIS (MYD29) Sea Ice Extent 5-Min L2 Swath 1 km datasets [29,30] from August 2012 to July 2015 are used to compare with the VIIIRS IST EDR.These data sets are in swath format at 1 km resolution at nadir for both daytime and nighttime under clear sky conditions, with swath coverage of 2330 km (cross track) by 2030 km (along track).To collocate the VIIRS IST EDR and MODIS IST, every overpass over the Arctic and Antarctic of the VIIRS IST EDR and MODIS IST are re-gridded to a 1 km Equal-Area Earth Grid (EASE-Grid), with 9025 by 9025 grid cells extending from 48.4 to 90 ˝N over the Arctic and 8025 by 8025 grid cells extending from 53.2 to 90 ˝S over the Antarctic.Only grids cells with VIIRS-MODIS observation time differences of less than 5 min, and indicated as cloud free from both VIIRS and MODIS cloud masks, are included in the matchup comparison.An additional comparison of VIIRS and MYD29 IST values based on S-NPP and Aqua Simultaneous Nadir Overpass (SNO) orbit swaths within a 5 min time window covering the Arctic is also performed.The matchup IST values for the SNO comparisons are binned and averaged to a 1/20th degree equal angle, Climate Model Grid (CMG).The SNO IST comparison matchups include only pixels indicated as confidently clear by the VIIRS cloud mask with satellite view zenith angles less than 45 degrees and ice concentrations greater than or equal to 95%, as estimated by the VIIRS Ice Concentration Intermediate Product (IP) [31].For this comparison, 20 SNO orbit swaths from March 2013 to April 2014 are used, thus capturing a relatively broad range of temperatures.
IST measured with an airborne Heitronics KT-19.85Series II infrared radiation pyrometer (KT-19) on a National Aeronautics and Space Administration (NASA) P3 during the NASA Operation IceBridge campaigns [32] is used to validate the VIIRS IST EDR.One of the NASA IceBridge mission goals is to better understand the processes that connect Polar Regions with the global climate system [33].All KT-19 observations in the Antarctic in 2012 and 2013, and in the Arctic in 2014, are used to validate collocated VIIRS IST EDR.All available KT-19 temperature measurements that are within 375 m of the center of the VIIRS pixel are averaged.Only KT-19 temperature samples with standard deviations less than 1 K are used in matchup comparisons in order to eliminate cases with small-scale IST outliers that VIIRS would not be able to resolve.In addition, the closest KT-19 point to the VIIRS pixel must be within 15 min and 100 meters to be considered.Finally, a rigorous quality control is done to make sure that the VIIRS IST value is not contaminated by cloud.
If a daylight (SZA less than 85 ˝) scene is used, a false-color image is created with VIIRS bands M5, M7, and M10 (0.67, 0.85, and 1.6 µm).These images are visually inspected to confirm that the daytime scene is clear sky.For the few available scenes that are night (SZA 90 ˝or more), a visual inspection of the M15 and M16 (10.8, 12.0 µm) images is done.It should be noted that no cases with SZA between 80 ˝and 90 ˝are used in the statistics because of the large increase in random error observed in ice and cloud products for such low-sun conditions.
The NCEP/NCAR Reanalysis dataset [34] provides surface air temperature at the 0.995 sigma level from 1948 to the present.The daily surface air temperature data is available at 00Z, 06Z, 12Z, and 18Z with a spatial resolution of 2.5 ˆ2.5 degrees latitude/longitude.The following approach is used to match the NCEP/NCAR reanalysis surface temperature with VIIRST IST EDR.For each VIIRS IST EDR overpass, the IST values in the 1 km EASE-Grid projection are remapped to a 2.5 ˆ2.5 degree latitude/longitude grid.If the ratio of available remapped samples to the maximum possible sample number in a grid cell is larger than a threshold of 0.9, the mean of all the remapped VIIRS IST EDRs is assigned to that grid cell.Using the median of all the remapped VIIRS IST EDRs produces the same results, as shown in the following section.
Surface air temperatures (SAT) measured by drifting ice buoys from the International Arctic Buoy Programme (IABP) [35] are also used to assess the accuracy of the VIIRS IST EDR.Unfortunately, issues such as solar heating, frost and snow can lead to inaccurate SAT buoy measurements.Sensors with the most accurate measurements are the 2 m, shielded thermistors (personal communication, I. Rigor, 2015).In this study, SAT measurements from four drifting buoys from August 2012 to June 2014 are used, with buoy IDs 300025010128510, 300025010125530, 300025010123530, and 300234011045700, which are all ice mass balance buoys.More details of these buoys are available at http://iabp.apl.washington.edu/maps_daily_table.html.To collocate the buoy SAT observations and the VIIRS IST EDR, both are remapped to 1 km EASE-Grid, and only cases of collocated buoy SAT measurements and the VIIRS IST EDR within 2 h of each other are included.

Definition of Statistics
In the following sections three quantities are used to quantify the differences between two datasets.The bias, , is the average difference between observations in two datasets: where e i is the difference between individual observations in the two samples-the "errors"-and n is the sample size.The root mean square error (RMSE) is the square root of the mean squared difference between two datasets: The RMSE with the bias removed is the square root of the average squared deviation of the errors from the mean error, which is the standard deviation of the errors: For the JPSS program, the absolute value of the bias in Equation ( 2) is termed the measurement accuracy.RMSE in Equation ( 3) is the measurement uncertainty.The RMSE with the bias removed, or the standard deviation, in Equation ( 4) is the measurement precision.If the bias is zero, RMSE and RMSE nobias are equal.

Comparison with MODIS IST
The MODIS IST product has been validated with surface air temperatures from weather stations and from buoy air temperatures, with a reported RMSE in the range of 1.2-1.3K [12].The MODIS IST product has been used in climate-related studies [29,30].Comparisons of VIIRS IST EDR with MODIS IST can therefore demonstrate both consistency and accuracy.Histograms of the differences between VIIRS IST and MODIS IST match-ups are shown in Figure 2 for all cases, and in Figure 3 for cases with MODIS IST in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K.For over 2 billion match-up pairs, the VIIRS IST EDR shows a bias of ´0.12 K with MODIS IST, and the RMSE nobias of 1.01 K.The majority of the IST differences have absolute values less than 1 K.This shows that in general the VIIRS IST EDR can be considered to be as accurate as the MODIS IST product.The good agreement can be attributed to similar split-window retrieval approach of both products, and good quality clear/cloudy identification using both the VIIRS and MODIS cloud masks.For the temperature bin with the MODIS IST between 270 and 275 K the bias is positive (0.22 K).The bias becomes negative (−0.10 K) for temperatures below 270 K, and becomes more negative (−0.47) for temperatures less than 230 K.The RMSEnobias values are relatively stable, with smaller For the temperature bin with the MODIS IST between 270 and 275 K the bias is positive (0.22 K).The bias becomes negative (´0.10 K) for temperatures below 270 K, and becomes more negative (´0.47) for temperatures less than 230 K.The RMSE nobias values are relatively stable, with smaller values for lower temperatures.RMSE is close to RMSE nobias because of the small bias, with RMSE of 1.02 K for all cases, and 1.00 K, 0.96 K, 1.05 K, 1.31 K, 1.31 K, and 1.22 K for MODIS IST in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K.The bias characteristics can be attributed to the differences between how the retrieval coefficients are derived for VIIRS and MODIS ISTs.For the VIIRS IST EDR, different sets of coefficients are used for daytime and nighttime.For MODIS, different retrieval coefficients are derived for three temperature ranges of the 11 µm brightness temperature: less than 240 K, 240-260 K, and greater than 260 K. Thus, the MODIS IST retrievals may benefit from this approach, and have a similar performance for all three temperature intervals.The VIIRS IST EDR thus tends to have close agreement with MODIS IST for the middle temperature range, and have a higher bias when temperatures are low or high.
Comparisons of VIIRS IST EDR and MODIS IST for the Arctic only are shown in Figures 4 and 5.The bias and RMSE nobias for all cases, and for different temperature bins, have similar values and characteristics as those from all cases in the Antarctic.The RMSE is close to RMSE nobias because of the small bias, with RMSE of 0.98 K for all cases, and 1.04 K, 0.88 K, 0.95 K, 1.25 K, 1.31 K, and 1.18 K for MODIS IST in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K.The more rigorous matchups of the VIIRS IST EDR and MODIS IST with SNO for the Arctic yield an absolute bias of 0.16 K, RMSE nobias of 0.94 K, and RMSE of 0.95 Ks. (Figure 6).

Validation with KT-19 IST
Though the inter-comparison of the VIIRS IST and the MODIS IST provides valuable information about the VIIRS IST EDR quality, the uncertainty of each needs to be determined with suitable "truth" data.The uncertainty in the MODIS IST has previously been evaluated with limited in situ validation [36].The validation of the VIIRS IST EDR with KT-19 IST measurements provides the most accurate assessment of VIIRS IST EDR quality.The KT-19 sensor has been shown to accurately measure sea ice lead (fracture) temperatures of for leads 40 m wide or larger.With the rigorous procedures to collocate the KT-19 measurements with the VIIRS IST EDR, and with the elimination of the cloud contamination as detailed above, the comparison results give a high confidence assessment of the VIIRS IST EDR.With over 200 pairs of collocated VIIRS IST EDR and KT-19 IST over the Arctic and the Antarctic, the VIIRS minus KT-19 bias is ´0.34 K, with an RMSE of 0.92 K (0.98 K) without (with) the bias included (Figure 7).With the bias removed, this RMSE is smaller than the 1.2-1.3K value for the MODIS surface air temperature comparison in [12].However, it should be noted that the IceBridge KT-19 measurements are concentrated in the springtime, when the surface temperatures are in the moderate range and the VIIRS IST EDR performs best.More measurements from other seasons, when the IST is near the melting point or is very low, are desirable.

Comparison with Drifting Buoy 2 m Air Temperature
The VIIRS IST EDR was also compared to the 2 m air temperature measured by drifting ice buoys in the Arctic (Figure 8).The results show a bias of ´3.2 K, RMSE nobias of 3.7, and RMSE of 4.9 K, which means the VIIRS IST EDR is lower (colder) than the 2 m air temperature by over 3 K.The bias absolute values increase with increasing VIIRS IST.

Comparison with NCEP/NCAR Surface Air Temperature
Surface air temperature is assimilated in most of reanalyses and is also an output parameter [34,37].However, the surface skin temperature over land and sea ice is not assimilated.The VIIRS IST EDR was compared to the surface air temperature from the NCEP/NCAR reanalysis [34].Note that an additional comparison between the VIIRS IST and the NCEP-DOE Reanalysis 2 dataset [38] was also performed, but in our study the differences between the VIIRS IST and the NCEP/NCAR reanalysis were smaller than those from the VIIRS IST and NCEP-DOE comparisons.The results of the NCEP/NCAR comparison (Figures 9 and 10) show an overall bias of ´3.41 K and an RMSE nobias of 6.04 K, so the VIIRS IST EDR is 3.41 K lower (colder) than the surface air temperature from the NCEP/NCAR reanalysis.On average, the RMSE nobias is larger than that of the drifting ice buoy comparison.The VIIRS IST has a minimum cold bias and RMSE nobias when VIIRS IST is between 270 and 275 K, opposite to the results from the comparison with drifting ice buoys.This cold bias generally increases with decreasing IST except for the IST range between 260 and 270 K.The maximum cold bias and RMSE nobias occur for IST between 213 and 230 K.The RMSE is 6.94 K for all cases, and 9.34 K, 7.35 K, 6.83 K, 6.87 K, 6.25 K, and 3.01 K for NCEP/NCAR surface air temperature in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K.   Radionov et al. [39] showed that the average difference between the snow surface temperature and the 2 m air temperature is ´0.7 K (surface is colder) during the June-August warm season, and ´0.6 K in cold season from October to May.The annual mean difference in their measurements was ´0.1 K.The largest differences are between ´4.5 and ´6.9 K from November to May, usually under clear sky conditions, and +2.1 to +5.5 in the warmer part of the year.Given that the VIIRS IST EDR is for clear-sky conditions only, differences larger than the averages reported in the Radionov et al., study would be expected.Nevertheless, with thousands of cases in the comparisons of VIIRS, buoy, and reanalysis temperatures, a VIIRS cold bias of over 3 K may not be solely the result of air-skin temperature differences.
The large VIIRS IST cold bias, or reanalysis and buoy warm bias, may also be attributable to differences in the spatial coverage of VIIRS, drifting buoys, and the NCEP/NCAR reanalysis, and therefore the heterogeneity of the observed IST field.Due to the coarser resolution, the reanalysis is unable to resolve some of the complex IST features.The reanalysis could be smoothing out extreme values caused by features such as leads and polynyas that create small-scale temperature variations.A factor that could also explain the bias with drifting buoys is that some drifting buoys experience heating from sunlight, and as a result their temperatures may be biased high (warm).This would also cause a warm bias in the reanalysis, as buoy temperatures and pressures are assimilated [34].Another third possibility is that cloud contamination, especially from high cirrus clouds, results in a low retrieved VIIRS IST.

Summary and Conclusions
Monitoring the surface temperature is important for climate studies globally, but it deserves special attention in the Polar Regions because of sparse surface observations, an extreme temperature range, and the unique surface and atmospheric characteristics.The VIIRS sensor onboard the Suomi NPP and future Joint Polar Satellite System (JPSS) satellites provides an IST product with higher spatial resolution and wider coverage than heritage sensors such as AVHRR and MODIS with excellent radiometric quality.This study provided an assessment of the VIIRS IST EDR accuracy through comparisons with in situ and aircraft measurements, another satellite product, and a climate reanalysis.
Validation with KT-19 IST measurements from NASA Operation IceBridge flights over the Arctic and Antarctic in 2012-2014 yields a VIIRS IST cold bias of ´0.34 K and an RMSE nobias of 0.92 K.An intercomparison with MODIS IST from 2012 to 2015 gives a VIIRS IST bias of ´0.12 K and RMSE nobias of 1.01 K.The bias and RMSE nobias values show a dependence on temperature, with a positive bias for the warmer part of the range (270-275 K) and negative biases for lower temperatures.However, the bias and RMSE nobias values are relatively small in all temperature bins.A more rigorous simultaneous nadir overpass comparison shows an absolute bias of 0.16 K and an RMSE nobias of 0.94 K.
Comparisons of the VIIRS IST with 2 m air temperatures show a larger cold bias.Validation with drifting ice buoy 2 m air temperature yields a VIIRS IST bias of ´3.2 K and RMSE nobias of 3.7 K. Comparisons to NCEP/NCAR Reanalysis surface air temperature produce a similar bias of ´3.41 K and RMSE nobias of 6.04 K.
These validation and evaluation results demonstrate the high quality of the VIIRS IST EDR as indicated by the absolute bias and RMSE values for the IceBridge KT-19 measurement comparisons, where the KT-19 provide the best "truth" data for skin temperature.The results also show that the science quality of the VIIRS IST is similar to that of the MODIS IST.Based on the IceBridge and MODIS results, the VIIRS IST uncertainty (RMSE) meets or exceeds the JPSS system requirement of 1.0 K.The product can therefore be considered useful for meteorological and climatological applications.Furthermore, while we are not aware of any numerical weather prediction centers or climate reanalysis projects that assimilate satellite-derived ice surface temperature, the demonstrated accuracy of the IST product and the relatively large difference found in the comparison with the reanalysis surface temperatures provide motivation for research into its use in numerical models.

Figure 1 .
Figure 1.Ice surface temperature (IST) composite of the Visible Infrared Imaging Radiometer Suite (VIIRS) IST Environmental Data Record (EDR) from all overpasses over the Arctic on 1 March 2015.

Figure 2 . 6 Figure 2 .
Figure 2. Histogram of ice surface temperature differences between Suomi National Polar-orbiting Partnership (S-NPP) VIIRS and MODIS (Aqua and Terra) in the Arctic and Antarctic for all cases from August 2012 to July 2015.

Figure 3 .
Figure 3. Histogram of ice surface temperature differences between NPP VIIRS and MODIS (Aqua and Terra) in the Arctic and Antarctic from August 2012 to July 2015 for cases with MODIS IST in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K. Measurement bias and RMSEnobias are indicated for each bin.

Figure 3 .
Figure 3. Histogram of ice surface temperature differences between NPP VIIRS and MODIS (Aqua and Terra) in the Arctic and Antarctic from August 2012 to July 2015 for cases with MODIS IST in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K. Measurement bias and RMSE nobias are indicated for each bin.

Figure 4 .
Figure 4. Histogram of ice surface temperature differences between NPP VIIRS and MODIS (Aqua and Terra) in the Arctic only for all cases from August 2012 to July 2015.

Figure 5 .
Figure 5. Histogram of ice surface temperature differences between NPP VIIRS and MODIS (Aqua and Terra) in the Arctic only from August 2012 to July 2015 for cases with MODIS ice surface temperature in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K. Measurement bias and RMSE nobias are indicated for each bin.

Figure 6 .
Figure 6.Scatter plot of IST from 20 NPP VIIRS and MODIS Aqua simultaneous nadir overpasses (SNO) over the Arctic from March 2013 to April 2014, with an overall bias of 0.16 K, RMSE nobias of 0.94 K, and RMSE of 0.95 K.

Figure 8 .
Figure 8. Scatterplot of surface air temperature from Arctic drifting ice buoys and VIIRS IST from August 2012 to June 2014.The thick line is the 1-to-1 ratio line; the thin line as the linear regression.

Figure 9 .
Figure 9. Histogram of ice surface temperature differences between VIIRS IST and National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) surface air temperature in the Arctic for all cases from August 2012 to July 2015.

Figure 10 .
Figure 10.Histogram of ice surface temperature differences between VIIRS IST and NCEP/NCAR surface air temperature in the Arctic from August 2012 to July 2015 for cases with NCEP/NCAR surface air temperature in the ranges 213-230 K, 230-240 K, 240-250 K, 250-260 K, 260-270 K, and 270-275 K. Measurement bias and RMSE nobias are indicated for each bin.

Table 1 .
Validation datasets used in this study.