Investigating net primary production in climate regions of Khuzestan Province, Iran using CASA model

This study aimed to investigate the vegetation production changes in Khuzestan province, Iran using MODIS data production, meteorological data, vegetation maps as well as topographic and field monitoring data in CASA model. The study area was divided into different climatic classes based on multivariate statistical method, so the vegetation of each climatic region was examined separately for changes in NPP values. Production changes due to degradation were calculated using the Miami model and subsequently, the rain use efficiency (RUE) and the light use efficiency (LUE) and correlation indices between the CASA model and ground data were determined. The results of this study (R2) showed that the accuracy of this model varies depending on the type of climatic regions (R2 = 0.80 to R2 = 0.15). In different climatic regions, the rate of NPP changes (very humid 68 gC/m2 to ultra-dry 15 gC/m2) varies in rangeland types. The highest rate of vegetation production is observed seasonally in May. Degradation conditions also reduced RUE and LUE. However, in hyper-arid regions, adaptations of plants in some different species (Hammada Spp.) increase their efficiency compared to other vegetation types. The results showed the importance of vegetation and climate classification in vegetation production studies.

A global dataset of the fraction of absorbed photosynthetically active radiation for 1982-2022

The fraction of absorbed photosynthetically active radiation (FPAR) is an essential biophysical parameter that characterizes the structure and function of terrestrial ecosystems. Despite the extensive utilization of several satellite-derived FPAR products, notable temporal inconsistencies within each product have been underscored. Here, the new generation of the GIMMS FPAR product, GIMMS FPAR4g, was developed using a combination of a machine learning algorithm and a pixel-wise multi-sensor records integration approach. PKU GIMMS NDVI, which eliminates the orbital drift and sensor degradation issues, was used as the data source. Comparisons with ground-based measurements indicate root mean square errors ranging from 0.10 to 0.14 with R-squared ranging from 0.73 to 0.87. More importantly, our product demonstrates remarkable spatiotemporal coherence and continuity, revealing a persistent terrestrial darkening over the past four decades (0.0004 yr-1, p < 0.001). The GIMMS FPAR4g, available for half-month intervals at a spatial resolution of 1/12° from 1982 to 2022, promises to be a valuable asset for in-depth analyses of vegetation structures and functions spanning the last 40 years.

Analysis of change process of NPP dominated by human activities in Northwest Hubei, China, from 2000 to 2020

Clarifying the spatial distribution of the impact of different human disturbance activities on the net primary productivity (NPP) in regions with single climatic conditions is of considerable importance to ecological protection. Time-series NPP from 2000 to 2020 was simulated in Northwest Hubei, China, and the effects of the climate and human activities on the NPP changes were separated. Research results showed that from 2000 to 2020, the NPP change with an area of 10,166.63 km2 in Northwest Hubei is influenced by climate and human activities. Among them, human activities account for as high as 84.53%. From 2000 to 2020, the NPP in Northwest Hubei showed a slight upward trend at a rate of 1.61 g C m-2 year-1. The significantly increased NPP accounted for 21.4% of the total, which was mainly distributed in north of Northwest Hubei. And the farming of cultivated land led to the increase of NPP in west as well as the reduced human distribution in cultivated land, which was scattered in forests. Only 6.67% of the total area demonstrated a significantly decreased NPP, which was distributed mainly in the central affected by the expansion of rural-urban land and change of broad-leaved forests to shrubs and in southeast regions of Northwest Hubei caused by the increase in potential evapotranspiration. This study refined the driving factors of spatial heterogeneity of NPP changes in Northwest Hubei, which is conducive to rational planning of terrestrial ecosystem protection measures.

Change and relationship between growing season metrics and net primary productivity in forestland and grassland in China

BACKGROUND

Vegetation phenology can characterize ecosystem functions and plays a key role in the dynamics of plant productivity. Here we investigated the changes in growing season metrics (start of growing season, SOS; end of growing season, EOS; length of growing season, LOS) and their relationships with net primary productivity (NPP) in forestland and grassland in China during 1981-2016.

RESULTS

SOS advanced, EOS delayed, LOS prolonged and NPP increased significantly in 23.7%, 21.0%, 40.5% and 19.9% of the study areas, with an average rate of 3.9 days decade-1, 3.3 days·decade-1, 6.7 days·decade-1 and 10.7 gC m-2·decade-1, respectively. The changes in growing season metrics were obvious in Northwest China (NWC) and North China (NC), but the least in Northeast China (NEC). NPP was negatively correlated with SOS and positively correlated with EOS and LOS in 22.0%, 16.3% and 22.8% of the study areas, respectively, and the correlation between NPP and growing season metrics was strong in NWC, NC and Southwest China (SWC), but weak in NEC and South China (SC).

CONCLUSION

The advanced SOS, delayed EOS and prolonged LOS all contribute to the increased NPP in forestland and grassland in China, especially in NWC, NC and SWC. This study also highlights the need to further study the response of NPP to growing season changes in different regions and under the influence of multiple factors.

Practical Guide to Measuring Wetland Carbon Pools and Fluxes

Wetlands cover a small portion of the world, but have disproportionate influence on global carbon (C) sequestration, carbon dioxide and methane emissions, and aquatic C fluxes. However, the underlying biogeochemical processes that affect wetland C pools and fluxes are complex and dynamic, making measurements of wetland C challenging. Over decades of research, many observational, experimental, and analytical approaches have been developed to understand and quantify pools and fluxes of wetland C. Sampling approaches range in their representation of wetland C from short to long timeframes and local to landscape spatial scales. This review summarizes common and cutting-edge methodological approaches for quantifying wetland C pools and fluxes. We first define each of the major C pools and fluxes and provide rationale for their importance to wetland C dynamics. For each approach, we clarify what component of wetland C is measured and its spatial and temporal representativeness and constraints. We describe practical considerations for each approach, such as where and when an approach is typically used, who can conduct the measurements (expertise, training requirements), and how approaches are conducted, including considerations on equipment complexity and costs. Finally, we review key covariates and ancillary measurements that enhance the interpretation of findings and facilitate model development. The protocols that we describe to measure soil, water, vegetation, and gases are also relevant for related disciplines such as ecology. Improved quality and consistency of data collection and reporting across studies will help reduce global uncertainties and develop management strategies to use wetlands as nature-based climate solutions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13157-023-01722-2.

Terrestrial CO2 exchange diagnosis using a peatland-optimized vegetation photosynthesis and respiration model (VPRM) for the Hudson Bay Lowlands

Satellite-based light use efficiency (LUE) models have been widely used to estimate gross primary production in various terrestrial ecosystems such as forests and croplands, but northern peatlands have received less attention. In particular, the Hudson Bay Lowlands (HBL) which is a massive peatland-rich region in Canada has been largely ignored in previous LUE-based studies. These peatland ecosystems have accumulated large stocks of organic carbon over many millennia, and play a vital role in the global carbon cycle. In this study, we used the satellite data-driven Vegetation Photosynthesis and Respiration Model (VPRM) to examine the suitability of LUE models for carbon flux diagnosis in the HBL. VPRM was driven alternately with the satellite-derived enhanced vegetation index (EVI) and solar-induced chlorophyll fluorescence (SIF). The model parameter values were constrained by eddy covariance (EC) tower observations from the Churchill fen and Attawapiskat River bog sites. The main objectives of the study were to (i) investigate if site-specific parameter optimization improved NEE estimates, (ii) determine which satellite-based proxy of photosynthesis produced more reliable estimates of peatland net carbon exchange, and (iii) examine how LUE and other model parameters vary within and between the study sites. The results indicate that the VPRM mean diurnal and monthly estimates of NEE had significant strong agreements with EC tower fluxes at the two study sites. A comparison of the site-optimized VPRM against a generic peatland-optimized version of the model revealed that the site-optimized VPRM provided better estimates of NEE only during the calibration period at the Churchill fen. The diurnal and seasonal cycles of peatland carbon exchange were better captured by the SIF-driven VPRM, demonstrating that SIF is a more accurate proxy for photosynthesis compared to EVI. Our study suggests that satellite-based LUE models have the potential to be applied on a larger scale to the HBL region.

