Effects of climate change and ozone on vegetation phenology on the Tibetan Plateau

The vegetation phenology, encompassing the start (SOS) and end (EOS) of the growing season on the Tibetan Plateau, has been significantly impacted by global climate change. Furthermore, ozone (O3) has gradually become the main pollutant in this region, substantially influencing carbon cycle and ecosystems on Earth. While ongoing studies have focused mainly on the implications of climate parameters, including temperature, precipitation, and radiation, the effects of O3 on the SOS and EOS remain unclear. Here, we compared the responses and sensitivities of the SOS and EOS to both climatic factors and O3 in this region. With the use of partial correlation analysis, we found that increased precipitation was the most important factor influencing the SOS and caused earlier occurrence (4.8 % vs. 21.9 %) for most plant functional types. In comparison, temperature only dominated in shrublands. In particular, we found that the EOS responded comparably to climatic factors with similar proportions between advancing and delaying patterns. However, higher O3 levels consistently advanced the EOS for almost all plant functional types and was the main factor controlling EOS variations based on the sensitivity analysis. Our results emphasized that O3 pollution should be considered for obtaining better phenological forecasts and determining the impacts of the environment and atmospheric composition on carbon sequestration in terrestrial ecosystems.

Parsimonious estimation of hourly surface ozone concentration across China during 2015-2020

Surface ozone is an important air pollutant detrimental to human health and vegetation productivity, particularly in China. However, high resolution surface ozone concentration data is still lacking, largely hindering accurate assessment of associated environmental impacts. Here, we collected hourly ground ozone observations (over 6 million records), remote sensing products, meteorological data, and social-economic information, and applied recurrent neural networks to map hourly surface ozone data (HrSOD) at a 0.1° × 0.1° resolution across China during 2015-2020. The coefficient of determination (R2) values in sample-based, site-based, and by-year cross-validations were 0.72, 0.65 and 0.71, respectively, with the root mean square error (RMSE) values being 11.71 ppb (mean = 30.89 ppb), 12.81 ppb (mean = 30.96 ppb) and 11.14 ppb (mean = 31.26 ppb). Moreover, it exhibits high spatiotemporal consistency with ground-level observations at different time scales (diurnal, seasonal, annual), and at various spatial levels (individual sites and regional scales). Meanwhile, the HrSOD provides critical information for fine-resolution assessment of surface ozone impacts on environmental and human benefits.

A single-point modeling approach for the intercomparison and evaluation of ozone dry deposition across chemical transport models (Activity 2 of AQMEII4)

A primary sink of air pollutants and their precursors is dry deposition. Dry deposition estimates differ across chemical transport models, yet an understanding of the model spread is incomplete. Here, we introduce Activity 2 of the Air Quality Model Evaluation International Initiative Phase 4 (AQMEII4). We examine 18 dry deposition schemes from regional and global chemical transport models as well as standalone models used for impact assessments or process understanding. We configure the schemes as single-point models at eight Northern Hemisphere locations with observed ozone fluxes. Single-point models are driven by a common set of site-specific meteorological and environmental conditions. Five of eight sites have at least 3 years and up to 12 years of ozone fluxes. The interquartile range across models in multiyear mean ozone deposition velocities ranges from a factor of 1.2 to 1.9 annually across sites and tends to be highest during winter compared with summer. No model is within 50 % of observed multiyear averages across all sites and seasons, but some models perform well for some sites and seasons. For the first time, we demonstrate how contributions from depositional pathways vary across models. Models can disagree with respect to relative contributions from the pathways, even when they predict similar deposition velocities, or agree with respect to the relative contributions but predict different deposition velocities. Both stomatal and nonstomatal uptake contribute to the large model spread across sites. Our findings are the beginning of results from AQMEII4 Activity 2, which brings scientists who model air quality and dry deposition together with scientists who measure ozone fluxes to evaluate and improve dry deposition schemes in the chemical transport models used for research, planning, and regulatory purposes.

Legislative and functional aspects of different metrics used for ozone risk assessment to forests

Surface ozone (O₃) is a threat to forests by decreasing photosynthesis and, consequently, influencing the strength of land carbon sink. However, due to the lack of continuous surface O₃ measurements, observational-based assessments of O₃ impacts on forests are largely missing at hemispheric to global scales. Currently, some metrics are used for regulatory purposes by governments or national agencies to protect forests against the negative impacts of ozone: in particular, both Europe and United States (US) makes use of two different exposure-based metrics, i.e. AOT40 and W126, respectively. However, because of some limitations in these metrics, a new standard is under consideration by the European Union (EU) to replace the current exposure metric. We analyse here the different air quality standards set or proposed for use in Europe and in the US to protect forests from O₃ and to evaluate their spatial and temporal consistency while assessing their effectiveness in protecting northern-hemisphere forests. Then, we compare their results with the information obtained from a complex land surface model (ORCHIDEE). We find that present O₃ uptake decreases gross primary production (GPP) in 37.7% of the NH forested area of northern hemisphere with a mean loss of 2.4% year^-1. We show how the proposed US (W126) and the currently used European (AOT40) air quality standards substantially overestimate the extension of potential vulnerable regions, predicting that 46% and 61% of the Northern Hemisphere (NH) forested area are at risk of O₃ pollution. Conversely, the new proposed European standard (POD1) identifies lower extension of vulnerability regions (39.6%).

The Influence of Ozone on Net Ecosystem Production of a Ryegrass–Clover Mixture under Field Conditions

In order to understand the effect of phytotoxic tropospheric ozone (O3) on terrestrial vegetation, we quantified the impact of current O3 concentration ([O3]) on net ecosystem production (NEP) when compared to the conditions of the pre-industrial era. We compared and tested linear mixed-effects models based on [O3] and stomatal O3 flux (Fsto). The managed ryegrass–clover (Lolium perenne and Trifolium pratense) mixture was grown on arable land in the Czech Republic, Central Europe. Values of [O3] and Fsto were measured and calculated based on resistance analogy, respectively, while NEP was calculated from eddy covariance CO2 fluxes. We found the Fsto-based model more precise when compared to measured NEP. High Fsto was found even at low [O3], while broad summer maximum of [O3] was not necessarily followed by significant NEP decline, due to low soil water content leading to a low stomatal conductivity and Fsto. Comparing to low pre-industrial O3 conditions, current levels of O3 resulted in the reduction of cumulative NEP over the entire growing season, up to 29.7 and 13.5% when the [O3]-based and Fsto-based model was applied, respectively. During the growing season, an O3-induced reduction of NEP ranged between 13.1% in May and 26.2% in July when compared to pre-industrial Fsto levels. Looking to the future, high [O3] and Fsto may lead to the reduction of current NEP by approximately 13.3% on average during the growing season, but may increase by up to 61–86.6% in autumn, indicating further O3-induced acceleration of the senescence. These findings indicate the importance of Fsto and its inclusion into the models estimating O3 effects on terrestrial vegetation. The interaction between environmental factors and stomatal conductance is therefore discussed in detail.

Increased greenhouse gas emissions intensity of major croplands in China: Implications for food security and climate change mitigation

Balancing crop production and greenhouse gas (GHG) emissions from agriculture soil requires a better understanding and quantification of crop GHG emissions intensity, a measure of GHG emissions per unit crop production. Here we conduct a state-of-the-art estimate of the spatial-temporal variability of GHG emissions intensities for wheat, maize, and rice in China from 1949 to 2012 using an improved agricultural ecosystem model (Dynamic Land Ecosystem Model-Agriculture Version 2.0) and meta-analysis covering 172 field-GHG emissions experiments. The results show that the GHG emissions intensities of these croplands from 1949 to 2012, on average, were 0.10-1.31 kg CO₂ -eq/kg, with a significant increase rate of 1.84-3.58 × 10^-3  kg CO₂ -eq kg^-1  year^-1 . Nitrogen fertilizer was the dominant factor contributing to the increase in GHG emissions intensity in northern China and increased its impact in southern China in the 2000s. Increasing GHG emissions intensity implies that excessive fertilizer failed to markedly stimulate crop yield increase in China but still exacerbated soil GHG emissions. This study found that overfertilization of more than 60% was mainly located in the winter wheat-summer maize rotation systems in the North China Plain, the winter wheat-rice rotation systems in the middle and lower reaches of the Yangtze River and southwest China, and most of the double rice systems in the South. Our simulations suggest that roughly a one-third reduction in the current N fertilizer application level over these "overfertilization" regions would not significantly influence crop yield but decrease soil GHG emissions by 29.60%-32.50% and GHG emissions intensity by 0.13-0.25 kg CO₂ -eq/kg. This reduction is about 29% and 5% of total agricultural soil GHG emissions in China and the world, respectively. This study suggests that improving nitrogen use efficiency would be an effective strategy to mitigate GHG emissions and sustain China's food security.

Effects of ozone on agriculture, forests and grasslands

The damage and injury that ground level ozone (O3) causes vegetation has become increasingly evident over the past half century with a large body of observational and experimental evidence demonstrating a variety of effects at ambient concentrations on crop, forest and grassland species and ecosystems. This paper explores the use of experimental data to develop exposure-response relationships for use in risk assessment studies. These studies have typically identified the USA mid-West, much of Europe, the Indo Gangetic Plain in South Asia and the Eastern coastal region of China as global regions where O3 is likely to threaten food supply and other ecosystems. Global risk assessment modelling estimates yield losses of staple crops between 3 to 16% causing economic losses of between US$14 to 26 billion in the year 2000. Changes in anthropogenic emissions of O3 precursors in recent decades have modified O3 concentration profiles (peaks versus background O3) and global distributions with the Northern Hemisphere seeing increases in O3 levels of between 1 and 5 ppb/decade since the 1950s and the emergence of Asia as the region with the highest O3 concentrations. In the future, O3 mitigation could focus on methane (CH4) and nitrogen oxide (NOx) emissions; these will differentially influence global and local/regional O3 concentrations and influence daily and seasonal profiles. The consequent effects on vegetation will in part depend on how these changes in O3 profile alter the exceedance of detoxification thresholds for plant damage. Adaptation options may play an important role in enhancing food supply while mitigation strategies are being implemented. An improved understanding of the mechanisms by which O3 affects plants, and how this might influence detoxification thresholds and interactions with other environmental variables such as water stress and nutrients, would help develop O3 deposition and impact models to support the development of crop, land-surface exchange and ultimately earth system models for holistic assessments of global change. This article is part of a discussion meeting issue 'Air quality, past present and future'.

Country-Scale Analysis of Methane Emissions with a High-Resolution Inverse Model Using GOSAT and Surface Observations

We employed a global high-resolution inverse model to optimize the CH4 emission using Greenhouse gas Observing Satellite (GOSAT) and surface observation data for a period from 2011–2017 for the two main source categories of anthropogenic and natural emissions. We used the Emission Database for Global Atmospheric Research (EDGAR v4.3.2) for anthropogenic methane emission and scaled them by country to match the national inventories reported to the United Nations Framework Convention on Climate Change (UNFCCC). Wetland and soil sink prior fluxes were simulated using the Vegetation Integrative Simulator of Trace gases (VISIT) model. Biomass burning prior fluxes were provided by the Global Fire Assimilation System (GFAS). We estimated a global total anthropogenic and natural methane emissions of 340.9 Tg CH4 yr−1 and 232.5 Tg CH4 yr−1, respectively. Country-scale analysis of the estimated anthropogenic emissions showed that all the top-emitting countries showed differences with their respective inventories to be within the uncertainty range of the inventories, confirming that the posterior anthropogenic emissions did not deviate from nationally reported values. Large countries, such as China, Russia, and the United States, had the mean estimated emission of 45.7 ± 8.6, 31.9 ± 7.8, and 29.8 ± 7.8 Tg CH4 yr−1, respectively. For natural wetland emissions, we estimated large emissions for Brazil (39.8 ± 12.4 Tg CH4 yr−1), the United States (25.9 ± 8.3 Tg CH4 yr−1), Russia (13.2 ± 9.3 Tg CH4 yr−1), India (12.3 ± 6.4 Tg CH4 yr−1), and Canada (12.2 ± 5.1 Tg CH4 yr−1). In both emission categories, the major emitting countries all had the model corrections to emissions within the uncertainty range of inventories. The advantages of the approach used in this study were: (1) use of high-resolution transport, useful for simulations near emission hotspots, (2) prior anthropogenic emissions adjusted to the UNFCCC reports, (3) combining surface and satellite observations, which improves the estimation of both natural and anthropogenic methane emissions over spatial scale of countries.

Southward propagation of Nazca subduction along the Andes

The Andean margin is the plate-tectonic paradigm for long-lived, continuous subduction, yet its geology since the late Mesozoic era (the past 100 million years or so) has been far from steady state. The episodic deformation and magmatism have been attributed to cyclic changes in the dip angle of the subducting slab, slab break-off and the penetration of the slab into the lower mantle; the role of plate tectonics remains unclear, owing to the extensive subduction of the Nazca-Farallon plate (which has resulted in more than 5,500 kilometres of lithosphere being lost to the mantle). Here, using tomographic data, we recreate the plate-tectonic geometry of the subducted Nazca slab, which enables us to reconstruct Andean plate tectonics since the late Mesozoic. Our model suggests that the current phase of Nazca subduction began at the northern Andes (5° S) during the late Cretaceous period (around 80 million years ago) and propagated southwards, reaching the southern Andes (40° S) by the early Cenozoic era (around 55 million year ago). Thus, contrary to the current paradigm, Nazca subduction has not been fully continuous since the Mesozoic but instead included episodic divergent phases. In addition, we find that foredeep sedimentation and the initiation of Andean compression are both linked to interactions between the Nazca slab and the lower mantle, consistent with previous modelling.

The role of plant phenology in stomatal ozone flux modeling

Plant phenology plays a pivotal role in the climate system as it regulates the gas exchange between the biosphere and the atmosphere. The uptake of ozone by forest is estimated through several meteorological variables and a specific function describing the beginning and the termination of plant growing season; actually, in many risk assessment studies, this function is based on a simple latitude and topography model. In this study, using two satellite datasets, we apply and compare six methods to estimate the start and the end dates of the growing season across a large region covering all Europe for the year 2011. Results show a large variability between the green-up and dormancy dates estimated using the six different methods, with differences greater than one month. However, interestingly, all the methods display a common spatial pattern in the uptake of ozone by forests with a marked change in the magnitude, up to 1.9 TgO₃ /year, and corresponding to a difference of 25% in the amount of ozone that enters the leaves. Our results indicate that improved estimates of ozone fluxes require a better representation of plant phenology in the models used for O₃ risk assessment.

Comparing concentration-based (AOT40) and stomatal uptake (PODY) metrics for ozone risk assessment to European forests

Tropospheric ozone (O3) produces harmful effects to forests and crops, leading to a reduction of land carbon assimilation that, consequently, influences the land sink and the crop yield production. To assess the potential negative O3 impacts to vegetation, the European Union uses the Accumulated Ozone over Threshold of 40 ppb (AOT40). This index has been chosen for its simplicity and flexibility in handling different ecosystems as well as for its linear relationships with yield or biomass loss. However, AOT40 does not give any information on the physiological O3 uptake into the leaves since it does not include any environmental constraints to O3 uptake through stomata. Therefore, an index based on stomatal O3 uptake (i.e. PODY), which describes the amount of O3 entering into the leaves, would be more appropriate. Specifically, the PODY metric considers the effects of multiple climatic factors, vegetation characteristics and local and phenological inputs rather than the only atmospheric O3 concentration. For this reason, the use of PODY in the O3 risk assessment for vegetation is becoming recommended. We compare different potential O3 risk assessments based on two methodologies (i.e. AOT40 and stomatal O3 uptake) using a framework of mesoscale models that produces hourly meteorological and O3 data at high spatial resolution (12 km) over Europe for the time period 2000-2005. Results indicate a remarkable spatial and temporal inconsistency between the two indices, suggesting that a new definition of European legislative standard is needed in the near future. Besides, our risk assessment based on AOT40 shows a good consistency compared to both in-situ data and other model-based datasets. Conversely, risk assessment based on stomatal O3 uptake shows different spatial patterns compared to other model-based datasets. This strong inconsistency can be likely related to a different vegetation cover and its associated parameterizations.

Complex spatiotemporal responses of global terrestrial primary production to climate change and increasing atmospheric CO2 in the 21st century

Quantitative information on the response of global terrestrial net primary production (NPP) to climate change and increasing atmospheric CO₂ is essential for climate change adaptation and mitigation in the 21^st century. Using a process-based ecosystem model (the Dynamic Land Ecosystem Model, DLEM), we quantified the magnitude and spatiotemporal variations of contemporary (2000s) global NPP, and projected its potential responses to climate and CO₂ changes in the 21^st century under the Special Report on Emission Scenarios (SRES) A2 and B1 of Intergovernmental Panel on Climate Change (IPCC). We estimated a global terrestrial NPP of 54.6 (52.8–56.4) PgC yr⁻¹ as a result of multiple factors during 2000–2009. Climate change would either reduce global NPP (4.6%) under the A2 scenario or slightly enhance NPP (2.2%) under the B1 scenario during 2010–2099. In response to climate change, global NPP would first increase until surface air temperature increases by 1.5°C (until the 2030s) and then level-off or decline after it increases by more than 1.5°C (after the 2030s). This result supports the Copenhagen Accord Acknowledgement, which states that staying below 2°C may not be sufficient and the need to potentially aim for staying below 1.5°C. The CO₂ fertilization effect would result in a 12%–13.9% increase in global NPP during the 21^st century. The relative CO₂ fertilization effect, i.e. change in NPP on per CO₂ (ppm) bases, is projected to first increase quickly then level off in the 2070s and even decline by the end of the 2080s, possibly due to CO₂ saturation and nutrient limitation. Terrestrial NPP responses to climate change and elevated atmospheric CO₂ largely varied among biomes, with the largest increases in the tundra and boreal needleleaf deciduous forest. Compared to the low emission scenario (B1), the high emission scenario (A2) would lead to larger spatiotemporal variations in NPP, and more dramatic and counteracting impacts from climate and increasing atmospheric CO₂.

Elevated carbon dioxide and ozone alter productivity and ecosystem carbon content in northern temperate forests

Three young northern temperate forest communities in the north-central United States were exposed to factorial combinations of elevated carbon dioxide (CO2 ) and tropospheric ozone (O3 ) for 11 years. Here, we report results from an extensive sampling of plant biomass and soil conducted at the conclusion of the experiment that enabled us to estimate ecosystem carbon (C) content and cumulative net primary productivity (NPP). Elevated CO2 enhanced ecosystem C content by 11%, whereas elevated O3 decreased ecosystem C content by 9%. There was little variation in treatment effects on C content across communities and no meaningful interactions between CO2 and O3 . Treatment effects on ecosystem C content resulted primarily from changes in the near-surface mineral soil and tree C, particularly differences in woody tissues. Excluding the mineral soil, cumulative NPP was a strong predictor of ecosystem C content (r(2) = 0.96). Elevated CO2 enhanced cumulative NPP by 39%, a consequence of a 28% increase in canopy nitrogen (N) content (g N m(-2) ) and a 28% increase in N productivity (NPP/canopy N). In contrast, elevated O3 lowered NPP by 10% because of a 21% decrease in canopy N, but did not impact N productivity. Consequently, as the marginal impact of canopy N on NPP (∆NPP/∆N) decreased through time with further canopy development, the O3 effect on NPP dissipated. Within the mineral soil, there was less C in the top 0.1 m of soil under elevated O3 and less soil C from 0.1 to 0.2 m in depth under elevated CO2 . Overall, these results suggest that elevated CO2 may create a sustained increase in NPP, whereas the long-term effect of elevated O3 on NPP will be smaller than expected. However, changes in soil C are not well-understood and limit our ability to predict changes in ecosystem C content.

Climate and land use controls on soil organic carbon in the loess plateau region of China

The Loess Plateau of China has the highest soil erosion rate in the world where billion tons of soil is annually washed into Yellow River. In recent decades this region has experienced significant climate change and policy-driven land conversion. However, it has not yet been well investigated how these changes in climate and land use have affected soil organic carbon (SOC) storage on the Loess Plateau. By using the Dynamic Land Ecosystem Model (DLEM), we quantified the effects of climate and land use on SOC storage on the Loess Plateau in the context of multiple environmental factors during the period of 1961–2005. Our results show that SOC storage increased by 0.27 Pg C on the Loess Plateau as a result of multiple environmental factors during the study period. About 55% (0.14 Pg C) of the SOC increase was caused by land conversion from cropland to grassland/forest owing to the government efforts to reduce soil erosion and improve the ecological conditions in the region. Historical climate change reduced SOC by 0.05 Pg C (approximately 19% of the total change) primarily due to a significant climate warming and a slight reduction in precipitation. Our results imply that the implementation of “Grain for Green” policy may effectively enhance regional soil carbon storage and hence starve off further soil erosion on the Loess Plateau.

Net greenhouse gas balance in response to nitrogen enrichment: perspectives from a coupled biogeochemical model

Increasing reactive nitrogen (N) input has been recognized as one of the important factors influencing climate system through affecting the uptake and emission of greenhouse gases (GHG). However, the magnitude and spatiotemporal variations of N-induced GHG fluxes at regional and global scales remain far from certain. Here we selected China as an example, and used a coupled biogeochemical model in conjunction with spatially explicit data sets (including climate, atmospheric CO2 , O3 , N deposition, land use, and land cover changes, and N fertilizer application) to simulate the concurrent impacts of increasing atmospheric and fertilized N inputs on balance of three major GHGs (CO2 , CH4 , and N2 O). Our simulations showed that these two N enrichment sources in China decreased global warming potential (GWP) through stimulating CO2 sink and suppressing CH4 emission. However, direct N2 O emission was estimated to offset 39% of N-induced carbon (C) benefit, with a net GWP of three GHGs averaging -376.3 ± 146.4 Tg CO2  eq yr(-1) (the standard deviation is interannual variability of GWP) during 2000-2008. The chemical N fertilizer uses were estimated to increase GWP by 45.6 ± 34.3 Tg CO2  eq yr(-1) in the same period, and C sink was offset by 136%. The largest C sink offset ratio due to increasing N input was found in Southeast and Central mainland of China, where rapid industrial development and intensively managed crop system are located. Although exposed to the rapidly increasing N deposition, most of the natural vegetation covers were still showing decreasing GWP. However, due to extensive overuse of N fertilizer, China's cropland was found to show the least negative GWP, or even positive GWP in recent decade. From both scientific and policy perspectives, it is essential to incorporate multiple GHGs into a coupled biogeochemical framework for fully assessing N impacts on climate changes.

Effects of elevated O₃ concentration on winter wheat and rice yields in the Yangtze River Delta, China

The effects of a continuing rise of ambient ozone on crop yield will seriously threaten food security in China. In the Yangtze River Delta, a rapidly developing and seriously air polluted region in China, innovative open-top chambers have been established to fumigate winter wheat and rice in situ with elevated O(3). Five years of study have shown that the yields of wheat and rice decreased with increasing O(3) concentration. There were significant relationships between the relative yield and AOT40 (accumulated hourly O(3) concentration over 40 ppb) for both winter wheat and rice. Winter wheat was more sensitive to O(3) than rice. O(3)-induced yield declines were attributed primarily to 1000-grain weight and harvest index for winter wheat, and attributed primarily to grain number per panicle and harvest index for rice. Control of ambient O(3) pollution and breeding of O(3) tolerant crops are urgent to guarantee food security in China.

China's crop productivity and soil carbon storage as influenced by multifactor global change

