Next Article in Journal
PM2.5 Air Pollution Prediction through Deep Learning Using Multisource Meteorological, Wildfire, and Heat Data
Previous Article in Journal
Interannual Variability of Summer Hotness in China: Synergistic Effect of Frequency and Intensity of High Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Particulate Matter and Ammonia Pollution in the Animal Agricultural-Producing Regions of North Carolina: Integrated Ground-Based Measurements and Satellite Analysis

Department of Marine, Earth, and Atmospheric Sciences, North Carolina State University, Raleigh, NC 27606, USA
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(5), 821; https://doi.org/10.3390/atmos13050821
Submission received: 22 March 2022 / Revised: 3 May 2022 / Accepted: 11 May 2022 / Published: 17 May 2022
(This article belongs to the Special Issue Improving Air Quality Predictions and Assessment across Scales)

Abstract

:
Intensive animal agriculture is an important part of the US and North Carolina’s (NC’s) economy. Large emissions of ammonia (NH3) gas emanate from the handling of animal wastes at these operations contributing to the formation of fine particulate matter (PM2.5) around the state causing a variety of human health and environmental effects. The objective of this research is to provide the relationship between ammonia, aerosol optical depth and meteorology and its effect on PM2.5 concentrations using satellite observations (column ammonia and aerosol optical depth (AOD)) and ground-based meteorological observations. An observational-based multiple linear regression model was derived to predict ground-level PM2.5 during the summer months (JJA) from 2008–2017 in New Hanover County, Catawba County and Sampson County. A combination of the Cumberland and Johnston County models for the summer was chosen and validated for Duplin County, NC, then used to predict Sampson County, NC, PM2.5 concentrations. The model predicted a total of six 24 h exceedances over the nine-year period. This indicates that there are rural areas of the state that may have air quality issues that are not captured for a lack of measurements. Moreover, PM2.5 chemical composition analysis suggests that ammonium is a major component of the PM2.5 aerosol.

1. Introduction

Animal feeding operations (AFOs) are agricultural facilities where animals are kept and raised in a small area. While animal agriculture is a large driver of the economy, these farms also create a strain on the environment. The eastern region of North Carolina (NC) has become highly populated with thousands of AFO farms, primarily housing poultry and hogs [1]. NC is ranked 4th in the United States for hog production and 1st in the United States for all poultry production due to the strong turkey production [2]. During 2016, Duplin and Sampson Counties, which are in Eastern NC, contained over 4 million hogs and 100 million chickens [3].
Ammonia (NH3) emissions from AFOs account for more than half of the reactive nitrogen released into the environment [4,5]. The largest emission source of NH3 into the atmosphere is from agricultural sources [6,7]. NH3 plays a significant role in the biogeochemical nitrogen cycle and causes health and environmental effects [8,9,10,11,12,13,14,15,16,17]. The manure that is produced by the animals is frequently used as fertilizer and spread on crops throughout the farm, further releasing the pollutant into the atmosphere. NH3 also has several human health effects including eye, nose and throat irritation, dizziness and headaches [18,19]. Once NH3 is released into the atmosphere it can then deposit by either rainfall or dry deposition to regional waterbodies or it can react with other compounds in the area to create other pollutants and cause further harm [14,15,16,17,20,21,22]. Gaseous NH3 in the atmosphere contributes to the formation of airborne fine particulate matter (PM2.5) through reactions with water vapor and other air pollutants, including oxidation products of sulfur dioxide (SO2) or nitrogen oxides (NO and NO2, or NOX) [23]. Ammonium compounds, including ammonium sulfates (NH4HSO4 and (NH4)2SO4) and ammonium nitrate (NH4NO3), make up a large fraction of PM2.5 [24] (defined as particulate matter with an aerodynamic diameter of 2.5 microns or smaller). Elevated concentrations of PM2.5 are associated with respiratory issues, heart problems and has been linked to premature death [25,26,27]. Ammonia can be important in the nucleation of new particles [28,29]. These ammoniated particles scatter light, attenuating visibility, and can result in some atmospheric cooling [30]. Particulate matter emissions are released through primary sources (e.g., power plants, transportation) [31] as well as secondary sources (e.g., gas to particle conversion). However, in the United States, regulatory strategies to reduce PM2.5 have focused on primary emissions of PM2.5 and on reductions of SO2 and NOX emissions; however, the control of NH3 emissions can also be effective for reducing concentrations of PM2.5.
The North Carolina Clean Smokestacks Act has reduced sulfur dioxide emissions by 89% compared to 1998 emissions [32] and the improved regulations on vehicle manufacturers has helped to reduce NOx emissions. Because these reductions have already been made, any additional effort to reduce both SO2 and NOx emissions would become less cost effective. Thus, some studies have started to indicate that a reduction in ammonia emissions would be a more cost-effective solution in terms of long-term reductions in particulate matter concentrations [33,34,35,36,37,38,39,40]. However, ammonia is not a criteria pollutant. Moreover, any reduction in particulate matter must be measured at monitoring stations, leading to economic, political and logistical problems.
Recent studies in Europe using detailed chemical transport models and time resolved NH3 emissions illustrate the strong nonlinearity between PM2.5 and NH3 emissions, and the reduction in NH3 emissions significantly reduces PM2.5 levels in both summer and winter periods [33]. The impact of NH3 emissions on PM2.5 depends on meteorological parameters (e.g., temperature, relative humidity), the magnitude of the perturbation to NH3 emissions and the abundance of particulate nitrate (NO3), gaseous nitric acid (HNO3) and particulate sulfate (SO4= and HSO4), which are the products of the oxidation of SO2 and NOx, two byproducts of combustion [31,33,35]. Utilizing this foundational knowledge of key processes and an integrated approach of using satellite measurements (ammonia and aerosol optical depth) and ground-based measurements (meteorology and PM2.5), an observational-based statistical model is developed to predict PM2.5 concentrations in these rural regions in an effort to advance Earth system predictability. The abundance of NO3T and SO4T was not measured.
Measurements of PM2.5 in rural regions (especially in intensively agricultural locations) of the US are in general not conducted and are therefore not available. Recent studies [21,41] have shown that increased ammonia from agriculture leads to increased PM2.5. Therefore, the purpose of this work is to estimate PM2.5 concentrations in ammonia-rich environments of Eastern NC using a combination of satellite and in situ data. Using satellite data will allow us to develop a method for predicting ground-level PM2.5 concentrations in areas of high agricultural influence which normally do not have a ground-based measurement site. By utilizing satellite-derived ammonia concentrations and aerosol optical depth (AOD) along with meteorological parameters, we can predict, with reasonable certainty, PM2.5 concentrations across NC. Moreover, a thorough review of the scientific literature did not provide any studies offering predictions of PM2.5 in rural regions based on our process-based approach (i.e., coupling of satellite and ground-based measurements). However, a complex numerical air quality model such as the Weather Research and Forecasting-Community Multiscale Air Quality (WRF-CMAQ) or WRF-Chem Modeling System may also be utilized to predict PM2.5 and NH3 emission from agricultural land using these air quality models across various regions of the world [42,43,44,45].

