基于简化参数方法的蒙古干旱区土壤湿度被动微波遥感

卫星被动微波遥感土壤湿度,是准确分析大空间尺度上陆表水分变化信息的有效手段.美国航天局(NASA)发布的基于AMSR-E观测亮温资料的全球土壤湿度反演产品,在蒙古干旱区的实际精度并不令人满意.本文基于对地表微波辐射传输中地表粗糙度和植被层影响的简化处理方法,采用AMSR-E的6.9 GHz,10.7 GHz和18.7 GHz之V极化亮温资料,应用多频率反演算法,并以国际能量和水循环协同观测计划(The Coordinated Energy and Water Cycle Observations Project)即CEOP实验在蒙古国东部荒漠地区的地面实验资料作为先验知识,获取被动微波遥感模型的优化参数,以期获得蒙古干旱区精度更高的土壤湿度遥感估算结果.分析表明,本文方法反演的白天和夜间土壤湿度结果与地面验证值之间的均方根误差(RMSE)接近0.030 cm³/cm³, 证明所用方法在不需要其他辅助资料或参数帮助下,可较精确地反演干旱区表层土壤湿度信息,能够全天候、动态监测大空间尺度的土壤湿度变化,可为干旱区气候变化研究及陆面过程模拟和数据同化研究提供高精度的表层土壤湿度初始场资料.     

A soil moisture assimilation scheme based on the ensemble Kalman filter using microwave brightness temperature

This study presents a soil moisture assimilation scheme, which could assimilate microwave brightness temperature directly, based on the ensemble Kalman filter and the shuffled complex evolution method (SCE-UA). It uses the soil water model of the land surface model CLM3.0 as the forecast operator, and a radiative transfer model (RTM) as the observation operator in the assimilation system. The assimilation scheme is implemented in two phases: the parameter calibration phase and the pure soil moisture assimilation phase. The vegetation optical thickness and surface roughness parameters in the RTM are calibrated by SCE-UA method and the optimal parameters are used as the final model parameters of the observation operator in the assimilation phase. The ideal experiments with synthetic data indicate that this scheme could significantly improve the simulation of soil moisture at the surface layer. Furthermore, the estimation of soil moisture in the deeper layers could also be improved to a certain extent. The real assimilation experiments with AMSR-E brightness temperature at 10.65 GHz (vertical polarization) show that the root mean square error (RMSE) of soil moisture in the top layer (0–10 cm) by assimilation is 0.03355 m³ · m⁻³, which is reduced by 33.6% compared with that by simulation (0.05052 m³ · m⁻³). The mean RMSE by assimilation for the deeper layers (10–50 cm) is also reduced by 20.9%. All these experiments demonstrate the reasonability of the assimilation scheme developed in this study.     

Calibration of an integrated land surface process and radiobrightness (LSP/R) model during summertime

In this study, a soil vegetation and atmosphere transfer (SVAT) model was linked with a microwave emission model to simulate microwave signatures for different terrain during summertime, when the energy and moisture fluxes at the land surface are strong. The integrated model, land surface process/radiobrightness (LSP/R), was forced with weather and initial conditions observed during a field experiment. It simulated the fluxes and brightness temperatures for bare soil and brome grass in the Northern Great Plains. The model estimates of soil temperature and moisture profiles and terrain brightness temperatures were compared with the observed values. Overall, the LSP model provides realistic estimates of soil moisture and temperature profiles to be used with a microwave model. The maximum mean differences and standard deviations between the modeled and the observed temperatures (canopy and soil) were 2.6 K and 6.8 K, respectively; those for the volumetric soil moisture were 0.9% and 1.5%, respectively. Brightness temperatures at 19 GHz matched well with the observations for bare soil, when a rough surface model was incorporated indicating reduced dielectric sensitivity to soil moisture by surface roughness. The brightness temperatures of the brome grass matched well with the observations indicating that a simple emission model was sufficient to simulate accurate brightness temperatures for grass typical of that region and surface roughness was not a significant issue for grass-covered soil at 19 GHz. Such integrated SVAT-microwave models allow for direct assimilation of microwave observations and can also be used to understand sensitivity of microwave signatures to changes in weather forcings and soil conditions for different terrain types.     

