Fertilizer management for global ammonia emission reduction.

Crop production is a large source of atmospheric ammonia (NH3), which poses risks to air quality, human health and ecosystems1-5. However, estimating global NH3 emissions from croplands is subject to uncertainties because of data limitations, thereby limiting the accurate identification of mitigation options and efficacy4,5. Here we develop a machine learning model for generating crop-specific and spatially explicit NH3 emission factors globally (5-arcmin resolution) based on a compiled dataset of field observations. We show that global NH3 emissions from rice, wheat and maize fields in 2018 were 4.3 ± 1.0 Tg N yr-1, lower than previous estimates that did not fully consider fertilizer management practices6-9. Furthermore, spatially optimizing fertilizer management, as guided by the machine learning model, has the potential to reduce the NH3 emissions by about 38% (1.6 ± 0.4 Tg N yr-1) without altering total fertilizer nitrogen inputs. Specifically, we estimate potential NH3 emissions reductions of 47% (44-56%) for rice, 27% (24-28%) for maize and 26% (20-28%) for wheat cultivation, respectively. Under future climate change scenarios, we estimate that NH3 emissions could increase by 4.0 ± 2.7% under SSP1-2.6 and 5.5 ± 5.7% under SSP5-8.5 by 2030-2060. However, targeted fertilizer management has the potential to mitigate these increases.

Climate change increases carbon allocation to leaves in early leaf green-up.

Global greening, characterized by an increase in leaf area index (LAI), implies an increase in foliar carbon (C). Whether this increase in foliar C under climate change is due to higher photosynthesis or to higher allocation of C to leaves remains unknown. Here, we explored the trends in foliar C accumulation and allocation during leaf green-up from 2000 to 2017 using satellite-derived LAI and solar-induced chlorophyll fluorescence (SIF) across the Northern Hemisphere. The accumulation of foliar C accelerated in the early green-up period due to both increased photosynthesis and higher foliar C allocation driven by climate change. In the late stage of green-up, however, we detected decreasing trends in foliar C accumulation and foliar C allocation. Such stage-dependent trends in the accumulation and allocation of foliar C are not represented in current terrestrial biosphere models. Our results highlight that a better representation of C allocation should be incorporated into models.

Urban warming advances spring phenology but reduces the response of phenology to temperature in the conterminous United States.

Urbanization has caused environmental changes, such as urban heat islands (UHIs), that affect terrestrial ecosystems. However, how and to what extent urbanization affects plant phenology remains relatively unexplored. Here, we investigated the changes in the satellite-derived start of season (SOS) and the covariation between SOS and temperature (RT ) in 85 large cities across the conterminous United States for the period 2001-2014. We found that 1) the SOS came significantly earlier (6.1 ± 6.3 d) in 74 cities and RT was significantly weaker (0.03 ± 0.07) in 43 cities when compared with their surrounding rural areas (P < 0.05); 2) the decreased magnitude in RT mainly occurred in cities in relatively cold regions with an annual mean temperature <17.3 °C (e.g., Minnesota, Michigan, and Pennsylvania); and 3) the magnitude of urban-rural difference in both SOS and RT was primarily correlated with the intensity of UHI. Simulations of two phenology models further suggested that more and faster heat accumulation contributed to the earlier SOS, while a decrease in required chilling led to a decline in RT magnitude in urban areas. These findings provide observational evidence of a reduced covariation between temperature and SOS in major US cities, implying the response of spring phenology to warming conditions in nonurban environments may decline in the warming future.

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.

Unusual characteristics of the carbon cycle during the 2015-2016 El Niño.

The 2015-2016 El Niño was one of the strongest on record, but its influence on the carbon balance is less clear. Using Northern Hemisphere atmospheric CO2 observations, we found both detrended atmospheric CO2 growth rate (CGR) and CO2 seasonal-cycle amplitude (SCA) of 2015-2016 were much higher than that of other El Niño events. The simultaneous high CGR and SCA were unusual, because our analysis of long-term CO2 observations at Mauna Loa revealed a significantly negative correlation between CGR and SCA. Atmospheric inversions and terrestrial ecosystem models indicate strong northern land carbon uptake during spring but substantially reduced carbon uptake (or high emissions) during early autumn, which amplified SCA but also resulted in a small anomaly in annual carbon uptake of northern ecosystems in 2015-2016. This negative ecosystem carbon uptake anomaly in early autumn was primarily due to soil water deficits and more litter decomposition caused by enhanced spring productivity. Our study demonstrates a decoupling between seasonality and annual carbon cycle balance in northern ecosystems over 2015-2016, which is unprecedented in the past five decades of El Niño events.

Photoperiod decelerates the advance of spring phenology of six deciduous tree species under climate warming.

Vegetation phenology in spring has substantially advanced under climate warming, consequently shifting the seasonality of ecosystem process and altering biosphere-atmosphere feedbacks. However, whether and to what extent photoperiod (i.e., daylength) affects the phenological advancement is unclear, leading to large uncertainties in projecting future phenological changes. Here we examined the photoperiod effect on spring phenology at a regional scale using in situ observation of six deciduous tree species from the Pan European Phenological Network during 1980-2016. We disentangled the photoperiod effect from the temperature effect (i.e., forcing and chilling) by utilizing the unique topography of the northern Alps of Europe (i.e., varying daylength but uniform temperature distribution across latitudes) and examining phenological changes across latitudes. We found prominent photoperiod-induced shifts in spring leaf-out across latitudes (up to 1.7 days per latitudinal degree). Photoperiod regulates spring phenology by delaying early leaf-out and advancing late leaf-out caused by temperature variations. Based on these findings, we proposed two phenological models that consider the photoperiod effect through different mechanisms and compared them with a chilling model. We found that photoperiod regulation would slow down the advance in spring leaf-out under projected climate warming and thus mitigate the increasing frost risk in spring that deciduous forests will face in the future. Our findings identify photoperiod as a critical but understudied factor influencing spring phenology, suggesting that the responses of terrestrial ecosystem processes to climate warming are likely to be overestimated without adequately considering the photoperiod effect.

Global vegetation biomass production efficiency constrained by models and observations.

Plants use only a fraction of their photosynthetically derived carbon for biomass production (BP). The biomass production efficiency (BPE), defined as the ratio of BP to photosynthesis, and its variation across and within vegetation types is poorly understood, which hinders our capacity to accurately estimate carbon turnover times and carbon sinks. Here, we present a new global estimation of BPE obtained by combining field measurements from 113 sites with 14 carbon cycle models. Our best estimate of global BPE is 0.41 ± 0.05, excluding cropland. The largest BPE is found in boreal forests (0.48 ± 0.06) and the lowest in tropical forests (0.40 ± 0.04). Carbon cycle models overestimate BPE, although models with carbon-nitrogen interactions tend to be more realistic. Using observation-based estimates of global photosynthesis, we quantify the global BP of non-cropland ecosystems of 41 ± 6 Pg C/year. This flux is less than net primary production as it does not contain carbon allocated to symbionts, used for exudates or volatile carbon compound emissions to the atmosphere. Our study reveals a positive bias of 24 ± 11% in the model-estimated BP (10 of 14 models). When correcting models for this bias while leaving modeled carbon turnover times unchanged, we found that the global ecosystem carbon storage change during the last century is decreased by 67% (or 58 Pg C).

Cryptic phenology in plants: Case studies, implications, and recommendations.

Plant phenology-the timing of cyclic or recurrent biological events in plants-offers insight into the ecology, evolution, and seasonality of plant-mediated ecosystem processes. Traditionally studied phenologies are readily apparent, such as flowering events, germination timing, and season-initiating budbreak. However, a broad range of phenologies that are fundamental to the ecology and evolution of plants, and to global biogeochemical cycles and climate change predictions, have been neglected because they are "cryptic"-that is, hidden from view (e.g., root production) or difficult to distinguish and interpret based on common measurements at typical scales of examination (e.g., leaf turnover in evergreen forests). We illustrate how capturing cryptic phenology can advance scientific understanding with two case studies: wood phenology in a deciduous forest of the northeastern USA and leaf phenology in tropical evergreen forests of Amazonia. Drawing on these case studies and other literature, we argue that conceptualizing and characterizing cryptic plant phenology is needed for understanding and accurate prediction at many scales from organisms to ecosystems. We recommend avenues of empirical and modeling research to accelerate discovery of cryptic phenological patterns, to understand their causes and consequences, and to represent these processes in terrestrial biosphere models.

Response of vegetation phenology to urbanization in the conterminous United States.

The influence of urbanization on vegetation phenology is gaining considerable attention due to its implications for human health, cycling of carbon and other nutrients in Earth system. In this study, we examined the relationship between change in vegetation phenology and urban size, an indicator of urbanization, for the conterminous United States. We studied more than 4500 urban clusters of varying size to determine the impact of urbanization on plant phenology, with the aids of remotely sensed observations since 2003-2012. We found that phenology cycle (changes in vegetation greenness) in urban areas starts earlier (start of season, SOS) and ends later (end of season, EOS), resulting in a longer growing season length (GSL), when compared to the respective surrounding urban areas. The average difference of GSL between urban and rural areas over all vegetation types, considered in this study, is about 9 days. Also, the extended GSL in urban area is consistent among different climate zones in the United States, whereas their magnitudes are varying across regions. We found that a tenfold increase in urban size could result in an earlier SOS of about 1.3 days and a later EOS of around 2.4 days. As a result, the GSL could be extended by approximately 3.6 days with a range of 1.6-6.5 days for 25th ~ 75th quantiles, with a median value of about 2.1 days. For different vegetation types, the phenology response to urbanization, as defined by GSL, ranges from 1 to 4 days. The quantitative relationship between phenology and urbanization is of great use for developing improved models of vegetation phenology dynamics under future urbanization, and for developing change indicators to assess the impacts of urbanization on vegetation phenology.

Informing climate models with rapid chamber measurements of forest carbon uptake.

Models predicting ecosystem carbon dioxide (CO2 ) exchange under future climate change rely on relatively few real-world tests of their assumptions and outputs. Here, we demonstrate a rapid and cost-effective method to estimate CO2 exchange from intact vegetation patches under varying atmospheric CO2 concentrations. We find that net ecosystem CO2 uptake (NEE) in a boreal forest rose linearly by 4.7 ± 0.2% of the current ambient rate for every 10 ppm CO2 increase, with no detectable influence of foliar biomass, season, or nitrogen (N) fertilization. The lack of any clear short-term NEE response to fertilization in such an N-limited system is inconsistent with the instantaneous downregulation of photosynthesis formalized in many global models. Incorporating an alternative mechanism with considerable empirical support - diversion of excess carbon to storage compounds - into an existing earth system model brings the model output into closer agreement with our field measurements. A global simulation incorporating this modified model reduces a long-standing mismatch between the modeled and observed seasonal amplitude of atmospheric CO2 . Wider application of this chamber approach would provide critical data needed to further improve modeled projections of biosphere-atmosphere CO2 exchange in a changing climate.

Seasonal responses of terrestrial ecosystem water-use efficiency to climate change.

Ecosystem water-use efficiency (EWUE) is an indicator of carbon-water interactions and is defined as the ratio of carbon assimilation (GPP) to evapotranspiration (ET). Previous research suggests an increasing long-term trend in annual EWUE over many regions and is largely attributed to the physiological effects of rising CO2 . The seasonal trends in EWUE, however, have not yet been analyzed. In this study, we investigate seasonal EWUE trends and responses to various drivers during 1982-2008. The seasonal cycle for two variants of EWUE, water-use efficiency (WUE, GPP/ET), and transpiration-based WUE (WUEt , the ratio of GPP and transpiration), is analyzed from 0.5° gridded fields from four process-based models and satellite-based products, as well as a network of 63 local flux tower observations. WUE derived from flux tower observations shows moderate seasonal variation for most latitude bands, which is in agreement with satellite-based products. In contrast, the seasonal EWUE trends are not well captured by the same satellite-based products. Trend analysis, based on process-model factorial simulations separating effects of climate, CO2 , and nitrogen deposition (NDEP), further suggests that the seasonal EWUE trends are mainly associated with seasonal trends of climate, whereas CO2 and NDEP do not show obvious seasonal difference in EWUE trends. About 66% grid cells show positive annual WUE trends, mainly over mid- and high northern latitudes. In these regions, spring climate change has amplified the effect of CO2 in increasing WUE by more than 0.005 gC m(-2)  mm(-1)  yr(-1) for 41% pixels. Multiple regression analysis further shows that the increase in springtime WUE in the northern hemisphere is the result of GPP increasing faster than ET because of the higher temperature sensitivity of GPP relative to ET. The partitioning of annual EWUE to seasonal components provides new insight into the relative sensitivities of GPP and ET to climate, CO2, and NDEP.

Seasonally different response of photosynthetic activity to daytime and night-time warming in the Northern Hemisphere.

Over the last century the Northern Hemisphere has experienced rapid climate warming, but this warming has not been evenly distributed seasonally, as well as diurnally. The implications of such seasonal and diurnal heterogeneous warming on regional and global vegetation photosynthetic activity, however, are still poorly understood. Here, we investigated for different seasons how photosynthetic activity of vegetation correlates with changes in seasonal daytime and night-time temperature across the Northern Hemisphere (>30°N), using Normalized Difference Vegetation Index (NDVI) data from 1982 to 2011 obtained from the Advanced Very High Resolution Radiometer (AVHRR). Our analysis revealed some striking seasonal differences in the response of NDVI to changes in day- vs. night-time temperatures. For instance, while higher daytime temperature (Tmax) is generally associated with higher NDVI values across the boreal zone, the area exhibiting a statistically significant positive correlation between Tmax and NDVI is much larger in spring (41% of area in boreal zone--total area 12.6×10(6) km2) than in summer and autumn (14% and 9%, respectively). In contrast to the predominantly positive response of boreal ecosystems to changes in Tmax, increases in Tmax tended to negatively influence vegetation growth in temperate dry regions, particularly during summer. Changes in night-time temperature (Tmin) correlated negatively with autumnal NDVI in most of the Northern Hemisphere, but had a positive effect on spring and summer NDVI in most temperate regions (e.g., Central North America and Central Asia). Such divergent covariance between the photosynthetic activity of Northern Hemispheric vegetation and day- and night-time temperature changes among different seasons and climate zones suggests a changing dominance of ecophysiological processes across time and space. Understanding the seasonally different responses of vegetation photosynthetic activity to diurnal temperature changes, which have not been captured by current land surface models, is important for improving the performance of next generation regional and global coupled vegetation-climate models.

Change in terrestrial ecosystem water-use efficiency over the last three decades.

Defined as the ratio between gross primary productivity (GPP) and evapotranspiration (ET), ecosystem-scale water-use efficiency (EWUE) is an indicator of the adjustment of vegetation photosynthesis to water loss. The processes controlling EWUE are complex and reflect both a slow evolution of plants and plant communities as well as fast adjustments of ecosystem functioning to changes of limiting resources. In this study, we investigated EWUE trends from 1982 to 2008 using data-driven models derived from satellite observations and process-oriented carbon cycle models. Our findings suggest positive EWUE trends of 0.0056, 0.0007 and 0.0001 g C m(-2)  mm(-1)  yr(-1) under the single effect of rising CO2 ('CO2 '), climate change ('CLIM') and nitrogen deposition ('NDEP'), respectively. Global patterns of EWUE trends under different scenarios suggest that (i) EWUE-CO2 shows global increases, (ii) EWUE-CLIM increases in mainly high latitudes and decreases at middle and low latitudes, (iii) EWUE-NDEP displays slight increasing trends except in west Siberia, eastern Europe, parts of North America and central Amazonia. The data-driven MTE model, however, shows a slight decline of EWUE during the same period (-0.0005 g C m(-2)  mm(-1)  yr(-1) ), which differs from process-model (0.0064 g C m(-2)  mm(-1)  yr(-1) ) simulations with all drivers taken into account. We attribute this discrepancy to the fact that the nonmodeled physiological effects of elevated CO2 reducing stomatal conductance and transpiration (TR) in the MTE model. Partial correlation analysis between EWUE and climate drivers shows similar responses to climatic variables with the data-driven model and the process-oriented models across different ecosystems. Change in water-use efficiency defined from transpiration-based WUEt (GPP/TR) and inherent water-use efficiency (IWUEt , GPP×VPD/TR) in response to rising CO2 , climate change, and nitrogen deposition are also discussed. Our analyses will facilitate mechanistic understanding of the carbon-water interactions over terrestrial ecosystems under global change.

Detection and attribution of vegetation greening trend in China over the last 30 years.

