This study evaluates the land-atmosphere interaction using various land surface models (LSMs). Selected schemes are the Oregon State University (OSU) LSM, Noah LSM, and Variable Infiltration Capacity (VIC) LSM. To investigate the characteristics of the land surface process, two types of comparison methods are employed in this study: 1) preliminary evaluation using a single-column model (SCM), 2) examination of the simulated results using a global climate model (GCM). Both models are a version of National Centers for Environmental Prediction (NCEP) Medium-Range Forecast (MRF) model. The SCM is forced by the Atmospheric Radiation Measurement (ARM) over the Southern Great Plains data and Coordinated Energy and Water Cycle Observations Project (CEOP) reference sites observations. The results with the GCM based on the NCEP global spectral model with the Land Data Assimilation System (LDAS) boundary forcing are compared to those from the SCM data.