Microplastics and hydrocarbons in soils: Quantification as an anthropic carbon source

The literature on the presence of microplastics (MPs) and their potential impact on terrestrial ecosystems is still scarce. Interestingly, soil MPs are detected as organic carbon (SOC) using traditional quantification methods (e.g., loss on ignition [LOI]), although its dynamics in the environment will be different. The objective of this study was to quantify the carbon (C) contribution of MPs to the SOC in superficial soil samples from a coastal urban wetland (Avellaneda, Buenos Aires, Argentina) with the features of a humid subtropical forest and compare with hydrocarbon contribution. Soil samples were split for analysis of moisture content; texture (sieve and pipet method); organic matter as a LOI (8 h at 450 °C); total hydrocarbons (THCs; gravimetry of solvent extractable matter); n-alkanes (solvent extraction and gas chromatography-flame ionization detection analysis); and extraction of MPs (floatation in NaCl_aq , filtration, H₂ O₂ digestion, and visual sorting under a stereomicroscope). The superficial soil was a sandy clay loam with a large organic matter content (19%-30%). The THC averaged 2.5 ± 1.9 g kg and the marked predominance of odd-numbered carbon n-alkanes maximizing at C29 and C31 show the contribution of the terrestrial plant waxes. The average number of MPs was 587 ± 277 items kg of dry soil, predominantly fibers. Taking account of the C content, THCs and MPs add to the soil 1.23 ± 1.10 ton C ha and 0.10-0.97 ton C ha, respectively. Therefore, in this system with humid forest characteristics, the MPs represent between 0.12% and 1.25% of soil estimated carbon, in a magnitude similar to the C contribution of THCs (0.6%-4.2%). This preliminary study shows the relevance of discriminating MPs from other carbon sources and presents a description of their impact on soils to advance future research or tools for decision-makers. Integr Environ Assess Manag 2023;19:698-705. © 2022 SETAC.

No Signs of Long-term Greening Trend in Western Mongolian Grasslands

Trends for increased vegetation greenness based on satellite-derived data have been repeatedly published for the temperate grassland biome (including forest steppes) of eastern Inner Asia since 1982. Although this greening trend has been attenuated or partially reversed by drought in the early twenty-first century, linear increases in the Normalized Difference Vegetation Index (NDVI) or other parameters of vegetation greenness are nevertheless evident when the period since 1982 is regarded. However, the question arises whether these trends are part of a long-term trend driven by climate change, as simultaneously forests in the region show widespread drought-induced growth reductions and mortality outbreaks. Therefore, we hypothesized that the post-1982 greening trend was neither part of a long-term trend nor unprecedented. To test this hypothesis, we analyzed monthly maximum NDVI data from AVHRR time series and correlated these data with standardized tree-ring data of Larix sibirica from two regions of western Mongolia. We used linear regression to model the NDVI from tree-ring anomalies and to reconstruct the NDVI since 1940. These reconstructions show that the availability of satellite-based NDVI data coincidentally began during a dry period of low vegetation greenness in the early 1980s and was followed by a wet phase in the 1990s, producing the linear greening trend. No positive long-term trend in the reconstructed NDVI was observed from 1940 to 2010. This result rules out a recent climate change-driven greening trend for the grasslands and forest steppes of western Mongolia and calls into question its existence for all of eastern Inner Asia.

Spatiotemporal Variability in Water-Use Efficiency in Tianshan Mountains (Xinjiang, China) and the Influencing Factors

Water-use efficiency (WUE) is a crucial physiological index in carbon–water interactions and is defined as the ratio of vegetation productivity to water loss. The variation in climatic variables and drought have the most significant effects on WUE and evapotranspiration (ET). Nevertheless, how WUE varies with climate factors and drought processes in the Tianshan Mountains (TMS) is still poorly understood. In the present work, we analyzed the spatiotemporal variations in WUE, and investigated the correlations between WUE, climate factors, and drought, in the study area. The results showed that, in the TMS during 2000–2020, annual net primary productivity (NPP) ranged from 147.9 to 189.4 gC·m−2, annual ET was in the range of 212.5–285.8 mm, and annual WUE ranged from 0.66 to 0.78 gC·kg−1·H2O. Both NPP and ET exhibited an increasing trend with some fluctuation, whereas WUE showed the opposite tendency during the study period. The obtained results demonstrated that the decrease in WUE was primarily because of the increase in ET. There were obvious differences in WUE, under different land-use types, caused by NPP and ET. However, the interannual variation in WUE showed small fluctuations and the dynamic process of WUE in each land-use type showed good consistency. Temperature and wind speed had a positive influence on WUE in the middle and eastern regions of the TMS. Precipitation also played a mainly positive role in enhancing WUE, especially on the northern slope of the TMS. There was strong spatial heterogeneity of the correlation coefficient (0.68, p < 0.05) between WUE and the temperature vegetation drought index (TVDI). Moreover, the slopes of WUE and TVDI showed good consistency in terms of spatial distribution, suggesting that drought had a significant impact on ecosystem WUE. This work will enhance the understanding of WUE variation, and provide scientific evidence for water resource management and sustainable utilization in the study area.

The Different Impacts of Climate Variability and Human Activities on NPP in the Guangdong–Hong Kong–Macao Greater Bay Area

As two main drivers of vegetation dynamics, climate variability and human activities greatly influence net primary productivity (NPP) variability by altering the hydrothermal conditions and biogeochemical cycles. Therefore, studying NPP variability and its drivers is crucial to understanding the patterns and mechanisms that sustain regional ecosystem structures and functions under ongoing climate variability and human activities. In this study, three indexes, namely the potential NPP (NPPp), actual NPP (NPPa), and human-induced NPP (NPPh), and their variability from 2000 to 2020 in the Guangdong–Hong Kong–Macao Greater Bay Area (GBA) were estimated and analyzed. Six main scenarios were generated based on change trends in the three indexes over the past 21 years, and the different relative impacts of climate variability and human activities on NPPa variability were quantitatively analyzed and identified. The results showed that the NPPp, NPPa, and NPPh had heterogeneous spatial distributions, and the average NPPp and NPPa values over the whole study area increased at rates of 3.63 and 6.94 gC·m−2·yr−1 from 2000 to 2020, respectively, while the NPPh decreased at a rate of −4.43 gC·m−2·yr−1. Climate variability and the combined effects of climate variability and human activities were the major driving factors of the NPPa increases, accounting for more than 72% of the total pixels, while the combined effects of the two factors caused the NPPa values to increase by 32–54% of the area in all cities expect Macao and across all vegetation ecosystems. Human activities often led to decreases in NPPa over more than 16% of the total pixels, and were mainly concentrated in the central cities of the GBA. The results can provide a reference for understanding NPP changes and can offer a theoretical basis for implementing ecosystem restoration, ecological construction, and conservation practices in the GBA.

Estimation of Terrestrial Net Primary Productivity in the Yellow River Basin of China Using Light Use Efficiency Model

