Abstract. We present airborne measurements of carbon dioxide (CO 2 ), carbon monoxide (CO), ozone (O 3 ), equivalent black carbon (EBC) and ultra fine particles over NorthEastern Siberia in July 2008 performed during the YAK-AEROSIB/POLARCAT experiment. During a "golden day" (11 July 2008) a number of biomass burning plumes were encountered with CO mixing ratio enhancements of up to 500 ppb relative to a background of 90 ppb. Number concentrations of aerosols in the size range 3.5-200 nm peaked at 4000 cm −3 and the EBC content reached 1.4 µg m −3 . These high concentrations were caused by forest fires in the vicinity of the landing airport in Yakutsk where measurements in fresh smoke could be made during the descent. We estimate a combustion efficiency of 90 ± 3% based on CO and CO 2 measurements and a CO emission factor of 65.5±10.8 g CO per kilogram of dry matter burned. This suggests a potential increase in the average northern hemispheric CO mixing ratio of 3.0-7.2 ppb per million hectares of Siberian forest burned. For BC, we estimate an emission factor of 0.52 ± 0.07 g BC kg −1 , comparable to values reported in the literature. The emission ratio of ultra-fine particles (3.5-200 nm) was 26 cm −3 (ppb CO) −1 , consistent with other airborne studies.The transport of identified biomass burning plumes was investigated using the FLEXPART Lagrangian model. Based on sampling of wildfire plumes from the same source but with different atmospheric ages derived from FLEXPART, Correspondence to: J.-D. Paris (jean-daniel.paris@lsce.ipsl.fr) we estimate that the e-folding lifetimes of EBC and ultra fine particles (between 3.5 and 200 nm in size) against removal and growth processes are 5.1 and 5.5 days respectively, supporting lifetime estimates used in various modelling studies.

Two airborne campaigns were carried out to measure the tropospheric concentrations and variability of CO2, CO and O3 over Siberia. In order to quantify the influence of remote and regional natural and anthropogenic sources, we analysed a total of 52 vertical profiles of these species collected in April and September 2006, every ∼200 km and up to 7 km altitude. CO2 and CO concentrations were high in April 2006 (respectively 385–390 ppm CO2 and 160–200 ppb CO) compared to background values. CO concentrations up to 220 ppb were recorded above 3.5 km over eastern Siberia, with enhancements in 500–1000 m thick layers. The presence of CO enriched air masses resulted from a quick frontal uplift of a polluted air mass exposed to northern China anthropogenic emissions and to fire emissions in northern Mongolia. A dominant Asian origin for CO above 4 km (71.0%) contrasted with a dominant European origin below this altitude (70.9%) was deduced both from a transport model analysis, and from the contrasted ΔCO/ΔCO2 ratio vertical distribution. In September 2006, a significant O3 depletion (∼–30 ppb) was repeatedly observed in the boundary layer, as diagnosed from virtual potential temperature profiles and CO2 gradients, compared to the free troposphere aloft, suggestive of a strong O3 deposition over Siberian forests.

Abstract. Eight surface observation sites providing quasicontinuous measurements of atmospheric methane mixing ratios have been operated since the mid-2000's in Siberia. For the first time in a single work, we assimilate 1 year of these in situ observations in an atmospheric inversion. Our objective is to quantify methane surface fluxes from anthropogenic and wetland sources at the mesoscale in the Siberian lowlands for the year 2010. To do so, we first inquire about the way the inversion uses the observations and the way the fluxes are constrained by the observation sites. As atmospheric inversions at the mesoscale suffer from mis-quantified sources of uncertainties, we follow recent innovations in inversion techniques and use a new inversion approach which quantifies the uncertainties more objectively than the previous inversion systems. We find that, due to errors in the representation of the atmospheric transport and redundant pieces of information, only one observation every few days is found valuable by the inversion. The remaining high-resolution quasicontinuous signal is representative of very local emission patterns difficult to analyse with a mesoscale system. An analysis of the use of information by the inversion also reveals that the observation sites constrain methane emissions within a radius of 500 km. More observation sites than the ones currently in operation are then necessary to constrain the whole Siberian lowlands. Still, the fluxes within the constrained areas are quantified with objectified uncertainties. Finally, the tolerance intervals for posterior methane fluxes are of roughly 20 % (resp. 50 %) of the fluxes for anthropogenic (resp. wetland) sources. About 50-70 % of Siberian lowlands emissions are constrained by the inversion on average on an annual basis. Extrapolating the figures on the constrained areas to the whole Siberian lowlands, we find a regional methane budget of 5-28 TgCH 4 for the year 2010, i.e. 1-5 % of the global methane emissions. As very few in situ observations are available in the region of interest, observations of methane total columns from the Greenhouse Gas Observing SATellite (GOSAT) are tentatively used for the evaluation of the inversion results, but they exhibit only a marginal signal from the fluxes within the region of interest.

