North American Land Data Assimilation System (NLDAS) land surface models have been run for a retrospective period forced by atmospheric observations from the Eta analysis and actual precipitation and downward solar radiation to calculate land hydrology. We evaluated these simulations using in situ observations over the southern Great Plains for the periods of May--September of 1998 and 1999 by comparing the model outputs with surface latent, sensible, and ground heat fluxes at 24 Atmospheric Radiation Measurement/Cloud and Radiation Testbed stations and with soil temperature and soil moisture observations at 72 Oklahoma Mesonet stations. The standard NLDAS models do a fairly good job but with differences in the surface energy partition and in soil moisture between models and observations and among models during the summer, while they agree quite well on the soil temperature simulations. To investigate why, we performed a series of experiments accounting for differences between model-specified soil types and vegetation and those observed at the stations, and differences in model treatment of different soil types, vegetation properties, canopy resistance, soil column depth, rooting depth, root density, snow-free albedo, infiltration, aerodynamic resistance, and soil thermal diffusivity. The diagnosis and model enhancements demonstrate how the models can be improved so that they can be used in actual data assimilation mode.