APPLICATION OF CMAQ ON HEMISPHERIC SCALES

Similar documents
CMAQ Modeling of Atmospheric Mercury

The North American Mercury Model Intercomparison Study (NAMMIS): Study description and model-to-model comparisons

Chemical Transport of Atmospheric Mercury over the Pacific

New Science Implementation in CMAQ-Hg: Test over a Continental United States Domain

ATMOSPHERIC MERCURY SIMULATION WITH CMAQ VERSION 4.5.1

Supplement of Photochemical grid model implementation and application of VOC, NO x, and O 3 source apportionment

Global Mercury Modeling at Environment Canada. Ashu Dastoor &Didier Davignon. Air Quality Research Division Environment Canada

Jesper H. Christensen NERI-ATMI, Frederiksborgvej Roskilde

PROBLEMS Sources of CO Sources of tropospheric ozone

Know and Respond AQ Alert Service. Paul Willis SCOTTISH AIR QUALITY DATABASE AND WEBSITE ANNUAL SEMINAR Stirling 30 th March 2011

HTAP-2 analysis for the Arctic

Atmospheric Composition Matters: To Air Quality, Weather, Climate and More

AIRQUEST Annual Report and State of the Model

Reactive mercury in the troposphere: Model formation and results for Florida, the northeastern U.S. and the Atlantic Ocean

MULTI-MODEL AIR QUALITY FORECASTING OVER NEW YORK STATE FOR SUMMER 2008

HIGH RESOLUTION CMAQ APPLICATION FOR THE REGIONAL MUNICIPALITY OF PEEL, ONTARIO, CANADA

Plume-in-grid Modeling of Atmospheric Mercury Deposition in the United States

Implications of Sulfate Aerosols on Clouds, Precipitation and Hydrological Cycle

Improvement of Meteorological Inputs for Air Quality Study

Vertical Structure of Atmosphere

Science Results Based on Aura OMI-MLS Measurements of Tropospheric Ozone and Other Trace Gases

Tananyag fejlesztés idegen nyelven

MCIPv3: Using WRF-EM Output with CMAQ

Updated Dust-Iron Dissolution Mechanism: Effects Of Organic Acids, Photolysis, and Dust Mineralogy

Ammonia Emissions and Nitrogen Deposition in the United States and China

TM4-ECPL model : Oceanic Sources for Oxygenated VOC and Aerosols

Accuracy and Cost Considerations in Choosing a Chemical Mechanism for Operational Use in AQ Models

Atmospheric mercury chemistry: Sensitivity of global model simulations to chemical reactions

CONTENTS 1 MEASURES OF ATMOSPHERIC COMPOSITION

3. Carbon Dioxide (CO 2 )

The Importance of Ammonia in Modeling Atmospheric Transport and Deposition of Air Pollution. Organization of Talk:

Creating Meteorology for CMAQ

surrounds Earth and protects it somewhat from solar radiation. Like all other matter, air has weight,

NO X emissions, isoprene oxidation pathways, and implications for surface ozone in the Southeast United States

Chemical sources and sinks of Hg(II) in the remote atmospheric marine boundary layer

CHAPTER 1. MEASURES OF ATMOSPHERIC COMPOSITION

TIME SERIES COMPARISONS OF MIPAS LEVEL 2 NEAR REAL TIME PRODUCTS WITH CLIMATOLOGY

The Challenge of. Guy Brasseur

The North Atlantic Oscillation: Climatic Significance and Environmental Impact

Incorporating Space-borne Observations to Improve Biogenic Emission Estimates in Texas (Project )

Transport of Asian ozone pollution into surface air over the western U.S. in spring. Meiyun Lin

PRELIMINARY EXPERIENCES WITH THE MULTI-MODEL AIR QUALITY FORECASTING SYSTEM FOR NEW YORK STATE

Air Quality Modelling under a Future Climate

CMAQ and CAMx Overview and Recent Updates 2014 Midwest and Central States Air Quality Workshop

14.4 NUMERICAL SIMULATION OF AIR POLLUTION OVER KANTO AREA IN JAPAN USING THE MM5/CMAQ MODEL

Baseline Ozone in Western North America: Measurements and Models. David Parrish

Responsibilities of Harvard Atmospheric Chemistry Modeling Group

CHAPTER 8. AEROSOLS 8.1 SOURCES AND SINKS OF AEROSOLS

OVERVIEW OF CMAQ 5.0 AND CAMX 5.4 5/17/2012 1

