Next Article in Journal
Distinguishing Photosynthetic and Non-Photosynthetic Vegetation: How Do Traditional Observations and Spectral Classification Compare?
Next Article in Special Issue
Comparison of Harmonic Analysis of Time Series (HANTS) and Multi-Singular Spectrum Analysis (M-SSA) in Reconstruction of Long-Gap Missing Data in NDVI Time Series
Previous Article in Journal
Precise Orbit Determination for BeiDou GEO/IGSO Satellites during Orbit Maneuvering with Pseudo-Stochastic Pulses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Climatologies of Average Monthly Air Temperature over Mongolia Using MODIS Land Surface Temperature (LST) Time Series and Machine Learning Techniques

by
Munkhdulam Otgonbayar
1,
Clement Atzberger
2,*,
Matteo Mattiuzzi
3 and
Avirmed Erdenedalai
1
1
Institute of Geography and Geoecology, Mongolian Academy of Sciences, Ulaanbaatar 15170, Mongolia
2
Institute of Surveying, Remote Sensing and Land Information (IVFL), University of Natural Resources and Life Sciences (BOKU), 1190 Vienna, Austria
3
European Environmental Agency, 1050 Copenhagen, Denmark
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(21), 2588; https://doi.org/10.3390/rs11212588
Submission received: 2 September 2019 / Revised: 28 October 2019 / Accepted: 30 October 2019 / Published: 4 November 2019

Abstract

:
The objective of this research was to develop a robust statistical model to estimate climatologies (2002–2017) of monthly average near-surface air temperature (Ta) over Mongolia using Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) time series products and terrain parameters. Two regression models were analyzed in this study linking automatic weather station data (Ta) with Earth observation (EO) images: Partial least squares (PLS) and random forest (RF). Both models were trained to predict Ta climatologies for each of the twelve months, using up to 17 variables as predictors. The models were applied to the entire land surface of Mongolia, the eighteenth largest but most sparsely populated country in the world. Twelve of the predictor variables were derived from the LST time series products of the Terra MODIS satellite. The LST MOD11A2 (collection 6) products provided thermal information at a spatial resolution of 1 km and with 8-day temporal resolution from 2002 to 2017. Three terrain variables, namely, elevation, slope, and aspect, were extracted using a Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM), and two variables describing the geographical location of weather stations were extracted from vector data. For training, a total of 8544 meteorological data points from 63 automatic weather stations were used covering the same period as MODIS LST products. The PLS regression resulted in a coefficient of determination (R2) between 0.74 and 0.87 and a root-mean-square error (RMSE) from 1.20 °C to 2.19 °C between measured and estimated monthly Ta. The non-linear RF regression yielded even more accurate results with R2 in the range from 0.82 to 0.95 and RMSE from 0.84 °C to 1.93 °C. Using RF, the two best modeled months were July and August and the two worst months were January and February. The four most predictive variables were day/nighttime LST, elevation, and latitude. Using the developed RF models, spatial maps of the monthly average Ta at a spatial resolution of 1 km were generated for Mongolia (~1566 × 106 km2). This spatial dataset might be useful for various environmental applications. The method is transparent and relatively easy to implement.

1. Introduction

Near-surface air temperature (Ta) is a key descriptor of the climate [1]. Ta is a critical variable to the effective understanding of the many physical and biological processes between the atmosphere and land systems [2,3,4] because it regulates many land surface processes such as photosynthesis, respiration, and evaporation [5]. As air temperature influences nearly all biotic processes [6], climatologies of Ta also permit a good characterization of terrestrial environmental conditions [5,7]. As this variable can change quickly over space, cost-efficient mapping procedures are needed that can depict Ta using high spatial resolution.
Since the early 1980s, various interpolation methods have been used to estimate Ta given adequate sample points [8,9]. The literature shows that the most common interpolation techniques are global interpolators, thin plate smoothing splines and different forms of kriging [10,11], inverse distance weighting [12], and climatologically-aided interpolation [13]. In the comparative study of [14], most interpolation methods gave similar results. However, interpolation errors typically range between 1 and 3 °C [15,16] depending on the spatial and temporal resolution of recorded Ta data and the density of the station network [17].
Direct measurements of Ta at a height of 2 m above ground are only available from a limited number of meteorological stations. In many cases such as Mongolia, the spatial coverage of these measurements is inadequate; in addition, typical Ta time series come with many missing values [6]. On the contrary, satellite-derived land surface temperature (LST) data are continuous in both spatial–temporal coverages and are relatively inexpensive. However, the satellite does not directly measure Ta but only the LST.
Based on the physical linkage between LST and Ta, several authors have offered methods to estimate Ta using remote sensing satellite data [5,18,19]. During past decades, a large body of research has been collected regarding the retrieval of LST from satellite-based thermal infrared (TIR) data [20], in particular that related to the better understanding of emissivity and atmospheric effects [18,21,22]. As a result, LST can be retrieved nowadays relatively accurately from remotely sensed TIR data [23]. Several studies have demonstrated that Ta and LST data are highly correlated [24]. However, as expected, large differences have been noticed [25] which are related, for example, to physical properties and atmospheric conditions [26].
Three major approaches have been used to estimate Ta from LST data [27]:
(1)
energy-balance parameterization based on thermodynamic approaches [15,19,28],
(2)
contextual approaches based on temperature–vegetation index relations (TVX) [5,29,30], and
(3)
statistical approaches using various form of regression techniques [15,16,23,31,32,33,34].
Good exemplary studies estimating air temperatures with MODIS LST products using the aforementioned methods can be found for example in Bartkowiak et al. [35], Lu et al. [36], Zhou et al. [37], Janatian et al. [38], Ho et al. [39], Duan et al. [40], and Benali et al. [3].
Within the last two decades, statistical approaches, including simple and advanced regression, (e.g., linear and multiple regression, and machine learning techniques) have been developed to estimate Ta from Moderate Resolution Imaging Spectroradiometer (MODIS) LST products with varying levels of success. More recently, several studies have investigated more complex and advanced approaches to estimate Ta from MODIS LST products, such as geographically weighted regression (GWR) and climate space weighted regression (CSWR) [6], spatiotemporal regression-kriging (STRK) [31], stepwise [38,41], random forest (RF) [39,41,42,43,44,45], generalized boosted model (GBM) [45], cubist [41,45], support vector machine (SVM) [39], ordinary least squares (OLS) [39] and M5 model tree [46]. To ensure high modeling accuracy, several papers have highlighted the usefulness of multivariate and non-parametric algorithms such as RF and STRK. For instance, Kilibarda et al. [31] estimated mean, maximum, and minimum daily Ta with a spatial resolution of 1 km at a global scale using STRK with MODIS 8-day time-series LST products along with elevation, wetness index, and geographical location. The performance of STRK to predict Ta from MODIS LST products was compared with the performance of the linear regression model. The results indicated that the root-mean-square errors (RMSEs) for predicting mean, maximum, and minimum daily Ta are ±2 °C for areas with a high density of stations and from ±2 °C to ±4 °C for areas with a coarse station density. The lowest accuracy was 6 °C in Antarctica and at locations with high altitudes. Yoo et al. (2018) [43] estimated maximum and minimum daily Ta in two megacities using LST data from MODIS Terra/Aqua and seven auxiliary variables based on the RF machine learning method, resulting in an RMSE of 1.1 °C and 1.2 °C for maximum and minimum Ta, respectively, in Seoul, and an RMSE of 1.7 °C and 1.2 °C for maximum and minimum Ta, respectively, in Los Angeles. Several authors have concluded that machine learning techniques perform better than more conventional methods which provide multi-variables and nonlinear and nonparametric regression and classification [38,39,43,45,47,48]. Machine learning algorithms are particularly useful for cases where no deterministic model is available to solve the problem. Our research objective was to develop a robust empirical model to estimate climatologies of average monthly Ta across Mongolia at 1 km spatial resolution using time-series of MODIS Terra LST products, terrain parameters (elevation, slope, and aspect), and other ancillary information.

2. Study Area

The area studied in this work covers the entirety of Mongolia with a total area of approximately 1566 × 106 km2. Mongolia is the eighteenth largest and most sparsely populated country in the world. Mongolia extends between the latitudes 41°35′ N–52°09′ N and the longitudes 87°44′ E–119°56′ E, with an average land surface elevation of 1580 m above sea level [49]. The elevation ranges between 524 m to 4320 m above sea level (Figure 1a). Its continentality increases from east to west. The climate conditions are extreme continental with semiarid and arid regions (Figure 1b). The “blue sky” country counts on average 260 sunny days per year and is characterized by a long cold winter and a short dry-hot summer with generally low precipitation [50]. The country’s average annual temperature varies between −8 °C and 6 °C [49] with strong temperature gradients. Annual total precipitation ranges between 50 mm in the desert steppe and desert regions to 500 mm in the high mountain regions. Generally, precipitation gradually increases from south to north.
According to the Köppen climate classification, Mongolia can be divided into seven climatic regions (Figure 1b). The most prevalent climates are BWk, BSk, and Dwc. By contrast, Dfc, Dsc, Dwb, and Dsb are less prevalent. Exemplary temperature and precipitation charts for the three most prevalent climates are shown in Figure 2, demonstrating a large variation in precipitation and temperature variation.

3. Data and Methods

3.1. Remote Sensing Data

MODIS LST products are distributed by the Land Processes Distributed Active Archive Center (LP DAAC) in a hierarchical data format or HDF file. We used observations from MOD11 from the Terra satellite. MODIS generates two daily observations which involve one for daytime (LSTd) and one for nighttime (LSTn) at approximately 10:30 and 22:30 local time, respectively. The new collection 6 (c006) of MODIS LST products have been used to estimate Ta [53]. This dataset was made available in 2016. It covers the entire period (2002–2017) and data are of higher quality compared to the earlier collection(s), which had been used for previous studies such as [3,17,46,47,54,55]. The LST accuracy of the c006 products is reported as being approximately twice as good as collection 5 (c005) due to the incorporation of the emissivity adjustment model in the MODIS split-window algorithm [53]. For instance, the c006 LST product reduced the RMSE of bare soil sites of the c005 LST product by 1.24 °C during the day and 0.58 °C at night ([53], p. 88).
To cover the land surface of Mongolia, seven tiles of granules with horizontal (h) and vertical (v) title numbers h23v03, h23v04, h24v03, h24v04, h25v03, h25v04, and h26v04 had been used. The MODIS MOD11A2 c006 data were obtained through the online Data Pool at the National Aeronautics and Space Administration (NASA), the LP DAAC, and the United States Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota. The retrieved MODIS LST used the generalized split-window algorithm [56] to derive surface temperature from the recorded at-satellite radiances.
For this study, we used the MODIS Terra 8-day LST product (MOD11A2) at a spatial resolution of 1 km, gridded in the Sinusoidal projection intervals, and covering the period 2002–2017. The HDF file for the MOD11A2 product includes 12 different scientific data sets (SDSs), as shown in Table 1. A detailed description of SDSs is given in [57,58]. This MOD11A2 product includes daytime and nighttime LST data (LSTd and LSTn), quality information (QCd and QCn), observation information (DvA, NvA, DvT, and NvT), emissivity data (Em31 and Em32), and clear sky coverage (CsD and CsN). The HDF file for this product also contains associated quality science dataset layers which provide users with information regarding the usability and usefulness of the data products. The MODIS LST quality science dataset layers are binary encoded and bit packed. The quality assurance (QA) layer containing integer values had been converted to a bit binary value for interpretation [59,60] (p. 19). The quality controls (QC) are defined by bit flags such as mandatory quality assessment (QA) flags, data quality flags, emissivity quality flags, and cloud error flags.
To retrieve and pre-process the products, the MODIS R-package (MODIS acquisition and processing package v1.1.4) was used [61]. The package is run in the R software system and environment for statistical computing and graphics [62]. The MODIS R-package allows automatic downloading of data and processing such as changing file format, mosaicking, subsetting, and time-series filtering [63]. Using the package, digital numbers (DN) of MODIS Terra LST products were converted into LST (Table 1). Additionally, three terrain parameters (elevation, slope, and aspect) originating from SRTM-DEM [64] were retrieved. All raster data were re-projected to MODIS sinusoidal projection.

3.2. In Situ Meteorological Data

Sixty-three synoptic weather stations are present in Mongolia. Their geographical locations are indicated in Figure 1a. The weather stations provide Ta every 3 h, i.e., eight times a day. Air temperature data between 2002 and 2004 were obtained from the Mongolian Information Research Institute of Meteorology, Hydrology, and the Environment (IRIMHE). Data from 2004 to 2017 was downloaded from the “Reliable Prognosis (RP5)” website (https://rp5.ru). From the three-hourly meteorological data, the average air temperature was calculated for every 8 days of MODIS LST, taking into account the eight daily observations. This led to a total of 8544 meteorological data-points from 63 automatic weather stations covering the same period as the MODIS LST products, allowing for the development of prediction models between the remotely sensed data and Ta. The frequency distribution of the measured Ta reference data from 63 weather stations for the period 2002–2017 is shown in Figure 3 (n = 8544). The value of measured Ta ranged from −36.6 °C to 27.2 °C with a mean value of 0.7 °C and a standard deviation of 14.6 °C.

3.3. Random Forest and Partial Least Square Regression

RF and PLS models were trained to predict Ta using up to 17 predictor variables. The use of two competing approaches permits the evaluation of the benefits of using non-linear machine learning approaches (e.g., RF) compared to classical linear regression models (e.g., PLS). Twelve of the seventeen variables were derived from LST time-series products of the Terra MODIS for the period 2002–2017 (Table 1). The five remaining variables were elevation, slope, and aspect (extracted from SRTM-DEM), and geographical location (latitude and longitude) of weather stations (extracted from vector data). Summary descriptive statistics of the response and the 17 predictor variables are reported in Table 2. Longitude was included as this indirectly depicts (for Mongolia) distance to sea [65].
The correlation matrix (Figure 4) reveals a strong correlation between Ta and daytime/nighttime LST of MODIS, as well as several other correlations and redundancies. Based on these intercorrelations and taking into account that the number of variable sets should be relatively small, the predictor variables were grouped into seven different groupings (Table 3).

3.3.1. RF Regression

The well-known random forest regression method [42,44] was chosen as the main approach to model the relation between our response variable (Ta) and the predictor variables listed in Table 2 and Table 3 (LST MODIS products plus elevation, slope, aspect, latitude, and longitude). RF is a non-linear statistical ensemble method that leverages uncorrelated decision trees for regression. Developed by Breiman [66], it is capable of modeling discrete and/or continuous data sets [67,68]. RF predictions are obtained by aggregating a large number of individual regression decision trees where each decision tree is built from bootstrapped training samples (as in bagging) and variables are randomly selected at each decision node. The algorithm then randomly selects a subset of the predictors as candidates for splitting [66,69]. To obtain the final regression model, the results of all the individual trees are averaged. Good examples of the benefits and drawbacks of RF are given in [70,71,72].
The RF algorithm provides out-of-bag error (OOB) estimates and variable importance rankings [73,74], as not all observations are included in the respective bootstraps of the individual trees. In each tree at each split, the enhancement in the split-criterion importance measure is characterized by the splitting variable and aggregates individually all the trees in the forest for each variable [75]. Variable importance is measured by computing the increase in mean square error (MSE) when the OOB data for each variable are again computed but without the left-out variable [66,76]. The variable importance measures can assist in defining which variables are most important in the reduction of prediction error [71]. Two kinds of variable importance measures widely use the “randomForest” package in R [69,77]: (1) percent increase in the mean square error (%IncMSE) and (2) increase in node purity (IncNodePurity). From these, our analysis computed and analysed %IncMSE. However, we also checked the IncNodePurity indicator, but found similar results (not shown). In our research, the basic algorithm shown in Equation (1) was used to build the RF predictor for regression [75]:
f ^   ( x ) = 1 B b = 1 B T b ( x )
A new bootstrap sample for each decision tree Tb that includes X = x1…… xi with responses Y = y1…… yi bagging repeatedly (B times) selects a random sample from training data and each unpruned decision tree is increased in the sample. To increase each individual tree Tb, the following steps are repeated at each terminal node of the tree:
  • Randomly select m variables from p variables
  • Pick the variable that best splits and the corresponding split point
  • Split the node into two nodes.
As mentioned above, to implement the RF regression model, two parameters must be set: The number of decision trees (ntree) and the number of variables to select for the best split (mtry). For both hyperparameters, standard settings have been chosen. Each decision tree is independently increased to its maximum size, focusing on a new bootstrap sample from the training data (2/3 of samples). The remaining 1/3 of the samples, not used to fit the given decision tree, are referred to as the out-of-bag sample. The OOB sample is used to calculate the OOB error rate and variable importance. For quantifying the OOB error (prediction error) for each RF decision tree, we used Equation (2), i.e.,
E r r o r O O B =   1 n i = 1 n ( y i y ^ i ) 2
where y ^ i is the estimated output of OOB samples, y i is the actual output, and n is the total number of OOB samples. RF regression is flexible and easy to use in comparison to other machine learning algorithms, even without hyper-parameter tuning.

3.3.2. PLS Regression

For comparison—and to assess the differences between linear and non-linear models—a prominent linear modeling technique was used: Partial least square regression. PLS is widely used by the remote sensing community for vegetation analysis [78,79,80], soil related studies [81,82,83], and climate and ecological studies [84,85,86,87] amongst others.
PLSR is a multivariate linear regression method used to predict a set of dependent variables from a set of independent variables or predictors [88]. PLSR was originally developed for econometrics and chemometrics [89], where commonly a large number of strongly correlated predictor variables exist [90]. PLSR reduces the variables to a smaller set of uncorrelated components and performs least squares regression on these components instead of on the original data. Compared to other techniques, PLSR is more robust and less susceptible to data redundancy and over-fitting [91].
PLSR extracts a set of latent variables that explain the correlation between dependent and independent variables. The optimum number of latent variables for each generated model are implemented using the minimum value of residual mean squared error and the leave-one-out-cross-validation (LOO-CV) methods, e.g., jackknife and bootstrap [92]. To assess which variables are most contributing to the PLS model, we used the variable importance in projection (VIP) method [90], as seen in Equation (3), i.e.,
V I P j = f = 1 F w j f 2   S S Y f   J S S Y t o t a l   F
where V I P j   is a measure of the contribution of the j variable in the PLSR model, W j f   is the weight value for the j variable and f latent variables (components), S S Y f is the sum of squares of explained variance for the f latent variable and J number of the predictor (independent) variables, S S Y t o t a l is the total sum of squares explained as the response (dependent) variables, and F is the total number of latent variables. The VIP values determine the contribution of the predictor variables to the PLSR latent variables. A VIP value greater than 0.80 ensures that only relevant variables are considered [90]. In [93] the VIP threshold of predictor variables that were identified as the most relevant variables ranged between 0.83 and 1.21. Predictor variables with ≤0.80 VIP values were classified as less important while variables with VIP values ≥1.20 were considered the most influential.

3.3.3. Model Evaluation and Statistics

Two widely used statistics were calculated to assess the accuracy of the models [94], including the R2 and the RMSE. The R2 describes the percentage of explained variance whereas the RMSE summarizes the deviations of predictions from the one-to-one line.
As both models provide quantitative information about the importance of different variables, we also report these findings. For the RF regression model, we assessed the importance of the individual predictors in Ta estimates focused on the %IncMSE [39]. For the PLS regression model, we used the VIP method [90,93,95].

4. Results

4.1. Comparison of RF and PLS Models: Variable Importance and Prediction Accuracy

The estimated importance of the 17 predictor variables in the RF regression model is shown in Figure 5 for each of the twelve months. Under the top three in each month, LSTn appeared 12 times, LSTd 11 times, elevation 7 times, latitude 5 times and aspect once. The remaining variables were never found under the top three in those rankings. Strong seasonality in the ranking can also be observed. For example, the warm season Ta (April to October) was strongly dependent on elevation. Conversely, the cold season Ta was more heavily affected by latitude. Using all 17 predictor variables led to R2 values in the range from 0.83 (April) to 0.96 (August), while RMSE were between 0.91 °C (September) to 1.92 °C (February) (see Table 4, column G7).
Similar results were obtained for the PLS models (Appendix A Figure A1). Using PLS regression, the variables most often listed under the top three were LSTn (12 times), LSTd (12 times), elevation (4 times), latitude (4 times), and emissivity (4 times). Again, the ranking was season-dependent. The variables LSTn, LSTd, and latitude were the most important variables for estimating Ta in autumn and winter (September to February). The Ta for spring (March–May) was strongly dependent on LSTn, LSTd, and emissivity. For the summer months (June–August), LSTn, LSTd, and elevation were strongly influenced by the estimation of Ta for summer.
Using the entire set of 17 predictor variables for estimating the monthly average air temperatures, accuracies of PLS models were constantly lower compared to the RF models. The PLS models gave R2 a measured and estimated monthly Ta between 0.74 and 0.86 and RMSE from 1.20 °C to 2.19 °C (Appendix A Table A2, column G7).
Concerning variable importance, the three variables LSTn, LSTd, and elevation were identical for PLS and RF regression models. This shows that LSTn, LSTd, and elevation play a key role in modeling Ta, with all other variables having a significantly smaller impact.
To further study the impact of the different predictor variables, the seven variable groupings highlighted in Table 3 were analyzed in more detail. Results for each month and the annual average air temperature are shown in Table 4 for the RF models. Compared to the full set of 17 variables (column G7), the reduced set with only three predictor variables LSTn, LSTd, and elevation (G2) achieved comparable results, again highlighting and confirming the importance of these three predictor variables. None of the other five groupings (G1 and G3 to G5) were able to yield similar model performances. The same findings also hold for the PLS models (Appendix A Table A2 for details) but with constantly lower accuracies compared to the RF models. PLS models and groupings G1 to G7 were, therefore, skipped for the remainder of the study.

4.2. Maps of Predicted Air Temperatures Using RF Models with the Reduced Feature Set

Both the results of the variable importance rankings (Figure 5) and the grouping of variables (Table 4) indicate that relatively simple RF prediction models can be built to estimate Ta using only daytime/nighttime LST and elevation information. Scatterplots between measured and estimated monthly average air temperatures using only these three predictor variables are shown in Figure 6. Corresponding maps of modeled air temperatures at 1 km spatial resolution and covering the entire land-mass of Mongolia are shown in Figure 7a (see Figure A2 and Table A3 in the Appendix A for corresponding scatterplots and maps generated using PLS models).
The scatterplots in Figure 6 reveal that the RF-predicted Ta is well distributed around the 1-to-1 line, with no apparent systematic deviations. In particular, we do not see any autocorrelation in the errors, nor saturation effects. The errors are generally low and the explained variance (R2) mostly above 0.85. Generally, however, the RMSE increases slightly during the colder months (Figure 7b).
The maps in Figure 7a depict in high spatial detail the model predictions. As expected, the predicted air temperatures decrease with elevation (Figure 1) but reveal additional detail and information. Monthly analyses of the coefficient of determination (R2; in blue) and root mean square error (RMSE; in red) for the period 2002–2017 are shown in Figure 7b. Overall, a good agreement between observed and estimated Ta values was found but reflected again the afore-mentioned seasonal pattern. Large discrepancies were found to occur in transition months, such as the start or end of seasons.
The PLS-generated maps of monthly average air temperatures are generally similar to the maps derived from RF models (see Appendix A Figure A2b). However, a more detailed analysis reveals sometimes larger differences, even if modeled air temperatures are averaged by season. For example, Figure 8 clearly shows that large method-specific differences occur (maps in the third column). The differences show large seasonal fluctuations. The deviations moreover show a clear north–south gradient with generally lower Ta estimated using PLS compared to RF (reddish colors). The deviations are usually strongest during the warmer months. As the RF model outperformed the PLS model when evaluated against the observed air temperatures (Table 4 for RF and Appendix A Table A2 for PLS), we interpret these findings as mainly being the result of a systematic underestimation of Ta by the PLS model.

5. Discussion

Using relatively simple RF models driven by a few predictor variables, climatologies of monthly air temperatures in Mongolia could be obtained in this study with high accuracy (RMSE of about 0.84–1.93°C). Without any hyperparameter tuning, the non-linear RF models outperformed linear PLS models in accordance with other studies [43,55,96,97]. Amongst the variables studied, the MODIS-derived land surface temperatures (day and night) together with elevation were the three most important predictors. The studies of Noi et al. [41] and Kilibarda et al. [31] have also reported the high importance of day- and nighttime LST observations as well as elevation.
As LST (both day and night) and elevation can be readily produced at 1 km spatial resolution, the models calibrated against weather station data permitted the creation of maps of average air temperature for each of the twelve months in unprecedented detail and accuracy. Although the RF-generated maps often follow elevation, the inclusion of remotely sensed land surface temperature from MODIS improved the accuracy and spatial detail.
The results of the importance analysis indicated that nighttime LST was slightly more important compared to LSTd. The same result has been noted in China [42] and in Portugal [3]. In accordance with these studies, we argue that nighttime observations are probably more predictive because LSTn is not affected by reflected solar radiation when using TIR sensors [17]. The daytime land surface temperature was nevertheless found to be important, as it reveals the strength of the latent heat flux and the energy available for generating sensible heat [42].
Other studies have confirmed that Ta predictions are possible using satellite observations and that there is a strong relationship between Ta, LSTn, and LSTd [32]. Several studies have produced Ta estimations using MODIS LST data using multivariate linear and non-parametric regression methods [39,41,42,43,44,45,93]. These already published studies have showed different levels of success. The performance of multivariate and non-parametric regression models has been strongly dependent on environmental parameters such as vegetation cover, slope, aspect, elevation, quality of MODIS LST products, and applied filter techniques. For instance, the accuracy of the MODIS LST has been found to depend on the employed split-window algorithm, cloud cover, and terrain parameters [32]. Nonetheless, estimation of Ta derived from MODIS LST studies using multivariate and non-parametric algorithms is suitable for generating results at high accuracy. For instance, [42] estimated monthly average Ta for the territory of China at a spatial resolution of 1 km using RF regression with MODIS LST, normalized difference vegetation index (NDVI), nighttime light and elevation. Using this dataset, the RMSE of the monthly average Ta ranged between 1.57 °C to 1.99 °C. Our study has showed that monthly average Ta can be accurately estimated using LSTn, LSTd, and elevation with similar RMSE ranging from 0.91 °C to 1.93 °C. The method is relatively easy to implement provided that there is a sufficient amount of training data with corresponding EO time-series observations [96].

6. Conclusions

In this study, PLS and RF regression models were applied to estimate monthly average Ta in Mongolia for the period 2002–2017 using MODIS LST time-series products and terrain parameters. Meteorological data from 63 automatic weather stations were used to calibrate and validate the PLSR and RF models. Both models were trained to predict Ta using up to 17 variables as predictor variables. Twelve variables were derived from LST time-series products of Terra MODIS and three variables were extracted from an SRTM DEM (elevation, slope, and aspect). The geographical location (longitude and latitude) was used as an additional variable. For training, a total of 8544 meteorological data points from 63 automatic weather stations and corresponding MODIS LST were used. Both datasets covered the period 2002–2017. Using only day/nighttime LST and elevation as predictor variables, the correlation between measured and estimated monthly average Ta RMSE ranged from 1.20 °C to 2.19 °C for the PLSR and 0.84 °C to 1.93 °C for the RF. The significantly lower errors of the RF models confirm the benefits of this machine learning approach compared to traditional (linear) modeling techniques (e.g., PLSR). We therefore recommend the use of RF models for similar studies.
Concerning the MODIS land surface temperature data, we found that this information contributed significantly to the modeling of air temperature. For example, it was not possible to obtain similarly low errors in the modeled air temperature using only terrain parameters as predictors. It is recommended that day- and nighttime LST be used simultaneously as both variables scored high in the feature importance metric.
Both machine learning models (RF and PLSR) represented well seasonal and spatial variations in Ta when time-series of LST were included as predictor variables. Using the models, maps of the monthly average Ta of Mongolia were developed at a spatial resolution of 1 km which were representative for the period 2002–2017. Although errors in the predicted Ta were generally low, the residual errors showed a significant seasonality; the warmer months were generally better modeled compared to the extremely cold winter months. Probably, the increased errors during the winter months reflect a lower accuracy in the input (LST) data. Further research is warranted to better understand the seasonality of the model quality.
Despite these trends, we firmly believe that this spatial dataset may be useful for various environmental applications; for instance, it may be useful for better assessing bioclimatic variations within the huge land-mass of Mongolia. The developed methodology is relatively easy and transparent and can be applied in different geographic regions, provided that enough weather stations are available to permit a model calibration. The spatial resolution of the final map product mainly depends on the ground sampling distance of the employed satellite sensors. As sensor technology advances at a rapid pace, the current 1 km spatial resolution can be further improved in the near future.

Author Contributions

The first author M.O. analyzed the data and performed the experiments; she also runs the statistical analysis in R software and drafted the first version of the manuscript. The second author C.A. contributed to data analysis and the revision on an earlier and on the final version of the manuscript. The third author M.M. contributed MODIS data analysis preprocessing in R software. All authors together developed and discussed the manuscript and wrote the paper.

Funding

This article was funded by the Eurasia-Pacific Uninet, Ernst Mach Grant (reference no. ICM-2018-10238).

Acknowledgments

We are grateful to Anja Klisch and Valentin Pesendorfer, who are researchers at the Natural Resources and Life Sciences University (BOKU), for their assistance in writing R code. Munkhdulam Otgonbayar received funding from Eurasia-Pacific Uninet on behalf of the Austrian Federal Ministry of Science, Research and Economy (BMWFW) during the preparation of the manuscript. We thank all colleagues at the Institute of Geography and Geoecology, Mongolian Academy of Sciences who helped in this study. We would like to thank our reviewers for providing valuable feedback and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Variable importance in projections (VIPs), R2, and RMSE for the 17 predictor variables in the twelve-monthly PLSR models.
Figure A1. Variable importance in projections (VIPs), R2, and RMSE for the 17 predictor variables in the twelve-monthly PLSR models.
Remotesensing 11 02588 g0a1
Table A1. Model coefficients for the twelve-monthly PLS regression models including all 17 predictor variables. In the last column, the model of the annual average Ta is also included.
Table A1. Model coefficients for the twelve-monthly PLS regression models including all 17 predictor variables. In the last column, the model of the annual average Ta is also included.
VariableJanFebMarAprMayJunJulAugSepOctNovDecYear
Intercept−52.83−93.6226.1755.5293.6889.3797.2745.1168.7662.9253.5422.1782.62
LSTd0.2550.2190.2080.1850.1590.1460.1370.1780.1500.1280.1970.2340.319
LSTn0.4200.4570.4250.2650.2230.2320.2270.2840.2100.1960.3440.4070.464
CsD0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*
CsN−0.014−0.0030.0070.000*0.0010.0010.0020.000*0.0040.000*0.000*−0.0110.000*
DvA0.0150.0150.007−0.0020.0050.0080.0110.0210.0060.0220.0030.004−0.007
DvT0.0940.109−0.015−0.130−0.072−0.041−0.0190.003−0.1050.061−0.044−0.0210.095
Em310.0690.1040.021−0.152−0.133−0.126−0.131−0.050−0.127−0.125−0.0470.006−0.079
Em320.1010.1580.003−0.012−0.208−0.198−0.209−0.096−0.199−0.195−0.087−0.0090.007
NvA−0.019−0.034−0.034−0.012−0.009−0.002−0.007−0.0170.007−0.004−0.017−0.028−0.038
NvT0.0310.100−0.1140.0450.0530.0500.0340.0290.1690.083−0.030−0.038−0.413
QCd0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*
QCn0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*0.000*
Elevation0.000*0.000*−0.001−0.001−0.001−0.001−0.001−0.002−0.001−0.001−0.0010.000*−0.001
Slope−0.0210.028−0.009−0.061−0.059−0.062−0.065−0.047−0.051−0.037−0.025−0.0310.026
Aspect−0.002−0.001−0.001−0.002−0.002−0.001−0.002−0.001−0.002−0.002−0.003−0.0030.003
Latitude−0.411−0.401−0.248−0.268−0.236−0.214−0.246−0.244−0.260−0.248−0.241−0.3930.336
Longitude0.0990.1060.0390.0230.0420.0460.0520.0650.0480.0320.0160.0620.012
* <0.0001.
Table A2. PLS regression results. Summary statistics (R2 and RMSE) for the monthly Ta prediction models, including six groups of variables.
Table A2. PLS regression results. Summary statistics (R2 and RMSE) for the monthly Ta prediction models, including six groups of variables.
nG1G2G3G4G5G6G7
R2RMSER2RMSER2RMSER2RMSER2RMSER2RMSER2RMSE
January7120.871.960.871.950.344.340.822.260.394.170.862.010.852.10
February7120.832.190.832.190.324.400.792.450.374.260.852.050.842.19
March7120.801.890.741.940.363.400.772.040.393.320.811.840.811.83
April7120.791.510.791.530.482.390.771.580.322.730.751.670.741.69
May7120.761.480.791.410.741.540.801.370.312.520.761.500.771.47
June7120.781.440.801.380.751.520.791.410.262.630.781.450.761.51
July7120.831.330.861.200.791.480.841.280.282.710.821.340.801.44
August7120.841.360.871.230.811.440.851.280.302.800.831.370.861.25
September7120.811.350.841.240.821.300.831.260.322.530.801.370.791.42
October7120.831.270.831.260.681.730.821.310.412.360.811.320.771.47
November7120.831.540.831.570.372.970.791.700.402.920.821.590.811.65
December7120.861.680.861.670.363.530.821.890.343.580.851.740.831.82
Figure A2. Estimated monthly average Ta based on the PLS regression model and using LSTd, LSTn, and elevation as predictor variables. (a) Spatial maps of estimated monthly average Ta over Mongolia at 1 km spatial resolution. (b) Monthly statistics of R2 (blue) and RMSE (red) between observed and predicted air temperature.
Figure A2. Estimated monthly average Ta based on the PLS regression model and using LSTd, LSTn, and elevation as predictor variables. (a) Spatial maps of estimated monthly average Ta over Mongolia at 1 km spatial resolution. (b) Monthly statistics of R2 (blue) and RMSE (red) between observed and predicted air temperature.
Remotesensing 11 02588 g0a2
Table A3. Model equations obtained from the PLS regression models using three variables: LSTd, LSTn, and elevation. The months are numbered from 01 to 12.
Table A3. Model equations obtained from the PLS regression models using three variables: LSTd, LSTn, and elevation. The months are numbered from 01 to 12.
Regression ModelsR2RMSE
Ta01 = –2.137 + 0.347 × LSTd + 0.497 × LSTn + 0.00033 × elevation0.871.95
Ta02= –3.037 + 0.297 × LSTd + 0.493 × LSTn + 0.00019 × elevation0.832.19
Ta03 = –1.986 + 0.252 × LSTd + 0.477 × LSTn – 0.001 × elevation0.741.94
Ta04 = 0.516 + 0.296 × LSTd + 0.424 × LSTn − 0.002 × elevation0.791.53
Ta05 = 3.863 + 0.272 × LSTd + 0.383 × LSTn − 0.002 × elevation0.791.41
Ta06 = 7.060 + 0.242 × LSTd + 0.384 × LSTn − 0.002 × elevation0.801.38
Ta07 = 8.440 + 0.241 × LSTd + 0.398 × LSTn − 0.002 × elevation0.861.20
Ta08 = 7.644 + 0.253 × LSTd + 0.419 × LSTn − 0.002 × elevation0.871.23
Ta09 = 5.294 + 0.291 × LSTd + 0.407 × LSTn − 0.002 × elevation0.841.24
Ta10 = 3.418 + 0.266 × LSTd + 0.406 × LSTn − 0.002 × elevation0.831.26
Ta11 = –0.912 + 0.271 × LSTd + 0.411 × LSTn − 0.001 × elevation0.831.57
Ta12 = –2.560 + 0.314 × LSTd + 0.463 × LSTn + 0.00037 × elevation0.861.67

References

  1. Nieto, H.; Sandholt, I.; Aguado, I.; Chuvieco, E.; Stisen, S. Air temperature estimation with MSG-SEVIRI data: Calibration and validation of the TVX algorithm for the Iberian Peninsula. Remote Sens. Environ. 2011, 115, 107–116. [Google Scholar] [CrossRef] [Green Version]
  2. Shamir, E.; Georgakakos, K.P. MODIS Land Surface Temperature as an index of surface air temperature for operational snowpack estimation. Remote Sens. Environ. 2014, 152, 83–98. [Google Scholar] [CrossRef]
  3. Benali, A.; Carvalho, A.C.; Nunes, J.P.; Carvalhais, N.; Santos, A. Estimating air surface temperature in Portugal using MODIS LST data. Remote Sens. Environ. 2012, 124, 108–121. [Google Scholar] [CrossRef]
  4. Stisen, S.; Sandholt, I.; Nørgaard, A.; Fensholt, R.; Eklundh, L. Estimation of diurnal air temperature using MSG SEVIRI data in West Africa. Remote Sens. Environ. 2007, 110, 262–274. [Google Scholar] [CrossRef]
  5. Prihodko, L.; Goward, S.N. Estimation of air temperature from remotely sensed surface observations. Remote Sens. Environ. 1997, 60, 335–346. [Google Scholar] [CrossRef]
  6. Hooker, J.; Duveiller, G.; Cescatti, A. A global dataset of air temperature derived from satellite remote sensing and weather stations. Sci. Data 2018, 5, 180246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Peón, J.; Recondo, C.; Calleja, J.F. Improvements in the estimation of daily minimum air temperature in peninsular Spain using MODIS land surface temperature. Int. J. Remote Sens. 2014, 35, 5148–5166. [Google Scholar] [CrossRef]
  8. Lam, N.S.N. Spatial interpolation methods: A review. Am. Cartogr. 1983, 10, 129–150. [Google Scholar] [CrossRef]
  9. Thiébaux, H.J. Statistics and the environment: The analysis of large-scale earth-oriented systems. Environmetrics 1991, 2, 5–24. [Google Scholar] [CrossRef]
  10. Ishida, T.; Kawashima, S. Use of cokriging to estimate surface air temperature from elevation. Theor. Appl. Clim. 1993, 47, 147–157. [Google Scholar] [CrossRef]
  11. Hutchinson, M.F. A New Objective Method for Spatial Interpolation of Meteorological Variables from Irregular Networks Applied to the Estimation of Monthly Mean Solar Radiation, Temperature, Precipitation and Windrun; Need for climatic and hydrologic data in agriculture in South East Asia; UN University Workshop: Canberra, Australia, 1983; pp. 95–104. [Google Scholar]
  12. Willmott, C.J.; Rowe, C.M.; Philpot, W.D. Small-Scale Climate Maps: A Sensitivity Analysis of Some Common Assumptions Associated with Grid-Point Interpolation and Contouring. Am. Cartogr. 1985, 12, 5–16. [Google Scholar] [CrossRef]
  13. Willmott, C.J.; Robeson, S.M. Climatologically aided interpolation (CAI) of terrestrial air temperature. Int. J. Clim. 1995, 15, 221–229. [Google Scholar] [CrossRef]
  14. Burrough, P.A.; McDonnell, R.A. Principles of Geographical Information Systems; Oxford University Press: Oxford, UK, 1998. [Google Scholar]
  15. Mostovoy, G.V.; King, R.L.; Reddy, K.R.; Kakani, V.G.; Filippova, M.G. Statistical Estimation of Daily Maximum and Minimum Air Temperatures from MODIS LST Data over the State of Mississippi. GIScience Remote Sens. 2006, 43, 78–110. [Google Scholar] [CrossRef] [Green Version]
  16. Vogt, J.V.; Viau, A.A.; Paquet, F. Mapping regional air temperature fields using satellite-derived surface skin temperatures. Int. J. Clim. 1997, 17, 1559–1579. [Google Scholar] [CrossRef]
  17. Vancutsem, C.; Ceccato, P.; Dinku, T.; Connor, S.J. Evaluation of MODIS land surface temperature data to estimate air temperature in different ecosystems over Africa. Remote Sens. Environ. 2010, 114, 449–465. [Google Scholar] [CrossRef]
  18. Dash, P.; Göttsche, F.-M.; Olesen, F.-S.; Fischer, H. Land surface temperature and emissivity estimation from passive sensor data: Theory and practice-current trends. Int. J. Remote Sens. 2002, 23, 2563–2594. [Google Scholar] [CrossRef]
  19. Oke, T. The urban energy balance. Prog. Phys. Geogr. Earth Environ. 1988, 12, 471–508. [Google Scholar] [CrossRef]
  20. Chatterjee, R.; Singh, N.; Thapa, S.; Sharma, D.; Kumar, D. Retrieval of land surface temperature (LST) from landsat TM6 and TIRS data by single channel radiative transfer algorithm using satellite and ground-based inputs. Int. J. Appl. Earth Obs. Geoinf. 2017, 58, 264–277. [Google Scholar] [CrossRef]
  21. Li, Z.-L.; Tang, B.-H.; Wu, H.; Ren, H.; Yan, G.; Wan, Z.; Trigo, I.F.; Sobrino, J.A. Satellite-derived land surface temperature: Current status and perspectives. Remote Sens. Environ. 2013, 131, 14–37. [Google Scholar] [CrossRef] [Green Version]
  22. Prata, A.J.; Caselles, V.; Coll, C.; Sobrino, J.A.; Ottlé, C. Thermal remote sensing of land surface temperature from satellites: Current status and future prospects. Remote Sens. Rev. 1995, 12, 175–224. [Google Scholar] [CrossRef]
  23. Cresswell, M.P.; Morse, A.P.; Thomson, M.C.; Connor, S.J. Estimating surface air temperatures, from Meteosat land surface temperatures, using an empirical solar zenith angle model. Int. J. Remote Sens. 1999, 20, 1125–1132. [Google Scholar] [CrossRef]
  24. Stoll, M.J.; Brazel, A.J. Surface-air temperature relationships in the urban environment of phoenix, Arizona. Phys. Geogr. 1992, 13, 160–179. [Google Scholar] [CrossRef]
  25. Garand, L.; Buehner, M.; Wagneur, N. Background Error Correlation between Surface Skin and Air Temperatures: Estimation and Impact on the Assimilation of Infrared Window Radiances. J. Appl. Meteorol. 2004, 43, 1853–1863. [Google Scholar] [CrossRef]
  26. Jin, M.; Dickinson, R.E. Land surface skin temperature climatology: Benefitting from the strengths of satellite observations. Environ. Res. Lett. 2010, 5, 044004. [Google Scholar] [CrossRef]
  27. Zakšek, K.; Schroedter-Homscheidt, M. Parameterization of air temperature in high temporal and spatial resolution from a combination of the SEVIRI and MODIS instruments. ISPRS J. Photogramm. Remote Sens. 2009, 64, 414–421. [Google Scholar] [CrossRef]
  28. Sun, Y.J.; Wang, J.F.; Zhang, R.H.; Gillies, R.R.; Xue, Y.Y.C.B.; Bo, Y.C. Air temperature retrieval from remote sensing data based on thermodynamics. Theor. Appl. Climatol. 2005, 80, 37–48. [Google Scholar] [CrossRef]
  29. Zhu, W.; Lű, A.; Jia, S. Estimation of daily maximum and minimum air temperature using MODIS land surface temperature products. Remote Sens. Environ. 2013, 130, 62–73. [Google Scholar] [CrossRef]
  30. Czajkowski, K.; Goward, S.; Stadler, S.; Walz, A. Thermal Remote Sensing of Near Surface Environmental Variables: Application Over the Oklahoma Mesonet. Prof. Geogr. 2000, 52, 345–357. [Google Scholar] [CrossRef]
  31. Kilibarda, M.; Hengl, T.; Heuvelink, G.B.M.; Gräler, B.; Pebesma, E.; Tadić, M.P.; Bajat, B. Spatio-temporal interpolation of daily temperatures for global land areas at 1 km resolution. J. Geophys. Res. Atmos. 2014, 119, 2294–2313. [Google Scholar] [CrossRef]
  32. Chen, F.; Liu, Y.; Liu, Q.; Qin, F. A statistical method based on remote sensing for the estimation of air temperature in China. Int. J. Climatol. 2015, 35, 2131–2143. [Google Scholar] [CrossRef]
  33. Xu, Y.; Qin, Z.; Shen, Y. Study on the estimation of near-surface air temperature from MODIS data by statistical methods. Int. J. Remote Sens. 2012, 33, 7629–7643. [Google Scholar] [CrossRef]
  34. Yan, H.; Zhang, J.; Hou, Y.; He, Y. Estimation of air temperature from MODIS data in east China. Int. J. Remote Sens. 2009, 30, 6261–6275. [Google Scholar] [CrossRef]
  35. Bartkowiak, P.; Castelli, M.; Notarnicola, C. Downscaling Land Surface Temperature from MODIS Dataset with Random Forest Approach over Alpine Vegetated Areas. Remote Sens. 2019, 11, 1319. [Google Scholar] [CrossRef]
  36. Lu, L.; Zhang, T.; Wang, T.; Zhou, X. Evaluation of Collection-6 MODIS Land Surface Temperature Product Using Multi-Year Ground Measurements in an Arid Area of Northwest China. Remote Sens. 2018, 10, 1852. [Google Scholar] [CrossRef]
  37. Zhou, W.; Peng, B.; Shi, J.; Wang, T.; Dhital, Y.P.; Yao, R.; Yu, Y.; Lei, Z.; Zhao, R. Estimating High Resolution Daily Air Temperature Based on Remote Sensing Products and Climate Reanalysis Datasets over Glacierized Basins: A Case Study in the Langtang Valley, Nepal. Remote Sens. 2017, 9, 959. [Google Scholar] [CrossRef]
  38. Janatian, N.; Sadeghi, M.; Sanaeinejad, S.H.; Bakhshian, E.; Farid, A.; Hasheminia, S.M.; Ghazanfari, S. A statistical framework for estimating air temperature using MODIS land surface temperature data. Int. J. Climatol. 2017, 37, 1181–1194. [Google Scholar] [CrossRef]
  39. Ho, H.C.; Knudby, A.; Sirovyak, P.; Xu, Y.; Hodul, M.; Henderson, S.B. Mapping maximum urban air temperature on hot summer days. Remote Sens. Environ. 2014, 154, 38–45. [Google Scholar] [CrossRef]
  40. Duan, S.-B.; Li, Z.-L.; Tang, B.-H.; Wu, H.; Tang, R.; Bi, Y.; Zhou, G. Estimation of Diurnal Cycle of Land Surface Temperature at High Temporal and Spatial Resolution from Clear-Sky MODIS Data. Remote Sens. 2014, 6, 3247–3262. [Google Scholar] [CrossRef] [Green Version]
  41. Noi, P.T.; Kappas, M.; Degener, J. Estimating Daily Maximum and Minimum Land Air Surface Temperature Using MODIS Land Surface Temperature Data and Ground Truth Data in Northern Vietnam. Remote Sens. 2016, 8, 1002. [Google Scholar] [CrossRef]
  42. Li, L.; Zha, Y. Estimating monthly average temperature by remote sensing in China. Adv. Space Res. 2019, 63, 2345–2357. [Google Scholar] [CrossRef]
  43. Yoo, C.; Im, J.; Park, S.; Quackenbush, L.J. Estimation of daily maximum and minimum air temperatures in urban landscapes using MODIS time series satellite data. ISPRS J. Photogramm. Remote Sens. 2018, 137, 149–162. [Google Scholar] [CrossRef]
  44. Yang, Y.Z.; Cai, W.H.; Yang, J. Evaluation of MODIS Land Surface Temperature Data to Estimate Near-Surface Air Temperature in Northeast China. Remote Sens. 2017, 9, 410. [Google Scholar] [CrossRef]
  45. Meyer, H.; Katurji, M.; Appelhans, T.; Müller, M.U.; Nauss, T.; Roudier, P.; Zawar-Reza, P. Mapping Daily Air Temperature for Antarctica Based on MODIS LST. Remote Sens. 2016, 8, 732. [Google Scholar] [CrossRef]
  46. Emamifar, S.; Rahimikhoob, A.; Noroozi, A.A. Daily mean air temperature estimation from MODIS land surface temperature products based on M5 model tree. Int. J. Clim. 2013, 33, 3174–3181. [Google Scholar] [CrossRef]
  47. Zhang, H.; Zhang, F.; Zhang, G.; He, X.; Tian, L. Evaluation of cloud effects on air temperature estimation using MODIS LST based on ground measurements over the Tibetan Plateau. Atmos. Chem. Phys. Discuss. 2016, 16, 13681–13696. [Google Scholar] [CrossRef] [Green Version]
  48. Keramitsoglou, I.; Kiranoudis, C.T.; Sismanidis, P.; Zakšek, K. An Online System for Nowcasting Satellite Derived Temperatures for Urban Areas. Remote Sens. 2016, 8, 306. [Google Scholar] [CrossRef]
  49. Nanzad, L.; Zhang, J.; Tuvdendorj, B.; Nabil, M.; Zhang, S.; Bai, Y. NDVI anomaly for drought monitoring and its correlation with climate factors over Mongolia from 2000 to 2016. J. Arid. Environ. 2019, 164, 69–77. [Google Scholar] [CrossRef]
  50. Nandintsetseg, B.; Greene, J.S.; Goulden, C.E. Trends in extreme daily precipitation and temperature near lake Hövsgöl, Mongolia. Int. J. Clim. 2007, 27, 341–347. [Google Scholar] [CrossRef]
  51. World Map of the Köppen Climate Classification. Available online: https://commons.wikimedia.org/wiki/File:World_K%C3%B6ppen_Classification_(with_authors).svg (accessed on 18 October 2019).
  52. Climate Data Org Webpage. Available online: https://en.climate-data.org/asia/mongolia-180/ (accessed on 17 October 2019).
  53. Duan, S.-B.; Li, Z.-L.; Wu, H.; Leng, P.; Gao, M.; Wang, C. Radiance-based validation of land surface temperature products derived from Collection 6 MODIS thermal infrared data. Int. J. Appl. Earth Obs. Geoinf. 2018, 70, 84–92. [Google Scholar] [CrossRef]
  54. Oyler, J.W.; Dobrowski, S.Z.; Holden, Z.A.; Running, S.W. Remotely Sensed Land Skin Temperature as a Spatial Predictor of Air Temperature across the Conterminous United States. J. Appl. Meteorol. Clim. 2016, 55, 1441–1457. [Google Scholar] [CrossRef]
  55. Xu, Y.; Knudby, A.; Ho, H.C. Estimating daily maximum air temperature from MODIS in British Columbia, Canada. Int. J. Remote Sens. 2014, 35, 8108–8121. [Google Scholar] [CrossRef]
  56. Wan, Z.; Dozier, J. A generalized split-window algorithm for retrieving land-surface temperature from space. IEEE Trans. Geosci. Remote Sens. 1996, 34, 892–905. [Google Scholar] [Green Version]
  57. Wan, Z. New refinements and validation of the MODIS Land-Surface Temperature/Emissivity products. Remote Sens. Environ. 2008, 112, 59–74. [Google Scholar] [CrossRef]
  58. Wan, Z. New refinements and validation of the collection-6 MODIS land-surface temperature/emissivity product. Remote Sens. Environ. 2014, 140, 36–45. [Google Scholar] [CrossRef]
  59. Wan, Z.; Zhang, Y.; Zhang, Q.; Li, Z.-L. Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data. Remote Sens. Environ. 2002, 83, 163–180. [Google Scholar] [CrossRef]
  60. Wan, Z. Collection-6 MODIS Land Surface Temperature Products Users’ Guide; ICESS: Santa Barbara, CA, USA, 2007. [Google Scholar]
  61. Mattiuzzi, M.; Verbesselt, J.; Klisch, A.; Stevens, F.; Mosher, S.; Evans, B.; Lobo, A.; Detsch, F. MODIS: MODIS Acquisition and Processing Package. R Package Version 1.1.4. Available online: https://cran.r-project.org/src/contrib/Archive/MODIS/MODIS.pdf (accessed on 25 October 2018).
  62. R Core Team. R: A Language and Environment for Statistical Computing. Available online: https://www.R-project.org/ (accessed on 31 October 2019).
  63. Mattiuzzi, M.; Verbesselt, J.; Hengel, T.; Klisch, A.; Stevens, F.; Mosher, S.; Evans, B.; Lobo, A.; Hufkens, K.; Detsch, F. MODIS: MODIS Acquisition and Processing package. R Package Version 1.1.5. Available online: https://cran.r-project.org/web/packages/MODIS/MODIS.pdf (accessed on 8 March 2019).
  64. Rabus, B.; Eineder, M.; Roth, A.; Bamler, R. The shuttle radar topography mission—A new class of digital elevation models acquired by spaceborne radar. ISPRS J. Photogramm. Remote Sens. 2003, 57, 241–262. [Google Scholar] [CrossRef]
  65. Bulgan, D.; Mandakh, N.; Odbayar, M.; Otgontugs, M.; Tsogtbaatar, J.; Elbegjargal, N.; Enerel, T.; Erdenetuya, M. Desertification Atlas of Mongolia; Tsogtbaatar, J., Khudulmur, S., Eds.; Admon Press: Ulaanbaatar, Mongolia, 2013; pp. 17–29. ISBN 978-9973-0-065-2. [Google Scholar]
  66. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  67. Mutanga, O.; Adam, E.; Cho, M.A. High density biomass estimation for wetland vegetation using WorldView-2 imagery and random forest regression algorithm. Int. J. Appl. Earth Obs. Geoinf. 2012, 18, 399–406. [Google Scholar] [CrossRef]
  68. Pal, M. Random forest classifier for remote sensing classification. Int. J. Remote Sens. 2005, 26, 217–222. [Google Scholar] [CrossRef]
  69. Liaw, A.; Wiener, M. Classification and regression by randomForest. R News 2002, 2, 18–22. [Google Scholar]
  70. Immitzer, M.; Stepper, C.; Bock, S.; Straub, C.; Atzberger, C. Use of WorldView-2 stereo imagery and National Forest Inventory data for wall-to-wall mapping of growing stock. For. Ecol. Manag. 2016, 359, 232–246. [Google Scholar] [CrossRef]
  71. Belgiu, M.; Drăguț, L. Random forest in remote sensing: A review of applications and future directions. ISPRS J. Photogramm. Remote Sens. 2016, 114, 24–31. [Google Scholar] [CrossRef]
  72. Hudak, A.T.; Crookston, N.L.; Evans, J.S.; Hall, D.E.; Falkowski, M.J. Nearest neighbor imputation of species-level, plot-scale forest structure attributes from LiDAR data. Remote Sens. Environ. 2008, 112, 2232–2245. [Google Scholar] [CrossRef] [Green Version]
  73. Rodriguez-Galiano, V.F.; Ghimire, B.; Rogan, J.; Olmo, M.C.; Rigol-Sanchez, J.P. An assessment of the effectiveness of a random forest classifier for land-cover classification. ISPRS J. Photogramm. Remote Sens. 2012, 67, 93–104. [Google Scholar] [CrossRef]
  74. Genuer, R.; Poggi, J.-M.; Tuleau-Malot, C. Variable selection using random forests. Pattern Recognit. Lett. 2010, 31, 2225–2236. [Google Scholar] [CrossRef] [Green Version]
  75. Trevor, H.; Tibshirani, R.; Friedman, J. The Elements of Statistical Learning; Springer: Stanford, CA, USA, 2011. [Google Scholar]
  76. Prasad, A.M.; Iverson, L.R.; Liaw, A. Newer Classification and Regression Tree Techniques: Bagging and Random Forests for Ecological Prediction. Ecosystems 2006, 9, 181–199. [Google Scholar] [CrossRef]
  77. Ishwaran, H. Variable importance in binary regression trees and forests. Electron. J. Stat. 2007, 1, 519–537. [Google Scholar] [CrossRef]
  78. Atzberger, C.; Guérif, M.; Baret, F.; Werner, W. Comparative analysis of three chemometric techniques for the spectroradiometric assessment of canopy chlorophyll content in winter wheat. Comput. Electron. Agric. 2010, 73, 165–173. [Google Scholar] [CrossRef]
  79. Laurin, G.V.; Chen, Q.; Lindsell, J.A.; Coomes, D.A.; Del Frate, F.; Guerriero, L.; Pirotti, F.; Valentini, R. Above ground biomass estimation in an African tropical forest with lidar and hyperspectral data. ISPRS J. Photogramm. Remote Sens. 2014, 89, 49–58. [Google Scholar] [CrossRef]
  80. Inoue, Y.; Sakaiya, E.; Zhu, Y.; Takahashi, W. Diagnostic mapping of canopy nitrogen content in rice based on hyperspectral measurements. Remote Sens. Environ. 2012, 126, 210–221. [Google Scholar] [CrossRef]
  81. Gomez, C.; Rossel, R.A.V.; McBratney, A.B. Soil organic carbon prediction by hyperspectral remote sensing and field vis-NIR spectroscopy: An Australian case study. Geoderma 2008, 146, 403–411. [Google Scholar] [CrossRef]
  82. Farifteh, J.; Van Der Meer, F.; Atzberger, C.; Carranza, E.; Carranza, E.J. Quantitative analysis of salt-affected soil reflectance spectra: A comparison of two adaptive methods (PLSR and ANN). Remote Sens. Environ. 2007, 110, 59–78. [Google Scholar] [CrossRef]
  83. Rossel, R.V.; Walvoort, D.; McBratney, A.; Janik, L.; Skjemstad, J.; Rossel, R.V. Visible, near infrared, mid infrared or combined diffuse reflectance spectroscopy for simultaneous assessment of various soil properties. Geoderma 2006, 131, 59–75. [Google Scholar] [CrossRef]
  84. Fernandes, L.S.; Fernandes, A.; Ferreira, A.; Cortes, R.; Pacheco, F. A partial least squares—Path modeling analysis for the understanding of biodiversity loss in rural and urban watersheds in Portugal. Sci. Total Environ. 2018, 626, 1069–1085. [Google Scholar] [CrossRef] [PubMed]
  85. Liu, Y.; Peng, J.; Wang, Y. Application of partial least squares regression in detecting the important landscape indicators determining urban land surface temperature variation. Landsc. Ecol. 2018, 33, 1133–1145. [Google Scholar] [CrossRef]
  86. Ceglar, A.; Toreti, A.; Lecerf, R.; Van Der Velde, M.; Dentener, F. Impact of meteorological drivers on regional inter-annual crop yield variability in France. Agric. For. Meteorol. 2016, 216, 58–67. [Google Scholar] [CrossRef]
  87. Carrascal, L.M.; Galván, I.; Gordo, O. Partial least squares regression as an alternative to current regression methods used in ecology. Oikos 2009, 118, 681–690. [Google Scholar] [CrossRef]
  88. Abdi, H. Partial least square regression (PLS regression). Encycl. Res. Methods Soc. Sci. 2003, 6, 1–7. [Google Scholar]
  89. Wold, H. Soft modeling: The basic design and some extensions. Syst. Indirect Obs. 1982, 2, 343. [Google Scholar]
  90. Wold, S.; Sjöström, M.; Eriksson, L. PLS-regression: A basic tool of chemometrics. Chemom. Intell. Lab Syst. 2001, 58, 109–130. [Google Scholar] [CrossRef]
  91. Geladi, P.; Kowalski, B.R. Partial least-squares regression: A tutorial. Anal. Chim. Acta 1986, 185, 1–17. [Google Scholar] [CrossRef]
  92. Mevik, B.; Cederkvist, H.R. Mean squared error of prediction (MSEP) estimates for principal component regression (PCR) and partial least squares regression (PLSR). J. Chemom. 2004, 18, 422–429. [Google Scholar] [CrossRef]
  93. Mkhabela, M.S.; Bullock, P.R.; Sapirstein, H.D. Characterising the most critical climatic parameters that impact the quality of spring-wheat (Triticum aestivum L.) on the Canadian Prairies using partial least squares (PLS) analysis. J. Cereal Sci. 2018, 81, 44–51. [Google Scholar] [CrossRef]
  94. Richter, K.; Atzberger, C.; Hank, T.B.; Mauser, W. Derivation of biophysical variables from Earth observation data: Validation and statistical measures. J. Appl. Remote Sens. 2012, 6, 63557. [Google Scholar] [CrossRef]
  95. Farres, M.; Platikanov, S.; Tsakovski, S.; Tauler, R. Comparison of the variable importance in projection (VIP) and of the selectivity ratio (SR) methods for variable selection and interpretation. J. Chemom. 2015, 29, 528–536. [Google Scholar] [CrossRef]
  96. Otgonbayar, M.; Atzberger, C.; Chambers, J.; Damdinsuren, A. Mapping pasture biomass in Mongolia using Partial Least Squares, Random Forest regression and Landsat 8 imagery. Int. J. Remote Sens. 2019, 40, 3204–3226. [Google Scholar] [CrossRef]
  97. Sun, J.; Passi, K.; Jain, C.K. Improved microarray data analysis using feature selection methods with machine learning methods. In Proceedings of the 2016 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Shenzhen, China, 15–18 December 2016; pp. 1527–1534. [Google Scholar]
Figure 1. Mongolia: (a) Topography and location of meteorological stations (n = 63) for which reference air temperature information from automatic weather stations was available. The digital elevation model (DEM) was derived from the Shuttle Radar Topography Mission (SRTM-DEM) with a resolution of 90 m. (b) Köppen climate classification of Mongolia based on long-term meteorological conditions (>30 years) [51].
Figure 1. Mongolia: (a) Topography and location of meteorological stations (n = 63) for which reference air temperature information from automatic weather stations was available. The digital elevation model (DEM) was derived from the Shuttle Radar Topography Mission (SRTM-DEM) with a resolution of 90 m. (b) Köppen climate classification of Mongolia based on long-term meteorological conditions (>30 years) [51].
Remotesensing 11 02588 g001
Figure 2. Climate graphs of the three most prevalent climatic zones in Mongolia [52]. (a) Dalanzadgad with cold desert climate (BWk), (b) Ulaangom with cold semi-arid climate (BSk), and (c) Sukhbaatar with cool continental climate/subarctic climate (Dwc). Information about altitude, average annual air temperature, and total precipitation is provided in each sub-plot. The data reflect long-term averages (>30 years).
Figure 2. Climate graphs of the three most prevalent climatic zones in Mongolia [52]. (a) Dalanzadgad with cold desert climate (BWk), (b) Ulaangom with cold semi-arid climate (BSk), and (c) Sukhbaatar with cool continental climate/subarctic climate (Dwc). Information about altitude, average annual air temperature, and total precipitation is provided in each sub-plot. The data reflect long-term averages (>30 years).
Remotesensing 11 02588 g002
Figure 3. Frequency distribution of measured average 8-day air temperature reference data (n = 8544) from the 63 automatic weather stations for the years 2002–2017. Monthly statistics are depicted in Table 2.
Figure 3. Frequency distribution of measured average 8-day air temperature reference data (n = 8544) from the 63 automatic weather stations for the years 2002–2017. Monthly statistics are depicted in Table 2.
Remotesensing 11 02588 g003
Figure 4. Correlation matrix between response and predictor variables (n = 8544). The saturation of the colors indicates the strength of the correlations. Positive correlations are shown in blue and negative correlations in red. In this graph, the air temperature data has been pooled across the twelve months. For the abbreviations, see Table 1.
Figure 4. Correlation matrix between response and predictor variables (n = 8544). The saturation of the colors indicates the strength of the correlations. Positive correlations are shown in blue and negative correlations in red. In this graph, the air temperature data has been pooled across the twelve months. For the abbreviations, see Table 1.
Remotesensing 11 02588 g004
Figure 5. Random forest (RF) variable importance for each month. The importance is here given as the percentage increase in mean square error (%IncMSE).
Figure 5. Random forest (RF) variable importance for each month. The importance is here given as the percentage increase in mean square error (%IncMSE).
Remotesensing 11 02588 g005
Figure 6. Comparison between measured and estimated monthly average Ta using LSTd, LSTn, and elevation for the RF regression model.
Figure 6. Comparison between measured and estimated monthly average Ta using LSTd, LSTn, and elevation for the RF regression model.
Remotesensing 11 02588 g006
Figure 7. Estimated monthly average Ta based on RF regression model using LSTd, LSTn, and elevation as predictor variables. (a) Spatial maps of estimated monthly average Ta over Mongolia at 1 km spatial resolution. (b) Monthly statistics of R2 (blue) and RMSE (red) between observed and predicted air temperature.
Figure 7. Estimated monthly average Ta based on RF regression model using LSTd, LSTn, and elevation as predictor variables. (a) Spatial maps of estimated monthly average Ta over Mongolia at 1 km spatial resolution. (b) Monthly statistics of R2 (blue) and RMSE (red) between observed and predicted air temperature.
Remotesensing 11 02588 g007
Figure 8. Estimated average Ta per season using LSTd, LSTn, and elevation as predictor variables. (a) Spatial maps of seasonal-average Ta over Mongolia at 1 km spatial resolution using the RF (first column) and partial least squares (PLS) regression models (second column). In the last column, the difference between the two model outputs is shown.
Figure 8. Estimated average Ta per season using LSTd, LSTn, and elevation as predictor variables. (a) Spatial maps of seasonal-average Ta over Mongolia at 1 km spatial resolution using the RF (first column) and partial least squares (PLS) regression models (second column). In the last column, the difference between the two model outputs is shown.
Remotesensing 11 02588 g008
Table 1. Description of Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) products used in this study (source: Land Processes Distributed Active Archive Center (LP DAAC), 2019).
Table 1. Description of Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) products used in this study (source: Land Processes Distributed Active Archive Center (LP DAAC), 2019).
Variable TypeAcronymUnitsData TypeFill ValueValid Range (VR)Scale Factor (SF)Additional Offset (AO)
Daytime LSTLSTdKelvin16 bit07500 to 65,5350.02N/A
Nighttime LSTLSTnKelvin16 bit07500 to 65,5350.02N/A
Day clear sky coverage CsDN/A16 bit01 to 65,5350.0005N/A
Night clear sky coverageCsNN/A16 bit01 to 65,5350.0005N/A
View zenith angle of daytimeDvADegree8 bit2550 to 1301.0−65
View zenith angle of nighttimeNvADegree8 bit2550 to 1301.0−65
Time of daytime (local solar)DvTHours8 bit2550 to 2400.1N/A
Time of nighttime (local solar)NvTHours8 bit2550 to 2400.1N/A
Emissivity band 31Em31None8 bit01 to 2550.0020.49
Emissivity band 32Em32None8 bit01 to 2550.0020.49
Quality controls of day LSTQCdBit8 bitN/A0 to 255N/AN/A
Quality controls of night LSTQCnBit8 bitN/A0 to 255N/AN/A
Note: Digital number (DN) = VR × SF − AO.
Table 2. List of response/predictor variables and corresponding descriptive statistics (period 2002 to 2017). The list includes the measured air temperature (Ta) reference data at the weather station level (n = 712 for each of the twelve months) as well as the corresponding seventeen predictor variables extracted from satellite and other geo-data. For the acronyms of the variables, see Table 1.
Table 2. List of response/predictor variables and corresponding descriptive statistics (period 2002 to 2017). The list includes the measured air temperature (Ta) reference data at the weather station level (n = 712 for each of the twelve months) as well as the corresponding seventeen predictor variables extracted from satellite and other geo-data. For the acronyms of the variables, see Table 1.
VariableNo. of Samples (n)MinimumMaximumMeanStandard Deviation
Ta01712−36.60−6.60−20.835.36
Ta02712−35.10−0.60−16.605.38
Ta03712−20.504.80−6.734.22
Ta04712−7.6012.403.923.32
Ta057122.8019.3010.753.05
Ta067129.2024.9016.983.06
Ta0771211.6027.2019.493.21
Ta087128.7025.6017.163.34
Ta097122.6019.7010.483.07
Ta10712−8.2011.701.493.08
Ta11712−22.700.30−9.623.76
Ta12712−31.50−6.10−17.854.40
LSTd8544−36.9048.6013.4320.08
LSTn8544−42.5024.40−5.7514.47
CsD85440.000.130.060.02
CsN85440.000.130.070.02
DvA8544−55.0062.005.04−52.35
DvT854410.4012.1011.820.71
Em3185440.960.990.980.50
Em3285440.970.990.980.49
NvA8544−65.0056.00−0.48−56.47
NvT854420.8022.7021.901.29
QCd85442.00133.0061.9816.58
QCn85442.00145.0055.1419.70
Elevation63667.002255.001369.10411.70
Slope630.0819.60N/AN/A
Aspect636.34358.10N/AN/A
Latitude6342.9751.11N/AN/A
Longitude6389.93118.67N/AN/A
Table 3. Seven model subsets studied. The seven groups were generated to study the relations between responses and up to 17 predictor variables. Nvar indicates the number of predictor variables in each group.
Table 3. Seven model subsets studied. The seven groups were generated to study the relations between responses and up to 17 predictor variables. Nvar indicates the number of predictor variables in each group.
AcronymVariablesNvar
Group 1G1LSTd and LSTn2
Group 2G2LSTd, LSTn, and elevation3
Group 3G3Elevation, slope, aspect, latitude, and longitude5
Group 4G4Combined G1 and G37
Group 5G5CsD, CsN, DvA, DvT, Em31, Em32, NvA, NvT, QCd, and QCn10
Group 6G6Combined G1 and G512
Group 7G7Combined G1, G3, and G517
Table 4. Modeling results obtained using the RF regression. Reported are the monthly summary statistics (coefficient of determination (R2) and root-mean-square error (RMSE)) for Ta prediction models for each of the seven groups of variables. For details of groupings G1 to G7, see Table 3.
Table 4. Modeling results obtained using the RF regression. Reported are the monthly summary statistics (coefficient of determination (R2) and root-mean-square error (RMSE)) for Ta prediction models for each of the seven groups of variables. For details of groupings G1 to G7, see Table 3.
nG1G2G3G4G5G6G7
R2RMSER2RMSER2RMSER2RMSER2RMSER2RMSER2RMSE
January7120.852.070.911.590.523.710.891.770.404.140.881.860.901.71
February7120.832.240.871.930.493.830.862.050.463.950.852.060.881.92
March7120.762.090.851.640.552.840.831.770.463.120.811.830.851.67
April7120.771.590.821.410.472.420.801.470.382.620.791.510.831.40
May7120.761.480.881.050.751.520.841.070.392.370.811.340.881.08
June7120.771.450.881.040.761.500.881.070.312.530.811.330.881.07
July7120.841.290.930.840.811.410.910.970.322.650.871.170.870.93
August7120.831.380.950.920.831.370.900.970.352.680.861.270.960.92
September7120.811.320.910.910.841.220.881.050.402.370.851.190.910.91
October7120.831.280.900.990.701.680.881.080.472.230.851.220.891.03
November7120.811.650.871.340.542.550.851.470.422.860.841.520.861.39
December7120.841.760.891.440.602.790.881.550.403.430.861.640.891.49

