The parameters of the GR4J-CemaNeige coupling model (GR4neige) are typically treated as constants. However, the maximum capacity of the production store (parX1) exhibits time-varying characteristics due to climate variability and vegetation coverage change. This study employed differentiable parameter learning (dPL) to identify the time-varying parX1 in the GR4neige across 671 catchments within the United States. We built two types of dPL, including static and dynamic parameter networks, to assess the advantages of the time-varying parameter. In the dynamic parameter network, we evaluated the impact of potential evapotranspiration (PET), precipitation (P), temperature (T), soil moisture (SM), and normalized difference vegetation index (NDVI) datasets on the performance of dPL. We then compared dPL with the empirical functional method (fm). The results demonstrated that the dynamic parameter network outperformed the static parameter network in streamflow estimation. There were differences in streamflow estimation among the dynamic parameter network driven by various input features. In humid catchments, simultaneously incorporating all five factors, including PET, P, T, SM, and the NDVI, achieved optimal streamflow simulation accuracy. In arid catchments, it was preferable to introduce PET, T, and the NDVI separately for improved performance. dPL significantly outperformed the empirical fm in estimating streamflow and uncalibrated intermediate variables, like evapotranspiration (ET). Both the derived parX1 from dPL and the empirical fm exhibited significant spatiotemporal variation across 671 catchments. Notably, compared to parX1 obtained through the empirical fm, parX1 derived from dPL exhibited a distinct spatial clustering pattern. This study highlights the potential of dPL in enhancing model accuracy and contributes to understanding the spatiotemporal variation characteristics of parX1 under the influence of climate factors, soil conditions, and vegetation change.

The demand for freshwater is increasing with population growth and rapid socio-economic development. It is more and more important for refined irrigation water management to conduct research on crop evapotranspiration (ET) data with a high spatiotemporal resolution in agricultural regions. We propose the unmixing–weight ET image fusion model (UWET), which integrates the advantages of the unmixing method in spatial downscaling and the weight-based method in temporal prediction to produce daily ET maps with a high spatial resolution. The Landsat-ET and MODIS-ET datasets for the UWET fusion data are retrieved from Landsat and MODIS images based on the surface energy balance model. The UWET model considers the effects of crop phenology, precipitation, and land cover in the process of the ET image fusion. The precision evaluation is conducted on the UWET results, and the measured ET values are monitored by eddy covariance at the Luancheng station, with average MAE values of 0.57 mm/day. The image results of UWET show fine spatial details and capture the dynamic ET changes. The seasonal ET values of winter wheat from the ET map mainly range from 350 to 660 mm in 2019–2020 and from 300 to 620 mm in 2020–2021. The average seasonal ET in 2019–2020 is 499.89 mm, and in 2020–2021, it is 459.44 mm. The performance of UWET is compared with two other fusion models: the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Spatial and Temporal Reflectance Unmixing Model (STRUM). UWET performs better in the spatial details than the STARFM and is better in the temporal characteristics than the STRUM. The results indicate that UWET is suitable for generating ET products with a high spatial–temporal resolution in agricultural regions.

Accurate estimation of gross primary production (GPP) of paddy rice fields is essential for understanding cropland carbon cycles, yet remains challenging due to spatial heterogeneity. In this study, we integrated high-resolution unmanned aerial vehicle (UAV) imagery into a leaf biochemical properties-based model for improving GPP estimation. The key parameter, maximum carboxylation rate at the top of the canopy ($V_{cmax,0}^{25}$), was quantified using various spatial information representation methods, including mean (μ_ref) and standard deviation (σ_ref) of reflectance, gray-level co-occurrence matrix (GLCM)-based features, local binary pattern histogram (LBPH), and convolutional neural networks (CNNs). Our models were evaluated using a two-year eddy covariance (EC) system and UAV measurements. The result shows that incorporating spatial information can vastly improve the accuracy of $V_{cmax,0}^{25}$ and GPP estimation. CNN methods achieved the best $V_{cmax,0}^{25}$ estimation, with an R of 0.94, an RMSE of 19.44 μmol m⁻² s⁻¹, and an MdAPE of 11%, and further produced highly accurate GPP estimates, with an R of 0.92, an RMSE of 6.5 μmol m⁻² s⁻¹, and an MdAPE of 23%. The μ_ref-GLCM texture features and μ_ref-LBPH joint-driven models also gave promising results. However, σ_ref contributed less to $V_{cmax,0}^{25}$ estimation. The Shapley value analysis revealed that the contribution of input features varied considerably across different models. The CNN model focused on nir and red-edge bands and paid much attention to the subregion with high spatial heterogeneity. The μ_ref-LBPH joint-driven model mainly prioritized reflectance information. The μ_ref-GLCM-based features joint-driven model emphasized the role of GLCM texture indices. As the first study to leverage the spatial information from high-resolution UAV imagery for GPP estimation, our work underscores the critical role of spatial information and provides new insight into monitoring the carbon cycle.