Time series are presented of radiocarbon and 13C contents in atmospheric carbon dioxide over eastern Europe (southern Poland), covering the periods 1983–1994 and 2000–2004. The carbon isotope composition was measured in biweekly composite samples of atmospheric CO2, collected about 20 m above the local ground level. The data for 2 observational sites are presented: i) city of Kraków (50°04′N, 19°55′E; 220 m asl; for 1983–1994 and 2000–2004); and ii) Kasprowy Wierch, Tatra Mountains (49°14′N, 19°56′E; 1989 m asl; for 2000–2004). The latter site is considered a regional reference station, relatively free of anthropogenic influences. During the period 1983–1994, observations in the Kraków area revealed a gradual decrease of 14C content with a broad minimum around 1991 and a small increase by about 10% in the subsequent years. δ13C also changes with time, showing a decreasing trend from approximately −9.6% in 1983, with a slope of −0.02%/yr. The observed trends for both isotopes coincide well with a substantial reduction of coal consumption in Poland and partial replacement of coal by natural gas, especially in urban regions. After 2000, the δ13C slightly increases, reaching a mean value of −10% in 2004, while Δ14C is below the reference level by ∼3.5%. Observations at Kasprowy Wierch (regional reference station) also reflect a diminishing input of fossil carbon into the regional atmosphere. The fossil component in atmospheric CO2, calculated with the aid of 14C data available for the 2 study periods, shows a reduction of anthropogenic input by a factor of 2, which is confirmed by annual statistics of coal consumption.

Regular measurements of atmospheric CO (2) mixing ratios and their carbon isotope composition ((13)C/(12)C and (14)C/(12)C ratios) performed between 2005 and 2009 at two sites of contrasting characteristics (Krakow and the remote mountain site Kasprowy Wierch) located in southern Poland were used to derive fossil fuel-related and biogenic contributions to the total CO (2) load measured at both sites. Carbon dioxide present in the atmosphere, not coming from fossil fuel and biogenic sources, was considered 'background' CO (2). In Krakow, the average contribution of fossil fuel CO (2) was approximately 3.4%. The biogenic component was of the same magnitude. Both components revealed a distinct seasonality, with the fossil fuel component reaching maximum values during winter months and the biogenic component shifted in phase by approximately 6 months. The partitioning of the local CO (2) budget for the Kasprowy Wierch site revealed large differences in the derived components: the fossil fuel component was approximately five times lower than that derived for Krakow, whereas the biogenic component was negative in summer, pointing to the importance of photosynthetic sink associated with extensive forests in the neighbourhood of the station. While the presented study has demonstrated the strength of combined measurements of CO (2) mixing ratios and their carbon isotope signature as efficient tools for elucidating the partitioning of local atmospheric CO (2) loads, it also showed the important role of the land cover and the presence of the soil in the footprint of the measurement location, which control the net biogenic surface CO (2) fluxes.

During the summer of 2018, a widespread drought developed over Northern and Central Europe. The increase in temperature and the reduction of soil moisture have influenced carbon dioxide (CO 2 ) exchange between the atmosphere and terrestrial ecosystems in various ways, such as a reduction of photosynthesis, changes in ecosystem respiration, or allowing more frequent fires. In this study, we characterize the resulting perturbation of the atmospheric CO 2 seasonal cycles. 2018 has a good coverage of European regions affected by drought, allowing the investigation of how ecosystem flux anomalies impacted spatial CO 2 gradients between stations. This density of stations is unprecedented compared to previous drought events in 2003 and 2015, particularly thanks to the deployment of the Integrated Carbon Observation System (ICOS) network of atmospheric greenhouse gas monitoring stations in recent years. Seasonal CO 2 cycles from 48 European stations were available for 2017 and 2018. Earlier data were retrieved for comparison from international databases or national networks. Here, we show that the usual summer minimum in CO 2 due to the surface carbon uptake was reduced by 1.4 ppm in 2018 for the 10 stations located in the area most affected by the temperature anomaly, mostly in Northern Europe. Notwithstanding, the CO 2 transition phases before and after July were slower in 2018 compared to 2017, suggesting an extension of the growing season, with either continued CO 2 uptake by photosynthesis and/or a reduction in respiration driven by the depletion of substrate for respiration inherited from the previous months due to the drought. For stations with sufficiently long time series, the CO 2 anomaly observed in 2018 was compared to previous European droughts in 2003 and 2015. Considering the areas most affected by the temperature anomalies, we found a higher CO 2 anomaly in 2003 (+3 ppm averaged over 4 sites), and a smaller anomaly in 2015 (+1 ppm averaged over 11 sites) compared to 2018. This article is part of the theme issue ‘Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Abstract. Concentration of radon (222Rn) in the near-ground atmosphere has been measured quasi-continuously from January 2005 to December 2009 at two continental sites in Europe: Heidelberg (south-west Germany) and Krakow (southern Poland). The atmosphere was sampled at ca. 30 and 20 m above the local ground. Both stations were equipped with identical instruments. Regular observations of 222Rn were supplemented by measurements of surface fluxes of this gas in the Krakow urban area, using two different approaches. The measured concentrations of 222Rn varied at both sites in a wide range, from less than 2.0 Bq m−3 to approximately 40 Bq m−3 in Krakow and 35 Bq m−3 in Heidelberg. The mean 222Rn content in Krakow, when averaged over the entire observation period, was 30% higher than in Heidelberg (5.86 ± 0.09 and 4.50 ± 0.07 Bq m−3, respectively). Distinct seasonality of 222Rn signal is visible in the obtained time series of 222Rn concentration, with higher values recorded generally during late summer and autumn. The surface 222Rn fluxes measured in Krakow also revealed a distinct seasonality, with broad maximum observed during summer and early autumn and minimum during the winter. When averaged over a 5-year observation period, the night-time surface 222Rn flux was equal to 46.8 ± 2.4 Bq m−2 h−1. Although the atmospheric 222Rn levels at Heidelberg and Krakow appeared to be controlled primarily by local factors, it was possible to evaluate the "continental effect" in atmospheric 222Rn content between both sites, related to gradual build-up of 222Rn concentration in the air masses travelling between Heidelberg and Krakow. The mean value of this build-up was equal to 0.78 ± 0.12 Bq m−3. The measured minimum 222Rn concentrations at both sites and the difference between them was interpreted in the framework of a simple box model coupled with HYSPLIT (Hybrid Single Particle Lagrangian Integrated Trajectory) analysis of air mass trajectories. The best fit of experimental data was obtained for the mean 222Rn flux over the European continent equal to 52 Bq m−2 h−1, the mean transport velocity of the air masses within the convective mixed layer of the planetary boundary layer (PBL) on their route from the Atlantic coast to Heidelberg and Krakow equal to 3.5 m s−1, the mean rate constant of 222Rn removal across the top of the PBL equal to the 222Rn decay constant and the mean height of the convective mixed layer equal to 1600 m.

Anthropogenic emissions of carbon dioxide (CO 2) and methane (CH 4) in the atmosphere constitute an important component of the related carbon budget. The main source of anthropogenic CO 2 is burning of fossil fuels, especially in densely populated areas. Similar emissions of CH 4 are associated with the agricultural sector, coal mining, and other human activities, such as waste management and storage and natural gas networks supplying methane to large urban, industrial centers. We discuss several methods aimed at characterizing and quantifying atmospheric loads and fluxes of CO 2 and CH 4 in Krakow, the second largest city in Poland. The methods are based on atmospheric observations of mixing ratios as well as isotopic composition of the investigated gases. Atmospheric mixing ratios of CO 2 and CH 4 were measured using gas chromatography (GC) and cavity ring-down spectroscopy (CRDS). The isotopic composition of CO 2 and CH 4 was analyzed using isotope ratio mass spectrometry (IRMS), accelerator mass spectrometry (AMS), and CRDS techniques. These data, combined with auxiliary information characterizing the intensity of vertical mixing in the lower atmosphere (height of the nocturnal boundary layer [NBL] and atmospheric 222 Rn concentration), were further used to quantify emission rates of CO 2 and CH 4 in the urban atmosphere of Krakow. These methods provide an efficient way of quantifying surface emissions of major greenhouse gases originating from distributed sources, thus complementing the widely used bottom-up methodology based on emission statistics.

Abstract. We developed a high-resolution surface flux inversion system based on the global Eulerian–Lagrangian coupled tracer transport model composed of the National Institute for Environmental Studies (NIES) transport model (TM; collectively NIES-TM) and the FLEXible PARTicle dispersion model (FLEXPART). The inversion system is named NTFVAR (NIES-TM–FLEXPART-variational) as it applies a variational optimization to estimate surface fluxes. We tested the system by estimating optimized corrections to natural surface CO2 fluxes to achieve the best fit to atmospheric CO2 data collected by the global in situ network as a necessary step towards the capability of estimating anthropogenic CO2 emissions. We employed the Lagrangian particle dispersion model (LPDM) FLEXPART to calculate surface flux footprints of CO2 observations at a spatial resolution of 0.1∘×0.1∘. The LPDM is coupled with a global atmospheric tracer transport model (NIES-TM). Our inversion technique uses an adjoint of the coupled transport model in an iterative optimization procedure. The flux error covariance operator was implemented via implicit diffusion. Biweekly flux corrections to prior flux fields were estimated for the years 2010–2012 from in situ CO2 data included in the Observation Package (ObsPack) data set. High-resolution prior flux fields were prepared using the Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) for fossil fuel combustion, the Global Fire Assimilation System (GFAS) for biomass burning, the Vegetation Integrative SImulator for Trace gases (VISIT) model for terrestrial biosphere exchange, and the Ocean Tracer Transport Model (OTTM) for oceanic exchange. The terrestrial biospheric flux field was constructed using a vegetation mosaic map and a separate simulation of CO2 fluxes at a daily time step by the VISIT model for each vegetation type. The prior flux uncertainty for the terrestrial biosphere was scaled proportionally to the monthly mean gross primary production (GPP) by the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD17 product. The inverse system calculates flux corrections to the prior fluxes in the form of a relatively smooth field multiplied by high-resolution patterns of the prior flux uncertainties for land and ocean, following the coastlines and fine-scale vegetation productivity gradients. The resulting flux estimates improved the fit to the observations taken at continuous observation sites, reproducing both the seasonal and short-term concentration variabilities including high CO2 concentration events associated with anthropogenic emissions. The use of a high-resolution atmospheric transport in global CO2 flux inversions has the advantage of better resolving the transported mixed signals from the anthropogenic and biospheric sources in densely populated continental regions. Thus, it has the potential to achieve better separation between fluxes from terrestrial ecosystems and strong localized sources, such as anthropogenic emissions and forest fires. Further improvements in the modelling system are needed as our posterior fit was better than that of the National Oceanic and Atmospheric Administration (NOAA)'s CarbonTracker for only a fraction of the monitoring sites, i.e. mostly at coastal and island locations where background and local flux signals are mixed.