Improvement of snow depth retrieval for FY3B-MWRI in China

The primary objective of this work is to develop an operational snow depth retrieval algorithm for the FengYun3B Microwave Radiation Imager(FY3B-MWRI)in China.Based on 7-year(2002–2009)observations of brightness temperature by the Advanced Microwave Scanning Radiometer-EOS(AMSR-E)and snow depth from Chinese meteorological stations,we develop a semi-empirical snow depth retrieval algorithm.When its land cover fraction is larger than 85%,we regard a pixel as pure at the satellite passive microwave remote-sensing scale.A 1-km resolution land use/land cover(LULC)map from the Data Center for Resources and Environmental Sciences,Chinese Academy of Sciences,is used to determine fractions of four main land cover types(grass,farmland,bare soil,and forest).Land cover sensitivity snow depth retrieval algorithms are initially developed using AMSR-E brightness temperature data.Each grid-cell snow depth was estimated as the sum of snow depths from each land cover algorithm weighted by percentages of land cover types within each grid cell.Through evaluation of this algorithm using station measurements from 2006,the root mean square error(RMSE)of snow depth retrieval is about 5.6 cm.In forest regions,snow depth is underestimated relative to ground observation,because stem volume and canopy closure are ignored in current algorithms.In addition,comparison between snow cover derived from AMSR-E and FY3B-MWRI with Moderate-resolution Imaging Spectroradiometer(MODIS)snow cover products(MYD10C1)in January 2010 showed that algorithm accuracy in snow cover monitoring can reach 84%.Finally,we compared snow water equivalence(SWE)derived using FY3B-MWRI with AMSR-E SWE products in the Northern Hemisphere.The results show that AMSR-E overestimated SWE in China,which agrees with other validations.     

Comparison of ensemble-based state and parameter estimation methods for soil moisture data assimilation

Model parameters are a source of uncertainty that can easily cause systematic deviation and significantly affect the accuracy of soil moisture generation in assimilation systems. This study addresses the issue of retrieving model parameters related to soil moisture via the simultaneous estimation of states and parameters based on the Common Land Model (CoLM). The state-parameter estimation algorithms AEnKF (Augmented Ensemble Kalman Filter), DEnKF (Dual Ensemble Kalman Filter) and SODA (Simultaneous optimization and data assimilation) are entirely implemented within an EnKF framework to investigate how the three algorithms can correct model parameters and improve the accuracy of soil moisture estimation. The analysis is illustrated by assimilating the surface soil moisture levels from varying observation intervals using data from Mongolian plateau sites. Furthermore, a radiation transfer model is introduced as an observation operator to analyze the influence of brightness temperature assimilation on states and parameters that are estimated at different microwave signal frequencies. Three cases were analyzed for both soil moisture and brightness temperature assimilation, focusing on the progressive incorporation of parameter uncertainty, forcing data uncertainty and model uncertainty. It has been demonstrated that EnKF is outperformed by all other methods, as it consistently maintains a bias. State-parameter estimation algorithms can provide a more accurate estimation of soil moisture than EnKF. AEnKF is the most robust method, with the lowest RMSE values for retrieving states and parameters dealing only with parameter uncertainty, but it possesses disadvantages related to increasing sources of uncertainty and decreasing numbers of observations. SODA performs well under the complex situations in which DEnKF shows slight disadvantages in terms of statistical indicators; however, the former consumes far more memory and time than the latter.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

Simultaneous estimation of soil moisture and hydraulic parameters using residual resampling particle filter

