In situ measurement of CO2 and CH4 from aircraft over northeast China and comparison with OCO-2 data

Several satellites have been launched to monitor the increasing concentrations of greenhouse gases, especially CO2 and CH4 in the atmosphere, through backscattered hyperspectral radiance in the shortwave infrared (SWIR) band. The vertical profiles of greenhouse gases and aerosols could strongly affect the results from these instruments. To investigate the effects of the vertical distribution of CO2 on the uncertainty of SWIR satellite retrieval results, we conducted observations of the vertical profiles of CO2, CH4 and aerosol particles at 0.6–7 km above sea level using a Beechcraft King Air 350ER in Jiansanjiang (46.77 N, 131.99 E), Heilongjiang Province, northeast China, on 7–12 August 2018. The profiles from this aircraft captured a decrease in CO2 from 2 km to the minimum altitude due to the absorption of vegetation at the surface in summer. CH4 measurements showed about a 0.2 ppm increase from 2.0 to 0.6 km on 10 August, which may result from emissions from the large area of paddy fields below, and a constant mole fraction between 1.951 and 1.976 ppm was recorded at 2 km and above. Comparison of CO2 profiles from a new version of the carbon cycle data assimilation system Tan-Tracker (v1), retrievals from OCO-2 and aircraft measurements was conducted. The results from OCO-2 and the assimilation model system Tan-Tracker captured the vertical structure of CO2 above 3 km, whereas below 3 km the values from OCO-2 and the Tan-Tracker model were lower than those from in situ measurements. Column-averaged CO2 volume mole fractions calculated from in situ measurements showed biases of −4.68±0.44 ppm (−1.18±0.11%) compared to OCO-2 retrievals.


