Sea surface Fresnel reflections difference driven de-glint algorithm for airborne optical images.

This paper presents a glint correction algorithm for high spatial resolution optical remote sensing imagery captured by the ER-2 Airborne Visual Infrared Imaging Spectrometer (AVIRIS). The algorithm employs linear and differential techniques to mitigate sun glint and sky glint effects, encompassing statistical glint reflections resulting from variations in imaging angles within strips and inter-strip variations due to Fresnel reflectance disparities. It aims to diminish Fresnel reflectance diversity on water surfaces and mitigate the distortions induced by glint reflectance during spectral and ocean color inversion. A comparative analysis of spectral and ocean color information in AVIRIS images before and after correction reveals enhanced accuracy following the glint correction. By systematically addressing multiple glint reflections and their ramifications, this method offers a valuable framework for correcting water surface glint in diverse high spatial resolution optical imagery.

China's Fossil Fuel CO2 Emissions Estimated Using Surface Observations of Coemitted NO2.

Accurate estimates of fossil fuel CO2 (FFCO2) emissions are of great importance for climate prediction and mitigation regulations but remain a significant challenge for accounting methods relying on economic statistics and emission factors. In this study, we employed a regional data assimilation framework to assimilate in situ NO2 observations, allowing us to combine observation-constrained NOx emissions coemitted with FFCO2 and grid-specific CO2-to-NOx emission ratios to infer the daily FFCO2 emissions over China. The estimated national total for 2016 was 11.4 PgCO2·yr-1, with an uncertainty (1σ) of 1.5 PgCO2·yr-1 that accounted for errors associated with atmospheric transport, inversion framework parameters, and CO2-to-NOx emission ratios. Our findings indicated that widely used "bottom-up" emission inventories generally ignore numerous activity level statistics of FFCO2 related to energy industries and power plants in western China, whereas the inventories are significantly overestimated in developed regions and key urban areas owing to exaggerated emission factors and inexact spatial disaggregation. The optimized FFCO2 estimate exhibited more distinct seasonality with a significant increase in emissions in winter. These findings advance our understanding of the spatiotemporal regime of FFCO2 emissions in China.

Improving sea surface floating matter identification from Sentinel-2 MSI imagery using optical radiative simulation of neighborhood difference.

The reflectance difference (ΔR) between a floating matter pixel and a nearby water reference pixel is a method of atmospheric radiation unmixing. This technique unveils target signals by referencing the background within the horizontal neighborhood. ΔR is effective for removing the mixed-pixel effect and partial atmospheric path radiance. However, other atmospheric interference sources in the difference pixel, including atmospheric extinction and sunglint, need to be clarified. To address these challenges, we combined in situ floating matter endmember spectra for simulation and Sentinel-2 Multispectral Instrument (MSI) sensors for validation. We focused on radiative transfer simulation of horizontal neighborhood and vertical atmospheric column, investigating the bilateral conversion of ΔR between bottom-of-atmosphere (BOA) and top-of-atmosphere (TOA) signals, and clarifying how the atmosphere affects the difference pixel (ΔR) and floating matter identification. Results showed that direct use of TOA ΔR works in discriminating algae from non-algae floating matters under weak sunglint, and is a suitable candidate for no bother with atmospheric correction, least uncertain, and wider coverage. And then, sunglint interference is also inevitable, whether serious or not.

Sunglint reflection facilitates performance of spaceborne UV sensor in oil spill detection.

Ultraviolet Imager (UVI) onboard Haiyang-1C/D (HY-1C/D) satellites has been providing ultraviolet (UV) data to detect marine oil spills since 2018. Although the scale effect of UV remote sensing has been preliminarily interpreted, the application characteristics of spaceborne UV sensors with medium spatial resolution in oil spill detection deserve further investigation, especially the role of sunglint in the process of detection. In this study, the performance of the UVI is thoroughly assessed by the following aspects: image features of oils under sunglint, sunglint requirement for spaceborne UV detection of oils, and the stability of the UVI signal. The results indicate that in UVI images, it is sunglint reflection that determines the image features of spilled oils, and the appearance of sunglint can strengthen the contrast between oils and seawater. Besides, the required sunglint strength in spaceborne UV detection has been deduced to be 10-3 - 10-4 sr-1, which is higher than that in the VNIR wavelengths. Moreover, uncertainties in the UVI signal can meet the demand to discriminate between oils and seawater. The above results can confirm the capability of the UVI and the critical role of sunglint in spaceborne UV detection of marine oil spills, and provide new reference for spaceborne UV remote sensing.

Correction of multi-scale sunglint reflections from the water surface in airborne high-spatial resolution optical images.

Airborne optical images (AOI) are often with complex sunglint reflections, which brings a certain influence to watercolor retrieval. This includes the sunglint reflection with water surface statistical distribution characteristics caused by imaging viewing angles differences, with high spatial resolution surface discrete characteristics sharing similar viewing angles, and the surface Fresnel reflection sunglint differences caused by the skylight difference during the flight of unmanned aerial vehicles. Aiming at the multiscale optical characteristics of sunglint reflection in high spatial resolution AOI, based on multi-path optical radiation transmission, the sunglint reflection interference from three different imaging processes is clarified. We developed a correction method to eliminate these different sunglint reflections on water surfaces and improve the reflectivity accuracy. The comparison with the in situ measured remote sensing reflectance of water indicated that the root mean square error (RMSE) was reduced from 0.0009 sr-1 to 0.0004 sr-1, and the mean relative error (MRE) decreased from 21.8% to 15.7%. This method has also been applied to correct the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) images, showing good applicability. The method is fast, effective, and without auxiliary parameters, which provides a correction reference for different surface sunglint corrections of various AOI.

Revisiting the cumulative effects of drought on global gross primary productivity based on new long-term series data (1982-2018).

Drought has broad and deep impacts on vegetation. Studies on the effects of drought on vegetation have been conducted over years. Recently, the cumulative effect of drought is recognized as another key factor affecting plant growth. However, global-scale studies on this phenomenon are still lacking. Thus, based on new satellite based gross primary productivity (GPP) and multi-temporal scale Standardized Precipitation Evapotranspiration Index data sets, we explored the cumulative effect duration (CED) of drought on global vegetation GPP and analyzed its variability across elevations and climatic zones. The main findings were as follows: (1) The cumulative effect of drought on GPP was widespread, with an average CED of 4.89 months. (2) CED of drought on GPP varied among vegetation types. Specifically, grasslands showed the longest duration, with an average value of 5.28 months, followed by shrublands (5.09 months), wetlands (5.03 months), croplands (4.85 months), savannas (4.58 months), and forestlands (4.57 months). (3) CED of drought on GPP changes with climate conditions. It decreased with the decrease of precipitation in the driest month (Pdry ) and mean annual precipitation in tropical and arid climate zones, respectively. In both temperate and cold climate zones, CED of drought on GPP was shorter in areas with dry winter than that in areas with dry summer. It increased with the decrease of mean annual air temperature in tropical climate zones and decreased with the increase of summer temperature in temperate and cold climatic zones. (4) With increasing elevation, CED of drought on GPP showed a pattern of increasing (0-3000 m), then decreasing (3000-5000 m), and increasing again (>5000 m). Our findings highlight the heterogeneity of CED of drought on GPP, owing to differences in vegetation types, long-term hydrothermal conditions, elevation, etc. The results could deepen our understanding of the effects of drought on global vegetation.

NO x Emission Changes Over China During the COVID-19 Epidemic Inferred From Surface NO2 Observations.

The COVID-19 epidemic has substantially limited human activities and affected anthropogenic emissions. In this work, daily NO x emissions are inferred using a regional data assimilation system and hourly surface NO2 measurement over China. The results show that because of the coronavirus outbreak, NO x emissions across the whole mainland China dropped sharply after 31 January, began to rise slightly in certain areas after 10 February, and gradually recover across the country after 20 February. Compared with the emissions before the outbreak, NO x emissions fell by more than 60% and ~30% in many large cities and most small to medium cities, respectively. Overall, NO x emissions were reduced by 36% over China, which were mainly contributed by transportation. Evaluations show that the inverted changes over eastern China are credible, whereas those in western China might be underestimated. These findings are of great significance for exploring the reduction potential of NO x emissions in China.

Interannual variation in grassland net ecosystem productivity and its coupling relation to climatic factors in China.

Grassland, as an important part of land cover, plays an important role in the global carbon cycle and carbon balance. Net ecosystem productivity (NEP) is a key indicator of the carbon cycle process and an important factor in assessing ecosystem security and maintaining ecosystem balance. In this paper, Boreal Ecosystem Productivity Simulator (BEPS) combining meteorological data, leaf area index, and land cover type data were used to simulate the grassland NEP of China from 1979 to 2008. This model was also used to analyze the responses to changes in climate factors, interannual variation in carbon conversion efficiency, drought stress coefficient, and water use efficiency of grassland in China. Results showed that from 1979 to 2008, the mean annual grassland NEP was 13.6 g C/m2 with weak carbon sinks. The grassland NEP distribution increased from northwest to southeast across China. Regions with NEP of > 0 (C sink) accounted for 73.1% of the total grassland area of China. The total C sequestration reached 26.6 Tg yearly, and grassland NEP was positive from 1979 to 2008. The annual changing characteristics were analyzed. Grassland NEP was positive with carbon sink from June to September, which was negative with carbon source in the remaining months. The carbon conversion efficiency and water use efficiency of the grassland increased significantly within 30 years. NEP showed positive correlation with precipitation (accounting for 74.2% of the total grassland area was positively correlated) but weakly positive correlation with temperature (50.2% of the case). Furthermore, significant positive correlation was found between grassland NEP and precipitation, especially in northeastern and central Inner Mongolia, northern Tianshan of Xinjiang, southwestern Tibet, and southern Qinghai Lake.

Satellite sun-induced chlorophyll fluorescence detects early response of winter wheat to heat stress in the Indian Indo-Gangetic Plains.

Extremely high temperatures represent one of the most severe abiotic stresses limiting crop productivity. However, understanding crop responses to heat stress is still limited considering the increases in both the frequency and severity of heat wave events under climate change. This limited understanding is partly due to the lack of studies or tools for the timely and accurate monitoring of crop responses to extreme heat over broad spatial scales. In this work, we use novel spaceborne data of sun-induced chlorophyll fluorescence (SIF), which is a new proxy for photosynthetic activity, along with traditional vegetation indices (Normalized Difference Vegetation Index NDVI and Enhanced Vegetation Index EVI) to investigate the impacts of heat stress on winter wheat in northwestern India, one of the world's major wheat production areas. In 2010, an abrupt rise in temperature that began in March adversely affected the productivity of wheat and caused yield losses of 6% compared to previous year. The yield predicted by satellite observations of SIF decreased by approximately 13.9%, compared to the 1.2% and 0.4% changes in NDVI and EVI, respectively. During early stage of this heat wave event in early March 2010, the SIF observations showed a significant reduction and earlier response, while NDVI and EVI showed no changes and could not capture the heat stress until late March. The spatial patterns of SIF anomalies closely tracked the temporal evolution of the heat stress over the study area. Furthermore, our results show that SIF can provide large-scale, physiology-related wheat stress response as indicated by the larger reduction in fluorescence yield (SIFyield ) than fraction of photosynthetically active radiation during the grain-filling phase, which may have eventually led to the reduction in wheat yield in 2010. This study implies that satellite observations of SIF have great potential to detect heat stress conditions in wheat in a timely manner and assess their impacts on wheat yields at large scales.

[A inversion model for remote sensing of leaf water content based on the leaf optical property].

Leaf water content is a fundamental physiological characteristic parameter of crops, and plays an important role in the study of the ecological environment. The aim of the work reported in this paper was to focus upon the retrieval of leaf water content from leaf-scale reflectance spectra by developing a physical inversion model based on the radiative transfer theory and wavelet analysis techniques. A continuous wavelet transform was performed on each of leaf component specific absorption coefficients to pick wavelet coefficients that were identified as highly sensitive to leaf water content and insensitive to other components. In the present study, for identifying the most appropriate wavelet, the six frequently used wavelet functions available within MATLAB were tested. Two biorl. 5 wavelet coefficients observed at the scale of 200 nm are provided with good performance, their wave-length positions are located at 1 405 and 1 488 nm, respectively. Two factors (α and Δ) of the predictive theoretical models based on the biorl. 5 wavelet coefficients of the leaf-scale reflectance spectra were determined by leaf structure parameter N. We built a database composed of thousands of simulated leaf reflectance spectra with the PROSPECT model. The entire dataset was split into two parts, with 60% the calibration subset assigned to calibrating two factors (α and Δ) of the predictive theoretical model. The remaining 40% the validation subset combined with the LOPEX93 experimental dataset used for validating the models. The results showed that the accuracy of the models compare to the statistical regression models derived from the traditional vegetation indices has improved with the highest predictive coefficient of determination (R2) of 0. 987, and the model becomes more robust. This study presented that wavelet analysis has the potential to capture much more of the information contained with reflectance spectra than previous analytical approaches which have tended to focus on using a small number of optimal wavebands while discarding the majority of the spectrum.

Satellite-derived subsurface urban heat island.

The subsurface urban heat island (SubUHI) is one part of the overall UHI specifying the relative warmth of urban ground temperatures against the rural background. To combat the challenge on measuring extensive underground temperatures with in situ instruments, we utilized satellite-based moderate-resolution imaging spectroradiometer data to reconstruct the subsurface thermal field over the Beijing metropolis through a three-time-scale model. The results show the SubUHI's high spatial heterogeneity. Within the depths shallower than 0.5 m, the SubUHI dominates along the depth profiles and analyses imply the moments for the SubUHI intensity reaching first and second extremes during a diurnal temperature cycle are delayed about 3.25 and 1.97 h per 0.1 m, respectively. At depths shallower than 0.05 m in particular, there is a subsurface urban cool island (UCI) in spring daytime, mainly owing to the surface UCI that occurs in this period. At depths between 0.5 and 10 m, the time for the SubUHI intensity getting to its extremes during an annual temperature cycle is lagged 26.2 days per meter. Within these depths, the SubUHI prevails without exception, with an average intensity of 4.3 K, varying from 3.2 to 5.3 K.

The carbon budget of China: 1980-2021.

As one of the world's largest emitters of greenhouse gases, China has set itself the ambitious goal of achieving carbon peaking and carbon neutrality. Therefore, it is crucial to quantify the magnitude and trend of sources and sinks of atmospheric carbon dioxide (CO2), and to monitor China's progress toward these goals. Using state-of-the-art datasets and models, this study comprehensively estimated the anthropogenic CO2 emissions from energy, industrial processes and product use, and waste along with natural sources and sinks of CO2 for all of China during 1980-2021. To recognize the differences among various methods of estimating greenhouse emissions, the estimates are compared with China's National Greenhouse Gas Inventories (NGHGIs) for 1994, 2005, 2010, 2012, and 2014. Anthropogenic CO2 emissions in China have increased by 7.39 times from 1980 to 12.77 Gt CO2 a-1 in 2021. While benefiting from ecological projects (e.g., Three Norths Shelter Forest System Project), the land carbon sink in China has reached 1.65 Gt CO2 a-1 averaged through 2010-2021, which is almost 15.81 times that of the carbon sink in the 1980s. On average, China's terrestrial ecosystems offset 14.69% ± 2.49% of anthropogenic CO2 emissions through 2010-2021. Two provincial-level administrative regions of China, Xizang and Qinghai, have achieved carbon neutrality according to our estimates, but nearly half of the administrative regions of China have terrestrial carbon sink offsets of less than 10% of anthropogenic CO2 emissions. This study indicated a high level of consistency between NGHGIs and various datasets used for estimating fossil CO2 emissions, but found notable differences for land carbon sinks. Future estimates of the terrestrial carbon sinks of NGHGIs urgently need to be verified with process-based models which integrate the comprehensive carbon cycle processes.

