Carbon fluxes of China's coastal wetlands and impacts of reclamation and restoration.

Coastal wetlands play an important role in regulating atmospheric carbon dioxide (CO2) concentrations and contribute significantly to climate change mitigation. However, climate change, reclamation, and restoration have been causing substantial changes in coastal wetland areas and carbon exchange in China during recent decades. Here we compiled a carbon flux database consisting of 15 coastal wetland sites to assess the magnitude, patterns, and drivers of carbon fluxes and to compare fluxes among contrasting natural, disturbed, and restored wetlands. The natural coastal wetlands have the average net ecosystem exchange of CO2 (NEE) of -577 g C m-2 year-1, with -821 g C m-2 year-1 for mangrove forests and -430 g C m-2 year-1 for salt marshes. There are pronounced latitudinal patterns for carbon dioxide exchange of natural coastal wetlands: NEE increased whereas gross primary production (GPP) and respiration of ecosystem decreased with increasing latitude. Distinct environmental factors drive annual variations of GPP between mangroves and salt marshes; temperature was the dominant controlling factor in salt marshes, while temperature, precipitation, and solar radiation were co-dominant in mangroves. Meanwhile, both anthropogenic reclamation and restoration had substantial effects on coastal wetland carbon fluxes, and the effect of the anthropogenic perturbation in mangroves was more extensive than that in salt marshes. Furthermore, from 1980 to 2020, anthropogenic reclamation of China's coastal wetlands caused a carbon loss of ~3720 Gg C, while the mangrove restoration project during the period of 2021-2025 may switch restored coastal wetlands from a carbon source to carbon sink with a net carbon gain of 73 Gg C. The comparison of carbon fluxes among these coastal wetlands can improve our understanding of how anthropogenic perturbation can affect the potentials of coastal blue carbon in China, which has implications for informing conservation and restoration strategies and efforts of coastal wetlands.

Deeper topsoils enhance ecosystem productivity and climate resilience in arid regions, but not in humid regions.

Understanding the controlling mechanisms of soil properties on ecosystem productivity is essential for sustaining productivity and increasing resilience under a changing climate. Here we investigate the control of topsoil depth (e.g., A horizons) on long-term ecosystem productivity. We used nationwide observations (n = 2401) of topsoil depth and multiple scaled datasets of gross primary productivity (GPP) for five ecosystems (cropland, forest, grassland, pasture, shrubland) over 36 years (1986-2021) across the conterminous USA. The relationship between topsoil depth and GPP is primarily associated with water availability, which is particularly significant in arid regions under grassland, shrubland, and cropland (r = .37, .32, .15, respectively, p < .0001). For every 10 cm increase in topsoil depth, the GPP increased by 114 to 128 g C m-2 year-1 in arid regions (r = .33 and .45, p < .0001). Paired comparison of relatively shallow and deep topsoils while holding other variables (climate, vegetation, parent material, soil type) constant showed that the positive control of topsoil depth on GPP occurred primarily in cropland (0.73, confidence interval of 0.57-0.84) and shrubland (0.75, confidence interval of 0.40-0.94). The GPP difference between deep and shallow topsoils was small and not statistically significant. Despite the positive control of topsoil depth on productivity in arid regions, its contribution (coefficients: .09-.33) was similar to that of heat (coefficients: .06-.39) but less than that of water (coefficients: .07-.87). The resilience of ecosystem productivity to climate extremes varied in different ecosystems and climatic regions. Deeper topsoils increased stability and decreased the variability of GPP under climate extremes in most ecosystems, especially in shrubland and grassland. The conservation of topsoil in arid regions and improvements of soil depth representation and moisture-retention mechanisms are critical for carbon-sequestration ecosystem services under a changing climate. These findings and relationships should also be included in Earth system models.

Decipher soil organic carbon dynamics and driving forces across China using machine learning.