2. Data and Methodology

2.1. Location Description

AFOs in Eastern NC, while good for the economy, create a strain on the environment. For example, many residences of these high agriculturally active counties have started filling nuisance lawsuits against large hog farms’ odor issues and a general decrease in quality of life for those in the area. For these reasons, our study focuses on various areas in NC, all of which, due to proximity and wind direction, will have a direct influence from these farms. This ultimately will allow for further understanding of how these agricultural areas affect ground-level PM2.5 concentrations.
The most recent measurements within the region of our study are for Kinston, NC, in 2007, as part of the U.S. EPA Chemical Speciation Network for PM2.5. Ammonium compounds formed a substantial portion of the mass of PM2.5 even during the summer months [46], as shown in Table 1 and Figure 1.
Ammonia is directly measured in North Carolina by passive monitors as part of the Ammonia Monitoring Network (AMoN); however, the locations of these monitors is sparse, with only 4 active sites across the state of North Carolina. Moreover, particulate matter is measured at various locations throughout the state with monitors placed and maintained by the Environmental Protection Agency (EPA) and the North Carolina Department of Environmental Quality (NCDEQ). The location of these monitors (Figure 2) is widespread throughout the state and lacks the ability to monitor at a high spatial resolution. To achieve the spatial resolution needed to accurately monitor PM2.5 reductions throughout the state satellite data may be utilized. Satellites allow us to obtain both the spatial and temporal data resolution needed to accurately determine if there is a relationship between ammonia emissions from agricultural farms in Eastern North Carolina and PM2.5 concentrations. Satellite data can also be used to predict PM2.5 concentrations in areas with a high agricultural trace gas emission.
Satellite data are a relatively new data source in air quality monitoring [47,48,49,50]. These satellites can be extremely useful in capturing a relatively high-resolution (such as 0.25° by 0.25° or 1 km by 1 km) air quality image over the entire world for a single day or multiple days. For this reason, satellite data can be especially useful in rural areas around the state, where monitoring networks are scarce. Using satellite data will allow us to develop a method for predicting ground-level PM2.5 concentrations in areas of high agricultural influence which normally do not have a ground-based measurement site. By utilizing satellite-derived ammonia concentrations and aerosol optical depth (AOD) along with meteorological parameters, we can predict, with reasonable certainty, PM2.5 concentrations across North Carolina.

2.2. Satellite-Derived Ammonia Data

Satellite-derived ammonia retrievals from the Infrared Atmospheric Sounding Interferometer (IASI, version 2.1, EUMETSAT) were used with <25% cloud cover for 2008 through to 2017. IASI was used because it provides a longer-term repository of measurement data than other satellite platforms, such as the Cross-track Infrared Sounder (CrIS). IASI is a collaboration between the European Organization for the Exploitation of Meteorological Satellites and the National Centre for Space Studies (CNES) or the French government space agency. IASI was launched in 2006 on the MetOp-A satellite, in 2012 on the MetOp-B satellite and finally on the MetOp-C satellite launched in 2018. They are polar orbiting sun synchronous satellites with an orbit altitude at around 817 km above the Earth’s surface. IASI has an orbital period of about 90 min, crossing the equator at around 9:30 a.m. and p.m. local time, with each subsequent pass displaced by about 22.5 degrees of longitude. The satellite measures detailed infrared spectra over a broad angular swath with a spatial resolution of about 12 km [51]. IASI researchers retrieve estimated total atmospheric column loadings of NH3 and other pollutants based on patterns of infrared absorption. Specific information on the development and comparison testing of the algorithm used to derive ammonia measurements is seen in both [51,52].
At the latitude of NC, morning and evening measurements are collected at about 10 a.m. and 10 p.m. local time. For this study, we used the morning pass to perform the analysis because NH3 concentrations are typically higher during this time. This is because during the morning hour, the mixing layer height is lower in the atmosphere, which allows for NH3 concentrations to remain closer to the surface. The data were filtered by the error estimate reported with the data. Only the points that had an error less than 100% were used in this study.

2.3. Satellite-Derived Aerosol Optical Depth

Satellite-based 550 nm AOD retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) from 2008–2017 were used in this study. MODIS is onboard National Aeronautics and Space Administration’s (NASA) Aqua and Terra satellites which, similar to IASI, are polar orbiting satellites. Aqua and Terra orbit at around 705 km above the Earth’s surface and pass NC once per day. For this study, MODIS level 2 collection 6.1 high confidence AOD at 10 km × 10 km measurements from Aqua were used. Retrievals onboard the Aqua satellite were used specifically due to deterioration issues with the MODIS sensor on the Terra satellite. A visual example of these data can be seen in Figure 3. The blank location in the image can be attributed to cloud cover over the area [51]. Another reason for the blank locations in the image is the filtering process of the data. The data were filtered to only include the best quality data, which NASA indicates is a quality assurance flag of 3.

2.4. Ground-Based PM2.5

