Estimating surface CO2 fluxes from space-borne dry air mole fraction observations using an Ensemble Kalman filter

We have developed an ensemble Kalman Filter (EnKF) to estimate 8-day regional surface fluxes of CO2 from space-borne CO2 dry-air mole fraction observations (XCO2) and evaluate the approach using a series of synthetic experiments, in preparation for data from the NASA Orbiting Carbon Observatory (OCO). The 32-day duty cycle of OCO alternates every 16 days between nadir and glint measurements of backscattered solar radiation at short-wave infrared wavelengths. The EnKF uses an ensemble of states to represent the error covariances to estimate 8-day CO2 surface fluxes over 144 geographical regions. We use a 12×8-day lag window, recognising that XCO2 measurements include surface flux information from prior time windows. The observation operator that relates surface CO2 fluxes to atmospheric distributions of XCO2 includes: a) the GEOS-Chem transport model that relates surface fluxes to global 3-D distributions of CO2 concentrations, which are sampled at the time and location of OCO measurements that are cloud-free and have aerosol optical depths <0.3; and b) scene-dependent averaging kernels that relate the CO2 profiles to XCO2, accounting for differences between nadir and glint measurements, and the associated scene-dependent observation errors. We show that OCO XCO2 measurements significantly reduce the uncertainties of surface CO2 flux estimates. Glint measurements are generally better at constraining ocean CO2 flux estimates. Nadir XCO2 measurements over the terrestrial tropics are sparse throughout the year because of either clouds or smoke. Glint measurements provide the most effective constraint for estimating tropical terrestrial CO2 fluxes by accurately sampling fresh continental outflow over neighbouring oceans. We also present results from sensitivity experiments that investigate how flux estimates change with 1) bias and unbiased errors, 2) alternative duty cycles, 3) measurement density and correlations, 4) the spatial resolution of estimated flux estimates, and 5) reducing the length of the lag window and the size of the ensemble. At the revision stage of this manuscript, the OCO instrument failed to reach its orbit after it was launched on 24 February 2009. The EnKF formulation presented here is also applicable to GOSAT measurements of CO2 and CH4.

Short-term temporal changes of bare soil CO2 fluxes after tillage described by first-order decay models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Hurricane disturbance and recovery of energy balance, CO2 fluxes and canopy structure in a mangrove forest of the Florida Everglades

Eddy covariance (EC) estimates of carbon dioxide (CO2) fluxes and energy balance are examined to investigate the functional responses of a mature mangrove forest to a disturbance generated by Hurricane Wilma on October 24, 2005 in the Florida Everglades. At the EC site, high winds from the hurricane caused nearly 100% defoliation in the upper canopy and widespread tree mortality. Soil temperatures down to -50 cm increased, and air temperature lapse rates within the forest canopy switched from statically stable to statically unstable conditions following the disturbance. Unstable conditions allowed more efficient transport of water vapor and CO2 from the surface up to the upper canopy layer. Significant increases in latent heat fluxes (LE) and nighttime net ecosystem exchange (NEE) were also observed and sensible heat fluxes (H) as a proportion of net radiation decreased significantly in response to the disturbance. Many of these impacts persisted through much of the study period through 2009. However, local albedo and MODIS (Moderate Resolution Imaging Spectro-radiometer) data (the Enhanced Vegetation Index) indicated a substantial proportion of active leaf area recovered before the EC measurements began 1 year after the storm. Observed changes in the vertical distribution and the degree of clumping in newly emerged leaves may have affected the energy balance. Direct comparisons of daytime NEE values from before the storm and after our measurements resumed did not show substantial or consistent differences that could be attributed to the disturbance. Regression analyses on seasonal time scales were required to differentiate the storm's impact on monthly average daytime NEE from the changes caused by interannual variability in other environmental drivers. The effects of the storm were apparent on annual time scales, and CO2 uptake remained approximately 250 g C m-2 yr-1 lower in 2009 compared to the average annual values measured in 2004-2005. Dry season CO2 uptake was relatively more affected by the disturbance than wet season values. Complex leaf regeneration dynamics on damaged trees during ecosystem recovery are hypothesized to lead to the variable dry versus wet season impacts on daytime NEE. In contrast, nighttime CO2 release (i.e., nighttime respiration) was consistently and significantly greater, possibly as a result of the enhanced decomposition of litter and coarse woody debris generated by the storm, and this effect was most apparent in the wet seasons compared to the dry seasons. The largest pre- and post-storm differences in NEE coincided roughly with the delayed peak in cumulative mortality of stems in 2007-2008. Across the hurricane-impacted region, cumulative tree mortality rates were also closely correlated with declines in peat surface elevation. Mangrove forest-atmosphere interactions are interpreted with respect to the damage and recovery of stand dynamics and soil accretion processes following the hurricane.

