1. Introduction
Many attempts have been made to estimate land surface temperature (LST) from satellite observations at various temporal and spatial scales (e.g., Becker and Li 1990, 1995; McFarland et al. 1990; Sobrino et al. 1994; Wan and Li 1997; Schmugge et al. 1998; Snyder et al. 1998; Sun and Pinker 2003, 2005; Pinker et al. 2009; Yu et al. 2005, 2009; Duan et al. 2014; Ren et al. 2018; Zheng et al. 2019). The objective of these studies varied from validation of climate models to modeling evapotranspiration of land surface, energy and water balance and draft detection (e.g., Fuchs 1990; Anderson et al. 2007, 2011; Kustas and Anderson 2009; Li et al. 2009; Karnieli et al. 2010; Tomlinson et al. 2011; Wang et al. 2014). Only in few previous studies attention has been paid to the diurnal cycle. In one such early attempt, Ignatov and Gutman (1999) used the 3 hourly International Satellite Cloud Climatology Project (ISCCP) data (at 280 km resolution) (Rossow and Schiffer 1991) in combination with ground observations, to derive the monthly mean diurnal cycle in surface temperature over land, suitable for global climate studies. Duan et al. (2014) tried the same using high-spatial-resolution data from clear-sky MODIS observations while Inamdar et al. (2008) disaggregated the diurnal cycle of LST at the Geostationary Operational Environment Satellite (GOES) pixel scale to that at the MODIS pixel scale. Wang and Prigent (2020) compared the diurnal variations of LSTs for a full season in 2010 using reanalysis results from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA-Interim and ERA5), infrared satellite observations from an updated ISCCP product (Young et al. 2018), the Spinning Enhanced Visible and Infrared Imager (SEVIRI), and in situ measurements. They found that SEVIRI had closer agreement with the in situ measurements than the other products with a bias often less than ±2 K. Over snow or in arid areas ISCCP had more systematic errors than the other products; both reanalyses had higher (lower) estimations at nighttime (daytime) than the in situ measurements. To explore such variability, there is a need to establish a consistent and seamless long-term global record of land surface properties, which requires to homogenize satellite observations from several sources (e.g., Pinheiro et al. 2006; Susskind and Blaisdell 2008; Seemann et al. 2008; Anderson et al. 2011; Hulley and Hook 2011; Hulley et al. 2018; Ermida et al. 2017, 2018, 2020). For instance, Scarino et al. (2013, 2017) used a single-channel thermal-infrared (TIR) method to retrieve surface LST under clear-sky from geostationary-Earth-orbit (GEO) and low-Earth-orbit (LEO) satellite imagers. They used an empirically adjusted theoretical model of angular anisotropy (Vinnikov et al. 2012) to improve the satellite LST retrievals. They demonstrated that the application of the anisotropic correction yields reduced mean bias and improved precision for GOES-13 LST relative to independent Moderate Resolution Imaging Spectroradiometer (MYD11_L2) retrievals and also against the Atmospheric Radiation Measurement (ARM) program ground observations (Ackerman et al. 2016). Geostationary satellite data at high temporal and spatial resolution could provide a detailed depiction of the LST diurnal cycle.
In the present study we use satellite observations of the highest practical spatial (5-km) and temporal (hourly) resolutions and formulate a framework for deriving LST from both GOES-East (GOES-E) and GOES-West (GOES-W) utilizing a recently developed inference scheme (Pinker et al. 2019) that inherently accounts for atmospheric anisotropy and utilizes recent auxiliary information from the Modern-Era Retrospective Analysis for Research and Applications (MERRA-2) (Gelaro et al. 2017). We report on results obtained during the period (2004–09) over the United States and we investigate the potential of the satellite products to reproduce the diurnal cycle and diurnal temperature range (DTR) as observed by ground observations. The methodology is described in section 2; results are discussed in section 3; summary and discussion are provided in section 4.
2. Approach
a. Review of retrieval methodology
The GOES systems (East and West) provide readily available data at high temporal frequency (every half hour) at continental-scale coverage (https://www.bou.class.noaa.gov/release/index.htm). The imager scan sectors in routine mode is shown in Table 1. The present study is based on observations at 15 and 30 min after each hour from GOES-12 (GOES-E) and GOES-10/GOES-11 (GOES-W), respectively.
Selected GOES Imager scan sectors in routine mode.


Typically, the GOES imager includes five spectral channels (one visible, four infrared, Table 2). Before deriving the LST, all channels except 3 and 6 are calibrated (Gunshor et al. 2009; Weinreb et al. 2007) and used in a cloud screening algorithm (Pinker et al. 2019).
Summary of GOES imager channels.


b. Issue of GOES-E and GOES-W overlap
Figures 1a and 1b show the GOES-W coverage used and a calibrated image. For GOES-W, the images are of different size and dimension than those for GOES-E. They need to be regenerated to allow the use of existing screening code and ancillary information in the formats used for GOES-E. This is done using “nearest neighbor” interpolation method. The reprojected GOES-W image to GOES-E domain is shown in Fig. 1c. To run the LST retrieval algorithm for GOES-W, there is also a need to prepare ancillary data, such as surface type, snow-free channel-1 radiance and a pixel position data information for gridding in the same dimension and location as the satellite images of GOES-W.

(a) Band 2 image for 0030 UTC 1 May 2000 (raw counts). The X–Y axis shows pixel lines and elements; (b) calibrated band-2 brightness temperature for 0030 UTC 1 May 2000; (c) reprojected GOES-W to GOES-E.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

(a) Band 2 image for 0030 UTC 1 May 2000 (raw counts). The X–Y axis shows pixel lines and elements; (b) calibrated band-2 brightness temperature for 0030 UTC 1 May 2000; (c) reprojected GOES-W to GOES-E.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
(a) Band 2 image for 0030 UTC 1 May 2000 (raw counts). The X–Y axis shows pixel lines and elements; (b) calibrated band-2 brightness temperature for 0030 UTC 1 May 2000; (c) reprojected GOES-W to GOES-E.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
c. Data used for evaluation
MOD11
The LST retrievals from MOD11 version 6 land surface temperature and emissivity product (Wan 2014) are used to evaluate the GOES-based LST estimates at two time scales (instantaneous and daily). There are two algorithms used to process the MODIS LST: the generalized split-window LST algorithm and the day–night LST algorithm (Wan and Dozier 1996; Wan and Li 1997; Wang et al. 2014). The data are produced as a series of nine products. We use both MOD11_L2 data (DOI: 10.5067/MODIS/MOD11_L2.006) and MOD11C3 data (DOI: 10.5067/MODIS/MOD11C3.006). The MOD11_L2 version 6 swath product (MOD11L2 hereafter) provides LST and emissivity with pixel size of 1 km in 5-min temporal increments of satellite acquisition (https://lpdaac.usgs.gov/products/mod11_l2v006/). The MOD11C3 version 6 product (hereafter MOD11C3v6) provides monthly LST and emissivity values in a 0.05° (5.6 km at the equator) latitude–longitude climate modeling grid (CMG), which has a geographic grid with 7200 columns and 3600 rows representing the entire globe (https://lpdaac.usgs.gov/products/mod11c3v006/). The present study is restricted to the continental United States. These two products are also used in Pinker et al. (2019) for a comprehensive comparisons with the GOES-E LST.
d. In situ measurements
1) SURFRAD/BSRN
The Surface Radiation Budget Network (SURFRAD) (http://www.esrl.noaa.gov/gmd/grad/surfrad/) provides the best available, continuous, long-term measurements of surface radiation budget over the United States, which became the continental U.S. contingent of the international Baseline Surface Radiation Network (BSRN) (Ohmura et al. 1998; Augustine et al. 2005). The sites provide upwelling (
2) ARM Southern Great Plains
The Southern Great Plains (SGP) atmospheric observatory was the first field measurement site established by the ARM user facility (https://www.arm.gov/capabilities/observatories/sgp). This observatory is the world’s largest and most extensive climate research facility. The SGP site offers high-quality data and simulations at two tower levels (10 and 25 m). The instruments are a Heitronics GmbH KT19.85 Infrared Radiation Pyrometer (https://www.arm.gov/capabilities/instruments/irt) that measure radiances between 9.6 and 11.5 μm. The temperature resolution is ±0.45 K at 293 K. The uncertainty is 0.0244 K for a 0–1 V output range and 100 K span. The instruments are checked annually and the data quality is monitored by the ARM Data Quality Office. The surface skin temperature used in this paper is observed at 60-s intervals at the Central Facility (36.60°N, 97.48°W).
3) Oklahoma Mesonet
As detailed in Pinker et al. (2019) the Oklahoma Mesonet is a world-class network of environmental monitoring stations (http://www.mesonet.org). The network was designed and implemented by scientists from the University of Oklahoma (OU) and Oklahoma State University (OSU). It is an automated network of 120 stations across Oklahoma. About eighty nine of the Mesonet sites are installed with infrared temperature (IRT) sensors (8–14 μm) (Apogee Instruments, Inc.). This sensor is water resistant, designed for continuous outdoor use. Sensor accuracy is approximately ±0.2 K from 288 to 308 K and ±0.3 K from 278 to 318 K. The sensors are installed at a height of 1.5 m and have a field of view of a diameter circle of 0.5 m. A combination of automated and manual tests are applied using simultaneous soil and atmospheric measurements to intercompare observations and ensure that the skin temperature observations are of research quality (Fiebrich et al. 2003).
e. Uncertainties
The primary uncertainties of ground-based LST retrievals depend on the accuracy of the radiometric measurements and the emissivity estimates used in Eq. (2) (e.g., Hook et al. 2004, 2007, 2020; Augustine and Dutton 2013; Heidinger et al. 2013; Guillevic et al. 2014, 2017; Sobrino and Skoković 2016; Göttsche and Hulley 2012; Göttsche et al. 2013, 2016; Martin et al. 2019). The error in the in situ LST caused by the uncertainty of upwelling (±5 W m−2) and downwelling (±5 W m−2) radiometric measurement is less than 0.2 K (Martin et al. 2019). Error introduced by uncertainty in the broadband emissivity (<0.1) is less than 0.25 K (Heidinger et al. 2013). It should be noted that the above referenced errors are not negligible but neither are they a major source of uncertainty in LST estimation. Sobrino and Skoković (2016) claim that “the biggest part of uncertainty is due to inhomogeneity, which varies for each station and season. The other components have less influence, especially in seasons where the inhomogeneity is high.” The homogeneity of each site was previously assessed in Pinker et al. (2019). It was found that these sites showed a high degree of homogeneity for the period of 2000–08 over a 0.05° × 0.05° box.
f. Matchup issues
The method of spatial and temporal matching is described in Pinker et al. (2019) for GOES-E and also used for GOES-W. Briefly, since the spatial scale of the ground-based LST and the satellite LST product is different, we averaged the pixel values that fall within a 0.05° × 0.05° box, with the target site as a center. To reduce time difference between the in situ and satellite-based LST, we took the averages of the in situ LST observations that fall in ±15 min interval around the start of the scanning time for each GOES satellite. This selection is based on the scan duration (10–15 min) of the GOES satellite operation.
3. Results
a. Evaluation of LST from GOES-E and GOES-W
As discussed in sections 2a and 2b, required changes have been implemented in the retrieval schemes for applicability to GOES-W. The RTTOV algorithm is run for GOES-E and GOES-W for the years of 2004–09. The retrieved results from GOES-E for July 2004 at 0615 UTC and from GOES-W for July 2004 at 0630 UTC were compared against each other as shown in Fig. 2. Their differences and the frequency distribution of these differences are also illustrated. The frequency distribution of LST differences between the monthly means is also given. As seen, in the overlap area, the LST has similar spatial distribution patterns. The mean and median value of differences are small, 0.01 and 0.11 K, respectively.

(a) Monthly mean LST for July 2004 at 0615 UTC retrieved from GOES-E; (b) monthly mean LST for July at 0630 UTC retrieved from GOES-W; (c) LST differences between (a) and (b); (d) frequency distribution of LST differences.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

(a) Monthly mean LST for July 2004 at 0615 UTC retrieved from GOES-E; (b) monthly mean LST for July at 0630 UTC retrieved from GOES-W; (c) LST differences between (a) and (b); (d) frequency distribution of LST differences.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
(a) Monthly mean LST for July 2004 at 0615 UTC retrieved from GOES-E; (b) monthly mean LST for July at 0630 UTC retrieved from GOES-W; (c) LST differences between (a) and (b); (d) frequency distribution of LST differences.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
In Fig. 3 we show monthly mean LST for July 2004 at 1815 and 1830 UTC as retrieved from GOES-E and GOES-W, respectively, and their difference. The frequency distribution of LST differences is also given. As seen, the differences are higher than at 0600 UTC, the mean and median value of differences are 0.15 and 1.33 K, respectively. The objective of this evaluation was to provide information on the range of differences during extreme times of the day.

(a) Monthly mean LST for July 2004 at 1815 UTC retrieved from GOES-E; (b) monthly mean LST for July at 1830 UTC retrieved from GOES-W; (c) LST differences between (a) and (b); (d) frequency distribution of LST differences. Only grids/points that have at least 3 days of values have been used.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

(a) Monthly mean LST for July 2004 at 1815 UTC retrieved from GOES-E; (b) monthly mean LST for July at 1830 UTC retrieved from GOES-W; (c) LST differences between (a) and (b); (d) frequency distribution of LST differences. Only grids/points that have at least 3 days of values have been used.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
(a) Monthly mean LST for July 2004 at 1815 UTC retrieved from GOES-E; (b) monthly mean LST for July at 1830 UTC retrieved from GOES-W; (c) LST differences between (a) and (b); (d) frequency distribution of LST differences. Only grids/points that have at least 3 days of values have been used.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
b. Evaluation against SURFRAD/BSRN at DRA
We have evaluated GOES-E and GOES-W LST estimates at DRA for the period 2004–09, independently for daytime and nighttime. The scatterplots of the instantaneous satellite-based LST for GOES-E and GOES-W against DRA for both daytime and nighttime from 2004 to 2009 are shown in Fig. 4. As evident, the satellite estimates and the ground observations have high correlation (Corr ≈ 0.99). Overall, GOES-W has smaller bias than GOES-E for both daytime and nighttime. For daytime (upper panel, Fig. 4) the bias of GOES-W is about 0.31 K (0.1%) while GOES-E is about −1.16 K (0.4%). For nighttime, the bias is −1.25 K (0.4%) and −2.18 K (0.8%) for GOES-W and GOES-E, respectively. The standard deviation (std) of the two products are comparable to each other for both daytime and nighttime (std < 2 K). The root-mean-square error (RMSE) is less than 1%. One needs to note that the GOES-E instantaneous LST we used here are 15 min after hour while the GOES-W instantaneous data are 30 min after the hour. However, the averaging of the ground observations are around each satellite overpass time, respectively.

Evaluation of instantaneous (left) GOES-E- and (right) GOES-W-based LST estimates against the DRA station, independently for (top) daytime and (bottom) nighttime from 2004 to 2009.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Evaluation of instantaneous (left) GOES-E- and (right) GOES-W-based LST estimates against the DRA station, independently for (top) daytime and (bottom) nighttime from 2004 to 2009.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Evaluation of instantaneous (left) GOES-E- and (right) GOES-W-based LST estimates against the DRA station, independently for (top) daytime and (bottom) nighttime from 2004 to 2009.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
c. Comparison between GOES-E, GOES-W at instantaneous time scale over the ARM/SGPC1 site
The probability distribution of differences between GOES-E and GOES-W-based LST retrievals and ARM/SGPC1 in situ LST is shown in Fig. 5 and statistics are provided in Table 3. Most of the differences (>80%) fall within the interval of 1 std (~4 K). The correlation between GOES-E and ARM/SGPC1 and GOES-W and ARM/SGPC1 is high for both observational levels (Corr > 0.95). The bias for GOES-E at 10 and 25 m is about −0.74 K (0.25%) and −0.92 K (0.32%), respectively, while GOES-W is about 1.3 K lower.

Evaluation of GOES-E and GOES-W against ARM/SGPC1 at (left) 10 and (right) 25 m during 2004–09. Black solid: GOES-E; black dashed: GOES-W; red dotted: ±std for GOES-E; blue dotted: ±3 std for GOES-E; red dash–dotted: ±1 std for GOES-W; blue dash–dotted: ±3 std for GOES-W.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Evaluation of GOES-E and GOES-W against ARM/SGPC1 at (left) 10 and (right) 25 m during 2004–09. Black solid: GOES-E; black dashed: GOES-W; red dotted: ±std for GOES-E; blue dotted: ±3 std for GOES-E; red dash–dotted: ±1 std for GOES-W; blue dash–dotted: ±3 std for GOES-W.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Evaluation of GOES-E and GOES-W against ARM/SGPC1 at (left) 10 and (right) 25 m during 2004–09. Black solid: GOES-E; black dashed: GOES-W; red dotted: ±std for GOES-E; blue dotted: ±3 std for GOES-E; red dash–dotted: ±1 std for GOES-W; blue dash–dotted: ±3 std for GOES-W.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Statistics from intercomparison of LST retrieved from GOES-E and GOES-W at the ARM/SGPC1 site during 2004–09 as shown in Fig. 5.


d. Comparison between GOES-E, GOES-W at instantaneous time scale at Oklahoma Mesonet
The distribution of Oklahoma Mesonet sites used in current evaluation is shown (Fig. 7). The evaluations are carried out against all the stations for both daytime and nighttime in January and July during 2004–07. Results are presented in Figs. 6 and 7 . Both of the retrieved GOES LSTs have high correlation with the in situ data (Corr > 0.8); see Table 4. Results for GOES-E at both daytime and nighttime for January and July are comparable. GOES-W has smaller bias for July than January at both daytime and nighttime. For January, the GOES-W data have comparable bias with GOES-E, but relatively have a larger spread than GOES-E. For July, the bias of GOES-E data is much smaller than GOES-W and its spread is also slightly smaller than GOES-W.

Evaluation against Mesonet sites for daytime and nighttime in January and July over 2004–07.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Evaluation against Mesonet sites for daytime and nighttime in January and July over 2004–07.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Evaluation against Mesonet sites for daytime and nighttime in January and July over 2004–07.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Distribution of RMSE from intercomparison of LST retrieved from (left) GOES-E and (right) GOES-W at Mesonet sites during 2004–07.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Distribution of RMSE from intercomparison of LST retrieved from (left) GOES-E and (right) GOES-W at Mesonet sites during 2004–07.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Distribution of RMSE from intercomparison of LST retrieved from (left) GOES-E and (right) GOES-W at Mesonet sites during 2004–07.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Statistics from intercomparison of LST retrieved from GOES-E, GOES-W at Mesonet sites for daytime and nighttime in January and July over 2004–07.


Figure 7 shows the distribution of RMSE from intercomparison of LST retrieved from GOES-E, GOES-W at Mesonet sites during 2004–07. As seen for GOES-E, the RMSE of most of the sites is less than 1% of their mean LST (282.03–298.66 K). Only at two sites the RMSE is greater than 1.5% but less than 2%. The minimum value is 1.76 K and the maximum is 6.53 K. For GOES-W, at most of the sites the RMSE is greater than 1% but less than 1.5%. The minimum value is 1.92 K and the maximum is 6.76 K.
e. Comparison between GOES-E, GOES-W, and MODIS at instantaneous time scale
We have conducted an intercomparison between MOD11L2 and the two GOES LST estimates. The intercomparison approach requires accounting for differences in spatial resolution, view angle and overpass time between the satellites (Guillevic et al. 2014). Matchups represent coincident pairs of granules with respect to satellite overpass times and view angles. The “near miss” time spans are usually referred to as simultaneous nadir overpasses (SNOs) when nadir view angles are considered (Cao et al. 2004). First, the MOD11L2 data are regridded to 0.05° × 0.05°. The GOES-W LST has a 15 min time difference with GOES-E and MODIS. Figure 8 shows the LST retrievals from GOES-E and MOD11L2 at 1815 UTC and GOES-W at 1830 UTC 29 August 2004. The differences of their distributions are presented in Fig. 9. The statistics are based on all available grid points for this case as shown in Table 5.

(left) LST retrievals and (right) their distributions from (bottom) MOD11L2 and (top) GOES-E for 1815 UTC and (middle) GOES-W for 1830 UTC 29 Aug 2004.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

(left) LST retrievals and (right) their distributions from (bottom) MOD11L2 and (top) GOES-E for 1815 UTC and (middle) GOES-W for 1830 UTC 29 Aug 2004.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
(left) LST retrievals and (right) their distributions from (bottom) MOD11L2 and (top) GOES-E for 1815 UTC and (middle) GOES-W for 1830 UTC 29 Aug 2004.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Comparison of LST retrieval differences between (top) GOES-E and MOD11L2, (middle) GOES-W and MOD11, and (bottom) GOES-E and GOES-W at 1815 or 1830 UTC 29 Aug 2004. The statistical summaries are shown in Table 5.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Comparison of LST retrieval differences between (top) GOES-E and MOD11L2, (middle) GOES-W and MOD11, and (bottom) GOES-E and GOES-W at 1815 or 1830 UTC 29 Aug 2004. The statistical summaries are shown in Table 5.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Comparison of LST retrieval differences between (top) GOES-E and MOD11L2, (middle) GOES-W and MOD11, and (bottom) GOES-E and GOES-W at 1815 or 1830 UTC 29 Aug 2004. The statistical summaries are shown in Table 5.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Statistics from intercomparison among LST retrieved from GOES-E and GOES-W and from MOD11L2 product at 1815 or 1830 UTC 29 Aug 2004 as shown in Fig. 9.


As seen in the left panel of Fig. 8, the highest LST is close to the southeast coast for all of the three products. The maximum values for GOES-E, GOES-W, and MOD11L2 are 332.2, 334.7, and 332.7 K, respectively. The position of lowest LST for all three products are similarly located. The minimum values are 278.8, 269.0 and 280.9 K. As evident from the right panel of Fig. 8, most of the grids are in the range of 307.5–313.0 K for all products.
From Fig. 9 and Table 5, it is evident that both GOES LST products have high correlation with MOD11L2 (coefficients are over 0.85 using more than 40 000 grid points) while GOES LST products are higher than MOD11L2. The largest differences occur in Washington State. The bias for GOES-E and GOES-W is about 2.54 K (0.8%) and 2.03 K (0.7%), respectively. Correlation coefficient between GOES-W and GOES-E LST is 0.9 and the bias is about −0.51 K (0.2%).
f. Comparison between GOES-E, GOES-W, and MODIS at seasonal time scale
We also evaluate the ability of capturing the LST seasonal changes at DRA including the MOD11C3v6 product. Figure 10 shows the LST distribution for GOES_E, GOES-W, and MOD11C3v6 over 2004–09 for each month and their monthly mean values compared with DRA (Table 6).

Comparison of (left) GOES-E and (right) GOES-W LST at monthly time scale for 2004–09 with MOD11C3v6 and the DRA site. The mean values (stars) are shown in Table 6.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Comparison of (left) GOES-E and (right) GOES-W LST at monthly time scale for 2004–09 with MOD11C3v6 and the DRA site. The mean values (stars) are shown in Table 6.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Comparison of (left) GOES-E and (right) GOES-W LST at monthly time scale for 2004–09 with MOD11C3v6 and the DRA site. The mean values (stars) are shown in Table 6.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Both of the retrieved GOES LSTs have similar patterns of variability to those from the DRA site. MOD11C3v6 LSTs are higher than the DRA observations for most seasons; the annual mean LST for all years used is 294.5 K, which is 2.3 K higher than the value measured at DRA (292.2 K), but GOES estimates are much closer to the site value than the others. The annul mean LST from GOES-E is 290.6 K, which is 1.6 K lower than DRA; the value from GOES-W is 291.5 K, which is 0.7 K lower than DRA.
g. Climatology of LST
A 6 yr (2004–09) mean LST at 0.05° spatial resolution for January and July is shown in Fig. 11; statistics is presented in Table 7. As shown, for July and January, the LST distribution pattern of the two products are similar to each other. They have high correlation (Corr > 0.9) and very small bias (0.19 K for July and −0.76 K for January). The std and RMSE are larger in July than in January.

Climatology of LST for (left) July and (right) January averaged for 2004–09 as derived from (top) GOES-EAST and (middle) GOES-WEST and (bottom) their differences.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Climatology of LST for (left) July and (right) January averaged for 2004–09 as derived from (top) GOES-EAST and (middle) GOES-WEST and (bottom) their differences.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Climatology of LST for (left) July and (right) January averaged for 2004–09 as derived from (top) GOES-EAST and (middle) GOES-WEST and (bottom) their differences.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Intercomparison statistics of LST from GOES-E and GOES-W climatology during 2004–09 over the United States in areas of overlap in January and July.


Figure 12 shows the diurnal cycle of LST at DRA and at the ARM/SGPC1 sites from ground and satellite-based estimates. Closer agreement between the satellite and ground observations is seen from about noon to late afternoon. The difference between DRA observations and GOES estimates are around 0–2 K while at ARM/SGPC1 the differences are about 0.5–3 K; before noon, the differences are larger around 0–6 K. The GOES-E LST has shown a good agreement in depicting the diurnal cycle at other SURFRAD sites (Pinker et al. 2019). These findings indicate that the GOES-E and GOES-W estimates have the ability to represent well the diurnal variations of the LST. This is of extreme interest since most satellite-based estimates of LST use polar orbiters unable to depict the true diurnal cycle. A comprehensive analysis over the entire United States will be conducted independently.

Diurnal cycle of average LST based on data from 2004 to 2009 as observed at DRA, ARM/SGPC1 and derived from GOES-E and GOES-W.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1

Diurnal cycle of average LST based on data from 2004 to 2009 as observed at DRA, ARM/SGPC1 and derived from GOES-E and GOES-W.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
Diurnal cycle of average LST based on data from 2004 to 2009 as observed at DRA, ARM/SGPC1 and derived from GOES-E and GOES-W.
Citation: Journal of Atmospheric and Oceanic Technology 38, 4; 10.1175/JTECH-D-20-0086.1
4. Summary and discussion
Prospects for achieving global information from GEO satellites are improving. New development in geostationary systems aims to make them more similar to each other in terms of spatial resolution and the spectral characteristics of relevant channels. For instance, the GOES-R Advanced Baseline Imager (ABI) Sensor (Schmit et al. 2017) is similar to the Advanced Himawari Imager (AHI) 8 and 9 (Okuyama et al. 2015). Such improvements make it possible to establish a consistent and seamless long-term global record of land surface properties, which requires homogenization of satellite observations from several sources. However, the issue of merging of the observations in areas of overlap due to differences in viewing geometry and the impact of the anisotropy of surface emissivity (Minnis et al. 2004; Cuenca and Sobrino 2004; Pinheiro et al. 2006) still exists. As illustrated in Evan et al. (2007) in respect to the International Satellite Cloud Climatology Project (ISCCP) cloud data (Rossow and Schiffer 1999), the abrupt changes in global cloud amounts result from sudden changes in the geostationary viewing angles. They also point out that these observations are consistent with the theory that changes in the number of geostationary satellites are altering information on the global mean cloud amounts (Campbell 2004). As such it is critical to evaluate and optimize observations from different sources, a formidable task and as yet, not comprehensively addressed. An early attempt to address the merging issue from geostationary satellites in respect to surface radiative fluxes is described in Zhang et al. (2007). They introduced a semiempirical orthogonal function (EOF) iteration scheme for homogenizing the fluxes. On the average, the latter reduced the RMSE in the fluxes as compared to ground observations by about 2–3 W m−2. The newly revived interest in the topic of merging GEO satellite observations can be attributed to progress made to improve satellite observations. The new capabilities of geostationary satellites (the European SEVIRI, ABI on GOES-R, Himawari-8/9) provide new capabilities to derive high-resolution (temporal and spatial) climatic parameters such as LST. The satellites have similar channels (GOES-R and Himawari-8/9), which reduces some of the problems in merging data from different platforms.
In this paper, we describe an approach to derive and evaluate high-temporal- and high-spatial-resolution information on LST from the GOES satellites across multiple missions and multiple sensors with overlapping coverage. Specifically, we have implemented the RTTOV radiative transfer model adjusted for channel 4 (10.7 μm) of GOES-E and GOES-W with the MERRA-2 atmospheric profiles and the CAMEL emissivity product to derive a 6 yr record (2004–09) of LST. The data are produced at 0.05° spatial resolution at hourly time scale and evaluated for the period of 2004–09 against best available ground observations and an independent well established product from MODIS. We report results of evaluation at instantaneous time scale as well as averaged over a month. We found that monthly mean differences in LST for July at 0600 UTC retrieved from GOES-E and GOES-W in terms of mean and median values were 0.01 and 0.11 K, respectively, but were higher at 1800 UTC, with respective mean and median differences of 0.15 and 1.33 K. The performance of GOES-based LST is comparable to the MODIS product and is in good agreement with in situ data. As such they are of sufficient quality to represent seasonal and diurnal variability and climatological characteristics of LST over the United States.
Another important application of the DTR as provided by the data that generated in this study is related to issues of climate change that currently are investigated by changes in mean temperature. According to several studies (Easterling et al. 1997; Dai et al. 1999; Davy et al. 2017) the DTR has been decreasing worldwide since the 1950s. According to Easterling et al. (1997) the decrease is primarily due to the increase in minimum temperature. The various possible causes for the increase in the minimum temperature have been discussed in numerous papers (Karl et al. 1993; Braganza et al. 2004; Stone and Weaver 2002, 2003; Dai et al. 1999). While the above studies used information on air temperature only, the LST data generated allow us to take a fresh look at this issue. As demonstrated, our results are of sufficient quality to represent seasonal and diurnal variability and climatological characteristics of LST over the United States as well as DTR.
Our findings can serve as guidelines for developing a strategy for global coverage from geostationary satellites. The methodology we use minimizes differences between satellites in areas of overlap. The proposed semiempirical corrections for viewing geometry that have been used in a previous study (Scarino et al. 2017) had only small impact on our results. Scarino et al. (2017) found a larger impact of the angular corrections on the retrieved LSTs; however, their experiment is based on independent observations and methodology. Possibly, findings can be also impacted by the limited database used for developing the angular corrections. A controlled experiment using same observations but independent methodology could shed more light on this issue in future.
Acknowledgments
We acknowledge the MERRA-2 data as provided by the Global Modeling and Assimilation Office (GMAO). Information from the MOD11_L2: MODIS/Terra land surface temperature and emissivity 5-min L2 swath 1 km V006 database (Zhengming Wan, PI), (https://search.earthdata.nasa.gov/search/granules/collectiondetails?p=C194001236-LPDAAC_ECS&m=-26.4375!136.6875!0!1!0!0%2C2&tl=1515438649!4!!&q=MOD11_L2%20V006) was provided under the courtesy of the NASA EOSDIS Land Processes Distributed Active Archive Center (LP DAAC), USGS/Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota. GOES data were obtained from the NOAA Comprehensive Large Array data Stewardship System (CLASS) (https://www.bou.class.noaa.gov/saa/products/search?sub_id=0&datatype_family=GVAR_IMG&submit.x=25&submit.y=8). The BSRN/SURFRAD data were provided by the NOAA Earth System Research Laboratory, Global Monitoring Division (https://www.esrl.noaa.gov/gmd/grad/surfrad/). The ARM data were obtained from the Radiation Measurement (ARM) Climate Research 10.5067/MEaSUREs/LSTE/CAM5K30EM.001. The Mesonet data were provided by https://weather.ok.gov/index.php/site/about/data_access_and_pricing (McPherson et al. 2007). The team efforts in generating and providing all the required data is greatly appreciated. This research was funded by Grant NNH12ZDA001N-MEASURES from NASA to JPL. Part of the research described in this paper was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. Government sponsorship acknowledged. We thank the anonymous reviewers for constructive comments that helped to improve the manuscript and to the editor for overseeing the review process.
Data availability statement
The data are being prepared for publication on a website. Prior to formal publication can be obtained upon request.
REFERENCES
Ackerman, T. P., T. S. Cress, W. Ferrell, J. H. Mather, and D. D. Turner, 2016: The programmatic maturation of the ARM program. The Atmospheric Radiation Measurement (ARM) Program: The First 20 Years, Meteor. Monogr., No. 57, Amer. Meteor. Soc., https://doi.org/10.1175/AMSMONOGRAPHS-D-15-0054.1.
Aires, F., C. Prigent, and W. B. Rossow, 2004: Temporal interpolation of global surface skin temperature diurnal cycle over land under clear and cloudy conditions. J. Geophys. Res., 109, D04313, https://doi.org/10.1029/2003JD003527.
Anderson, M. C., J. M. Norman, J. R. Mecikalski, J. A. Otkin, and W. P. Kustas, 2007: Climatological study of evapotranspiration and moisture stress across the continental United States: 1. Model formulation. J. Geophys. Res., 112, D10117, https://doi.org/10.1029/2006JD007506.
Anderson, M. C., and Coauthors, 2011: Mapping daily evapotranspiration at field to continental scales using geostationary and polar orbiting satellite imagery. Hydrol. Earth Syst. Sci., 15, 223–239, https://doi.org/10.5194/hess-15-223-2011.
Augustine, J. A., and E. G. Dutton, 2013: Variability of the surface radiation budget over the United States from 1996 through 2011 from high-quality measurements. J. Geophys. Res. Atmos., 118, 43–53, https://doi.org/10.1029/2012JD018551.
Augustine, J. A., G. B. Hodges, C. R. Cornwall, J. J. Michalsky, and C. I. Medina, 2005: An update on SURFRAD: The GCOS Surface Radiation Budget network for the continental United States. J. Atmos. Oceanic Technol., 22, 1460–1472, https://doi.org/10.1175/JTECH1806.1.
Barkan, J., H. Shafir, and P. Alpert, 2020: Multi-factor analysis of DTR variability over Israel in the sea/desert border. Theor. Appl. Climatol., 139, 287–295, https://doi.org/10.1007/s00704-019-02958-x.
Becker, F., and Z. Li, 1990: Towards a local split window method over land surfaces. Int. J. Remote Sens., 11, 369–393, https://doi.org/10.1080/01431169008955028.
Becker, F., and Z. Li, 1995: Surface temperature and emissivity at various scales: Definition, measurement and related problems. Remote Sens. Rev., 12, 225–253, https://doi.org/10.1080/02757259509532286.
Borbas, E., G. Hulley, M. Feltz, R. Knuteson, and S. Hook, 2018: The Combined ASTER MODIS Emissivity over Land (CAMEL) part 1: Methodology and high spectral resolution application. Remote Sens., 10, 643, https://doi.org/10.3390/rs10040643.
Braganza, K., D. J. Karoly, and J. M. Arblaster, 2004: Diurnal temperature range as an index of global climate change during the twentieth century. Geophys. Res. Lett., 31, L13217, https://doi.org/10.1029/2004GL019998.
Campbell, G., 2004: View angle dependence of cloudiness and the trend in ISCCP cloudiness. 13th Conf. on Satellite Meteorology and Oceanography, Norfolk, VA, Amer. Meteor. Soc., P6.7, https://ams.confex.com/ams/13SATMET/techprogram/paper_79041.htm.
Cao, C., M. Weinreb, and H. Xu, 2004: Predicting simultaneous nadir overpasses among polar-orbiting meteorological satellites for the intersatellite calibration of radiometers. J. Atmos. Oceanic Technol., 21, 537–542, https://doi.org/10.1175/1520-0426(2004)021<0537:PSNOAP>2.0.CO;2.
Carlson, T. N., 1986: Regional scale estimates of surface moisture availability and thermal inertia using remote thermal measurements. Remote Sens. Rev., 1, 197–247, https://doi.org/10.1080/02757258609532069.
Cuenca, J., and J. Sobrino, 2004: Experimental measurements for studying angular and spectral variation of thermal infrared emissivity. Appl. Opt., 43, 4598–4602, https://doi.org/10.1364/AO.43.004598.
Dai, A., and K. E. Trenberth, 2004: The diurnal cycle and its depiction in the Community Climate System Model. J. Climate, 17, 930–951, https://doi.org/10.1175/1520-0442(2004)017<0930:TDCAID>2.0.CO;2.
Dai, A., K. E. Trenberth, and T. R. Karl, 1999: Effects of clouds, soil moisture, precipitation, and water vapor on diurnal temperature range. J. Climate, 12, 2451–2473, https://doi.org/10.1175/1520-0442(1999)012<2451:EOCSMP>2.0.CO;2.
Davy, R., I. Esau, A. Chernokulsky, S. Outten, and S. Zilitinkevich, 2017: Diurnal asymmetry to the observed global warming. Int. J. Climatol., 37, 79–93, https://doi.org/10.1002/joc.4688.
Duan, S.-B., Z.-L. Li, B.-H. Tang, H. Wu, R. Tang, Y. Bi, and G. Zhou, 2014: Estimation of diurnal cycle of land surface temperature at high temporal and spatial resolution from clear-sky MODIS data. Remote Sens., 6, 3247–3262, https://doi.org/10.3390/rs6043247.
Easterling, D. R., and Coauthors, 1997: Maximum and minimum temperature trends for the globe. Science, 277, 364–367, https://doi.org/10.1126/science.277.5324.364.
Ermida, S. L., C. C. DaCamara, I. F. Trigo, A. C. Pires, D. Ghent, and J. Remedios, 2017: Modelling directional effects on remotely sensed land surface temperature. Remote Sens. Environ., 190, 56–69, https://doi.org/10.1016/J.RSE.2016.12.008.
Ermida, S. L., I. F. Trigo, C. C. DaCamara, and A. C. Pires, 2018: A methodology to simulate LST directional effects based on parametric models and landscape properties. Remote Sens., 10, 1114, https://doi.org/10.3390/rs10071114.
Ermida, S. L., I. F. Trigo, G. Hulley, and C. C. DaCamara, 2020: A multi-sensor approach to retrieve emissivity angular dependence over desert regions. Remote Sens. Environ., 237, 111559, https://doi.org/10.1016/j.rse.2019.111559.
Evan, A. T., A. K. Heidinger, and D. J. Vimont, 2007: Arguments against a physical long-term trend in global ISCCP cloud amounts. Geophys. Res. Lett., 34, L04701, https://doi.org/10.1029/2006GL028083.
Feltz, M., E. Borbas, R. Knuteson, G. Hulley, and S. Hook, 2018: The Combined ASTER MODIS Emissivity over Land (CAMEL) part 2: Uncertainty and validation. Remote Sens., 10, 664, https://doi.org/10.3390/rs10050664.
Fiebrich, C. A., J. E. Martinez, J. A. Brotzge, and J. B. Basara, 2003: The Oklahoma Mesonet’s skin temperature network. J. Atmos. Oceanic Technol., 20, 1496–1504, https://doi.org/10.1175/1520-0426(2003)020<1496:TOMSTN>2.0.CO;2.
Fuchs, M., 1990: Infrared measurement of canopy temperature and detection of plant water stress. Theor. Appl. Climatol., 42, 253–261, https://doi.org/10.1007/BF00865986.
Gelaro, R., and Coauthors, 2017: The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). J. Climate, 30, 5419–5454, https://doi.org/10.1175/JCLI-D-16-0758.1.
Göttsche, F.-M., and G. C. Hulley, 2012: Validation of six satellite-retrieved land surface emissivity products over two land cover types in a hyper-arid region. Remote Sens. Environ., 124, 149–158, https://doi.org/10.1016/j.rse.2012.05.010.
Göttsche, F.-M., F.-S. Olesen, and A. Bork-Unkelbach, 2013: Validation of land surface temperature derived from MSG/SEVIRI with in situ measurements at Gobabeb, Namibia. Int. J. Remote Sens., 34, 3069–3083, https://doi.org/10.1080/01431161.2012.716539.
Göttsche, F.-M., F.-S. Olesen, I. F. Trigo, A. Bork-Unkelbach, and M. A. Martin, 2016: Long term validation of land surface temperature retrieved from MSG/SEVIRI with continuous in-situ measurements in Africa. Remote Sens., 8, 410, https://doi.org/10.3390/rs8050410.
Guillevic, P. C., J. L. Privette, B. Coudert, M. A. Palecki, J. Demarty, C. Ottlé, and J. A. Augustine, 2012: Land surface temperature product validation using NOAA’s surface climate observation networks—Scaling methodology for the Visible Infrared Imager Radiometer Suite (VIIRS). Remote Sens. Environ., 124, 282–298, https://doi.org/10.1016/j.rse.2012.05.004.
Guillevic, P. C., and Coauthors, 2014: Validation of land surface temperature products derived from the Visible Infrared Imaging Radiometer Suite (VIIRS) using ground-based and heritage satellite measurements. Remote Sens. Environ., 154, 19–37, https://doi.org/10.1016/j.rse.2014.08.013.
Guillevic, P. C., and Coauthors, 2017: Land surface temperature product validation best practice protocol, version 1.0. Land Product Validation Subgroup Publ., 60 pp.
Gunshor, M. M., T. J. Schmit, W. P. Menzel, and D. C. Tobin, 2009: Intercalibration of broadband geostationary imagers using AIRS. J. Atmos. Oceanic Technol., 26, 746–758, https://doi.org/10.1175/2008JTECHA1155.1.
Heidinger, A. K., I. Laszlo, C. C. Molling, and D. Tarpley, 2013: Using SURFRAD to verify the NOAA single-channel land surface temperature algorithm. J. Atmos. Oceanic Technol., 30, 2868–2884, https://doi.org/10.1175/JTECH-D-13-00051.1.
Hook, S. J., and Coauthors, 2004: In-flight validation and recovery of water surface temperature with Landsat-5 thermal infrared data using an automated high-altitude lake validation site at Lake Tahoe. IEEE Trans. Geosci. Remote Sens., 42, 2767–2776, https://doi.org/10.1109/TGRS.2004.839092.
Hook, S. J., R. G. Vaughan, H. Tonooka, and S. G. Schladow, 2007: Absolute radiometric in-flight validation of mid infrared and thermal infrared data from ASTER and MODIS on the Terra spacecraft using the Lake Tahoe, CA/NV, USA, automated validation site. IEEE Trans. Geosci. Remote Sens., 45, 1798–1807, https://doi.org/10.1109/TGRS.2007.894564.
Hook, S. J., K. Cawse-Nicholson, J. Barsi, R. Radocinski, G. C. Hulley, W. R. Johnson, G. Rivera, and B. Markham, 2020: In-flight validation of the ECOSTRESS, Landsats 7 and 8 thermal infrared spectral channels using the Lake Tahoe CA/NV and Salton Sea CA automated validation sites. IEEE Trans. Geosci. Remote Sens., 58, 1294–1302, https://doi.org/10.1109/TGRS.2019.2945701.
Hulley, G. C., and S. Hook, 2011: Generating consistent land surface temperature and emissivity products between ASTER and MODIS data for Earth science research. IEEE Trans. Geosci. Remote Sens., 49, 1304–1315, https://doi.org/10.1109/TGRS.2010.2063034.
Hulley, G. C., S. Veraverbeke, and S. Hook, 2014: Thermal-based techniques for land cover change detection using a new dynamic MODIS multispectral emissivity product (MOD21). Remote Sens. Environ., 140, 755–765, https://doi.org/10.1016/j.rse.2013.10.014.
Hulley, G. C., N. K. Malakar, T. Islam, and R. J. Freepartner, 2018: NASA’s MODIS and VIIRS land surface temperature and emissivity products: A long-term and consistent Earth system data record. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens., 11, 522–535, https://doi.org/10.1109/JSTARS.2017.2779330.
Idso, S. B., R. D. Jackson, and R. J. Reginato, 1976: Compensating for environmental variability in the thermal inertia approach to remote sensing of soil moisture. J. Appl. Meteor., 15, 811–817, https://doi.org/10.1175/1520-0450(1976)015<0811:CFEVIT>2.0.CO;2.
Ignatov, A., and G. Gutman, 1999: Monthly mean diurnal cycles in surface temperature over and land for global climate studies. J. Climate, 12, 1900–1910, https://doi.org/10.1175/1520-0442(1999)012<1900:MMDCIS>2.0.CO;2.
Inamdar, A. K., A. French, S. Hook, G. Vaughan, and W. Luckett, 2008: Land surface temperature retrieval at high spatial and temporal resolutions over the southwestern United States. J. Geophys. Res., 113, D07107, https://doi.org/10.1029/2007JD009048.
Karl, T. R., and Coauthors, 1993: A new perspective on recent global warming: Asymmetric trends of daily maximum and minimum temperature. Bull. Amer. Meteor. Soc., 74, 1007–1023, https://doi.org/10.1175/1520-0477(1993)074<1007:ANPORG>2.0.CO;2.
Karnieli, A., N. Agam, R. T. Pinker, M. Anderson, M. L. Imhoff, G. G. Gutman, N. Panov, and A. Goldberg, 2010: Use of NDVI and land surface temperature for drought assessment: Merits and limitations. J. Climate, 23, 618–633, https://doi.org/10.1175/2009JCLI2900.1.
Kustas, W., and M. Anderson, 2009: Advances in thermal infrared remote sensing for land surface modeling. Agric. For. Meteor., 149, 2071–2081, https://doi.org/10.1016/j.agrformet.2009.05.016.
Li, Z.-L., R. Tang, Z. Wan, Y. Bi, C. Zhou, B. Tang, G. Yan, and X. Zhang, 2009: A review of current methodologies for regional evapotranspiration estimation from remotely sensed data. Sensors, 9, 3801–3853, https://doi.org/10.3390/s90503801.
Malakar, N. K., G. C. Hulley, S. J. Hook, K. Laraby, M. Cook, and J. R. Schott, 2018: An operational land surface temperature product for Landsat thermal data: Methodology and validation. IEEE Trans. Geosci. Remote Sens., 56, 5717–5735, https://doi.org/10.1109/TGRS.2018.2824828.
Martin, M., and Coauthors, 2019: Comprehensive in situ validation of five satellite land surface temperature data sets over multiple stations and years. Remote Sens., 11, 479, https://doi.org/10.3390/rs11050479.
Matricardi, M., and R. Saunders, 1999: Fast radiative transfer model for simulation of Infrared Atmospheric Sounding Interferometer radiances. Appl. Opt., 38, 5679–5691, https://doi.org/10.1364/AO.38.005679.
McFarland, M. J., R. L. Miller, and C. M. U. Neale, 1990: Land surface temperature derived from the SSM/I passive microwave brightness temperatures. IEEE Trans. Geosci. Remote Sens., 28, 839–845, https://doi.org/10.1109/36.58971.
McPherson, R. A., and Coauthors, 2007: Statewide monitoring of the mesoscale environment: A technical update on the Oklahoma Mesonet. J. Atmos. Oceanic Technol., 24, 301–321, https://doi.org/10.1175/JTECH1976.1.
Minnis, P., A. V. Gambheer, and D. R. Doelling, 2004: Azimuthal anisotropy of longwave and infrared window radiances from CERES TRMM and Terra data. J. Geophys. Res., 109, D08202, https://doi.org/10.1029/2003JD004471.
Ohmura, A., and Coauthors, 1998: Baseline Surface Radiation Network (BSRN/WCRP): New precision radiometry for climate research. Bull. Amer. Meteor. Soc., 79, 2115–2136, https://doi.org/10.1175/1520-0477(1998)079<2115:BSRNBW>2.0.CO;2.
Okuyama, A., and Coauthors, 2015: Preliminary validation of Himawari-8/AHI navigation and calibration. Proc. SPIE, 9607, 96072E, https://doi.org/10.1117/12.2188978.
Pinheiro, A. C. T., J. L. Privette, and P. Guillevic, 2006: Modeling the observed angular anisotropy of land surface temperature in a Savanna. IEEE Trans. Geosci. Remote Sens., 44, 1036–1047, https://doi.org/10.1109/TGRS.2005.863827.
Pinker, R. T., D. Sun, M.-P. Hung, C. Li, and J. B. Basara, 2009: Evaluation of satellite estimates of land surface temperature from GOES over the United States. J. Appl. Meteor. Climatol., 48, 167–180, https://doi.org/10.1175/2008JAMC1781.1.
Pinker, R. T, and Coauthors, 2019: Towards a unified and coherent land surface temperature Earth system data record from geostationary satellites. Remote Sens., 11, 1399, https://doi.org/10.3390/rs11121399.
Price, J. C., 1977: Thermal inertia mapping: A new view of the Earth. J. Geophys. Res., 82, 2582–2590, https://doi.org/10.1029/JC082i018p02582.
Ren, H., X. Ye, R. Liu, J. Dong, and Q. Qin, 2018: Improving land surface temperature and emissivity retrieval from the Chinese Gaofen-5 satellite using a hybrid algorithm. IEEE Trans. Geosci. Remote Sens., 56, 1080–1090, https://doi.org/10.1109/TGRS.2017.2758804.
Rossow, W. B., and R. A. Schiffer, 1991: ISCCP cloud data products. Bull. Amer. Meteor. Soc., 72, 2–20, https://doi.org/10.1175/1520-0477(1991)072<0002:ICDP>2.0.CO;2.
Saunders, R. W., M. Matricardi, and P. Brunel, 1999: An improved fast radiative transfer model for assimilation of satellite radiance observations. Quart. J. Roy. Meteor. Soc., 125, 1407–1425, https://doi.org/10.1002/qj.1999.49712555615.
Scarino, B., P. Minnis, R. Palikonda, R. H. Reichle, D. Morstad, C. Yost, B. Shan, and Q. Liu, 2013: Retrieving clear-sky surface skin temperature for numerical weather prediction applications from geostationary satellite data. Remote Sens., 5, 342–366, https://doi.org/10.3390/rs5010342.
Scarino, B., P. Minnis, T. Chee, K. M. Bedka, C. R. Yost, and R. Palikonda, 2017: Global clear-sky surface skin temperature from multiple satellites using a single-channel algorithm with angular anisotropy corrections. Atmos. Meas. Tech., 10, 351–371, https://doi.org/10.5194/amt-10-351-2017.
Schmit, T. J., P. Griffith, M. M. Gunshor, J. M. Daniels, S. J. Goodman, and W. J. Lebair, 2017: A closer look at the ABI on the GOES-R series. Bull. Amer. Meteor. Soc., 98, 681–698, https://doi.org/10.1175/BAMS-D-15-00230.1.
Schmugge, T., S. J. Hook, and C. Coll, 1998: Recovering surface temperature and emissivity from thermal infrared multispectral data. Remote Sens. Environ., 65, 121–131, https://doi.org/10.1016/S0034-4257(98)00023-6.
Seemann, S. W., E. E. Borbas, R. O. Knuteson, G. R. Stephenson, and H.-L. Huang, 2008: Development of a global infrared land surface emissivity database for application to clear sky sounding retrievals from multispectral satellite radiance measurements. J. Appl. Meteor. Climatol., 47, 108–123, https://doi.org/10.1175/2007JAMC1590.1.
Snyder, W., Z. Wan, Y. Zhang, and Y. Feng, 1998: Classification based emissivity for land surface temperature measurement from space. Int. J. Remote Sens., 19, 2753–2774, https://doi.org/10.1080/014311698214497.
Sobrino, J. A., and D. Skoković, 2016: Permanent stations for calibration/validation of thermal sensors over Spain. Data,