Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling vegetation change”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Variability in Surface BRDF at Different Spatial Scales (30 m-500 m) Over a Mixed Agricultural Landscape as Retrieved from Airborne and Satellite Spectral Measurements

Over the past decade, the role of multiangle remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities represented by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained approx.1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75 off-nadir, and at spatial resolutions ranging from 3 m - 500 m. This unique dataset was used to examine the interaction of the spatial and angular characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertain ties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure.

Roman, Miguel O.↗
NASA NTRS2012

Variability in Surface BRDF at Different Spatial Scales (30m-500m) Over a Mixed Agricultural Landscape as Retrieved from Airborne and Satellite Spectral Measurements

Over the past decade, the role of multiangle 1 remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities represented by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained approx.1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75deg off-nadir, and at spatial resolutions ranging from 3 m - 500 m. This unique dataset was used to examine the interaction of the spatial and angular 18 characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertainties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure.

Roman, Miguel O.↗
NASA NTRS1992

Botswana water and surface energy balance research program. Part 2: Large scale moisture and passive microwaves

The Botswana water and surface energy balance research program was developed to study and evaluate the integrated use of multispectral satellite remote sensing for monitoring the hydrological status of the Earth's surface. The research program consisted of two major, mutually related components: a surface energy balance modeling component, built around an extensive field campaign; and a passive microwave research component which consisted of a retrospective study of large scale moisture conditions and Nimbus scanning multichannel microwave radiometer microwave signatures. The integrated approach of both components are explained in general and activities performed within the passive microwave research component are summarized. The microwave theory is discussed taking into account: soil dielectric constant, emissivity, soil roughness effects, vegetation effects, optical depth, single scattering albedo, and wavelength effects. The study site is described. The soil moisture data and its processing are considered. The relation between observed large scale soil moisture and normalized brightness temperatures is discussed. Vegetation characteristics and inverse modeling of soil emissivity is considered.

Vandegriend, A. A.↗
DOE OSTI2025

Cooling outweighs warming across phenological transitions in the Northern Hemisphere

Vegetation phenology, i.e., seasonal biological events such as leaf-out and leaf-fall, regulates local climate through biophysical processes like evapotranspiration (ET) and albedo. However, the net surface temperature impact of these processes—whether ET cooling or albedo-induced warming predominates—and how the dominance changes across phenological transitions and regions remains poorly understood. Here, we investigated the effects of vegetation foliage on daytime land surface temperature (LST) following six phenological transitions, spanning from the start of season to end of season, in deciduous and mixed forests across the mid- to high-latitude Northern Hemisphere during 2013–2021 using multiple satellite products and ground observations. We quantified vegetation effect as the difference between observed LST and LST estimates from the Annual Temperature Cycle (ATC) model, representing a no-foliage scenario. We found that vegetation-induced cooling consistently outweighs warming following all phenological transitions except for the end of the season. Cooling intensity increased with vegetation greenness, ranging from 1.0 ± 0.5 °C (mean ± 0.15 SD) in 59% of forests after the start of the season (SOS) to 6.1 ± 0.8 °C in 89% of forests following the onset of maturity, before declining toward the end of the season. Over half of the regions experiencing cooling showed intensification of surface cooling with climate warming, suggesting an amplified vegetation-mediated cooling under future climate change. The findings provide a more precise understanding of the role of vegetation in modulating climate at the intraseasonal scale, highlighting the importance of integrating phenological impacts into climate adaptation strategies and Earth system modeling.

albedo↗
DOE OSTI2020

Crop Modeling Application to Improve Irrigation Efficiency in Year-Round Vegetable Production in the Texas Winter Garden Region

Given a rising demand for quality assurance, rather than solely yield, supplemental irrigation plays an important role to ensure the viability and profitability of vegetable crops from unpredictable changes in weather. However, under drought conditions, agricultural irrigation is often given low priority for water allocation. This reduced water availability for agriculture calls for techniques with greater irrigation efficiency, that do not compromise crop quality and yield, and that provide economic benefit for producers. This study developed vegetable growing models for eight different vegetable crops (bush bean, green bean, cabbage, peppermint, spearmint, yellow straight neck squash, zucchini, and bell pepper) based on data from several years of field research. The ALMANAC model accurately simulated yields and water use efficiency (WUE) of all eight vegetables. The developed vegetable models were used to evaluate the effects of various irrigation regimes on vegetable growth and production in several locations in the Winter Garden Region of Texas, under variable weather conditions. Based on our simulation results from 960 scenarios, optimal irrigation amounts that produce high yield as well as reasonable economic profit to producers were determined for each vegetable crop. Overall, yields for all vegetables increased as irrigation amounts increased. However, irrigation amounts did not have a sustainable impact on vegetable yield at high irrigation treatments, and the WUEs of most vegetables were not significantly different among various irrigation regimes. When vegetable yields were compared with water cost, the rate decreased as irrigation amounts increased. Thus, producers will not receive economic benefits when vegetable irrigation water demand is too high.

Kim, Sumin↗
DOE OSTI2023

Machine learning models inaccurately predict current and future high-latitude C balances

The high-latitude carbon (C) cycle is a key feedback to the global climate system, yet because of system complexity and data limitations, there is currently disagreement over whether the region is a source or sink of C. Recent advances in big data analytics and computing power have popularized the use of machine learning (ML) algorithms to upscale site measurements of ecosystem processes, and in some cases forecast the response of these processes to climate change. Due to data limitations, however, ML model predictions of these processes are almost never validated with independent datasets. To better understand and characterize the limitations of these methods, we develop an approach to independently evaluate ML upscaling and forecasting. We mimic data-driven upscaling and forecasting efforts by applying ML algorithms to different subsets of regional process-model simulation gridcells, and then test ML performance using the remaining gridcells. In this study, we simulate C fluxes and environmental data across Alaska using ecosys, a process-rich terrestrial ecosystem model, and then apply boosted regression tree ML algorithms to training data configurations that mirror and expand upon existing AmeriFLUX eddy-covariance data availability. We first show that a ML model trained using ecosys outputs from currently-available Alaska AmeriFLUX sites incorrectly predicts that Alaska is presently a modeled net C source. Increased spatial coverage of the training dataset improves ML predictions, halving the bias when 240 modeled sites are used instead of 15. However, even this more accurate ML model incorrectly predicts Alaska C fluxes under 21st century climate change because of changes in atmospheric CO 2 , litter inputs, and vegetation composition that have impacts on C fluxes which cannot be inferred from the training data. Our results provide key insights to future C flux upscaling efforts and expose the potential for inaccurate ML upscaling and forecasting of high-latitude C cycle dynamics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

The prospects for detecting spectral shifts due to satellite sensor aging

Along with responsivity changes due to sensor aging, there may be concurrent spectral changes. A field-measurement approach for detecting postlaunch spectral changes is described. For illustration, one hypothetical model of change (spectral band shift) is explored through simulation for five satellite sensors. Two different types of natural terrain - vegetation and bare soil - are used as test targets.

Suits, G.↗
NASA NTRS2008

Generating Vegetation Leaf Area Index Earth System Data Record from Multiple Sensors: Theory - Part 1

The generation of multi-decade long Earth System Data Records (ESDRs) of Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation absorbed by vegetation (FPAR) from remote sensing measurements of multiple sensors is key to monitoring long-term changes in vegetation due to natural and anthropogenic influences. Challenges in developing such ESDRs include problems in remote sensing science (modeling of variability in global vegetation, scaling, atmospheric correction) and sensor hardware (differences in spatial resolution, spectral bands, calibration, and information content). In this paper, we develop a physically based approach for deriving LAI and FPAR products from the Advanced Very High Resolution Radiometer (AVHRR) data that are of comparable quality to the Moderate resolution Imaging Spectroradiometer (MODIS) LAI and FPAR products, thus realizing the objective of producing a long (multi-decadal) time series of these products. The approach is based on the radiative transfer theory of canopy spectral invariants which facilitates parameterization of the canopy spectral bidirectional reflectance factor (BRF). The methodology permits decoupling of the structural and radiometric components and obeys the energy conservation law. The approach is applicable to any optical sensor, however, it requires selection of sensor-specific values of configurable parameters, namely, the single scattering albedo and data uncertainty. According to the theory of spectral invariants, the single scattering albedo is a function of the spatial scale, and thus, accounts for the variation in BRF with sensor spatial resolution. Likewise, the single scattering albedo accounts for the variation in spectral BRF with sensor bandwidths. The second adjustable parameter is data uncertainty, which accounts for varying information content of the remote sensing measurements, i.e., Normalized Difference Vegetation Index (NDVI, low information content), vs. spectral BRF (higher information content). Implementation of this approach indicates good consistency in LAI values retrieved from NDVI (AVHRRmode) and spectral BRF (MODIS-mode). Specific details of the implementation and evaluation of the derived products are detailed in the second part of this two-paper series.

Ganguly, Sangram↗
DOE OSTI2020

Pervasive shifts in forest dynamics in a changing world

Forest dynamics arise from the interplay of environmental drivers and disturbances with the demographic processes of recruitment, growth, and mortality, subsequently driving biomass and species composition. However, forest disturbances and subsequent recovery are shifting with global changes in climate and land use, altering these dynamics. Changes in environmental drivers, land use, and disturbance regimes are forcing forests toward younger, shorter stands. Rising carbon dioxide, acclimation, adaptation, and migration can influence these impacts. Recent developments in Earth system models support increasingly realistic simulations of vegetation dynamics. In parallel, emerging remote sensing datasets promise qualitatively new and more abundant data on the underlying processes and consequences for vegetation structure. In combination, these advances hold promise for improving the scientific understanding of changes in vegetation demographics and disturbances.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Dynamics of Aboveground Phytomass of the Circumpolar Arctic Tundra During the Past Three Decades

Numerous studies have evaluated the dynamics of Arctic tundra vegetation throughout the past few decades, using remotely sensed proxies of vegetation, such as the normalized difference vegetation index (NDVI). While extremely useful, these coarse-scale satellite-derived measurements give us minimal information with regard to how these changes are being expressed on the ground, in terms of tundra structure and function. In this analysis, we used a strong regression model between NDVI and aboveground tundra phytomass, developed from extensive field-harvested measurements of vegetation biomass, to estimate the biomass dynamics of the circumpolar Arctic tundra over the period of continuous satellite records (1982-2010). We found that the southernmost tundra subzones (C-E) dominate the increases in biomass, ranging from 20 to 26%, although there was a high degree of heterogeneity across regions, floristic provinces, and vegetation types. The estimated increase in carbon of the aboveground live vegetation of 0.40 Pg C over the past three decades is substantial, although quite small relative to anthropogenic C emissions. However, a 19.8% average increase in aboveground biomass has major implications for nearly all aspects of tundra ecosystems including hydrology, active layer depths, permafrost regimes, wildlife and human use of Arctic landscapes. While spatially extensive on-the-ground measurements of tundra biomass were conducted in the development of this analysis, validation is still impossible without more repeated, long-term monitoring of Arctic tundra biomass in the field.

Arctic tundra↗
NASA NTRS1982

Variation of directional reflectance factors with structural changes of a developing alfalfa canopy

Directional reflectance factors of an alfalfa canopy were determined and related to canopy structure, agronomic variables, and irradiance conditions at four periods during a cutting cycle. Nadir and off-nadir reflectance factors decreased with increasing biomass in Thematic Mapper band 3(0.63-0.69 micrometer) and increased with increasing biomass in band 4(0.76-0.90 micrometer). The sensor view angle had less impact on perceived reflectance as the alfalfa progressed from an erectophile canopy of stems after harvest to a near planophile canopy of leaves at maturity. Studies of directional reflectance are needed for testing and upgrading vegetation canopy models and to aid in the complex interpretation problems presented by aircraft scanners and pointable satellites where illumination and viewing geometries may vary widely. Distinct changes in the patterns of radiance observed by a sensor as structural and biomass changes occur are keys to monitoring the growth and condition of crops.

Kirchner, J. A.↗
NASA NTRS1991

Relationships of soil, grass, and bedrock over the Kaweah serpentine melange through spectral mixture analysis of AVIRIS data

A linear mixing model is used to model the spectral variability of an AVIRIS scene from the western foothills of the Sierra Nevada and calibrate these radiance data to reflectance. Five spectral endmembers from the AVIRIS data, plus an ideal 'shade' endmember were required to model the continuum reflectance of each pixel in the image. Three of the endmembers were interpreted to model the surface constituents green vegetation, dry grass, and illumination. These are the main transient surface constituents that are expected to change with shifts in land use or climatic influences and viewing conditions ('shade' only). The spectral distinction between the other three endmembers is very small, yet the spatial distributions are coherent and interpretable. These distributions cross anthropogenic and vegetation boundaries and are best interpreted as different soil types. Comparison of the fraction images to the bedrock geology maps indicates that substrate composition must be a factor contributing to the spectral properties of these endmembers. Detailed examination of the reflectance spectra of the three soil endmembers reveals that differences in the amount of ferric and ferrous iron and/or organic constituents in the soils is largely responsible for the differences in spectral properties of these endmembers.

Mustard, John F.↗
NASA NTRSDate not supplied

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale for post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of fire impacts on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify wildfire impacts on runoff generation. We consider the vegetation and soil disturbances as the primary contributors to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty and consider model sensitivity, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We evaluated the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration. Analysis of two events in central California shows that 1) the LAI assimilation has a greater impact on water balance; and 2) the soil scheme captures a range of outcomes, with substantial impacts on surface runoff for some ensemble members, which could affect flood potential.

Timothy M Lahmers↗
DOE OSTI2022

4D electrical resistivity tomography for assessing the influence of vegetation and subsurface moisture on railway cutting condition

Instability of slopes, embankments, and cuttings on the railway network is increasingly prevalent globally. Monitoring vulnerable infrastructure aids in geotechnical asset management, and improvements to transport safety and efficiency. Here, we examine the use of a novel, near-real-time Electrical Resistivity Tomography (ERT) monitoring system for assessing the stability of a railway cutting in Leicestershire, United Kingdom. In 2015, an ERT monitoring system was installed across a relict landslide (grassed) and an area of more stable ground on either side (wooded), to monitor changes in electrical resistivity through time and space, and to assess the influence of different types of vegetation on the stability of transportation infrastructure. Two years of 4-Dimensional ERT monitoring results are presented here, and petrophysical relationships developed in the laboratory are applied to calibrate the resistivity models in order to provide an insight into hydrogeological pathways within a railway cutting. The influence of vegetation type on subsurface moisture pathways and on slope stability is also assessed – here we find that seasonal subsurface changes in moisture content and soil suction are exacerbated by the presence of trees (wooded area). This results in shrink-swell behaviour of the clays comprising the railway cutting, resulting in fissuring and a reduction in shear strength, leading to instability. As such, it is proposed that on slopes comprised of expansive soils, grassed slopes are beneficial for stability. Insights into the use of 4-D ERT for monitoring railway infrastructure gained from this study may be applied to the monitoring of critical geotechnical assets elsewhere.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

The Land Surface in Earth System Reanalysis

Current global atmospheric reanalysis products such as ERA-5, MERRA-2, and JRA-55 provide estimates of land surface states and fluxes, including soil moisture, soil temperature, snow mass, latent and sensible heat fluxes, and runoff. These land surface estimates are based on land surface process models and, depending on the reanalysis product, on precipitation observations or the assimilation of land surface observations of soil moisture and snow conditions from satellites and in situ networks. The land models used in the current global atmospheric reanalysis systems are generally focused on water and energy balance processes, with prescribed static land cover and prescribed climatological vegetation conditions. That is, current global atmospheric reanalysis products do not reflect changes in land use or land cover, nor do they provide dynamic estimates of the land carbon cycle. These and other gaps, including the omission of anthropogenic processes such as irrigation in the land model, need to be addressed in the development of the next generation Earth system reanalysis systems. This presentation provides an overview of the state-of-the-art and future developments of land surface estimates in global reanalysis systems.

Rolf Reichle↗
DOE OSTI2022

A Process‐Model Perspective on Recent Changes in the Carbon Cycle of North America