Energy and CO2 fluxes in temperate deciduous broadleaf forests: using local add-covariance data for a global optimization of ORCHIDEE model

Non peer reviewed

The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

Daytime CO2 urban surface fluxes from airborne measurements, eddy-covariance observations and emissions inventory in Greater London

Airborne measurements within the urban mixing layer (360 m) over Greater London are used to quantify CO2 emissions at the meso-scale. Daytime CO2 fluxes, calculated by the Integrative Mass Boundary Layer (IMBL) method, ranged from 46 to 104 μmol CO2 m−2 s−1 for four days in October 2011. The day-to-day variability of IMBL fluxes is at the same order of magnitude as for surface eddy-covariance fluxes observed in central London. Compared to fluxes derived from emissions inventory, the IMBL method gives both lower (by −37%) and higher (by 19%) estimates. The sources of uncertainty of applying the IMBL method in urban areas are discussed and guidance for future studies is given.

Fluxes of CO2 above a sugarcane plantation in Brazil

Fluxes of CO2 were measured above a sugarcane plantation using the eddy-covariance method covering two growth cycles, representing the second and third re-growth (ratoons) harvested with stubble burning. The total net ecosystem exchange (NEE) in the first cycle (second ratoon, 393 days long) was −1964 ± 44 g C m−2; the gross ecosystem productivity (GEP) was 3612 ± 46 g C m−2 and the ecosystem respiration (RE) was 1648 ± 14 g C m−2. The NEE and GEP totals in the second cycle (third ratoon, 374 days long) decreased 51% and 25%, respectively and RE increased 7%. Accounting for the carbon emitted during biomass burning and the removal of stalks at harvest, net ecosystem carbon balance (NECB) totals were 102 ± 130 g C m−2 and 403 ± 84 g C m−2 in each cycle respectively. Thus the sugarcane agrosystem was approximately carbon neutral in the second ratoon. Yield in stalks fresh weight (SFW) attained the regional average (8.3 kg SFW m−2). Although it was a carbon source to the atmosphere, observed productivity (6.2 kg SFW m−2) of the third ratoon was 19% lower than the regional average due to the lower water availability observed during the initial 120 days of re-growth. However, the overall water use efficiency (WUE) achieved in the first cycle (4.3 g C kg−1 H2O) decreased only 5% in the second cycle. © 2013 Elsevier B.V. All rights reserved

Air-sea fluxes of CO2 and CH4 from the Penlee Point Atmospheric Observatory on the south-west coast of the UK

We present air–sea fluxes of carbon dioxide (CO2), methane (CH4), momentum, and sensible heat measured by the eddy covariance method from the recently established Penlee Point Atmospheric Observatory (PPAO) on the south-west coast of the United Kingdom. Measurements from the south-westerly direction (open water sector) were made at three different sampling heights (approximately 15, 18, and 27m above mean sea level, a.m.s.l.), each from a different period during 2014–2015. At sampling heights ≥18ma.m.s.l., measured fluxes of momentum and sensible heat demonstrate reasonable (≤ ±20% in the mean) agreement with transfer rates over the open ocean. This confirms the suitability of PPAO for air–sea exchange measurements in shelf regions. Covariance air–sea CO2 fluxes demonstrate high temporal variability. Air-to-sea transport of CO2 declined from spring to summer in both years, coinciding with the breakdown of the spring phytoplankton bloom. We report, to the best of our knowledge, the first successful eddy covariance measurements of CH4 emissions from a marine environment. Higher sea-to-air CH4 fluxes were observed during rising tides (20±3; 38±3; 29±6 μmolem-2 d-1 at 15, 18, 27ma.m.s.l.) than during falling tides (14±2; 22±2; 21±5 μmolem-2 d-1), consistent with an elevated CH4 source from an estuarine outflow driven by local tidal circulation. These fluxes are a few times higher than the predicted CH4 emissions over the open ocean and are significantly lower than estimates from other aquatic CH4 hotspots (e.g. polar regions, freshwater). Finally, we found the detection limit of the air–sea CH4 flux by eddy covariance to be 20 μmolem-2 d-1 over hourly timescales (4 μmolem-2 d-1 over 24 h).