The dynamics of soil organic carbon (SOC) play a critical role in modulating global warming. However, the long-term spatiotemporal changes of SOC at large scale, and the impacts of driving forces remain unclear. In this study, we investigated the dynamics of SOC in different soil layers across China through the1980s to 2010s using a machine learning approach and quantified the impacts of the key factors based on factorial simulation experiments.Our results showed that the latest (2000-2014) SOC stock in the first meter soil (SOC100 ) was 80.68 ± 3.49 Pg C, of which 42.6% was stored in the top 20 cm, sequestrating carbon with a rate of 30.80 ± 12.37 g C m-2  yr-1 since the 1980s. Our experiments focusing on the recent two periods (2000s and 2010s) revealed that climate change exerted the largest relative contributions to SOC dynamics in both layers and warming or drying can result in SOC loss. However, the influence of climate change weakened with soil depth, while the opposite for vegetation growth. Relationships between SOC and forest canopy height further confirmed this strengthened impact of vegetation with soil depth and highlighted the carbon sink function of deep soil in mature forest. Moreover, our estimates suggested that SOC dynamics in 71% of topsoil were controlled by climate change and its coupled influence with environmental variation (CE). Meanwhile, CE and the combined influence of climate change and vegetation growth dominated the SOC dynamics in 82.05% of the first meter soil. Additionally, the national cropland topsoil organic carbon increased with a rate of 23.6 ± 7.6 g C m-2  yr-1 since the 1980s, and the widely applied nitrogenous fertilizer was a key stimulus. Overall, our study extended the knowledge about the dynamics of SOC and deepened our understanding about the impacts of the primary factors.

Soil moisture regulates warming responses of autumn photosynthetic transition dates in subtropical forests.

Autumn phenology plays a key role in regulating the terrestrial carbon and water balance and their feedbacks to the climate. However, the mechanisms underlying autumn phenology are still poorly understood, especially in subtropical forests. In this study, we extracted the autumn photosynthetic transition dates (APTD) in subtropical China over the period 2003-2017 based on a global, fine-resolution solar-induced chlorophyll fluorescence (SIF) dataset (GOSIF) using four fitting methods, and then explored the temporal-spatial variations of APTD and its underlying mechanisms using partial correlation analysis and machine learning methods. We further predicted the APTD shifts under future climate warming conditions by applying process-based and machine learning-based models. We found that the APTD was significantly delayed, with an average rate of 7.7 days per decade, in subtropical China during 2003-2017. Both partial correlation analysis and machine learning methods revealed that soil moisture was the primary driver responsible for the APTD changes in southern subtropical monsoon evergreen forest (SEF) and middle subtropical evergreen forest (MEF), whereas solar radiation controlled the APTD variations in the northern evergreen-broadleaf deciduous mixed forest (NMF). Combining the effects of temperature, soil moisture and radiation, we found a significantly delayed trend in APTD during the 2030-2100 period, but the trend amplitude (0.8 days per decade) was much weaker than that over 2003-2017. In addition, we found that machine learning methods outperformed process-based models in projecting APTD. Our findings generate from different methods highlight that soil moisture is one of the key players in determining autumn photosynthetic phenological processes in subtropical forests. To comprehensively understand autumn phenological processes, in-situ manipulative experiments are urgently needed to quantify the contributions of different environmental and physiological factors in regulating plants' response to ongoing climate change.

Disentangling the role of photosynthesis and stomatal conductance on rising forest water-use efficiency.

Multiple lines of evidence suggest that plant water-use efficiency (WUE)-the ratio of carbon assimilation to water loss-has increased in recent decades. Although rising atmospheric CO2 has been proposed as the principal cause, the underlying physiological mechanisms are still being debated, and implications for the global water cycle remain uncertain. Here, we addressed this gap using 30-y tree ring records of carbon and oxygen isotope measurements and basal area increment from 12 species in 8 North American mature temperate forests. Our goal was to separate the contributions of enhanced photosynthesis and reduced stomatal conductance to WUE trends and to assess consistency between multiple commonly used methods for estimating WUE. Our results show that tree ring-derived estimates of increases in WUE are consistent with estimates from atmospheric measurements and predictions based on an optimal balancing of carbon gains and water costs, but are lower than those based on ecosystem-scale flux observations. Although both physiological mechanisms contributed to rising WUE, enhanced photosynthesis was widespread, while reductions in stomatal conductance were modest and restricted to species that experienced moisture limitations. This finding challenges the hypothesis that rising WUE in forests is primarily the result of widespread, CO2-induced reductions in stomatal conductance.

Seasonal changes in GPP/SIF ratios and their climatic determinants across the Northern Hemisphere.

