Prioritisation of pesticides and target organ systems for dietary cumulative risk assessment based on the 2019-2021 monitoring cycle.

Aiming at accelerating the implementation of cumulative risk assessment to pesticide residues, this report describes a two-step prioritisation analysis, on individual pesticides and on target organ systems, that allows to identify (i) low-priority substances expected to have a marginal contribution to cumulative risk, and (ii) high priority organ systems to be addressed in future cumulative risk assessments. The analysis encompassed 350 substances and 36 raw primary commodities of plant origin surveyed in the monitoring cycle 2019-2021, carried out in 30 population groups, covering 3 age classes, and 17 EU countries. Probabilistic exposure calculations, for chronic and acute effects, were executed on the occurrence and consumption data by a two-dimensional procedure, modelling variability and uncertainty. In the first step, the prioritisation method adopted allowed to reduce the number of substances by about 80%. These substances were in turn grouped based on their capacity to cause toxicological effects on common organ systems and, as second step, probabilistic combined exposure calculations were carried out for 16 target organ systems. This step allowed to identify the organ systems that need further assessment, reducing their initial number by about 70%. The organ systems would need to be prioritised as follows: reproductive and developmental toxicity, liver, kidney, male reproductive system, and haematopoietic system and haematology. The sources of uncertainty stemming from the modelling procedure and from methodological assumptions were discussed and their impact qualitatively assessed. Overall, it was concluded that the risk estimates for the different organ systems were more likely to be overestimated than underestimated.

Near-real-time global gridded daily CO2 emissions 2021.

We present a near-real-time global gridded daily CO2 emissions dataset (GRACED) throughout 2021. GRACED provides gridded CO2 emissions at a 0.1° × 0.1° spatial resolution and 1-day temporal resolution from cement production and fossil fuel combustion over seven sectors, including industry, power, residential consumption, ground transportation, international aviation, domestic aviation, and international shipping. GRACED is prepared from the near-real-time daily national CO2 emissions estimates (Carbon Monitor), multi-source spatial activity data emissions and satellite NO2 data for time variations of those spatial activity data. GRACED provides the most timely overview of emissions distribution changes, which enables more accurate and timely identification of when and where fossil CO2 emissions have rebounded and decreased. Uncertainty analysis of GRACED gives a grid-level two-sigma uncertainty of value of ±19.9% in 2021, indicating the reliability of GRACED was not sacrificed for the sake of higher spatiotemporal resolution that GRACED provides. Continuing to update GRACED in a timely manner could help policymakers monitor energy and climate policies' effectiveness and make adjustments quickly.

Near-real-time global gridded daily CO2 emissions.

Precise and high-resolution carbon dioxide (CO2) emission data is of great importance in achieving carbon neutrality around the world. Here we present for the first time the near-real-time Global Gridded Daily CO2 Emissions Dataset (GRACED) from fossil fuel and cement production with a global spatial resolution of 0.1° by 0.1° and a temporal resolution of 1 day. Gridded fossil emissions are computed for different sectors based on the daily national CO2 emissions from near-real-time dataset (Carbon Monitor), the spatial patterns of point source emission dataset Global Energy Infrastructure Emissions Database (GID), Emission Database for Global Atmospheric Research (EDGAR), and spatiotemporal patters of satellite nitrogen dioxide (NO2) retrievals. Our study on the global CO2 emissions responds to the growing and urgent need for high-quality, fine-grained, near-real-time CO2 emissions estimates to support global emissions monitoring across various spatial scales. We show the spatial patterns of emission changes for power, industry, residential consumption, ground transportation, domestic and international aviation, and international shipping sectors from January 1, 2019, to December 31, 2020. This gives thorough insights into the relative contributions from each sector. Furthermore, it provides the most up-to-date and fine-grained overview of where and when fossil CO2 emissions have decreased and rebounded in response to emergencies (e.g., coronavirus disease 2019 [COVID-19]) and other disturbances of human activities of any previously published dataset. As the world recovers from the pandemic and decarbonizes its energy systems, regular updates of this dataset will enable policymakers to more closely monitor the effectiveness of climate and energy policies and quickly adapt.

High resolution temporal profiles in the Emissions Database for Global Atmospheric Research.

Emissions into the atmosphere from human activities show marked temporal variations, from inter-annual to hourly levels. The consolidated practice of calculating yearly emissions follows the same temporal allocation of the underlying annual statistics. However, yearly emissions might not reflect heavy pollution episodes, seasonal trends, or any time-dependant atmospheric process. This study develops high-time resolution profiles for air pollutants and greenhouse gases co- emitted by anthropogenic sources in support of atmospheric modelling, Earth observation communities and decision makers. The key novelties of the Emissions Database for Global Atmospheric Research (EDGAR) temporal profiles are the development of (i) country/region- and sector- specific yearly profiles for all sources, (ii) time dependent yearly profiles for sources with inter-annual variability of their seasonal pattern, (iii) country- specific weekly and daily profiles to represent hourly emissions, (iv) a flexible system to compute hourly emissions including input from different users. This work creates a harmonized emission temporal distribution to be applied to any emission database as input for atmospheric models, thus promoting homogeneity in inter-comparison exercises.

Modelling black carbon absorption of solar radiation: combining external and internal mixing assumptions.

An accurate simulation of the absorption properties is key for assessing the radiative effects of aerosol on meteorology and climate. The representation of how chemical species are mixed inside the particles (the mixing state) is one of the major uncertainty factors in the assessment of these effects. Here we compare aerosol optical properties simulations over Europe and North America, coordinated in the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII), to 1 year of AERONET sunphotometer retrievals, in an attempt to identify a mixing state representation that better reproduces the observed single scattering albedo and its spectral variation. We use a single post-processing tool (FlexAOD) to derive aerosol optical properties from simulated aerosol speciation profiles, and focus on the absorption enhancement of black carbon when it is internally mixed with more scattering material, discarding from the analysis scenes dominated by dust. We found that the single scattering albedo at 440 nm (ω 0,440) is on average overestimated (underestimated) by 3-5 % when external (core-shell internal) mixing of particles is assumed, a bias comparable in magnitude with the typical variability of the quantity. The (unphysical) homogeneous internal mixing assumption underestimates ω 0,440 by ~ 14 %. The combination of external and core-shell configurations (partial internal mixing), parameterized using a simplified function of air mass aging, reduces the ω 0,440 bias to -1/-3 %. The black carbon absorption enhancement (E abs) in core-shell with respect to the externally mixed state is in the range 1.8-2.5, which is above the currently most accepted upper limit of ~ 1.5. The partial internal mixing reduces E abs to values more consistent with this limit. However, the spectral dependence of the absorption is not well reproduced, and the absorption Ångström exponent AAE 675 440 is overestimated by 70-120 %. Further testing against more comprehensive campaign data, including a full characterization of the aerosol profile in terms of chemical speciation, mixing state, and related optical properties, would help in putting a better constraint on these calculations.

Modeled deposition of nitrogen and sulfur in Europe estimated by 14 air quality model systems: evaluation, effects of changes in emissions and implications for habitat protection.

