Search

Ocean ecosystems play an important role in the global carbon cycle. In this paper, we develop and calibrate a pelagic ocean ecosystem model by simultaneously fitting data from two sites in the North Atlantic: the site of the Bermuda Atlantic Time-series Study (BATS, 32°N, 64°W) and Ocean Weather Ship Station “I” (OWSI, 59°N, 20°W). These sites differ dramatically in the magnitude and timing of primary production, chlorophyll concentrations, and light and nutrient availabilities. We find that a relatively simple model can be simultaneously fit to ecosystem data at both sites. As in the model of Hurtt and Armstrong (1996, Deep-Sea Research II 43(2–3): 653–683), the fit largely depends on the inclusion of multiple size classes of phytoplankton and detritus and a variable chlorophyll-to-biomass (nitrogen) ratio for phyto-plankton. The inclusion of these features in the model enables it to adjust to predict important biological differences temporally within sites and spatially between sites. In addition, the model also includes a phytoplankton self-shading parameter to capture differences in water clarity between sites, and the parameterization of iron limitation at OWSI as a hypothesis for the moderately HNLC condition there. Finally, we suggest that the model performs reasonably well at the Hawaii Ocean Time-series (HOT, 23°N, 158°W), a site where no model calibration was done.

A wide variety of models have illustrated the potential importance of terrestrial biological feedbacks on climate and climate change, yet our ability to make precise predictions is severely limited, due to a high degree of uncertainty. In this paper, after briefly reviewing current models, we present challenges for new terrestrial models and introduce a simple mechanistic approach that may complement existing approaches.

Current trajectories of global change may lead to regime shifts at regional scales, driving coupled human-environment systems to highly degraded states in terms of biodiversity, ecosystem services, and human well-being. For business-as-usual socioeconomic development pathways, regime shifts are projected to occur within the next several decades, to be difficult to reverse, and to have regional- to global-scale impacts on human society. We provide an overview of ecosystem, socioeconomic, and biophysical mechanisms mediating regime shifts and illustrate how these interact at regional scales by aggregation, synergy, and spreading processes. We give detailed examples of interactions for terrestrial ecosystems of central South America and for marine and coastal ecosystems of Southeast Asia. This analysis suggests that degradation of biodiversity and, ecosystem services over the twenty-first century could be far greater than was previously predicted. We identify key policy and management opportunities at regional to global scales to avoid these shifts.

Previous studies have demonstrated the importance of enhanced vegetation growth under future elevated atmospheric CO2 for 21st century climate warming. Surprisingly no study has completed an analogous assessment for the historical period, during which emissions of greenhouse gases increased rapidly and land-use changes (LUC) dramatically altered terrestrial carbon sources and sinks. Using the Geophysical Fluid Dynamics Laboratory comprehensive Earth System Model ESM2G and a reconstruction of the LUC, we estimate that enhanced vegetation growth has lowered the historical atmospheric CO2 concentration by 85 ppm, avoiding an additional 0.31 ± 0.06 °C warming. We demonstrate that without enhanced vegetation growth the total residual terrestrial carbon flux (i.e., the net land flux minus LUC flux) would be a source of 65–82 Gt of carbon (GtC) to atmosphere instead of the historical residual carbon sink of 186–192 GtC, a carbon saving of 251–274 GtC.