Satellite-derived sun-induced chlorophyll fluorescence (SIF) has been increasingly used for estimating gross primary production (GPP). However, the relationship between SIF and GPP has not been well defined, impeding the translation of satellite observed SIF to GPP. Previous studies have generally assumed a linear relationship between SIF and GPP at daily and longer time scales, but support for this assumption is lacking. Here, we used the GPP/SIF ratio to investigate seasonal variations in the relationship between SIF and GPP over the Northern Hemisphere (NH). Based on multiple SIF products and MODIS and FLUXCOM GPP data, we found strong seasonal hump-shaped patterns for the GPP/SIF ratio over northern latitudes, with higher values in the summer than in the spring or autumn. This hump-shaped GPP/SIF seasonal variation was confirmed by examining different SIF products and was evident for most vegetation types except evergreen broadleaf forests. The seasonal amplitude of the GPP/SIF ratio decreased from the boreal/arctic region to drylands and the tropics. For most of the NH, the lowest GPP/SIF values occurred in October or September, while the maximum GPP/SIF values were evident in June and July. The most pronounced seasonal amplitude of GPP/SIF occurred in intermediate temperature and precipitation ranges. GPP/SIF was positively related to temperature in the early and late parts of the growing season, but not during the peak growing months. These shifting relationships between temperature and GPP/SIF across different months appeared to play a key role in the seasonal dynamics of GPP/SIF. Several mechanisms may explain the patterns we observed, and future research encompassing a broad range of climate and vegetation settings is needed to improve our understanding of the spatial and temporal relationships between SIF and GPP. Nonetheless, the strong seasonal variation in GPP/SIF we identified highlights the importance of incorporating this behavior into SIF-based GPP estimations.

Moisture availability mediates the relationship between terrestrial gross primary production and solar-induced chlorophyll fluorescence: Insights from global-scale variations.

Effective use of solar-induced chlorophyll fluorescence (SIF) to estimate and monitor gross primary production (GPP) in terrestrial ecosystems requires a comprehensive understanding and quantification of the relationship between SIF and GPP. To date, this understanding is incomplete and somewhat controversial in the literature. Here we derived the GPP/SIF ratio from multiple data sources as a diagnostic metric to explore its global-scale patterns of spatial variation and potential climatic dependence. We found that the growing season GPP/SIF ratio varied substantially across global land surfaces, with the highest ratios consistently found in boreal regions. Spatial variation in GPP/SIF was strongly modulated by climate variables. The most striking pattern was a consistent decrease in GPP/SIF from cold-and-wet climates to hot-and-dry climates. We propose that the reduction in GPP/SIF with decreasing moisture availability may be related to stomatal responses to aridity. Furthermore, we show that GPP/SIF can be empirically modeled from climate variables using a machine learning (random forest) framework, which can improve the modeling of ecosystem production and quantify its uncertainty in global terrestrial biosphere models. Our results point to the need for targeted field and experimental studies to better understand the patterns observed and to improve the modeling of the relationship between SIF and GPP over broad scales.

Potential of hotspot solar-induced chlorophyll fluorescence for better tracking terrestrial photosynthesis.

Remote sensing of solar-induced fluorescence (SIF) opens a new window for quantifying a key ecological variable, the terrestrial ecosystem gross primary production (GPP), because of the revealed strong SIF-GPP correlation. However, similar to many other remotely sensed metrics, SIF observations suffer from the sun-sensor geometry effects, which may have important impacts on the SIF-GPP relationship but remain poorly understood. Here we used remotely sensed SIF, globally distributed tower GPP data, and a mechanistic model to provide a systematic analysis. Our results reveal that leaf physiology, canopy structure, and sun-sensor geometries all affect the SIF-GPP relationship. In particular, we found that SIF observations in the sun-tracking hotspot direction can be a better proxy of GPP due to the similar responses of light use efficiency and SIF escaping probability in the hotspot direction to the increasing incoming solar radiation. Such conclusions are supported by a variety of modeling simulations and satellite observations over various plant function types, at different time scales and with satellite observational modes. This study demonstrates the potential and advantage of normalizing SIF observations to the hotspot direction for better global GPP estimations. This study also demonstrates the great potentials of current and future spaceborne sun-tracking satellite missions for a significant improvement in measuring and monitoring, at a wide range of spatial and temporal scales, the changes in terrestrial ecosystem GPP in response to anticipated changes in the Earth's environmental conditions.

