Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scale

Conforming finite element methods for the clamped plate problem

Finite element methods for solving biharmonic boundary value problems are considered. The particular problem discussed is that of a clamped thin plate. This problem is reformulated in a weak, form in the Sobolev space Techniques for setting up conforming trial Functions are utilized in a Galerkin technique to produce finite element solutions. The shortcomings of various trial function formulations are discussed, and a macro—element approach to local mesh refinement using rectangular elements is given

Optimization of Nordsieck's Method for the Numerical Integration of Ordinary Differential Equations

Stability and accuracy of Nordsieck's integration method can be improved by choosing the zero-positions of the extraneous roots of the characteristic equation in a suitable way. Optimum zero-positions have been found by minimizing the lower bound of the interval of absolute stability and the coefficient of the truncation error. Various properties of the improved methods have been analysed, such as the behaviour of the high-order terms, the equivalence with multistep methods and the damping of perturbations

Carbon Sequestration in Forests: Current Rates and Attribution

No abstract availabl

Fire-induced Carbon Emissions and Regrowth Uptake in Western U.S. Forests: Documenting Variation Across Forest Types, Fire Severity, and Climate Regions

The forest area in the western United States that burns annually is increasing with warmer temperatures, more frequent droughts, and higher fuel densities. Studies that examine fire effects for regional carbon balances have tended to either focus on individual fires as examples or adopt generalizations without considering how forest type, fire severity, and regional climate influence carbon legacies. This study provides a more detailed characterization of fire effects and quantifies the full carbon impacts in relation to direct emissions, slow release of fire-killed biomass, and net carbon uptake from forest regrowth. We find important variations in fire-induced mortality and combustion across carbon pools (leaf, live wood, dead wood, litter, and duff) and across low- to high-severity classes. This corresponds to fire-induced direct emissions from 1984 to 2008 averaging 4 TgC/yr and biomass killed averaging 10.5 TgC/yr, with average burn area of 2723 sq km/yr across the western United States. These direct emission and biomass killed rates were 1.4 and 3.7 times higher, respectively, for high-severity fires than those for low-severity fires. The results show that forest regrowth varies greatly by forest type and with severity and that these factors impose a sustained carbon uptake legacy. The western U.S. fires between 1984 and 2008 imposed a net source of 12.3 TgC/yr in 2008, accounting for both direct fire emissions (9.5 TgC/yr) and heterotrophic decomposition of fire-killed biomass (6.1 TgC yr1) as well as contemporary regrowth sinks (3.3 TgC/yr). A sizeable trend exists toward increasing emissions as a larger area burns annually

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements

Hydroclimatic Controls over Global Variations in Phenology and Carbon Flux

The connection between phenological and hydroclimatological variations are quantified through joint analyses of global NDVI, LAI, and precipitation datasets. The global distributions of both NDVI and LAI in the warm season are strongly controlled by three quantities: mean annual precipitation, the standard deviation of annual precipitation, and Budyko's index of dryness. Upon demonstrating that these same basic (if biased) relationships are produced by a dynamic vegetation model (the dynamic vegetation and carbon storage components of the NCAR Community Land Model version 4 combined with the water and energy balance framework of the Catchment Land Surface Model of the NASA Global Modeling and Assimilation Office), we use the model to perform a sensitivity study focusing on how phenology and carbon flux might respond to climatic change. The offline (decoupled from the atmosphere) simulations show us, for example, where on the globe a given small increment in precipitation mean or variability would have the greatest impact on carbon uptake. The analysis framework allows us in addition to quantify the degree to which climatic biases in a free-running GCM are manifested as biases in simulated phenology

Global estimation of burned area using MODIS active fire observations

We present a method for estimating monthly burned area globally at 1° spatial resolution using Terra MODIS data and ancillary vegetation cover information. Using regression trees constructed for 14 different global regions, MODIS active fire observations were calibrated to burned area estimates derived from 500-m MODIS imagery based on the assumption that burned area is proportional to counts of fire pixels. Unlike earlier methods, we allow the constant of proportionality to vary as a function of tree and herbaceous vegetation cover, and the mean size of monthly cumulative fire-pixel clusters. In areas undergoing active deforestation, we implemented a subsequent correction based on tree cover information and a simple measure of fire persistence. Regions showing good agreement between predicted and observed burned area included Boreal Asia, Central Asia, Europe, and Temperate North America, where the estimates produced by the regression trees were relatively accurate and precise. Poorest agreement was found for southern-hemisphere South America, where predicted values of burned area are both inaccurate and imprecise; this is most likely a consequence of multiple factors that include extremely persistent cloud cover, and lower quality of the 500-m burned area maps used for calibration. Application of our approach to the nine remaining regions yielded comparatively accurate, but less precise, estimates of monthly burned area. We applied the regional regression trees to the entire archive of Terra MODIS fire data to produce a monthly global burned area data set spanning late 2000 through mid-2005. Annual totals derived from this approach showed good agreement with independent annual estimates available for nine Canadian provinces, the United States, and Russia. With our data set we estimate the global annual burned area for the years 2001-2004 to vary between 2.97 million and 3.74 million km², with the maximum occurring in 2001. These coarse-resolution burned area estimates may serve as a useful interim product until long-term burned area data sets from multiple sensors and retrieval approaches become available