Surface layer SM simulations from the WRF model for a period of four months (July 31 to December 01, 2010) were used in this paper. We used the Noah land surface model (LSM) as a land parameterization scheme.20 Noah LSM is one of the multiple land surface schemes coupled with the WRF modeling system.21 The initial conditions for SM simulations in the model are defined by the National Center for Environmental Prediction (NCEP) global final analysis (FNL) data with a spatial resolution of ( in our study area) and 6-h frequency.12 The input variables from the NCEP FNL data include (but are not limited to) soil moisture/water content, soil temperature, precipitation, heat flux, humidity, surface winds, and land cover (see http://rda.ucar.edu/datasets/ds083.2). The Noah LSM only provides surface heat and moisture fluxes as lower boundary conditions to the coupled atmospheric (i.e., the WRF) model.22 These fluxes are then transported throughout the boundary layer, and interact with other model physics involving cloud, radiation, and precipitation processes.22 In the Noah LSM, the green vegetation fraction (Fg) parameter, which is defined as the coverage of vegetation over an area unit or a pixel,12 is used to differentiate moisture flux from vegetation. As a result, the WRF model coupled with the Noah LSM is able to separate soil and vegetation moisture fluxes and also accounts for precipitation effects in SM approximation. To calculate SM, a prognostic diffusion equation for the volumetric SM content is used.23 Sensitivity of the WRF model to different land surface parameterization schemes, such as the Noah LSM, is discussed in the literature,12,21 which includes SM simulations. Noah LSM simulates the near-surface SM in four predefined layers (with a top-down sequence of 0.05, 0.25, 0.70, and 1.5 m depth). Outputs for the top most layer (5-cm depth) were used in our analysis. The first 24 h of the simulations were discarded from the analysis as the spin-up period required by the model to reach a balanced state with the boundary conditions. Grid spacing of the simulations was set to 5 km for the first and 1 km for the second domain of the model to match the spatial resolution of the MODIS LST. The time interval of the model outputs was set to a 30-min rate for all simulations.