The evaluation and intercomparison of air quality models is key to reducing model errors and uncertainty. The projects AQMEII3 and EURODELTA-Trends, in the framework of the Task Force on Hemispheric Transport of Air Pollutants and the Task Force on Measurements and Modelling, respectively (both task forces under the UNECE Convention on the Long Range Transport of Air Pollution, LTRAP), have brought together various regional air quality models to analyze their performance in terms of air concentrations and wet deposition, as well as to address other specific objectives. This paper jointly examines the results from both project communities by intercomparing and evaluating the deposition estimates of reduced and oxidized nitrogen (N) and sulfur (S) in Europe simulated by 14 air quality model systems for the year 2010. An accurate estimate of deposition is key to an accurate simulation of atmospheric concentrations. In addition, deposition fluxes are increasingly being used to estimate ecological impacts. It is therefore important to know by how much model results differ and how well they agree with observed values, at least when comparison with observations is possible, such as in the case of wet deposition. This study reveals a large variability between the wet deposition estimates of the models, with some performing acceptably (according to previously defined criteria) and others underestimating wet deposition rates. For dry deposition, there are also considerable differences between the model estimates. An ensemble of the models with the best performance for N wet deposition was made and used to explore the implications of N deposition in the conservation of protected European habitats. Exceedances of empirical critical loads were calculated for the most common habitats at a resolution of 100 × 100 m2 within the Natura 2000 network, and the habitats with the largest areas showing exceedances are determined. Moreover, simulations with reduced emissions in selected source areas indicated a fairly linear relationship between reductions in emissions and changes in the deposition rates of N and S. An approximate 20 % reduction in N and S deposition in Europe is found when emissions at a global scale are reduced by the same amount. European emissions are by far the main contributor to deposition in Europe, whereas the reduction in deposition due to a decrease in emissions in North America is very small and confined to the western part of the domain. Reductions in European emissions led to substantial decreases in the protected habitat areas with critical load exceedances (halving the exceeded area for certain habitats), whereas no change was found, on average, when reducing North American emissions in terms of average values per habitat.

Influence of anthropogenic emissions and boundary conditions on multi-model simulations of major air pollutants over Europe and North America in the framework of AQMEII3.

In the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3), and as contribution to the second phase of the Hemispheric Transport of Air Pollution (HTAP2) activities for Europe and North America, the impacts of a 20 % decrease of global and regional anthropogenic emissions on surface air pollutant levels in 2010 are simulated by an international community of regional-scale air quality modeling groups, using different state-of-the-art chemistry and transport models (CTMs). The emission perturbations at the global level, as well as over the HTAP2-defined regions of Europe, North America and East Asia, are first simulated by the global Composition Integrated Forecasting System (C-IFS) model from European Centre for Medium-Range Weather Forecasts (ECMWF), which provides boundary conditions to the various regional CTMs participating in AQMEII3. On top of the perturbed boundary conditions, the regional CTMs used the same set of perturbed emissions within the regional domain for the different perturbation scenarios that introduce a 20 % reduction of anthropogenic emissions globally as well as over the HTAP2-defined regions of Europe, North America and East Asia. Results show that the largest impacts over both domains are simulated in response to the global emission perturbation, mainly due to the impact of domestic emission reductions. The responses of NO2, SO2 and PM concentrations to a 20 % anthropogenic emission reduction are almost linear (~ 20 % decrease) within the global perturbation scenario with, however, large differences in the geographical distribution of the effect. NO2, CO and SO2 levels are strongly affected over the emission hot spots. O3 levels generally decrease in all scenarios by up to ~ 1 % over Europe, with increases over the hot spot regions, in particular in the Benelux region, by an increase up to ~ 6 % due to the reduced effect of NOx titration. O3 daily maximum of 8 h running average decreases in all scenarios over Europe, by up to ~ 1 %. Over the North American domain, the central-to-eastern part and the western coast of the US experience the largest response to emission perturbations. Similar but slightly smaller responses are found when domestic emissions are reduced. The impact of intercontinental transport is relatively small over both domains, however, still noticeable particularly close to the boundaries. The impact is noticeable up to a few percent, for the western parts of the North American domain in response to the emission reductions over East Asia. O3 daily maximum of 8 h running average decreases in all scenarios over north Europe by up to ~ 5 %. Much larger reductions are calculated over North America compared to Europe. In addition, values of the Response to Extra-Regional Emission Reductions (RERER) metric have been calculated in order to quantify the differences in the strengths of nonlocal source contributions to different species among the different models. We found large RERER values for O3 (~ 0.8) over both Europe and North America, indicating a large contribution from non-local sources, while for other pollutants including particles, low RERER values reflect a predominant control by local sources. A distinct seasonal variation in the local vs. non-local contributions has been found for both O3 and PM2.5, particularly reflecting the springtime long-range transport to both continents.

Assessment and economic valuation of air pollution impacts on human health over Europe and the United States as calculated by a multi-model ensemble in the framework of AQMEII3.

