Forest carbon flux is critical to climate change, and the accurate modeling of forest carbon flux is an extremely challenging task. The remote sensing model (the MODIS MOD_17 gross primary productivity (GPP) model (MOD_17)) has strong practicability and is widely used around the world. The ecological process (the Biome-BioGeochemical Cycles Multilayer Soil Module model (Biome-BGCMuSo)) model can describe most of the vegetation’s environmental and physiological processes on fine time scales. Nevertheless, complex parameters and calibrations pose challenges to the application and development of models. In this study, we optimized all the input parameters of the MOD_17 model for the calibration of the Biome-BGCMuSo model to obtain GPP with continuous spatiality. To determine the contribution of input parameters to the GPP of different forest types, an Extended Fourier Amplitude Sensitivity Test (EFAST) was performed on the Biome-BGCMuSo model firstly. Then, we selected the sample points of each forest type and its different ecological gradients (30 for each type), using the GPP simulation value of the optimized MOD_17 model corresponding to the time and space scale to calibrate the Biome-BGCMuSo model, to drive the calibrated Biome-BGCMuSo, and we simulated the different forest types’ net primary productivity (NPP). According to dendrochronological measurements, the NPP simulation results were verified on the whole regional scale. The results showed that the GPP values of different forest types were highly sensitive to C:N_leaf (C:N of leaf), SLA1 (canopy average specific leaf area in phenological phase 1), and FLNR (fraction of leaf N in Rubisco). The coefficient of determination (R²) between the simulated forest NPP and the measured NPP was 0.64, and the root-mean-square (RMSE) was 26.55 g/C/m²/year. Our study aims to reduce uncertainty in forest carbon fluxes simulated by the Biome-BGCMuSo model, providing feedback for understanding forest ecosystem carbon cycling, vegetation productivity, and climate change.

The leaf area index (LAI) is a crucial parameter for analyzing terrestrial ecosystem carbon cycles and global climate change. Obtaining high spatiotemporal resolution forest stand vegetation LAI products over large areas is essential for an accurate understanding of forest ecosystems. This study takes the northwestern part of the Inner Mongolia Autonomous Region (the northern section of the Greater Khingan Mountains) in northern China as the research area. It also generates the LAI time series product of the 8-day and 30 m forest stand vegetation growth period from 2013 to 2017 (from the 121st to the 305th day of each year). The Simulated Annealing-Back Propagation Neural Network (SA-BPNN) model was used to estimate LAI from Landsat8 OLI, and the multi-period GaoFen-1 WideField-View satellite images (GF-1 WFV) and the spatiotemporal adaptive reflectance fusion mode (STARFM) was used to predict high spatiotemporal resolution LAI by combining inversion LAI and Global LAnd Surface Satellite-derived vegetation LAI (GLASS LAI) products. The results showed the following: (1) The SA-BPNN estimation model has relatively high accuracy, with R² = 0.75 and RMSE = 0.38 for the 2013 LAI estimation model, and R² = 0.74 and RMSE = 0.17 for the 2016 LAI estimation model. (2) The fused 30 m LAI product has a good correlation with the LAI verification of the measured sample site (R² = 0.8775) and a high similarity with the GLASS LAI product. (3) The fused 30 m LAI product has a high similarity with the GLASS LAI product, and compared with the GLASS LAI interannual trend line, it accords with the growth trend of plants in the seasons. This study provides a theoretical and technical reference for forest stand vegetation growth period LAI spatiotemporal fusion research based on high-score data, and has an important role in exploring vegetation primary productivity and carbon cycle changes in the future.