Leaf area index (LAI) is often used to quantify plant production and evapotranspiration with terrestrial ecosystem models (TEMs). This study evaluated theLAIsimulation in North America using a data assimilation technique and a process‐based TEM as well as in situ and satellite data. We first optimized the parameters related toLAIin theTEMusing a Markov Chain Monte Carlo method, and AmeriFlux site‐level and regionalLAIdata from advanced very high‐resolution radiometer. The parameterized model was then verified with the observed monthlyLAIof major ecosystem types at site level. SimulatedLAIwas compared well with the observed data at sites of Harvard Forest (R2 = 0.96), University of Michigan Biological Station (R2 = 0.87), Howland Forest (R2 = 0.96), Morgan Monroe State Forest (R2 = 0.85), Shidler Tallgrass Prairie (R2 = 0.82), and Donaldson (R2 = 0.75). The root‐mean‐square error (RMSE) between modeled and satellite‐based monthlyLAIin North America is 1.4 m2/m2for the period of 1985–2010. The simulated average monthlyLAIin recent three decades increased by (3 ± 0.5)% in the region, with 1.24, 1.46, and 2.21 m2/m2on average, in Alaska, Canada, and the conterminous United States, respectively, which is consistent with satellite data. The model performed well for wet tundra, boreal forest, temperate coniferous forests, temperate deciduous forests, grasslands, and xeric shrublands (RMSE < 1.5 m2/m2), but not for alpine tundra and xeric woodlands (RMSE > 1.5 m2/m2). Both the spring and fallLAIin the 2000s are higher than that in the 1980s in the region, suggesting that the leaf phenology has an earlier onset and later senescence in the 2000s. The averageLAIincreased in April and September by 0.03 and 0.24 m2/m2, respectively. This study provides a way to quantifyLAIwith ecosystem models, which will improve future carbon and water cycling studies.