The impact of air pollution on human health and the associated external costs in Europe and the United States (US) for the year 2010 are modeled by a multi-model ensemble of regional models in the frame of the third phase of the Air Quality Modelling Evaluation International Initiative (AQMEII3). The modeled surface concentrations of O3, CO, SO2 and PM2.5 are used as input to the Economic Valuation of Air Pollution (EVA) system to calculate the resulting health impacts and the associated external costs from each individual model. Along with a base case simulation, additional runs were performed introducing 20 % anthropogenic emission reductions both globally and regionally in Europe, North America and east Asia, as defined by the second phase of the Task Force on Hemispheric Transport of Air Pollution (TF-HTAP2). Health impacts estimated by using concentration inputs from different chemistry-transport models (CTMs) to the EVA system can vary up to a factor of 3 in Europe (12 models) and the United States (3 models). In Europe, the multi-model mean total number of premature deaths (acute and chronic) is calculated to be 414 000, while in the US, it is estimated to be 160 000, in agreement with previous global and regional studies. The economic valuation of these health impacts is calculated to be EUR 300 billion and 145 billion in Europe and the US, respectively. A subset of models that produce the smallest error compared to the surface observations at each time step against an all-model mean ensemble results in increase of health impacts by up to 30 % in Europe, while in the US, the optimal ensemble mean led to a decrease in the calculated health impacts by ~ 11 %. A total of 54 000 and 27 500 premature deaths can be avoided by a 20 % reduction of global anthropogenic emissions in Europe and the US, respectively. A 20 % reduction of North American anthropogenic emissions avoids a total of ~ 1000 premature deaths in Europe and 25 000 total premature deaths in the US. A 20 % decrease of anthropogenic emissions within the European source region avoids a total of 47 000 premature deaths in Europe. Reducing the east Asian anthropogenic emissions by 20 % avoids ~ 2000 total premature deaths in the US. These results show that the domestic anthropogenic emissions make the largest impacts on premature deaths on a continental scale, while foreign sources make a minor contribution to adverse impacts of air pollution.

Evaluation and uncertainty estimation of the impact of air quality modelling on crop yields and premature deaths using a multi-model ensemble.

This study promotes the critical use of air pollution modelling results for health and agriculture impacts, with the primary goal of providing more reliable estimates to decision makers. To date, the accuracy of air quality (AQ) models and the effects of model-to-model result variability (which we will refer to as model uncertainty) on impact assessment studies have been often ignored, thus undermining the robustness of the information used in the decision making process and the confidence in the results obtained. A suite of twelve PM2.5 and ozone concentration fields produced by regional-scale chemistry transport Air Quality (AQ) models during the third phase of the Air Quality Model Evaluation International Initiative (AQMEII) has been used to calculate the impact of air pollution on premature deaths and crop yields. An innovative technique is applied to bias-adjust the models to available observations. The model results for ozone and PM2.5 are combined in a multi-model (MM) ensemble, which is used to estimate the damage and economic cost to human health and crop yields, as well as the associated uncertainties. The MM ensemble quantifies directly the uncertainty introduced by AQ models into the air pollution impact assessment chain, while the indirect use of experimental information through a bias adjustment, reduces the uncertainty in the ozone and PM2.5 fields and subsequently the uncertainty of the final impact assessment and cost valuation. The analysis over the European countries analysed in this study shows a mean number of premature deaths due to exposure to PM2.5 and ozone of approximately 370,000 (inter-quantile range between 260,000 and 415,000) and a relative yield loss of approximately 7% to 9% (depending on the exposure metrics used, for wheat and maize together). Furthermore, the results indicate that a reduction in the uncertainty of the modelled ozone by 61% and by 80% (depending on the aggregation metric used) and by 46% for PM2.5, produces a reduction in the uncertainty in premature mortality and crop loss of >60%, and of an equivalent percentage in the final uncertainty of cost valuation, providing decision makers with more accurate estimations for more targeted interventions.