Air-sea fluxes of CO2 and CH4 from the Penlee Point Atmospheric Observatory on the south-west coast of the UK

We present air–sea fluxes of carbon dioxide (CO2), methane (CH4), momentum, and sensible heat measured by the eddy covariance method from the recently established Penlee Point Atmospheric Observatory (PPAO) on the south-west coast of the United Kingdom. Measurements from the south-westerly direction (open water sector) were made at three different sampling heights (approximately 15, 18, and 27m above mean sea level, a.m.s.l.), each from a different period during 2014–2015. At sampling heights ≥18ma.m.s.l., measured fluxes of momentum and sensible heat demonstrate reasonable (≤ ±20% in the mean) agreement with transfer rates over the open ocean. This confirms the suitability of PPAO for air–sea exchange measurements in shelf regions. Covariance air–sea CO2 fluxes demonstrate high temporal variability. Air-to-sea transport of CO2 declined from spring to summer in both years, coinciding with the breakdown of the spring phytoplankton bloom. We report, to the best of our knowledge, the first successful eddy covariance measurements of CH4 emissions from a marine environment. Higher sea-to-air CH4 fluxes were observed during rising tides (20±3; 38±3; 29±6 μmolem-2 d-1 at 15, 18, 27ma.m.s.l.) than during falling tides (14±2; 22±2; 21±5 μmolem-2 d-1), consistent with an elevated CH4 source from an estuarine outflow driven by local tidal circulation. These fluxes are a few times higher than the predicted CH4 emissions over the open ocean and are significantly lower than estimates from other aquatic CH4 hotspots (e.g. polar regions, freshwater). Finally, we found the detection limit of the air–sea CH4 flux by eddy covariance to be 20 μmolem-2 d-1 over hourly timescales (4 μmolem-2 d-1 over 24 h).

Soil CO2 flux: a method comparison of closed static chambers in a sugarcane field

A large variety of techniques have been used to measure soil CO2 released from the soil surface, and much of the variability observed between locations must be attributed to the different methods used by the investigators. Therefore, a minimum protocol of measurement procedures should be established. The objectives of this study were (a) to compare different absorption areas, concentrations and volumes of the alkali trapping solution used in closed static chambers (CSC), and (b) to compare both, the optimized alkali trapping solution and the soda-lime trapping using CSC to measure soil respiration in sugarcane areas. Three CO2 absorption areas were evaluated (7; 15 and 20 % of the soil emission area or chamber); two volumes of NaOH (40 and 80 mL) at three concentrations (0.1, 0.25 and 0.5 mol L-1). Three different types of alkaline traps were tested: (a), 80 mL of 0.5 mol L-1 NaOH in glass containers, absorption area 15 % (V0.5); (b) 40 mL of 2 mol L-1 NaOH retained in a sponge, absorption area 80 % (S2) and (c) 40 g soda lime, absorption area 15 % (SL). NaOH concentrations of 0.5 mol L-1 or lower underestimated the soil CO2-C flux or CO2 flux. The lower limit of the alkali trap absorption area should be a minimum of 20 % of the area covered by the chamber. The 2 mol L-1 NaOH solution trap (S2) was the most efficient (highest accuracy and highest CO2 fluxes) in measuring soil respiration.