ASSESSMENT OF PM2.5 RETRIEVALS USING A COMBINATION OF SATELLITE AOD AND WRF PBL HEIGHTS IN COMPARISON TO WRF/CMAQ BIAS CORRECTED OUTPUTS

A B C D PROBLEMS Dilution of power plant plumes. z z z z

Regional modelling and assessment of atmospheric particulate matter concentrations at rural background locations in Europe

Direct radiative forcing due to aerosols in Asia during March 2002

Arctic Chemistry And Climate

Review of the SAPRC-16 Chemical Mechanism and Comparison with the Regional Atmospheric Chemistry Mechanism, Version-2

Eurasian Snow Cover Variability and Links with Stratosphere-Troposphere Coupling and Their Potential Use in Seasonal to Decadal Climate Predictions

Surface Ozone Problem in Two Polluted Regions in China and VOGA-NCP 2013 Summer Campaign

Estimation of Speciated and Total Mercury Dry Deposition. Leiming Zhang. Pierrette Blanchard, David Gay, and many others

SCICHEM: A Puff Model with Chemistry. Part 2: Ozone and Particulate Matter

Extremes of Weather and the Latest Climate Change Science. Prof. Richard Allan, Department of Meteorology University of Reading

Atmospheric models. Chem Christoph Knote - NCAR Spring 2013

Using GOME and SCIAMACHY NO 2 measurements to constrain emission inventories potential and limitations

Chapter outline. Reference 12/13/2016

Atmospheric composition modeling over the Arabian Peninsula for Solar Energy applications

Supporting Online Material for

Data Assimilation for Tropospheric CO. Avelino F. Arellano, Jr. Atmospheric Chemistry Division National Center for Atmospheric Research

Supplementary Information for:

CALIOPE EU: Air Quality

Mae Gustin University of Nevada

A STUDY OF THE IMPACT OF SYNOPTIC WEATHER CONDITIONS AND WATER VAPOUR ON AEROSOL-CLOUD RELATIONSHIPS

Arctic Oxidation Chemistry

Revisi&ng the GEOS- Chem atmospheric Hg simula&on: chemistry and emissions

Aircraft Observations of the Lower Atmosphere and Surface Exchange Processes. Jennifer Murphy Department of Chemistry, University of Toronto

Regional methane emissions estimates in northern Pennsylvania gas fields using a mesoscale atmospheric inversion system

Introduction to atmospheric dispersion modelling. M.Sofiev

Temperature (T) degrees Celsius ( o C) arbitrary scale from 0 o C at melting point of ice to 100 o C at boiling point of water Also (Kelvin, K) = o C

REGIONAL AIR QUALITY FORECASTING OVER GREECE WITHIN PROMOTE

Climate Change 2007: The Physical Science Basis

Source apportionment of fine particulate matter over the Eastern U.S. Part I. Source sensitivity simulations using CMAQ with the Brute Force method

NOAA s Air Quality Forecasting Activities. Steve Fine NOAA Air Quality Program

7.5 ASSESSMENT OF BASELINE AIR QUALITY IN WYOMING S GREEN RIVER BASIN PRIOR TO ACCELERATED GAS AND OIL WELL DEVELOPMENT

1. The frequency of an electromagnetic wave is proportional to its wavelength. a. directly *b. inversely

2. Fargo, North Dakota receives more snow than Charleston, South Carolina.

Comparing Modal and Sectional Approaches in Modeling Particulate Matter in Northern California

Air Quality Modeling from the Offshore Energy Sector in the Gulf of Mexico: An Overview for the Oil and Gas Industry

Biogenic aerosols and their interactions with climate. Yuzhong Zhang

SUPPLEMENTARY INFORMATION

Impact of Solar and Sulfate Geoengineering on Surface Ozone

Approach Estimating Mercury Dry Deposition for AMNeT Leiming Zhang

Rela%ve Merit of MODIS AOD and Surface PM2.5 for Aerosol Analysis and Forecast

Transport of stratospheric aerosols in the field of averaged vertical wind

Introduction. 1 Division of Engineering and Applied Sciences, Harvard University, Cambridge, MA

Ambient Ammonium Contribution to total Nitrogen Deposition


Challenges in the Use of Satellite Measurements for Studies of Tropospheric Ozone

Aerosol/Transport issues in VOCALS region

Long-lead prediction of the 2015 fire and haze episode in Indonesia

Application of an Adaptive Nudging Scheme in Air Quality Forecasting in China

Preliminary Experiences with the Multi Model Air Quality Forecasting System for New York State

ATOC 3500/CHEM 3151 Air Pollution Chemistry Lecture 1

Transcription:

APPLICATION OF CMAQ ON HEMISPHERIC SCALES O. Russell Bullock, Jr.* and Rohit Mathur US Environmental Protection Agency, Research Triangle Park, NC, USA Francis S. Binkowski and Neil N. Davis Carolina Environmental Program, The University of North Carolina, Chapel Hill, NC, USA 1. INTRODUCTION The Community Multi-scale Air Quality (CMAQ) model has been applied to simulate a number of different air pollutants by environmental institutions and agencies all around the world. These applications have generally been on limited-area domains where estimation of the air concentration of pertinent pollutant species at the lateral boundaries is required. Results from the recent North American Mercury Model Intercomparison Study (NAMMIS) (Bullock et al., 2008; 2009) suggest that prescribed lateral boundary air concentrations can lead to unrealistic simulation results if they are not consistent with the meteorological and chemical conditions of the CMAQ simulation. For example, if water-soluble pollutants are prescribed to exist at a boundary location where the CMAQ simulation indicates heavy rainfall, unrealistically large deposition fluxes will be simulated. Mercury (Hg) is recognized as a global pollutant with a typical atmospheric lifetime on the order of months or years. However, certain forms of atmospheric Hg (e.g., oxidized gases) are known to deposit rather rapidly once they are emitted or formed in atmospheric chemical reactions. Thus, atmospheric Hg modeling has both long-range (global) and short-range (local/regional) components and provides an opportunity to investigate the effect of prescribed lateral boundary values that may be inconsistent with internal model processes. To avoid inconsistent lateral boundary concentrations, this work has applied the CMAQ model on a domain covering the Northern Hemisphere to the extent that boundary fluxes are much less of a concern, at least for Hg species. We plan to eventually use hemispheric simulations and domain nesting to regional scales to study the effect of inconsistent boundary concentrations for a variety of pollutants. This paper describes our *Corresponding author: O. Russell Bullock, Jr., U.S. Environmental Protection Agency, Mail Drop E243-01, Research Triangle Park, NC 27711; e-mail: bullock.russell@epa.gov initial 15-day hemispheric application of CMAQ for atmospheric Hg, which also involved simulation of ozone, acid deposition and particulate matter. 2. MODELING DOMAIN The CMAQ was designed to work on any conformal map projection (i.e., where geometric shape is conserved). A polar stereographic projection was chosen to cover most of the Northern Hemisphere. We wanted to get far enough south so that boundary fluxes of Hg in CMAQ simulations would be significantly reduced. There was some concern that the meteorology of the Inter-Tropical Convergence Zone (ITCZ) near the Equator might be difficult to simulate. The US EPA worked with the UNC Institute for the Environment (UNC-IE) to develop a meteorological modeling domain covering all of the Northern Hemisphere north of 10 N. This was accomplished using a 149x149 array of 108-km grid cells on a secant latitude of 45 N. To isolate boundary effects in the meteorological modeling, the CMAQ domain was set back by 6 grid cells on all sides, resulting in a 137x137 grid cell array with corner cells near the equator and mid-lateral cells near 20 N. The CMAQ vertical domain was resolved using 34 layers defined in sigmapressure coordinates with a model top at 100 hpa. 3. METEOROLOGICAL DATA The meteorological input data for this study were produced by UNC-IE using version 3.0 of the Weather Research and Forecasting (WRF) - Advanced Research WRF (WRF-ARW) model. The WRF-ARW simulations used the 149x149 grid cell array described above with the same 34-layer vertical resolution as in CMAQ. The Pleim-Xiu (PX) surface layer, land cover and PBL schemes were used along with the RRTMG radiative transfer model. The WRF-ARW outputs were processed using MCIP version 3.4.1. 1

UNC-IE encountered some difficulties when it first tried to process WRF-ARW outputs using MCIP and display the results with PAVE and VERDI. The problem was traced to the way MCIP was converting WRF header data to the I/O API format. MCIP has been modified to correct this problem. When a major bug-fix to WRF-ARW was published in late 2008, EPA and UNC-IE decided it was prudent to restart the meteorological modeling with the corrected WRF- ARW. UNC-IE eventually provided EPA with CMAQ-ready hemispheric meteorological data for the period from March 16 through 30, 2006. 4. POLLUTANT EMISSION DATA Emission data for anthropogenic SOx, NOx, VOC, PM, and CO were obtained from an on-line data base developed by NASA s Arctic Research of the Composition of the Troposphere from Aircraft and Satellite (ARCTAS) project. This data base is global and provides 1 x1 spatial resolution. The data are in terms of annual fluxes, so no temporal resolution is provided and none was estimated for this work. The base year for this inventory is 2005 or 2006 depending on the pollutant. Hg emissions are also available from the ARCTAS data base, but without any physical or chemical speciation. Emission data for biogenic isoprene and terpene were obtained from an on-line database provided by the Precursors of Ozone and their Effects on the Troposphere (POET) project. These data also have a 1 x1 spatial resolution. Biogenic emissions from POET are specified on a monthly basis and the base year for the data is 2000. Emission data for state-of-the-science atmospheric Hg modeling must provide information about the physical form (gas vs. aerosol) and the oxidation state for all significant Hg emissions. The first-ever global Hg emission inventory with this kind of speciation was recently developed by the Norwegian Institute for Air Research (NILU), the Arctic Monitoring and Assessment Programme (AMAP), and the Arctic Center, University of Groningen. It is available online providing 0.5 x0.5 spatial resolution and speciation in terms of gaseous elemental Hg (GEM), gaseous oxidized Hg (GOM) and particulate Hg (PHg). Some resolution of emission height based on industrial sector is also provided, but that was not used in this work. The base year for this inventory is 2000. 5. SIMULATION RESULTS Fig. 1. Vertical concentration profiles of gaseous elemental mercury (GEM) and gaseous oxidized mercury (GOM) used for initial and lateral boundary conditions. We applied a special multi-pollutant configuration of CMAQ version 4.7 where Hg was the only air toxic pollutant added to the base model for photo-oxidants, sulfur, nitrogen and particulate matter. As previously mentioned, CMAQ could only be applied on the hemispheric domain for a two-week period in late March of 2006. However, this simulation period was sufficient to demonstrate intercontinental transport of Hg and also highlight some additional technical obstacles that need to be overcome in the CMAQ treatment of atmospheric Hg cycling on hemispheric and global scales. 15-day simulations have been conducted using two different vertical profiles for initial and boundary concentrations for Hg species. Fig. 1 shows the concentration profiles used for GEM and GOM. The first set of vertical profiles had a discontinuous and modest decrease of GEM with height and a corresponding increase of GOM with height. The second set had profiles that were continuous with height, but the decrease of GEM and increase of GOM with height was more pronounced. In both cases, the initial and boundary concentrations of PHg were set to 17 pg m -3 at all levels. Two additional hemispheric simulations have been conducted using the second set of vertical profiles to test model sensitivity to changes in gaseous and aqueous Hg chemistry. The purpose for these tests was to determine the source of high GOM concentrations simulated near the surface in remote oceanic locations. 2

Fig. 2. Simulated concentration of gaseous oxidized mercury in layer 26 at 24-hour intervals showing trans-pacific transport from Asia to North America. 5.1 Intercontinental Transport Atmospheric Hg emissions in Asia have often been cited as a significant source of Hg deposition to the U.S. It has long been known that emissions of GEM all around the globe contribute to a nearly ubiquitous air concentration of 1.0 to 1.5 ng m -3 for that particular Hg species. However, atmospheric deposition of GEM is rather slow compared to the other species of Hg. The main concern is the possibility for emissions of GOM (and possibly PHg) to be transported across the Pacific Ocean with subsequent deposition to North America. Because GOM is readily deposited from air to any underlying surface, its transport near the surface should be terminated rather quickly. Thus, long-range transport of GOM at altitudes above the planetary boundary layer has always been the main concern. Fig. 2 shows results from the first CMAQ simulation which demonstrates transport across the Pacific of an air mass laden with ROM. Other less obvious examples of long-range GOM transport can be found in a variety of locations. The results shown here are from layer 26 of the model, which corresponds to an altitude of around 6000 m above sea level. 3

