Mercury Modeling
EPA has been working on the development of atmospheric mercury models since the early 1990's when the Regional Lagrangian Model of Air Pollution (RELMAP) was adapted to simulate mercury in support of EPA’s Mercury Study Report to Congress. As the scientific understanding of atmospheric mercury continued to develop in the late 1990's, it became apparent that Lagrangian-type models, also known as “puff” models, would have difficulties simulating the complex chemical and physical interactions of mercury with other pollutants that were being discovered. Thus, EPA’s focus for atmospheric mercury model development was moved to the Community Multi-scale Air Quality model (CMAQ). The CMAQ simulates atmospheric processes within a 3-dimensional array of pre-defined finite volume elements and can model complex interactions between all of the pollutants that might exist within each volume element. The CMAQ was previously developed to simulate photochemical oxidants, acidic and nutrient pollutants, and aerosol particulate matter, all of which have been shown to interact with mercury in air and in cloud water and influence its deposition to sensitive aquatic ecosystems. The “multi-pollutant” approach of CMAQ where all pollutants are simulated together just as they exist in the real atmosphere is applied in atmospheric mercury modeling at EPA.
A number of modification were made to the standard CMAQ model to allow it to simulate atmospheric mercury which are described in detail in Bullock and Brehme (2002). Because new information about chemical and physical processes affecting atmospheric mercury is continually being published, refinement of the model code is an ongoing process. Further modification of the CMAQ-Hg chemical mechanisms for mercury in both the gaseous and aqueous phases is expected as additional chemical reactions are identified and studied. The latest public release of CMAQ provides the ability to simulate atmospheric mercury in the “multi-pollutant” version of the model. We found this to be the most efficient way to maintain and disseminate mercury simulation capabilities in CMAQ due to the increasing number of pollutants with which mercury is known to react.
EPA has participated in two major model intercomparison studies for atmospheric mercury. The first was the Intercomparison of Numerical Models for Long-Range Atmospheric Transport of Mercury sponsored by the European Monitoring and Evaluation Programme (EMEP) and organized by EMEP’s Meteorological Synthesizing Center - East in Moscow, Russia. The first phase of this EMEP study involved the simulation of mercury chemistry in a closed cloud volume given a variety of initial conditions and results were published in Ryaboshapko et al. (2002). The second phase of the study involved full-scale model simulations of the emission, transport, transformation, and deposition of mercury over Europe for two short periods of 10 to 14 days each. Model simulations were compared to field measurements of elemental mercury gas, reactive gaseous mercury and particulate mercury in air. Results from this phase of the study were reported in Ryaboshapko et al. (2007a). The third and final phase of the EMEP intercomparison involved model simulations for longer periods of time (up to one year) and comparisons to observations of the wet deposition of mercury. Results from this phase of the study were reported in Ryaboshapko et al. (2007b).
As the EMEP study was nearing completion, EPA organized a second mercury model intercomparison study, this time with a focus on North America. The North American Mercury Model Intercomparison Study (NAMMIS) took advantage of standardized weekly wet deposition samples taken by the Mercury Deposition Network (MDN) as described in Vermette et al. (1995) and separate event-based precipitation samples taken at Underhill, Vermont (Keeler et al., 2005). In addition to CMAQ, two other regional models were tested in the NAMMIS; the Regional Modeling System for Aerosols and Deposition (REMSAD) and the Trace Element Analysis Model (TEAM). All three models were each applied to simulate the entire year of 2001 three times, each time using initial condition and boundary conditions developed from a different global model. The NAMMIS provided not only a comparison between regional atmospheric mercury models, but also a measure of the sensitivity of each regional model to uncertainties regarding intercontinental transport. The NAMMIS evaluated each regional model for its agreement to observations of wet deposition of mercury from the MDN and Underhill observations.
Results from the NAMMIS statistical model evaluation are shown in the table below. For most of the evaluation metrics, CMAQ was found to have superior agreement to the observations. Results from all three of the regional models tested (CMAQ, REMSAD and TEAM) varied depending on the global model used to define lateral boundary concentrations for mercury (CTM, GEOS-Chem or GRAHM). All of the regional models used meteorological data provided by the MM5 model. Statistics for the precipitation data obtained from MM5 are also shown in the table. Obviously, the level of accuracy one can expect from the regional air quality models is limited by the accuracy of the input precipitation data. It does appear that the best performing air-quality simulations have about the same level of accuracy as the precipitation data provided to those simulations. Thus, it can be reasoned that significant improvements in the simulation of mercury wet deposition are contingent upon improvements in the modeling of physical meteorology. Complete descriptions of the NAMMIS study design, participating models, and modeling results are available in two articles published in the Journal of Geophysical Research (Bullock et al., 2008; Bullock et al., 2009).
CMAQ mercury modeling capabilities have been applied to support various EPA regulatory actions for mercury. They have also been used to provide information regarding mercury deposition from global background concentrations to tribal, state and regional environmental authorities in the development of their water quality protection strategies. EPA is currently working with the United Nations Environment Program towards the development of international treaties to reduce mercury emissions to the environment. EPA is working to expand the CMAQ modeling domain to cover the Northern Hemisphere. This will allow CMAQ to provide modeling assessments of the intercontinental transport of mercury and its importance as a global pollutant.
|
Table 1. Evaluation statistics from the North American Mercury Model Intercomparison Study
|
Contacts: Russell Bullock
Relevant Publications & Presentations:
-
Bullock, O.R., Jr. and Brehme, K.A. Atmospheric mercury simulation using the CMAQ model: formulation description and analysis of wet deposition results. Atmospheric Environment 36, 2135-2146 (2002).
-
Bullock, O.R., Jr., Atkinson, D., Braverman, T., Civerolo, K., Dastoor, A., Davignon, D., Ku, J-.Y., Lohman, K., Myers, T.C., Park, R.J., Seigneur, C., Selin, N.E., Sistla, G., and Vijayaraghavan, K. (2008) The North American Mercury Model Intercomparison Study (NAMMIS): Study description and model-to-model comparisons. Journal of Geophysical Research 113, D17310, doi:10.1029/2008JD009803 (2008).
-
Bullock, O.R., Jr., Atkinson, D., Braverman, T., Civerolo, K., Dastoor, A., Davignon, D., Ku, J-.Y., Lohman, K., Myers, T.C., Park, R.J., Seigneur, C., Selin, N.E., Sistla, G., and Vijayaraghavan, K. An Analysis of Simulated Wet Deposition of Mercury from the North American Mercury Model Intercomparison Study (NAMMIS). Journal of Geophysical Research 114, D08301, doi:10.1029/2008JD011224 (2009).
-
Keeler, G. J., Gratz, L., and Al-Wali, K. Influences on the long-term atmospheric mercury wet deposition at Underhill, Vermont, Ecotoxicology 14, 71-83 (2005).
Ryaboshapko, A., Bullock, R., Ebinghaus, R., Ilyin, I., Lohman, K., Munthe, J., Petersen, G., Seigneur, C., and Wängberg, I. Comparison of mercury chemistry models. Atmospheric Environment 36, 3881-3898 (2002). -
Ryaboshapko, A., Bullock, O.R., Christensen, J., Cohen, M., Dastoor, A., Ilyin, I., Petersen, G., Syrakov, D., Artz, R.S., Davignon, D., Draxler, R.R., and Munthe, J. Intercomparison study of atmospheric mercury models: 1. Comparison of models with short-term measurements, Science of the Total Environment 376(1-3), 228-240 (2007a).
-
Ryaboshapko, A., Bullock, O.R., Christensen, J., Cohen, M., Dastoor, A., Ilyin, I., Petersen, G., Syrakov, D., Travnikov, O., Artz, R.S., Davignon, D., Draxler, R.R., Munthe, J., and Pacyna, J. Intercomparison study of atmospheric mercury models: 2. Modelling results vs. long-term observations and comparison of country atmospheric balances, Science of the Total Environment 377(2-3), 319-333 (2007b).
-
Vermette, S., Lindberg, S., and Bloom, N. Field tests for a regional mercury deposition network—sampling design and preliminary test results, Atmospheric Environment 29, 1247–1251 (1995).
![[logo] US EPA](http://www.epa.gov/epafiles/images/logo_epaseal.gif)