Large diurnal compensatory effects mitigate the response of Amazonian forests to atmospheric warming and drying.

Photosynthesis and evapotranspiration in Amazonian forests are major contributors to the global carbon and water cycles. However, their diurnal patterns and responses to atmospheric warming and drying at regional scale remain unclear, hindering the understanding of global carbon and water cycles. Here, we used proxies of photosynthesis and evapotranspiration from the International Space Station to reveal a strong depression of dry season afternoon photosynthesis (by 6.7 ± 2.4%) and evapotranspiration (by 6.1 ± 3.1%). Photosynthesis positively responds to vapor pressure deficit (VPD) in the morning, but negatively in the afternoon. Furthermore, we projected that the regionally depressed afternoon photosynthesis will be compensated by their increases in the morning in future dry seasons. These results shed new light on the complex interplay of climate with carbon and water fluxes in Amazonian forests and provide evidence on the emerging environmental constraints of primary productivity that may improve the robustness of future projections.

Improved estimation of global gross primary productivity during 1981-2020 using the optimized P model.

Accurate estimation of terrestrial gross primary productivity (GPP) is essential for quantifying the net carbon exchange between the atmosphere and biosphere. Light use efficiency (LUE) models are widely used to estimate GPP at different spatial scales. However, difficulties in proper determination of maximum LUE (LUEmax) and downregulation of LUEmax into actual LUE result in uncertainties in GPP estimated by LUE models. The recently developed P model, as a LUE-like model, captures the deep mechanism of photosynthesis and simplifies parameterization. Site level studies have proved the outperformance of P model over LUE models. However, the global application of the P model is still lacking. Thus, the effectiveness of 5 water stress factors integrated into the P model was compared. The optimal P model was used to generate a new long-term (1981-2020) global monthly GPP dataset at a spatial resolution of 0.1° × 0.1°, called PGPP. Validation at globally distributed 109 FLUXNET sites indicated that PGPP is better than three widely-used GPP products. R2 between PGPP and observed GPP equals to 0.75, the corresponding root mean squared error (RMSE) and mean absolute error (MAE) equal to 1.77 g C m-2 d-1 and 1.28 g C m-2 d-1. During the period from 1981 to 2020, PGPP significantly increased in 69.02% of global vegetated regions (p < 0.05). Overall, PGPP provides a new GPP product choice for global ecology studies and the comparison of various water stress factors provides a new idea for the improvement of GPP model in the future.

Anthropogenic emissions estimated using surface observations and their impacts on PM2.5 source apportionment over the Yangtze River Delta, China.

Source-tagged source apportionment (SA) has advantages for quantifying the contribution of various source regions and categories to PM2.5; however, it is highly affected by uncertainty in the emission inventory. In this study, we used a Regional multi-Air Pollutant Assimilation System (RAPAS) to optimize daily SO2, NOx and primary PM2.5 (PPM2.5) emissions in the Yangtze River Delta (YRD) in December 2016 by assimilating hourly in-situ measurements. The CMAQ-ISAM model was implemented with prior and posterior emissions respectively to investigate the impacts of optimizing emissions on PM2.5 SA in the YRD megalopolis (YRDM) and three megacities of Shanghai, Nanjing, and Hangzhou in the YRDM. The results showed that RAPAS significantly improved the simulations and reduced the emission uncertainties of the different pollutants. Compared with prior emissions, the posterior emissions in the YRD decreased by 13% and 11% for SO2 and NOx respectively, and increased by 24% for PPM2.5. Compared with SA using prior emissions, the contributions from Hangzhou, northern Zhejiang, and areas outside of the YRD to the YRDM increased. The local contributions from the YRDM, Nanjing and Shanghai decreased by 1.8%, 9.7%, and 2.3%, respectively, whereas that of Hangzhou increased by 5.6%. The changes in the daily local contributions caused by optimizing emissions ranged from -18.0% to 23.6%. Generally, under stable weather conditions, the local contribution changed the most, whereas under unstable weather conditions, the contribution from upwind areas changed significantly. Overall, with optimized emissions, we found in Nanjing, Shanghai, and Hangzhou, local emissions contributed 18.2%, 39.6% and 36.8% of their PM2.5 concentrations, respectively; long-range transport from outside the YRDM contributed 59.2%, 48.1%, and 48.2%, respectively. This study emphasizes the importance of improving emission estimations for source-tagged SA and provides more reliable SA results for the main cities in the YRD, which will contribute to pollution control in these regions.

A global 0.05° dataset for gross primary production of sunlit and shaded vegetation canopies from 1992 to 2020.

Distinguishing gross primary production of sunlit and shaded leaves (GPPsun and GPPshade) is crucial for improving our understanding of the underlying mechanisms regulating long-term GPP variations. Here we produce a global 0.05°, 8-day dataset for GPP, GPPshade and GPPsun over 1992-2020 using an updated two-leaf light use efficiency model (TL-LUE), which is driven by the GLOBMAP leaf area index, CRUJRA meteorology, and ESA-CCI land cover. Our products estimate the mean annual totals of global GPP, GPPsun, and GPPshade over 1992-2020 at 125.0 ± 3.8 (mean ± std) Pg C a-1, 50.5 ± 1.2 Pg C a-1, and 74.5 ± 2.6 Pg C a-1, respectively, in which EBF (evergreen broadleaf forest) and CRO (crops) contribute more than half of the totals. They show clear increasing trends over time, in which the trend of GPP (also GPPsun and GPPshade) for CRO is distinctively greatest, and that for DBF (deciduous broadleaf forest) is relatively large and GPPshade overwhelmingly outweighs GPPsun. This new dataset advances our in-depth understanding of large-scale carbon cycle processes and dynamics.

The status of carbon neutrality of the world's top 5 CO2 emitters as seen by carbon satellites.

China, the Unite States (US), the European Union (EU), India, and Russia are the world's top 5 fossil fuel and cement CO2 (FFC) emitting countries or regions (CRs). It is very important to understand their status of carbon neutrality, and to monitor their future changes of net carbon fluxes (NCFs). In this study, we implemented a well-established global carbon assimilation system (GCAS, Version 2) to infer global surface carbon fluxes from May 2009 to December 2019 using both GOSAT and OCO-2 XCO2 retrievals. The reductions of flux uncertainty and XCO2 bias, and the evaluation of posterior flux show that GCAS has comparable and good performance in the 5 CRs. The results suggest that Russia has achieved carbon neutrality, but the other 4 are still far from being carbon neutral, especially China. The mean annual NCFs in China, the US, the EU, India, and Russia are 2.33 ± 0.29, 0.82 ± 0.20, 0.42 ± 0.16, 0.50 ± 0.12, and -0.33 ± 0.23 PgC yr-1, respectively. From 2010 to 2019, the NCFs showed an increasing trend in the US and India, a slight downward trend after 2013 in China, and were stable in the EU. The changes of land sinks in China and the US might be the main reason for their trends. India's trend was mainly due to the increase of FFC emission. The relative contributions of NCFs to the global land net carbon emission of China and the EU have decreased, while those of the US and India have increased, implying the US and India must take more active measures to control carbon emissions or increase their sinks. This study indicates that satellite XCO2 could be successfully used to monitor the changes of regional NCFs, which is of great significance for major countries to achieve greenhouse gas control goals.

The Nord Stream pipeline gas leaks released approximately 220,000 tonnes of methane into the atmosphere.

Sudden mega natural gas leaks of two Nord Stream pipelines in the Baltic Sea (Denmark) occurred from late September to early October 2022, releasing large amounts of methane into the atmosphere. We inferred the methane emissions of this event based on surface in situ observations using two inversion methods and two meteorological reanalysis datasets, supplemented with satellite-based observations. We conclude that approximately 220 ± 30 Gg of methane was released from September 26 to October 1, 2022.

