Climate-limited vegetation change in the conterminous United States of America

The effects of climate change on vegetation composition and distribution are evident in different ecosystems around the world. Although some climate-derived alterations on vegetation are expected to result in changes in lifeform fractional cover, disentangling the direct effects of climate change from different non-climate factors, such as land-use change, is challenging. By applying "Liebig's law of the minimum" in a geospatial context, we determined the climate-limited potential for tree, shrub, herbaceous, and non-vegetation fractional cover change for the conterminous United States and compared these potential rates to observed change rates for the period 1986 to 2018. We found that 10% of the land area of the conterminous United States appears to have climate limitations on the change in fractional cover, with a high proportion of these sites located in arid and semiarid ecosystems in the Southwest part of the country. The rates of change in lifeform fractional cover for the remaining area of the country are likely limited by non-climate factors such as the disturbance regime, land management, land-use history, soil conditions, and species interactions and adaptations.

Global net biome CO2 exchange predicted comparably well using parameter-environment relationships and plant functional types

Accurate estimation and forecasts of net biome CO₂ exchange (NBE) are vital for understanding the role of terrestrial ecosystems in a changing climate. Prior efforts to improve NBE predictions have predominantly focused on increasing models' structural realism (and thus complexity), but parametric error and uncertainty are also key determinants of model skill. Here, we investigate how different parameterization assumptions propagate into NBE prediction errors across the globe, pitting the traditional plant functional type (PFT)-based approach against a novel top-down, machine learning-based "environmental filtering" (EF) approach. To do so, we simulate these contrasting methods for parameter assignment within a flexible model-data fusion framework of the terrestrial carbon cycle (CARDAMOM) at a global scale. In the PFT-based approach, model parameters from a small number of select locations are applied uniformly within regions sharing similar land cover characteristics. In the EF-based approach, a pixel's parameters are predicted based on underlying relationships with climate, soil, and canopy properties. To isolate the role of parametric from structural uncertainty in our analysis, we benchmark the resulting PFT-based and EF-based NBE predictions with estimates from CARDAMOM's Bayesian optimization approach (whereby "true" parameters consistent with a suite of data constraints are retrieved on a pixel-by-pixel basis). When considering the mean absolute error of NBE predictions across time, we find that the EF-based approach matches or outperforms the PFT-based approach at 55% of pixels-a narrow majority. However, NBE estimates from the EF-based approach are susceptible to compensation between errors in component flux predictions and predicted parameters can align poorly with the assumed "true" values. Overall, though, the EF-based approach is comparable to conventional approaches and merits further investigation to better understand and resolve these limitations. This work provides insight into the relationship between terrestrial biosphere model performance and parametric uncertainty, informing efforts to improve model parameterization via PFT-free and trait-based approaches.

Representing plant diversity in land models: An evolutionary approach to make "Functional Types" more functional

Plants are critical mediators of terrestrial mass and energy fluxes, and their structural and functional traits have profound impacts on local and global climate, biogeochemistry, biodiversity, and hydrology. Yet, Earth System Models (ESMs), our most powerful tools for predicting the effects of humans on the coupled biosphere-atmosphere system, simplify the incredible diversity of land plants into a handful of coarse categories of "Plant Functional Types" (PFTs) that often fail to capture ecological dynamics such as biome distributions. The inclusion of more realistic functional diversity is a recognized goal for ESMs, yet there is currently no consistent, widely accepted way to add diversity to models, that is, to determine what new PFTs to add and with what data to constrain their parameters. We review approaches to representing plant diversity in ESMs and draw on recent ecological and evolutionary findings to present an evolution-based functional type approach for further disaggregating functional diversity. Specifically, the prevalence of niche conservatism, or the tendency of closely related taxa to retain similar ecological and functional attributes through evolutionary time, reveals that evolutionary relatedness is a powerful framework for summarizing functional similarities and differences among plant types. We advocate that Plant Functional Types based on dominant evolutionary lineages ("Lineage Functional Types") will provide an ecologically defensible, tractable, and scalable framework for representing plant diversity in next-generation ESMs, with the potential to improve parameterization, process representation, and model benchmarking. We highlight how the importance of evolutionary history for plant function can unify the work of disparate fields to improve predictive modeling of the Earth system.

Development of Land Cover Naturalness in Lithuania on the Edge of the 21st Century: Trends and Driving Factors

Landscape naturalness is an important indicator for supporting sustainable development-driven policies and suggesting associated decisions in land management. This study used CORINE Land Cover data to estimate the changes in land cover naturalness in Lithuania since 1995. All the land cover types were ranked according to naturalness level, ranging from purely anthropogenic to natural landscapes. Spatial patterns of the increase or decline in landscape naturalness were investigated at the level of municipalities. Then, publicly available geographic data were mobilised to explain the reasons behind the trends observed. A minor increase in land cover naturalness in the whole area of Lithuania was observed; however, this increase was statistically insignificant. Nevertheless, statistically significant clusters with both increasing and decreasing levels of land cover naturalness were identified when moving to the level of municipalities. The trends in the development of landscape naturalness were associated with the specificity of agricultural and forestry activities in the municipalities. The suitability of lands for agriculture due to soil, terrain, current land use specifics, and related drivers, such as the availability of land reclamation installations and the intensity of land use, were the main drivers for the declining level of land cover naturalness, usually concentrated in northern and central Lithuania. The land cover naturalness did increase in less suitable areas for agriculture, i.e., in the more forested southeastern municipalities. The study emphasised the need for a systematic and spatially explicit monitoring of the land cover patterns and their changes as well as elaborated proposals for land management policies over the next decade, which were mostly in the line with current European Union and national strategies.

Towards an open and synergistic framework for mapping global land cover

Global land-cover datasets are key sources of information for understanding the complex inter-actions between human activities and global change. They are also among the most critical variables for climate change studies. Over time, the spatial resolution of land cover maps has increased from the kilometer scale to 10-m scale. Single-type historical land cover datasets, including for forests, water, and impervious surfaces, have also been developed in recent years. In this study, we present an open and synergy framework to produce a global land cover dataset that combines supervised land cover classification and aggregation of existing multiple thematic land cover maps with the Google Earth Engine (GEE) cloud computing platform. On the basis of this method of classification and mosaicking, we derived a global land cover dataset for 6 years over a time span of 25 years. The overall accuracies of the six maps were around 75% and the accuracy for change area detection was over 70%. Our product also showed good similarity with the FAO and existing land cover maps.

Priority list of biodiversity metrics to observe from space

Monitoring global biodiversity from space through remotely sensing geospatial patterns has high potential to add to our knowledge acquired by field observation. Although a framework of essential biodiversity variables (EBVs) is emerging for monitoring biodiversity, its poor alignment with remote sensing products hinders interpolation between field observations. This study compiles a comprehensive, prioritized list of remote sensing biodiversity products that can further improve the monitoring of geospatial biodiversity patterns, enhancing the EBV framework and its applicability. The ecosystem structure and ecosystem function EBV classes, which capture the biological effects of disturbance as well as habitat structure, are shown by an expert review process to be the most relevant, feasible, accurate and mature for direct monitoring of biodiversity from satellites. Biodiversity products that require satellite remote sensing of a finer resolution that is still under development are given lower priority (for example, for the EBV class species traits). Some EBVs are not directly measurable by remote sensing from space, specifically the EBV class genetic composition. Linking remote sensing products to EBVs will accelerate product generation, improving reporting on the state of biodiversity from local to global scales.

Global land use for 2015-2100 at 0.05° resolution under diverse socioeconomic and climate scenarios