Land data assimilation (DA) has gradually developed into an important earth science research method because of its ability to combine model simulations and observations. Integrating new observations into a land surface model by the DA method can correct the predicted trajectory of the model and thus, improve the accuracy of state variables. It can also reduce uncertainties in the model by estimating some model parameters simultaneously. Among the various DA methods, the particle filter is free from the constraints of linear models and Gaussian error distributions, and can be applicable to any nonlinear and non-Gaussian state-space model; therefore, its importance in land data assimilation research has increased. In this study, a DA scheme was developed based on the residual resampling particle filter. Microwave brightness temperatures were assimilated into the macro-scale semi-distributed variance infiltration capacity model to estimate the surface soil moisture and three hydraulic parameters simultaneously. Finally, to verify the scheme, a series of comparative experiments was performed with experimental data obtained during the Soil Moisture Experiment of 2004 in Arizona. The results show that the scheme can improve the accuracy of soil moisture estimations significantly. In addition, the three hydraulic parameters were also well estimated, demonstrating the effectiveness of the DA scheme.     

The effects of scene heterogeneity on soil moisture retrieval from passive microwave data

The τ–ω model of microwave emission from soil and vegetation layers is widely used to estimate soil moisture content from passive microwave observations. Its application to prospective satellite-based observations aggregating several thousand square kilometres requires understanding of the effects of scene heterogeneity. The effects of heterogeneity in soil surface roughness, soil moisture, water area and vegetation density on the retrieval of soil moisture from simulated single- and multi-angle observing systems were tested. Uncertainty in water area proved the most serious problem for both systems, causing errors of a few percent in soil moisture retrieval. Single-angle retrieval was largely unaffected by the other factors studied here. Multiple-angle retrievals errors around one percent arose from heterogeneity in either soil roughness or soil moisture. Errors of a few percent were caused by vegetation heterogeneity. A simple extension of the model vegetation representation was shown to reduce this error substantially for scenes containing a range of vegetation types.     

Evaluation of the SMOS and SMAP soil moisture products under different vegetation types against two sparse in situ networks over arid mountainous watersheds,Northwest China

Assessment of the suitability of satellite soil moisture products at large scales is urgently needed for numerous climatic and hydrological researches, particularly in arid mountainous watersheds where soil moisture plays a key role in landatmosphere exchanges. This study presents evaluation of the SMOS(L2) and SMAP(L2_P_E and L2_P) products against ground-based observations from the Upstream of the Heihe River Watershed in situ Soil Moisture Network(UHRWSMN) and the Ecological and Hydrological Wireless Sensor Network(EHWSN) over arid high mountainous watersheds, Northwest China.Results show that all the three products are reliable in catching the temporal trend of the in situ observations at both point and watershed scales in the study area. Due to the uncertainty in brightness temperature and the underestimation of effective temperature, the SMOS L2 product and both the SMAP L2 products show "dry bias" in the high, cold mountainous area. Because of the more accurate brightness temperature observations viewing at a constant angle and more suitable estimations of single scattering albedo and optical depth, both the SMAP L2 products performed significantly better than the SMOS product.Moreover, comparing with station density of in situ network, station representation is much more important in the evaluation of the satellite soil moisture products. Based on our analysis, we propose the following suggestions for improvement of the SMOS and SMAP product suitability in the mountainous areas: further optimization of effective temperature; revision of the retrieval algorithm of the SMOS mission to reduce the topographic impacts; and, careful selection of in situ observation stations for better representation of in situ network in future evaluations. All these improvements would lead to better applicability of the SMOS and SMAP products for soil moisture estimation to the high elevation and topographically complex mountainous areas in arid regions.     

Improving the estimation of hydrothermal state variables in the active layer of frozen ground by assimilating <Emphasis Type="Italic">in situ</Emphasis> observations and SSM/I data