Share and Cite

MDPI and ACS Style

Otgonbayar, M.; Atzberger, C.; Mattiuzzi, M.; Erdenedalai, A. Estimation of Climatologies of Average Monthly Air Temperature over Mongolia Using MODIS Land Surface Temperature (LST) Time Series and Machine Learning Techniques. Remote Sens. 2019, 11, 2588. https://doi.org/10.3390/rs11212588

AMA Style

Otgonbayar M, Atzberger C, Mattiuzzi M, Erdenedalai A. Estimation of Climatologies of Average Monthly Air Temperature over Mongolia Using MODIS Land Surface Temperature (LST) Time Series and Machine Learning Techniques. Remote Sensing. 2019; 11(21):2588. https://doi.org/10.3390/rs11212588

Chicago/Turabian Style

Otgonbayar, Munkhdulam, Clement Atzberger, Matteo Mattiuzzi, and Avirmed Erdenedalai. 2019. "Estimation of Climatologies of Average Monthly Air Temperature over Mongolia Using MODIS Land Surface Temperature (LST) Time Series and Machine Learning Techniques" Remote Sensing 11, no. 21: 2588. https://doi.org/10.3390/rs11212588

APA Style

Otgonbayar, M., Atzberger, C., Mattiuzzi, M., & Erdenedalai, A. (2019). Estimation of Climatologies of Average Monthly Air Temperature over Mongolia Using MODIS Land Surface Temperature (LST) Time Series and Machine Learning Techniques. Remote Sensing, 11(21), 2588. https://doi.org/10.3390/rs11212588

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop