scholarly journals Study of the footprints of short-term variation in XCO<sub>2</sub> observed by TCCON sites using NIES and FLEXPART atmospheric transport models

Author(s):  
D. A. Belikov ◽  
S. Maksyutov ◽  
A. Ganshin ◽  
R. Zhuravlev ◽  
N. M. Deutscher ◽  
...  

Abstract. The Total Carbon Column Observing Network (TCCON) is a network of ground-based Fourier Transform Spectrometers (FTS) that record near-infrared (NIR) spectra of the Sun. From these spectra, accurate and precise observations of CO2 column-averaged dry-air mole fraction (denoted XCO2) are retrieved. TCCON FTS observations have previously been used to validate satellite estimations of XCO2; however, our knowledge of the short-term spatial and temporal variations in XCO2 surrounding the TCCON sites is limited. In this work, we use the National Institute for Environmental Studies (NIES) Eulerian three-dimensional transport model and the FLEXPART (FLEXible PARTicle) Lagrangian Particle Dispersion Model (LPDM) to determine the footprints of short-term variations in XCO2 observed by operational, past, future, and possible TCCON sites. We propose a footprint-based method for the colocation of satellite and TCCON XCO2 observations, and estimate the performance of the method using the NIES model and five GOSAT XCO2 product datasets. Comparison of the proposed approach with a standard geographic method shows higher number of colocation points and average bias reduction up to 0.15 ppm for a subset of 16 stations for the period from January 2010 to January 2014. Case studies of the Darwin and La Réunion sites reveal that when the footprint area is rather curved, non-uniform and significantly different from a geographical rectangular area, the differences between these approaches are more noticeable. This emphasizes that the colocation is sensitive to local meteorological conditions and flux distributions.

2017 ◽  
Vol 17 (1) ◽  
pp. 143-157 ◽  
Author(s):  
Dmitry A. Belikov ◽  
Shamil Maksyutov ◽  
Alexander Ganshin ◽  
Ruslan Zhuravlev ◽  
Nicholas M. Deutscher ◽  
...  

Abstract. The Total Carbon Column Observing Network (TCCON) is a network of ground-based Fourier transform spectrometers (FTSs) that record near-infrared (NIR) spectra of the sun. From these spectra, accurate and precise observations of CO2 column-averaged dry-air mole fractions (denoted XCO2) are retrieved. TCCON FTS observations have previously been used to validate satellite estimations of XCO2; however, our knowledge of the short-term spatial and temporal variations in XCO2 surrounding the TCCON sites is limited. In this work, we use the National Institute for Environmental Studies (NIES) Eulerian three-dimensional transport model and the FLEXPART (FLEXible PARTicle dispersion model) Lagrangian particle dispersion model (LPDM) to determine the footprints of short-term variations in XCO2 observed by operational, past, future and possible TCCON sites. We propose a footprint-based method for the collocation of satellite and TCCON XCO2 observations and estimate the performance of the method using the NIES model and five GOSAT (Greenhouse Gases Observing Satellite) XCO2 product data sets. Comparison of the proposed approach with a standard geographic method shows a higher number of collocation points and an average bias reduction up to 0.15 ppm for a subset of 16 stations for the period from January 2010 to January 2014. Case studies of the Darwin and Reunion Island sites reveal that when the footprint area is rather curved, non-uniform and significantly different from a geographical rectangular area, the differences between these approaches are more noticeable. This emphasises that the collocation is sensitive to local meteorological conditions and flux distributions.


2016 ◽  
Vol 16 (14) ◽  
pp. 9149-9161 ◽  
Author(s):  
Misa Ishizawa ◽  
Osamu Uchino ◽  
Isamu Morino ◽  
Makoto Inoue ◽  
Yukio Yoshida ◽  
...  

Abstract. Extremely high levels of column-averaged dry-air mole fractions of atmospheric methane (XCH4) were detected in August and September 2013 over northeast Asia (∼  20 ppb above the averaged summertime XCH4 over 2009–2012, after removing a long-term trend), as being retrieved from the Short-Wavelength InfraRed (SWIR) spectral data observed with the Thermal And Near-infrared Sensor for carbon Observation – Fourier Transform Spectrometer (TANSO-FTS) onboard Greenhouse Gases Observing Satellite (GOSAT). Similar enhancements of XCH4 were also observed by the ground-based measurements at two Total Carbon Column Observing Network (TCCON) sites in Japan. The analysis of surface CH4 concentrations observed at three monitoring sites around the Japan archipelago suggest that the extreme increase of XCH4 has occurred in a limited area. The model analysis was conducted to investigate this anomalously high XCH4 event, using an atmospheric transport model. The results indicate that the extreme increase of XCH4 is attributed to the anomalous atmospheric pressure pattern over East Asia during the summer of 2013, which effectively transported the CH4-rich air to Japan from the strong CH4 source areas in east China. The two Japanese TCCON sites, ∼  1000 km east–west apart each other, coincidentally located along the substantially CH4-rich air flow from east China. This analysis demonstrates the capability of GOSAT to monitor an XCH4 event on a synoptic scale. We anticipate that the synoptic information of XCH4 from GOSAT data contributes to improve our understanding of regional carbon cycle and the regional flux estimation.


2015 ◽  
Vol 15 (19) ◽  
pp. 11147-11164 ◽  
Author(s):  
B. Oney ◽  
S. Henne ◽  
N. Gruber ◽  
M. Leuenberger ◽  
I. Bamberger ◽  
...  

Abstract. We describe a new rural network of four densely placed (< 100 km apart), continuous atmospheric carbon (CO2, CH4, and CO) measurement sites in north-central Switzerland and analyze its suitability for regional-scale (~ 100–500 km) carbon flux studies. We characterize each site for the period from March 2013 to February 2014 by analyzing surrounding land cover, observed local meteorology, and sensitivity to surface fluxes, as simulated with the Lagrangian particle dispersion model FLEXPART-COSMO (FLEXible PARTicle dispersion model-Consortium for Small-Scale Modeling). The Beromünster measurements are made on a tall tower (212 m) located on a gentle hill. At Beromünster, regional CO2 signals (measurement minus background) vary diurnally from −4 to +4 ppmv, on average, and are simulated to come from nearly the entire Swiss Plateau, where 50 % of surface influence is simulated to be within 130–260 km distance. The Früebüel site measurements are made 4 m above ground on the flank of a gently sloping mountain. Nearby (< 50 km) pasture and forest fluxes exert the most simulated surface influence, except during convective summertime days when the site is mainly influenced by the eastern Swiss Plateau, which results in summertime regional CO2 signals varying diurnally from −5 to +12 ppmv and elevated summer daytime CH4 signals (+30 ppbv above other sites). The Gimmiz site measurements are made on a small tower (32 m) in flat terrain. Here, strong summertime regional signals (−5 to +60 ppmv CO2) stem from large, nearby (< 50 km) crop and anthropogenic fluxes of the Seeland region, except during warm or windy days when simulated surface influence is of regional scale (< 250 km). The Lägern-Hochwacht measurements are made on a small tower (32 m) on top of the steep Lägern crest, where simulated surface influence is typically of regional scale (130–300 km) causing summertime regional signals to vary from −5 to +8 ppmv CO2. Here, considerable anthropogenic influence from the nearby industrialized region near Zurich causes the average wintertime regional CO2 signals to be 5 ppmv above the regional signals simultaneously measured at the Früebüel site. We find that the suitability of the data sets from our current observation network for regional carbon budgeting studies largely depends on the ability of the high-resolution (2 km) atmospheric transport model to correctly capture the temporal dynamics of the stratification of the lower atmosphere at the different sites. The current version of the atmospheric transport model captures these dynamics well, but it clearly reaches its limits at the sites in steep topography and at the sites that generally remain in the surface layer. Trace gas transport and inverse modeling studies will be necessary to determine the impact of these limitations on our ability to derive reliable regional-scale carbon flux estimates in the complex Swiss landscape.


2007 ◽  
Vol 7 (7) ◽  
pp. 1851-1868 ◽  
Author(s):  
G. Pérez-Landa ◽  
P. Ciais ◽  
G. Gangoiti ◽  
J. L. Palau ◽  
A. Carrara ◽  
...  

Abstract. Vertical profiles of CO2 concentration were collected during an intensive summer campaign in a coastal complex-terrain region within the frame of the European Project RECAB (Regional Assessment and Modelling of the Carbon Balance in Europe). The region presents marked diurnal mesoscale circulation patterns. These circulations result in a specific coupling between atmospherically transported CO2 and its surface fluxes. To understand the effects of this interaction on the spatial variability of the observed CO2 concentrations, we applied a high-resolution transport simulation to an idealized model of land biotic fluxes. The regional Net Ecosystem Exchange fluxes were extrapolated for different land-use classes by using a set of eddy-covariance measurements. The atmospheric transport model is a Lagrangian particle dispersion model, driven by a simulation of the RAMS mesoscale model. Our simulations were able to successfully reproduce some of the processes controlling the mesoscale transport of CO2. A semi-quantitative comparison between simulations and data allowed us to characterize how the coupling between mesoscale transport and surface fluxes produced CO2 spatial gradients in the domain. Temporal averages in the simulated CO2 field show a covariance between flux and transport consisting of: 1) horizontally, a CO2 deficit over land, mirrored by a CO2 excess over the sea and 2) vertically, the prevalence of a mean CO2 depletion between 500 and 2000 m, and a permanent build-up of CO2 in the lower levels.


