The co-evolution of multi-cities has emerged as the primary form of urbanization in China in recent years. However, the processes, patterns, and coordination are not well characterized and understood, which hinders the understanding of the driving forces, consequences, and management of polycentric urban development. We used the Continuous Change Detection and Classification (CCDC) algorithm to integrate all available Landsat 5, 7, and 8 images and map annual land use and land cover (LULC) from 2001 to 2017 in the Chang–Zhu–Tan urban agglomeration (CZTUA), a typical urban agglomeration in China. Results showed that the impervious surface in the study area expanded by 371 km² with an annual growth rate of 2.25%, primarily at the cost of cropland (169 km²) and forest (206 km²) during the study period. Urban growth has evolved from infilling being the dominant type in the earlier period to mainly edge-expansion and leapfrogging in the core cities, and from no dominant type to mainly leapfrogging in the satellite cities. The unfolding of the “cool center and hot edge” urban growth pattern in CZTUA, characterized by higher expansion rates in the peripheral than in the core cities, may signify a new form of the co-evolution of multi-cities in the process of urbanization. Detailed urban management and planning policies in CZTUA were analyzed. The co-evolution of multi-cities principles need to be studied in more extensive regions, which could help policymakers to promote sustainable and livable development in the future.

The impacts of ice-storms on forests have received growing attention in recent years. Although there is a wide agreement that ice-storms significantly affect forest structure and functions, how frequent ice-storms and subsequent salvage logging impact productivity of subtropical coniferous forests in the future still remains poorly understood. In this study, we used the Ecosystem Demography model, Version 2.2 (ED-2.2), to project the impact of salvage logging of ice-storm-damaged trees on the productivity of Cunninghamia lanceolata-dominated coniferous forest and C. lanceolata-dominated mixed coniferous and broadleaved forests. The results show that forest productivity recovery is delayed in coniferous forests when there is no shade-tolerant broadleaved species invasion after ice-storms, and C. lanceolata could continue to dominate the canopy in the mixed coniferous and broadleaved forests under high-frequency ice-storms and subsequent salvage logging. The resistance and resilience of the mixed coniferous and broadleaved forests to high-frequency ice-storms and subsequent salvage logging were stronger compared to coniferous forests. Although conifers could continue to dominate the canopy under shade-tolerant broadleaved species invasion, we could not rule out the possibility of a future forest community dominated by shade-tolerant broadleaf trees because there were few coniferous saplings and shade-tolerant broadleaf species dominated the understory. Our results highlight that post-disaster forest management should be continued after high-frequency ice-storms and subsequent salvage logging in C. lanceolata forests to prevent possible shade-tolerant, late successional broadleaf trees from dominating the canopy in the future.

The impacts of drought and/or warming on forests have received great attention in recent decades. Although the extreme drought and/or warming events significantly changed the forest demography and regional carbon cycle, the seasonality quantifying the impacts of these climate extremes with different severities on the productivity of subtropical coniferous forests remains poorly understood. This study evaluated the effects of seasonal drought and/or warming on the net primary productivity (NPP) of subtropical coniferous forests (i.e., Cunninghamia lanceolata and Pinus massoniana forests) from Hengyang–Shaoyang Basin in southern China using the Ecosystem Demography model, Version 2.2 (ED-2.2) and based on the datasets from forest inventory, meteorological reanalysis, and remotely sensed products. The results showed that the goodness of fit of the DBH-height allometric equations was better than that of the default in ED-2.2 after model calibration; the ED-2.2 model qualitatively captured the seasonality of NPP in the subtropical coniferous forests; and the mismatch between simulated annual NPP and MODIS-NPP (MOD17A3HGF) became smaller over time. The effect of seasonal drought on NPP was greater than that of warming; the decline rate of NPP gradually increased and decreased with time (from July to October) under the seasonal drought and warming scenarios, respectively; NPP decreased more seriously under the combined drought-warming scenario in October, with an average decrease of 31.72%, than the drought-only and warming-only scenarios; seasonal drought had an obvious legacy impact on productivity recovery of subtropical coniferous forests, but it was not the case for warming. With the increase in drought severity, the average values of soil available water and NPP together showed a downward trend. With the increase in warming severity, the average values of canopy air space temperature increased, but NPP decreased. Seasonal drought and/or warming limit forest production through decreasing soil moisture and/or increasing canopy air space temperature, which impact on plant photosynthesis and productivity, respectively. Our results highlight the significance of taking into account the impacts of seasonal warming and drought when evaluating the productivity of subtropical coniferous forests, as well as the significance of enhancing the resistance and resilience of forests to future, more severe global climate change.