The net primary productivity (NPP) of vegetation is an essential factor of ecosystem functions, including the biological geochemical carbon cycle, which is often impacted by climate change and human activities. It plays a significant role in comprehending the nature of carbon balance in an ecosystem and demonstrates the global and regional carbon cycle dynamics. The present study used an upgraded CASA model to calculate the NPP in the Yellow River Basin (YRB), China. The model’s simulation ability was improved by changing the model parameters. Further, the CASA model was validated by comparing with MODIS-NPP and in situ observed NPP, wherein the accuracy of the CASA model estimation was found satisfactory to estimate NPP changes in the study area. The simulated results of the improved CASA model showed that the mean annual NPP value of vegetation in the YRB was 283.4 gC m–2 a–1 from 2001 to 2020, with a declining trend in spatial distribution from south to north. In contrast, the NPP appeared as an increasing trend in the YRB temporally from 212 gC m–2 a–1 in 2001 to 342 gC m–2 a–1 in 2020, with a mean annual growth rate of 4.6 gC m–2 a–1. The total NPP in the YRB increased by 40,088.3 GgC between 2001 and 2020, from 226.06 TgC to 266.15 TgC. This rise can be attributed to the increase in forests. The average grassland area has reduced by 4651 km2 during the last two decades, significantly impacting the total NPP of grasslands. Although the increase in NPP in wetlands was minimal, accounting for 815.53 GgC, the highest change percentage of 79.78%, could be observed among the six vegetation types due to the anthropogenic influences and climate change. The conditions favorable for vegetation growth and a sustained environment were enhanced by the increased precipitation and temperature and the reinforced ecological protection by the government.

Assessing the Net Primary Productivity Dynamics of the Desert Steppe in Northern China during the Past 20 Years and Its Response to Climate Change

The net primary productivity (NPP) dynamics in arid and semi-arid ecosystems are critical for regional carbon management. Our study applied a light-utilization-efficiency model (CASA: Carnegie–Ames–Stanford Approach) to evaluate the vegetation NPP dynamics of a desert steppe in northern China over the past 20 years, and its response to climate change. Our results show that the annual average NPP of the desert steppe was 132 g C m−2 y−1, of which the grass- and shrub-dominated biome values were 142 and 91 g C m−2 y−1, respectively. The average change rate of NPP was 1.13 g C m−2 y−1, and in the grassland biome 1.31 g C m−2 y−1, a value which was significantly higher than that in shrubland, at 0.84 g C m−2 y−1. The precipitation and temperature at different time scales in the desert steppe showed a slow upward trend, and the degree of aridity tended to weaken. The correlation analysis shows that NPP changes were significantly positively and negatively correlated with precipitation and temperature, respectively. In terms of temperature, 43% of the area was significantly correlated during the growing season, which decreased to 12% on the annual scale. In 31% of the changed areas, the average NPP was 148.1 g C m−2 y−1, which was higher than the remaining significant areas. This suggests that higher NPP levels help to attenuate the negative effects of high temperature during the growing season on plant productivity in the desert steppe. This improves the understanding of the carbon cycle mechanism of arid and semi-arid ecosystems, which is beneficial to improving sustainable grassland development strategies.

Evaluation of ecosystem stability against climate changes via satellite data in the eastern sandy area of northern China

Intensive and frequent climate change events (e.g., droughts or extreme weather) significantly affect vulnerable water-limited ecosystems. Until now, the ecosystem stability against climate changes in regional scale sandy lands remain unclear. In this study, the AutoRegression (ARx) model was combined with time-series Net Primary Productivity (NPP) data to extract stability metrics (e.g., temporal stability, resilience, drought-resistance, and temperature-resistance) to evaluate the stability of the main sandy land regions of Northern China. Strong correlations among ecosystem stability metrics were found in the study area, such as the significant negative correlation between resilience and resistance (r = -0.49, p < 0.01), the strong positive correlation between drought-resistance and temperature-resistance, (r = 0.81, p < 0.01), except for the uncorrelation between resilience and temporal stability. Meanwhile, more unstable regions were found in the western low- or moderate-cover sandy grassland. Due to the differences of factors (e.g. hydrothermal conditions, vegetation species composition, and other disturbances or anthropogenic impacts), the unstable grasslands and barren regions, Otindag and Hulun Buir sandy lands, and slightly desertified area (SL) presented more resilience but less resistance and variance than the forest and cropland, Horqin Sandy Land, and Moderate (M) or Severe desertified areas (S), respectively. Thus, the unstable low-or moderate-cover grassland and SL area should be paid much more attention to meet the challenges of more intense climate extremes in the future.

Remote Sensing Estimation and Spatiotemporal Pattern Analysis of Terrestrial Net Ecosystem Productivity in China

Net ecosystem productivity (NEP) plays an important role in understanding ecosystem function and the global carbon cycle. In this paper, the key parameters of the Carnegie Ames Stanford Approach (CASA) model, maximum light use efficiency (εmax), was optimized by using vegetation classification data. Then, the NEP was estimated by coupling the optimized CASA model, geostatistical model of soil respiration (GSMSR) and the soil respiration–soil heterotrophic respiration (Rs-Rh) relationship model. The ground observations from ChinaFLUX were used to verify the NEP estimation accuracy. The results showed that the R2 of the optimized CASA model increased from 0.411 to 0.774, and RMSE decreased from 21.425 gC·m−2·month−1 to 12.045 gC·m−2·month−1, indicating that optimizing CASA model by vegetation classification data was an effective method to improve the estimation accuracy of NEP. On this basis, the spatial and temporal distribution of NEP in China was analyzed. The research indicated that the monthly variation of NEP in China was a single peak curve with summer as the peak, which generally presented the pattern of southern region > northern region > Qinghai–Tibet region > northwest region. Furthermore, from 2001 to 2016, most regions of China showed a non-significant level upward trend, but main cropland (e.g., North China Plain and Northeast Plain) and some grassland (e.g., Ngari in Qinghai–Tibet Plateau and Xilin Gol League in Inner Mongolia) showed a non-significant-level downward trend. The study can deepen the understanding of the distribution of carbon sources/sinks in China, and provide a reference for regional carbon cycle research.

Estimating mangrove forest gross primary production by quantifying environmental stressors in the coastal area

Mangrove ecosystems play an important role in global carbon budget, however, the quantitative relationships between environmental drivers and productivity in these forests remain poorly understood. This study presented a remote sensing (RS)-based productivity model to estimate the light use efficiency (LUE) and gross primary production (GPP) of mangrove forests in China. Firstly, LUE model considered the effects of tidal inundation and therefore involved sea surface temperature (SST) and salinity as environmental scalars. Secondly, the downscaling effect of photosynthetic active radiation (PAR) on the mangrove LUE was quantified according to different PAR values. Thirdly, the maximum LUE varied with temperature and was therefore determined based on the response of daytime net ecosystem exchange and PAR at different temperatures. Lastly, GPP was estimated by combining the LUE model with the fraction of absorbed photosynthetically active radiation from Sentinel-2 images. The results showed that the LUE model developed for mangrove forests has higher overall accuracy (RMSE = 0.0051, R² = 0.64) than the terrestrial model (RMSE = 0.0220, R² = 0.24). The main environmental stressor for the photosynthesis of mangrove forests in China was PAR. The estimated GPP was, in general, in agreement with the in-situ measurement from the two carbon flux towers. Compared to the MODIS GPP product, the derived GPP had higher accuracy, with RMSE improving from 39.09 to 19.05 g C/m²/8 days in 2012, and from 33.76 to 19.51 g C/m²/8 days in 2015. The spatiotemporal distributions of the mangrove GPP revealed that GPP was most strongly controlled by environmental conditions, especially temperature and PAR, as well as the distribution of mangroves. These results demonstrate the potential of the RS-based productivity model for scaling up GPP in mangrove forests, a key to explore the carbon cycle of mangrove ecosystems at national and global scales.

Estimating aboveground net primary productivity of reforested trees in an urban landscape using biophysical variables and remotely sensed data

Recently, urban reforestation programs have emerged as potential carbon sinks and climate mitigates in urban landscapes. Thus, spatially explicit information on net primary productivity (NPP) of reforested trees in urban environments is central to understanding the value of reforestation initiatives in the global carbon budget and climate regulation potential. To date, numerous studies have mainly focused on natural and commercial forests NPP at a regional scale based on coarse spatial resolution remotely sensed data. Generally, local scale NPP studies based on fine spatial resolution data are limited. Therefore, this study sought to estimate aboveground NPP of an urban reforested landscape using biophysical and Sentinel-2 Multispectral Imager data derived variables. Using the MOD17 model, results showed that mean NPP ranged between 6.24 Mg C ha^-1 with high coefficient of determination (R²: 0.92) and low RMSE (0.82 Mg ha^-1) across all reforested trees within the study area. Results also showed a considerable variation in NPP among the reforested trees, with deciduous Acacia and Dalbergia obovate species showing the highest NPP (7.62 Mg C ha^-1 and 7.58 Mg C ha^-1, respectively), while the evergreen Syzygium cordatum and shrub Artemisia afra had the lowest NPP (4.54 Mg C ha^-1 and 5.26 Mg C ha^-1). Furthermore, the multiple linear regression analysis showed that vegetation specific biophysical variables (i.e. leaf area index, Normalized Difference Vegetation Index and Fraction of Photosynthetically Active Radiation) significantly improved the estimation of reforested aboveground NPP at a fine-scale resolution. These findings demonstrate the effectiveness of biophysical and remotely sensed variables in determining NPP (as carbon sequestration surrogate) at fine-scaled reforested urban landscape. Furthermore, the utility of species biometric measurements and MOD17 model offers unprecedented opportunity for improved local scale reforestation assessment and monitoring schedules.

Nonlinear Characteristics of NPP Based on Ensemble Empirical Mode Decomposition from 1982 to 2015—A Case Study of Six Coastal Provinces in Southeast China

Monitoring vegetation net primary productivity (NPP) is very important for evaluating ecosystem health. However, the nonlinear characteristics of the vegetation NPP remain unclear in the six provinces along the Maritime Silk Road in China. In this study, using NDVI and meteorological data from 1982 to 2015, NPP was estimated with the Carnegie-Ames-Stanford Approach (CASA) model based on vegetation type dynamics, and its nonlinear characteristics were explored through the ensemble empirical mode decomposition (EEMD) method. The results showed that: (1) The total NPP in the changed vegetation types caused by ecological engineering and urbanization increased but decreased in those caused by agricultural reclamation and vegetation destruction, (2) the vegetation NPP was dominated by interannual variations, mainly in the middle of the study area, while by long-term trends, mainly in the southwest and northeast, (3) for most of the vegetation types, NPP was dominated by the monotonically increasing trend. Although vegetation NPP in the urban land mainly showed a decreasing trend (monotonic decrease and decrease from increase), there were large areas in which NPP increased from decreasing. Although vegetation NPP in the farmland mainly showed increasing trends, there were large areas that faced the risk of NPP decreasing; (4) dynamical changes of vegetation type by agricultural reclamation and vegetation destruction made the NPP trend monotonically decrease in large areas, leading to ecosystem degradation, while those caused by urbanization and ecological engineering mainly made the NPP increase from decreasing, leading to later recovery from early degradation. Our results highlighted the importance of vegetation type dynamics for accurately estimating vegetation NPP, as well as for assessing their impacts, and the importance of nonlinear analysis for deepening our understanding of vegetation NPP changes.

BERM: a Belowground Ecosystem Resiliency Model for estimating Spartina alterniflora belowground biomass

Spatiotemporal patterns of Spartina alterniflora belowground biomass (BGB) are important for evaluating salt marsh resiliency. To solve this, we created the BERM (Belowground Ecosystem Resiliency Model), which estimates monthly BGB (30-m spatial resolution) from freely available data such as Landsat-8 and Daymet climate summaries. Our modeling framework relied on extreme gradient boosting, and used field observations from four Georgia salt marshes as ground-truth data. Model predictors included estimated tidal inundation, elevation, leaf area index, foliar nitrogen, chlorophyll, surface temperature, phenology, and climate data. The final model included 33 variables, and the most important variables were elevation, vapor pressure from the previous four months, Normalized Difference Vegetation Index (NDVI) from the previous five months, and inundation. Root mean squared error for BGB from testing data was 313 g m^-2 (11% of the field data range), explained variance (R² ) was 0.62-0.77. Testing data results were unbiased across BGB values and were positively correlated with ground-truth data across all sites and years (r = 0.56-0.82 and 0.45-0.95, respectively). BERM can estimate BGB within Spartina alterniflora salt marshes where environmental parameters are within the training data range, and can be readily extended through a reproducible workflow. This provides a powerful approach for evaluating spatiotemporal BGB and associated ecosystem function.

Aboveground Biomass Mapping of Crops Supported by Improved CASA Model and Sentinel-2 Multispectral Imagery

The net primary productivity (NPP) and aboveground biomass mapping of crops based on remote sensing technology are not only conducive to understanding the growth and development of crops but can also be used to monitor timely agricultural information, thereby providing effective decision making for agricultural production management. To solve the saturation problem of the NDVI in the aboveground biomass mapping of crops, the original CASA model was improved using narrow-band red-edge information, which is sensitive to vegetation chlorophyll variation, and the fraction of photosynthetically active radiation (FPAR), NPP, and aboveground biomass of winter wheat and maize were mapped in the main growing seasons. Moreover, in this study, we deeply analyzed the seasonal change trends of crops’ biophysical parameters in terms of the NDVI, FPAR, actual light use efficiency (LUE), and their influence on aboveground biomass. Finally, to analyze the uncertainty of the aboveground biomass mapping of crops, we further discussed the inversion differences of FPAR with different vegetation indices. The results demonstrated that the inversion accuracies of the FPAR of the red-edge normalized vegetation index (NDVIred-edge) and red-edge simple ratio vegetation index (SRred-edge) were higher than those of the original CASA model. Compared with the reference data, the accuracy of aboveground biomass estimated by the improved CASA model was 0.73 and 0.70, respectively, which was 0.21 and 0.13 higher than that of the original CASA model. In addition, the analysis of the FPAR inversions of different vegetation indices showed that the inversion accuracies of the red-edge vegetation indices NDVIred-edge and SRred-edge were higher than those of the other vegetation indices, which confirmed that the vegetation indices involving red-edge information can more effectively retrieve FPAR and aboveground biomass of crops.

Ecosystem Services in Working Lands of the Southeastern USA

Agriculture and natural systems interweave in the southeastern US, including Florida, Georgia, and Alabama, where topographic, edaphic, hydrologic, and climatic gradients form nuanced landscapes. These are largely working lands under private control, comprising mosaics of timberlands, grazinglands, and croplands. According to the “ecosystem services” framework, these landscapes are multifunctional. Generally, working lands are highly valued for their provisioning services, and to some degree cultural services, while regulating and supporting services are harder to quantify and less appreciated. Trade-offs and synergies exist among these services. Regional ecological assessments tend to broadly paint working lands as low value for regulating and supporting services. But this generalization fails to consider the complexity and tight spatial coupling of land uses and land covers evident in such regions. The challenge of evaluating multifunctionality and ecosystem services is that they are not spatially concordant. While there are significant acreages of natural systems embedded in southeastern working lands, their spatial characteristics influence the balance of tradeoffs between ecosystem services at differing scales. To better understand this, we examined the configuration of working lands in the southeastern US by comparing indicators of ecosystem services at multiple scales. Indicators included measurements of net primary production (provisioning), agricultural Nitrogen runoff (regulating), habitat measured at three levels of land use intensity, and biodiversity (supporting). We utilized a hydrographic and ecoregional framework to partition the study region. We compared indicators aggregated at differing scales, ranging from broad ecoregions to local landscapes focused on the USDA Long-Term Agroecosystem Research (LTAR) Network sites in Florida and Georgia. Subregions of the southeastern US differ markedly in contributions to overall ecosystem services. Provisioning services, characterized by production indicators, were very high in northern subregions of Georgia, while supporting services, characterized by habitat and biodiversity indicators, were notably higher in smaller subregions of Florida. For supporting services, the combined contributions of low intensity working lands with embedded natural systems made a critical difference in their regional evaluation. This analysis demonstrated how the inclusion of working lands combined with examining these at different scales shifted our understanding of ecosystem services trade-offs and synergies in the southeastern United States.

A physiology-based Earth observation model indicates stagnation in the global gross primary production during recent decades

Earth observation-based estimates of global gross primary production (GPP) are essential for understanding the response of the terrestrial biosphere to climatic change and other anthropogenic forcing. In this study, we attempt an ecosystem-level physiological approach of estimating GPP using an asymptotic light response function (LRF) between GPP and incoming photosynthetically active radiation (PAR) that better represents the response observed at high spatiotemporal resolutions than the conventional light use efficiency approach. Modelled GPP is thereafter constrained with meteorological and hydrological variables. The variability in field-observed GPP, net primary productivity and solar-induced fluorescence was better or equally well captured by our LRF-based GPP when compared with six state-of-the-art Earth observation-based GPP products. Over the period 1982-2015, the LRF-based average annual global terrestrial GPP budget was 121.8 ± 3.5 Pg C, with a detrended inter-annual variability of 0.74 ± 0.13 Pg C. The strongest inter-annual variability was observed in semi-arid regions, but croplands in China and India also showed strong inter-annual variations. The trend in global terrestrial GPP during 1982-2015 was 0.27 ± 0.02 Pg C year-1 , and was generally larger in the northern than the southern hemisphere. Most positive GPP trends were seen in areas with croplands whereas negative trends were observed for large non-cropped parts of the tropics. Trends were strong during the eighties and nineties but levelled off around year 2000. Other GPP products either showed no trends or continuous increase throughout the study period. This study benchmarks a first global Earth observation-based model using an asymptotic light response function, improving simulations of GPP, and reveals a stagnation in the global GPP after the year 2000.

Niche differentiation between deeply divergent phylogenetic lineages of an endemic newt: implications for Species Distribution Models

Species distribution models (SDMs) treat species as a single unit, neglecting intraspecific variation. Few studies address the impact of intraspecific variation on SDM performance, and none of them account for the well-known inter-algorithm variability in prediction performance. The endemic Iberian amphibian Lissotriton boscai comprises two geographically highly structured phylogenetic lineages, which allowed us to explore how intraspecific variation affects the overall performance of SDMs and the predicted ecological niche. We built species and lineage distribution models using three different presence-only algorithms. We also tested for niche overlap, niche equivalency and niche similarity, using an ordination technique. We found differences in the predicted potential distribution of the two lineages and the underlying environmental factors. Moreover, intraspecific differences in model predictive capacity existed irrespective of which algorithm was used to build the distribution models. This was coupled with lineages showing a low degree of niche overlap and occurring in relatively different environmental niches spaces. The intraspecific variation observed in L. boscai led to an improved intraspecific predictivity of the lineage level based-distribution models. There was partial spatial agreement between the niche overlap and independently reported secondary contact zones. Thus SDMs built only at the species level may be too naive to predict impacts of global change on species distributions.

Forest and Crop Leaf Area Index Estimation Using Remote Sensing: Research Trends and Future Directions

Leaf area index (LAI) is an important vegetation leaf structure parameter in forest and agricultural ecosystems. Remote sensing techniques can provide an effective alternative to field-based observation of LAI. Differences in canopy structure result in different sensor types (active or passive), platforms (terrestrial, airborne, or satellite), and models being appropriate for the LAI estimation of forest and agricultural systems. This study reviews the application of remote sensing-based approaches across different system configurations (passive, active, and multisource sensors on different collection platforms) that are used to estimate forest and crop LAI and explores uncertainty analysis in LAI estimation. A comparison of the difference in LAI estimation for forest and agricultural applications given the different structure of these ecosystems is presented, particularly as this relates to spatial scale. The ease of use of empirical models supports these as the preferred choice for forest and crop LAI estimation. However, performance variation among different empirical models for forest and crop LAI estimation limits the broad application of specific models. The development of models that facilitate the strategic incorporation of local physiology and biochemistry parameters for specific forests and crop growth stages from various temperature zones could improve the accuracy of LAI estimation models and help develop models that can be applied more broadly. In terms of scale issues, both spectral and spatial scales impact the estimation of LAI. Exploration of the quantitative relationship between scales of data from different sensors could help forest and crop managers more appropriately and effectively apply different data sources. Uncertainty coming from various sources results in reduced accuracy in estimating LAI. While Bayesian approaches have proven effective to quantify LAI estimation uncertainty based on the uncertainty of model inputs, there is still a need to quantify uncertainty from remote sensing data source, ground measurements and related environmental factors to mitigate the impacts of model uncertainty and improve LAI estimation.

The effects of climate factors and human activities on net primary productivity in Xinjiang

Net primary productivity (NPP) is an index of the increase in plant biomass. Plant biomass is an important component of the global carbon cycle that indicates the health of an ecosystem. Environmental restoration has recently received much attention in Xinjiang, and it is thus important to quantify the dynamic effects of the drivers of NPP in the region. NPP was calculated for the annual growing season from 1982 to 2013 using the Carnegie-Ames-Stanford Approach (CASA) model. The effects of climate factors on NPP were analyzed, and the relationships between NPP and climate factors as well as human activity were quantified. Additionally, an innovative method based on partial derivatives and residual error was proposed to calculate the contributions of climate factors and human activities. The results show that average annual NPP in Xinjiang was 57.45 g C m^-2 from 1982 to 2013 and that the average increase in annual NPP was 0.23 g C m^-2 year^-1. The average increases in annual NPP due to temperature, precipitation, and solar radiation were 0.0095, 0.2679, and 0.2541 g C m^-2 year^-1; the average decreases were respectively - 0.0133, - 0.0521, and - 0.0725 g C m^-2 year^-1. Precipitation and solar radiation influence NPP more than temperature. Precipitation had the greatest effect on NPP in the first 19 years, but solar radiation became more influential after 2000. Climate conditions were favorable for increase in NPP before 2000. The environmental restoration also occurred in Xinjiang during that period, and human activity slightly decreased NPP. Human activity increased and had a greater effect on NPP after 2000.

Assessing the spatiotemporal variation of NPP and its response to driving factors in Anhui province, China

