Search

Ammonia (NH3), a significant precursor of particulate matter, affects not only biodiversity, ecosystems, and soil acidification but also climate and human health. In addition, its concentrations are constantly rising due to increasing feeding needs and the large use of fertilization and animal farming. Despite the significance of ammonia, its emissions are associated with large uncertainties, while its atmospheric abundance is difficult to measure. Nowadays, satellite products can effectively measure ammonia with low uncertainty and a global coverage. Here, we use satellite observations of column ammonia in combination with an inversion algorithm to derive ammonia emissions with a high resolution over Europe for the period 2013–2020. Ammonia emissions peak in northern Europe due to agricultural application and livestock management, in western Europe (industrial activity), and over Spain (pig farming). Emissions have decreased by −26 % since 2013 (from 5431 Gg in 2013 to 3994 Gg in 2020), showing that the abatement strategies adopted by the European Union have been very efficient. The slight increase (+4.4 %) in 2015 is also reproduced here and is attributed to some European countries exceeding annual emission targets. Ammonia emissions are low in winter (286 Gg) and peak in summer (563 Gg) and are dominated by the temperature-dependent volatilization of ammonia from the soil. The largest emission decreases were observed in central and eastern Europe (−38 %) and in western Europe (−37 %), while smaller decreases were recorded in northern (−17 %) and southern Europe (−7.6 %). When complemented with ground observations, modelled concentrations using the posterior emissions showed improved statistics, also following the observed seasonal trends. The posterior emissions presented here also agree well with respective estimates reported in the literature and inferred from bottom-up and top-down methodologies. These results indicate that satellite measurements combined with inverse algorithms constitute a robust tool for emission estimates and can infer the evolution of ammonia emissions over large timescales.

Low concentrations of 106Ru were detected across Europe at the turn of September and October 2017. The origin of 106Ru has still not been confirmed; however, current studies agree that the release occurred probably near Mayak in the southern Urals. The source reconstructions are mostly based on an analysis of concentration measurements coupled with an atmospheric transport model. Since reasonable temporal resolution of concentration measurements is crucial for proper source term reconstruction, the standard 1-week sampling interval could be limiting. In this paper, we present an investigation of the usability of the newly developed AMARA (Autonomous Monitor of Atmospheric Radioactive Aerosol) and CEGAM (carousel gamma spectrometry) real-time monitoring systems, which are based on the gamma-ray counting of aerosol filters and allow for determining the moment when 106Ru arrived at the monitoring site within approx. 1 h and detecting activity concentrations as low as several mBq m−3 in 4 h intervals. These high-resolution data were used for inverse modeling of the 106Ru release. We perform backward runs of the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) atmospheric transport model driven with meteorological data from the Global Forecast System (GFS), and we construct a source–receptor sensitivity (SRS) matrix for each grid cell of our domain. Then, we use our least squares with adaptive prior covariance (LS-APC) method to estimate possible locations of the release and the source term of the release. With Czech monitoring data, the use of concentration measurements from the standard regime and from the real-time regime is compared, and a better source reconstruction for the real-time data is demonstrated in the sense of the location of the source and also the temporal resolution of the source. The estimated release location, Mayak, and the total estimated source term, 237±107 TBq, are in agreement with previous studies. Finally, the results based on the Czech monitoring data are validated with the IAEA-reported (International Atomic Energy Agency) dataset with a much better spatial resolution, and the agreement between the IAEA dataset and our reconstruction is demonstrated. In addition, we validated our findings also using the FLEXPART (FLEXible PARTicle dispersion) model coupled with meteorological analyses from the European Centre for Medium-Range Weather Forecasts (ECMWF).