2016 ◽  
Vol 9 (2) ◽  
pp. 749-764 ◽  
Author(s):  
Dmitry A. Belikov ◽  
Shamil Maksyutov ◽  
Alexey Yaremchuk ◽  
Alexander Ganshin ◽  
Thomas Kaminski ◽  
...  

Abstract. We present the development of the Adjoint of the Global Eulerian–Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian Particle Dispersion Model (LPDM). The forward tangent linear and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com), with additional manual pre- and post-processing aimed at improving transparency and clarity of the code and optimizing the performance of the computing, including MPI (Message Passing Interface). The Lagrangian component did not require any code modification, as LPDMs are self-adjoint and track a significant number of particles backward in time in order to calculate the sensitivity of the observations to the neighboring emission areas. The constructed Eulerian adjoint was coupled with the Lagrangian component at a time boundary in the global domain. The simulations presented in this work were performed using the A-GELCA model in forward and adjoint modes. The forward simulation shows that the coupled model improves reproduction of the seasonal cycle and short-term variability of CO2. Mean bias and standard deviation for five of the six Siberian sites considered decrease roughly by 1 ppm when using the coupled model. The adjoint of the Eulerian model was shown, through several numerical tests, to be very accurate (within machine epsilon with mismatch around to ±6 e−14) compared to direct forward sensitivity calculations. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation. A-GELCA will be incorporated into a variational inversion system designed to optimize surface fluxes of greenhouse gases.


2015 ◽  
Vol 15 (17) ◽  
pp. 24995-25020 ◽  
Author(s):  
M. Ishizawa ◽  
O. Uchino ◽  
I. Morino ◽  
M. Inoue ◽  
Y. Yoshida ◽  
...  

Abstract. Extremely high levels of column-averaged dry-air mole fractions of atmospheric methane (XCH4) were detected in August and September 2013 over Northeast Asia (~ 20 ppb above the averaged summertime XCH4 over 2009–2012, after removing a long-term trend), as being retrieved from the Short-Wavelength InfraRed (SWIR) spectral data observed with the Thermal And Near-infrared Sensor for carbon Observation-Fourier Transform Spectrometer (TANSO-FTS) onboard Greenhouse Gases Observing Satellite (GOSAT). Similar enhancements of XCH4 were also observed by the ground-based measurements at two Total Carbon Column Observing Network (TCCON) sites in Japan. The analysis of surface CH4 concentrations observed at three monitoring sites around the Japan islands suggest that the extreme increase of XCH4 has occurred in a limited area. The model analysis was conducted to investigate this anomalously high XCH4 event, using an atmospheric transport model. The results indicate that the extreme increase of XCH4 is attributed to the anomalous atmospheric pressure pattern over East Asia during the summer of 2013, which effectively transported the CH4-rich air to Japan from the strong CH4 source areas in East China. The two Japanese TCCON sites, ~ 1000 km east-west apart each other, coincidentally located along the substantially CH4-rich air flow from East China. The GOSAT orbiting with three-day recurrence successfully observed the synoptic-scale XCH4 enhancement in the comparable accuracy to the TCCON data. This analysis demonstrates the capability of GOSAT to monitor an XCH4 event on a synoptic scale.


2015 ◽  
Vol 8 (7) ◽  
pp. 5983-6019
Author(s):  
D. A. Belikov ◽  
S. Maksyutov ◽  
A. Yaremchuk ◽  
A. Ganshin ◽  
T. Kaminski ◽  
...  

Abstract. We present the development of the Adjoint of the Global Eulerian–Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian plume diffusion model (LPDM). The tangent and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com), with additional manual pre- and post-processing aimed at improving the performance of the computing, including MPI (Message Passing Interface). As results, the adjoint of Eulerian model is discrete. Construction of the adjoint of the Lagrangian component did not require any code modification, as LPDMs are able to track a significant number of particles back in time and thereby calculate the sensitivity of observations to the neighboring emissions areas. Eulerian and Lagrangian adjoint components were coupled at the time boundary in the global domain.The results are verified using a series of test experiments. The forward simulation shown the coupled model is effective in reproducing the seasonal cycle and short-term variability of CO2 even in the case of multiple limiting factors, such as high uncertainty of fluxes and the low resolution of the Eulerian model. The adjoint model demonstrates the high accuracy compared to direct forward sensitivity calculations and fast performance. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation.


2007 ◽  
Vol 4 (4) ◽  
pp. 2877-2914 ◽  
Author(s):  
C. D. Nevison ◽  
N. M. Mahowald ◽  
S. C. Doney ◽  
I. D. Lima