The reliable detection and attribution of changes in vegetation growth is a prerequisite for the development of strategies for the sustainable management of ecosystems. This is an extraordinary challenge. To our knowledge, this study is the first to comprehensively detect and attribute a greening trend in China over the last three decades. We use three different satellite-derived Leaf Area Index (LAI) datasets for detection as well as five different process-based ecosystem models for attribution. Rising atmospheric CO2 concentration and nitrogen deposition are identified as the most likely causes of the greening trend in China, explaining 85% and 41% of the average growing-season LAI trend (LAIGS) estimated by satellite datasets (average trend of 0.0070 yr(-1), ranging from 0.0035 yr(-1) to 0.0127 yr(-1)), respectively. The contribution of nitrogen deposition is more clearly seen in southern China than in the north of the country. Models disagree about the contribution of climate change alone to the trend in LAIGS at the country scale (one model shows a significant increasing trend, whereas two others show significant decreasing trends). However, the models generally agree on the negative impacts of climate change in north China and Inner Mongolia and the positive impact in the Qinghai-Xizang plateau. Provincial forest area change tends to be significantly correlated with the trend of LAIGS (P < 0.05), and marginally significantly (P = 0.07) correlated with the residual of LAIGS trend, calculated as the trend observed by satellite minus that estimated by models through considering the effects of climate change, rising CO2 concentration and nitrogen deposition, across different provinces. This result highlights the important role of China's afforestation program in explaining the spatial patterns of trend in vegetation growth.

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models: Current status and future directions.

Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and C loss from soil accounts for a large proportion of land-atmosphere C exchange. Therefore, a small change in soil organic C (SOC) can affect atmospheric carbon dioxide (CO2) concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil C exchange with the atmosphere through site measurements, inventories, and empirical/process-based modeling. However, these estimates are highly uncertain, and identifying major driving forces controlling soil C dynamics remains a key research challenge. This study has compiled century-long (1901-2010) estimates of SOC storage and heterotrophic respiration (Rh) from 10 terrestrial biosphere models (TBMs) in the Multi-scale Synthesis and Terrestrial Model Intercomparison Project and two observation-based data sets. The 10 TBM ensemble shows that global SOC estimate ranges from 425 to 2111 Pg C (1 Pg = 1015 g) with a median value of 1158 Pg C in 2010. The models estimate a broad range of Rh from 35 to 69 Pg C yr-1 with a median value of 51 Pg C yr-1 during 2001-2010. The largest uncertainty in SOC stocks exists in the 40-65°N latitude whereas the largest cross-model divergence in Rh are in the tropics. The modeled SOC change during 1901-2010 ranges from -70 Pg C to 86 Pg C, but in some models the SOC change has a different sign from the change of total C stock, implying very different contribution of vegetation and soil pools in determining the terrestrial C budget among models. The model ensemble-estimated mean residence time of SOC shows a reduction of 3.4 years over the past century, which accelerate C cycling through the land biosphere. All the models agreed that climate and land use changes decreased SOC stocks, while elevated atmospheric CO2 and nitrogen deposition over intact ecosystems increased SOC stocks-even though the responses varied significantly among models. Model representations of temperature and moisture sensitivity, nutrient limitation, and land use partially explain the divergent estimates of global SOC stocks and soil C fluxes in this study. In addition, a major source of systematic error in model estimations relates to nonmodeled SOC storage in wetlands and peatlands, as well as to old C storage in deep soil layers.

Toward verifying fossil fuel CO2 emissions with the CMAQ model: motivation, model description and initial simulation.

UNLABELLED: Motivated by the question of whether and how a state-of-the-art regional chemical transport model (CTM) can facilitate characterization of CO2 spatiotemporal variability and verify CO2 fossil-fuel emissions, we for the first time applied the Community Multiscale Air Quality (CMAQ) model to simulate CO2. This paper presents methods, input data, and initial results for CO2 simulation using CMAQ over the contiguous United States in October 2007. Modeling experiments have been performed to understand the roles of fossil-fuel emissions, biosphere-atmosphere exchange, and meteorology in regulating the spatial distribution of CO2 near the surface over the contiguous United States. Three sets of net ecosystem exchange (NEE) fluxes were used as input to assess the impact of uncertainty of NEE on CO2 concentrations simulated by CMAQ. Observational data from six tall tower sites across the country were used to evaluate model performance. In particular, at the Boulder Atmospheric Observatory (BAO), a tall tower site that receives urban emissions from Denver CO, the CMAQ model using hourly varying, high-resolution CO2 fossil-fuel emissions from the Vulcan inventory and Carbon Tracker optimized NEE reproduced the observed diurnal profile of CO2 reasonably well but with a low bias in the early morning. The spatial distribution of CO2 was found to correlate with NO(x), SO2, and CO, because of their similar fossil-fuel emission sources and common transport processes. These initial results from CMAQ demonstrate the potential of using a regional CTM to help interpret CO2 observations and understand CO2 variability in space and time. The ability to simulate a full suite of air pollutants in CMAQ will also facilitate investigations of their use as tracers for CO2 source attribution. This work serves as a proof of concept and the foundation for more comprehensive examinations of CO2 spatiotemporal variability and various uncertainties in the future. IMPLICATIONS: Atmospheric CO2 has long been modeled and studied on continental to global scales to understand the global carbon cycle. This work demonstrates the potential of modeling and studying CO2 variability at fine spatiotemporal scales with CMAQ, which has been applied extensively, to study traditionally regulated air pollutants. The abundant observational records of these air pollutants and successful experience in studying and reducing their emissions may be useful for verifying CO2 emissions. Although there remains much more to further investigate, this work opens up a discussion on whether and how to study CO2 as an air pollutant.

