Fine Particulate Matter Concentrations during Independence Day Fireworks Display in the Lower Rio Grande Valley Region, South Texas, USA

Fireworks are typically discharged as a mark of celebration and joy in many societies spanning various cultures. In the United States of America, 4th July is celebrated as the Independence Day when the nation overthrew the British colonial yoke in 1776. While this day instills a sense of patriotism in every American's heart, it is also a major PM2.5 air pollution concern. This study is first of its type in the Lower Rio Grande Valley (RGV) Region of South Texas, USA, that characterizes fine particulate matter pollution. Using a low-cost sensor (TSI BlueSky Air Quality Monitor), real-time PM2.5 measurements were assessed at eleven different locations in four different towns and cities of Lower RGV Region: Brownsville, Edinburg, Weslaco, and Port Isabel. Hourly PM2.5 concentrations from July 03–06, 2021 are presented in this research work. Intraurban PM2.5 spatial and temporal variations provide an insight on the general population's exposure burden during the festive period. Results indicate an increase in fine particulate matter pollution across the region, but the levels do not exceed the U.S. National Ambient Air Quality Standards (NAAQS). Findings from this study would possibly help in the formulation of effective firework policies to minimize the pollution impact.


Introduction
Fireworks display can be a major source of air pollution and has the potential to cause deleterious health effects [1][2][3][4]. Of all the major criteria air pollutants, fine particulate matter (aerosol particles with an aerodynamic diameter of less than 2.5 μm), i.e., PM 2.5 is measured and studied due to its documented health effects [5][6][7]. PM 2.5 particles can deposit deep in the lungs thereby resulting in damage to the lower thoracic region [8]. PM 2.5 exposures have both acute and chronic impacts on human health [2,3,[9][10][11]. Continued exposure to fine PM results in the development of relevant cardiovascular and or respiratory diseases [1,4,10]. e major sources of fine particulate matter pollution are typically traffic emissions including both tailpipe and break and tire tear and wear, power plant emissions, emissions from smokestacks, and agricultural burning [4,12,13].
A study published in 2015 reported that PM 2.5 concentrations are elevated on 4 th July and continue to remain high until the morning of 5 th July [20]. e authors analyzed PM 2.5 data from 315 US air quality monitory sites and showed a 42% increase (5 µg/m ) of 24 hr PM 2.5 concentrations in the United States during Independence Day [4,20]. During the Chinese Lantern Festival in Beijing, Wang and colleagues showed a six-fold increase in PM 2.5 concentrations on the Lantern Day as compared to normal days [17]. Another study from Beijing, China quantified PM 2.5 average concentrations and documented highest levels during firework days at 248.9 µg·m −3 in the 2015 Spring Festival. [4,21,22]. During the Montreal International Fireworks Competition in Quebec, Canada, researchers showed that PM 2.5 levels can be as high as 1000 µg/ m during the display period of about 45 minutes [23]. In New Delhi, India, during the day of the Diwali (Festival of Lights) festival, PM 2.5 levels were 588 µg/m in 2007, and 389 µg/m in 2008 [13].
Till date, no study has been conducted till date in the Lower RGV Region of South Texas, USA, to assess an increase in fine particulate matter pollution levels during the US Independence Day celebrations. e present study is first of its kind to use TSI BlueSky low-cost sensor to characterize fine particulate matter pollution in this region. Recently, the usage of low-cost air quality sensors to measure PM pollution has been gaining credence [9,24]. ese sensors are easy to operate, very convenient in terms of price, energy, and mobility, and record data in real time and can store data for future downloading purposes [9,25].

Site Selection and Study Period.
e study area is in the Lower RGV Region of South Texas, the USA, with a population of 1,402,512 persons [26]. Eleven BlueSky sensors were deployed in four towns and cities and labeled as follows: five monitors in Brownsville (B1-B5), three monitors in Edinburg (E1-E3), two monitors in Weslaco (W1 and W2), and one monitor in Port Isabel (PI). ese intraurban sensor locations are illustrated in Figure 1 with color-coordinated indicators to differentiate the cities. Brownsville locations are marked with blue indicators, Edinburg has orange, Weslaco has green, and lastly, Port Isabel is purple. Each marker signifies a deployed BlueSky monitor with tactical outdoor locations to portray significant areas such as schools, residential area, and other neighborhood locations.
Sensor B1 in Brownsville is right across from Dean Porter Park and adjacent to the Gladys Porter Zoo in a semiresidential area. B2 is deployed in the University of Texas Rio Grande Valley (UTRGV) Brownsville Campus police department, near the campus student dormitories. B3 is located at the Music and Science Learning Center on the University Boulevard at the UTRGV campus, and B4 is situated adjacent to Texas State Highway 69E only 0.04 miles away from the U.S.-Mexico International Port of Entry. e entire university is located on the main highway 69E on the edge of the US-Mexican border, as seen in Figure 1, so exposure to traffic pollutants from cross-border vehicular traffic is a health concern. B5 is in a neighborhood  Figure 1 also shows the locations of the five Continuous Ambient Monitoring Station (CAMS) sites in the region. ese sites are maintained by the Texas Commission on Environmental Quality (TCEQ) and they monitor air pollutant such as PM 2.5 and meteorological parameters such as temperature, resultant wind speed, and solar radiation. e CAMS sites are shown as black triangular symbol in the study map. e five CAMS sites are: C323 (Port Isabel), C80 (Brownsville), C1023 (Harlingen), C1046 (Edinburg), and C43 (Mission). Out of the five CAMS sites, only three record data for PM 2.5 . ese are CAMS sites C43, C80, and C323. e other two C1046 and C 1023 do not collect PM 2.5 data. All five sites collect various meteorological parameters such as temperature, resultant wind speed, and two sites (C43 and C80) provide data on solar radiation. e UTRGV Edinburg campus has two monitors, E1 and E2 facing commonly used roads (Schunior St. and 107 Texas, respectively). E3 is deployed in a closed gated residential community. In Weslaco, W1 is located further from the city in a residence right off Farm to Market Road 88 (FM88). W2 is deployed at one of the buildings of the Weslaco Police department on the frontage of Texas State Expressway 83. Lastly, the PI sensor is deployed at the UTRGV coastal labs between two neighborhoods near the popular tourist destination of South Padre Island.
Sensor deployment locations were considered such that they are an accurate representation of the daily exposure levels of PM 2.5 in the neighborhood. e RGV Region has previously hosted multiple fireworks in their cities, as well as supporting many fireworks stands, so residents can easily purchase firecrackers for their celebrations. In 2021, the cities and towns of Brownsville, Edinburg, Harlingen, McAllen, Weslaco, and South Padre Island celebrated the Independence Day with much fanfare by grand firework display shows. e study period starts the evening of July 3 rd , 7 : 00 pm through July 6 th at 6 : 59 pm. BlueSky monitors were configured to log PM concentrations every 5 minutes, and these data were collected by a built-in microSD card. ose 5minute concentrations were formulated into hourly data to be further assessed.
e day prior and after 4 th July are considered primarily as control days to help understand the temporal variation in the PM 2.5 concentrations.

Instrumentation. In this study, BlueSky Low-Cost
Sensors were used (Model: 8143 by TSI Incorporated, Minnesota, U.S.). is instrumentation is easy to install and weighs about 0.35 lb. e sensor does not require much power, approximately less than 5 W (5 VDC @ 1 Amp). e PM sensor included is precalibrated similarly to other highquality TSI equipment like the DustTrak ™ models (TSI Incorporated, Minnesota, U.S). Self-diagnostic tests are configured to daily cleaning intervals to attain high-quality data. e PM sensor measures from the range 0 to 1000 μg/ m 3 with the measurement resolution of 1 μg/m 3 and a response time being 1 second [27]. It is prudent to mention here that, albeit the usage of low-cost sensors for air quality monitoring has increased substantially in the last few years, issues such as sensor baseline drift, sensitivity to variations in meteorological parameters such as ambient temperature and relative humidity, instrument measurement artifact needs to be addressed and accounted for as outlined succinctly by Morawska and colleagues [28]. Extensive quality assurance and quality control (QA and QC) studies conducted by the South Coast Air Quality Management District (South Coast AQMD) on BlueSky sensors have shown a good record of performance and evaluations with the sensor showing a moderate to strong PM 2.5 (0.66 < R 2 < 0.78) correlation with other Federal Equivalent Method (FEM) instruments like FEM GRIMM and FEM Teledyne API T640 in the field [29].

Statistical Data Analysis.
Data from these low-cost sensors were inputted in Microsoft Excel (2021) for calculating hourly concentrations and time series. e time series depicts PM 2.5 concentrations through the study period hours, while also comparing data from each site. Spearman's Rho correlations were calculated with SPSS for MacOS (SPSS, Inc., Chicago, IL) to estimate temporal variability in PM 2.5 concentrations across the eleven locations. Visual data analysis from R programming (RStudio, Inc., Boston, MA) was used to demonstrate box plot hourly variability. e boxes are the interquartile ranges (75 th and 25 th ), the whiskers show the minimum and maximum values, the outliers are shown in asterisk. e median is indicated by the black line inside the boxes and the diamond in the boxes indicate the mean.
Spatial variation in PM 2.5 concentrations in each site are analyzed with the performance of Coefficient of Divergence (COD) analysis was performed to study the spatial variation between the various study sites [30][31][32]. e COD provides a degree of uniformity between two simultaneously sampled sites, j and k by the following equation: where the number of observations is indicated by p and x ij is the i th concentration measured at site j over the sampling period. COD values less than 0.20 indicate the two observed sites are spatially homogeneous in terms of the pollutant concentration. COD equal to greater than 0.20 establishes spatial heterogeneity in the pollutant concentrations or significant differences between the two simultaneously sampled sites [31].     In the city of Brownsville, across the five sampling sites one can observe a slight increase in mean PM 2.5 concentrations on the day of Independence. For example, site B1 showed the mean concentration of 13.6 ± 4.9 µg/m for between July 04, 19 : 00 hours-July 05, 18 : 59 hours. For the same time frame, site B3 had the lowest concentration (8.8 ± 2.9) µg/m and site B5 had the highest (14.8 ± 7.1) µg/ m . In the city of Edinburg, hourly concentrations ranged from 11.1 ± 2.7 µg/m at E3 to 14.1 ± 4.1 µg/m at E2 for the same study period. Similarly, the sensor at Port Isabel (PI) showed a mean of 10.8 ± 4.7 µg/m . e PM 2.5 concentrations reverted to control values starting July 05. As is obvious from Tables 1 and 2, the mean concentrations varied from the lowest at W2 (5.3 ± 2.1 µg/m ) to the highest at B1 (9.7 ± 5.1 µg/m ) and E2 (9.7 ± 4.1 µg/m ). It is also important to mention here a striking similarity in the pollutant concentrations a day prior and after the 4 th of July Independence Day celebrations as is obvious from Table 1.

Results and Discussion
Site W1 is located in a residential area as mentioned before and it is quite probable that the huge amount of fireworks lit in this neighborhood by the residents resulted in an hourly maximum of 51.2 µg/m on the Independence Day. For that same period, the maximum recorded at site W2 is only 12.3 µg/m . is makes sense because site W2 is located at one of the buildings of the Weslaco Police Department and there were no fireworks reported in this neighborhood. is same trend is observed in the city of Brownsville. Sites B1 and B5 were in residential neighborhoods where we expect people to light firecrackers. Hence the maximum concentration recorded at these two sites were 20.3 µg/m (B1) and 38.5 µg/m (B5). is is in sharp contrast to site B4 (adjacent to the interstate highway) where the maximum concentration was recorded as 15.5 µg/m , as well as site B3 which is the University campus. ese findings, therefore, suggest the role of fireworks in the high pollutant concentration observed during this festive period. e three CAMS sites also exhibited higher PM 2.5 concentrations between July 04 and July 05. e mean concentration at C43 (Mission) was 19.54 ± 5.96 µg/m , C80 (Brownsville), 16.54 ± 3.82 µg/m , and C1046 (Edinburg), 18.17 ± 5.38 µg/m . e maximum hourly concentration recorded at C43 and C1046 for this period is 30 µg/m . For the next 24-hour period, the PM 2.5 values at the three CAMS sites ranged from 9.67 ± 6.03 µg/m (C80) to 12.08 ± 4.83 µg/ m (C43) suggesting a decrement in the pollutant concentrations after the firework festive period. ese values between July 05 and July 06 (post celebrations) reflect the same low concentration trend as was observed between July 03 and July 04 (pre celebrations, 10.78 ± 6.48 µg/m at C80 to 14.83 ± 9.02 µg/m at C1046).
e United States Environmental Protection Agency is mandated by the United States Congress, under the auspices of the Clean Air Act of 1970 and its amendments thereof in 1990 and onwards, to set National Ambient Air Quality Standards (NAAQS) for the criteria air pollutants-one of which is PM 2.5. Currently, the 24-hour mean average for PM 2.5 is 35 µg/m . Based on our research findings, we can posit that none of the study sites during the three-day study period exceeded the PM 2.5 NAAQS standards. However, it would be prudent to mention here that the daily mean PM 2.5 standard by the World Health Organization is 15 µg/m and our study sites are just within this stricter standard as well during the study period. Figure 3 displays the box plots for the hourly PM 2.5 concentrations across the eleven sites and the three CAMS sites from July 03 to July 06. e color coordination for the box plots is as followed: blue for Brownsville, gray for the CAMS sites, green for Edinburg, gold for Port Isabel, and salmon for Weslaco. Box plots are helpful to detect any noticeable outliers, which are seen during for PM 2.5 concentrations on the Independence Day. e high values on this day and the outliers accentuate the role play by firework displays across most of the sites. Outliers are also observed for July 05 and July 06 for Port Isabel, close to South Padre Island-a major U.S. tourist destination that organizes massive firework displays.

Spatial and Temporal Variations in PM 2.5 Hourly
Concentrations. COD matrix for PM 2.5 concentrations across the eleven sites and the three CAMS sites are shown in Table 3. COD values equal to greater than 0.2 are in italic and bold. In Brownsville, spatial uniformity in PM 2. e two sampling sites at Weslaco exhibit a COD value 0.3 indicating spatial nonuniformity. In the city of Edinburg, spatial homogeneity was observed across the three sampled sites with COD values less than 0.2. However, the three Edinburg sites demonstrated spatial nonuniformity in the PM 2.5 concentrations when compared with the CAMS C43 site in Mission (E1: 0.27, E2: 0.23, and E3: 0.30). Port Isabel sampling site also exhibited spatial heterogeneity in PM 2.5 pollutant concentrations with all the sampled sites except for B1 (0.19) and B4 (0.17). e COD value (0.32) between PI site and C323 site in Port Isabel also exhibit slight spatial nonhomogeneity.
ese COD findings demonstrate the importance of intra-and interurban sampling of criteria air pollutants like fine particulate matter.
Spearman's rho correlation coefficients were computed between all the study and CAMS sites for fine particulate matter and are shown in Table 4

Conclusions
To the best of our knowledge, the findings from this study are the first in the lower RGV region of South Texas, characterizing fine particulate matter pollution across eleven sites before, during, and after the 4 th of July Independence Day celebrations. ese celebrations are accompanied by the display of massive fireworks display resulting in short-term increase in particulate matter pollution. Whereas, our study findings do not exceed the US NAAQS 24-hr mean PM 2.5 concentrations, even a short-term exposure of a few hours to high PM 2.5 levels can result in respiratory health issues such a wheezing, tightness of chest, and asthma aggravation and subsequent asthma attacks. Findings from this study may, therefore, contribute toward the future formulation of policies for firework display at the local level. is study demonstrates the usefulness of using low-cost sensors to characterize intraurban spatial variability in fine particulate matter concentrations as well as address fine-scale changes in PM 2.5 pollution, which is inadequately captured by federal and state fixed continuous ambient monitoring sites.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.