Net primary productivity (NPP) of terrestrial ecosystems is an important metric of ecosystem functioning; however, the understanding of response mechanism of NPP to influencing factors and driving mechanisms are still limited. In this study, taking Anhui province as an example, spatio-temporal changes of NPP and its response to influencing factors were investigated for evaluating the effects of climate change and land use and land cover change (LUCC) on regional NPP. The Carnegie-Ames-Stanford Approach (CASA) model was employed for NPP simulation by using the MODIS normalized difference vegetation index (NDVI) data and meteorological data over 2001-2016. Combined domestic LUCC, the spatiotemporal distribution pattern and dynamic change characteristics of NPP under a long time series and its response to climate factors and human activities were analyzed in the Anhui province. The results indicated that from 2001 to 2016, total NPP had a fluctuated and decreased trend with the variation range between 30.52 and 38.07 TgC in Anhui province. The multi-year average of total NPP was about 34.62 TgC. The highest value was in 2008 and the lowest value was in 2011. Among them, amount of forestland NPP was the most. The spatial distribution of NPP shows that the high value area was mainly distributed in southern Anhui mountain areas and western Anhui Dabie mountain areas; the lower value was distributed in the middle in the study area. The area of which the NPP showed a slight decrease and essentially unchanged accounted for 59.35% and 31.82%, respectively. In general, the correlation between vegetation NPP and temperature was greater than that between precipitation. The vegetations NPP of eight land use types were all positively correlated with temperature. However, the other seven types of land use were negatively correlated with precipitation except cultivated land. In the past 16 years, the decrease of cultivated land areas and the increase of urban and construction land areas contributed a lot to the decrease of vegetation NPP in Anhui province. The NPP changes of different land use types were closely related to climatic factors, land cover area, and vegetation types.

A Framework to Tackling the Synchrony between Social and Ecological Phases of the Annual Cyclic Movement of Transhumant Pastoralism

Transhumant pastoralism is mobile livelihood strategy in which families and their herds move seasonally from lowlands, where they settle during the winter, towards the highlands, located in mountainous areas, during the summer. We propose a framework, rooted in a socio-environmental coevolutionary perspective, for the transhumant annual cycle as comprised by the winter-phase, the summer-phase, and movement transitions between them. The aim was to assess the level of synchrony between ecological phases and social phases and the benefit of moving between pasturelands in selected study cases from Patagonia, Argentina. Ecological phases were addressed by the difference between vegetation productivity of winter- and summer-lands, with Fourier transform applied to data series of the Normalized Difference Vegetation Index (NDVI). Social phases were estimated by the proportion of annual time spent by pastoralists and their herds in each site and during transitions, respectively, obtained from interviews. The framework was sensitive to capturing differences across study cases. There was an observed tendency towards more synchronisation in the cases with closer distances and asynchrony in the cases with longer distances and longer movement transitions between pasturelands. Results are encouraging as a step towards the development of a monitoring system of both transhumant pastoralism activity and environmental changes.

An improved understanding of ungulate population dynamics using count data: Insights from western Montana

Understanding the dynamics of ungulate populations is critical given their ecological and economic importance. In particular, the ability to evaluate the evidence for potential drivers of variation in population trajectories is important for informed management. However, the use of age ratio data (e.g., juveniles:adult females) as an index of variation in population dynamics is hindered by a lack of statistical power and difficult interpretation. Here, we show that the use of a population model based on count, classification and harvest data can dramatically improve the understanding of ungulate population dynamics by: 1) providing estimates of vital rates (e.g., per capita recruitment and population growth) that are easier to interpret and more useful to managers than age ratios and 2) increasing the power to assess potential sources of variation in key vital rates. We used a time series of elk (Cervus canadensis) spring count and classification data (2004 to 2016) and fall harvest data from hunting districts in western Montana to construct a population model to estimate vital rates and assess evidence for an association between a series of environmental covariates and indices of predator abundance on per capita recruitment rates of elk calves. Our results suggest that per capita recruitment rates were negatively associated with cold and wet springs, and severe winters, and positively associated with summer precipitation. In contrast, an analysis of the raw age ratio data failed to detect these relationships. Our approach based on a population model provided estimates of the region-wide mean per capita recruitment rate (mean = 0.25, 90% CI = 0.21, 0.29), temporal variation in hunting-district-specific recruitment rates (minimum = 0.09; 90% CI = [0.07, 0.11], maximum = 0.43; 90% CI = [0.38, 0.48]), and annual population growth rates (minimum = 0.83; 90% CI = [0.78, 0.87], maximum = 1.20; 90% CI = [1.11, 1.29]). We recommend using routinely collected population count and classification data and a population modeling approach rather than interpreting estimated age ratios as a substantial improvement in understanding population dynamics.

Inversion of the Fraction of Absorbed Photosynthetically Active Radiation (FPAR) from FY-3C MERSI Data

An accurate inversion of the fraction of absorbed photosynthetically active radiation (FPAR) based on remote sensing data is particularly important for understanding global climate change. At present, there are relatively few studies focusing on the inversion of FPAR using Chinese autonomous satellites. This work intends to investigate the inversion of the FPAR obtained from the FengYun-3C (FY-3C) data of domestic satellites by using the PROSAIL model and the look-up table (LUT) algorithm for different vegetation types from various places in China. After analyzing the applicability of existing models using FY-3C data and MOD09GA data, an inversion strategy for FY-3C data is implemented. This strategy is applied to areas with various types of vegetation, such as grasslands, croplands, shrubs, broadleaf forests, and needleleaf forests, and produces FPAR products, which are cross-validated against the FPAR products from the Moderate Resolution Imaging Spectro Radiometer (MODIS), Geoland Version 1 (GEOV1), and Global Land Surface Satellite (GLASS). Accordingly, the results show that the FPAR retrieved from the FY-3C data has good spatial and temporal consistency and correlation with the three FPAR products. However, this technique does not favor all types of vegetation equally; the FY-FPAR is relatively more suitable for the inversion of grasslands and croplands during the lush period than for others. Therefore, the inversion strategy provides the potential to generate large-area and long-term sequence FPAR products from FY-3C data.

A Method for Estimating Annual Cumulative Soil/Ecosystem Respiration and CH4 Flux from Sporadic Data Collected Using the Chamber Method

Measurements of greenhouse gas fluxes over many ecosystems have been made as part of the attempt to quantify global carbon and nitrogen cycles. In particular, annual flux observations are of great value for regional flux assessments, as well as model development and optimization. The chamber method is a popular approach for soil/ecosystem respiration and CH4 flux observations of terrestrial ecosystems. However, in situ flux chamber measurements are usually made with non-continuous sampling. To date, efficient methods for the application of such sporadic data to upscale temporally and obtain annual cumulative fluxes have not yet been determined. To address this issue, we tested the adequacy of non-continuous sampling using multi-source data aggregation. We collected 330 site-years monthly soil/ecosystem respiration and 154 site-years monthly CH4 flux data in China, all obtained using the chamber method. The data were randomly divided into a training group and verification group. Fluxes of all possible sampling months of a year, i.e., 4094 different month combinations were used to obtain the annual cumulative flux. The results showed a good linear relationship between the monthly flux and the annual cumulative flux. The flux obtained during the warm season from May to October generally played a more important role in annual flux estimations, as compared to other months. An independent verification analysis showed that the monthly flux of 1 to 4 months explained up to 67%, 89%, 94%, and 97% of the variability of the annual cumulative soil/ecosystem respiration and 92%, 99%, 99%, and 99% of the variability of the annual cumulative CH4 flux. This study supports the use of chamber-observed sporadic flux data, which remains the most commonly-used method for annual flux estimating. The flux estimation method used in this study can be used as a guide for designing sampling programs with the intention of estimating the annual cumulative flux.

Land use management based on multi-scenario allocation and trade-offs of ecosystem services in Wafangdian County, Liaoning Province, China