Consistent time allocation fraction to vegetation green-up versus senescence across northern ecosystems despite recent climate change.

Extended growing season lengths under climatic warming suggest increased time for plant growth. However, research has focused on climatic impacts to the timing or duration of distinct phenological events. Comparatively little is known about impacts to the relative time allocation to distinct phenological events, for example, the proportion of time dedicated to leaf growth versus senescence. We use multiple satellite and ground-based observations to show that, despite recent climate change during 2001 to 2020, the ratio of time allocated to vegetation green-up over senescence has remained stable [1.27 (± 0.92)] across more than 83% of northern ecosystems. This stability is independent of changes in growing season lengths and is caused by widespread positive relationships among vegetation phenological events; longer vegetation green-up results in longer vegetation senescence. These empirical observations were also partly reproduced by 13 dynamic global vegetation models. Our work demonstrates an intrinsic biotic control to vegetation phenology that could explain the timing of vegetation senescence under climate change.

Thermal, water, and land cover factors led to contrasting urban and rural vegetation resilience to extreme hot months.

With continuing global warming and urbanization, it is increasingly important to understand the resilience of urban vegetation to extreme high temperatures, but few studies have examined urban vegetation at large scale or both concurrent and delayed responses. In this study, we performed an urban-rural comparison using the Enhanced Vegetation Index and months that exceed the historical 90th percentile in mean temperature (referred to as "hot months") across 85 major cities in the contiguous United States. We found that hot months initially enhanced vegetation greenness but could cause a decline afterwards, especially for persistent (≥4 months) and intense (≥+2 °C) episodes in summer. The urban responses were more positive than rural in the western United States or in winter, but more negative during spring-autumn in the eastern United States. The east-west difference can be attributed to the higher optimal growth temperatures and lower water stress levels of the western urban vegetation than the rural. The urban responses also had smaller magnitudes than the rural responses, especially in deciduous forest biomes, and least in evergreen forest biomes. Within each biome, analysis at 1 km pixel level showed that impervious fraction and vegetation cover, local urban heat island intensity, and water stress were the key drivers of urban-rural differences. These findings advance our understanding of how prolonged exposure to warm extremes, particularly within urban environments, affects vegetation greenness and vitality. Urban planners and ecosystem managers should prioritize the long and intense events and the key drivers in fostering urban vegetation resilience to heat waves.

A comprehensive analysis of the crop effect on the urban-rural differences in land surface phenology.

The response of land surface phenology (LSP) to the urban heat island effect (UHI) is a useful biological indicator for understanding how vegetated ecosystems will be affected by future climate warming. However, vegetation cover in rural areas is often dominated by cultivated land, whose phenological timing is considerably influenced by agricultural managements (e.g., timing of sowing and harvesting), leading to biased conclusions derived from the urban-rural LSP differences. To demonstrate this problem, we investigated the crop influence on the phenological response to a warmer environment resulting from the UHI effect. We partitioned cities in the United States into cultivated and non-cultivated categories according to the proportion of crops in rural areas. We then built continuous buffer zones starting from the urban boundary to explore the urban-rural LSP differences considering the UHI effect on them. The results suggest crop inclusion is likely to lead to >14 days of urban-rural differences at both the start of the season (SOS) and the end of the season (EOS) between cultivated and non-cultivated cities. The temperature sensitivity (ST) of SOS is overestimated by approximately 2.7 days/°C, whereas the EOS is underestimated by 3.6 days/°C. Removing crop-dominated pixels (i.e., above 50 %) can minimize the influence of crop planting/harvesting on LSP and derive reliable results. We, therefore, suggest explicit consideration of crop impacts in future studies of phenological differences between urban and rural areas and the UHI effect on LSP in urban domains, as presented by this comprehensive study.