Two-scale multi-model ensemble: is a hybrid ensemble of opportunity telling us more?

In this study we introduce a hybrid ensemble consisting of air quality models operating at both the global and regional scale. The work is motivated by the fact that these different types of models treat specific portions of the atmospheric spectrum with different levels of detail, and it is hypothesized that their combination can generate an ensemble that performs better than mono-scale ensembles. A detailed analysis of the hybrid ensemble is carried out in the attempt to investigate this hypothesis and determine the real benefit it produces compared to ensembles constructed from only global-scale or only regional-scale models. The study utilizes 13 regional and 7 global models participating in the Hemispheric Transport of Air Pollutants phase 2 (HTAP2)-Air Quality Model Evaluation International Initiative phase 3 (AQMEII3) activity and focuses on surface ozone concentrations over Europe for the year 2010. Observations from 405 monitoring rural stations are used for the evaluation of the ensemble performance. The analysis first compares the modelled and measured power spectra of all models and then assesses the properties of the mono-scale ensembles, particularly their level of redundancy, in order to inform the process of constructing the hybrid ensemble. This study has been conducted in the attempt to identify that the improvements obtained by the hybrid ensemble relative to the mono-scale ensembles can be attributed to its hybrid nature. The improvements are visible in a slight increase of the diversity (4 % for the hourly time series, 10 % for the daily maximum time series) and a smaller improvement of the accuracy compared to diversity. Root mean square error (RMSE) improved by 13-16 % compared to G and by 2-3 % compared to R. Probability of detection (POD) and false-alarm rate (FAR) show a remarkable improvement, with a steep increase in the largest POD values and smallest values of FAR across the concentration ranges. The results show that the optimal set is constructed from an equal number of global and regional models at only 15 % of the stations. This implies that for the majority of the cases the regional-scale set of models governs the ensemble. However given the high degree of redundancy that characterizes the regional-scale models, no further improvement could be expected in the ensemble performance by adding yet more regional models to it. Therefore the improvement obtained with the hybrid set can confidently be attributed to the different nature of the global models. The study strongly reaffirms the importance of an in-depth inspection of any ensemble of opportunity in order to extract the maximum amount of information and to have full control over the data used in the construction of the ensemble.

Seasonal ozone vertical profiles over North America using the AQMEII3 group of air quality models: model inter-comparison and stratospheric intrusions.

This study evaluates simulated vertical ozone profiles produced in the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) against ozonesonde observations in North America for the year 2010. Four research groups from the United States (US) and Europe have provided modeled ozone vertical profiles to conduct this analysis. Because some of the modeling systems differ in their meteorological drivers, wind speed and temperature are also included in the analysis. In addition to the seasonal ozone profile evaluation for 2010, we also analyze chemically inert tracers designed to track the influence of lateral boundary conditions on simulated ozone profiles within the modeling domain. Finally, cases of stratospheric ozone intrusions during May-June 2010 are investigated by analyzing ozonesonde measurements and the corresponding model simulations at Intercontinental Chemical Transport Experiment Ozonesonde Network Study (IONS) experiment sites in the western United States. The evaluation of the seasonal ozone profiles reveals that, at a majority of the stations, ozone mixing ratios are underestimated in the 1-6 km range. The seasonal change noted in the errors follows the one seen in the variance of ozone mixing ratios, with the majority of the models exhibiting less variability than the observations. The analysis of chemically inert tracers highlights the importance of lateral boundary conditions up to 250 hPa for the lower-tropospheric ozone mixing ratios (0-2 km). Finally, for the stratospheric intrusions, the models are generally able to reproduce the location and timing of most intrusions but underestimate the magnitude of the maximum mixing ratios in the 2-6 km range and overestimate ozone up to the first kilometer possibly due to marine air influences that are not accurately described by the models. The choice of meteorological driver appears to be a greater predictor of model skill in this altitude range than the choice of air quality model.

Technical note: Coordination and harmonization of the multi-scale, multi-model activities HTAP2, AQMEII3, and MICS-Asia3: simulations, emission inventories, boundary conditions, and model output formats.

