1. Introduction
The East Antarctic Ice Sheet (EAIS) is always a region of interest due to its critical role in the climate system of the Earth (Studinger and others, Reference Studinger2005; Seo and others, Reference Seo2016) and its possession of the largest ice volume, equivalent to ~52 m of global sea level (Pattyn and Morlighem, Reference Pattyn and Morlighem2020). In recent decades, the EAIS was considered to be less susceptible to climate warming (Remy and Frezzotti, Reference Remy and Frezzotti2006) than the Arctic (Graversen and others, Reference Graversen, Mauritsen, Tjernström, Källén and Svensson2008), West Antarctica (Steig and others, Reference Steig2009) and Antarctic Peninsula (Vaughan and others, Reference Vaughan2003). However, mass loss has recently been observed in the EAIS coastal regions, including Wilkes Land and Northern Victoria Land (Rignot and others, Reference Rignot2008; Miles and others, Reference Miles, Stokes and Vieli2013; Mengel and others, Reference Mengel and Levermann2014), thereby causing the global sea level rise. The regional mass loss in these regions was associated with Circumpolar Deep Water modified by warming (Miles and others, Reference Miles, Stokes and Vieli2013), suggesting that the Earth's largest ice sheet was likely more sensitive to climate change than previously recognized. Various climate projections indicated that the Antarctic Ice sheet melting could contribute ~1 m to global sea level rise by 2100 (DeConto and Pollard, Reference DeConto and Pollard2016). However, the extent of the EAIS's contribution to future sea level rise under future climate warming remained highly uncertain (Wong and others, Reference Wong, Bakker and Keller2017). In this context, it is imperative to make accurate and extensive temperature measurements, and subsequently to re-evaluate the sensitivity of the EAIS to temperature changes.
The land surface temperature (LST) reflects the energy budget at the land surface, and its long-term variations across the EAIS have the potential to significantly alter the surface mass balance and thus affect the ice-sheet stability (Remy and Frezzotti, Reference Remy and Frezzotti2006; Tuckett and others, Reference Tuckett2019; Baumhoer and others, Reference Baumhoer, Dietz, Kneisel, Paeth and Kuenzer2021). It also impacts changes in atmospheric circulation (Bromwich and Wang, Reference Bromwich and Wang2008; Clem and Fogt, Reference Clem and Fogt2015; Screen and others, Reference Screen, Bracegirdle and Simmonds2018) and polar terrestrial ecosystems (Amesbury and others, Reference Amesbury2017). Although considerable efforts have been made to establish staffed weather stations for meteorological observations on the Antarctic continent since the International Geophysical Year (IGY) of 1957/1958, among them only three (Dome C, South Pole and Syowa Stations also included in the Baseline Surface Radiation Network [BSRN]) can provide continuous longwave radiation measurements for the estimation of surface skin temperature (Driemel and others, Reference Driemel2018). Since the 1980s, more than 260 automatic weather stations (AWSs) have been set up in Antarctica (Allison and others, Reference Allison, Wendler and Radok1993; Turner and others, Reference Turner2004; Lazzara and others, Reference Lazzara, Weidner, Keller, Thom and Cassano2012; Wang and others, Reference Wang2023), some of which alternatively offer skin temperature measurements. However, until recently, in situ surface temperature observations have been extremely sparse in space and time largely because of the harsh environmental conditions.
Benefiting from the short revisiting cycle and extensive spatial coverage, satellites have enabled faster retrieval of skin temperature across the entire EAIS (Comiso, Reference Comiso2000; Retamales and others, Reference Retamales, Durán-Alarcón and Mattar2019). Numerous spaceborne sensors including Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) (Wang and Liang, Reference Wang and Liang2009; Duan and others, Reference Duan, Li, Cheng and Leng2017), Thematic Mapper (TM) (Nguemhe Fils and others, Reference Nguemhe Fils2018), Enhanced Thematic Mapper Plus (ETM + ) (Kumari and others, Reference Kumari2018), Moderate-Resolution Imaging Spectroradiometer (MODIS) (Wan and others, Reference Wan, Zhang, Zhang and Li2002, Reference Wan, Zhang, Zhang and Li2004) and Thermal Infrared Sensor (TIRS) (Tan and others, Reference Tan, Liao, Du and Wu2017), are currently used to generate thermal infrared (IR) LST products. Although ASTER (90 m) and Landsat (TM/ETM + /TIRS, 30 m) products have higher spatial resolution than MODIS, their longer revisit cycle (~16 d) constrains their application in detecting daily LST variations. Owing to high temporal resolution, MODIS-derived LST products, covering the period from 2000 onwards, have been widely applied in many fields (e.g. Xu and Shen, Reference Xu and Shen2013; Kilpatrick and others, Reference Kilpatrick2015; Qie and others, Reference Qie, Wang, Wu and Chen2020), especially in inaccessible polar regions, such as reconstruction of near-surface air temperature (Zhang and others, Reference Zhang2022), investigation of LST patterns (Kang and others, Reference Kang, Im, Lee and Quackenbush2014; Retamales and others, Reference Retamales, Durán-Alarcón and Mattar2019) and determination of surface melting of the ice sheets (Nghiem and others, Reference Nghiem2012; Hall and others, Reference Hall2013).
Various validations or comparisons of the MODIS LST C5 and C6 products have already been conducted in non-polar regions, such as the Savannah landscape (Jacob and others, Reference Jacob2004), city (Ayanlade, Reference Ayanlade2016), cropland (Coll and others, Reference Coll2005, Reference Coll, Wan and Galvem2009), inland water (Wan and others, Reference Wan2008; Crosman and Horel, Reference Crosman and Horel2009), permafrost (Muster and others, Reference Muster, Langer, Abnizova, Young and Boike2015), bare soil surfaces (Duan and others, Reference Duan, Li, Cheng and Leng2017), sand dune and forest (Li and others, Reference Li2014; Duan and others, Reference Duan2018). The offsets of MODIS LST values from in situ measurements were generally within 1°C for most surface types (Wan and others, Reference Wan, Zhang, Zhang and Li2004; Wan, Reference Wan2014). For the cryosphere, the validations of MODIS LST have been performed in the Arctic (Hall and others, Reference Hall2008; Koenig and Hall, Reference Koenig and Hall2010; Westermann and others, Reference Westermann, Langera and Boikea2012; Østby and others, Reference Østby, Schuler and Westermann2014) and alpine mountains (Williamson and others, Reference Williamson2017; Kindstedt and others, Reference Kindstedt2022). However, studies concerning this topic were relatively rare in Antarctica (Wang and others, Reference Wang, Wang and Zhao2013; Liu and others, Reference Liu2015), likely due to the limited availability of in situ skin temperature observations. Most recently, a preliminary assessment of accuracy of MODIS LST C5 products was conducted in East Antarctica (Freville and others, Reference Freville2014). Scambos and others (Reference Scambos2018) also compared MODIS C5 and C6 LST with air temperature data from weather stations during winter in the upper East Antarctic ice divide. Nevertheless, the latest version (C6.1) of the MODIS LST product remains unvalidated in East Antarctica, and the factors influencing its accuracy should be explored.
The main objective of this study was to assess the accuracy of the recently released MODIS C6.1 LST swath product across East Antarctica, using the in situ LSTs calculated by longwave radiation measurements at 12 weather stations. The factors affecting MODIS LST accuracy were also examined to identify potential enhancements for the ongoing algorithm development.
2. Data and methods
2.1. MODIS LST products
The MODIS instruments onboard the Terra and Aqua spacecraft on the NASA Earth Observation System offer important observations of atmospheric water vapor, land surface conditions, cloud properties, aerosols and surface temperature, etc., over the entire Earth's surface every 1–2 d (Justice and others, Reference Justice1998). This study utilized the LST products generated by a generalized split-window LST algorithm using the emissivity in band 31 (≈11 μm) and band 32 (≈12 μm) developed by Wan and Dozier (Reference Wan and Dozier1996). The MODIS LST family has undergone several updates since the release of the first version. Here, we assessed the MOD/MYD11 collection 6.1 (C6.1) Level 2 swath products (here called MxD11_L2) from the Terra and Aqua satellites, which provided per-pixel LST values with a pixel size of 1 km and severed as the basis for the MODIS LST Level 3 products. On average, the available LSTs from the MxD11_L2 product were approximately five times per day at the coastal regions and more than nine times per day on the Antarctic Plateau. This most recent LST product made several calibration enhancements, including modification to the response-versus-scan angle (RVS) approach affecting reflectance bands for MODIS instruments, and adjustments of the optical crosstalk in IR bands of Terra MODIS (Wan and others, Reference Wan, Hook and Hulley2021). We also assessed the gridded MxD11C1 dataset, a Level 3 LST product at a spatial resolution of 0.05 × 0.05° in an equal-angle geographic projection, which had four LST values at four distinct overpass times in a day (twice each at the Aqua and Terra satellites). Quality Control Science Dataset (QC SDS) were included in the two products. We selected only those LST values at each gridcell flagged as ‘good quality’ in the QC layer, which only contained the pixels under clear-sky conditions at a confidence level (defined in the MxD35 cloud mask product) of ≥95% above the land surface.
2.2. In situ measurements
The downwelling and upwelling longwave radiation (LWdown and LWup) records at 12 stations (Fig. 1) were collected to determine snow skin temperatures, and subsequently to validate MODIS LST products. These stations comprised nine AWSs (Jakobs and others, Reference Jakobs, Carleen, Michiel, Smeets and Gert2020) and three BSRN stations (Driemel and others, Reference Driemel2018), each with records at varying timespans (Table 1). AWSs and BSRN stations provided hourly and minute-by-minute radiation records, respectively. Of these, two BSRN stations (DOM and SPO) and two AWSs (AWS12 and AWS13) were situated on the East Antarctic Plateau, whereas three AWS stations (AWS5, AWS10 and AWS11) and the SYO BSRN station were distributed along the coast of the EAIS. AWS6, AWS8, AWS9 and AWS16 were located at the edge of the plateau (Fig. 1). If LWup is lower than LWdown of +5 W m−2, the corresponding radiation values were excluded to ensure the accuracy of the observations (Freville and others, Reference Freville2014). According to the Stefan–Boltzmann law, the snow skin temperature (T s) was calculated by LWdown and LWup at the stations via the following formula.
where T s represented the snow skin temperature (hereafter, in situ T s), σ was the Stefan–Boltzmann constant (5.67 × 10−8 W m−2 K−4) and ɛ denoted the snow surface emissivity. The surface emissivity of snow was estimated to be between 0.985 and 0.990 for all grain sizes (Dozier and Warren, Reference Dozier and Warren1982). Sensitivity tests by Freville and others (Reference Freville2014) showed that setting the snow surface emissivity to 1.0 and 0.98 in comparison to 0.99 resulted in an average surface temperature difference of only ~0.08°C (0.080 and 0.079°C, respectively). Thus, we used a constant snow surface emissivity value of 0.99 to calculate the in situ T s (Brun and others, Reference Brun2011).
Map of the study area and the weather station locations.
Summary of the station information
All stations were divided into three groups by their elevations. The first four stations (AWS5, AWS10, AWS11 and SYO) were classified as the coastal group (<1000 m) and the last four (AWS12, AWS13, DOM and SPO) as the plateau group (>3000 m). The middle four (AWS6, AWS9, AWS8 and AWS16) were classified as the group located on the edge of the plateau (1000–3000 m)
At all 12 stations, we recorded in situ T s within 5 min of the view time (hereafter, satellite overpass time) provided by the MODIS LST products. The mean bias (MB), root mean square error (RMSE) and determination coefficient (R 2) were used to quantify the uncertainties of the MODIS LST products.
where n, y i and t i were the number of in situ observations, the MODIS LST values and in situ T s, respectively. $\bar{y}$
and $\bar{t}$
were the average of satellite-derived and in situ surface temperatures, respectively.
To more effectively remove cloud contamination from MODIS LST pixels, and to identify the environmental factors affecting the accuracy of satellite-derived LSTs, the cloud fraction, 2 m air temperature (T 2m), solar radiation, wind speed and relative humidity observations from the nine AWSs were used. The cloud fraction is the physically meaningful cloud cover obtained from the longwave radiation balance measured by the AWS and was calculated based on the difference in emissivity between clear and cloudy atmospheres. Kuipers Munneke and others (Reference Kuipers Munneke, Reijmer and van den Broeke2011) plotted the hourly LWdown against the T 2m data in a scatterplot. The upper boundary of the scatterplot coincided with ${\rm L}{\rm W}_{{\rm down}} = \sigma T_{2{\rm m}}^4$
(atmospheric emissivity ɛ ≈ 1, completely cloudy-sky emitting as a blackbody radiator); the lower boundary of the scatterplot could be approximately fitted by a second-order polynomial to represent clear-sky conditions (ɛ ≈ 0). The cloud fraction (between 0 and 1) was produced using linear interpolation between the two boundaries (details in Kuipers Munneke and others, Reference Kuipers Munneke, Reijmer and van den Broeke2011). Using these cloud fraction data, we set up a filter to eliminate cloud-contaminated MODIS LSTs. MODIS LST grid cells with the corresponding AWS cloud fraction exceeding the threshold of 0.3 were classified as cloud-contaminated, and subsequently screened out. After removing the cloud-contaminated values, Pearson's correlation coefficient (ryt) was employed to quantify the relationships between meteorological factors and the biases in the two sets of skin temperature data.
3. Results
3.1. Accuracy of the MODIS LST products
3.1.1. Accuracy of the original MxD11_L2 and MxD11C1 LST products at satellite overpass times
At all the stations, the biases between the MxD11_L2 LST and in situ T s predominantly fell below 0°C (Figs 2, 3), with the highest frequency between −4 and 0°C, suggesting that the MODIS LST was generally lower than in situ T s in most of the cases. Specially, the MODIS LST exhibited mean cold biases ranging from −4.86 to −0.66°C relative to the near-synchronous in situ T s, with the exception of a minor warm bias (0.11°C) at the Syowa station. The RMSE values ranged from 3.7 to 7.32°C, and R 2 values exceeded 0.75 at the majority of stations. Spatially, cold biases were less pronounced at stations on the East Antarctic Plateau (−1.39 to −0.66°C) and more significant at those along the edge of plateau (−4.86 to −1.68°C), with coastal stations showing intermediate values (−3 to −1.76°C). In contrast, the R 2 values were relatively low for the stations at the edge of plateau (0.59–0.84), higher for those on plateau (0.76–0.91) and for coastal stations (0.64–0.83) still in between. However, these spatial patterns were not observed in the RMSEs.
Histograms of the biases between the MODIS LST and in situ T s at (a) all stations, (b) Antarctic coast, (c) edge of the plateau and (d) Antarctic Plateau. The red vertical line denoted the bias value equal to 0.
comparison between the MODIS LST and in situ T s at satellite overpass times at stations located on the Antarctic coast (top group, a–d), the edge of the plateau (middle group, e–h) and the East Antarctic Plateau (bottom group, i–l). The red line represented the one-to-one line. The green and blue points represented the samples in warm period (October to March of the following year) and cold period (April to September), respectively. N stood for the paired amount of two types of LST observations.
Similarly, compared with the in situ T s, the MxD11C1 LST exhibited mean cold biases (−5.27 to −0.48°C), except for a slight warm bias (0.17°C) at the Syowa station, with RMSE values ranging from 3.17 to 6.46°C and R 2 exceeding 0.7 at all stations (Fig. S1). The spatial pattern of the MxD11C1 LST biases was also the same as that of the L2 product, indicating smaller biases at stations on the East Antarctic Plateau (−1.07 to −0.48°C) and larger ones along the edge of plateau (−5.27 to −4.49°C), with the intermediate values at coastal stations (−1.84 to 0.17°C). The R 2 values were relatively low for the stations at the edge of plateau (0.7–0.73), higher for those on the plateau (0.87–0.94) and intermediate for coastal stations (0.76–0.89). These results basically coincided with those of the MxD11_L2 product in terms of magnitude.
A large numerous of sample points were observed well below the 1:1 line at almost all stations (Fig. 3), suggesting substantial cold biases in MODIS LST in comparison to in situ T s. In particular, some bias values dipped below −20°C (Fig. 2). Spatially, these large cold biases were more prevalent at stations near the coast, while they were comparatively infrequent on the Antarctic Plateau. A possible reason for these large biases was the imperfect cloud mask of MODIS LST products, as noted in earlier studies on previous MODIS LST versions in Greenland (Hall and others, Reference Hall, Key, Casey, Riggs and Cavalieri2004; Westermann and others, Reference Westermann, Langera and Boikea2012; Adolph and others, Reference Adolph, Albert and Hall2018). Since the temperature of clouds tend to be much lower than that of snow surface, misidentification of cloud and snow surfaces can limit the accuracy of surface temperature products in ice- and snow-covered regions. This could be further confirmed by the significant correlations between the MODIS LST biases and AWS cloud fraction (Table 2). Seasonally, this phenomenon occurred more frequently during warm periods than during cold periods (Fig. 3 and Table 2). Therefore, we initially implemented a more restrictive cloud mask based on AWS cloud fraction to further reduce the effects of cloud on the MODIS LST. Due to the absence of cloud records at the three BSRN stations, it was not possible to reduce cloud contamination using a unified filter, and thus these data were not included in the following analyses.
Correlation coefficients (r) between AWS cloud fraction and the biases of original MODIS LST at different AWSs
Bolded r values were statistically significant at p < 0.05 level.
3.1.2. Accuracy of the MxD11_L2 LST through filtering by in situ cloud data
As shown in Figure 4, the stricter cloud mask based on the AWS cloud fraction greatly reduced the number of cloud-affected LST measurements, and the filtered MODIS LST could be considered to be cloud-free. Even after filtering out the cloud effects by means of AWS cloud data, the MODIS LST still exhibited cold biases (−5.18 to −0.33°C) in comparison to in situ T s at all AWSs (Fig. 4). The RMSE values varied between 2.37 and 6.28°C, and R 2 values ranged from 0.72 to 0.95. Clearly, by reducing the impacts of cloud, accuracy of the MODIS LST was notably improved at all stations except AWS9. Specially, the MODIS LST biases decreased by ~1°C at these stations. The RMSE values reduced by 0.16–2.26°C, with an average of 1.39°C, and R 2 rose by 0.04–0.18 with an average of 0.09. Spatially, the improvement in MODIS LST accuracy was most pronounced at the stations located along the Antarctic coast and the edge of plateau. The RMSEs reduced by more than 1°C and R 2 values increased by more than 0.1 at almost all stations. This indicated that cloud cover significantly impacted the MODIS LST accuracy in these regions. However, at AWS9, the improvement in LST accuracy was minimal, with a slight increase in RMSE and bias by 0.09 and 0.32°C, respectively. This can be explained by no significant correlation between MODIS LST biases and cloud fraction (Table 2), suggesting that factors other than cloud conditions primarily affect the MODIS LST accuracy here.
Comparison between in situ T s and the MODIS LST after removing cloud-affected data at the satellite overpass times at stations located on the Antarctic coast (a–c), the edge of the plateau (d–g) and the East Antarctic Plateau (h–i). The red line represented the one-to-one line. The green and blue points represented the samples during warm period (October to March of the following year) and during cold period (April to September), respectively. N was the paired amount of two types of LST observations.
We calculated the daily averages of MODIS LST and in situ T s respectively to examine the performance of the satellite-derived LST on daily scale. Relative to the corresponding in situ T s (Fig. 5), the daily-averaged MODIS LST still presented cold biases at all stations (−5 to −0.39°C). The RMSE values ranged from 1.47 to 5.65°C, which were lower than those at satellite overpass times in magnitude. The R 2 at these stations also increased, with all above 0.73 and most above 0.9. Spatially, MODIS LST accuracy improvement was most prominent at stations on the Antarctic coast and on the plateau, with the R 2 values of above 0.87, the RMSE values of below 3.1°C and cold biases of smaller than −1.52°C. For the stations located at the edge of the plateau, the LST accuracy was relatively higher in AWS6 and AWS8, while there was still more room for improvement in AWS9 and AWS16.
Comparison between in situ T s and the MODIS LST after removing cloud-affected data at the daily-scale at the stations located on the Antarctic coast (a–c), the edge of the plateau (d–g) and the East Antarctic Plateau (h–i). The red line represented the one-to-one line. The green and blue points represented the samples during warm period (October to March of the following year) and during cold period (April to September), respectively. N was the paired amount of two types of LST observations.
3.1.3. Cloud-screened MODIS LST accuracy for the monthly periods
Aggregating the data into daily periods led to closer agreement between the MODIS LST and in situ T s (see the description above and Figs 4, 5). To determine if the monthly synthesis can result in more robust alignment, the mean monthly biases of the two datasets were calculated. As shown in Figure 6, the difference between the two datasets was not eliminated by the monthly synthesis. We found that the biases showed a distinct seasonality at stations along the coast and the edge of plateau, with larger biases during the cold period (April to September), peaking during May–July, and smaller biases during the warm period (October to March). Both the means and medians of biases were below 0°C for each month with the mean values obviously lower than the medians. This pattern was also evident across all stations (Fig. 6a). Unfortunately, the relatively low pair counts (no more than 25) between the two datasets limited their statistical significance during the cold period at plateau stations, leading to the omission of the results for those months (Fig. 6d). During the warm period, the monthly biases were maintained at a low level at these stations, indicating the high accuracy of the MODIS LST.
Temperature difference between cloud-screened MODIS LST and in situ T s in each month at all stations (a), the stations on the Antarctic coast (b), the edge of the plateau (c) and the East Antarctic Plateau (d). Boxplots showed the 25th and 75th percentiles (boxes), the averages (green triangles) and the medians (orange short horizontal lines). The numbers at the top of the boxes were the paired amount of two types of LST observations in each month, and gray boxes in subplot (d) represented the months with fewer than 25 paired amounts. The red horizontal line represented the bias value equal to 0°C. The results cannot be extrapolated to other places without an AWS (equipped with cloud fraction).
3.2. Factors influencing the MODIS LST accuracy in East Antarctica
Cloud fraction, near-surface temperature inversion intensity (T air–T surface) and original MODIS LST were potential causes of the discrepancies between the MODIS LST and in situ T s. As previously mentioned, cloud contaminations significantly affect the MxD11_L2 LST accuracy over EAIS. Given that the MxD11_L2 LST served as the input for the generation of MxD11C1, and that no additional cloud-reducing algorithms were applied, cloud issues still occurred in MxD11C1. This was further confirmed by the similar accuracies of MxD11_L2 LST and MxD11C1. Near-surface inversions are widespread on the EAIS (Hudson and Brandt, Reference Hudson and Brandt2005; Scambos and others, Reference Scambos, Haran and Massom2006, Reference Scambos2018; Pietroni and others, Reference Pietroni, Argentini and Petenko2014) and typically occur under clear-sky conditions, low incoming solar irradiance and low wind speeds (Table S2 and Fig. S2). In this study, the difference between T air and T surface was calculated to estimate the near-surface inversion intensity. The value (T air–T surface) is positive when the 2 m air temperature is higher than the snow surface temperature, indicating the occurrence of near-surface temperature inversion. Cloud-screened MODIS LST biases correlated significantly with near-surface temperature inversions at all stations (−0.34 to 0.69, p < 0.05) during satellite overpass times (Fig. 7a). By comparison, the correlations at plateau stations were higher than those at the coastal stations, peaking in summer or fall (0.6–0.88, p < 0.05). Conversely, correlations at coastal stations did not exhibit clear seasonality. For stations at the edge of plateau, the correlation pattern was similar to that of two plateau stations (AWS8 and AWS9), possibly due to their higher elevations and greater distance from the ocean. In particular, the MODIS LST biases at AWS9 appeared to be associated with temperature inversions rather than cloud variations (Table 2). In contrast, AWS6 resembled coastal stations because of its lower elevation and closer proximity to the coast.
Correlations between biases and (a) near-surface inversions (T air-T surface), and (b) MODIS LST at different time scales, respectively. From left to right, the values in each column represented the correlations at satellite overpass time (Sot), daily, austral summer (Sum), fall (Fal), winter (Win) and spring (Spr), respectively. The correlations labeled ‘*’ were significant at p < 0.05 level with the two-sided test. The winter at AWS 9 and 12 was filled in white due to the lack of data.
Significant correlations between biases and MODIS LST measurements were observed at the stations near the coast (AWS6, AWS16 and coastal stations; Fig. 7b). Seasonally, these correlations were highest in austral winter, a period characterized by more frequent occurrence of extreme low temperatures in MODIS LSTs. This difference was most pronounced at AWS16, where all correlation coefficients were above 0.52 (p < 0.05) at all scales. This means that when the original MODIS LST values were lower in the coastal region, their biases tended to be larger.
4. Discussion
To our knowledge, this study provided the first assessment of the accuracy of the latest MODIS C6.1 product over East Antarctica using available surface radiation-derived in situ T s. Our results suggested that this MODIS LST product had a good performance in fitting in situ T s, which was relatively accurate in aiding our understanding LST spatiotemporal patterns and ice–atmosphere interactions across the EAIS. However, cloud contamination emerged as a critical factor affecting MODIS LST accuracy, which was consistent with the findings from several previous studies on polar regions (Hall and others, Reference Hall2008; Østby and others, Reference Østby, Schuler and Westermann2014; Sun and others, Reference Sun2017; Zhang and others, Reference Zhang, Tan and Shi2018). In this study, we chose to use clear-sky MODIS LST values flagged as ‘good quality’. However, significant correlations (−0.56 to −0.39, p < 0.05) were observed between AWS cloud fraction and MODIS LST biases at all stations except AWS9, especially during the warm period characterized by higher cloud fraction (Fig. 3 and Table 2). Many cloud-contaminated data could also be found in the scatterplot (Fig. 3). This indicated that the current cloud mask algorithm of MODIS LST product was imperfect, and only partially effective over the EAIS. The pixels contaminated by clouds were not fully identified and excluded, which likely contributed much to the large cold biases (Fig. 3). Despite losing 29.4% of the data, the application of AWS cloud fraction as a stricter filter to remove cloud contamination resulted in the significant decreases in the biases and RMSEs between the two LST datasets, along with an increase in R 2. The improvement was more obvious at the stations located along the coast and the edge of plateau (Figs 3, 4), suggesting that cloud contamination occurred more frequently in places near the coast where air was relatively moist, and less on inland plateau where the sky was clearer. The MxD11_L2 LST was the original swath product of MODIS LST family and provided the input for the generation of Level 3 LST products. If a deficient cloud mask was employed, the accuracy of Level 3 product will likely be compromised. We verified this phenomenon using the MxD11C1 product at several stations with longer time series (Fig. S1). At these stations, the MxD11C1 LST exhibited cold biases (−5.27 to −0.48°C) with respective to the in situ T s, except for a slight warm bias (0.17°C) at the Syowa station. The spatial pattern of the MxD11C1 LST biases was also the same as that of the L2 product. These bias values were basically coincided with the magnitude of the original MxD11_L2 product (without using AWS data to remove cloud effects), and the impact of cloud on the MODIS LST was largely retained in the Level 3 product (Fig. 3 and Fig. S1).
Cloud-screened MODIS LSTs still exhibited significant cold biases relative to in situ T s at different timescales, which was in accordance with previous validations focusing on ice-covered polar regions (Scambos and others, Reference Scambos, Haran and Massom2006, Reference Scambos2018; Koenig and Hall, Reference Koenig and Hall2010; Westermann and others, Reference Westermann, Langera and Boikea2012; Østby and others, Reference Østby, Schuler and Westermann2014; Zikan and others, Reference Zikan, Adolph, Brown and Fausto2023). The bias and RMSE values were comparable in magnitude to those found in earlier estimation of the previous version of the MODIS LST at these sites. The MODIS Version C5 LST in East Antarctica showed cold biases (−2.66 to 0.08°C) in 2009, with worse accuracy at coastal stations (Freville and others, Reference Freville2014). This seemed to show that there was no distinct improvement of C6 MODIS LST over C5 version across the EAIS. We further explicitly presented the spatial differences in biases, with the smaller cold biases occurring at the stations on the Antarctic plateau, but the larger biases occurring at those along the edge of plateau, and biases at coastal stations roughly in between. Accurate calibration of these deviations was essential to faithfully represent meteorological and glaciological conditions.
Near-surface temperature inversions also strongly influenced the MODIS LST accuracy. Previous studies showed that larger cold biases of MODIS LST compared to 2 m air temperatures occur under stronger inversions in polar regions (Scambos and others, Reference Scambos, Haran and Massom2006, Reference Scambos2018; Adolph and others, Reference Adolph, Albert and Hall2018; Kindstedt and others, Reference Kindstedt2022). When comparing the MODIS LST with 2 m air temperatures in a test focusing on the EAIS, similar conclusions were drawn. When strong temperature inversions occurred, the in situ 2 m air temperature was distinctly higher than that at the surface. In this case, since MODIS measured the surface temperature, a greater cold difference is inevitable between the MODIS LST and in situ 2 m temperature. When we compared in situ skin temperatures with MODIS LST, high positive correlations demonstrated that inversions still exert a strong influence on the accuracy of MODIS LST, especially on the East Antarctic Plateau, where temperature inversions are more frequent and intense (Connolley, Reference Connolley1996; Scambos and others, Reference Scambos2018). This highly positive correlation suggested that MODIS LST tended to exhibit either warm biases or small cold biases compared to in situ surface temperatures when strong temperature inversions occurred. Some evidence of this can be found in Figure 4 (after removing cloud effects), with some points located far above the 1:1 line, representing the warm biases of MODIS LST in comparison to in situ surface temperatures. Yu and others (Reference Yu, Rothrock and Lindsay1995) showed that ice-crystal precipitation (or clear-sky precipitation) occurring during temperature inversions may cause warm biases of satellite-derived LST compared to ground-based measurements (below the ice-crystal precipitation layer). The East Antarctic Plateau was subjected to frequent temperature inversions and ice-crystal precipitation (Kikuchi and Hogan, Reference Kikuchi and Hogan1979; Walden and others, Reference Walden, Warren and Tuttle2003; Vignon and others, Reference Vignon, Roussel, Gorodetskaya, Genthon and Berne2021), which to some extent explained the occurrence of warm biases and positive correlations between temperature inversions and MODIS LST biases at the Antarctic Plateau stations. Therefore, the influence of atmospheric phenomena occurring simultaneously with inversion and high humidity on MODIS LST accuracy should be examined in the future.
In the coastal region, a ‘progressive’ cold bias was found that was the bias between two LST datasets gradually increased as the MODIS LST decreased. This was most obvious for the AWS16 station and could also be seen in the scatterplot (Fig. 4e). The MODIS LST exhibited greater cold biases at lower temperatures on the coastal area, indicating that MODIS LST needs to be corrected for low temperatures in the future (Xiong and others, Reference Xiong2015; Zikan and others, Reference Zikan, Adolph, Brown and Fausto2023).
We further explored the potential relationships between MODIS LST biases and other environmental factors. Several previous studies reported that biases of C5 and C6 MODIS LST products were associated with the scan angle of the Terra and Aqua satellites (Adolph and others, Reference Adolph, Albert and Hall2018). However, our study demonstrated that there was no correlation between MODIS LST biases and MODIS view angle in the latest version (C6.1) (Table S1). Furthermore, the latest C6.1 product has undergone improvements in the RVS approach, which affected reflectance bands for Aqua and Terra MODIS (Wan and others, Reference Wan, Hook and Hulley2021). These findings suggested that the impact of the scan angle of satellites have greatly decreased, and likely have been even eliminated. There was no relationship between relative humidity and MODIS LST biases at almost all stations (Table S1), which was inconsistent with the results of Zikan and others (Reference Zikan, Adolph, Brown and Fausto2023) that station-level atmospheric humidity may affect the biases of MODIS LST (C6). Zikan and others (Reference Zikan, Adolph, Brown and Fausto2023) proposed that the presence of hoarfrost growing on sensors under high humidity could cause inaccurate temperature measurements and contribute to errors, but this issue was not found in this study. Additionally, the station-level wind speeds are almost uncorrelated with MODIS LST biases (Table S1). However, a weak negative correlation (r = −0.22, p < 0.05) was observed between solar radiation and biases at stations on the Antarctic Plateau. We assumed that these relationships were likely linked through temperature inversions, which usually occur at low solar radiation (Table S2 and Fig. S2). Previous results also demonstrated strong and widespread temperature inversions in the interior of Antarctica, but relatively weak in coastal regions (Hudson and Brandt, Reference Hudson and Brandt2005; Pietroni and others, Reference Pietroni, Argentini and Petenko2014; Scambos and others, Reference Scambos2018).
Temperature inversions and relative humidity were analyzed on only a narrow vertical scale (at the station level) in this study, and hence it is necessary to explore the temperature inversion intensity and variation of atmospheric water vapor in the thicker atmospheric troposphere, and to investigate its impacts on satellite-derived products in the future. Moreover, the harsh environmental conditions over EAIS such as extreme weather (winds and extremely low temperatures) usually result in intermittent gaps in AWS meteorological observations (Zhang and others, Reference Zhang2022). This somewhat limited our study of factors affecting MODIS LST accuracy, especially in winter (Fig. 7).
5. Conclusions
MODIS has proved to be a valuable tool for monitoring surface temperatures across East Antarctica. The latest MODIS LST C6.1 product highly and significantly correlated with in situ T s calculated by radiation from 12 stations in East Antarctica, with the determination coefficient (R 2) of >0.75 at most stations. However, there was no distinct improvement in the C6.1 MODIS LST over C5 version on the EAIS. Significant cold biases occurred in MODIS LST at different timescales (transient, daily and monthly), and these biases were more pronounced during cold period than during warm period. Spatially, the biases were smaller at inland stations, but larger for the stations on coastal regions and the edge of plateau. The cloud fraction, low temperature and near-surface temperature inversion intensity contributed much to the MODIS LST biases. Cloud contamination was effectively removed by utilizing a stricter filter based on AWS cloud fraction, which thereby significantly improved MODIS LST accuracy. Therefore, there was a need to develop a more stringent cloud mask algorithm to further reduce cloud effects on MODIS LST product, especially at the Antarctic coastal regions. The application of AWS observations was a good choice for improving MODIS cloud mask algorithm. Low temperature corrections were also required in coastal regions. In addition, on the Antarctic Plateau, ice-crystal precipitation occurring during temperature inversions significantly impacted the MODIS LST biases. The energy balance model appears to be a promising approach for calibrating MODIS LST in the future research. Notably, validating satellite-generated LST observations in East Antarctica remains challenging due to the sparseness of in situ observations. Therefore, strengthening in situ LST monitoring particularly during winter is highly needed.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/jog.2024.17.
Data
The MODIS LST products used in this study are available at http://modis.gsfc.nasa.gov. The AWSs datasets are available at https://doi.org/10.1594/PANGAEA.910473 and the BSRN stations datasets are available at https://dataportals.pangaea.de/bsrn/3. The unpublished data used in this paper can be obtained by contacting the author if they are needed.
Acknowledgements
This work was financially supported by the National Key Research and Development Program of China (2020YFA0608202), the National Natural Science Foundation of China (41971081 and 41830644), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA19070103) and the Project for Outstanding Youth Innovation Team in the Universities of Shandong Province (2019KJH011). The authors are grateful for the comments from Ghislain Picard (reviewer) and the anonymous reviewer which have improved this manuscript.
Competing interest
None.
Open access