The ground-based PM2.5 data were obtained from the Environmental Protection Agency’s (EPA) Outdoor Air Data website, developed based on EPA’s Air Quality System (AQS) network (https://www.epa.gov/outdoor-air-quality-data/download-daily-data; accessed on 30 October 2019). The air quality stations around the state are a collaboration between the EPA office in Research Triangle Park and NCDEQ. Two NC sites were chosen (i.e., training data), one in Cumberland County (35.04 N, 78.95 W) and another in Johnston County (35.59 N, 78.46 W) as input data to develop the statistical model. These locations were chosen based on their proximity to high agricultural areas in NC. The PM2.5 data were available at each of those sites for 3-day intervals for 2008–2017. PM2.5 data were then collected (i.e., test data) from a site in Duplin County (34.95 N, 77.96 W), New Hanover County (34.36 N, 77.84 W6) and Catawba County (35.73 N, 81.36 W) to validate the model performance. The monitor used in NC is an active monitor that takes in a predetermined amount of air and pulls it through a filter within the sensor. After 24 h, the filter is then taken out of the sensor and gravimetrically measured and the mass concentration in μg/m3 is reported [52]. PM2.5 was the only particulate matter size incorporated into this study because PM10 only makes up around less than 20% of the measured PM concentrations in these locations.

2.5. Meteorology Data

Ground-based meteorology data were obtained from the NC State Climate Office for the entire study period from sites in Cumberland, Johnston, Duplin, Sampson, Catawba and New Hanover Counties. Daily averaged temperature in degrees Celsius (°C), pressure in millibar (mb), wind speed in meters per second (m/s) and relative humidity in percent (%) where available from each station and used in the development of the model. Days with reported rainfall were excluded from the final analysis because they would lead to washout events for both the ammonia and AOD data. Figure 4 shows the location of these sites, which for Cumberland, Johnston, Duplin, Catawba and New Hanover are collocated with the ground-level PM2.5 sites. The figure also shows the counties with the highest agricultural production as of 2017 and the predominant wind directions to illustrate the agricultural influence on each of these areas are shown in Figure 5.

2.6. Methodology

IASI ammonia retrievals are filtered so that only those retrievals with a relative error less than 100% are used. AOD retrievals are filtered to use only those observations with the best quality rating (quality assurance flag = 3). For the Cumberland County training case, a total of 135 observation days were identified for which NH3 retrievals and AOD retrievals meeting the above filtration criteria could be compared with available ground-level PM2.5 measurements. For the Johnston County training case, a total of 142 observation days were identified. The quality assurance flag is given by the algorithm team as their assessment of the data quality [53]. These retrievals are then collocated with the ground-level meteorology sites and the ground-level PM2.5 sites. This is carried out through the “average” approach [54,55]. In the averaging method, a 100 km radius of the PM2.5 ground site and the ground-level meteorology station are averaged together to receive a single value for the time period, which in this study is one day. This process is repeated every day for the months of June, July and August for the 9-year period. June, July and August were chosen because they would allow for the highest ammonia retrieval as, statistically, the summer months have the highest ammonia concentrations in NC. Days with reported rainfall are excluded from the analysis because this would result in a washout event for both the PM2.5 data and the ammonia data as they are typically scavenged by rainfall [41].
At the sites in Cumberland and Johnston Counties, the data were then compiled for each day using averaged meteorology parameters and subsequently averaged by month. A multivariate regression model was run in the Statistical Analysis System (SAS) using the IASI ammonia data, the MODIS AOD data and the meteorology data per month at both locations. We then tested the accuracy (i.e., validated) of our model using data from Duplin County. For the Duplin County test, a total of 175 observation days were identified, for which NH3 and AOD retrievals met the filtration criteria, and PM2.5 measurements were available for comparison. Once the model was tested, PM2.5 values were then predicted in Catawba County, New Hanover County and Sampson County to see the influence of agriculture on PM2.5 concentrations across the state and in different eco-regions. For Catawba County, 222 observation days were identified for which NH3 and AOD retrievals met the filtration criteria, 168 valid observation days were identified for New Hanover County and 490 observation days for Sampson County.

3. Results

Multiple linear regression is a classical and well-known statistical technique to quantify the association of a variable (called the dependent variable) with several other variables (called independent variables). ANOVA or analysis of variance is a statistical technique to gauge the statistical significance and practical usefulness of the linear regression towards explaining the dependent variable. In this analysis, we implemented multiple linear regression and ANOVA to quantify the association of PM2.5 (the dependent variable) with six independent variables: ammonia, AOD, T, P, WS and RH. Furthermore, we constructed the ANOVA table to study the statistical significance of this multiple linear regression model. Multiple linear regression models were created to assess the ability of remotely sensed data to predict ground-based PM2.5 concentrations in rural areas that are usually characterized as having high agricultural activity and no ground-based PM2.5 monitor. The models were created for Cumberland and Johnston Counties in NC (i.e., training data) and then validated from Duplin County NC (i.e., test data) against the ground-based measurements located in the county. Cumberland and Johnston Counties were chosen specifically because of their proximity to Sampson County, which has the highest agricultural activity in that state next to Duplin County, which was chosen as the validation site for this reason. Sampson County could not be used for validation because there is no ground-based PM2.5 monitor in the county. After reviewing and testing different model inputs (e.g., a variety of meteorological parameters with data from each county of interest), a combination of Cumberland County and Johnston County data for the entire summer period were chosen:
P M 2.5 = E X P ( 15.14 + ( 1.05 × 10 17 A m m o n i a ) + ( 1.51 A O D ) + ( 0.26 T ) + ( 0.013 P ) + ( 0.040 W S ) + ( 0.013 R H ) )
where PM2.5 is given in units of µg m−3, Ammonia is the total atmospheric column loading of ammonia as retrieved from IASI observations in molecules cm−2, AOD is the MODIS aerosol optical depth, T is temperature in °C, P is pressure in millibars, WS is wind speed in meters per second and RH is relative humidity in percent.
This model gave an r2 value of 0.43. The specific statistics can be seen in Table 2 (an analysis of variance or ANOVA table). Here, DF stands for “degrees of freedom”, i.e., the number of values that are free to vary as we estimate the parameters of the model. The total sum of squares (TSS) is the sum of the squared difference between each PM2.5 observation and the sample average, which is decomposed into two parts: the model sum of squares which is the part of the TSS that can be explained by the regression model, and the error sum of squares which is the remaining part of the TSS which cannot be explained by the model. The mean squares are the sum of squares divided by the respective degrees of freedom. The F value is the ratio of the model mean square to the error mean square, which is a classical statistical metric to quantify the explanatory power of the regression model and the final column is the p-value associated with the F value which indicates the statistical significance of the model (lower means more significant). The F value and the p-value seen in the table indicate that our model was able to explain a statistically significant portion of the data variation [54,55]. Table 3 is a parameter estimate table. The table indicates the significance of a parameter based on the t statistic used in the model. For our model, the table shows that AOD is the most significant variable in our model. This is to be expected, since AOD results from the presence of PM2.5 in the atmosphere. However, AOD is not directly proportional to the mass concentration of PM2.5, but is also dependent on the particle size distribution and other properties of the aerosol and the atmosphere. Relative humidity, temperature and ammonia are also significant variables in our model. The intercept is also relatively significant while pressure and wind speed are the least significant. We know this because Pr > F values have a significant range to them. Anything less than 0.01 has strong significance while anything between 0.01 to 0.05 has appreciable significances.
The model is dominated by data in the range of ~5 to ~25 micrograms per cubic meter. The model Pr < 0.0001 and F Value of ~39 (Table 2), coupled with r2 of ~0.5 (Figure 6, Duplin County), suggest that the model is offering a good prediction overall. Moreover, having Pr for each of the variable (Table 3) less than 0.05 suggests the robust contribution of each variable.
The results of the Duplin County prediction can be seen in Figure 6. The figure shows four different scatter plots each representing a summer month (June, July and August) and a total graph that included all three months in one plot. These plots act to illustrate a visual representation of the model performance at this location. Table 4 illustrates the normalized mean bias (NMB) and normalized mean error (NME) values, which are commonly used to assess the performance of models [54,55].
Once this model was verified for Duplin, it was then used to predict PM2.5 concentrations in New Hanover and Catawba Counties. These locations are also affected by high agricultural production areas; however, they have very different meteorological conditions and processes affecting them daily. The results of this analysis can be seen in Figure 7 and Figure 8. Figure 7 shows the results for New Hanover County and illustrates that despite the different meteorological mechanisms, the model can predict PM2.5 concentrations for this area at a relatively high degree of accuracy. Table 5 shows the model prediction parameters used for this study and clearly shows the model’s prediction struggle in July. This particularly high value is the result of the model’s dependence on AOD values. A high AOD value will sometimes result in an overpredicted PM2.5 value.
Similarly, Figure 8 shows the results for Catawba County; however, the model consistently under predicts concentrations at this location compared to both New Hanover and Duplin Counties, which can be seen in Table 6. This indicates that the different meteorological mechanisms that dominate in Catawba County are not well understood by the model. Catawba County is also at a different elevation than the other two locations, which could explain some of the inconsistent predictions. Overall, these results suggest that the multiple regression model can predict (at a relatively high certainty) for the eastern portion of NC and loses some capabilities in the western portion of the state likely due to topography.
We then used the combination model to predict PM2.5 values in Sampson County. As previously mentioned, Sampson County has one of the highest concentrations of agricultural activity (both animal and crop) in the state. Due to its rural landscape, there is no PM2.5 monitor available and thus no PM2.5 data for this location. Our model results showed, also seen in Figure 9, that the PM2.5 concentrations are low; however, the model did indicate that six days out of the nine-year period were over the EPA NAAQS of 35 μg m−3. In order to further investigate how many exceedances were predicted by the model, the normalized mean error values calculated for New Hanover County were used to see how the model’s error would affect the number of exceedances calculated. New Hanover County errors were used because they were the highest of the three errors calculated for the different counties. Given the errors in the model, the number of exceedances could range from ten to three during the ten-year period. This could indicate that areas near agricultural activity could see higher PM2.5 values, most of which are not being captured due to monitoring locations. This could also suggest that a reduction in ammonia emissions could have a positive impact on PM2.5 concentrations in these high agricultural areas.

4. Conclusions

We have created a multiple regression model that can predict PM2.5 mass concentrations in counties with the highest agricultural activity in NC. High agricultural activity results in high concentration of ammonia being released from agricultural farms. Ammonia is also a precursor to PM2.5 formation and may cause an increase in PM2.5 concentrations in these areas which could lead to poor air quality and quality of life for those living in the area.
This model was developed using remotely sensed data (column ammonia and AOD) and ground-based meteorology, making it resistant to issues in the scarcity of ground-based PM2.5 monitors around the state. Many of the areas around these agricultural farms are rural and have no monitor in place to track these concentrations daily. Satellite data, however, can introduce other limitations. For example, the instrument cannot directly measure ammonia concentrations at the surface, the concentrations are calculated through an algorithm developed by researchers which introduces limitation on the accuracy of the concentrations from the satellite. The same issues can be said for the AOD data. Despite the limitations, satellite algorithms are improving over time and as they improve, utilizing these data will become more important. Satellite data can also be utilized for the meteorology data if ground-based data are not available. Modern-era retrospective analysis for research and applications (MERRA) data is globally modeled meteorological data based on GEOS-5 atmospheric data assimilation system. The data models ground-based temperature, pressure and wind speed well. However, the data do not contain relative humidity, thus it must be obtained experimentally or mathematically. Moreover, PM2.5 chemical composition analysis (Figure 1) suggests that ammonium is a major component of the PM2.5 aerosol.
The model was developed (i.e., training data) for the Cumberland and Johnston Counties. The model was validated (i.e., tested) for Duplin County and applied to predict PM2.5 in New Hanover, Catawba and Sampson Counties, all of which have a predominate influence from these agricultural farms. The model predicted PM2.5 values in Sampson County that were above the National Ambient Air Quality Standard (NAAQS) limit of 35 μg m−3. The need to investigate an ammonia reduction strategy due to its effects on PM2.5 concentrations in high agricultural areas is becoming more prevalent. If the PM2.5 reduction strategies seen in the past, such as the Clean Smokestacks Act and the Clean Air Act, have been as successful as they can be, then reducing ammonia emissions will not only provide air quality improvements, but also a reduction in PM2.5-related issues.

Author Contributions

Conceptualization, V.P.A.; Formal analysis, R.W., W.H.B. and C.B.M.; Funding acquisition, V.P.A.; Investigation, R.W., W.H.B., C.B.M. and V.P.A.; Methodology, R.W., W.H.B. and V.P.A.; Project administration, W.H.B. and V.P.A.; Supervision, V.P.A.; Writing —original draft, R.W.; Writing —review & editing, R.W., W.H.B., C.B.M. and V.P.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Acknowledgments

We would like to first acknowledge the scientist in the Atmospheric Spectroscopy group at the Université Libre de Bruxelles, Belgium (ULB-LATMOS) for access to the IASI retrievals and advice on using the data. Next, we would like to acknowledge various people at the NC Department of Environmental Quality, Elliot Tardif, Bradley McLamb, Nicolas Witcraft and Joette Steger for their advice and guidance on the PM2.5 monitoring data. We thank David Dickey, William Neal Reynolds Distinguished Professor and Srijan Sengupta, Department of Statistics, NC State University, Raleigh, NC 27695, for reviewing and discussing the statistical development in the paper. We thank Swarnali Sanyal, University of Illinois, for her assistance in air quality modeling discussions. We would also like to acknowledge the air quality research group at NC State University for their assistance.

Conflicts of Interest

There are no conflict of interest.

References

  1. Aneja, V.P.; Arya, S.P.; Kim, D.S.; Rumsey, I.C.; Arkinson, H.L.; Semunegus, H.; Bajwa, K.S.; Dickey, D.A.; Stefanski, L.A.; Todd, L.; et al. Characterizing ammonia emissions from swine farms in eastern North Carolina: Part 1—Conventional lagoon and spray technology for waste treatment. J. Air Waste Manag. Assoc. 2008, 58, 1130–1144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ross, M. Poultry by the Numbers. Available online: https://sampson.ces.ncsu.edu/2017/11/poultry-by-the-numbers/ (accessed on 7 March 2019).
  3. United States Department of Agriculture National Agricultural Statistics Service, North Carolina Department of Agriculture and Consumer Science. North Carolina Agricultural Statistics. (Publication No. 218); 2017. Available online: https://www.nass.usda.gov/Statistics_by_State/North_Carolina/Publications/Annual_Statistical_Bulletin/AgStat/NCAgStatBook.pdf (accessed on 30 October 2019).
  4. Galloway, J.N.; Aber, J.D.; Erisman, J.W.; Seitzinger, S.P.; Howarth, R.W.; Cowling, E.B.; Cosby, B.J. The nitrogen cascade. Bioscience 2003, 53, 341–356. [Google Scholar] [CrossRef]
  5. Behera, S.N.; Sharma, M.; Aneja, V.P.; Balasubramanian, R. Ammonia in the atmosphere: A review on emission sources, atmospheric chemistry and deposition on terrestrial bodies. Environ. Sci. Pollut. Res. 2013, 20, 8092–8131. [Google Scholar] [CrossRef] [PubMed]
  6. Faulkner, W.B.; Shaw, B.W. Review of ammonia emission factors for United States animal agriculture. Atmos. Environ. 2008, 42, 6567–6574. [Google Scholar] [CrossRef]
  7. Aneja, V.P.; Schlesinger, W.H.; Li, Q.; Nahas, A.; Battye, W.H. Characterization of the Global Sources of Atmospheric Ammonia from Agricultural Soils. J. Geophys. Res. Atmos. 2020, 125, e2019JD031684. [Google Scholar] [CrossRef]
  8. Erisman, J.W.; Bleeker, A.; Galloway, J.; Sutton, M.S. Reduced nitrogen in ecology and the environment. Environ. Pollut. 2007, 150, 140–149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Wing, S.; Wolf, S. Intensive livestock operations, health, and quality of life among eastern North Carolina residents. Environ. Health Perspect. 2000, 108, 233–238. [Google Scholar] [CrossRef]
  10. Aneja, V.P.; Schlesinger, W.H.; Erisman, J.W. Effects of agriculture upon the air quality and climate: Research, policy, and regulations. Environ. Sci. Technol. 2009, 43, 4234–4240. [Google Scholar] [CrossRef] [Green Version]
  11. Sutton, M.A.; Reis, S.; Riddick, S.N.; Dragosits, U.; Nemitz, E.; Theobald, M.R.; Tang, Y.S.; Braban, C.F.; Vieno, M.; Dore, A.J.; et al. Towards a climate-dependent paradigm of ammonia emission and deposition. Philos. Trans. R. Soc. B Biol. Sci. 2013, 368, 20130166. [Google Scholar] [CrossRef]
  12. Battye, W.; Aneja, V.P.; Schlesinger, W.H. Is nitrogen the next carbon? Earth’s Future 2017, 5, 894–904. [Google Scholar] [CrossRef] [Green Version]
  13. Baker, J.; Battye, W.H.; Robarge, W.P.; Arya, S.P.; Aneja, V.P. Characterization and Fate of Ammonia from Poultry Operations: Their Emissions, Transport, and Deposition in the Chesapeake Bay. Sci. Total Environ. 2019, 706, 135290. [Google Scholar] [CrossRef] [PubMed]
  14. Aneja, V.P.; Arya, S.P.; Rumsey, I.C.; Kim, D.S.; Bajwa, K.S.; Williams, C.M. Characterizing ammonia emissions from swine farms in eastern North Carolina: Reduction of emissions from water-holding structures at two candidate superior technologies for waste treatment. Atmos. Environ. 2008, 42, 3291–3300. [Google Scholar] [CrossRef]
  15. Aneja, V.P.; Schlesinger, W.H.; Erisman, J.W. Farming pollution. Nat. Geosci. 2008, 1, 409–411. [Google Scholar] [CrossRef]
  16. Stephen, K.; Aneja, V.P. Trends in agricultural ammonia emissions and ammonium concentrations in precipitation over the Southeast and Midwest United States. Atmos. Environ. 2008, 42, 3238–3252. [Google Scholar] [CrossRef] [Green Version]
  17. Linker, L.C.; Dennis, R.; Shenk, G.W.; Batiuk, R.A.; Grimm, J.; Wang, P. Computing atmospheric nutrient loads to the Chesapeake Bay watershed and tidal waters. JAWRA J. Am. Water Resour. Assoc. 2013, 49, 1025–1041. [Google Scholar] [CrossRef]
  18. Schiffman, S.S. Livestock odors: Implications for human health and well-being. J. Anim. Sci. 1998, 76, 1343–1355. [Google Scholar] [CrossRef] [Green Version]
  19. Wilson, S.M.; Serre, M.L. Examination of atmospheric ammonia levels near hog CAFOs, homes, and schools in Eastern North Carolina. Atmos. Environ. 2007, 41, 4977–4987. [Google Scholar] [CrossRef]
  20. Bittman, S.; Mikkelsen, R. Ammonia emissions from agricultural operations: Livestock. Better Crops 2009, 93, 28–31. [Google Scholar]
  21. Aneja, V.P.; Arya, S.P.; Rumsey, I.C.; Kim, D.S.; Arkinson, H.L.; Semunegus, H.; Bajwa, K.S.; Dickey, D.A.; Stefanski, L.A.; Todd, L.; et al. Characterizing Ammonia Emissions from Swine Farms in Eastern North Carolina: Part 2—Potential Environmentally Superior Technologies for Waste Treatment. J. Air Waste Manag. Assoc. 2008, 58, 1145–1157. [Google Scholar] [CrossRef]
  22. Baker, J.L. Characterization and Fate of Ammonia from Poultry Operations: Their Emissions, Transport, and Deposition in the Chesapeake Bay. Ph.D. Thesis, North Carolina State University, Raleigh, NC, USA, 2020. [Google Scholar]
  23. Li, H.; Xin, H.; Burns, R.T.; Hoff, S.J.; Harmon, J.D.; Jacobson, L.D.; Noll, S.; Koziel, J.A. Ammonia and PM emissions from a tom turkey barn in Iowa. In 2008 Providence, Rhode Island, June 29–July 2, 2008; American Society of Agricultural and Biological Engineers: Saint Joseph, MI, USA, 2008; p. 1. [Google Scholar]
  24. Kwok, R.H.F.; Napelenok, S.L.; Baker, K.R. Implementation and evaluation of PM2.5 source contribution analysis in a photochemical model. Atmos. Environ. 2013, 80, 398–407. [Google Scholar] [CrossRef]
  25. Kampa, M.; Castanas, E. Human health effects of air pollution. Environ. Pollut. 2008, 151, 362–367. [Google Scholar] [CrossRef] [PubMed]
  26. Pui, D.Y.; Chen, S.C.; Zuo, Z. PM2.5 in China: Measurements, sources, visibility and health effects, and mitigation. Particuology 2014, 13, 1–26. [Google Scholar] [CrossRef]
  27. Kravchenko, J.; Rhew, S.H.; Akushevich, I.; Agarwal, P.; Lyerly, H.K. Mortality and health outcomes in North Carolina communities located in close proximity to hog concentrated animal feeding operations. N. C. Med. J. 2018, 79, 278–288. [Google Scholar] [CrossRef] [PubMed]
  28. Holmes, N.S. A review of particle formation events and growth in the atmosphere in the various environments and discussion of mechanistic implications. Atmos. Environ. 2007, 41, 2183–2201. [Google Scholar] [CrossRef] [Green Version]
  29. Herb, J.; Nadykto, A.B.; Yu, F. Large ternary hydrogen-bonded pre-nucleation clusters in the Earth’s atmosphere. Chem. Phys. Lett. 2011, 518, 7–14. [Google Scholar] [CrossRef]
  30. Pinder, W.; Davidson, E.A.; Goodale, C.L.; Greaver, T.L.; Herrick, J.D.; Liu, L. Climate change impacts of US reactive nitrogen. Proc. Natl. Acad. Sci. USA 2012, 109, 7671–7675. [Google Scholar] [CrossRef] [Green Version]
  31. Seinfeld, J.H.; Pandis, S.N. Atmospheric Chemistry and Physics: From Air Pollution to Climate Change; John Wiley & Sons: New York, NY, USA, 2016. [Google Scholar]
  32. North Carolina Department of Environmental and Natural Resources, North Carolina Utilities Commission. Implementation of the “Clean Smokestacks Act” (Report No. XII). 2014. Available online: https://files.nc.gov/ncdeq/Air%20Quality/news/leg/2014_Clean_Smokestacks_Act_Report.pdf (accessed on 30 October 2019).
  33. Makar, P.A.; Moran, M.D.; Zheng, Q.; Cousineau, S.; Sassi, M.; Duhamel, A.; Besner, M.; Davignon, D.; Crevier, L.P.; Bouchet, V.S. Modelling the impacts of ammonia emissions reductions on North American air quality. Atmos. Chem. Phys. 2009, 9, 7183–7212. [Google Scholar] [CrossRef] [Green Version]
  34. Paulot, F.; Jacob, D.J. Hidden Cost of U.S. Agricultural Exports: Particulate Matter from Ammonia Emissions. Environ. Sci. Technol. 2014, 48, 903–908. [Google Scholar] [CrossRef]
  35. Gong, L.; Lewicki, R.; Griffin, R.J.; Tittel, F.K.; Lonsdale, C.R.; Stevens, R.G.; Pierce, J.R.; Malloy, Q.G.; Travis, S.A.; Bobmanuel, L.M.; et al. Role of atmospheric ammonia in particulate matter formation in Houston during summertime. Atmos. Environ. 2013, 77, 893–900. [Google Scholar] [CrossRef]
  36. Stokstad, E. Ammonia pollution from farming may exact hefty health costs. Science 2014, 343, 238. [Google Scholar] [CrossRef]
  37. Wu, Y.; Gu, B.; Erisman, J.W.; Reis, S.; Fang, Y.; Lu, X.; Zhang, X. PM2.5 pollution is substantially affected by ammonia emissions in China. Environ. Pollut. 2016, 218, 86–94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Zhao, Z.Q.; Bai, Z.H.; Winiwarter, W.; Kiesewetter, G.; Heyes, C.; Ma, L. Mitigating ammonia emission from agriculture reduces PM2.5 pollution in the Hai River Basin in China. Sci. Total Environ. 2017, 609, 1152–1160. [Google Scholar] [CrossRef] [PubMed]
  39. Guo, H.; Otjes, R.; Schlag, P.; Kiendler-Scharr, A.; Nenes, A.; Weber, R.J. Effectiveness of ammonia reduction on control of fine particle nitrate. Atmos. Chem. Phys. 2018, 18, 12241–12256. [Google Scholar] [CrossRef] [Green Version]
  40. Cheng, T.; Chen, H.; Gu, X.; Yu, T.; Guo, J.; Guo, H. The inter-comparison of MODIS, MISR and GOCART aerosol products against AERONET data over China. J. Quant. Spectrosc. Radiat. Transf. 2012, 113, 2135–2145. [Google Scholar] [CrossRef]
  41. Boylan, J.W.; Russell, A.G. PM and light extinction model performance metrics, goals, and criteria for three-dimensional air quality models. Atmos. Environ. 2006, 40, 4946–4959. [Google Scholar] [CrossRef]
  42. Mazzeo, A.; Zhong, J.; Hood, C.; Smith, S.; Stocker, J.; Cai, X.; Bloss, W.J. Modelling the impact of national vs. local emission reduction on PM2.5 in the West Midlands, UK using WRF-CMAQ. Atmosphere 2022, 13, 377. [Google Scholar] [CrossRef]
  43. Kim, E.; Kim, B.U.; Kim, H.C.; Kim, S. Sensitivity of fine particulate matter concentrations in South Korea to regional ammonia emissions in Northeast Asia. Environ. Pollut. 2021, 273, 116428. [Google Scholar] [CrossRef]
  44. Liu, Z.; Zhou, M.; Chen, Y.; Pan, Y.; Song, T.; Ji, D.; Chen, Q.; Zhang, L. The nonlinear response of fine particulate matter pollution to ammonia emission reductions in North China. Environ. Res. Lett. 2021, 16, 034014. [Google Scholar] [CrossRef]
  45. Choi, H.; Sunwoo, Y. Environmental Benefits of Ammonia Reduction in an Agriculture-Dominated Area in South Korea. Atmosphere 2022, 13, 384. [Google Scholar] [CrossRef]
  46. Megaritis, A.; Fountoukis, C.; Charalampidis, P.; Pilinis, C.; Pandis, S. Response of fine particulate matter concentrations to changes of emissions and temperature in Europe. Atmos. Chem. Phys. 2013, 13, 3423–3443. [Google Scholar] [CrossRef] [Green Version]
  47. Battye, W.H.; Bray, C.D.; Aneja, V.P.; Tong, D.; Lee, P.; Tang, Y. Evaluating ammonia (NH3) predictions in the NOAA National Air Quality Forecast Capability (NAQFC) using in situ aircraft, ground-level, and satellite measurements from the DISCOVER-AQ Colorado campaign. Atmos. Environ. 2016, 140, 342–351. [Google Scholar] [CrossRef]
  48. Bray, C.D.; Battye, W.; Aneja, V.P.; Tong, D.; Lee, P.; Tang, Y.; Nowak, J.B. Evaluating ammonia (NH3) predictions in the NOAA National Air Quality Forecast Capability (NAQFC) using in-situ aircraft and satellite measurements from the CalNex. Atmos. Environ. 2010, 163, 65–76. [Google Scholar] [CrossRef]
  49. Aneja, V.P.; Pillai, P.R.; Isherwood, A.; Morgan, P.; Aneja, S.P. Particulate matter pollution in the coal-producing regions of the Appalachian Mountains: Integrated ground-based measurements and satellite analysis. J. Air Waste Manag. Assoc. 2017, 67, 421–430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. House Bill 515: An Act to Enact the Clean Water Responsibility and Environmentally Sound Policy Act, a Comprehensive and Balanced Program to Protect Water Quality, Public Health, and the Environment (HB 515); 1997 Session; General Assembly of North Carolina: Raleigh, NC, USA, 1997.
  51. Van Damme, M.; Wichink Kruit, R.J.; Schaap, M.; Clarisse, L.; Clerbaux, C.; Coheur, P.F.; Dammers, E.; Dolman, A.J.; Erisman, J.W. Evaluating 4 years of atmospheric ammonia (NH3) over Europe using IASI satellite observations and LOTOS-EUROS model results. J. Geophys. Res. Atmos. 2014, 119, 9549–9566. [Google Scholar] [CrossRef] [Green Version]
  52. Levy, R.C.; Remer, L.A.; Kleidman, R.G.; Mattoo, S.; Ichoku, C.; Kahn, R.; Eck, T.F. Global evaluation of the Collection 5 MODIS dark-target aerosol products over land. Atmos. Chem. Phys. 2010, 10, 10399. [Google Scholar] [CrossRef] [Green Version]
  53. Shimshock, J.P.; De Pena, R.G. Below-cloud scavenging of tropospheric ammonia. Tellus B Chem. Phys. Meteorol. 1989, 41, 296–304. [Google Scholar] [CrossRef]
  54. Devore, J.L. Probability and Statistics for Engineering and the Sciences; Brooks/Cole: Belmont, CA, USA, 2011. [Google Scholar]
  55. Teetor, P. R Cookbook: Proven Recipes for Data Analysis, Statistics, and Graphics; O’Reilly Media, Inc.: Newton, MA, USA, 2011. [Google Scholar]
Figure 1. Speciation of PM2.5 (µg m−3) during the summer months of 2007 in Kinston, NC, USA.
Figure 1. Speciation of PM2.5 (µg m−3) during the summer months of 2007 in Kinston, NC, USA.
Atmosphere 13 00821 g001
Figure 2. PM2.5 site locations. These sites are set up and maintained by the Environmental Protection Agency (EPA) and the North Carolina Department of Environmental Quality (NCDEQ).
Figure 2. PM2.5 site locations. These sites are set up and maintained by the Environmental Protection Agency (EPA) and the North Carolina Department of Environmental Quality (NCDEQ).
Atmosphere 13 00821 g002
Figure 3. Aerosol optical depth (AOD) retrievals from the MODIS Aqua satellite for 30 July 2011.
Figure 3. Aerosol optical depth (AOD) retrievals from the MODIS Aqua satellite for 30 July 2011.
Atmosphere 13 00821 g003
Figure 4. The initial locations that we built the model from are seen in blue while the locations we modeled can be seen in black. Total ammonia emissions from agricultural farms are indicated by the shading.
Figure 4. The initial locations that we built the model from are seen in blue while the locations we modeled can be seen in black. Total ammonia emissions from agricultural farms are indicated by the shading.
Atmosphere 13 00821 g004
Figure 5. The initial locations that we built the model from are seen in blue while the locations we modeled can be seen in black. Animal density is indicated by the shading and the wind roses are used to indicate the predominate wind direction and speed from three sites, New Hanover, Sampson and Catawba Counties.
Figure 5. The initial locations that we built the model from are seen in blue while the locations we modeled can be seen in black. Animal density is indicated by the shading and the wind roses are used to indicate the predominate wind direction and speed from three sites, New Hanover, Sampson and Catawba Counties.
Atmosphere 13 00821 g005
Figure 6. Scatter plots for Duplin County indicating model performance. The black points indicate the modeled versus predicted point and the grey bars indicate the ±1 standard deviation (s.d.) of the modeled dataset. The red line indicated the one-to-one line and the blue line is the best fit line for the dataset. The slope of the best fit line is indicated in blue in the corner of each figure.
Figure 6. Scatter plots for Duplin County indicating model performance. The black points indicate the modeled versus predicted point and the grey bars indicate the ±1 standard deviation (s.d.) of the modeled dataset. The red line indicated the one-to-one line and the blue line is the best fit line for the dataset. The slope of the best fit line is indicated in blue in the corner of each figure.
Atmosphere 13 00821 g006
Figure 7. Scatter plots for New Hanover County indicating model performance. The black points indicate the modeled versus predicted point, and the grey bars indicate the ±1 standard deviation of the modeled dataset. The red line indicated the one-to-one line and the blue line is the best fit line for the dataset. The slope of the best fit line is indicated in blue in the corner of each figure.
Figure 7. Scatter plots for New Hanover County indicating model performance. The black points indicate the modeled versus predicted point, and the grey bars indicate the ±1 standard deviation of the modeled dataset. The red line indicated the one-to-one line and the blue line is the best fit line for the dataset. The slope of the best fit line is indicated in blue in the corner of each figure.
Atmosphere 13 00821 g007
Figure 8. Scatter plots for Catawba County indicating model performance. The black points indicate the modeled versus predicted point, and the grey bars indicate the ±1 standard deviation of the modeled dataset. The red line indicated the one-to-one line and the blue line is the best fit line for the dataset. The slope of the best fit line is indicated in blue in the corner of each figure.
Figure 8. Scatter plots for Catawba County indicating model performance. The black points indicate the modeled versus predicted point, and the grey bars indicate the ±1 standard deviation of the modeled dataset. The red line indicated the one-to-one line and the blue line is the best fit line for the dataset. The slope of the best fit line is indicated in blue in the corner of each figure.
Atmosphere 13 00821 g008
Figure 9. Histogram plot illustrating the model predictions in Sampson County for June (a), July (b), August (c) and total (d). The blue lines indicate the number of times (y-axis) the PM2.5 concentrations reached a certain amount (x-axis). The red line indicates the NAAQS exceedance value.
Figure 9. Histogram plot illustrating the model predictions in Sampson County for June (a), July (b), August (c) and total (d). The blue lines indicate the number of times (y-axis) the PM2.5 concentrations reached a certain amount (x-axis). The red line indicates the NAAQS exceedance value.
Atmosphere 13 00821 g009
Table 1. PM2.5 speciation in Kinston, NC, during the summer months of 2006 and 2007 *.
Table 1. PM2.5 speciation in Kinston, NC, during the summer months of 2006 and 2007 *.
PM2.5 Composition (µg m−3)
Summer 2006Summer 2007
Ammonium sulfates and ammonium nitrate7.087.34
Organic carbon4.445.10
Elemental carbon0.220.24
Other and unidentified2.282.45
Total14.0315.13
Fraction ammonium sulfates and nitrate50%49%
Fraction EC/OC33%35%
* Data source: PM2.5 Chemical Speciation Network (CSN) for the months of June, July and August. https://www.epa.gov/amtic/chemical-speciation-network-csn, https://www.epa.gov/outdoor-air-quality-data/interactive-map-air-quality-monitors, accessed on 30 October 2019.
Table 2. Analysis of variance (ANOVA) table.
Table 2. Analysis of variance (ANOVA) table.
SourceDFSum of SquaresMean SquareF ValuePr > F
Model618.843.1402738.99<0.0001
Error30924.880.08053
Corrected Total31543.73
Table 3. Parameter estimates for the model.
Table 3. Parameter estimates for the model.
VariableParameter EstimateStandard ErrorType II SSF ValuePr > F
Intercept15.145.070.728.920.0030
Ammonia1.05 × 10−172.99 × 10−180.9812.210.0005
AOD1.510.167.0587.53<0.0001
Temp0.260.0071.0613.140.0003
Pressure−0.0130.0050.546.650.0104
WS−0.040.020.384.750.0301
RH−0.0130.0032.1226.35<0.0001
Table 4. Normalized mean bias (NMB) and normalized mean error values (NME) for Duplin County by month.
Table 4. Normalized mean bias (NMB) and normalized mean error values (NME) for Duplin County by month.
MonthNormalized Mean BiasNormalized Mean Error
June0.47%25.31%
July13.19%24.66%
August14.59%25.23%
Table 5. Normalized mean bias (NMB) and normalized mean error values (NME) for New Hanover County by month.
Table 5. Normalized mean bias (NMB) and normalized mean error values (NME) for New Hanover County by month.
MonthNormalized Mean BiasNormalized Mean Error
June17.62%28.55%
July17.70%43.61%
August19.86%28.04
Table 6. Normalized mean bias (NMB) and normalized mean error values (NME) for Catawba County by month.
Table 6. Normalized mean bias (NMB) and normalized mean error values (NME) for Catawba County by month.
MonthNormalized Mean BiasNormalized Mean Error
June−9.68%23.92%
July−13.19%26.01%
August−14.85%27.10%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wiegand, R.; Battye, W.H.; Myers, C.B.; Aneja, V.P. Particulate Matter and Ammonia Pollution in the Animal Agricultural-Producing Regions of North Carolina: Integrated Ground-Based Measurements and Satellite Analysis. Atmosphere 2022, 13, 821. https://doi.org/10.3390/atmos13050821

AMA Style

Wiegand R, Battye WH, Myers CB, Aneja VP. Particulate Matter and Ammonia Pollution in the Animal Agricultural-Producing Regions of North Carolina: Integrated Ground-Based Measurements and Satellite Analysis. Atmosphere. 2022; 13(5):821. https://doi.org/10.3390/atmos13050821

Chicago/Turabian Style

Wiegand, Rebecca, William H. Battye, Casey Bray Myers, and Viney P. Aneja. 2022. "Particulate Matter and Ammonia Pollution in the Animal Agricultural-Producing Regions of North Carolina: Integrated Ground-Based Measurements and Satellite Analysis" Atmosphere 13, no. 5: 821. https://doi.org/10.3390/atmos13050821

APA Style

Wiegand, R., Battye, W. H., Myers, C. B., & Aneja, V. P. (2022). Particulate Matter and Ammonia Pollution in the Animal Agricultural-Producing Regions of North Carolina: Integrated Ground-Based Measurements and Satellite Analysis. Atmosphere, 13(5), 821. https://doi.org/10.3390/atmos13050821

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop