Above-ground biomass (AGB) directly relates to the productivity of forests. Precisely, AGB mapping for regional forests based on very high resolution (VHR) imagery is widely needed for evaluation of productivity. However, the diversity of variables and algorithms and the difficulties inherent in high resolution optical imagery make it complex. In this paper, we explored the potentials of the state-of-art algorithm convolutional neural networks (CNNs), which are widely used for its high-level representation, but rarely applied for AGB estimation. Four experiments were carried out to compare the performance of CNNs and other state-of-art Machine Learning (ML) algorithms: (1) performance of CNN using bands, (2) performance of Random Forest (RF), support vector regression (SVR), artificial neural network (ANN) on bands, and vegetation indices (VIs). (3) Performance of RF, SVR, and ANN on gray-level co-occurrence matrices (GLCM), and exploratory spatial data analysis (ESDA), and (4) performance of RF, SVR, and ANN based on all combined data and ESDA+VIs. CNNs reached satisfactory results (with R² = 0.943) even with limited input variables (i.e., only bands). In comparison, RF and SVR with elaborately designed data obtained slightly better accuracy than CNN. For examples, RF based on GLCM textures reached an R² of 0.979 and RF based on all combined data reached a close R² of 0.974. However, the results of ANN were much worse (with the best R² of 0.885).

The subtropical vegetation plays an important role in maintaining the structure and function of global ecosystems, and its contribution to the global carbon balance are receiving increasing attention. The fractional vegetation cover (FVC) as an important indicator for monitoring environment change, is widely used to analyze the spatiotemporal pattern of regional and even global vegetation. China is an important distribution area of subtropical vegetation. Therefore, we first used the dimidiate pixel model to extract the subtropical FVC of China during 2001–2018 based on MODIS land surface reflectance data, and then used the linear regression analysis and the variation coefficient to explore its spatiotemporal variations characteristics. Finally, the partial correlation analysis and the partial derivative model were used to analyze the influences and contributions of climate factors on FVC, respectively. The results showed that (1) the subtropical FVC had obvious spatiotemporal heterogeneity; the FVC high-coverage and medium-coverage zones were concentratedly and their combined area accounted for more than 70% of the total study area. (2) The interannual variation in the average subtropical FVC from 2001 to 2018 showed a significant growth trend. (3) In 76.28% of the study area, the regional FVC showed an increasing trend, and the remaining regional FVC showed a decreasing trend. However, the overall fluctuations in the FVC (increasing or decreasing) in the region were relatively stable. (4) The influences of climate factors to the FVC exhibited obvious spatial differences. More than half of all pixels exhibited the influence of the average annual minimum temperature and the annual precipitation had positive on FVC, while the average annual maximum temperature had negative on FVC. (5) The contributions of climate changes to FVC had obvious heterogeneity, and the average annual minimum temperature was the main contribution factor affecting the dynamic variations of FVC.

Bamboo forests are widespread in subtropical areas and are well known for their rapid growth and great carbon sequestration ability. To recognize the potential roles and functions of bamboo forests in regional ecosystems, forest aboveground biomass (AGB)—which is closely related to forest productivity, the forest carbon cycle, and, in particular, carbon sinks in forest ecosystems—is calculated and applied as an indicator. Among the existing studies considering AGB estimation, linear or nonlinear regression models are the most frequently used; however, these methods do not take the influence of spatial heterogeneity into consideration. A geographically weighted regression (GWR) model, as a spatial local model, can solve this problem to a certain extent. Based on Landsat 8 OLI images, we use the Random Forest (RF) method to screen six variables, including TM457, TM543, B7, NDWI, NDVI, and W7B6VAR. Then, we build the GWR model to estimate the bamboo forest AGB, and the results are compared with those of the cokriging (COK) and orthogonal least squares (OLS) models. The results show the following: (1) The GWR model had high precision and strong prediction ability. The prediction accuracy (R²) of the GWR model was 0.74, 9%, and 16% higher than the COK and OLS models, respectively, while the error (RMSE) was 7% and 12% lower than the errors of the COK and OLS models, respectively. (2) The bamboo forest AGB estimated by the GWR model in Zhejiang Province had a relatively dense spatial distribution in the northwestern, southwestern, and northeastern areas. This is in line with the actual bamboo forest AGB distribution in Zhejiang Province, indicating the potential practical value of our study. (3) The optimal bandwidth of the GWR model was 156 m. By calculating the variable parameters at different positions in the bandwidth, close attention is given to the local variation law in the estimation of the results in order to reduce the model error.

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⁻¹. 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⁻¹, and the annual mean NEP increased to 441.91 gC·m⁻²·yr⁻¹. 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.

Urbanization inevitably poses a threat to urban ecology by altering its external structure and internal attributes. Nighttime light (NTL) has become increasingly extensive and practical, offering a special perspective on the world in revealing urbanization. In this study, we applied the Normalized Impervious Surface Index (NISI) constructed by NTL and MODIS NDVI to examine the urbanization process in the Yangtze River Delta (YRD). Geographical detectors combined with factors involving human and natural influences were utilized to investigate the drive mechanism. Urban ecology stress was evaluated based on changes in urban morphological patterns and fractional vegetation cover (FVC). The results showed that the NISI can largely overcome the obstacle of directly coupling NTL data in performing urbanization and has efficient applicability in the long-term pixel scale. Built-up areas in the YRD increased by 2.83 times during the past two decades, from 2053.5 to 7872.5 km². Urbanization intensity has saturated the city center and is spilling over into the suburbs, which show a “cold to hot” spatial clustering distribution. Economic factors are the primary forces driving urbanization, and road network density is becoming essential as factor that reflects urban infrastructure. Urban geometry pattern changes in fractal dimension (FD) and compactness revealed the ecological stress from changing urban external structure, and internal ecological stress was clear from the negative effect on 63.4% FVC. This impact gradually increased in urban expanded area and synchronously decreased when urbanization saturated the core area. An analysis of ecological stress caused by urbanization from changing physical structure and social attributes can provide evidence for urban management and coordinated development.

Dynamic monitoring of carbon storage in forests resources is important for tracking ecosystem functionalities and climate change impacts. In this study, we used multi-year Landsat data combined with a Random Forest (RF) algorithm to estimate the forest aboveground carbon (AGC) in a forest area in China (Hang-Jia-Hu) and analyzed its spatiotemporal changes during the past two decades. Maximum likelihood classification was applied to make land-use maps. Remote sensing variables, such as the spectral band, vegetation indices, and derived texture features, were extracted from 20 Landsat TM and OLI images over five different years (2000, 2004, 2010, 2015, and 2018). These variables were subsequently selected according to their importance and subsequently used in the RF algorithm to build an estimation model of forest AGC. The results showed the following: (1) Verification of classification results showed maximum likelihood can extract land information effectively. Our land cover classification yielded overall accuracies between 86.86% and 89.47%. (2) Additionally, our RF models showed good performance in predicting forest AGC, with R² from 0.65 to 0.73 in the training and testing phase and a RMSE range between 3.18 and 6.66 Mg/ha. RMSEr in the testing phase ranged from 20.27 to 22.27 with a low model error. (3) The estimation results indicated that forest AGC in the past two decades increased with density at 10.14 Mg/ha, 21.63 Mg/ha, 26.39 Mg/ha, 29.25 Mg/ha, and 44.59 Mg/ha in 2000, 2004, 2010, 2015, and 2018. The total forest AGC storage had a growth rate of 285%. (4) Our study showed that, although forest area decreased in the study area during the time period under study, the total forest AGC increased due to an increment in forest AGC density. However, such an effect is overridden in the vicinity of cities by intense urbanization and the loss of forest covers. Our study demonstrated that the combined use of remote sensing data and machine learning techniques can improve our ability to track the forest changes in support of regional natural resource management practices.