Abstract. The results of regular observations of atmospheric CO2 mixing ratios and its carbon isotope composition (δ13C, Δ14C), carried out at two continental sites located in central Europe are presented and discussed. The sites (Kasprowy Wierch, 49°14' N, 19°59' E, 1989 m a.s.l.; Krakow, 50°04' N, 19°55' E, 220 m a.s.l.), are located in two contrasting environments: (i) high-altitude mountaneous area, relatively free of anthropogenic influences, and (ii) typical urban environment with numerous local sources of carbon dioxide. Despite of relative proximity of those sites (ca. 100 km), substantial differences in both the recorded CO2 levels and their isotopic composition were detected. The CO2 mixing ratios measured in the urban atmosphere revealed quasi-permanent excess concentration of this gas when compared with near-by background atmosphere. The annual mean CO2 concentration recorded in Krakow in 2004 was almost 10% higher than that recorded at high-altitude mountain site (Kasprowy Wierch). Such effect is occuring probably in all urban centers. Carbon isotopic composition of atmospheric CO2 proved to be efficient tool for identification the surface CO2 fluxes into the atmosphere related to fossil fuel burning and their influence on the recorded levels of this gas in the local atmosphere. The available records of Δ14C for Krakow and Kasprowy Wierch suggest gradual reduction of 14C-free CO2 fluxes into the urban atmosphere of Krakow in the past several years.

Abstract. We developed a high-resolution surface flux inversion system based on the global Lagrangian–Eulerian coupled tracer transport model composed of National Institute for Environmental Studies Transport Model (NIES-TM) and FLEXible PARTicle dispersion model (FLEXPART). The inversion system is named NTFVAR (NIES-TM-FLEXPART-variational) as it applies variational optimisation to estimate surface fluxes. We tested the system by estimating optimized corrections to natural surface CO2 fluxes to achieve best fit to atmospheric CO2 data collected by the global in-situ network, as a necessary step towards capability of estimating anthropogenic CO2 emissions. We employ the Lagrangian particle dispersion model (LPDM) FLEXPART to calculate the surface flux footprints of CO2 observations at a 0.1° × 0.1° spatial resolution. The LPDM is coupled to a global atmospheric tracer transport model (NIES-TM). Our inversion technique uses an adjoint of the coupled transport model in an iterative optimization procedure. The flux error covariance operator is being implemented via implicit diffusion. Biweekly flux corrections to prior flux fields were estimated for the years 2010–2012 from in-situ CO2 data included in the Observation Package (ObsPack) dataset. High-resolution prior flux fields were prepared using Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) for fossil fuel combustion, Global Fire Assimilation System (GFAS) for biomass burning, the Vegetation Integrative SImulator for Trace gases (VISIT) model for terrestrial biosphere exchange and Ocean Tracer Transport Model (OTTM) for oceanic exchange. The terrestrial biospheric flux field was constructed using a vegetation mosaic map and separate simulation of CO2 fluxes at daily time step by the VISIT model for each vegetation type. The prior flux uncertainty for terrestrial biosphere was scaled proportionally to the monthly mean Gross Primary Production (GPP) by the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD17 product. The inverse system calculates flux corrections to the prior fluxes in the form of a relatively smooth field multiplied by high-resolution patterns of the prior flux uncertainties for land and ocean, following the coastlines and vegetation productivity gradients. The resulting flux estimates improve fit to the observations at continuous observations sites, reproducing both the seasonal variation and short-term concentration variability, including high CO2 concentration events associated with anthropogenic emissions. The use of high-resolution atmospheric transport in global CO2 flux inversion has the advantage of better resolving the transport from the mix of the anthropogenic and biospheric sources in densely populated continental regions and shows potential for better separation between fluxes from terrestrial ecosystems and strong localised sources such as anthropogenic emissions and forest fires. Further improvements in the modelling system are needed as the posterior fit is better than that by the National Oceanic and Atmospheric Administration (NOAA) CarbonTracker only for a fraction of the monitoring sites, mostly at coastal and island locations experiencing mix of background and local flux signals.