Estimation of the temporal profile of an atmospheric release, also called the source term, is an important problem in environmental sciences. The problem can be formalized as a linear inverse problem wherein the unknown source term is optimized to minimize the difference between the measurements and the corresponding model predictions. The problem is typically ill-posed due to low sensor coverage of a release and due to uncertainties, e.g., in measurements or atmospheric transport modeling; hence, all state-of-the-art methods are based on some form of regularization of the problem using additional information. We consider two kinds of additional information: the prior source term, also known as the first guess, and regularization parameters for the shape of the source term. While the first guess is based on information independent of the measurements, such as the physics of the potential release or previous estimations, the regularization parameters are often selected by the designers of the optimization procedure. In this paper, we provide a sensitivity study of two inverse methodologies on the choice of the prior source term and regularization parameters of the methods. The sensitivity is studied in two cases: data from the European Tracer Experiment (ETEX) using FLEXPART v8.1 and the caesium-134 and caesium-137 dataset from the Chernobyl accident using FLEXPART v10.3.

In the fall of 2011, iodine-131 (131I) was detected at several radionuclide monitoring stations in central Europe. After investigation, the International Atomic Energy Agency (IAEA) was informed by Hungarian authorities that 131I was released from the Institute of Isotopes Ltd. in Budapest, Hungary. It was reported that a total activity of 342 GBq of 131I was emitted between 8 September and 16 November 2011. In this study, we use the ambient concentration measurements of 131I to determine the location of the release as well as its magnitude and temporal variation. As the location of the release and an estimate of the source strength became eventually known, this accident represents a realistic test case for inversion models. For our source reconstruction, we use no prior knowledge. Instead, we estimate the source location and emission variation using only the available 131I measurements. Subsequently, we use the partial information about the source term available from the Hungarian authorities for validation of our results. For the source determination, we first perform backward runs of atmospheric transport models and obtain source-receptor sensitivity (SRS) matrices for each grid cell of our study domain. We use two dispersion models, FLEXPART and Hysplit, driven with meteorological analysis data from the global forecast system (GFS) and from European Centre for Medium-range Weather Forecasts (ECMWF) weather forecast models. Second, we use a recently developed inverse method, least-squares with adaptive prior covariance (LS-APC), to determine the 131I emissions and their temporal variation from the measurements and computed SRS matrices. For each grid cell of our simulation domain, we evaluate the probability that the release was generated in that cell using Bayesian model selection. The model selection procedure also provides information about the most suitable dispersion model for the source term reconstruction. Third, we select the most probable location of the release with its associated source term and perform a forward model simulation to study the consequences of the iodine release. Results of these procedures are compared with the known release location and reported information about its time variation. We find that our algorithm could successfully locate the actual release site. The estimated release period is also in agreement with the values reported by IAEA and the reported total released activity of 342 GBq is within the 99 % confidence interval of the posterior distribution of our most likely model.

Estimation of pollutant releases into the atmosphere is an important problem in the environmental sciences. It is typically formalized as an inverse problem using a linear model that can explain observable quantities (e.g., concentrations or deposition values) as a product of the source-receptor sensitivity (SRS) matrix obtained from an atmospheric transport model multiplied by the unknown source-term vector. Since this problem is typically ill-posed, current state-of-the-art methods are based on regularization of the problem and solution of a formulated optimization problem. This procedure depends on manual settings of uncertainties that are often very poorly quantified, effectively making them tuning parameters. We formulate a probabilistic model, that has the same maximum likelihood solution as the conventional method using pre-specified uncertainties. Replacement of the maximum likelihood solution by full Bayesian estimation also allows estimation of all tuning parameters from the measurements. The estimation procedure is based on the variational Bayes approximation which is evaluated by an iterative algorithm. The resulting method is thus very similar to the conventional approach, but with the possibility to also estimate all tuning parameters from the observations. The proposed algorithm is tested and compared with the standard methods on data from the European Tracer Experiment (ETEX) where advantages of the new method are demonstrated. A MATLAB implementation of the proposed algorithm is available for download.

The source term of atmospheric emissions of hazardous materials is a crucial aspect of the analysis of unintended release. Motivated by wildfires of regions contaminated by radioactivity, the focus is placed on the case of airborne transmission of material from 5 dimensions: spatial location described by longitude and latitude in a given area with potentially many sources, time profiles, height above ground level, and the size of particles carrying the material. Since the atmospheric inverse problem is typically ill-posed and the number of measurements is usually too low to estimate the whole 5D tensor, some prior information is necessary. For the first time in this domain, a method based on deep image prior utilizing the structure of a deep neural network to regularize the inversion is proposed. The network is initialized randomly without the need to train it on any dataset first. In tandem with variational optimization, this approach not only introduces smoothness in the spatial estimate of the emissions but also reduces the number of unknowns by enforcing a prior covariance structure in the source term. The strengths of this method are demonstrated on the case of 137 Cs emissions during the Chernobyl wildfires in 2020., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2025 Elsevier B.V. All rights reserved.)

