Abstract. The processes involved in the exchange of water, energy and carbon in terrestrial ecosystems are strongly intertwined. To accurately represent the terrestrial biosphere in land surface models (LSM), the intrinsic coupling between these processes is required. Soil moisture and leaf area index are two key variables at the nexus of water, energy and vegetation. Here, we evaluated three LSM (ISBA, ORCHIDEE and a diagnostic model, based on the LSA SAF algorithms) in their ability to simulate the latent heat flux (LE) and gross primary production (GPP) coherently, and their interactions through leaf area index (LAI) and soil moisture. The models were validated using in situ eddy covariance observations, soil moisture measurements and remote sensed LAI. It was found that the diagnostic model performed consistently well, regardless land cover, whereas important shortcomings of the prognostic models were revealed for in herbaceous/dry sites. Despite their different architecture and parametrization, ISBA and ORCHIDEE shared some key weaknesses. In both models, LE and GPP were found to be oversensitive to drought stress. Though the simulated soil water dynamics could be improved, this was not the main cause of errors in the surface fluxes. Instead, these errors were strongly correlated to errors in LAI. The simulated phenological cycle in ISBA and ORCHIDEE was delayed compared to observations, and failed to capture the observed seasonal variability. The feedback mechanism between GPP and LAI (i.e. the biomass allocation scheme) was identified as a key element to improve the intricate coupling between energy, water and vegetation in LSM.