The active layer of frozen ground data assimilation system adopts the SHAW (Simulteneous Heat and Water) model as the model operator. It employs an ensemble kalman filter to fuse state variables predicted by the SHAW model with in situ observation and the SSM/I 19 GHz brightness temperature for the purpose of optimizing model hydrothermal state variables. When there is little water movement in the frozen soil during the winter season, the unfrozen water content depends primarily on soil temperature. Thus, soil temperature is the crucial state variable to be improved. In contrast, soil moisture is heavily influenced by precipitation during the summer season. The simulation accuracy of soil moisture has a strong and direct impact on the soil temperature. In this case, the crucial state variable to be improved is soil moisture. One-dimensional assimilation experiments that have been carried out at AMDO station show that land data assimilation method can improve the estimation of hydrothermal state variables in the soil by fusing model information and observation information. The reasonable model error covariance matrix plays a key role in transferring the optimized surface state information to the deep soil, and it provides improved estimations of whole soil state profiles. After assimilating the 4-cm soil temperature by in situ observation, the soil temperature RMSE (Root Mean Square Error) of each soil layer decreased by 0.96°C on average relative to the SHAW simulation. After assimilating the 4-cm soil moisture in situ observation, the soil moisture RMSE of each soil layer decreased by 0.020 m³·m⁻³. When assimilating the SSM/I 19 GHz brightness temperature, the soil temperature RMSE of each soil layer during the winter decreased by 0.76°C, while the soil moisture RMSE of each soil layer during the summer decreased by 0.018 m³·m⁻³.     

A mathematical model of soil moisture spatial distribution on the hill slopes of the Loess Plateau

Based on important factors that affect soil moisture spatial distribution, such as the slope gradients, land use, vegetation cover, and surface water diffusion characteristics together with field measurements of soil moisture data obtained from the surface soil under different land use structures, a soil moisture spatial distribution model was established. The diffusion degree coefficient of surface water for different vegetations was estimated from soil moisture values obtained from field measurements. The model can be solved using the finite unit method. The soil moisture spatial distribution on the hill slopes in the Loess Plateau were simulated by the model. A comparison of the simulated values with measurement data shows that the model is a good fit.     

Remote sensing in hydrology

Remote sensing provides a means of observing hydrological state variables over large areas. The ones which we will consider in this paper are land surface temperature from thermal infrared data, surface soil moisture from passive microwave data, snow cover using both visible and microwave data, water quality using visible and near-infrared data and estimating landscape surface roughness using lidar. Methods for estimating the hydrometeorlogical fluxes, evapotranspiration and snowmelt runoff, using these state variables are also described.     

Numerical Simulation of the Sensitivity of Summer Monsoon Circulation and Rainfall over India to Land Surface Processes

—The influence of soil moisture and vegetation variation on simulation of monsoon circulation and rainfall is investigated. For this purpose a simple land surface parameterization scheme is incorporated in a three-dimensional regional high resolution nested grid atmospheric model. Based on the land surface parameterization scheme, latent heat and sensible heat fluxes are explicitly estimated over the entire domain of the model. Two sensitivity studies are conducted; one with bare dry soil conditions (no latent heat flux from land surface) and the other with realistic representation of the land surface parameters such as soil moisture, vegetation cover and landuse patterns in the numerical simulation. The sensitivity of main monsoon features such as Somali jet, monsoon trough and tropical easterly jet to land surface processes are discussed.¶Results suggest the necessity of including a detailed land surface parameterization in the realistic short-range weather numerical predictions. An enhanced short-range prediction of hydrological cycle including precipitation was produced by the model, with land surface processes parameterized. This parameterization appears to simulate all the main circulation features associated with the summer monsoon in a realistic manner.     

Spatio-temporal evolution and time-stable characteristics of soil moisture within remote sensing footprints with varying soil, slope, and vegetation

Air-borne passive microwave remote sensors measure soil moisture at the footprint scale, a scale of several hundred square meters or kilometers that encompasses different characteristic combinations of soil, topography, vegetation, and climate. Studies of within-footprint variability of soil moisture are needed to determine the factors governing hydrologic processes and their relative importance, as well as to test the efficacy of remote sensors. Gridded ground-based impedance probe water content data and aircraft-mounted Electronically Scanned Thinned Array Radiometer (ESTAR) pixel-average soil moisture data were used to investigate the spatio-temporal evolution and time-stable characteristics of soil moisture in three selected (LW03, LW13, LW21) footprints from the Southern Great Plains 1997 (SGP97) Hydrology Experiment. Better time-stable features were observed within a footprint containing sandy loam soil than within two pixels containing silty loam soil. Additionally, flat topography with split wheat/grass land cover produced the largest spatio-temporal variability and the least time stability in soil moisture patterns. A comparison of ground-based and remote sensing data showed that ESTAR footprint-average soil moisture was well calibrated for the LW03 pixel with sandy loam soil, rolling topography, and pasture land cover, but improved calibration is warranted for the LW13 (silty loam soil, rolling topography, pasture land) and LW21 (silty loam soil, flat topography, split vegetation of wheat and grass land with tillage practice) pixels. Footprint-scale variability and associated nonlinear soil moisture dynamics may prove to be critical in the regional-scale hydroclimatic models.     

