Great advances have been made in the last decade in quantifying and understanding the spatiotemporal patterns of terrestrial gross primary production (GPP) with ground, atmospheric, and space observations. However, although global GPP estimates exist, each data set relies upon assumptions and none of the available data are based only on measurements. Consequently, there is no consensus on the global total GPP and large uncertainties exist in its benchmarking. The objective of this review is to assess how the different available data sets predict the spatiotemporal patterns of GPP, identify the differences among data sets, and highlight the main advantages/disadvantages of each data set. We compare GPP estimates for the historical period (1990-2009) from two observation-based data sets (Model Tree Ensemble and Moderate Resolution Imaging Spectroradiometer) to coupled carbon-climate models and terrestrial carbon cycle models from the Fifth Climate Model Intercomparison Project and TRENDY projects and to a new hybrid data set (CARBONES). Results show a large range in the mean global GPP estimates. The different data sets broadly agree on GPP seasonal cycle in terms of phasing, while there is still discrepancy on the amplitude. For interannual variability (IAV) and trends, there is a clear separation between the observation-based data that show little IAV and trend, while the process-based models have large GPP variability and significant trends. These results suggest that there is an urgent need to improve observation-based data sets and develop carbon cycle modeling with processes that are currently treated either very simplistically to correctly estimate present GPP and better quantify the future uptake of carbon dioxide by the world's vegetation.

High-latitude ecosystems have the capacity to release large amounts of carbon dioxide (CO2) to the atmosphere in response to increasing temperatures, representing a potentially significant positive feedback within the climate system. Here, we combine aircraft and tower observations of atmospheric CO2 with remote sensing data and meteorological products to derive temporally and spatially resolved year-round CO2 fluxes across Alaska during 2012–2014. We find that tundra ecosystems were a net source of CO2 to the atmosphere annually, with especially high rates of respiration during early winter (October through December). Long-term records at Barrow, AK, suggest that CO2 emission rates from North Slope tundra have increased during the October through December period by 73% ± 11% since 1975, and are correlated with rising summer temperatures. Together, these results imply increasing early winter respiration and net annual emission of CO2 in Alaska, in response to climate warming. Our results provide evidence that the decadal-scale increase in the amplitude of the CO2 seasonal cycle may be linked with increasing biogenic emissions in the Arctic, following the growing season. Early winter respiration was not well simulated by the Earth System Models used to forecast future carbon fluxes in recent climate assessments. Therefore, these assessments may underestimate the carbon release from Arctic soils in response to a warming climate.

[1] The ability to reliably estimate CO 2 fluxes from current in situ atmospheric CO 2 measurements and future satellite CO 2 measurements is dependent on transport model performance at synoptic and shorter timescales. The TransCom continuous experiment was designed to evaluate the performance of forward transport model simulations at hourly, daily, and synoptic timescales, and we focus on the latter two in this paper. Twenty-five transport models or model variants submitted hourly time series of nine predetermined tracers (seven for CO 2 ) at 280 locations. We extracted synoptic-scale variability from daily averaged CO 2 time series using a digital filter and analyzed the results by comparing them to atmospheric measurements at 35 locations. The correlations between modeled and observed synoptic CO 2 variabilities were almost always largest with zero time lag and statistically significant for most models and most locations. Generally, the model results using diurnally varying land fluxes were closer to the observations compared to those obtained using monthly mean or daily average fluxes, and winter was often better simulated than summer. Model results at higher spatial resolution compared better with observations, mostly because these models were able to sample closer to the measurement site location. The amplitude and correlation of model-data variability is strongly model and season dependent. Overall similarity in modeled synoptic CO 2 variability suggests that the first-order transport mechanisms are fairly well parameterized in the models, and no clear distinction was found between the meteorological analyses in capturing the synoptic-scale dynamics.

[1] A forward atmospheric transport modeling experiment has been coordinated by the TransCom group to investigate synoptic and diurnal variations in CO 2 . Model simulations were run for biospheric, fossil, and air-sea exchange of CO 2 and for SF 6 and radon for [2000][2001][2002][2003]. Twenty-five models or model variants participated in the comparison. Hourly concentration time series were submitted for 280 sites along with vertical profiles, fluxes, and meteorological variables at 100 sites. The submitted results have been analyzed for diurnal variations and are compared with observed CO 2 in 2002. Mean summer diurnal cycles vary widely in amplitude across models. The choice of sampling location and model level account for part of the spread suggesting that representation errors in these types of models are potentially large. Despite the model spread, most models simulate the relative variation in diurnal amplitude between sites reasonably well. The modeled diurnal amplitude only shows a weak relationship with vertical resolution across models; differences in near-surface transport simulation appear to play a major role. Examples are also presented where there is evidence that the models show useful skill in simulating seasonal and synoptic changes in diurnal amplitude.

Abstract.Resolving the discrepancies between NEE estimates based upon (1) ground studies and (2) atmospheric inversion results, demands increasingly sophisticated techniques. In this paper we present a high-resolution inversion based upon a regional meteorology model (RAMS) and an underlying biosphere (SiB3) model, both running on an identical 40 km grid over most of North America. Current operational systems like CarbonTracker as well as many previous global inversions including the Transcom suite of inversions have utilized inversion regions formed by collapsing biome-similar grid cells into larger aggregated regions. An extreme example of this might be where corrections to NEE imposed on forested regions on the east coast of the United States might be the same as that imposed on forests on the west coast of the United States while, in reality, there likely exist subtle differences in the two areas, both natural and anthropogenic. Our current inversion framework utilizes a combination of previously employed inversion techniques while allowing carbon flux corrections to be biome independent. Temporally and spatially high-resolution results utilizing biome-independent corrections provide insight into carbon dynamics in North America. In particular, we analyze hourly CO 2 mixing ratio data from a sparse network of eight towers in North America for 2004. A prior estimate of carbon fluxes due to Gross Primary Productivity (GPP) and Ecosystem Respiration (ER) is constructed from the SiB3 biosphere model on a 40 km grid. A combination of transport from Correspondence to: A. E. Schuh (aschuh@atmos.colostate.edu) the RAMS and the Parameterized Chemical Transport Model (PCTM) models is used to forge a connection between upwind biosphere fluxes and downwind observed CO 2 mixing ratio data. A Kalman filter procedure is used to estimate weekly corrections to biosphere fluxes based upon observed CO 2 . RMSE-weighted annual NEE estimates, over an ensemble of potential inversion parameter sets, show a mean estimate 0.57 Pg/yr sink in North America. We perform the inversion with two independently derived boundary inflow conditions and calculate jackknife-based statistics to test the robustness of the model results. We then compare final results to estimates obtained from the CarbonTracker inversion system and at the Southern Great Plains flux site. Results are promising, showing the ability to correct carbon fluxes from the biosphere models over annual and seasonal time scales, as well as over the different GPP and ER components. Additionally, the correlation of an estimated sink of carbon in the South Central United States with regional anomalously high precipitation in an area of managed agricultural and forest lands provides interesting hypotheses for future work.

Abstract. Inter-annual variations in the tropical land carbon (C) balance are a dominant component of the global atmospheric CO2 growth rate. Currently, the lack of quantitative knowledge on processes controlling net tropical ecosystem C balance on inter-annual timescales inhibits accurate understanding and projections of land–atmosphere C exchanges. In particular, uncertainty on the relative contribution of ecosystem C fluxes attributable to concurrent forcing anomalies (concurrent effects) and those attributable to the continuing influence of past phenomena (lagged effects) stifles efforts to explicitly understand the integrated sensitivity of a tropical ecosystem to climatic variability. Here we present a conceptual framework – applicable in principle to any land biosphere model – to explicitly quantify net biospheric exchange (NBE) as the sum of anomaly-induced concurrent changes and climatology-induced lagged changes to terrestrial ecosystem C states (NBE = NBECON+NBELAG). We apply this framework to an observation-constrained analysis of the 2001–2015 tropical C balance: we use a data–model integration approach (CARbon DAta-MOdel fraMework – CARDAMOM) to merge satellite-retrieved land-surface C observations (leaf area, biomass, solar-induced fluorescence), soil C inventory data and satellite-based atmospheric inversion estimates of CO2 and CO fluxes to produce a data-constrained analysis of the 2001–2015 tropical C cycle. We find that the inter-annual variability of both concurrent and lagged effects substantially contributes to the 2001–2015 NBE inter-annual variability throughout 2001–2015 across the tropics (NBECON IAV = 80 % of total NBE IAV, r =  0.76; NBELAG IAV = 64 % of NBE IAV, r = 0.61), and the prominence of NBELAG IAV persists across both wet and dry tropical ecosystems. The magnitude of lagged effect variations on NBE across the tropics is largely attributable to lagged effects on net primary productivity (NPP; NPPLAG IAV 113 % of NBELAG IAV, r = −0.93, p value < 0.05), which emerge due to the dependence of NPP on inter-annual variations in foliar C and plant-available H2O states. We conclude that concurrent and lagged effects need to be explicitly and jointly resolved to retrieve an accurate understanding of the processes regulating the present-day and future trajectory of the terrestrial land C sink.

Accurately quantifying the timing and magnitude of respiration and photosynthesis by high‐latitude ecosystems is important for understanding how a warming climate influences global carbon cycling. Data‐driven estimates of photosynthesis across Arctic regions often rely on satellite‐derived enhanced vegetation index (EVI); we find that satellite observations of solar‐induced chlorophyll fluorescence (SIF) provide a more direct proxy for photosynthesis. We model Alaskan tundra CO2 cycling (2012–2014) according to temperature and shortwave radiation and alternately input EVI or SIF to prescribe the annual seasonal cycle of photosynthesis. We find that EVI‐based seasonality indicates spring “green‐up” to occur 9 days prior to SIF‐based estimates, and that SIF‐based estimates agree with aircraft and tower measurements of CO2. Adopting SIF, instead of EVI, for modeling the seasonal cycle of tundra photosynthesis can result in more accurate estimates of growing season duration and net carbon uptake by arctic vegetation.

Abstract. Synoptic variations of atmospheric CO 2 produced by interactions between weather and surface fluxes are investigated mechanistically and quantitatively in midlatitude and tropical regions using continuous in-situ CO 2 observations in North America, South America and Europe and forward chemical transport model simulations with the Parameterized Chemistry Transport Model. Frontal CO 2 climatologies show consistently strong, characteristic frontal CO 2 signals throughout the midlatitudes of North America and Europe. Transitions between synoptically identifiable CO 2 air masses or transient spikes along the frontal boundary typically characterize these signals. One case study of a summer cold front shows CO 2 gradients organizing with deformational flow along weather fronts, producing strong and spatially coherent variations. In order to differentiate physical and biological controls on synoptic variations in midlatitudes and a site in Amazonia, a boundary layer budget equation is constructed to break down boundary layer CO 2 tendencies into components driven by advection, moist convection, and surface fluxes. This analysis suggests that, in midlatitudes, advection is dominant throughout the year and responsible for 60-70% of day-to-day variations on average, with moist convection contributing less than 5%. At a site in Amazonia, vertical mixing, in particular coupling between convective transport and surface CO 2 flux, is most important, with advection responsible for 26% of variations, moist convection 32%, and surface flux 42%. Transport model sensitivity experiments agree with budget analysis. These results imply the existence of a recharge-discharge mechanism in Amazonia important for controlling synoptic variations of boundary Correspondence to: N. C. Parazoo (nparazoo@atmos.colostate.edu) layer CO 2 , and that forward and inverse simulations should take care to represent moist convective transport. Due to the scarcity of tropical observations at the time of this study, results in Amazonia are not generalized for the tropics, and future work should extend analysis to additional tropical locations.