Spatial non-stationary characteristics between grass yield and its influencing factors in the Ningxia temperate grasslands based on a mixed geographically weighted regression model

Spatial models are effective in obtaining local details on grassland biomass, and their accuracy has important practical significance for the stable management of grasses and livestock. To this end, the present study utilized measured quadrat data of grass yield across different regions in the main growing season of temperate grasslands in Ningxia of China (August 2020), combined with hydrometeorology, elevation, net primary productivity (NPP), and other auxiliary data over the same period. Accordingly, non-stationary characteristics of the spatial scale, and the effects of influencing factors on grass yield were analyzed using a mixed geographically weighted regression (MGWR) model. The results showed that the model was suitable for correlation analysis. The spatial scale of ratio resident-area index (PRI) was the largest, followed by the digital elevation model, NPP, distance from gully, distance from river, average July rainfall, and daily temperature range; whereas the spatial scales of night light, distance from roads, and relative humidity (RH) were the most limited. All influencing factors maintained positive and negative effects on grass yield, save for the strictly negative effect of RH. The regression results revealed a multiscale differential spatial response regularity of different influencing factors on grass yield. Regression parameters revealed that the results of Ordinary least squares (OLS) (Adjusted R² = 0.642) and geographically weighted regression (GWR) (Adjusted R² = 0.797) models were worse than those of MGWR (Adjusted R² = 0.889) models. Based on the results of the RMSE and radius index, the simulation effect also was MGWR > GWR > OLS models. Ultimately, the MGWR model held the strongest prediction performance (R² = 0.8306). Spatially, the grass yield was high in the south and west, and low in the north and east of the study area. The results of this study provide a new technical support for rapid and accurate estimation of grassland yield to dynamically adjust grazing decision in the semi-arid loess hilly region.