In the sense of 17-year climatology, the most significant differences between LAI/FASIR and LAI/GIMMS were shown in tropical regions of Amazon basin and central Africa not only in the magnitude but also the interannual variability. In boreal summer season, there also exist large differences in mid- and high latitudinal North American and Eurasian continents. The seasonal and interannual variations of the data sets were examined using EOF analysis for different regions to find anomalous variations. The singular value decomposition (SVD) analysis was applied to examine the precipitation-LAI coupled variability. The most-dominant coupled mode explains about 30% of the total variability and correlation coefficients between principle components of precipitation and LAI are greater than 0.6 over Africa and South America. But, it was hard to recognize significant differences of coupled variability between LAI/FASIR and LAI/GIMMS datasets. Although it is not easy to identify merits and limitations of the currently available LAI dataset due to the lack of ground measurements, it is necessary for the climate simulation to take account for interannual variability of vegetation intensity, which can be provided by only satellite-based LAI datasets. In order to investigate feasibility to use these datasets as a surface boundary forcing for climate simulation, we are conducting 17-year climate simulation with the NCEP GSM coupled with the SSiB. In the preliminary results, two LAI datasets induced significant differences in the simulated climate especially for the monsoon evolution over some areas, such as East Asia. Further detail analysis of LAI datasets and simulation results will be discussed in the presentation.