The effects of land-use changes on climate are assessed using specified-concentration simulations complementary to the representative concentration pathway 2.6 (RCP2.6) and RCP8.5 scenarios performed for phase 5 of the Coupled Model Intercomparison Project (CMIP5). This analysis focuses on differences in climate and land–atmosphere fluxes between the ensemble averages of simulations with and without land-use changes by the end of the twenty-first century. Even though common land-use scenarios are used, the areas of crops and pastures are specific for each Earth system model (ESM). This is due to different interpretations of land-use classes. The analysis reveals that fossil fuel forcing dominates land-use forcing. In addition, the effects of land-use changes are globally not significant, whereas they are significant for regions with land-use changes exceeding 10%. For these regions, three out of six participating models—the Second Generation Canadian Earth System Model (CanESM2); Hadley Centre Global Environmental Model, version 2 (Earth System) (HadGEM2-ES); and Model for Interdisciplinary Research on Climate, Earth System Model (MIROC-ESM)—reveal statistically significant changes in mean annual surface air temperature. In addition, changes in land surface albedo, available energy, and latent heat fluxes are small but significant for most ESMs in regions affected by land-use changes. These climatic effects are relatively small, as land-use changes in the RCP2.6 and RCP8.5 scenarios are small in magnitude and mainly limited to tropical and subtropical regions. The relative importance of the climatic effects of land-use changes is higher for the RCP2.6 scenario, which considers an expansion of biofuel croplands as a climate mitigation option. The underlying similarity among all models is the loss in global land carbon storage due to land-use changes.

Recent drought events underscore the vulnerability of Amazon forests to understorey fires. The long-term impact of fires on biodiversity and forest carbon stocks depends on the frequency of fire damages and deforestation rates of burned forests. Here, we characterized the spatial and temporal dynamics of understorey fires (1999–2010) and deforestation (2001–2010) in southern Amazonia using new satellite-based estimates of annual fire activity (greater than 50 ha) and deforestation (greater than 10 ha). Understorey forest fires burned more than 85 500 km 2 between 1999 and 2010 (2.8% of all forests). Forests that burned more than once accounted for 16 per cent of all understorey fires. Repeated fire activity was concentrated in Mato Grosso and eastern Pará, whereas single fires were widespread across the arc of deforestation. Routine fire activity in Mato Grosso coincided with annual periods of low night-time relative humidity, suggesting a strong climate control on both single and repeated fires. Understorey fires occurred in regions with active deforestation, yet the interannual variability of fire and deforestation were uncorrelated, and only 2.6 per cent of forests that burned between 1999 and 2008 were deforested for agricultural use by 2010. Evidence from the past decade suggests that future projections of frontier landscapes in Amazonia should separately consider economic drivers to project future deforestation and climate to project fire risk.

Mechanistic models of ocean ecosystem dynamics are of fundamental importance to understanding and predicting the role of marine ecosystems in the oceanic uptake of carbon. In this paper, a new pelagic ecosystem model that is descended from the model of Fasham et al. (Journal of Marine Research, 99 (1990) 591–639) (FDM model) is presented. During model development, the FDM model was first simplified to reduce the number of variables unconstrained by data and to reduce the number of parameters to be estimated. Many alternative simplified model formulations were tested in an attempt to fit 1988–1991 Bermuda Atlantic Time-series Study (BATS) data. The model presented here incorporates the changes found to be important. (i) A feature of the FDM physics that gives rise to a troublesome fall bloom was replaced. (ii) A biodiversity effect was added: the addition of larger algal and detrital size classes as phytoplankton and detrital biomasses increase. (iii) A phytoplankton physiological effect was also added: the adjustment of the chlorophyll-to-nitrogen ratio by phytoplankton in response to light and nutrient availabilities. The new model has only four state variables and a total of 11 biological parameters; yet it fits the average annual cycle in BATS data better than the FDM model. The new model also responds reasonably well to interannual variability in physical forcing. Based on the justification for changes (i)--(iii) from empirical studies and the success of this simple model at fitting BATS data, it is argued that these changes may be generally important. It is also shown that two alternative assumptions about ammonium concentrations lead to very different model calibrations, emphasizing the need for time series data on ammonium.

Plant competition for space is studied using analytical and simulation models. Here, the interaction is viewed as a local competition between juveniles of different species for environmentally variable sites vacated by the random deaths of adults. Because plants are sedentary and have finite fecundity, often only a subset of species will compete for an available site. When a dominant species is recruitment limited, inferior competitors will win some sites by forfeit. It is shown that recruitment limitation allows “winning-by-forfeit” which lessens the effect of competitive asymmetries and slows population and community dynamics. Moreover, since recruitment limitation is likely to be most pronounced in highly diverse communities because of the rarity of many species, it is suggested that there is no conflict between the hypothesis that species-rich plant communities are more influenced by chance and history than regulated by competition, and observations of strong interspecific differences among plants.

This is the full workshop report for the modeling workshop we did for the National Climate Assessment, with DOE support.

Tropical forests are marked by high biological diversity and complex vegetation dynamics that result in a spatially diverse array of forest stand structures (Richards 1952, Denslow 1987, Salati and Vose 1984, Terborgh 1992, Ozanne et al. 2003). Knowledge of the forest structure is vital for estimation of carbon stocks and fluxes (Houghton et al. 2000, 2001), habitat and faunal distributions (Schwarzkopt and Rylands 1989), and interactions between the biosphere and atmosphere (Keller et al. 2004a). With deforestation and land use change occurring throughout the tropics, improved understanding of these dynamic and complex forests are vital for the development of regional and global carbon budgets (Nobre et al. 1991, Werth and Avissar 2002, Houghton et al. 2001, Davidson and Artaxo 2004).

G. C. Stevens showed that an equatorward increase in species richness is often paralleled by a decrease in the mean latitudinal range of species-a pattern he called Rapoport's rule. He reported a similar pattern for elevational gradients and suggested that both latitudinal and elevational species richness gradients may be a result of the Rapoport effect. Using null models that assume no environmental gradients, we show that "nonbiological" gradients in species richness arise inevitably from the assumption of a random latitudinal (or elevational) association between the size and placement of species' ranges. These models predict a peak in species richness at tropical latitudes and, at intermediate elevations (regardless of latitude), common patterns in empirical data. Using simulated sampling from a parametric distribution of ranges that incorporates a richness gradient but no Rapoport effect, we then show that a spurious Rapoport effect can be caused by sampling bias alone, even when total sampling effort...

[1] For more than a century, scientists have recognized the importance of vegetation structure in understanding forest dynamics. Now future satellite missions such as Deformation, Ecosystem Structure, and Dynamics of Ice (DESDynI) hold the potential to provide unprecedented global data on vegetation structure needed to reduce uncertainties in terrestrial carbon dynamics. Here, we briefly review the uses of data on vegetation structure in ecosystem models, develop and analyze theoretical models to quantify model-data requirements, and describe recent progress using a mechanistic modeling approach utilizing a formal scaling method and data on vegetation structure to improve model predictions. Generally, both limited sampling and coarse resolution averaging lead to model initialization error, which in turn is propagated in subsequent model prediction uncertainty and error. In cases with representative sampling, sufficient resolution, and linear dynamics, errors in initialization tend to compensate at larger spatial scales. However, with inadequate sampling, overly coarse resolution data or models, and nonlinear dynamics, errors in initialization lead to prediction error. A robust model-data framework will require both models and data on vegetation structure sufficient to resolve important environmental gradients and tree-level heterogeneity in forest structure globally.

[1] In this paper we present the results of an experiment to measure forest structure and biomass dynamics over the tropical forests of La Selva Biological Station in Costa Rica using a medium resolution lidar. Our main objective was to observe changes in forest canopy height, related height metrics, and biomass, and from these map sources and sinks of carbon across the landscape. The Laser Vegetation Imaging Sensor (LVIS) measured canopy structure over La Selva in 1998 and again in 2005. Changes in waveform metrics were related to field-derived changes in estimated aboveground biomass from a series of old growth and secondary forest plots. Pairwise comparisons of nearly coincident lidar footprints between years showed canopy top height changes that coincided with expected changes based on land cover types. Old growth forests had a net loss in height of −0.33 m, while secondary forests had net gain of 2.08 m. Multiple linear regression was used to relate lidar metrics with biomass changes for combined old growth and secondary forest plots, giving an r2 of 0.65 and an RSE of 10.5 Mg/ha, but both parametric and bootstrapped confidence intervals were wide, suggesting weaker model performance. The plot level relationships were then used to map biomass changes across La Selva using LVIS at a 1 ha scale. The spatial patterns of biomass changes matched expected patterns given the distribution of land cover types at La Selva, with secondary forests showing a gain of 25 Mg/ha and old growth forests showing little change (2 Mg/ha). Prediction intervals were calculated to assess uncertainty for each 1 ha cell to ascertain whether the data and methods used could confidently estimate the sign (source or sink) of the biomass changes. The resulting map showed most of the old growth areas as neutral (no net biomass change), with widely scattered and isolated sources and sinks. Secondary forests in contrast were mostly sinks or neutral, but were never sources. By quantifying both the magnitude of biomass changes and the sensitivity of lidar to detect them, this work will help inform the formulation of future space missions focused on biomass dynamics, such as NASA's Deformation Ecosystem Structure and Dynamics of Ice mission.

Assessing Biodiversity Declines Understanding human impact on biodiversity depends on sound quantitative projection. Pereira et al. (p. 1496 , published online 26 October) review quantitative scenarios that have been developed for four main areas of concern: species extinctions, species abundances and community structure, habitat loss and degradation, and shifts in the distribution of species and biomes. Declines in biodiversity are projected for the whole of the 21st century in all scenarios, but with a wide range of variation. Hoffmann et al. (p. 1503 , published online 26 October) draw on the results of five decades' worth of data collection, managed by the International Union for Conservation of Nature Species Survival Commission. A comprehensive synthesis of the conservation status of the world's vertebrates, based on an analysis of 25,780 species (approximately half of total vertebrate diversity), is presented: Approximately 20% of all vertebrate species are at risk of extinction in the wild, and 11% of threatened birds and 17% of threatened mammals have moved closer to extinction over time. Despite these trends, overall declines would have been significantly worse in the absence of conservation actions.

[1] Abrupt forest disturbances generating gaps >0.001 km 2 impact roughly 0.4–0.7 million km 2 a 1 . Fire, windstorms, logging, and shifting cultivation are dominant disturbances; minor contributors are land conversion, flooding, landslides, and avalanches. All can have substantial impacts on canopy biomass and structure. Quantifying disturbance location, extent, severity, and the fate of disturbed biomass will improve carbon budget estimates and lead to better initialization, parameterization, and/or testing of forest carbon cycle models. Spaceborne remote sensing maps large-scale forest disturbance occurrence, location, and extent, particularly with moderate- and fine-scale resolution passive optical/near-infrared (NIR) instruments. High-resolution remote sensing (e.g., 1 m passive optical/NIR, or small footprint lidar) can map crown geometry and gaps, but has rarely been systematically applied to study small-scale disturbance and natural mortality gap dynamics over large regions. Reducing uncertainty in disturbance and recovery impacts on global forest carbon balance requires quantification of (1) predisturbance forest biomass; (2) disturbance impact on standing biomass and its fate; and (3) rate of biomass accumulation during recovery. Active remote sensing data (e.g., lidar, radar) are more directly indicative of canopy biomass and many structural properties than passive instrument data; a new generation of instruments designed to generate global coverage/sampling of canopy biomass and structure can improve our ability to quantify the carbon balance of Earth’s forests. Generating a high-quality quantitative assessment of disturbance impacts on canopy biomass and structure with spaceborne remote sensing requires comprehensive, well designed, and well coordinated field programs collecting high-quality ground-based data and linkages to dynamical models that can use this information.

[1] We have developed a dynamic land model (LM3V) able to simulate ecosystem dynamics and exchanges of water, energy, and CO2 between land and atmosphere. LM3V is specifically designed to address the consequences of land use and land management changes including cropland and pasture dynamics, shifting cultivation, logging, fire, and resulting patterns of secondary regrowth. Here we analyze the behavior of LM3V, forced with the output from the Geophysical Fluid Dynamics Laboratory (GFDL) atmospheric model AM2, observed precipitation data, and four historic scenarios of land use change for 1700–2000. Our analysis suggests a net terrestrial carbon source due to land use activities from 1.1 to 1.3 GtC/a during the 1990s, where the range is due to the difference in the historic cropland distribution. This magnitude is substantially smaller than previous estimates from other models, largely due to our estimates of a secondary vegetation sink of 0.35 to 0.6 GtC/a in the 1990s and decelerating agricultural land clearing since the 1960s. For the 1990s, our estimates for the pastures' carbon flux vary from a source of 0.37 to a sink of 0.15 GtC/a, and for the croplands our model shows a carbon source of 0.6 to 0.9 GtC/a. Our process-based model suggests a smaller net deforestation source than earlier bookkeeping models because it accounts for decelerated net conversion of primary forest to agriculture and for stronger secondary vegetation regrowth in tropical regions. The overall uncertainty is likely to be higher than the range reported here because of uncertainty in the biomass recovery under changing ambient conditions, including atmospheric CO2 concentration, nutrients availability, and climate.

Assessments from field plots steer much of our current understanding of global change impacts on forest ecosystem structure and function. Recent widespread observations of net carbon accumulation in field plots have suggested that terrestrial ecosystems may be a carbon sink, possibly resulting from climate change and/or CO(2) fertilization. We hypothesize that field plots may inadequately sample inherently rare mortality events, leading to bias when plot level measurements are scaled up to larger domains. In this study, we constructed a simple computer simulation model of forest dynamics to investigate the effects of disturbance patterns on landscape-scale carbon balance estimates. The model was constructed to be a balanced biosphere at the landscape-scale with a uniform spatial pattern of forest growth rates. Disturbance gap-size distributions across the landscape were modelled with a power-law distribution. Small and frequent disturbances result in a well-mixed heterogeneous forest where even small sample plots represented domain-wide behaviour. However, with disturbances dominated by large and rare events, sample plots as large as 50 ha displayed significant bias towards growth. We suggest that the accuracy of domain level estimates of carbon balance from sample plots are highly sensitive to the distribution of disturbance events across the landscape, and to the number, size and distribution of field plots that comprise the estimate. Assumptions that small clusters of field plots may be representative of domain-wide conditions should only be made very cautiously, and warrant further investigation for verification.

dead vegetation, and surface litter), soil, and shade for seasonally matched Landsat 5 images captured before and after the storm. The fractional change in NPV (∆NPV) from 2003 to 2006 provided a quantitative measure of the change in dead vegetation associated with Katrina. A subset for the Pearl River basin was stratified by DNPV to generate disturbance classes, and forest inventory plots were randomly established across the entire DNPV disturbance gradient (Fig. 1A). In each plot, tree mortality and damage, species composition, and biomass loss were quantified. A strong correlation between Landsat-derived DNPVand field-measured tree mortality and damage (fig. S1) enabled development of tree mortality and damagemaps from the Landsat imagery. Next, a second scaling function was generated by comparing Landsat- and MODIS-derived DNPV. With the high temporal frequency and large spatial dimension of MODIS imagery, the Landsat-MODIS scaling provided an assessment of hurricane disturbance across the entire impact region (Fig. 1B). To carry out this scaling, we generated distribution functions for stem density and tree biomass from our forest inventory plots and additional U.S. Forest Service data. A Monte Carlo model was devel

Reducing the rate of biodiversity loss and averting dangerous biodiversity change are international goals, reasserted by the Aichi Targets for 2020 by Parties to the United Nations (UN) Convention on Biological Diversity (CBD) after failure to meet the 2010 target (1, 2). However, there is no global, harmonized observation system for delivering regular, timely data on biodiversity change (3). With the first plenary meeting of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES) soon under way, partners from the Group on Earth Observations Biodiversity Observation Network (GEO BON) (4) are developing—and seeking consensus around—Essential Biodiversity Variables (EBVs) that could form the basis of monitoring programs worldwide.

[1] We used the Regional Atmospheric Modeling System (RAMS) model to investigate the possible impact of land cover change on the July climate of the coterminous United States over the last 290 years. Vegetation data were estimated using the Ecosystem Demography model. The observed change in land cover leads to a weak warming along the Atlantic coast and a strong cooling of more than 1 K over the Midwest and the Great Plains region. The precipitation signal is weaker and shows some reduction in the Midwest because of changes in the patterns of large-scale moisture advection.

Atmospheric and ground-based methods agree on the presence of a carbon sink in the coterminous United States (the United States minus Alaska and Hawaii), and the primary causes for the sink recently have been identified. Projecting the future behavior of the sink is necessary for projecting future net emissions. Here we use two models, the Ecosystem Demography model and a second simpler empirically based model (Miami Land Use History), to estimate the spatio-temporal patterns of ecosystem carbon stocks and fluxes resulting from land-use changes and fire suppression from 1700 to 2100. Our results are compared with other historical reconstructions of ecosystem carbon fluxes and to a detailed carbon budget for the 1980s. Our projections indicate that the ecosystem recovery processes that are primarily responsible for the contemporary U.S. carbon sink will slow over the next century, resulting in a significant reduction of the sink. The projected rate of decrease depends strongly on scenarios of future land use and the long-term effectiveness of fire suppression.

For the period 1980–89, we estimate a carbon sink in the coterminous United States between 0.30 and 0.58 petagrams of carbon per year (petagrams of carbon = 10 15 grams of carbon). The net carbon flux from the atmosphere to the land was higher, 0.37 to 0.71 petagrams of carbon per year, because a net flux of 0.07 to 0.13 petagrams of carbon per year was exported by rivers and commerce and returned to the atmosphere elsewhere. These land-based estimates are larger than those from previous studies (0.08 to 0.35 petagrams of carbon per year) because of the inclusion of additional processes and revised estimates of some component fluxes. Although component estimates are uncertain, about one-half of the total is outside the forest sector. We also estimated the sink using atmospheric models and the atmospheric concentration of carbon dioxide (the tracer-transport inversion method). The range of results from the atmosphere-based inversions contains the land-based estimates. Atmosphere- and land-based estimates are thus consistent, within the large ranges of uncertainty for both methods. Atmosphere-based results for 1980–89 are similar to those for 1985–89 and 1990–94, indicating a relatively stable U.S. sink throughout the period.

Carbon accumulation in forests has been attributed to historical changes in land use and the enhancement of tree growth by CO 2 fertilization, N deposition, and climate change. The relative contribution of land use and growth enhancement is estimated by using inventory data from five states spanning a latitudinal gradient in the eastern United States. Land use is the dominant factor governing the rate of carbon accumulation in these states, with growth enhancement contributing far less than previously reported. The estimated fraction of aboveground net ecosystem production due to growth enhancement is 2.0 ± 4.4%, with the remainder due to land use.

In temperate forests of the eastern US, tropical cyclones are a principal agent of catastrophic wind damage, with dramatic impacts on the structure and functioning of forests. Substantial progress has been made to quantify forest damage and resulting gross carbon emissions from tropical cyclones. However, the net effect of storms on the carbon balance of forests depends not only on the biomass lost in single events, but also on the uptake during recovery from a mosaic of past events. This study estimates the net impacts of tropical cyclones on the carbon balance of US forests over the period 1851–2000. To track both disturbance and recovery and to isolate the effects of storms, a modeling framework is used combining gridded historical estimates of mortality and damage with a mechanistic model using an ensemble approach. The net effect of tropical cyclones on the carbon balance is shown to depend strongly on the spatial and temporal scales of analysis. On average, tropical cyclones contribute a net carbon source over latter half of the 19th century. However, throughout much of the 20th century a regional carbon sink is estimated resulting from periods of forest recovery exceeding damage. The large-scale net annual flux resulting from tropical cyclones varies by up to 50 Tg C yr−1, an amount equivalent to 17%–36% of the US forest carbon sink.

The present and future concentration of atmospheric carbon dioxide depends on both anthropogenic and natural sources and sinks of carbon. Most proposed climate mitigation strategies rely on a progressive transition to carbon-efficient technologies to reduce industrial emissions, substantially supported by policies to maintain or enhance the terrestrial carbon stock in forests and other ecosystems. This strategy may be challenged if terrestrial sequestration capacity is affected by future climate feedbacks, but how and to what extent is little understood. Here, we show that climate mitigation strategies are highly sensitive to future natural disturbance rates (e.g. fires, hurricanes, droughts), because of the potential effect of disturbances on the terrestrial carbon balance. Generally, altered disturbance rates affect the pace of societal and technological transitions required to achieve the mitigation target, with substantial consequences on the energy sector and the global economy. An understanding of the future dynamics and consequences of natural disturbances on terrestrial carbon balance is thus essential for developing robust climate mitigation strategies and policies.

44 pages, 9 tables, 9 figures.-- Corinne Le Quéré ... et al.-- This work is distributed under the Creative Commons Attribution 4.0 License, Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the "global carbon budget" – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on land-cover change data and bookkeeping models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2007–2016), EFF was 9.4 ± 0.5 GtC yr−1, ELUC 1.3 ± 0.7 GtC yr−1, GATM 4.7 ± 0.1 GtC yr−1, SOCEAN 2.4 ± 0.5 GtC yr−1, and SLAND 3.0 ± 0.8 GtC yr−1, with a budget imbalance BIM of 0.6 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For year 2016 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1. Also for 2016, ELUC was 1.3 ± 0.7 GtC yr−1, GATM was 6.1 ± 0.2 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 2.7 ± 1.0 GtC yr−1, with a small BIM of −0.3 GtC. GATM continued to be higher in 2016 compared to the past decade (2007–2016), reflecting in part the high fossil emissions and the small SLAND consistent with El Niño conditions. The global atmospheric CO2 concentration reached 402.8 ± 0.1 ppm averaged over 2016. For 2017, preliminary data for the first 6–9 months indicate a renewed growth in EFF of +2.0 % (range of 0.8 to 3.0 %) based on national emissions projections for China, USA, and India, and projections of gross domestic product (GDP) corrected for recent changes in the carbon intensity of the economy for the rest of the world. This living data update documents changes in the methods and data sets used in this new global carbon budget compared with previous publications of this data set (Le Quéré et al., 2016, 2015b, a, 2014, 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2017 (GCP, 2017)

James S. Clark, Jason Lynch, and Peter Wyckoff are in the Department of Botany, Duke University, Durham, NC 27708; Chris Fastie and Stephen T. Jackson are in the Department of Botany, University of Wyoming, Laramie, WY 82701; George Hurtt and Stephen Pacala are in the Department of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ 08544-1003; Carter Johnson is in the Department of Horticulture and Forestry, South Dakota State University, Brookings, SD 57007; George A. King is at Dynamic Corporation, US EPA National Health and Environmental Effects Research Laboratory, Corvallis, OR 97333; Mark Lewis is in the Math Department, University of Utah, Salt Lake City, UT 84112; Colin Prentice is at the School of Ecology, Lund University, Lund, Sweden; Eugene W. Schupp is in the Department of Rangeland Resources, Utah State University, Logan, UT 84322; and Thompson Webb III is in the Department of Geological Sciences, Brown University, Providence, RI 029121846. ? 1998 American Institute of Biological Sciences. A plausible explanation

Accurately predicting forest dynamics and associated carbon fluxes requires both knowledge of the current state of the ecosystem and an understanding of the underlying processes and environmental conditions that influence the ecosystem processes. Here, we apply a combination of light detection and ranging (lidar) remote sensing (LVIS), an individual-based height-structured ecosystem model (ED), and detailed topographic and climate data to address these requirements to predict carbon dynamics at the Hubbard Brook Experimental Forest (HBEF) in the White Mountains of New Hampshire. Lidar data provided substantial constraints on model estimates of carbon stocks and annual net ecosystem production (ANEP). Lidar-initialized model estimates of carbon stocks (10.77 kg C m−2) were within 5% of the field estimates over the domain and accounted for the 44% decrease in carbon stocks observed between minimum and maximum elevation at HBEF. Lidar-initialized model estimates of ANEP (0.023 kg C m−2 year−1) also compared ...