Search

Permafrost-affected ecosystems of the Arctic–boreal zone in northwestern North America are undergoing profound transformation due to rapid climate change. NASA's Arctic Boreal Vulnerability Experiment (ABoVE) is investigating characteristics that make these ecosystems vulnerable or resilient to this change. ABoVE employs airborne synthetic aperture radar (SAR) as a powerful tool to characterize tundra, taiga, peatlands, and fens. Here, we present an annotated guide to the L-band and P-band airborne SAR data acquired during the 2017, 2018, 2019, and 2022 ABoVE airborne campaigns. We summarize the ∼80 SAR flight lines and how they fit into the ABoVE experimental design (Miller et al., 2023; https://doi.org/10.3334/ORNLDAAC/2150). The Supplement provides hyperlinks to extensive maps, tables, and every flight plan as well as individual flight lines. We illustrate the interdisciplinary nature of airborne SAR data with examples of preliminary results from ABoVE studies including boreal forest canopy structure from TomoSAR data over Delta Junction, AK, and the Boreal Ecosystem Research and Monitoring Sites (BERMS) area in northern Saskatchewan and active layer thickness and soil moisture data product validation. This paper is presented as a guide to enable interested readers to fully explore the ABoVE L- and P-band airborne SAR data (https://uavsar.jpl.nasa.gov/cgi-bin/data.pl).

Climate change, rising CO2 concentration, and land use and land cover change (LULCC) are primary driving forces for terrestrial gross primary productivity (GPP), but their impacts on the temporal changes in GPP are uncertain. In this study, the effects of the three main factors on the interannual variation (IAV) and seasonal cycle amplitude (SCA) of GPP in China were investigated using 12 terrestrial biosphere models from the Multi-scale Synthesis and Terrestrial Model Intercomparison Project. The simulated ensemble mean value of China's GPP between 1981 and 2010, driven by common climate forcing, LULCC and CO2 data, was found to be 7.4±1.8 Pg C yr−1. In general, climate was the dominant control factor of the annual trends, IAV and seasonality of China's GPP. The overall rising CO2 led to enhanced plant photosynthesis, thus increasing annual mean and IAV of China's total GPP, especially in northeastern and southern China, where vegetation is dense. LULCC decreased the IAV of China's total GPP by ∼7 %, whereas rising CO2 induced an increase of 8 %. Compared to climate change and elevated CO2, LULCC showed less contributions to GPP's temporal variation, and its impact acted locally, mainly in southwestern China. Furthermore, this study also examined subregional contributions to the temporal changes in China's total GPP. Southern and southeastern China showed higher contributions to China's annual GPP, whereas southwestern and central parts of China explained larger fractions of the IAV in China's GPP.

This study investigates and compares soil moisture and hydrology projections of broadly used land models with permafrost processes and highlights the causes and impacts of permafrost zone soil moisture projections. Climate models project warmer temperatures and increases in precipitation (P) which will intensify evapotranspiration (ET) and runoff in land models. However, this study shows that most models project a long-term drying of the surface soil (0–20 cm) for the permafrost region despite increases in the net air–surface water flux (P-ET). Drying is generally explained by infiltration of moisture to deeper soil layers as the active layer deepens or permafrost thaws completely. Although most models agree on drying, the projections vary strongly in magnitude and spatial pattern. Land models tend to agree with decadal runoff trends but underestimate runoff volume when compared to gauge data across the major Arctic river basins, potentially indicating model structural limitations. Coordinated efforts to address the ongoing challenges presented in this study will help reduce uncertainty in our capability to predict the future Arctic hydrological state and associated land–atmosphere biogeochemical processes across spatial and temporal scales.

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.

Recent observations of near-surface soil temperatures over the circumpolar Arctic show accelerated warming of permafrost-affected soils. The availability of a comprehensive near-surface permafrost and active layer dataset is critical to better understanding climate impacts and to constraining permafrost thermal conditions and its spatial distribution in land system models. We compiled a soil temperature dataset from 72 monitoring stations in Alaska using data collected by the U.S. Geological Survey, the National Park Service, and the University of Alaska Fairbanks permafrost monitoring networks. The array of monitoring stations spans a large range of latitudes from 60.9 to 71.3° N and elevations from near sea level to ∼ 1300 m, comprising tundra and boreal forest regions. This dataset consists of monthly ground temperatures at depths up to 1 m, volumetric soil water content, snow depth, and air temperature during 1997–2016. These data have been quality controlled in collection and processing. Meanwhile, we implemented data harmonization evaluation for the processed dataset. The final product (PF-AK, v0.1) is available at the Arctic Data Center (https://doi.org/10.18739/A2KG55).

To improve our understanding of the global carbon balance and its representation in terrestrial biosphere models, we present here a first dual-species application of the CarbonTracker Data Assimilation System (CTDAS). The system's modular design allows for assimilating multiple atmospheric trace gases simultaneously to infer exchange fluxes at the Earth surface. In the prototype discussed here, we interpret signals recorded in observed carbon dioxide (CO2) along with observed ratios of its stable isotopologues 13CO2∕12CO2 (δ13C). The latter is in particular a valuable tracer to untangle CO2 exchange from land and oceans. Potentially, it can also be used as a proxy for continent-wide drought stress in plants, largely because the ratio of 13CO2 and 12CO2 molecules removed from the atmosphere by plants is dependent on moisture conditions.The dual-species CTDAS system varies the net exchange fluxes of both 13CO2 and CO2 in ocean and terrestrial biosphere models to create an ensemble of 13CO2 and CO2 fluxes that propagates through an atmospheric transport model. Based on differences between observed and simulated 13CO2 and CO2 mole fractions (and thus δ13C) our Bayesian minimization approach solves for weekly adjustments to both net fluxes and isotopic terrestrial discrimination that minimizes the difference between observed and estimated mole fractions.With this system, we are able to estimate changes in terrestrial δ13C exchange on seasonal and continental scales in the Northern Hemisphere where the observational network is most dense. Our results indicate a decrease in stomatal conductance on a continent-wide scale during a severe drought. These changes could only be detected after applying combined atmospheric CO2 and δ13C constraints as done in this work. The additional constraints on surface CO2 exchange from δ13C observations neither affected the estimated carbon fluxes nor compromised our ability to match observed CO2 variations. The prototype presented here can be of great benefit not only to study the global carbon balance but also to potentially function as a data-driven diagnostic to assess multiple leaf-level exchange parameterizations in carbon-climate models that influence the CO2, water, isotope, and energy balance.

Many ecosystem processes that influence Earth system feedbacks – vegetation growth, water and nutrient cycling, disturbance regimes – are strongly influenced by multidecadal- to millennial-scale climate variations that cannot be directly observed. Paleoclimate records provide information about these variations, forming the basis of our understanding and modeling of them. Fossil pollen records are abundant in the NE US, but cannot simultaneously provide information about paleoclimate and past vegetation in a modeling context because this leads to circular logic. If pollen data are used to constrain past vegetation changes, then the remaining paleoclimate archives in the northeastern US (NE US) are quite limited. Nonetheless, a growing number of diverse reconstructions have been developed but have not yet been examined together. Here we conduct a systematic review, assessment, and comparison of paleotemperature and paleohydrological proxies from the NE US for the last 3000 years. Regional temperature reconstructions (primarily summer) show a long-term cooling trend (1000 BCE–1700 CE) consistent with hemispheric-scale reconstructions, while hydroclimate data show gradually wetter conditions through the present day. Multiple proxies suggest that a prolonged, widespread drought occurred between 550 and 750 CE. Dry conditions are also evident during the Medieval Climate Anomaly, which was warmer and drier than the Little Ice Age and drier than today. There is some evidence for an acceleration of the longer-term wetting trend in the NE US during the past century; coupled with an abrupt shift from decreasing to increasing temperatures in the past century, these changes could have wide-ranging implications for species distributions, ecosystem dynamics, and extreme weather events. More work is needed to gather paleoclimate data in the NE US to make inter-proxy comparisons and to improve estimates of uncertainty in reconstructions.

Abstract Terrestrial ecosystems play a vital role in regulating the accumulation of carbon (C) in the atmosphere. Understanding the factors controlling land C uptake is critical for reducing uncertainties in projections of future climate. The relative importance of changing climate, rising atmospheric CO2, and other factors, however, remains unclear despite decades of research. Here, we use an ensemble of land models to show that models disagree on the primary driver of cumulative C uptake for 85% of vegetated land area. Disagreement is largest in model sensitivity to rising atmospheric CO2 which shows almost twice the variability in cumulative land uptake since 1901 (1 s.d. of 212.8 PgC vs. 138.5 PgC, respectively). We find that variability in CO2 and temperature sensitivity is attributable, in part, to their compensatory effects on C uptake, whereby comparable estimates of C uptake can arise by invoking different sensitivities to key environmental conditions. Conversely, divergent estimates of C uptake can occur despite being based on the same environmental sensitivities. Together, these findings imply an important limitation to the predictability of C cycling and climate under unprecedented environmental conditions. We suggest that the carbon modeling community prioritize a probabilistic multi-model approach to generate more robust C cycle projections.

Respiration in frozen soils is limited to thawed substrate within the thin water films surrounding soil particles. As temperatures decrease and the films become thinner, the available substrate also decreases, with respiration effectively ceasing at −8 °C. Traditional exponential scaling factors to model this effect do not account for substrate availability and do not work at the century to millennial timescales required to model the fate of the nearly 1100 Gt of carbon in permafrost regions. The exponential scaling factor produces a false, continuous loss of simulated permafrost carbon in the 20th century and biases in estimates of potential emissions as permafrost thaws in the future. Here we describe a new frozen biogeochemistry parameterization that separates the simulated carbon into frozen and thawed pools to represent the effects of substrate availability. We parameterized the liquid water fraction as a function of temperature based on observations and use this to transfer carbon between frozen pools and thawed carbon in the thin water films. The simulated volumetric water content (VWC) as a function of temperature is consistent with observed values and the simulated respiration fluxes as a function of temperature are consistent with results from incubation experiments. The amount of organic matter was the single largest influence on simulated VWC and respiration fluxes. Future versions of the parameterization should account for additional, non-linear effects of substrate diffusion in thin water films on simulated respiration. Controlling respiration in frozen soils based on substrate availability allows us to maintain a realistic permafrost carbon pool by eliminating the continuous loss caused by the original exponential scaling factors. The frozen biogeochemistry parameterization is a useful way to represent the effects of substrate availability on soil respiration in model applications that focus on century to millennial timescales in permafrost regions.

Permafrost-affected soils contain twice as much carbon as currently exists in the atmosphere. Studies show that warming of the perennially frozen ground could initiate significant release of the frozen soil carbon into the atmosphere. Initializing the frozen permafrost carbon with the observed soil carbon distribution from the Northern Circumpolar Soil Carbon Database reduces the uncertainty associated with the modeling of the permafrost carbon feedback. To improve permafrost thermal and carbon dynamics we implemented a dynamic surface organic layer with vertical carbon redistribution, and introduced dynamic root growth controlled by active layer thickness, which improved soil carbon exchange between frozen and thawed pools. These changes increased the initial amount of simulated frozen carbon from 313 to 560 Gt C, consistent with observed frozen carbon stocks, and increased the spatial correlation of the simulated and observed distribution of frozen carbon from 0.12 to 0.63.

We present an enhanced version of the SiBCASA terrestrial biosphere model that is extended with (a) biomass burning emissions from the SiBCASA carbon pools using remotely sensed burned area from the Global Fire Emissions Database (GFED), (b) an isotopic discrimination scheme that calculates 13C signatures of photosynthesis and autotrophic respiration, and (c) a separate set of 13C pools to carry isotope ratios into heterotrophic respiration. We quantify in this study the terrestrial exchange of CO2 and 13CO2 as a function of environmental changes in humidity and biomass burning. The implementation of biomass burning yields similar fluxes as CASA-GFED both in magnitude and spatial patterns. The implementation of isotope exchange gives a global mean discrimination value of 15.2‰, ranges between 4 and 20‰ depending on the photosynthetic pathway in the plant, and compares favorably (annually and seasonally) with other published values. Similarly, the isotopic disequilibrium is similar to other studies that include a small effect of biomass burning as it shortens the turnover of carbon. In comparison to measurements, a newly modified starch/sugar storage pool propagates the isotopic discrimination anomalies to respiration much better. In addition, the amplitude of the drought response by SiBCASA is lower than suggested by the measured isotope ratios. We show that a slight increase in the stomatal closure for large vapor pressure deficit would amplify the respired isotope ratio variability. Our study highlights the importance of isotope ratio observations of 13C to assess and improve biochemical models like SiBCASA, especially with regard to the allocation and turnover of carbon and the responses to drought.

Ecosystems are important and dynamic components of the global carbon cycle, and terrestrial biospheric models (TBMs) are crucial tools in further understanding of how terrestrial carbon is stored and exchanged with the atmosphere across a variety of spatial and temporal scales. Improving TBM skills, and quantifying and reducing their estimation uncertainties, pose significant challenges. The Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal multi-scale and multi-model intercomparison effort set up to tackle these challenges. The MsTMIP protocol prescribes standardized environmental driver data that are shared among model teams to facilitate model–model and model–observation comparisons. This paper describes the global and North American environmental driver data sets prepared for the MsTMIP activity to both support their use in MsTMIP and make these data, along with the processes used in selecting/processing these data, accessible to a broader audience. Based on project needs and lessons learned from past model intercomparison activities, we compiled climate, atmospheric CO2 concentrations, nitrogen deposition, land use and land cover change (LULCC), C3 / C4 grasses fractions, major crops, phenology and soil data into a standard format for global (0.5° × 0.5° resolution) and regional (North American: 0.25° × 0.25° resolution) simulations. In order to meet the needs of MsTMIP, improvements were made to several of the original environmental data sets, by improving the quality, and/or changing their spatial and temporal coverage, and resolution. The resulting standardized model driver data sets are being used by over 20 different models participating in MsTMIP. The data are archived at the Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC, http://daac.ornl.gov) to provide long-term data management and distribution.

Climate change is leading to a disproportionately large warming in the high northern latitudes, but the magnitude and sign of the future carbon balance of the Arctic are highly uncertain. Using 40 terrestrial biosphere models for the Alaskan Arctic from four recent model intercomparison projects – NACP (North American Carbon Program) site and regional syntheses, TRENDY (Trends in net land atmosphere carbon exchanges), and WETCHIMP (Wetland and Wetland CH4 Inter-comparison of Models Project) – we provide a baseline of terrestrial carbon cycle uncertainty, defined as the multi-model standard deviation (σ) for each quantity that follows. Mean annual absolute uncertainty was largest for soil carbon (14.0 ± 9.2 kg C m−2), then gross primary production (GPP) (0.22 ± 0.50 kg C m−2 yr−1), ecosystem respiration (Re) (0.23 ± 0.38 kg C m−2 yr−1), net primary production (NPP) (0.14 ± 0.33 kg C m−2 yr−1), autotrophic respiration (Ra) (0.09 ± 0.20 kg C m−2 yr−1), heterotrophic respiration (Rh) (0.14 ± 0.20 kg C m−2 yr−1), net ecosystem exchange (NEE) (−0.01 ± 0.19 kg C m−2 yr−1), and CH4 flux (2.52 ± 4.02 g CH4 m−2 yr−1). There were no consistent spatial patterns in the larger Alaskan Arctic and boreal regional carbon stocks and fluxes, with some models showing NEE for Alaska as a strong carbon sink, others as a strong carbon source, while still others as carbon neutral. Finally, AmeriFlux data are used at two sites in the Alaskan Arctic to evaluate the regional patterns; observed seasonal NEE was captured within multi-model uncertainty. This assessment of carbon cycle uncertainties may be used as a baseline for the improvement of experimental and modeling activities, as well as a reference for future trajectories in carbon cycling with climate change in the Alaskan Arctic and larger boreal region.

Drained thermokarst lake basins (DTLBs) are ubiquitous landforms on Arctic tundra lowland. Their dynamic states are seldom investigated, despite their importance for landscape stability, hydrology, nutrient fluxes, and carbon cycling. Here we report results based on high-resolution Interferometric Synthetic Aperture Radar (InSAR) measurements using space-borne data for a study area located on the North Slope of Alaska near Prudhoe Bay, where we focus on the seasonal thaw settlement within DTLBs, averaged between 2006 and 2010. The majority (14) of the 18 DTLBs in the study area exhibited seasonal thaw settlement of 3–4 cm. However, four of the DTLBs examined exceeded 4 cm of thaw settlement, with one basin experiencing up to 12 cm. Combining the InSAR observations with the in situ active layer thickness measured using ground penetrating radar and mechanical probing, we calculated thaw strain, an index of thaw settlement strength along a transect across the basin that underwent large thaw settlement. We found thaw strains of 10–35% at the basin center, suggesting the seasonal melting of ground ice as a possible mechanism for the large settlement. These findings emphasize the dynamic nature of permafrost landforms, demonstrate the capability of the InSAR technique to remotely monitor surface deformation of individual DTLBs, and illustrate the combination of ground-based and remote sensing observations to estimate thaw strain. Our study highlights the need for better description of the spatial heterogeneity of landscape-scale processes for regional assessment of surface dynamics on Arctic coastal lowlands.

Terrestrial biosphere models (TBMs) have become an integral tool for extrapolating local observations and understanding of land–atmosphere carbon exchange to larger regions. The North American Carbon Program (NACP) Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal model intercomparison and evaluation effort focused on improving the diagnosis and attribution of carbon exchange at regional and global scales. MsTMIP builds upon current and past synthesis activities, and has a unique framework designed to isolate, interpret, and inform understanding of how model structural differences impact estimates of carbon uptake and release. Here we provide an overview of the MsTMIP effort and describe how the MsTMIP experimental design enables the assessment and quantification of TBM structural uncertainty. Model structure refers to the types of processes considered (e.g., nutrient cycling, disturbance, lateral transport of carbon), and how these processes are represented (e.g., photosynthetic formulation, temperature sensitivity, respiration) in the models. By prescribing a common experimental protocol with standard spin-up procedures and driver data sets, we isolate any biases and variability in TBM estimates of regional and global carbon budgets resulting from differences in the models themselves (i.e., model structure) and model-specific parameter values. An initial intercomparison of model structural differences is represented using hierarchical cluster diagrams (a.k.a. dendrograms), which highlight similarities and differences in how models account for carbon cycle, vegetation, energy, and nitrogen cycle dynamics. We show that, despite the standardized protocol used to derive initial conditions, models show a high degree of variation for GPP, total living biomass, and total soil carbon, underscoring the influence of differences in model structure and parameterization on model estimates.

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.

The Unified North American Soil Map (UNASM) was developed to provide more accurate regional soil information for terrestrial biosphere modeling. The UNASM combines information from state-of-the-art US STATSGO2 and Soil Landscape of Canada (SLCs) databases. The area not covered by these datasets is filled by using the Harmonized World Soil Database version 1.21 (HWSD1.21). The UNASM contains maximum soil depth derived from the data source as well as seven soil attributes (including sand, silt, and clay content, gravel content, organic carbon content, pH, and bulk density) for the topsoil layer (0–30 cm) and the subsoil layer (30–100 cm), respectively, of the spatial resolution of 0.25 degrees in latitude and longitude. There are pronounced differences in the spatial distributions of soil properties and soil organic carbon between UNASM and HWSD, but the UNASM overall provides more detailed and higher-quality information particularly in Alaska and central Canada. To provide more accurate and up-to-date estimate of soil organic carbon stock in North America, we incorporated Northern Circumpolar Soil Carbon Database (NCSCD) into the UNASM. The estimate of total soil organic carbon mass in the upper 100 cm soil profile based on the improved UNASM is 365.96 Pg, of which 23.1% is under trees, 14.1% is in shrubland, and 4.6% is in grassland and cropland. This UNASM data will provide a resource for use in terrestrial ecosystem modeling both for input of soil characteristics and for benchmarking model output.

The El Niño Southern Oscillation is the dominant year-to-year mode of global climate variability. El Niño effects on terrestrial carbon cycling are mediated by associated climate anomalies, primarily drought, influencing fire emissions and biotic net ecosystem exchange (NEE). Here we evaluate whether El Niño produces a consistent response from the global carbon cycle. We apply a novel bottom-up approach to estimating global NEE anomalies based on FLUXNET data using land cover maps and weather reanalysis. We analyze 13 years (1997–2009) of globally gridded observational NEE anomalies derived from eddy covariance flux data, remotely-sensed fire emissions at the monthly time step, and NEE estimated from an atmospheric transport inversion. We evaluate the overall consistency of biospheric response to El Niño and, more generally, the link between global CO2 flux anomalies and El Niño-induced drought. Our findings, which are robust relative to uncertainty in both methods and time-lags in response, indicate that each event has a different spatial signature with only limited spatial coherence in Amazônia, Australia and southern Africa. For most regions, the sign of response changed across El Niño events. Biotic NEE anomalies, across 5 El Niño events, ranged from –1.34 to +0.98 Pg C yr−1, whereas fire emissions anomalies were generally smaller in magnitude (ranging from –0.49 to +0.53 Pg C yr−1). Overall drought does not appear to impose consistent terrestrial CO2 flux anomalies during El Niños, finding large variation in globally integrated responses from –1.15 to +0.49 Pg C yr−1. Despite the significant correlation between the CO2 flux and El Niño indices, we find that El Niño events have, when globally integrated, both enhanced and weakened terrestrial sink strength, with no consistent response across events.

Croplands are man-made ecosystems that have high net primary productivity during the growing season of crops, thus impacting carbon and other exchanges with the atmosphere. These exchanges play a major role in nutrient cycling and climate change related issues. An accurate representation of crop phenology and physiology is important in land-atmosphere carbon models being used to predict these exchanges. To better estimate time-varying exchanges of carbon, water, and energy of croplands using the Simple Biosphere (SiB) model, we developed crop-specific phenology models and coupled them to SiB. The coupled SiB-phenology model (SiBcrop) replaces remotely-sensed NDVI information, on which SiB originally relied for deriving Leaf Area Index (LAI) and the fraction of Photosynthetically Active Radiation (fPAR) for estimating carbon dynamics. The use of the new phenology scheme within SiB substantially improved the prediction of LAI and carbon fluxes for maize, soybean, and wheat crops, as compared with the observed data at several AmeriFlux eddy covariance flux tower sites in the US mid continent region. SiBcrop better predicted the onset and end of the growing season, harvest, interannual variability associated with crop rotation, day time carbon uptake (especially for maize) and day to day variability in carbon exchange. Biomass predicted by SiBcrop had good agreement with the observed biomass at field sites. In the future, we will predict fine resolution regional scale carbon and other exchanges by coupling SiBcrop with RAMS (the Regional Atmospheric Modeling System).

No abstract available.