Net ecosystem productivity (NEP) is an important index for the quantitative evaluation of carbon sources and sinks in terrestrial ecosystems. Based on MOD17A3 and meteorological data, the vegetation NEP was estimated from 2000 to 2021 in the Loess Plateau (LP) and its six ecological subregions of the LP (loess sorghum gully subregions:A1, A2; loess hilly and gully subregions:B1, B2; sandy land and agricultural irrigation subregion:C; and earth-rock mountain and river valley plain subregion:D). Combined with the terrain, remote sensing, and human activity data, Theil-Sen Median trend analysis, correlation analysis, multiple regression residual analysis, and geographic detector were used, respectively, to explore the spatio-temporal characteristics of NEP and its response mechanism to climate, terrain, and human activity. The results showed that:① On the temporal scale, from 2000 to 2021 the annual mean NEP of the LP region (in terms of C) was 104.62 g·(m2·a)-1. The annual mean NEP for both the whole LP and each of the ecological subregions showed a significant increase trend, and the NEP of the LP increased by 6.10 g·(m2·a)-1 during the study period. The highest growth rate of the NEP was 9.04 g·(m2·a)-1, occurring in the A2 subregion of the loess sorghum gully subregions. The subregion C had the lowest growth rate of 2.74 g·(m2·a)-1. Except for the C subregion, all other ecological subregions (A1, A2, B1, B2, and D) were carbon sinks. ② On the spatial scale, the spatial distribution of annual NEP on the LP was significantly different, with the higher NEP distribution in the southeast of the LP and the lower in the northwest of the LP. The high carbon sink area was mainly distributed in the southern part of the loess sorghum gully subregions, and the carbon source area was mainly distributed in the northern part of the loess sorghum gully subregions and most of the C subregion. The high growth rate was mainly distributed in the central and the southern part of the A2 subregion and the southwest part of the B2 subregion. ③ Human activities had the greatest influence on the temporal variation in NEP in the LP and all the ecological subregions, with the correlation coefficient between human activity data and NEP being above 0.80, and the relative contribution rates of human factors was greater than 50%. The spatial distribution was greatly affected by meteorological factors, among which the precipitation and solar radiation were the main factors affecting the spatial changes in the NEP of the LP. The temporal and spatial variations in the NEP in the LP were influenced by natural and human social factors. To some extent, these results can provide a reference for the terrestrial ecosystem in the LP to reduce emissions and increase sinks and to achieve the goal of double carbon.
Keywords: Loess Plateau（LP）; correlation analysis; ecological subregion; geographic detector; net ecosystem productivity; residual analysis.