A swarm of semi-classical quantum mechanics/molecular mechanics molecular-dynamics simulations where OM2/MNDO is combined with the Gromacs program for consideration of explicit water is performed, solving the time-dependent Schrödinger equation in each step of the trajectories together with the Tully's fewest switches algorithm. Within this stochastic treatment, time dependent probabilities of the three lowest electronic states are determined. The fact that nucleobases are quickly deactivated is confirmed in the cytosine case where our best lifetime estimation is τ 1 =0.82 ps for the model with 100 water molecules with the SPCE force field and a time step of 0.1 fs. Lifetimes of the remaining molecules are visibly longer: 5-azacytosine, 2,4-diamino-1,3,5-triazine (DT), and 2,4,6-triamino-1,3,5-triazine (TT) molecules have an S 1 → S 0 de-excitation time of slightly above 10 ps. The lifetimes of the triazine family increases with the increasing number of exocyclic amino groups, that is, s -triazine < 2-amino-1,3,5-triazine < DT < TT. This can be explained by a higher mobility of the carbon-bonded hydrogen atoms in comparison with heavier amino groups since their movement is slowed down due to a substantially higher mass than hydrogen atoms, which can easier reach the out-of-plane positions required in the conical intersection structures. Moreover, bulkier NH 2 ligands suffer due to greater friction caused by the surrounding water environment. These mechanical aspects caused a change in the explored lifetime dependences in comparison with our previous gas-phase study.

We combine observations from Western USA and inverse modelling to constrain global atmospheric emissions of microplastics (MPs) and microfibers (MFs). The latter are used further to model their global atmospheric dynamics. Global annual MP emissions were calculated as 9.6 ± 3.6 Tg and MF emissions as 6.5 ± 2.9 Tg. Global average monthly MP concentrations were 47 ng m -3 and 33 ng m -3 for MFs, at maximum. The largest deposition of agricultural MPs occurred close to the world's largest agricultural regions. Road MPs mostly deposited in the East Coast of USA, Central Europe, and Southeastern Asia; MPs resuspended with mineral dust near Sahara and Middle East. Only 1.8% of the emitted mass of oceanic MPs was transferred to land, and 1.4% of land MPs to ocean; the rest were deposited in the same environment. Previous studies reported that 0.74-1.9 Tg y -1 of land-based atmospheric MPs/MFs (< 5 mm) are transported to the ocean, while riverine transport is between 3.3 and 14 Tg y -1 . We calculate that 0.418 ± 0.201 Tg y -1 MPs/MFs (size up to 250 and 2500 µm) were transported from the land to ocean (large particles were ignored). Model validation against observations showed that particle removal must be urgently updated in global models., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Estimation of a source term, i.e. release rate, of atmospheric radionuclide emissions is of key interest for nuclear emergency response and further accident analysis. The source term estimate is, however, often very inaccurate due to biases in atmospheric transport and used meteorological analysis. We propose a method for atmospheric plume bias correction which uses not only concentrations modeled at a measuring site but also the information on concentration gradient from the neighborhood of each measuring site, i.e. information already available from the atmospheric transport model. To properly regularize the model, we propose an elastic model of the plume bias correction based on regularization with the use of known topology of the measurement network. The proposed plume bias correction method can be coupled with an arbitrary source term estimation algorithm and can be instantly applied to any other atmospheric release of hazardous material. We demonstrate the method in two real cases. First, we use data from the European Tracer Experiment to validate the methodology. Second, we use data from the 106 Ru occurrence over Europe in 2017 to demonstrate the methodology in a more demanding case where agreement with state-of-the-art estimates is shown with much better reconstruction of measurements., (Copyright © 2021 Elsevier B.V. All rights reserved.)