Continental North America has been found to be a carbon (C) sink over recent decades by multiple studies employing a variety of estimation approaches. However, several key questions and uncertainties remain with these assessments. Here we used results from an ensemble of 19 state-of-the-art dynamic global vegetation models from the TRENDYv9 project to improve these estimates and study the drivers of its interannual variability. Our results show that North America has been a C sink with a magnitude of 0.37 ± 0.38 (mean and one standard deviation) PgC year —1 for the period 2000–2019 (0.31 and 0.44 PgC year —1 in each decade); split into 0.18 ± 0.12 PgC year —1 in Canada (0.15 and 0.20), 0.16 ± 0.17 in the United States (0.14 and 0.17), 0.02 ± 0.05 PgC year —1 in Mexico (0.02 and 0.02) and 0.01 ± 0.02 in Central America and the Caribbean (0.01 and 0.01). About 57% of the new C assimilated by terrestrial ecosystems is allocated into vegetation, 30% into soils, and 13% into litter. Losses of C due to fire account for 41% of the interannual variability of the mean net biome productivity for all North America in the model ensemble. Finally, we show that drought years (e.g., 2002) have the potential to shift the region to a small net C source in the simulations (—0.02 ± 0.46 PgC year —1 ). Our results highlight the importance of identifying the major drivers of the interannual variability of the continental-scale land C cycle along with the spatial distribution of local sink-source dynamics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

A Process-Model Perspective on Recent Changes in the Carbon Cycle of North America

Continental North America has been found to be a carbon (C) sink over recent decades by multiple studies employing a variety of estimation approaches. However, several key questions and uncertainties remain with these assessments. Here we used results from an ensemble of 19 state-of-the-art dynamic global vegetation models from the TRENDYv9 project to improve these estimates and study the drivers of its interannual variability. Our results show that North America has been a C sink with a magnitude of 0.37 ± 0.38 (mean and one standard deviation) PgC year−1 for the period 2000–2019 (0.31 and 0.44 PgC year−1 in each decade); split into 0.18 ± 0.12 PgC year−1 in Canada (0.15 and 0.20), 0.16 ± 0.17 in the United States (0.14 and 0.17), 0.02 ± 0.05 PgC year−1 in Mexico (0.02 and 0.02) and 0.01 ± 0.02 in Central America and the Caribbean (0.01 and 0.01). About 57% of the new C assimilated by terrestrial ecosystems is allocated into vegetation, 30% into soils, and 13% into litter. Losses of C due to fire account for 41% of the interannual variability of the mean net biome productivity for all North America in the model ensemble. Finally, we show that drought years (e.g., 2002) have the potential to shift the region to a small net C source in the simulations (−0.02 ± 0.46 PgC year−1). Our results highlight the importance of identifying the major drivers of the interannual variability of the continental-scale land C cycle along with the spatial distribution of local sink-source dynamics.

Guillermo Murray-Tortarolo↗
DOE OSTI2022

The global decline in the sensitivity of vegetation productivity to precipitation from 2001 to 2018

The sensitivity of vegetation productivity to precipitation (S ppt ) is a key metric for understanding the variations in vegetation productivity under changing precipitation and predicting future changes in ecosystem functions. However, a comprehensive assessment of S ppt over all the global land is lacking. Here, we investigated spatial patterns and temporal changes of S ppt across the global land from 2001 to 2018 with multiple streams of satellite observations. We found consistent spatial patterns of S ppt with different satellite products: S ppt was highest in dry regions while low in humid regions. Grassland and shrubland showed the highest S ppt , and evergreen needle-leaf forest and wetland showed the lowest. Temporally, S ppt showed a generally declining trend over the past two decades (p < .05), yet with clear spatial heterogeneities. The decline in S ppt was especially noticeable in North America and Europe, likely due to the increase in precipitation. In central Russia and Australia, however, S ppt showed an increasing trend. Biome-wise, most ecosystem types exhibited significant decrease in S ppt , while grassland, evergreen broadleaf forest, and mixed forest showed slight increases or non-significant changes in S ppt . Our finding of the overall decline in S ppt implies a potential stabilization mechanism for ecosystem productivity under climate change. However, the revealed S ppt increase for some regions and ecosystem types, in particular global grasslands, suggests that grasslands might be increasingly vulnerable to climatic variability with continuing global climate change.

54 ENVIRONMENTAL SCIENCES↗