Abstract. A three dimensional, time-evolving field of atmospheric potential oxygen (APO ~ O2/N2 + CO2) is estimated using surface O2, N2 and CO2 fluxes from the WHOI ocean ecosystem model to force the MATCH atmospheric transport model. Land and fossil carbon fluxes are also run in MATCH and translated into O2 tracers using assumed O2:CO2 stoichiometries. The model seasonal cycles in APO agree well with the observed cycles at 13 global monitoring stations, with agreement helped by the inclusion of oceanic CO2 in the APO calculation. The model latitudinal gradient in APO is strongly influenced by seasonal rectifier effects in atmospheric transport, which appear at least partly unrealistic based on comparison to observations. An analysis of the APO vs.~CO2 method for partitioning land and ocean carbon sinks is performed in the controlled context of the MATCH simulation, in which the true surface carbon and oxygen fluxes are known exactly. This analysis suggests uncertainty ranging up to ±0.2 PgC in the inferred sinks due to transport-induced variability. It also shows that interannual variability in oceanic O2 fluxes can cause increasingly large error in the sink partitioning when the method is applied over increasingly short timescales. However, when decadal or longer averages are used, the variability in the oceanic O2 flux is relatively small, allowing carbon sinks to be partitioned to within a standard deviation of 0.1 Pg C/yr of the true values, provided one has an accurate estimate of long-term mean O2 outgassing.


2008 ◽  
Vol 5 (3) ◽  
pp. 875-889 ◽  
Author(s):  
C. D. Nevison ◽  
N. M. Mahowald ◽  
S. C. Doney ◽  
I. D. Lima ◽  
N. Cassar

Abstract. A three dimensional, time-evolving field of atmospheric potential oxygen (APO ~O2/N2+CO2) was estimated using surface O2, N2 and CO2 fluxes from the WHOI ocean ecosystem model to force the MATCH atmospheric transport model. Land and fossil carbon fluxes were also run in MATCH and translated into O2 tracers using assumed O2:CO2 stoichiometries. The modeled seasonal cycles in APO agree well with the observed cycles at 13 global monitoring stations, with agreement helped by including oceanic CO2 in the APO calculation. The modeled latitudinal gradient in APO is strongly influenced by seasonal rectifier effects in atmospheric transport. An analysis of the APO-vs.-CO2 mass-balance method for partitioning land and ocean carbon sinks was performed in the controlled context of the MATCH simulation, in which the true surface carbon and oxygen fluxes were known exactly. This analysis suggests uncertainty of up to ±0.2 PgC in the inferred sinks due to variability associated with sparse atmospheric sampling. It also shows that interannual variability in oceanic O2 fluxes can cause large errors in the sink partitioning when the method is applied over short timescales. However, when decadal or longer averages are used, the variability in the oceanic O2 flux is relatively small, allowing carbon sinks to be partitioned to within a standard deviation of 0.1 Pg C/yr of the true values, provided one has an accurate estimate of long-term mean O2 outgassing.


2020 ◽  
Author(s):  
Yuanhong Zhao ◽  
Marielle Saunois ◽  
Philippe Bousquet ◽  
Xin Lin ◽  
Antoine Berchet ◽  
...  

Abstract. The hydroxyl radical (OH), which is the dominant sink of methane (CH4), plays a key role to close the global methane budget. Previous research that assessed the impact of OH changes on the CH4 budget mostly relied on box modeling inversions with a very simplified atmospheric transport and no representation of the heterogeneous spatial distribution of OH radicals. Here using a variational Bayesian inversion framework and a 3D chemical transport model, LMDz, combined with 10 different OH fields derived from chemistry-climate models (CCMI experiment), we evaluate the influence of OH burden, spatial distribution, and temporal variations on the global CH4 budget. The global tropospheric mean CH4-reaction-weighted [OH] ([OH]GM-CH4) ranges 10.3–16.3 × 105 molec cm−3 across 10 OH fields during the early 2000s, resulting in inversion-based global CH4 emissions between 518 and 757 Tg yr−1. The uncertainties in CH4 inversions induced by the different OH fields are comparable to, or even larger than the uncertainty typically given by bottom-up and top-down estimates. Based on the LMDz inversions, we estimate that a 1 %-increase in OH burden leads to an increase of 4 Tg yr−1 in the estimate of global methane emissions, which is about 25 % smaller than what is estimated by box-models. The uncertainties in emissions induced by OH are largest over South America, corresponding to large inter-model differences of [OH] in this region. From the early to the late 2000s, the optimized CH4 emissions increased by 21.9 ± 5.7 Tg yr−1 (16.6–30.0 Tg yr−1), of which ~ 25 % (on average) is contributed by −0.5 to +1.8 % increase in OH burden. If the CCMI models represent the OH trend properly over the 2000s, our results show that a higher increasing trend of CH4 emissions is needed to match the CH4 observations compared to the CH4 emission trend derived using constant OH. This study strengthens the importance to reach a better representation of OH burden and of OH spatial and temporal distributions to reduce the uncertainties on the global CH4 budget.


Sign in / Sign up

Export Citation Format

Share Document