567
Views
2
CrossRef citations to date
0
Altmetric
Visibility 2022

PM2.5 in Carlsbad Caverns National Park: Composition, sources, and visibility impacts

, , ORCID Icon, , , , , & ORCID Icon show all
Pages 1201-1218 | Received 25 Jan 2022, Accepted 15 May 2022, Published online: 01 Nov 2022

Figures & data

Figure 1. Active oil and natural gas wells near Carlsbad Caverns National Park (CAVE) (HIFLD, Citation2019). The color bar gives the number of active O/G wells in the given grid cell. The grid cell resolution is 0.1 degree by 0.1 degree. Cells with 500 or more active wells are shown with the same color. National Park boundaries are shown in green. Larger cities in the surrounding area are marked with black dots and labeled. The monitoring site, in CAVE is shown as a yellow star. Coal-fired power plants in the region are shown as yellow plus signs (U.S. EIA, Citation2022).

Figure 1. Active oil and natural gas wells near Carlsbad Caverns National Park (CAVE) (HIFLD, Citation2019). The color bar gives the number of active O/G wells in the given grid cell. The grid cell resolution is 0.1 degree by 0.1 degree. Cells with 500 or more active wells are shown with the same color. National Park boundaries are shown in green. Larger cities in the surrounding area are marked with black dots and labeled. The monitoring site, in CAVE is shown as a yellow star. Coal-fired power plants in the region are shown as yellow plus signs (U.S. EIA, Citation2022).

Table 1. Measurements made during the CarCavAQS study are listed with their respective time-resolutions and measured species. Only instruments used in this work are listed. Abbreviations used for the measured species: sodium (Na+), ammonium (NH4+), calcium (Ca2+), potassium (K+), magnesium (Mg2+), chloride (Cl), nitrite (NO2), nitrate (NO3), sulfate (SO42-), organic carbon (OC), black carbon (BC), elemental carbon (EC), ammonia (NH3), nitric acid (HNO3), sulfur dioxide (SO2).

Figure 2. Time series of major PM2.5 species and CH4 measured at the monitoring site. (a.) BC measured by the Aethalometer at 880 nm, WSOC measured by the PILS-TOC, and PM2.5 measured by the TEOM. Note that there is a gap in the PILS-IC and PILS-TOC data at the beginning of the study period associated with a power outage. The BC and WSOC are plotted at a time-resolution of 2-minutes. TEOM total PM2.5 is plotted at an hourly resolution. (b.) The five PILS-IC species with the largest contributions to total PM2.5 are plotted at 30-minute time-resolution. (c.) CH4 measured by the Picarro at 15-second time-resolution. Shaded regions represent case study transport dates.

Figure 2. Time series of major PM2.5 species and CH4 measured at the monitoring site. (a.) BC measured by the Aethalometer at 880 nm, WSOC measured by the PILS-TOC, and PM2.5 measured by the TEOM. Note that there is a gap in the PILS-IC and PILS-TOC data at the beginning of the study period associated with a power outage. The BC and WSOC are plotted at a time-resolution of 2-minutes. TEOM total PM2.5 is plotted at an hourly resolution. (b.) The five PILS-IC species with the largest contributions to total PM2.5 are plotted at 30-minute time-resolution. (c.) CH4 measured by the Picarro at 15-second time-resolution. Shaded regions represent case study transport dates.

Figure 3. URG annular denuder/filter pack sampling shows the particle/gas phase partitioning of (a.) NH3/NH4+, (b.) HNO3/ NO3, and (c.) SO2/ SO42-. Particle species are shown in Orange, and gas species are shown in blue. There is one day of missing data, 2019–08-04, due to a denuder extraction error.

Figure 3. URG annular denuder/filter pack sampling shows the particle/gas phase partitioning of (a.) NH3/NH4+, (b.) HNO3/ NO3−, and (c.) SO2/ SO42-. Particle species are shown in Orange, and gas species are shown in blue. There is one day of missing data, 2019–08-04, due to a denuder extraction error.

Figure 4. PILS-IC species (a.) NH4+ and SO42-, and (b.) NO3 and Na+ are shown as nano equivalents m−3. Shaded regions represent case study transport dates.

Figure 4. PILS-IC species (a.) NH4+ and SO42-, and (b.) NO3− and Na+ are shown as nano equivalents m−3. Shaded regions represent case study transport dates.

Figure 5. The ratio of cation to anion species measured by the PILS-IC is plotted in equiv. m−3. The dotted gray line shows a ratio of unity. The radius of the Orange point size changes proportionally with Ca2+ concentration. Shaded regions represent case study transport dates.

Figure 5. The ratio of cation to anion species measured by the PILS-IC is plotted in equiv. m−3. The dotted gray line shows a ratio of unity. The radius of the Orange point size changes proportionally with Ca2+ concentration. Shaded regions represent case study transport dates.

Figure 6. 48-hour HYSPLIT Back Trajectories are shown for three different transport patterns during the study. The bottom panel shows back trajectory altitudes above ground level. The dates were selected based on aerosol composition. For each date 12 daytime back trajectories of 48-hour duration are plotted.

Figure 6. 48-hour HYSPLIT Back Trajectories are shown for three different transport patterns during the study. The bottom panel shows back trajectory altitudes above ground level. The dates were selected based on aerosol composition. For each date 12 daytime back trajectories of 48-hour duration are plotted.

Figure 7. (a.) IMPROVE mass closure is plotted in the upper panel. Reconstructed PM2.5 mass is shown as the sum of OM, fine soil, (NH4)2SO4, NH4NO3, and EC. The IMPROVE measured PM2.5 mass is also shown by the black dots. (b.) The contribution to light extinction from each species plotted in (a.) as calculated using the second IMPROVE algorithm for light extinction. Site specific Rayleigh Scattering (10 Mm−1) and NO2 light extinction are also included. NO2 (g) light extinction is calculated using the second IMPROVE algorithm and a daily mean of the measured NO2.

Figure 7. (a.) IMPROVE mass closure is plotted in the upper panel. Reconstructed PM2.5 mass is shown as the sum of OM, fine soil, (NH4)2SO4, NH4NO3, and EC. The IMPROVE measured PM2.5 mass is also shown by the black dots. (b.) The contribution to light extinction from each species plotted in (a.) as calculated using the second IMPROVE algorithm for light extinction. Site specific Rayleigh Scattering (10 Mm−1) and NO2 light extinction are also included. NO2 (g) light extinction is calculated using the second IMPROVE algorithm and a daily mean of the measured NO2.

Figure 8. The second IMPROVE equation was applied to higher time resolution composition data sets (PILS-IC, PILS-TOC, and Aethalometer BC) to estimate hourly light extinction values. These values include the Rayleigh background contribution as well as extinction from NO2. Fractional light extinction contributions, excluding Rayleigh background, from different species are shown for three data points (shown as blue triangles) as examples. Extinction peaks highlighting maximum contributions are represented by points a and c (30 July 2019 09:00 MST and 28 August 2019 12:00 MST). Point b represents a low light extinction period on 18 August 2019 12:00 MST.

Figure 8. The second IMPROVE equation was applied to higher time resolution composition data sets (PILS-IC, PILS-TOC, and Aethalometer BC) to estimate hourly light extinction values. These values include the Rayleigh background contribution as well as extinction from NO2. Fractional light extinction contributions, excluding Rayleigh background, from different species are shown for three data points (shown as blue triangles) as examples. Extinction peaks highlighting maximum contributions are represented by points a and c (30 July 2019 09:00 MST and 28 August 2019 12:00 MST). Point b represents a low light extinction period on 18 August 2019 12:00 MST.

Figure 9. Incremental probability distributions are plotted for 7 species measured directly during the field campaign: SO42-, NH4+, WSOC, Ca2+, NO3, BC, and CH4 and 1 calculated value: particle light extinction (bext). The IPs are based on 24-hr back trajectories and are distance normalized. This is done by multiplying each grid cell value by the distance from the origin site, shown as a black dot. The grid cell resolution is 0.25 degrees square.

Figure 9. Incremental probability distributions are plotted for 7 species measured directly during the field campaign: SO42-, NH4+, WSOC, Ca2+, NO3−, BC, and CH4 and 1 calculated value: particle light extinction (bext). The IPs are based on 24-hr back trajectories and are distance normalized. This is done by multiplying each grid cell value by the distance from the origin site, shown as a black dot. The grid cell resolution is 0.25 degrees square.
Supplemental material

Supplementary_Materials_Naimie_et_al_CAVE.docx

Download MS Word (1.4 MB)

Data availability statement

The data that support the findings of this study are available from the corresponding author, Jeffrey L. Collett, upon reasonable request. https://doi.org/10.25675/10217/235481.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.