Introduction
Global warming due to greenhouse gases (GHGs) has become one of the most urgent and widely studied issues in recent years. The Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC) noted that the global average temperature has increased by 0.85 • over the period of 1880-2012. GHGs, especially the increasing CO 2 levels in the atmosphere related to anthropogenic activities, are blamed for global warming, because they absorb and emit radiant energy within the thermal infrared range. Emission of CO 2 from fossil fuel combustion and industrial processes has contributed about 78 % of the total GHG emissions increase from 1970 to 2010 (IPCC, 2014). Accurate measurement of CO 2 concentrations and their spatial and temporal variations in the atmosphere is essential for estimation of sources and sinks in regional and global models (Patra et al., 2005a, b;Zhang et al., 2008). The Global Atmospheric Watch program (https://community.wmo.int/activity-areas/gaw, last access: 19 June 2020) coordinates the systematic observation and analysis of GHGs and other trace substances, providing an important source of local and global GHG data. However, ground-based and in situ measurements near the surface can only provide information about the lower atmosphere and are insufficient for analysis of total-column GHGs, which exhibits variations in both the vertical and horizontal directions. Over the past few years, several satellites, including the Greenhouse Gases Observing Satellite (GOSAT, launched in January 2009), Second Orbiting Carbon Observatory, (OCO-2, launched in 2014) and TanSat (launched in 2016), have been launched into space to monitor CO 2 by observing backscattered hyperspectral radiance in the shortwave infrared (SWIR) wavelength, which can provide global coverage of the column-averaged dry-air mole fraction of CO 2 (X CO 2 ). Studies have shown that, given a 1-2 ppm accuracy of X CO 2 , the use of spaceborne instrument data can reduce the uncertainties in regional (8 • ×10 • footprint) estimation of CO 2 sources and sinks (Rayner and O'Brien, 2001). In addition, CO 2 vertical profiles in the 5-25 km altitude range can be obtained using limb viewing spaceborne sounders such as the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS, launched in August 2003). Foucher et al. (2009) reported the feasibility and difficulty of obtaining vertical CO 2 profiles using this method.
To validate and calibrate the X CO 2 data from satellite measurement products, the Total Carbon Column Observing Network (TCCON), a network of ground-based solar Fourier transform spectrometers operating in the SWIR spectral region, was established . Several studies have been conducted to determine the column-averaged volume mole fraction of CO 2 , CH 4 and other trace gases (X CO 2 , X CH 4 and X gas ) from TCCON data, which have shown good accuracy (Hedelius et al., 2017;Mendonca et al., 2019). In addition, commercial mobile solar-viewing near-infrared spectrometers of lower resolution than the TCCON instruments, such as the Bruker ™ EM27/SUN, show potential for measurement of X gas with an acceptable bias (Hedelius et al., 2016).
Retrieval accuracy is affected by knowledge of the vertical distribution of aerosols and CO 2 . Vertical profiles of CO 2 also affect the accuracy of estimation for regional carbon fluxes in the atmospheric transport model and can help elucidate the global carbon cycle and climate change. Many experiments have been conducted to measure the vertical profiles of CO 2 , CH 4 and other trace gases. The AirCore sampling system can be used to obtain vertical profiles of CO 2 and CH 4 from near the surface to 8-12 km with high accuracy (Karion et al., 2010;Membrive et al., 2017). Active remote sensing of atmospheric X CO 2 with the Raman lidar (light detection and ranging) technique has been developed and used to measure CO 2 vertically in the troposphere (Zhao et al., 2007;Gong et al., 2013;Han et al., 2017). CO 2 concentrations were measured at 8-12 km by Tohoku University (Sendai, Japan) through flask sampling on a commercial airliner operated by Japan Airlines (JAL) between Japan and Australia in 1984 and 1985 (Nakazawa et al., 1991). The Comprehensive Observation Network for TRace gases by AIrLiner (CONTRAIL) project installed continuous CO 2 measurement equipment on board aircraft operated by JAL for in situ measurement (Machida et al., 2008). The data for CONTRAIL are collected at altitudes between a few kilometers and 10 km, taking advantage of the frequent movement of commercial aircraft around the world. The Civil Aircraft for Remote Sensing and In Situ Measurements Based on the Instrumentation Container Concept (CARIBIC) project (Brenninkmeijer et al., 1999(Brenninkmeijer et al., , 2007 aimed to observe trace gases such as CO, O 3 and CO 2 by deploying measurement equipment in passenger aircraft. The HIAPER Pole-to-Pole Observation (HIPPO) project involved a sequence of five global aircraft measurement programs to sample the atmosphere from near the north pole to the coastal waters of Antarctica (Wofsy, 2011). Direct measurements that are independently collected from the aircraft provide validation information for satellite products. Several studies have shown that profile measurements of CO 2 and CH 4 obtained using aircraft and AirCore are useful for bias correction of both TCCON measurements Geibel et al., 2012;Hedelius et al., 2016;Messerschmidt et al., 2011) and satellite products (Araki et al., 2010;Inoue et al., 2013Inoue et al., , 2014Miyamoto et al., 2013;Frankenberg et al., 2016;Wunch et al., 2017).
Three satellites designed for CO 2 measurement, TanSat (Yang et al., 2018, GMI/GF-5 (Li et al., 2016) and GAS/FY-3D (Qi et al., 2020), were launched into space in 2016, 2017 and 2018, respectively. Measurement of profiles is crucial to further validate the retrieved hyperspectral measurements from these three satellites. Because the algorithm for satellite retrieval requires a priori profiles based on the model and in situ measurements, the lack of direct and independent airborne observations may increase the bias in the satellite results over China.
In this study, in situ aircraft-based measurements of CO 2 and CH 4 were conducted in Jiansanjiang, northeast China, in August 2018. An ultraportable greenhouse gas analyzer (UGGA; model 915-0011; Los Gatos Research, San Jose, CA, USA) was used on board the aircraft to measure the vertical mole fractions of CO 2 and CH 4 at altitudes of 0.6-7 km. Descriptions of the aircraft, the onboard instruments, and the OCO-2 and Tan-tracker data used in the paper are provided in Sect. 2. Details of the experimental site and the flight trajectory are provided in Sect. 3. The data processing method including the water vapor correction is given in Sect. 4. A comparison of the profiles obtained using aircraft with OCO-2 and the assimilation system Tan-Tracker (v1) is described in Sect. 5. The methods used to calculate X CO 2 and extrapolate in situ profiles, as well as error estimation, are discussed in Sect. 6.

Aircraft instrumentation
The aircraft used for this experiment was a Beechcraft King Air 350ER, which is a twin-turboprop aircraft designed for weather modification missions and measurement of trace gases and aerosols by the China Meteorological Administration (CMA). The cruising speed and maximum speed of the aircraft are 441 and 561 km h −1 , respectively. Temperature, wind speed, relative humidity and other meteorological data were detected and recorded by an aircraft-integrated meteorological measurement system (AIMMS-20AG) installed on the aircraft. The geolocation information including latitude, longitude, ambient pressure and height of the aircraft is also measured by the AIMMS-20AG. The relative humidity is calculated by temperature and dew point, measured by the total temperature sensor (model 102 non-de-iced, Rosemount Aerospace Inc.) and dew point hygrometer (model 137 Vigilant ™ , EdgeTech), respectively. The ultraportable greenhouse gas analyzer, UGGA (model 915-0011; Los Gatos Research), was connected to an aircraft-based impactor inlet system which consists of CVI (counterflow virtual impactor inlet system, model 1204; Brechtel Manufacturing Inc.) and ISO inlet (model 1200; Brechtel Manufacturing Inc.) in the pressurized cabin for continuous measurement of CO 2 and CH 4 . The CVI and/or ISO inlets were mounted on the top of the aircraft body as shown in Fig. 1, and the air flow rate of the inlets was kept constant by the automatic air flow controller of the inlets (aircraft-based counterflow virtual impactor inlet system CVI -model 1204, brochure; isokinetic inlet system ISO inlet -model 1200, brochure). The UGGA uses a laser absorption technology called off-axis integrated cavity output spectroscopy to determine the trace gas concentration with a high precision of < 300 ppb (CO 2 ) and < 2 ppb (CH 4 ) and a 10 s response time (UGGA user manual, model 915-0011; Los Gatos Research) and was tested and controlled in the laboratory. As shown in the in-flight schematic diagram ( Fig. 1), the external oilless diaphragm vacuum pump (F-9A 08-03, GAST) was mounted between the CVI inlet and/or the ISO inlet, and it was well-designed to keep a stable airflow with the maximum flow rate of 31.15 L min −1 . The ISO inlet was used as the aircraft passed through clouds, and the CVI inlet was used the other times. A similar system for airborne GHG measurement has been reported by O'Shea et al. (2013) and Palmer et al. (2013).
During the flight, the pressure of the sample cavity was kept constant by a small pump inside the instrument with the airflow about 0.3 L min −1 . The sample cavity temperature was also kept stable and constant by the temperature controller of the instrument. The instrument automatically recorded and saved the temperature and pressure in the cavity during operation. According to the records, the standard deviation of the cell pressure during three flights is 0.029, 0.029 and 0.033 on 7, 9 and 10 August, and the range of the cell pressure on each flight is below 0.16 hPa. For the cell temperature, the standard deviation is 0.46, 1.55 and 1.18 on each day and the range is below 3.11 • . The UGGA was calibrated against standard GHGs (provided by the National Institute of Metrology, China) before takeoff and after landing of each flight to ensure the accuracy of the data measured with the UGGA. Before this study, the GHG standard gases have been used by the CMA, Chinese Academy of Sciences, and other scientific research institutions for calibration and validation, showing that these standard gases have good performance and reliability. The standard gas we used is based on dry and clean air with known concentration values, filled in a 29.5 L aluminum alloy cylinder with salinization and other special treatment on the inner wall, traceable to the World Meteorological Organization Global Atmosphere Watch (WMO GAW) level 1 standard gas. The concentration of the CO 2 is 400.13 ppm and CH 4 is 1.867 ppm. The standard gas has been measured in the laboratory for the proportion of δ 13 C in CO 2 , and the proportion is between −8.0 ‰ and −8.2 ‰, close to the natural content, so it will not cause a significant isotopic effect on the measurement of CO 2 by the optical method and meet the requirements of standard gas (Yao et al., 2013). Just before taking off, UGGA was calibrated against standard gas, and the stability of the instrument was checked and tested again using the same standard gas of CO 2 and CH 4 immediately after landing. As shown in Fig. 2, the concentration of CO 2 and CH 4 before and after landing is stable around the values of standard gas concentration, and there was almost no drift after the flight. The precision and repeatability of the instruments are also checked and tested multiple times in the laboratory, and the results show that they are stable and good for the measurements.

Tan-Tracker and OCO-2 data
Based on the nonlinear least-squares four-dimensional variational (NLS 4D-Var) data assimilation algorithm and the Goddard Earth Observing System atmospheric chemistry transport model (GEOS-Chem), Tan-Tracker provides surface flux inversion estimates and profiles of CO 2 with 47 levels of vertical resolution from the surface to 0.03 hPa and horizontal resolution of 2.5 • × 2 • . The NLS 4D-Var assimilation model Tan-Tracker (v1) and OCO-2 X CO 2 (v9r) retrievals are used to optimize surface terrestrial ecosystem CO 2 flux and ocean CO 2 flux, while prior fossil fuel emission and fire emission remain unchanged (details of model setting and prior flux information can be found in Han and Tian, 2019).
The Orbiting Carbon Observatory-2 (OCO-2), successfully launched on 2 July 2014, obtained global measurement of CO 2 through hyperspectral measurement of reflected sunlight from Earth's atmosphere in one near-infrared (NIR) and two SWIR bands centered at 0.76, 1.61 and 2.06 µm; more details about the mission, retrieving algorithm and data char-  acteristics can be found in Crisp et al. (2008) and O'Dell et al. (2012). The uncertainty and bias of the X CO 2 products related to surface properties, aerosol and cloud, and the retrieving algorithm have been reported by Butz et al. (2009), Jung et al. (2016) and Connor et al. (2016). The OCO-2 data (V9r) including X CO 2 , the CO 2 profile and the a priori profile were used in this study.  Figure 3 shows the geolocation of the Jiansanjiang aircraft and the flight paths. The area is mostly covered with large tracts of farmland. Rice cultivation is carried out primarily in summer, and crop growth is vigorous during this period. Due to the influence of plant photosynthesis, a large amount of CO 2 uptake occurs near the surface. Three profiles were obtained between around 08:00 and 11:00 LT (GMT+8) on 7, 9 and 10 August 2018. The aircraft is designed for weather modification by the China Meteorological Administration (CMA), so the infrastructure of the aircraft and the gas flow system is also designed and completed in the USA by the team of the weather modification agency. The CMA is in charge of the flight route, and there is a chance (more flights are planned for the future) that it can carry the greenhouse gas analyzer to measure the profiles of CO 2 and CH 4 . The greenhouse gas analyzer was loaded on the aircraft and some parts of air flow arrangements were modified to better fit the requirement for profile measurement. Due to the logistical problem and the air traffic control restriction, we must fly in the morning from around 07:30 to 11:00 LT on these days to avoid obstructing civil aviation. The details of the three flights are listed in Table 1.
The flight trajectory on 7 August is shown in Fig. 4. The aircraft climbed up quickly and directly to the maximum height of about 7.1 km 30 min after taking off and then descended down step by step at about every 300 m. Since the 3-D figures on these 3 d look identical, the flight trajectories of the other 2 d (9 and 10 August) are not shown in Fig. 4 but available in the Supplement. Considering the sensitivity of the UGGA response, measurements during the ascent were discarded due to the rapid changes in air pressure, and only data collected while spiralling downward were regarded as valid and analyzed further. Data recorded below 0.6 km were also rejected because samples were easily contaminated with exhaust emissions during the slowing and descent of the aircraft before landing. The spiral descent of the aircraft lasted about 2.5 h on the 3 d, and the number of effective layers of measurements are 17, 21 and 20, respectively, on 7, 9 and 10 August.

Water vapor correction
The mole fraction of CO 2 or CH 4 measured during flight is the volume in proportion to the air containing water vapor, which cannot be directly compared with values from other data sources due to different water vapor contents of the sampled air. Therefore, the effect of water vapor is corrected and the mole fractions of CO 2 and CH 4 to dry air are given by where f gas_dry (mol mol −1 ) is the mole fraction of a gas in dry air, and f gas (mol mol −1 ) is the measured mole fraction of a gas under the real air conditions with water vapor. P H 2 O is water vapor pressure in hectopascals, which can be calculated as where e s (hPa) is the saturated water vapor pressure at the temperature T (K) at aircraft altitudes, which can be derived from the Clausius-Clapeyron equation (Wallace and Hobbs, 2006): where L v = 2.500 × 10 6 J kg −1 ; M w = 18.016, which is the molecular weight of water; R = 8.3145 J K −1 mol −1 ; and e s (hPa) is the saturated water vapor pressure at temperature T (K). Pressure p (hPa) of the ambient atmosphere is measured by the aircraft meteorology system, AIMMS-20AG, and the temperature T (K) was measured by total temperature sensor (model 102 type non-de-iced). The relative humidity (RH, %) was calculated by the dew point and temperature. The dew point data are obtained by a dew point hygrometer (model 137 Vigilant ™ , EdgeTech).

Accuracy and precision
Before the aircraft takeoff, the clocks of the UGGA, AIMMS-20AG, total temperature sensor and other instruments were adjusted to match those of the CO 2 , CH 4 and weather system measurements, synchronizing these data to the altitude and geolocation of the aircraft. The data from UGGA and synchronous meteorology measurements, including temperature, pressure and humidity of ambient atmosphere, are recorded every second and then smoothed with a 10 s running average to further remove errors caused by temporal mismatch considering the response time of the UGGA. Because the flights followed the spiral trajectories that descended approximately every 300 m, only data collected during level flight were retained and analyzed, whereas data from the descent periods were removed to avoid the effects  from vertical variations in sampling during rapid descent. The time points at the beginning and end of level flight are determined according to the altitude and its variation of the aircraft. Considering the residual time of the GHG measurement system, the data obtained 220 s from the start of the level flight are considered to be observed when the aircraft is descending rather than level, which may cause uncertainty of the measurement. Therefore, the data were kept after the level flight starting for 220 s. If the duration time of certain level flight lasted less than 220 s, the data observed during that level flight were also discarded. The instrument was calibrated against the standard gas before and after each flight. All of the measurements during the calibration process, including the standard gas used for calibration, can trace back to the WMO scale. The maximum and the average values of the difference between the standard gas and the measurement of the instrument of each day were considered for the accuracy of the aircraft data. For the precision, note that the instrument was not continuously calibrated against the standard gas during the flight. We calculated the 1 standard deviation of the data in each level flight, and the maximum of the average value of 1σ on each day is considered to be the precision of the aircraft measurement. The accuracy of CO 2 and CH 4 is below 0.66 and 0.002 ppm, 0.16 % and 0.10 % of the CO 2 and CH 4 concentration in standard gas, respectively. For precision, the 1σ value is below 0.71 and 0.0062 ppm for CO 2 and CH 4 , respectively.
5 Results and discussion 5.1 CO 2 and CH 4 profiles Figure 5 shows vertical profiles of the CO 2 and CH 4 mole fractions measured with the UGGA during the flight over Jiansanjiang, which is an agricultural area that produces a large amount of rice. The CO 2 concentration increased with height in the troposphere (Fig. 5a), which may result from CO 2 uptake by rice plants near the surface during the summer growth season. The greater increase rate of CO 2 in the lower troposphere on 7 August compared to the other 2 d was probably attributed to differing weather conditions on the 3 sampling days. It was sunny on 7 August, but it was overcast on 9 and 10 August, which may have weakened photosynthesis in rice and reduced CO 2 uptake. During all three flights, the maximum mole fraction of CO 2 reached about 418 ppm in the free troposphere at the top of the profile.
The mole fraction of CH 4 (Fig. 5b) showed a consistent decrease with increasing altitude, ranging from 1.95 to 2.10 ppm from about 2 km to near the surface, possibly as a result of CH 4 emissions from agricultural activity at the surface. CH 4 showed low variability of less than 0.5 ppm at higher altitudes, from above 2 to 7 km, indicating a wellmixed vertical structure of CH 4 in the free troposphere.
Comparing the CO 2 and CH 4 observation data, the mole fraction of CH 4 varied less than that of CO 2 from 1.5 to 2 km up to the free troposphere, with a stable value of about 1.925 ppm, indicating that CH 4 was evenly mixed at these heights and there were no obvious sources or sinks of CH 4 . CO 2 increased with altitude in the free troposphere from about 400 to 418 ppm. This increase may have been due to photosynthesis by vegetation and the large number of crops planted locally, creating a CO 2 sink at the surface and causing the CO 2 concentration to rise with height in the free troposphere. The results show that the vertical profile of CO 2 in summer increases with height in the upper troposphere, whereas that of CH 4 changed little with height and was relatively stable over the Jiansanjiang area during the experiment.

Comparison of profiles from the model and satellite product
Aircraft measurements were compared with CO 2 data obtained from OCO-2 (v9r) retrievals and the recently developed data assimilation system for the global carbon cycle, Tan-Tracker (v1) (Han and Tian, 2019). The assimilation data are collected and linearly interpolated spatially and temporally based on the geolocation of the observation site and time. Because no data were obtained from OCO-2 (v9r) over Jiansanjiang during the flight, the results of OCO-2 within 1 • ×1 • spatially at the closest time to the flight were used for comparison and were collected on 5 August. The height of the profile is available on the satellite product.
The structure of CO 2 varying with height could be roughly divided into three segments: the surface to 2, 2 to 3, and 3 to 8 km (Fig. 6). Below 2 km, CO 2 of the Tan-Tracker model is assumed to be well-mixed and uniformly distributed with height, with values ranging from 385 to 395 ppm. Therefore, the model could not reproduce the strong decrease in CO 2 from 2 km to the surface due to uptake by vegetation. From 2 to 3 km, CO 2 increased to about 400 ppm with altitude. The averaged satellite retrieval profiles correctly reproduced the decrease in CO 2 from 2 km to the surface, but the decrease rate was lower than those of in situ profiles, decreasing from 393 ppm at 2 km to 390 ppm near the surface. Flight data showed a significant CO 2 sink in this region, most notably on 7 August when it decreased from 400 ppm at 2 km to 380 ppm at 0.6 km. The impact of ground sinks was more pronounced and apparent than that from satellite inversion and model simulations, indicating that the strong variations in the lower atmosphere and planetary boundary layer (PBL) should be more carefully considered in model and retrieval algorithms. Between 2 and 4 km, aircraft profiles showed a relatively uniform mixing level of CO 2 , with roughly stable concentrations around 400 ppm.
In general, all profiles from the aircraft, satellite retrieval and model showed a similar vertical distribution trend in the troposphere above 2 km, but with differences in values. The average of the difference between OCO-2 and the aircraft profiles above 2 km is 4.22, 8.16 and −2.57 ppm on 7, 9 and 10 August, respectively. The volume mole fraction of CO 2 from both satellite and aircraft measurements indicated a CO 2 sink. GHGs profiles have rarely been observed before near the experiment site, or over the northeast of China as far as we know. The model simulations are based on data of a regional emission inventory. The accuracy of simulated profiles and concentrations near the surface over the experiment site still remains unknown. Thus continuous and regular observation of the GHGs profiles are necessary for better understanding of the regional emission amounts and the variation in the GHGs.

Comparison of X CO 2 products
The total column amount of CO 2 can be derived by integrating the CO 2 concentrations from the surface to the top of the atmosphere under the assumption of hydrostatic conditions: where VC CO 2 is the total column amount of CO 2 ; f dry CO 2 is the dry-air mole fraction (DMF) of CO 2 (mol mol −1 ); f H 2 O (p) is the aircraft profile of H 2 O (mol mol −1 ), which is measured by the onboard AIMMS system; m(p) is the mean molecular mass of wet air; g(p) is gravitational acceleration; m H 2 O = 18.02×10 −3 /N A kg molecule −1 ; m dry air = 28.964 × 10 −3 /N A kg molecule −1 ; and N A is Avogadro's constant. Data beyond the flight limits are taken from the National Centers for Environmental Prediction (NCEP) reanalysis data interpolated to the time of flight.
The column-averaged DMF ofCO 2 (X CO 2 ) from aircraft measurements was calculated based on the method of Wunch et al. (2010), which considers the average kernel in OCO-2 satellite retrievals: where a is the average kernel (Rodgers and Connor, 2003), X a CO 2 is the column-averaged DMF for the a priori profile t a , h j is the pressure weighting function of OCO-2 and t in situ is the in situ profile from aircraft measurement. Because in situ measurements available from aircraft are limited, values outside the aircraft's vertical observation range must be estimated to calculate X CO 2 . Two extrapolation methods were used to extend the profile of the aircraft measurements and then estimate the X CO 2 value of the in situ measurement respectively. (1) The unknown part of the aircraft profile was directly from the OCO-2 a priori profile.
(2) A well-mixed and constant mixing ratio of CO 2 is assumed from the surface to the lower limit of flight and from the upper limit of flight to the tropopause. The CO 2 concentrations above the tropopause were calculated with an empirical model (Toon and Wunch, 2014) which considers tropopause height as well as realistic latitude and time dependencies through curve fitting of data from high-altitude balloons, AirCore, Observations of the Middle Stratosphere balloon and aircraft. In general, the mole fraction of CO 2 decreased exponentially with height from the tropopause to upper stratosphere, and the tropopause height was obtained from NCEP reanalysis data with a 2.5 • × 2.5 • resolution, which was linearly interpolated to the geographic coordinates of Jiansanjiang. Figure 7 shows the extrapolated CO 2 profiles using method (2). X CO 2 calculated from the aircraft measurements and differences with that from OCO-2 are listed in Tables 2 and 3. The results showed that X CO 2 values of OCO-2 were lower, with an average difference of −4.68 ± 0.44 ppm (−1.18 % ± 0.11 %) and −5.09 ± 1.28 ppm (−1.28 % ± 0.32 %) by method (1) and by method (2). Table 2. X CO 2 derived from aircraft on each observation day (7, 9 and 10 August) supplemented the aircraft profile by method (1). OCO-2 (V9r) X CO 2 values were from 5 August, which was the closest time point of X CO 2 data from OCO-2 over Jiansanjiang to the observation period. Differences between aircraft X CO 2 and OCO-2 are shown in the fourth (ppm) and fifth (%) columns. The average difference and standard deviation are shown in the fifth row.

Date
Aircraft  Uncertainties induced by extrapolation of profiles outside the height limits of aircraft and errors in tropopause estimation were analyzed. Errors in extrapolation of the profile below the lower limit and above the upper limit of flight were estimated by recalculating X CO 2 after a 1 ppm positive shift in the CO 2 concentrations at these altitudes. For method (1), since the values of CO 2 mole factions of the unknown part are the same as those of the OCO-2 a priori profile, as Eq. (5) shows, no extra uncertainty would be introduced by extrapolation. But for method (2), as the profile is assumed to decrease exponentially with height above the tropopause, the height of the tropopause also introduces uncertainties for X CO 2 . Table 4 lists the errors resulting from three sources: (1) uncertainties from in situ measurement, (2) extrapolation of the profile in the PBL where no in situ measurements were collected and (3) profile assumptions above the upper limit of flight observations. Errors due to uncertainty in tropopause height were analyzed by shifting the tropopause height upward by 1 km, and the results are also listed in Table 4. These results indicated that the extrapolation method and assumptions used to construct profiles where no measurements were made were the primary source of errors, among which the greatest error was from the profile above the upper limit of the flight (0.323 ppm). Errors due to uncertainty in the tropopause height were also non-negligible. Because of the lack of observation data near the surface, the missing mea-surements were directly replaced by the data at the lowest altitude measured by the aircraft. The error caused by this practice is shown in Table 4, with an average of 0.079 ppm for X CO 2 . This is also the impact of the lack of near-surface observations on X CO 2 estimation. Therefore, observations from near the surface to about 1 km from other methods, such as in situ GHG measurements by tethered balloon and high tower, are necessary for accurate estimation of X CO 2 .

Conclusion
The vertical distributions of CO 2 and CH 4 were measured by using a Beechcraft King Air 350ER over Jiansanjiang, an extensive paddy area in northeast China, and three vertical profiles from 0.6 to 7.5 km were obtained on 7, 9 and 10 August. Measurements of the mole fraction of CO 2 showed an increase with height, whereas CH 4 decreased with height. These results are reasonable, because paddies are sinks for CO 2 and sources of CH 4 during the summer growing season. Comparing the observed profiles from aircraft with those from the carbon cycle data assimilation system Tan-Tracker (v1) and OCO-2 retrievals showed that the general vertical structure was consistent, but the values of mole fraction of CO 2 from Tan-Tracker and OCO-2 had negative bias estimates. The average bias between aircraft and OCO-2 is Figure 7. Extrapolated CO 2 profiles observed on 7, 9 and 10 August 2018, over Jiansanjiang by method (2). Red, blue and yellow solid lines show the aircraft-based (in situ) data collected on 7, 9 and 10 August, respectively, averaged for each flat stage of the flight. Dotted lines show the extrapolated parts of the profiles, with colors corresponding to sampling dates in accordance with the solid lines. Black horizontal lines show the tropopause height from NCEP reanalysis data. Table 4. Aircraft integration error budget of X CO 2 estimation for method (2). Errors in the three profiles from multiple error sources contributed to the calculation results of the integrated total column. There are four sources of error, similar to previously described error budgets : the contribution from the aircraft profile itself, the contribution from the unknown surface to the bottom of the profile, the contribution from the upper troposphere and stratosphere, and error from the tropopause height. Supplement. The supplement related to this article is available online at: https://doi.org/10.5194/amt-13-3595-2020-supplement.
Author contributions. MD and XS determined the main goal of this study. XS carried it out, analyzed the data and prepared the paper with contributions from all co-authors. YG and NC provided technical guidance for related instruments.