We present an overview of the coordinated global numerical modelling experiments performed during 2012-2016 by the Task Force on Hemispheric Transport of Air Pollution (TF HTAP), the regional experiments by the Air Quality Model Evaluation International Initiative (AQMEII) over Europe and North America, and the Model Intercomparison Study for Asia (MICS-Asia). To improve model estimates of the impacts of intercontinental transport of air pollution on climate, ecosystems, and human health and to answer a set of policy-relevant questions, these three initiatives performed emission perturbation modelling experiments consistent across the global, hemispheric, and continental/regional scales. In all three initiatives, model results are extensively compared against monitoring data for a range of variables (meteorological, trace gas concentrations, and aerosol mass and composition) from different measurement platforms (ground measurements, vertical profiles, airborne measurements) collected from a number of sources. Approximately 10 to 25 modelling groups have contributed to each initiative, and model results have been managed centrally through three data hubs maintained by each initiative. Given the organizational complexity of bringing together these three initiatives to address a common set of policy-relevant questions, this publication provides the motivation for the modelling activity, the rationale for specific choices made in the model experiments, and an overview of the organizational structures for both the modelling and the measurements used and analysed in a number of modelling studies in this special issue.

Advanced error diagnostics of the CMAQ and Chimere modelling systems within the AQMEII3 model evaluation framework.

The work here complements the overview analysis of the modelling systems participating in the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) by focusing on the performance for hourly surface ozone by two modelling systems, Chimere for Europe and CMAQ for North America. The evaluation strategy outlined in the course of the three phases of the AQMEII activity, aimed to build up a diagnostic methodology for model evaluation, is pursued here and novel diagnostic methods are proposed. In addition to evaluating the "base case" simulation in which all model components are configured in their standard mode, the analysis also makes use of sensitivity simulations in which the models have been applied by altering and/or zeroing lateral boundary conditions, emissions of anthropogenic precursors, and ozone dry deposition. To help understand of the causes of model deficiencies, the error components (bias, variance, and covariance) of the base case and of the sensitivity runs are analysed in conjunction with timescale considerations and error modelling using the available error fields of temperature, wind speed, and NO x concentration. The results reveal the effectiveness and diagnostic power of the methods devised (which remains the main scope of this study), allowing the detection of the timescale and the fields that the two models are most sensitive to. The representation of planetary boundary layer (PBL) dynamics is pivotal to both models. In particular, (i) the fluctuations slower than ~ 1.5 days account for 70-85 % of the mean square error of the full (undecomposed) ozone time series; (ii) a recursive, systematic error with daily periodicity is detected, responsible for 10-20 % of the quadratic total error; (iii) errors in representing the timing of the daily transition between stability regimes in the PBL are responsible for a covariance error as large as 9 ppb (as much as the standard deviation of the network-average ozone observations in summer in both Europe and North America); (iv) the CMAQ ozone error has a weak/negligible dependence on the errors in NO2, while the error in NO2 significantly impacts the ozone error produced by Chimere; (v) the response of the models to variations of anthropogenic emissions and boundary conditions show a pronounced spatial heterogeneity, while the seasonal variability of the response is found to be less marked. Only during the winter season does the zeroing of boundary values for North America produce a spatially uniform deterioration of the model accuracy across the majority of the continent.

Evaluation and error apportionment of an ensemble of atmospheric chemistry transport modeling systems: multivariable temporal and spatial breakdown.