Developing effective methods to coordinate the trade-offs among ecosystem services (ES) is important for achieving inclusive growth and sustainable development, and has been the focus of scholars and ecosystem managers globally. Using remote sensing and geographic information system (GIS) data, our study examined Wafangdian County of Liaoning Province as a case study to reveal the spatiotemporal evolution of four ES (food supply [FS], net primary productivity [NPP], water yield [WY], and soil conservation [SC]) and changes among their interactions. Then, an ordered weighted averaging model was introduced to simulate the optimal scenario of ES allocation. Results showed that: (1) the spatial and temporal changes in ES were significant over 14 years. All ES presented an inverted U-shaped growth curve from 2000–2014. (2) Synergies were observed within provisioning services, and there were trade-offs between provisioning services and regulating services, as well as provisioning services and supporting services. (3) The optimal scenario for Wafangdian was scenario 5 (trade-off coefficient, 0.68). The allocation of FS, NPP, WY, and SC in scenario 5 were 0.187, 0.427, 0.131, and 0.063, respectively. Implementing each ES weight of optimal scenario in land use management contributed to achieving intercoordination of ES. We propose to coordinate land and sea management to restore natural habitats that were expanded into in the high ES area. It is our anticipation that this study could provide a scientific basis for optimizing the allocation of ES and improving land use structure of coastal zones in the future.

Spatiotemporal Simulation of Net Ecosystem Productivity and Its Response to Climate Change in Subtropical Forests

Subtropical forests have great potential as carbon sinks; however, the relationship between net ecosystem productivity (NEP) and climate change is still unclear. This study took Zhejiang Province, a subtropical region, as an example. Based on remote sensing classification data of forest resources, the integrated terrestrial ecosystem carbon cycle (InTEC) model was used to simulate the spatiotemporal dynamics of the forest NEP in Zhejiang Province during 1985–2015 and analyze its response to meteorological factors such as temperature, precipitation, relative humidity, and radiation. Three patterns emerged: (1) The optimized InTEC model can better simulate the forest NEP in Zhejiang Province, and the correlation coefficient between the simulated NEP and observed NEP was up to 0.75. (2) From 1985 to 2015, the increase in the total NEP was rapid, with an average annual growth rate of 1.52 Tg·C·yr−1. During 1985–1988, the forests in Zhejiang Province were carbon sources. After 1988, the forests turned into carbon sinks and this continued to increase. During 2000–2015, more than 97% of the forests in Zhejiang Province were carbon sinks. The total NEP reached 32.02 Tg·C·yr−1, and the annual mean NEP increased to 441.91 gC·m−2·yr−1. The carbon sequestration capacity of forests in the east and southwest of Zhejiang Province is higher than that in the northeast of Zhejiang Province. (3) From 2000 to 2015, there was an extremely significant correlation between forest NEP and precipitation, with a correlation coefficient of 0.85. Simultaneously, the forest NEP showed a negative correlation with temperature and radiation, with a correlation coefficient of −0.56 for both, and the forest NEP was slightly negatively correlated with relative humidity. The relative contribution rates of temperature, precipitation, relative humidity, and radiation data to NEP showed that the contribution of precipitation to NEP is the largest, reaching 61%, followed by temperature and radiation at 18% and 17%, respectively. The relative contribution rate of relative humidity is the smallest at only 4%. During the period of 1985–1999, due to significant man-made disturbances, the NEP had a weak correlation with temperature, precipitation, relative humidity, and radiation. The results of this study are important for addressing climate change and illustrating the response mechanism between subtropical forest NEP and climate change.

Monitoring Land Cover Change and Disturbance of the Mount Wutai World Cultural Landscape Heritage Protected Area, Based on Remote Sensing Time-Series Images from 1987 to 2018

The contextual-based multi-source time-series remote sensing and proposed Comprehensive Heritage Area Threats Index (CHATI) index are used to analyze the spatiotemporal land use/land cover (LULC) and threats to the Mount Wutai World Heritage Area. The results show disturbances, such as forest coverage, vegetation conditions, mining area, and built-up area, in the research area changed dramatically. According to the CHATI, although different disturbances have positive or negative influences on environment, as an integrated system it kept stable from 1987 to 2018. Finally, this research uses linear regression and the F-test to mark the remarkable spatial-temporal variation. In consequence, the threats on Mount Wutai be addressed from the macro level and the micro level. Although there still have some drawbacks, the effectiveness of threat identification has been tested using field validation and the results are a reliable tool to raise the public awareness of WHA protection and governance.

Comparison of Normalized Difference Vegetation Index Derived from Landsat, MODIS, and AVHRR for the Mesopotamian Marshes Between 2002 and 2018

The Mesopotamian marshes are a group of water bodies located in southern Iraq, in the shape of a triangle, with the cities Amarah, Nasiriyah, and Basra located at its corners. The marshes are appropriate habitats for a variety of birds and most of the commercial fisheries in the region. The normalized difference vegetation index (NDVI) has been derived using observations from various satellite sensors, such as the Moderate Resolution Imaging Spectroradiometer (MODIS), Advanced Very-High-Resolution Radiometer (AVHRR), and Landsat over the Mesopotamian marshlands for the 17-year period between 2002 and 2018. We have chosen this time series (2002–2018) to monitor the change in vegetation of the study area since it is considered as a period of rehabilitation for the marshes (following a period when there was little to no water flowing into the marshes). Statistical analyses were performed to monitor the variability of the maximum biomass time (month of June). The results illustrated a strong positive correlation between the NDVI derived from Landsat, MODIS, and AVHRR. The statistical correlations were 0.79, 0.77, and 0.96 between Landsat and AVHRR, MODIS and AVHRR, and Landsat and MODIS, respectively. The linear slope of NDVI (Landsat, MODIS, and AVHRR) for each pixel over the period 2002–2018 displays a long-term trend of green biomass (NDVI) change in the study area, and the slope is slightly negative over most of the area. Slope values (−0.002 to −0.05) denote a slight decrease in the observed vegetation index over 17 years. The green biomass of the marshlands increased by 33.2% of the total area over 17 years. The areas of negative and positive slopes correspond to the same areas in slope map when calculated from Landsat, MODIS, and AVHRR, although they are different in spatial resolution (30 m, 1 km, and 5 km, respectively). The time series of the average NDVI (2002–2018) for three different sensors shows the highest and lowest NDVI values during the same years (for the month of June each year). The highest values were 0.19, 0.22, and 0.22 for Landsat, MODIS, and AVHRR, respectively, in 2006, and the lowest values were 0.09, 0.14, and 0.09 for Landsat, MODIS, and AVHRR, respectively, in 2003.

An Improved CASA Model for Estimating Winter Wheat Yield from Remote Sensing Images

The accurate and timely monitoring and evaluation of the regional grain crop yield is more significant for formulating import and export plans of agricultural products, regulating grain markets and adjusting the planting structure. In this study, an improved Carnegie–Ames–Stanford approach (CASA) model was coupled with time-series satellite remote sensing images to estimate winter wheat yield. Firstly, in 2009 the entire growing season of winter wheat in the two districts of Tongzhou and Shunyi of Beijing was divided into 54 stages at five-day intervals. Net Primary Production (NPP) of winter wheat was estimated by the improved CASA model with HJ-1A/B satellite images from 39 transits. For the 15 stages without HJ-1A/B transit, MOD17A2H data products were interpolated to obtain the spatial distribution of winter wheat NPP at 5-day intervals over the entire growing season of winter wheat. Then, an NPP-yield conversion model was utilized to estimate winter wheat yield in the study area. Finally, the accuracy of the method to estimate winter wheat yield with remote sensing images was verified by comparing its results to the ground-measured yield. The results showed that the estimated yield of winter wheat based on remote sensing images is consistent with the ground-measured yield, with R2 of 0.56, RMSE of 1.22 t ha−1, and an average relative error of −6.01%. Based on time-series satellite remote sensing images, the improved CASA model can be used to estimate the NPP and thereby the yield of regional winter wheat. This approach satisfies the accuracy requirements for estimating regional winter wheat yield and thus may be used in actual applications. It also provides a technical reference for estimating large-scale crop yield.

