Atmospheric Simulations of Total Column CO2 Mole Fractions from Global to Mesoscale within the Carbon Monitoring System Flux Inversion Framework

Quantifying the uncertainty of inversion-derived CO2 surface fluxes and attributing the uncertainty to errors in either flux or atmospheric transport simulations continue to be challenges in the characterization of surface sources and sinks of carbon dioxide (CO2). Despite recent studies inferring fluxes while using higher-resolution modeling systems, the utility of regional-scale models remains unclear when compared to existing coarse-resolution global systems. Here, we present an off-line coupling of the mesoscale Weather Research and Forecasting (WRF) model to optimized biogenic CO2 fluxes and mole fractions from the global Carbon Monitoring System inversion system (CMS-Flux). The coupling framework consists of methods to constrain the mass of CO2 introduced into WRF, effectively nesting our regional domain covering most of North America (except the northern half of Canada) within the CMS global model. We test the coupling by simulating Greenhouse gases Observing SATellite (GOSAT) column-averaged dry-air mole fractions (XCO2) over North America for 2010. We find mean model-model differences in summer of ∼0.12 ppm, significantly lower than the original coupling scheme (from 0.5 to 1.5 ppm, depending on the boundary). While 85% of the XCO2 values are due to long-range transport from outside our North American domain, most of the model-model differences appear to be due to transport differences in the fraction of the troposphere below 850 hPa. Satellite data from GOSAT and tower and aircraft data are used to show that vertical transport above the Planetary Boundary Layer is responsible for significant model-model differences in the horizontal distribution of column XCO2 across North America.

Atmospheric observation-based estimation of fossil fuel CO2 emissions from regions of central and southern California

Combustion of fossil fuel is the dominant source of greenhouse gas emissions to the atmosphere in California. Here, we describe radiocarbon (¹⁴CO₂) measurements and atmospheric inverse modeling to estimate fossil fuel CO₂ (ffCO₂) emissions for 2009-2012 from a site in central California, and for June 2013-May 2014 from two sites in southern California. A priori predicted ffCO₂ mixing ratios are computed based on regional atmospheric transport model (WRF-STILT) footprints and an hourly ffCO₂ prior emission map (Vulcan 2.2). Regional inversions using observations from the central California site suggest that emissions from the San Francisco Bay Area (SFBA) are higher in winter and lower in summer. Taking all years together, the average of a total of fifteen 3-month inversions from 2009 to 2012 suggests ffCO₂ emissions from SFBA were within 6 ± 35% of the a priori estimate for that region, where posterior emission uncertainties are reported as 95% confidence intervals. Results for four 3-month inversions using measurements in Los Angeles South Coast Air Basin (SoCAB) during June 2013-May 2014 suggest that emissions in SoCAB are within 13 ± 28% of the a priori estimate for that region, with marginal detection of any seasonality. While emissions from the SFBA and SoCAB urban regions (containing ~50% of prior emissions from California) are constrained by the observations, emissions from the remaining regions are less constrained, suggesting that additional observations will be valuable to more accurately estimate total ffCO₂ emissions from California as a whole.

Interactive effects of environmental change and management strategies on regional forest carbon emissions

Climate mitigation activities in forests need to be quantified in terms of the long-term effects on forest carbon stocks, accumulation, and emissions. The impacts of future environmental change and bioenergy harvests on regional forest carbon storage have not been quantified. We conducted a comprehensive modeling study and life-cycle assessment of the impacts of projected changes in climate, CO2 concentration, and N deposition, and region-wide forest management policies on regional forest carbon fluxes. By 2100, if current management strategies continue, then the warming and CO2 fertilization effect in the given projections result in a 32-68% increase in net carbon uptake, overshadowing increased carbon emissions from projected increases in fire activity and other forest disturbance factors. To test the response to new harvesting strategies, repeated thinnings were applied in areas susceptible to fire to reduce mortality, and two clear-cut rotations were applied in productive forests to provide biomass for wood products and bioenergy. The management strategies examined here lead to long-term increased carbon emissions over current harvesting practices, although semiarid regions contribute little to the increase. The harvest rates were unsustainable. This comprehensive approach could serve as a foundation for regional place-based assessments of management effects on future carbon sequestration by forests in other locations.

Evaluating atmospheric CO2 inversions at multiple scales over a highly inventoried agricultural landscape

An intensive regional research campaign was conducted by the North American Carbon Program (NACP) in 2007 to study the carbon cycle of the highly productive agricultural regions of the Midwestern United States. Forty-five different associated projects were conducted across five US agencies over the course of nearly a decade involving hundreds of researchers. One of the primary objectives of the intensive campaign was to investigate the ability of atmospheric inversion techniques to use highly calibrated CO2 mixing ratio data to estimate CO2 flux over the major croplands of the United States by comparing the results to an inventory of CO2 fluxes. Statistics from densely monitored crop production, consisting primarily of corn and soybeans, provided the backbone of a well studied bottom-up inventory flux estimate that was used to evaluate the atmospheric inversion results. Estimates were compared to the inventory from three different inversion systems, representing spatial scales varying from high resolution mesoscale (PSU), to continental (CSU) and global (CarbonTracker), coupled to different transport models and optimization techniques. The inversion-based mean CO2 -C sink estimates were generally slightly larger, 8-20% for PSU, 10-20% for CSU, and 21% for CarbonTracker, but statistically indistinguishable, from the inventory estimate of 135 TgC. While the comparisons show that the MCI region-wide C sink is robust across inversion system and spatial scale, only the continental and mesoscale inversions were able to reproduce the spatial patterns within the region. In general, the results demonstrate that inversions can recover CO2 fluxes at sub-regional scales with a relatively high density of CO2 observations and adequate information on atmospheric transport in the region.