Relationships between net primary productivity and stand age for several forest types and their influence on China's carbon balance.

Affected by natural and anthropogenic disturbances such as forest fires, insect-induced mortality and harvesting, forest stand age plays an important role in determining the distribution of carbon pools and fluxes in a variety of forest ecosystems. An improved understanding of the relationship between net primary productivity (NPP) and stand age (i.e., age-related increase and decline in forest productivity) is essential for the simulation and prediction of the global carbon cycle at annual, decadal, centurial, or even longer temporal scales. In this paper, we developed functions describing the relationship between national mean NPP and stand age using stand age information derived from forest inventory data and NPP simulated by the BEPS (Boreal Ecosystem Productivity Simulator) model in 2001. Due to differences in ecobiophysical characteristics of different forest types, NPP-age equations were developed for five typical forest ecosystems in China (deciduous needleleaf forest (DNF), evergreen needleleaf forest in tropic and subtropical zones (ENF-S), deciduous broadleaf forest (DBF), evergreen broadleaf forest (EBF), and mixed broadleaf forest (MBF)). For DNF, ENF-S, EBF, and MBF, changes in NPP with age were well fitted with a common non-linear function, with R(2) values equal to 0.90, 0.75, 0.66, and 0.67, respectively. In contrast, a second order polynomial was best suitable for simulating the change of NPP for DBF, with an R(2) value of 0.79. The timing and magnitude of the maximum NPP varied with forest types. DNF, EBF, and MBF reached the peak NPP at the age of 54, 40, and 32 years, respectively, while the NPP of ENF-S maximizes at the age of 13 years. The highest NPP of DBF appeared at 122 years. NPP was generally lower in older stands with the exception of DBF, and this particular finding runs counter to the paradigm of age-related decline in forest growth. Evaluation based on measurements of NPP and stand age at the plot-level demonstrates the reliability and applicability of the fitted NPP-age relationships. These relationships were used to replace the normalized NPP-age relationship used in the original InTEC (Integrated Terrestrial Ecosystem Carbon) model, to improve the accuracy of estimated carbon balance for China's forest ecosystems. With the revised NPP-age relationship, the InTEC model simulated a larger carbon source from 1950-1980 and a larger carbon sink from 1985-2001 for China's forests than the original InTEC model did because of the modification to the age-related carbon dynamics in forests. This finding confirms the importance of considering the dynamics of NPP related to forest age in estimating regional and global terrestrial carbon budgets.

The effect of water stress on net primary productivity in northwest China.

Net primary productivity (NPP) has been widely used as the indicator of vegetation function and exhibits large spatial and temporal variations caused by numerous factors. Northwest China (NWC) is one of the driest regions in China, and water supply is the key determinant of NPP here. However, studies on the effects of water stress on NPP in NWC at the regional scale are still relatively lacking. Thus, in this study, based on a set of Moderate-Resolution Imaging Spectroradiometer (MODIS) NPP and evapotranspiration (ET) datasets, we quantified the response of NPP to water stress, which is indicated by crop water stress index (CWSI). Regional average of annual NPP in NWC showed an increasing trend during the study period, at a rate of 0.84 g C m-2 yr-1. At the province level, the NPP increase rates increased in the order of Ningxia (7.7%), Shaanxi (6.5%), Gansu (4.5%), Qinghai (3.8%), and Xinjiang (1.7%). NPP was negatively correlated with CWSI (p<0.05) in 73% of areas, indicating the key role of water stress in constraining NPP over this arid region. The effect of water stress on NPP changes with elevation. Water stress has the strongest negative impact on NPP in areas with elevations around 2000 m. In elevations above 5000 m, NPP is not limited by water stress, mostly positively correlated with CWSI. Our findings further clarify the importance of water stress in dryland ecosystems, while highlighting that elevation gradients can significantly affect the correlation between NPP and water stress.