Abstract. Black carbon (BC) emissions from open biomass burning (BB) are known to have a considerable impact on the radiative budget of the atmosphere at both global and regional scales; however, these emissions are poorly constrained in models by atmospheric observations, especially in remote regions. Here, we investigate the feasibility of constraining BC emissions from BB using satellite observations of the aerosol absorption optical depth (AAOD) and the aerosol extinction optical depth (AOD) retrieved from OMI (Ozone Monitoring Instrument) and MODIS (Moderate Resolution Imaging Spectroradiometer) measurements, respectively. We consider the case of Siberian BB BC emissions, which have the strong potential to impact the Arctic climate system. Using aerosol remote sensing data collected at Siberian sites of the AErosol RObotic NETwork (AERONET) along with the results of the fourth Fire Lab at Missoula Experiment (FLAME-4), we establish an empirical parameterization relating the ratio of the elemental carbon (EC) and organic carbon (OC) contents in BB aerosol to the ratio of AAOD and AOD at the wavelengths of the satellite observations. Applying this parameterization to the BC and OC column amounts simulated using the CHIMERE chemistry transport model, we optimize the parameters of the BB emission model based on MODIS measurements of the fire radiative power (FRP); we then obtain top-down optimized estimates of the total monthly BB BC amounts emitted from intense Siberian fires that occurred from May to September 2012. The top-down estimates are compared to the corresponding values obtained using the Global Fire Emissions Database (GFED4) and the Fire Emission Inventory–northern Eurasia (FEI-NE). Our simulations using the optimized BB aerosol emissions are verified against AAOD and AOD data that were withheld from the estimation procedure. The simulations are further evaluated against in situ EC and OC measurements at the Zotino Tall Tower Observatory (ZOTTO) and also against aircraft aerosol measurement data collected in the framework of the Airborne Extensive Regional Observations in SIBeria (YAK-AEROSIB) experiments. We conclude that our BC and OC emission estimates, considered with their confidence intervals, are consistent with the ensemble of the measurement data analyzed in this study. Siberian fires are found to emit 0.41±0.14 Tg of BC over the whole 5-month period considered; this estimate is a factor of 2 larger and a factor of 1.5 smaller than the corresponding estimates based on the GFED4 (0.20 Tg) and FEI-NE (0.61 Tg) data, respectively. Our estimates of monthly BC emissions are also found to be larger than the BC amounts calculated using the GFED4 data and smaller than those calculated using the FEI-NE data for any of the 5 months. Particularly large positive differences of our monthly BC emission estimates with respect to the GFED4 data are found in May and September. This finding indicates that the GFED4 database is likely to strongly underestimate BC emissions from agricultural burns and grass fires in Siberia. All of these differences have important implications for climate change in the Arctic, as it is found that about a quarter of the huge BB BC mass emitted in Siberia during the fire season of 2012 was transported across the polar circle into the Arctic. Overall, the results of our analysis indicate that a combination of the available satellite observations of AAOD and AOD can provide the necessary constraints on BB BC emissions.

Abstract. A four-dimensional variational (4D-Var) method is a popular algorithm for inverting atmospheric greenhouse gas (GHG) measurements. In order to meet the computationally intense 4D-Var iterative calculation, offline forward and adjoint transport models are developed based on the Nonhydrostatic ICosahedral Atmospheric Model (NICAM). By introducing flexibility into the temporal resolution of the input meteorological data, the forward model developed in this study is not only computationally efficient, it is also found to nearly match the transport performance of the online model. In a transport simulation of atmospheric carbon dioxide (CO 2 ), the data-thinning error (error resulting from reduction in the time resolution of the meteorological data used to drive the offline transport model) is minimized by employing high temporal resolution data of the vertical diffusion coefficient; with a low 6-hourly temporal resolution, significant concentration biases near the surface are introduced. The new adjoint model can be run in discrete or continuous adjoint mode for the advection process. The discrete adjoint is characterized by perfect adjoint relationship with the forward model that switches off the flux limiter, while the continuous adjoint is characterized by an imperfect but reasonable adjoint relationship with its corresponding forward model. In the latter case, both the forward and adjoint models use the flux limiter to ensure the monotonicity of tracer concentrations and sensitivities. Trajectory analysis for high CO 2 concentration events are performed to test adjoint sensitivities. We also demonstrate the potential usefulness of our adjoint model for diagnosing tracer transport. Both the offline forward and adjoint models have computational efficiency about 10 times higher than the online model. A description of our new 4D-Var system that includes an optimization method, along with its application in an atmospheric CO 2 inversion and the effects of using either the discrete or continuous adjoint method, is presented in an accompanying paper (Niwa et al., 2016).