This study improved simulation of forest carbon fluxes in the Changbai Mountains with a process-based model (Biome-BGC) using incorporation and data assimilation. Firstly, the original remote sensing-based MODIS MOD_17 GPP (MOD_17) model was optimized using refined input data and biome-specific parameters. The key ecophysiological parameters of the Biome-BGC model were determined through the Extended Fourier Amplitude Sensitivity Test (EFAST) sensitivity analysis. Then the optimized MOD_17 model was used to calibrate the Biome-BGC model by adjusting the sensitive ecophysiological parameters. Once the best match was found for the 10 selected forest plots for the 8-day GPP estimates from the optimized MOD_17 and from the Biome-BGC, the values of sensitive ecophysiological parameters were determined. The calibrated Biome-BGC model agreed better with the eddy covariance (EC) measurements (R² = 0.87, RMSE = 1.583 gC·m⁻²·d⁻¹) than the original model did (R² = 0.72, RMSE = 2.419 gC·m⁻²·d⁻¹). To provide a best estimate of the true state of the model, the Ensemble Kalman Filter (EnKF) was used to assimilate five years (of eight-day periods between 2003 and 2007) of Global LAnd Surface Satellite (GLASS) LAI products into the calibrated Biome-BGC model. The results indicated that LAI simulated through the assimilated Biome-BGC agreed well with GLASS LAI. GPP performances obtained from the assimilated Biome-BGC were further improved and verified by EC measurements at the Changbai Mountains forest flux site (R² = 0.92, RMSE = 1.261 gC·m⁻²·d⁻¹).

We propose a long-term parameterization scheme for two critical parameters, zero-plane displacement height (d) and aerodynamic roughness length (z_0m), that we further use in the Surface Energy Balance System (SEBS). A sensitivity analysis of SEBS indicated that these two parameters largely impact the estimated sensible heat and latent heat fluxes. First, we calibrated regression relationships between measured forest vertical parameters (Lorey’s height and the frontal area index (FAI)) and forest aboveground biomass (AGB). Next, we derived the interannual Lorey’s height and FAI values from our calibrated regression models and corresponding forest AGB dynamics that were converted from interannual carbon fluxes, as simulated from two incorporated ecological models and a 2009 forest basis map These dynamic forest vertical parameters, combined with refined eight-day Global LAnd Surface Satellite (GLASS) LAI products, were applied to estimate the eight-day d, z_0m, and, thus, the heat roughness length (z_0h). The obtained d, z_0m and z_0h were then used as forcing for the SEBS model in order to simulate long-term forest evapotranspiration (ET) from 2000 to 2012 within the Qilian Mountains (QMs). As compared with MODIS, MOD16 products at the eddy covariance (EC) site, ET estimates from the SEBS agreed much better with EC measurements (R² = 0.80 and RMSE = 0.21 mm·day⁻¹).