Much concern has been raised about how multifactor global change has affected food security and carbon sequestration capacity in China. By using a process-based ecosystem model, the Dynamic Land Ecosystem Model (DLEM), in conjunction with the newly developed driving information on multiple environmental factors (climate, atmospheric CO2 , tropospheric ozone, nitrogen deposition, and land cover/land use change), we quantified spatial and temporal patterns of net primary production (NPP) and soil organic carbon storage (SOC) across China's croplands during 1980-2005 and investigated the underlying mechanisms. Simulated results showed that both crop NPP and SOC increased from 1980 to 2005, and the highest annual NPP occurred in the Southeast (SE) region (0.32 Pg C yr(-1) , 35.4% of the total NPP) whereas the largest annual SOC (2.29 Pg C yr(-1) , 35.4% of the total SOC) was found in the Northeast (NE) region. Land management practices, particularly nitrogen fertilizer application, appear to be the most important factor in stimulating increase in NPP and SOC. However, tropospheric ozone pollution and climate change led to NPP reduction and SOC loss. Our results suggest that China's crop productivity and soil carbon storage could be enhanced through minimizing tropospheric ozone pollution and improving nitrogen fertilizer use efficiency.

Impacts of urbanization on carbon balance in terrestrial ecosystems of the Southern United States

Using a process-based Dynamic Land Ecosystem Model, we assessed carbon dynamics of urbanized/developed lands in the Southern United States during 1945-2007. The results indicated that approximately 1.72 (1.69-1.77) Pg (1P = 10(15)) carbon was stored in urban/developed lands, comparable to the storage of shrubland or cropland in the region. Urbanization resulted in a release of 0.21 Pg carbon to the atmosphere during 1945-2007. Pre-urbanization vegetation type and time since land conversion were two primary factors determining the extent of urbanization impacts on carbon dynamics. After a rapid decline of carbon storage during land conversion, an urban ecosystem gradually accumulates carbon and may compensate for the initial carbon loss in 70-100 years. The carbon sequestration rate of urban ecosystem diminishes with time, nearly disappearing in two centuries after land conversion. This study implied that it is important to take urbanization effect into account for assessing regional carbon balance.

Effect of nitrogen deposition on China's terrestrial carbon uptake in the context of multifactor environmental changes

The amount of atmospheric nitrogen (N) deposited on the land surface has increased globally and by nearly five times in China from 1901 to 2005. Little is known about how elevated reactive N input has affected the carbon (C) sequestration capability of China's terrestrial ecosystems, largely due to the lack of reliable data on N deposition. Here we have used a newly developed data set of historical N deposition at a spatial resolution of 10 km x 10 km in combination with other gridded historical information on climate, atmospheric composition, land use, and land management practices to drive a process-based ecosystem model, the dynamic land ecosystem model (DLEM) for examining how increasing N deposition and its interactions with other environmental changes have affected C fluxes and storage in China's terrestrial ecosystems during 1901-2005. Our model simulations indicate that increased N deposition has resulted in a net C sink of 62 Tg C/yr (1 Tg = 1012 g) in China's terrestrial ecosystems, totaling up to 6.51 Pg C (1 Pg = 10(15) g) in the past 105 years. During the study period, the N-induced C sequestration can compensate for more than 25% of fossil-fuel CO2 emission from China. The largest C sink was found in southeast China, a region that experienced the most significant increase of N deposition in the period 1901-2005. However, the net primary productivity induced by per-unit N deposition (referred to as ecosystem N use efficiency, ENUE, in this paper) has leveled off or declined since the 1980s. This indicates that part of the deposited N may not be invested to stimulate plant growth, but instead leave the ecosystem by various pathways. Except shrubland and northwest/southwest China, signs of N saturation are apparent in the rest major biome types and regions, with ENUE peaking in the 1980s and leveling off or declining thereafter. Therefore, to minimize the excessive N pollution while keeping the N-stimulated C uptake in China's terrestrial ecosystems, optimized management practices should be taken to increase N use efficiency rather than to keep raising N input level in the near future.