Through the comparison of several regional-scale chemistry transport modeling systems that simulate meteorology and air quality over the European and North American continents, this study aims at (i) apportioning error to the responsible processes using timescale analysis, (ii) helping to detect causes of model error, and (iii) identifying the processes and temporal scales most urgently requiring dedicated investigations. The analysis is conducted within the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII) and tackles model performance gauging through measurement-to-model comparison, error decomposition, and time series analysis of the models biases for several fields (ozone, CO, SO2, NO, NO2, PM10, PM2.5, wind speed, and temperature). The operational metrics (magnitude of the error, sign of the bias, associativity) provide an overallsense of model strengths and deficiencies, while apportioning the error to its constituent parts (bias, variance, and covariance) can help assess the nature and quality of the error. Each of the error components is analyzed independently and apportioned to specific processes based on the corresponding timescale (long scale, synoptic, diurnal, and intraday) using the error apportionment technique devised in the former phases of AQMEII. The application of the error apportionment method to the AQMEII Phase 3 simulations provides several key insights. In addition to reaffirming the strong impact of model inputs (emission and boundary conditions) and poor representation of the stable boundary layer on model bias, results also highlighted the high interdependencies among meteorological and chemical variables, as well as among their errors. This indicates that the evaluation of air quality model performance for individual pollutants needs to be supported by complementary analysis of meteorological fields and chemical precursors to provide results that are more insightful from a model development perspective. This will require evaluaion methods that are able to frame the impact on error of processes, conditions, and fluxes at the surface. For example, error due to emission and boundary conditions is dominant for primary species (CO, particulate matter (PM)), while errors due to meteorology and chemistry are most relevant to secondary species, such as ozone. Some further aspects emerged whose interpretation requires additional consideration, such as the uniformity of the synoptic error being region- and model-independent, observed for several pollutants; the source of unexplained variance for the diurnal component; and the type of error caused by deposition and at which scale.

World Meteorological Organization's model simulations of the radionuclide dispersion and deposition from the Fukushima Daiichi nuclear power plant accident.

Five different atmospheric transport and dispersion model's (ATDM) deposition and air concentration results for atmospheric releases from the Fukushima Daiichi nuclear power plant accident were evaluated over Japan using regional (137)Cs deposition measurements and (137)Cs and (131)I air concentration time series at one location about 110 km from the plant. Some of the ATDMs used the same and others different meteorological data consistent with their normal operating practices. There were four global meteorological analyses data sets available and two regional high-resolution analyses. Not all of the ATDMs were able to use all of the meteorological data combinations. The ATDMs were configured identically as much as possible with respect to the release duration, release height, concentration grid size, and averaging time. However, each ATDM retained its unique treatment of the vertical velocity field and the wet and dry deposition, one of the largest uncertainties in these calculations. There were 18 ATDM-meteorology combinations available for evaluation. The deposition results showed that even when using the same meteorological analysis, each ATDM can produce quite different deposition patterns. The better calculations in terms of both deposition and air concentration were associated with the smoother ATDM deposition patterns. The best model with respect to the deposition was not always the best model with respect to air concentrations. The use of high-resolution mesoscale analyses improved ATDM performance; however, high-resolution precipitation analyses did not improve ATDM predictions. Although some ATDMs could be identified as better performers for either deposition or air concentration calculations, overall, the ensemble mean of a subset of better performing members provided more consistent results for both types of calculations.

Comparison of the complex terrain algorithms incorporated into two commonly used local-scale air pollution dispersion models (ADMS and AERMOD) using a hybrid model.

ADMS and AERMOD are the two most widely used dispersion models for regulatory purposes. It is, therefore, important to understand the differences in the predictions of the models and the causes of these differences. The treatment by the models of flat terrain has been discussed previously; in this paper the focus is on their treatment of complex terrain. The paper includes a discussion of the impacts of complex terrain on airflow and dispersion and how these are treated in ADMS and AERMOD, followed by calculations for two distinct cases: (i) sources above a deep valley within a relatively flat plateau area (Clifty Creek power station, USA); (ii) sources in a valley in hilly terrain where the terrain rises well above the stack tops (Ribblesdale cement works, England). In both cases the model predictions are markedly different. At Clifty Creek, ADMS suggests that the terrain markedly increases maximum surface concentrations, whereas the AERMOD complex terrain module has little impact. At Ribblesdale, AERMOD predicts very large increases (a factor of 18) in the maximum hourly average surface concentrations due to plume impaction onto the neighboring hill; although plume impaction is predicted by ADMS, the increases in concentration are much less marked as the airflow model in ADMS predicts some lateral deviation of the streamlines around the hill.