Canopy photosynthetic capacity drives contrasting age dynamics of resource use efficiencies between mature temperate evergreen and deciduous forests.

Forest resource use efficiencies (RUEs) can vary with tree age, but the nature of these trends and their underlying mechanisms are not well understood. Understanding the age dynamics of forest RUEs and their drivers is vital for assessing the trade-offs between forest functions and resource consumption, making rational management policy, and projecting ecosystem carbon dynamics. Here we used the FLUXNET2015 and AmeriFlux datasets and published literature to explore the age-dependent variability of forest light use efficiency (LUE) and inherent water use efficiency as well as their main regulatory variables in temperate regions. Our results showed that evergreen forest RUEs initially increased before reaching the mature stage (i.e., around 90 years old), and then gradually declined; in contrast, RUEs continuously increased with age for mature deciduous forests. Changing canopy photosynthetic capacity (Amax) was the primary cause of age-related changes in RUEs across temperate forest sites. More importantly, soil nitrogen (N) increased in mature deciduous forests through time but decreased in older evergreen forests. The age-dependent changes in soil N were closely linked with the age dynamics of Amax for mature temperate forests. Additionally, soil nutrient conditions played a greater role in deciduous forest RUEs than evergreen forest RUEs. This study highlights the importance of stand age and forest type on temperate forest RUEs over the long term.

Solar-induced chlorophyll fluorescence exhibits a universal relationship with gross primary productivity across a wide variety of biomes.

In our recent study in Global Change Biology (Li et al., ), we examined the relationship between solar-induced chlorophyll fluorescence (SIF) measured from the Orbiting Carbon Observatory-2 (OCO-2) and gross primary productivity (GPP) derived from eddy covariance flux towers across the globe, and we discovered that there is a nearly universal relationship between SIF and GPP across a wide variety of biomes. This finding reveals the tremendous potential of SIF for accurately mapping terrestrial photosynthesis globally.

Solar-induced chlorophyll fluorescence is strongly correlated with terrestrial photosynthesis for a wide variety of biomes: First global analysis based on OCO-2 and flux tower observations.