Fig. 3. Map of GOM concentrations in layer 26 near the end of the first simulation showing the influx of low concentrations at the lateral boundaries. Fig. 4. Map of GOM concentrations in layer 26 near the end of the second simulation showing boundary influx concentrations similar to those in the interior. 5.2 Boundary effects Measurements and global model simulations have generally indicated higher concentrations of GOM at higher altitudes with an opposite relationship for GEM. The first set of initial and boundary values was designed to reflect this. However, towards the end of the first 15-day simulation, it became apparent that the CMAQ model was not simulating the same general vertical distribution of GEM vs. GOM as was defined in that first set of boundary values. Fig. 3 shows an obvious lateral boundary influx of low GOM concentrations in layer 26 (~6000 m ASL) compared to what CMAQ simulated for interior portions of the domain. Three separate global models used to define boundary conditions for a recent atmospheric mercury modeling study indicated GOM concentrations of up to 700 pg m -3 (8 10-8 ppmv) near the 100 hpa level (Bullock et al., 2008). The second set of initial and boundary concentrations was developed to better match these values. Fig. 4 shows the result of using the second set of boundary values. Note that the concentration scale has been doubled from that used in Fig. 3. Boundary influx concentrations are similar to those within the interior of the domain at the end of the 15-day simulation. However, in this case the boundary values were found to be too high in the top model layer. CMAQ did not simulate a continuing increase of GOM to the top of the model like the global models. This may be because CMAQ does not extend as high to simulate the effect of stratospheric ozone on mercury oxidation as the global models do. 5.3 Surface-level RGM in marine air In the first two simulations, some rather high concentrations of GOM were noted near the surface over remote oceanic areas. Attendant low GEM concentration suggested some sort of chemical oxidation was the cause. A test was conducted where the simulated reaction of GEM with atomic chlorine (added in CMAQ v4.7) was slowed by a factor of 100. This had only a minor effect. In a second test, the aqueous chemistry was changed to allow the reaction of all dissolved species of Hg 2+ with HO 2. With CMAQ v4.7, only the free Hg 2+ ion was involved in this mercuryreducing reaction. However, the method employed in the only rate determination for this reaction (Pehkonen and Lin, 1998) used the changing concentration of total Hg 2+ to estimate the kinetic rate (Lin, personal communication). Figure 5 shows that reacting all dissolved Hg 2+ instead of only the free ion results in a significant reduction in GOM concentrations in marine air near the surface. 4

concentration profiles applied at the boundary need to be chosen with care. Further work is underway to push these boundaries even closer to the equator where cross-boundary flows should be reduced further. The lack of a continued upward gradient with height for simulated GOM near the top of model suggests that stratospheric chemistry needs to be included. Further work is underway to raise the model top to 50 hpa. Finally, this work has shown that certain changes to the aqueous mercury chemistry in CMAQ v4.7 may need to be reconsidered. The reduction of aqueous Hg 2+ to Hg 0 by reaction with HO 2 has been deemed unlikely to be important under normal atmospheric conditions (Gårdfeldt and Jonsson, 2003) and this reaction was scaled down considerably in CMAQ v4.7. However, it seems obvious from the results shown here that some similar reaction is required to maintain a reasonable reduction-oxidation balance for Hg. 7. REFERENCES Bullock, O.R., Jr., D. Atkinson, T. Braverman, K. Civerolo, A. Dastoor, D. Davignon, J-.Y. Ku, K. Lohman, T.C. Myers, R.J. Park, C. Seigneur, N.E. Selin, G. Sistla, K. Vijayaraghavan, 2009: An Analysis of Simulated Wet Deposition of Mercury from the North American Mercury Model Intercomparison Study (NAMMIS). J. of Geophys. Res., 114, D08301, doi:10.1029/2008jd011224. Fig. 5. Simulated GOM concentration in layer 1 after 15 days where (a) aqueous HO 2 reacts only with free Hg 2+ ions and (b) aqueous HO 2 reacts with all dissolved Hg 2+. 6. SUMMARY This work has demonstrated that CMAQ can be applied for hemispheric-scale modeling using a horizontal domain based on a polar stereographic projection. Simulations on this domain show evidence of intercontinental transport of oxidized forms of Hg. However, evaluation of simulated deposition of Hg and the fraction attributable to this intercontinental transport has not yet been performed. This work has also shown that, even with lateral boundaries between 0 and 20 N, crossboundary flows can still be significant and vertical Bullock, O.R., Jr., D. Atkinson, T. Braverman, K. Civerolo, A. Dastoor, D. Davignon, J-.Y. Ku, K. Lohman, T.C. Myers, R.J. Park, C. Seigneur, N.E. Selin, G. Sistla, K. Vijayaraghavan, 2008: The North American Mercury Model Intercomparison Study (NAMMIS): Study description and model-tomodel comparisons. J. of Geophys. Res., 113, D17310, doi:10.1029/2008jd009803. Gårdfeldt, K. and M. Jonsson, 2003: Is bimolecular reduction of Hg(II)-complexes possible in aqueous systems of environmental importance? J. Phys. Chem., 107, 4478-4482. Pehkonen, S.O. and C.-J. Lin, 1998: Aqueous photochemistry of mercury with organic acids. J. Air & Waste Manage. Assoc., 48, 144-150. 5