Evaluation of the hydrology of the IBIS land surface model in a semi‐arid catchment

This paper evaluates the Integrated BIosphere Simulator (IBIS) land surface model using daily soil moisture data over a 3‐year period (2005–2007) at a semi‐arid site in southeastern Australia, the Stanley catchment, using the Monte Carlo generalized likelihood uncertainty estimation (GLUE) approach. The model was satisfactorily calibrated for both the surface 30 cm and full profile 90 cm. However, full‐profile calibration was not as good as that for the surface, which results from some deficiencies in the evapotranspiration component in IBIS. Relatively small differences in simulated soil moisture were associated with large discrepancies in the predictions of surface runoff, drainage and evapotranspiration. We conclude that while land surface schemes may be effective at simulating heat fluxes, they may be ineffective for prediction of hydrology unless the soil moisture is accurately estimated. Sensitivity analyses indicated that the soil moisture simulations were most sensitive to soil parameters, and the wilting point was the most identifiable parameter. Significant interactions existed between three soils parameters: porosity, saturated hydraulic conductivity and Campbell ‘b’ exponent, so they could not be identified independent of each other. There were no significant differences in parameter sensitivity and interaction for different hydroclimatic years. Even though the data record contained a very dry year and another year with a very large rainfall event, this indicated that the soil model could be calibrated without the data needing to explore the extreme range of dry and wet conditions. IBIS was much less sensitive to vegetation parameters. The leaf area index (LAI) could affect the mean of daily soil moisture time series when LAI < 1, while the variance of the soil moisture time series was sensitive to LAI > 1. IBIS was insensitive to the Jackson rooting parameter, suggesting that the effect of the rooting depth distribution on predictions of hydrology was insignificant. Copyright © 2014 John Wiley & Sons, Ltd.     

Radar altimetry backscattering signatures at Ka,Ku, C,and S bands over West Africa

This study presents a comprehensive comparison of radar altimetry signatures at Ka-, Ku-, C-, and S-bands using SARAL, ENVISAT and Jason-2 data over the major bioclimatic zones, soil and vegetation types encountered in West-Africa, with an emphasis on the new information at Ka-band provided by the recently launched SARAL–Altika mission. Spatio-temporal variations of the radar altimetry responses were related to changes in surface roughness, land cover and soil wetness. Analysis of time series of backscattering coefficients along the West African bioclimatic gradient shows that radar echoes at nadir incidence are well correlated to soil moisture in semi-arid savannah environments. Radar altimeters are able to detect the presence of water even under a dense canopy cover at all frequencies. But only measurements at Ka-band are able to penetrate underneath the canopy of non-inundated tropical evergreen forests.     

Dual state-parameter estimation of root zone soil moisture by optimal parameter estimation and extended Kalman filter data assimilation

With well-determined hydraulic parameters in a hydrologic model, a traditional data assimilation method (such as the Kalman filter and its extensions) can be used to retrieve root zone soil moisture under uncertain initial state variables (e.g., initial soil moisture content) and good simulated results can be achieved. However, when the key soil hydraulic parameters are incorrect, the error is non-Gaussian, as the Kalman filter will produce a persistent bias in its predictions. In this paper, we propose a method coupling optimal parameters and extended Kalman filter data assimilation (OP-EKF) by combining optimal parameter estimation, the extended Kalman filter (EKF) assimilation method, a particle swarm optimization (PSO) algorithm, and Richards’ equation. We examine the accuracy of estimating root zone soil moisture through the optimal parameters and extended Kalman filter data assimilation method by using observed in situ data at the Meiling experimental station, China. Results indicate that merely using EKF for assimilating surface soil moisture content to obtain soil moisture content in the root zone will produce a persistent bias between simulated and observed values. Using the OP-EKF assimilation method, estimates were clearly improved. If the soil profile is heterogeneous, soil moisture retrieval is accurate in the 0-50 cm soil profile and is inaccurate at 100 cm depth. Results indicate that the method is useful for retrieving root zone soil moisture over large areas and long timescales even when available soil moisture data are limited to the surface layer, and soil moisture content are uncertain and soil hydraulic parameters are incorrect.     