Solar-induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite-observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3 km × 2.25 km) measured from the Orbiting Carbon Observatory-2 (OCO-2) provides the first opportunity to examine the SIF-GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO-2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO-2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757 nm (R2  = 0.72, p < 0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes (R2  = 0.57-0.79, p < 0.0001) except evergreen broadleaf forests (R2  = 0.16, p < 0.05) at the daily timescale. A higher slope was found for C4 grasslands and croplands than for C3 ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome-specific SIF-GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO-2 SIF generally had a better performance for predicting GPP than satellite-derived vegetation indices and a light use efficiency model. The universal SIF-GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO-2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies.

Contrasting ecosystem CO2 fluxes of inland and coastal wetlands: a meta-analysis of eddy covariance data.

Wetlands play an important role in regulating the atmospheric carbon dioxide (CO2 ) concentrations and thus affecting the climate. However, there is still lack of quantitative evaluation of such a role across different wetland types, especially at the global scale. Here, we conducted a meta-analysis to compare ecosystem CO2 fluxes among various types of wetlands using a global database compiled from the literature. This database consists of 143 site-years of eddy covariance data from 22 inland wetland and 21 coastal wetland sites across the globe. Coastal wetlands had higher annual gross primary productivity (GPP), ecosystem respiration (Re ), and net ecosystem productivity (NEP) than inland wetlands. On a per unit area basis, coastal wetlands provided large CO2 sinks, while inland wetlands provided small CO2 sinks or were nearly CO2 neutral. The annual CO2 sink strength was 93.15 and 208.37 g C m-2 for inland and coastal wetlands, respectively. Annual CO2 fluxes were mainly regulated by mean annual temperature (MAT) and mean annual precipitation (MAP). For coastal and inland wetlands combined, MAT and MAP explained 71%, 54%, and 57% of the variations in GPP, Re , and NEP, respectively. The CO2 fluxes of wetlands were also related to leaf area index (LAI). The CO2 fluxes also varied with water table depth (WTD), although the effects of WTD were not statistically significant. NEP was jointly determined by GPP and Re for both inland and coastal wetlands. However, the NEP/Re and NEP/GPP ratios exhibited little variability for inland wetlands and decreased for coastal wetlands with increasing latitude. The contrasting of CO2 fluxes between inland and coastal wetlands globally can improve our understanding of the roles of wetlands in the global C cycle. Our results also have implications for informing wetland management and climate change policymaking, for example, the efforts being made by international organizations and enterprises to restore coastal wetlands for enhancing blue carbon sinks.

Habitat quality evaluation and pattern simulation of coastal salt marsh wetlands.

Coastal salt marsh wetlands not only sequester a large amount of organic carbon, mitigating the effect of climate change, but also nurture rich wetland resources and diverse ecological environments. In this study, habitat pattern and quality of the Jiangsu Yancheng Wetland Rare Birds National Nature Reserve were studied. The evolution of habitat patterns was analyzed using the U-Net model and Sentinel-2 data. The habitat quality was evaluated using the InVEST model, while the future habitat pattern in 2027 under different scenarios were simulated using the PLUS model. Our results showed that, during 2017-2022, the Suaeda salsa habitat showed a net decrease in area of 2077.61 ha, while Spartina alterniflora and Phragmites australis habitats manifested a net increase in different degrees. The overall habitat pattern was characterized by fragmentation decline and regularization enhancement. The habitat quality decreased from 0.75 to 0.72, mainly due to the loss of the S. salsa habitat and the expansion of the P. australis habitat. The simulation results indicated that, the habitat quality is expected to further decline to 0.71 under the natural development scenario, and 390.27 ha of S. salsa habitat will convert to P. australis. While in government control scenario, the habitat quality is expected to improve to 0.78, which was 0.07 higher than that in natural development scenario, and S. salsa habitat can be restored well. This study provides a scientific basis for the protection of suitable habitats for waterfowl and is crucial for the ecological conservation and management planning of nature reserves and coastal salt marsh wetlands.

Global convergence in terrestrial gross primary production response to atmospheric vapor pressure deficit.

Atmospheric vapor pressure deficit (VPD) increases with climate warming and may limit plant growth. However, gross primary production (GPP) responses to VPD remain a mystery, offering a significant source of uncertainty in the estimation of global terrestrial ecosystems carbon dynamics. In this study, in-situ measurements, satellite-derived data, and Earth System Models (ESMs) simulations were analysed to show that the GPP of most ecosystems has a similar threshold in response to VPD: first increasing and then declining. When VPD exceeds these thresholds, atmospheric drought stress reduces soil moisture and stomatal conductance, thereby decreasing the productivity of terrestrial ecosystems. Current ESMs underscore CO2 fertilization effects but predict significant GPP decline in low-latitude ecosystems when VPD exceeds the thresholds. These results emphasize the impacts of climate warming on VPD and propose limitations to future ecosystems productivity caused by increased atmospheric water demand. Incorporating VPD, soil moisture, and canopy conductance interactions into ESMs enhances the prediction of terrestrial ecosystem responses to climate change.

A daily gap-free normalized difference vegetation index dataset from 1981 to 2023 in China.

Long-term, daily, and gap-free Normalized Difference Vegetation Index (NDVI) is of great significance for a better Earth system observation. However, gaps and contamination are quite severe in current daily NDVI datasets. This study developed a daily 0.05° gap-free NDVI dataset from 1981-2023 in China by combining valid data identification and spatiotemporal sequence gap-filling techniques based on the National Oceanic and Atmospheric Administration daily NDVI dataset. The generated NDVI in more than 99.91% of the study area showed an absolute percent bias (|PB|) smaller than 1% compared with the original valid data, with an overall R2 and root mean square error (RMSE) of 0.79 and 0.05, respectively. PB and RMSE between our dataset and the MODIS daily gap-filled NDVI dataset (MCD19A3CMG) during 2000 to 2023 are 7.54% and 0.1, respectively. PB between our dataset and three monthly NDVI datasets (i.e., GIMMS3g, MODIS MOD13C2, and SPOT/PROBA) are only -5.79%, 4.82%, and 2.66%, respectively. To the best of our knowledge, this is the first long-term daily gap-free NDVI in China by far.

Air quality improvements can strengthen China's food security.

Air pollution exerts crucial influence on crop yields and impacts regional and global food supplies. Here we employ a statistical model using satellite-based observations and flexible functional forms to analyse the synergistic effects of reductions in ozone and aerosols on China's food security. The model consistently shows that ozone is detrimental to crops, whereas aerosol has variable effects. China's maize, rice and wheat yields are projected to increase by 7.84%, 4.10% and 3.43%, respectively, upon reaching two air quality targets (60 µg m-3 for peak-season ozone and 35 µg m-3 for annual fine particulate matter). Average calories produced from these crops would surge by 4.51%, potentially allowing China to attain grain self-sufficiency 2 years earlier than previously estimated. These results show that ozone pollution control should be a high priority to increase staple crop edible calories, and future stringent air pollution regulations would enhance China's food security.

The reduced net carbon uptake over Northern Hemisphere land causes the close-to-normal CO2 growth rate in 2021 La Niña.

La Niña climate anomalies have historically been associated with substantial reductions in the atmospheric CO2 growth rate. However, the 2021 La Niña exhibited a unique near-neutral impact on the CO2 growth rate. In this study, we investigate the underlying mechanisms by using an ensemble of net CO2 fluxes constrained by CO2 observations from the Orbiting Carbon Observatory-2 in conjunction with estimates of gross primary production and fire carbon emissions. Our analysis reveals that the close-to-normal atmospheric CO2 growth rate in 2021 was the result of the compensation between increased net carbon uptake over the tropics and reduced net carbon uptake over the Northern Hemisphere mid-latitudes. Specifically, we identify that the extreme drought and warm anomalies in Europe and Asia reduced the net carbon uptake and offset 72% of the increased net carbon uptake over the tropics in 2021. This study contributes to our broader understanding of how regional processes can shape the trajectory of atmospheric CO2 concentration under climate change.

Soil environmental anomalies dominate the responses of net ecosystem productivity to heatwaves in three Mongolian grasslands.

Climate change is causing more frequent and intense heatwaves. Therefore, it is important to understand how heatwaves affect the terrestrial carbon cycle, especially in grasslands, which are especially susceptible to climate extremes. This study assessed the impact of naturally occurring, simultaneous short-term heatwaves on CO2 fluxes in three ecosystems on the Mongolia Plateau: meadow steppe (MDW), typical steppe (TPL), and shrub-grassland (SHB). During three heatwaves, net ecosystem productivity (NEP) was reduced by 86 %, 178 %, and 172 % at MDW, TPL, and SHB, respectively. The changes in ecosystem respiration, gross primary production, evapotranspiration, and water use efficiency were divergent, indicating the mechanisms underlying the observed NEP decreases among the sites. The impact of the heatwave in MDW was mitigated by the high soil water content, which enhanced evapotranspiration and subsequent cooling effects. However, at TPL, insufficient soil water led to combined thermal and drought stress and low resilience. At SHB, the ecosystem's low tolerance to an August heatwave was heavily influenced by species phenology, as it coincided with the key phenological growing phase of plants. The potential key mechanism of divergent NEP response to heatwaves lies in the divergent stability and varying importance of environmental factors, combined with the specific sensitivity of NEP to each factor in ecosystems. Furthermore, our findings suggest that anomalies in soil environment, rather than atmospheric anomalies, are the primary determinants of NEP anomalies during heatwaves. This challenges the conventional understanding of heatwaves as a discrete and ephemeral periods of high air temperatures. Instead, heatwaves should be viewed as chronologically variable, compound, and time-sensitive environmental stressors. The ultimate impact of heatwaves on ecosystems is co-determined by a complex interplay of environmental, biological, and heatwave features.

A national-scale assessment of land subsidence in China's major cities.

China's massive wave of urbanization may be threatened by land subsidence. Using a spaceborne synthetic aperture radar interferometry technique, we provided a systematic assessment of land subsidence in all of China's major cities from 2015 to 2022. Of the examined urban lands, 45% are subsiding faster than 3 millimeters per year, and 16% are subsiding faster than 10 millimeters per year, affecting 29 and 7% of the urban population, respectively. The subsidence appears to be associated with a range of factors such as groundwater withdrawal and the weight of buildings. By 2120, 22 to 26% of China's coastal lands will have a relative elevation lower than sea level, hosting 9 to 11% of the coastal population, because of the combined effect of city subsidence and sea-level rise. Our results underscore the necessity of enhancing protective measures to mitigate potential damages from subsidence.