Global future land use (LU) is an important input for Earth system models for projecting Earth system dynamics and is critical for many modeling studies on future global change. Here we generated a new global gridded LU dataset using the Global Change Analysis Model (GCAM) and a land use spatial downscaling model, named Demeter, under the five Shared Socioeconomic Pathways (SSPs) and four Representative Concentration Pathways (RCPs) scenarios. Compared to existing similar datasets, the presented dataset has a higher spatial resolution (0.05° × 0.05°) and spreads under a more comprehensive set of SSP-RCP scenarios (in total 15 scenarios), and considers uncertainties from the forcing climates. We compared our dataset with the Land Use Harmonization version 2 (LUH2) dataset and found our results are in general spatially consistent with LUH2. The presented dataset will be useful for global Earth system modeling studies, especially for the analysis of the impacts of land use and land cover change and socioeconomics, as well as the characterizing the uncertainties associated with these impacts.

An Assessment of Global Forest Change Datasets for National Forest Monitoring and Reporting

Global Forest Change datasets have the potential to assist countries with national forest measuring, reporting and verification (MRV) requirements. This paper assesses the accuracy of the Global Forest Change data against nationally derived forest change data by comparing the forest loss estimates from the global data with the equivalent data from Guyana for the period 2001–2017. To perform a meaningful comparison between these two datasets, the initial year 2000 forest state needs first to be matched to the definition of forest land cover appropriate to a local national setting. In Guyana, the default definition of 30% tree cover overestimates forest area is by 483,000 ha (18.15%). However, by using a tree canopy cover (i.e., density of tree canopy coverage metric) threshold of 94%, a close match between the Guyana-MRV non-forest area and the Global Forest Change dataset is achieved with a difference of only 24,210 ha (0.91%) between the two maps. A complimentary analysis using a two-stage stratified random sampling design showed the 94% tree canopy cover threshold gave a close correspondence (R2 = 0.98) with the Guyana-MRV data, while the Global Forest Change default setting of 30% tree canopy cover threshold gave a poorer fit (R2 = 0.91). Having aligned the definitions of forest for the Global Forest Change and the Guyana-MRV products for the year 2000, we show that over the period 2001–2017 the Global Forest Change data yielded a 99.34% overall Correspondence with the reference data and a 94.35% Producer’s Accuracy. The Guyana-MRV data yielded a 99.36% overall Correspondence with the reference data and a 95.94% Producer’s Accuracy. A year-by-year analysis of change from 2001–2017 shows that in some years, the Global Forest Change dataset underestimates change, and in other years, such as 2016 and 2017, change is detected that is not forest loss or gain, hence the apparent overestimation. The conclusion is that, when suitably calibrated for percentage tree cover, the Global Forest Change datasets give a good first approximation of forest loss (and, probably, gains). However, in countries with large areas of forest cover and low levels of deforestation, these data should not be relied upon to provide a precise annual loss/gain or rate of change estimate for audit purposes without using independent high-quality reference data.

Fractional Woody Cover Mapping of Texas Savanna at Landsat Scale

Texas savanna experienced substantial woody plant encroachment during the past several decades, resulting in habitat fragmentation and species loss. A detailed map of woody plant abundance and distribution in this area is critically needed for management purpose. This study endeavors to map the fractional woody cover of Texas savanna at Landsat scale (30 m) in an affordable way. The top of atmosphere reflectance, thermal bands, and NDVI layer of Web-Enabled Landsat Data (WELD) of 2012 were used as predictors, together with mean annual precipitation. Classification and Regression Trees (CART) were calibrated against training data of a whole range of fractional woody cover, which were derived from 1-m resolution digital orthophotos of 2012. Validation indicates a reasonable pixel level accuracy of the result fractional woody cover map, with a R-squared value of 0.45. Moreover, the result map clearly depicts the distribution of woody plants across the study area, as reflected by the orthophotos. Furthermore, this new map proves an improvement over the existing Landsat Vegetation Continuous Fields (VCF) tree cover product. The method developed here, combining remote sensing and statistical techniques, can contribute to savanna management through revealing the abundance and distribution of woody plants.

Global land change from 1982 to 2016

Land change is a cause and consequence of global environmental change1,2. Changes in land use and land cover considerably alter the Earth's energy balance and biogeochemical cycles, which contributes to climate change and-in turn-affects land surface properties and the provision of ecosystem services1-4. However, quantification of global land change is lacking. Here we analyse 35 years' worth of satellite data and provide a comprehensive record of global land-change dynamics during the period 1982-2016. We show that-contrary to the prevailing view that forest area has declined globally5-tree cover has increased by 2.24 million km2 (+7.1% relative to the 1982 level). This overall net gain is the result of a net loss in the tropics being outweighed by a net gain in the extratropics. Global bare ground cover has decreased by 1.16 million km2 (-3.1%), most notably in agricultural regions in Asia. Of all land changes, 60% are associated with direct human activities and 40% with indirect drivers such as climate change. Land-use change exhibits regional dominance, including tropical deforestation and agricultural expansion, temperate reforestation or afforestation, cropland intensification and urbanization. Consistently across all climate domains, montane systems have gained tree cover and many arid and semi-arid ecosystems have lost vegetation cover. The mapped land changes and the driver attributions reflect a human-dominated Earth system. The dataset we developed may be used to improve the modelling of land-use changes, biogeochemical cycles and vegetation-climate interactions to advance our understanding of global environmental change1-4,6.

Mapping local and global variability in plant trait distributions

Our ability to understand and predict the response of ecosystems to a changing environment depends on quantifying vegetation functional diversity. However, representing this diversity at the global scale is challenging. Typically, in Earth system models, characterization of plant diversity has been limited to grouping related species into plant functional types (PFTs), with all trait variation in a PFT collapsed into a single mean value that is applied globally. Using the largest global plant trait database and state of the art Bayesian modeling, we created fine-grained global maps of plant trait distributions that can be applied to Earth system models. Focusing on a set of plant traits closely coupled to photosynthesis and foliar respiration-specific leaf area (SLA) and dry mass-based concentrations of leaf nitrogen ([Formula: see text]) and phosphorus ([Formula: see text]), we characterize how traits vary within and among over 50,000 [Formula: see text]-km cells across the entire vegetated land surface. We do this in several ways-without defining the PFT of each grid cell and using 4 or 14 PFTs; each model's predictions are evaluated against out-of-sample data. This endeavor advances prior trait mapping by generating global maps that preserve variability across scales by using modern Bayesian spatial statistical modeling in combination with a database over three times larger than that in previous analyses. Our maps reveal that the most diverse grid cells possess trait variability close to the range of global PFT means.

Representing agriculture in Earth System Models: approaches and priorities for development

Earth System Model (ESM) advances now enable improved representations of spatially and temporally varying anthropogenic climate forcings. One critical forcing is global agriculture, which is now extensive in land-use and intensive in management, owing to 20th century development trends. Agriculture and food systems now contribute nearly 30% of global greenhouse gas emissions and require copious inputs and resources, such as fertilizer, water, and land. Much uncertainty remains in quantifying important agriculture-climate interactions, including surface moisture and energy balances and biogeochemical cycling. Despite these externalities and uncertainties, agriculture is increasingly being leveraged to function as a net sink of anthropogenic carbon, and there is much emphasis on future sustainable intensification. Given its significance as a major environmental and climate forcing, there now exist a variety of approaches to represent agriculture in ESMs. These approaches are reviewed herein, and range from idealized representations of agricultural extent to the development of coupled climate-crop models that capture dynamic feedbacks. We highlight the robust agriculture-climate interactions and responses identified by these modeling efforts, as well as existing uncertainties and model limitations. To this end, coordinated and benchmarking assessments of land-use-climate feedbacks can be leveraged for further improvements in ESM's agricultural representations. We suggest key areas for continued model development, including incorporating irrigation and biogeochemical cycling in particular. Lastly, we pose several critical research questions to guide future work. Our review focuses on ESM representations of climate-surface interactions over managed agricultural lands, rather than on ESMs as an estimation tool for crop yields and productivity.

Remote sensing of plant functional types

Conceptually, plant functional types represent a classification scheme between species and broad vegetation types. Historically, these were based on physiological, structural and/or phenological properties, whereas recently, they have reflected plant responses to resources or environmental conditions. Often, an underlying assumption, based on an economic analogy, is that the functional role of vegetation can be identified by linked sets of morphological and physiological traits constrained by resources, based on the hypothesis of functional convergence. Using these concepts, ecologists have defined a variety of functional traits that are often context dependent, and the diversity of proposed traits demonstrates the lack of agreement on universal categories. Historically, remotely sensed data have been interpreted in ways that parallel these observations, often focused on the categorization of vegetation into discrete types, often dependent on the sampling scale. At the same time, current thinking in both ecology and remote sensing has moved towards viewing vegetation as a continuum rather than as discrete classes. The capabilities of new remote sensing instruments have led us to propose a new concept of optically distinguishable functional types ('optical types') as a unique way to address the scale dependence of this problem. This would ensure more direct relationships between ecological information and remote sensing observations.

The emergence of land change science for global environmental change and sustainability

Land change science has emerged as a fundamental component of global environmental change and sustainability research. This interdisciplinary field seeks to understand the dynamics of land cover and land use as a coupled human-environment system to address theory, concepts, models, and applications relevant to environmental and societal problems, including the intersection of the two. The major components and advances in land change are addressed: observation and monitoring; understanding the coupled system-causes, impacts, and consequences; modeling; and synthesis issues. The six articles of the special feature are introduced and situated within these components of study.

Satellite detection of bird communities in tropical countryside

The future of biodiversity hinges partly on realizing the potentially high conservation value of human-dominated countryside. The characteristics of the countryside that promote biodiversity preservation remain poorly understood, however, particularly at the fine scales at which individual farmers tend to make land use decisions. To address this problem, we explored the use of a rapid remote sensing method for estimating bird community composition in tropical countryside, using a two-step process. First, we asked how fine-grained variation in land cover affected community composition. Second, we determined whether the observed changes in community composition correlated with three easily accessible remote sensing metrics (wetness, greenness, and brightness), derived from performing a tasseled-cap transformation on a Landsat Enhanced Thematic Mapper Plus image. As a comparison, we also examined whether the most commonly used remote sensing indicator in ecology, the Normalized Difference Vegetation Index (NDVI), correlated with community composition. We worked within an agricultural landscape in southern Costa Rica, where the land comprised a complex and highly heterogeneous mosaic of remnant native vegetation, pasture, coffee cultivation, and other crops. In this region, we selected 12 study sites (each < 60 ha) that encompassed the range of available land cover possibilities in the countryside. Within each site, we surveyed bird communities within all major land cover types, and we conducted detailed field mapping of land cover. We found that the number of forest-affiliated species increased with forest cover and decreased with residential area across sites. Conversely, the number of agriculture-affiliated species using forest increased with land area devoted to agricultural and residential uses. Interestingly, we found that the wetness and brightness metrics predicted the number of forest- and agriculture-affiliated species within a site as well as did detailed field-generated maps of land cover. In contrast, NDVI and the closely correlated greenness metric did not correlate with land cover or with bird communities. Our study shows the strong potential of the tasseled-cap transformation as a tool for assessing the conservation value of countryside for biodiversity.

The relevance of phylogeny to studies of global change

Phylogenetic thinking has infiltrated many areas of biological research, but has had little impact on studies of global ecology or climate change. Here, we illustrate how phylogenetic information can be relevant to understanding vegetation-atmosphere dynamics at ecosystem or global scales by re-analyzing a data set of carbonic anhydrase (CA) activity in leaves that was used to estimate terrestrial gross primary productivity. The original calculations relied on what appeared to be low CA activity exclusively in C4 grasses, but our analyses indicate that such activity might instead characterize the PACCAD grass lineage, which includes many widespread C3 species. We outline how phylogenetics can guide better taxon sampling of key physiological traits, and discuss how the emerging field of phyloinformatics presents a promising new framework for scaling from organism physiology to global processes.

Carbon emissions from tropical deforestation and regrowth based on satellite observations for the 1980s and 1990s

Carbon fluxes from tropical deforestation and regrowth are highly uncertain components of the contemporary carbon budget, due in part to the lack of spatially explicit and consistent information on changes in forest area. We estimate fluxes for the 1980s and 1990s using subpixel estimates of percent tree cover derived from coarse (National Oceanic and Atmospheric Administration's Advanced Very High Resolution Radiometer) satellite data in combination with a terrestrial carbon model. The satellite-derived estimates of change in forest area are lower than national reports and remote-sensing surveys from the United Nations Food and Agriculture Organization Forest Resource Assessment (FRA) in all tropical regions, especially for the 1980s. However, our results indicate that the net rate of tropical forest clearing increased approximately 10% from the 1980s to 1990s, most notably in southeast Asia, in contrast to an 11% reduction reported by the FRA. We estimate net mean annual carbon fluxes from tropical deforestation and regrowth to average 0.6 (0.3-0.8) and 0.9 (0.5-1.4) petagrams (Pg).yr(-1) for the 1980s and 1990s, respectively. Compared with previous estimates of 1.9 (0.6-2.5) Pg.yr(-1) based on FRA national statistics of changes in forest area, this alternative estimate suggests less "missing" carbon from the global carbon budget but increasing emissions from tropical land-use change.

Recent patterns and mechanisms of carbon exchange by terrestrial ecosystems

Knowledge of carbon exchange between the atmosphere, land and the oceans is important, given that the terrestrial and marine environments are currently absorbing about half of the carbon dioxide that is emitted by fossil-fuel combustion. This carbon uptake is therefore limiting the extent of atmospheric and climatic change, but its long-term nature remains uncertain. Here we provide an overview of the current state of knowledge of global and regional patterns of carbon exchange by terrestrial ecosystems. Atmospheric carbon dioxide and oxygen data confirm that the terrestrial biosphere was largely neutral with respect to net carbon exchange during the 1980s, but became a net carbon sink in the 1990s. This recent sink can be largely attributed to northern extratropical areas, and is roughly split between North America and Eurasia. Tropical land areas, however, were approximately in balance with respect to carbon exchange, implying a carbon sink that offset emissions due to tropical deforestation. The evolution of the terrestrial carbon sink is largely the result of changes in land use over time, such as regrowth on abandoned agricultural land and fire prevention, in addition to responses to environmental changes, such as longer growing seasons, and fertilization by carbon dioxide and nitrogen. Nevertheless, there remain considerable uncertainties as to the magnitude of the sink in different regions and the contribution of different processes.

From tropics to tundra: global convergence in plant functioning

Despite striking differences in climate, soils, and evolutionary history among diverse biomes ranging from tropical and temperate forests to alpine tundra and desert, we found similar interspecific relationships among leaf structure and function and plant growth in all biomes. Our results thus demonstrate convergent evolution and global generality in plant functioning, despite the enormous diversity of plant species and biomes. For 280 plant species from two global data sets, we found that potential carbon gain (photosynthesis) and carbon loss (respiration) increase in similar proportion with decreasing leaf life-span, increasing leaf nitrogen concentration, and increasing leaf surface area-to-mass ratio. Productivity of individual plants and of leaves in vegetation canopies also changes in constant proportion to leaf life-span and surface area-to-mass ratio. These global plant functional relationships have significant implications for global scale modeling of vegetation-atmosphere CO2 exchange.