Trace and heavy metal and pesticide contents of different tobacco brands sold in Finland during 1920 to 1984 were determined by atomic absorption spectroscopy and gas chromatography-mass spectrometry with the purpose of investigating the variation in contents over time. In addition, the following were studied: the effect of commercial filters in cigarettes in preventing the inhalation of metals, the transfer of metals and pesticides to the mainstream smoke; and the correlation between the cadmium content of adipose tissue of smokers and their smoking habits. No differences were found in trace and heavy metal contents of different brands and packs from different decades. The pesticide residues in tobacco have followed the pattern of their use in the tobacco plantation. The dichlorodiphenyltrichloroethane (DDT) content of tobacco has decreased during the past 20 yr to about 1/200 of the peak value, i.e., from 34.5 to 0.17 micrograms/g. Cigarette filters significantly prevent the inhalation of cadmium, lead, magnesium, and iron. The mean content of cadmium in fat tissue of male smokers was four times that of non-smokers. The difference was statistically significant.

Abstract. Reliable quantification of the sources and sinks of greenhouse gases, together with trends and uncertainties, is essential to monitoring the progress in mitigating anthropogenic emissions under the Paris Agreement. This study provides a consolidated synthesis of CH4 and N2O emissions with consistently derived state-of-the-art bottom-up (BU) and top-down (TD) data sources for the European Union and UK (EU27 + UK). We integrate recent emission inventory data, ecosystem process-based model results and inverse modeling estimates over the period 1990–2017. BU and TD products are compared with European national greenhouse gas inventories (NGHGIs) reported to the UN climate convention UNFCCC secretariat in 2019. For uncertainties, we used for NGHGIs the standard deviation obtained by varying parameters of inventory calculations, reported by the member states (MSs) following the recommendations of the IPCC Guidelines. For atmospheric inversion models (TD) or other inventory datasets (BU), we defined uncertainties from the spread between different model estimates or model-specific uncertainties when reported. In comparing NGHGIs with other approaches, a key source of bias is the activities included, e.g., anthropogenic versus anthropogenic plus natural fluxes. In inversions, the separation between anthropogenic and natural emissions is sensitive to the geospatial prior distribution of emissions. Over the 2011–2015 period, which is the common denominator of data availability between all sources, the anthropogenic BU approaches are directly comparable, reporting mean emissions of 20.8 Tg CH4 yr−1 (EDGAR v5.0) and 19.0 Tg CH4 yr−1 (GAINS), consistent with the NGHGI estimates of 18.9 ± 1.7 Tg CH4 yr−1. The estimates of TD total inversions give higher emission estimates, as they also include natural emissions. Over the same period regional TD inversions with higher-resolution atmospheric transport models give a mean emission of 28.8 Tg CH4 yr−1. Coarser-resolution global TD inversions are consistent with regional TD inversions, for global inversions with GOSAT satellite data (23.3 Tg CH4 yr−1) and surface network (24.4 Tg CH4 yr−1). The magnitude of natural peatland emissions from the JSBACH–HIMMELI model, natural rivers and lakes emissions, and geological sources together account for the gap between NGHGIs and inversions and account for 5.2 Tg CH4 yr−1. For N2O emissions, over the 2011–2015 period, both BU approaches (EDGAR v5.0 and GAINS) give a mean value of anthropogenic emissions of 0.8 and 0.9 Tg N2O yr−1, respectively, agreeing with the NGHGI data (0.9 ± 0.6 Tg N2O yr−1). Over the same period, the average of the three total TD global and regional inversions was 1.3 ± 0.4 and 1.3 ± 0.1 Tg N2O yr−1, respectively. The TD and BU comparison method defined in this study can be operationalized for future yearly updates for the calculation of CH4 and N2O budgets both at the EU+UK scale and at the national scale. The referenced datasets related to figures are visualized at https://doi.org/10.5281/zenodo.4590875 (Petrescu et al., 2020b).