Response of net primary production to land use and climate changes in the middle-reaches of the Heihe River Basin

Net primary production (NPP) supplies matter, energy, and services to facilitate the sustainable development of human society and ecosystem. The response mechanism of NPP to land use and climate changes is essential for food security and biodiversity conservation but lacks a comprehensive understanding, especially in arid and semi-arid regions. To this end, taking the middle-reaches of the Heihe River Basin (MHRB) as an example, we uncovered the NPP responses to land use and climate changes by integrating multisource data (e.g., MOD17A3 NPP, land use, temperature, and precipitation) and multiple methods. The results showed that (a) land use intensity (LUI) increased, and climate warming and wetting promoted NPP. From 2000 to 2014, the LUI, temperature, and precipitation of MHRB increased by 1.46, 0.58°C, and 15.76 mm, respectively, resulting in an increase of 14.62 gC/m2 in annual average NPP. (b) The conversion of low-yield cropland to forest and grassland increased NPP. Although the widespread conversion of unused land and grassland to cropland boosted both LUI and NPP, it was not conducive to ecosystem sustainability and stability due to huge water consumption and human-appropriated NPP. Urban sprawl occupied cropland, forest, and grassland and reduced NPP. (c) Increase in temperature and precipitation generally improved NPP. The temperature decreasing <1.2°C also promoted the NPP of hardy vegetation due to the simultaneous precipitation increase. However, warming-induced water stress compromised the NPP in arid sparse grassland and deserts. Cropland had greater NPP and NPP increase than natural vegetation due to the irrigation, fertilizers, and other artificial inputs it received. The decrease in both temperature and precipitation generally reduced NPP, but the NPP in the well-protection or less-disturbance areas still increased slightly.

Improving Forest Aboveground Biomass (AGB) Estimation by Incorporating Crown Density and Using Landsat 8 OLI Images of a Subtropical Forest in Western Hunan in Central China

Forest aboveground biomass (AGB) estimation modeling based on remote sensing is an important method for large-scale biomass estimation; the accuracy of the estimation models has been a topic of broad and current interest. In this study, we used permanent sample plot data and Landsat 8 Operational Land Imager (OLI) images of western Hunan. Remote-sensing-based models were developed for different vegetation types, and different crown density classes were incorporated. The linear model, linear dummy variable model, and linear mixed-effects model were used to determine the most effective and accurate method for remote-sensing-based AGB estimation. The results show that the adjusted coefficient of determination (R2adj) and root mean square error (RMSE) of the linear dummy model and linear mixed-effects model were significantly better than those of the linear model; the R2adj increased more than 0.16 and the RMSE decreased more than 2.12 for each vegetation type, and the F-test also showed significant differences between the linear model and linear dummy variable model and between the linear model and linear mixed-effects model. The accuracies of the AGB estimations of the linear dummy variable model and the linear mixed-effects model were significantly better than those of linear model in the thin and dense crown density classes. There were no significant differences in the AGB estimation performance between the linear dummy variable model and linear mixed-effects model; these two models were more flexible and more suitable than the linear model for remote-sensing-based AGB estimation. The results of this study provide a new approach for solving the low-accuracy estimations of linear models.

Temporal-Spatial Evolution of the Urban Ecological Footprint Based on Net Primary Productivity: A Case Study of Xuzhou Central Area, China

The urbanization process all over the world has caused serious ecological and environmental problems which have recently become a focus for study. Ecological footprint analysis, which is widely used to assess the sustainability of regional development, can quantitatively measure the human occupation of natural capital. In this study, the ecological footprint based on net primary production (EF-NPP) and MODIS data were used to measure the ecological footprint in Xuzhou central area from 2005 to 2014. The results showed that from 2005 to 2014, the per capita ecological footprint increased from 1.06 to 1.17 hm2/person; the per capita ecological capacity decreased from 0.10 to 0.09 hm2/person; the per capita ecological deficit increased from −0.96 to −1.09 hm2/person; and the ecological pressure index increased from 6.87 to 11.97. The composition of the ecological footprint showed that grassland contributed most to the ecological footprint and deficit, and cultivated land contributed most to the ecological capacity. The spatial distribution of the ecological footprint changed significantly, especially in the expansion of the area of lower value. The ecological capacity and deficit changed little. The ecological situation in Xuzhou central area was unbalanced. Based on this study, Xuzhou city was recommended to control the increase of the ecological footprint, improve the ecological capacity and balance the ecological pattern for sustainable development.

Spatial Responses of Net Ecosystem Productivity of the Yellow River Basin under Diurnal Asymmetric Warming

The net ecosystem productivity (NEP) of drainage basins plays an important role in maintaining the carbon balance of those ecosystems. In this study, the modified CASA (Carnegie Ames Stanford Approach) model and a soil microbial respiration model were used to estimate net primary productivity (NPP) and NEP of the Yellow River Basin’s (YRB) vegetation in the terrestrial ecosystem (excluding rivers, floodplain lakes and other freshwater ecosystems) from 1982 to 2015. After analyzing the spatiotemporal variations in the NEP using slope analysis, the coefficient of variation, and the Hurst exponent, precipitation was identified as the main factor limiting vegetation growth in the YRB. Hence, precipitation was treated as the control variable and a second-order partial correlation method was used to determine the correlation between diurnal asymmetric warming and the YRB’s NEP. The results indicate that: (i) diurnal asymmetric warming occurred in the YRB from 1982 to 2015, with nighttime warming (Tmin) being 1.50 times that of daytime warming (Tmax). There is a significant correlation between variations in NPP and diurnal warming; (ii) the YRB’s NEP are characterized by upward fluctuations in terms of temporal variations, large differences between the various vegetation types, high values in the western and southeastern regions but low values in the northern region in terms of spatial distribution, overall relative stability in the YRB’s vegetation cover, and changes in the same direction being more dominant than those in the opposite direction (although the former is not sustained); and (iii) positive correlations between the NEP and nighttime and daytime warming are approximately 48.37% and 67.51% for the YRB, respectively, with variations in nighttime temperatures having more extensive impacts on vegetation cover.

Improving remote estimation of winter crops gross ecosystem production by inclusion of leaf area index in a spectral model

The hysteresis of the seasonal relationships between vegetation indices (VIs) and gross ecosystem production (GEP) results in differences between these relationships during vegetative and reproductive phases of plant development cycle and may limit their applicability for estimation of croplands productivity over the entire season. To mitigate this problem and to increase the accuracy of remote sensing-based models for GEP estimation we developed a simple empirical model where greenness-related VIs are multiplied by the leaf area index (LAI). The product of this multiplication has the same seasonality as GEP, and specifically for vegetative periods of winter crops, it allowed the accuracy of GEP estimations to increase and resulted in a significant reduction of the hysteresis of VIs vs. GEP. Our objective was to test the multiyear relationships between VIs and daily GEP in order to develop more general models maintaining reliable performance when applied to years characterized by different climatic conditions. The general model parametrized with NDVI and LAI product allowed to estimate daily GEP of winter and spring crops with an error smaller than 14%, and the rate of GEP over- (for spring barley) or underestimation (for winter crops and potato) was smaller than 25%. The proposed approach may increase the accuracy of crop productivity estimation when greenness VIs are saturating early in the growing season.