基于光学与被动微波遥感的青藏高原地区土壤水分反演

青藏高原地区高精度的土壤水分反演对高原能水循环、全球大气循环研究有着极大的影响.因此,获取青藏高原土壤水分时空布信息是一个迫切需要解决的问题.温度植被干旱指数(TVDI),是基于光学与热红外遥感通道数据反演土壤水分的重要方法,但在研究区域较大、地表覆盖格局差异显著时,TVDI模型反演精度会受到地表温度(Ts)等因素的影响.被动微波AMSR-E数据精确记录了像元内的土壤水分信息,但空间分辨率低.本文利用同时期的MODIS与被动微波数据,发展了针对青藏高原地区高精度土壤水分反演算法.首先,在TVDI模型中,利用修正型土壤调整植被指数(MSAVI)代替归一化植被指数(NDVI),以改进NDVI易饱和的缺点;其次,利用ASTER GDEM数据,对地形高程和纬度差异引起的地表温度变化进行了校正;然后,通过神经网络训练建立基于TVDI、被动微波以及辅助气象数据的土壤水分反演模型,并应用该模型反演了青藏高原地区三个观测网(CAMP/Tibet、玛曲和那曲)的土壤水分;最后,利用实测土壤水分数据对反演结果进行验证,结果表明该模型的精度均方根误差(RMSE)数值可达到0.031~0.041 m~3·m~(-3).本文还应用该算法反演了青藏高原连续的土壤水分的空间分布,并比较了土壤水分的变化趋势与实测降水变化趋势,结果表明二者变化量的正负关系一致.     

Inferring land surface parameters from the diurnal variability of microwave and infrared temperatures

This study investigates the properties of the diurnal cycle of microwave brightness temperatures (TB), namely the phase and the amplitude, and their variability in time and space over the globe to infer information on key land surface parameters like changes in soil texture spatial distribution, soil moisture conditions, and vegetation density. The phase corresponds to the lag between Land Surface Temperature (LST) and TB diurnal cycles. The amplitude is determined as the difference between the maximum and the minimum of TB diurnal cycle. The diurnal cycle of TB was constructed using observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and the Special Sensor Microwave/Imager (SSM/I). The latter offer a series of sensors, namely, F13, F14, and F15 that were used in this study for a higher temporal coverage and more accurate diurnal cycle determination. LST estimates, which are available every 3 h from the International Satellite Cloud Climatology Project (ISCCP) database were used to build the LST diurnal cycle. ISCCP LST data is an infrared-based temperature with almost no penetration and is the representative of top skin temperature.The analyses of the diurnal cycles showed that the diurnal amplitude of TB decreases as the vegetation density increases, especially in the case of low frequencies which penetrate deeper into the canopy which makes them more sensitive to changes in vegetation density. The interannual variations of TB diurnal amplitudes were also in agreement with the seasonality of the vegetation cover. Over desert and rain forest regions where surface conditions do not vary significantly throughout the year, the changes in diurnal amplitudes were the lowest. A relationship between phase and amplitude values was established. It was found that the amplitude of TB diurnal cycle decreases when the phase lag increases. The spatial distribution of the determined diurnal properties, namely, phase and amplitude of TB, showed an agreement with lithology maps in desert areas. Lower TB amplitudes were observed over regions with loose siliceous rocks. Phase lag values between 1.5 and 3 h corresponded to 83% of the class “loose siliceous rocks” in the Sahara Desert, which corroborates the potential of using the diurnal properties of TB as an indicator of land surface parameters.