The surface albedo time series, CLARA-A2 SAL, was used to study trends in the snowmelt start and end dates, the melting season length and the albedo value preceding the melt onset in Finland during 1982-2016. In addition, the melt onset from the JSBACH land surface model was compared with the timing of green-up estimated from Moderate Resolution Imaging Spectroradiometer (MODIS) data. Moreover, the melt onset was compared with the timing of the greening up based on MODIS data. Similarly, the end of snowmelt timing predicted by JSBACH was compared with the meltoff dates based on the Finnish Meteorological Institute (FMI) operational in situ measurements and the Fractional Snow Cover (FSC) time-series product provided by the EU FP7 CryoLand project. It was found that the snowmelt date estimated using the 20 % threshold of the albedo range during the melting period corresponded well to the melt estimate of the permanent snow layer. The longest period, during which the ground is continuously half or more covered by snow, defines the permanent snow layer (Solantie et al., 1996). The greening up followed within 5-13 days the date when the albedo reached the 1 % threshold of the albedo dynamic range during the melting period. The time difference between greening up and complete snowmelt was smaller in mountainous areas than in coastal areas. In two northern vegetation map areas (Northern Karelia-Kainuu and Southwestern Lapland), a clear trend towards earlier snowmelt onset (5-6 days per decade) and increasing melting season length (6-7 days per decade) was observed. In the forested part of northern Finland, a clear decreasing trend in albedo (2 %-3 % per decade in absolute albedo percentage) before the start of the melt onset was observed. The decreasing albedo trend was found to be due to the increased stem volume.

Abstract. Knowledge of the spatial distribution of the fluxes of greenhouse gases (GHGs) and their temporal variability as well as flux attribution to natural and anthropogenic processes is essential to monitoring the progress in mitigating anthropogenic emissions under the Paris Agreement and to inform its global stocktake. This study provides a consolidated synthesis of CH4 and N2O emissions using bottom-up (BU) and top-down (TD) approaches for the European Union and UK (EU27 + UK) and updates earlier syntheses (Petrescu et al., 2020, 2021). The work integrates updated emission inventory data, process-based model results, data-driven sector model results and inverse modeling estimates, and it extends the previous period of 1990–2017 to 2019. BU and TD products are compared with European national greenhouse gas inventories (NGHGIs) reported by parties under the United Nations Framework Convention on Climate Change (UNFCCC) in 2021. Uncertainties in NGHGIs, as reported to the UNFCCC by the EU and its member states, are also included in the synthesis. Variations in estimates produced with other methods, such as atmospheric inversion models (TD) or spatially disaggregated inventory datasets (BU), arise from diverse sources including within-model uncertainty related to parameterization as well as structural differences between models. By comparing NGHGIs with other approaches, the activities included are a key source of bias between estimates, e.g., anthropogenic and natural fluxes, which in atmospheric inversions are sensitive to the prior geospatial distribution of emissions. For CH4 emissions, over the updated 2015–2019 period, which covers a sufficiently robust number of overlapping estimates, and most importantly the NGHGIs, the anthropogenic BU approaches are directly comparable, accounting for mean emissions of 20.5 Tg CH4 yr−1 (EDGARv6.0, last year 2018) and 18.4 Tg CH4 yr−1 (GAINS, last year 2015), close to the NGHGI estimates of 17.5±2.1 Tg CH4 yr−1. TD inversion estimates give higher emission estimates, as they also detect natural emissions. Over the same period, high-resolution regional TD inversions report a mean emission of 34 Tg CH4 yr−1. Coarser-resolution global-scale TD inversions result in emission estimates of 23 and 24 Tg CH4 yr−1 inferred from GOSAT and surface (SURF) network atmospheric measurements, respectively. The magnitude of natural peatland and mineral soil emissions from the JSBACH–HIMMELI model, natural rivers, lake and reservoir emissions, geological sources, and biomass burning together could account for the gap between NGHGI and inversions and account for 8 Tg CH4 yr−1. For N2O emissions, over the 2015–2019 period, both BU products (EDGARv6.0 and GAINS) report a mean value of anthropogenic emissions of 0.9 Tg N2O yr−1, close to the NGHGI data (0.8±55 % Tg N2O yr−1). Over the same period, the mean of TD global and regional inversions was 1.4 Tg N2O yr−1 (excluding TOMCAT, which reported no data). The TD and BU comparison method defined in this study can be operationalized for future annual updates for the calculation of CH4 and N2O budgets at the national and EU27 + UK scales. Future comparability will be enhanced with further steps involving analysis at finer temporal resolutions and estimation of emissions over intra-annual timescales, which is of great importance for CH4 and N2O, and may help identify sector contributions to divergence between prior and posterior estimates at the annual and/or inter-annual scale. Even if currently comparison between CH4 and N2O inversion estimates and NGHGIs is highly uncertain because of the large spread in the inversion results, TD inversions inferred from atmospheric observations represent the most independent data against which inventory totals can be compared. With anticipated improvements in atmospheric modeling and observations, as well as modeling of natural fluxes, TD inversions may arguably emerge as the most powerful tool for verifying emission inventories for CH4, N2O and other GHGs. The referenced datasets related to figures are visualized at https://doi.org/10.5281/zenodo.7553800 (Petrescu et al., 2023).

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level