Search

Abstract Grasslands grow in a sequence of seasonal growth stages that respond to both climate and weather, and these relationships can be used to establish a strategy for predicting plant phenology. Current plant states (phenophase) can be represented as one of established growth stages that dictate carbon allocation and leaf photosynthetic capacity. Calculating daily phenophases from climate and environmental relationships allows for sequential growth stages (i.e., well‐defined seasonal cycles with a single growth period) or dynamic growth stages (i.e., multiple growth periods during a growing season). Senescence results from biomass mortality in response to environmental conditions. This approach uses a single mechanistic framework to represent grassland ecology, removing the dependence on satellite‐based vegetation indices and individual site tuning of parameters. Rather than being specified, a variety of properties emerge, from allometric relationships such as root‐shoot ratios, to behavior across moisture gradients, to interannual variability in growing season lengths, carbon stores, and land surface fluxes. Using dynamic phenology stages to link biophysical and biogeochemical processes provides a mechanism to predict self‐consistent land‐atmosphere exchanges of carbon, water, energy, radiation, and momentum, as well as carbon storage in cascading pools of biomass; and describing these processes in a mathematically determinate model makes them clear, testable, and usable for predictions. This paper describes this new phenology method as it is implemented in the Simple Biosphere Model Version 4 (SiB4), and a companion paper evaluates this method at grassland sites worldwide.

Abstract The interaction of land with the atmosphere is sensitive to soil moisture (W). Evapotranspiration (ET) reacts to soil moisture in a nonlinear way, f(W), as soils dry from saturation to wilt point. This nonlinear behavior and the fact that soil moisture varies on scales as small as 1–10 m in nature, while numerical general circulation models (GCMs) have grid cell sizes on the order of 1 to 100s of kilometers, makes the calculation of grid cell‐average ET problematic. It is impractical to simulate the land in GCMs on the small scales seen in nature, so techniques have been developed to represent subgrid scale heterogeneity, including: (1) statistical‐dynamical representations of grid subelements of varying wetness, (2) relaxation of f(W), (3) moderating f(W) with approximations of catchment hydrology, (4) “tiling” the landscape into vegetation types, and (5) hyperresolution. Here we present an alternative method for representing subgrid variability in W, one proven in a conceptual framework where landscape‐scale W is represented as a series of “Bins” of increasing wetness from dry to saturated. The grid cell‐level f(W) is defined by the integral of the fractional area of the wetness bins and the value of f(W) associated with each. This approach accounts for the spatiotemporal dynamics of W. We implemented this approach in the SiB3 land surface parameterization and then evaluated its performance against a control, which assumes a horizontally uniform field of W. We demonstrate that the Bins method, with a physical basis, attenuates unrealistic jumps in model state and ET seen in the control runs.

Earth system processes exhibit complex patterns across time, as do the models that seek to replicate these processes. Model output may or may not be significantly related to observations at different times and on different frequencies. Conventional model diagnostics provide an aggregate view of model–data agreement, but usually do not identify the time and frequency patterns of model–data disagreement, leaving unclear the steps required to improve model response to environmental drivers that vary on characteristic frequencies. Wavelet coherence can quantify the times and timescales at which two time series, for example time series of models and measurements, are significantly different. We applied wavelet coherence to interpret the predictions of 20 ecosystem models from the North American Carbon Program (NACP) Site-Level Interim Synthesis when confronted with eddy-covariance-measured net ecosystem exchange (NEE) from 10 ecosystems with multiple years of available data. Models were grouped into classes with similar approaches for incorporating phenology, the calculation of NEE, the inclusion of foliar nitrogen (N), and the use of model–data fusion. Models with prescribed, rather than prognostic, phenology often fit NEE observations better on annual to interannual timescales in grassland, wetland and agricultural ecosystems. Models that calculated NEE as net primary productivity (NPP) minus heterotrophic respiration (HR) rather than gross ecosystem productivity (GPP) minus ecosystem respiration (ER) fit better on annual timescales in grassland and wetland ecosystems, but models that calculated NEE as GPP minus ER were superior on monthly to seasonal timescales in two coniferous forests. Models that incorporated foliar nitrogen (N) data were successful at capturing NEE variability on interannual (multiple year) timescales at Howland Forest, Maine. The model that employed a model–data fusion approach often, but not always, resulted in improved fit to data, suggesting that improving model parameterization is important but not the only step for improving model performance. Combined with previous findings, our results suggest that the mechanisms driving daily and annual NEE variability tend to be correctly simulated, but the magnitude of these fluxes is often erroneous, suggesting that model parameterization must be improved. Few NACP models correctly predicted fluxes on seasonal and interannual timescales where spectral energy in NEE observations tends to be low, but where phenological events, multi-year oscillations in climatological drivers, and ecosystem succession are known to be important for determining ecosystem function. Mechanistic improvements to models must be made to replicate observed NEE variability on seasonal and interannual timescales.

Conventionally, measurements of carbon isotopes in atmospheric CO2 (δ13CO2) have been used to partition fluxes between terrestrial and ocean carbon pools. However, novel analytical approaches combined with an increase in the spatial extent and frequency of δ13CO2 measurements allow us to conduct a global analysis of δ13CO2 variability to infer the isotopic composition of source CO2 to the atmosphere (δs). This global analysis yields coherent seasonal patterns of isotopic enrichment. Our results indicate that seasonal values of δs are more highly correlated with vapor pressure deficit (r = 0.404) than relative humidity (r = 0.149). We then evaluate two widely used stomatal conductance models and determine that the Leuning Model, which is primarily driven by vapor pressure deficit is more effective globally at predicting δs (RMSE = 1.6‰) than the Ball-Woodrow-Berry model, which is driven by relative humidity (RMSE = 2.7‰). Thus stomatal conductance on a global scale may be more sensitive to changes in vapor pressure deficit than relative humidity. This approach highlights a new application of using δ13CO2 measurements to validate global models.

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 CO2 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 the RAMS and the Parameterized Chemical Transport Model (PCTM) models is used to forge a connection between upwind biosphere fluxes and downwind observed CO2 mixing ratio data. A Kalman filter procedure is used to estimate weekly corrections to biosphere fluxes based upon observed CO2. 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.

Synoptic variations of atmospheric CO2 produced by interactions between weather and surface fluxes are investigated mechanistically and quantitatively in midlatitude and tropical regions using continuous in-situ CO2 observations in North America, South America and Europe and forward chemical transport model simulations with the Parameterized Chemistry Transport Model. Frontal CO2 climatologies show consistently strong, characteristic frontal CO2 signals throughout the midlatitudes of North America and Europe. Transitions between synoptically identifiable CO2 air masses or transient spikes along the frontal boundary typically characterize these signals. One case study of a summer cold front shows CO2 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 CO2 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 CO2 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 layer CO2, 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.

No abstract available.