Next Article in Journal
Comparative Analysis, Use Recommendations, and Application Cases of Methods for Develop Ship Emission Inventories
Previous Article in Journal
Impact of NH3 Emissions on Particulate Matter Pollution in South Korea: A Case Study of the Seoul Metropolitan Area
Previous Article in Special Issue
A Review of Selected Types of Indoor Air Purifiers in Terms of Microbial Air Contamination Reduction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Shale Gas Exploration and Exploitation Activities on the Quality of Ambient Air—The Case Study of Wysin, Poland

by
Janusz Jarosławski
*,
Izabela Pawlak
,
Jakub Guzikowski
and
Aleksander Pietruczuk
Institute of Geophysics Polish Academy of Sciences, 01-452 Warszawa, Poland
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(8), 1228; https://doi.org/10.3390/atmos13081228
Submission received: 21 June 2022 / Revised: 20 July 2022 / Accepted: 25 July 2022 / Published: 2 August 2022
(This article belongs to the Special Issue Air Quality in Poland)

Abstract

:
The continuous two-year monitoring of a set of air pollutants, as well as gases directly related to shale gas exploration processes (methane, non-methane hydrocarbons, carbon dioxide), was carried out at Stary Wiec village in the vicinity (1100 m) of the shale gas wells area in Wysin (Pomeranian voivodeship, north of Poland), covering the stages of preparation, drilling, hydrofracturing and closing of wells. The results of analysis of air pollution data from Stary Wiec and nearby urban and rural stations, over the period 2012–2017 (starting three years before preparations for hydraulic fracturing) indicated that Stary Wiec represents a clean rural environment with an average concentration of nitrogen oxides, carbon monoxide and particulate matter that is one of the lowest in the Pomeranian region. The aim of this study was to explore the range of potential impact of shale gas exploration on local ambient air quality. Analysis of dependence of the concentration level of pollutants on the wind direction indicated that during the drilling period, when the air was coming directly from the area of the wells, nitrogen oxide concentration increased by 13%. Increases of concentration during the hydro-fracturing period, recorded at the Stary Wiec station, were equal to 108%, 21%, 18%, 12%, 7%, 4%, 1% for nitrogen oxide, non-methane hydrocarbons, carbon monoxide, nitrogen dioxide, particulate matter, carbon dioxide and methane. The results of one-minute concentration values for the period 1–4 September 2016 showed a series of short peaks up to 7.45 ppm for methane and up to 3.03 ppm for non-methane hydrocarbons, being probably the result of operations carried out at the area of the wells.

Graphical Abstract

1. Introduction

Natural gas development and exploration can have a significant impact on air quality [1,2]. Shale gas production contributes to emissions into the air, including methane (CH4), nonmethane hydrocarbons (NMHC), carbon dioxide (CO2), and hydrogen sulfide (H2S) [3], although individual stages of gas production generate different sets of air pollutants [4,5]. The first stage of shale gas production (preparation and drilling) engages a large number of diesel-powered engines (e.g., trucks and pumps) that emit mainly nitrogen oxides (NOx), volatile organic compounds (VOCs), CO2, and fine particulate matter (PM2.5 and PM10) [6,7]. Research by Purvis et al. [8] highlights that the preoperational phase, covering 20 weeks of intensive activities related to preparation for hydraulic fracturing, leads to a significant increase of the annual mean of nitrogen oxide (NO) (3-fold) and NOx (2-fold) equal to 5.3 ppb and 10.9 ppb, respectively. Hydraulic fracturing processes often result in uncontrolled CH4 leaks [9]. At the stage of completion of the well, the process of venting shale gas above the surface can result in increased emission of VOCs [10]. The problem of air pollution is linked to direct emission because drilling and hydraulic fracturing may concern not only the borehole area, but the transport of these pollutants to downwind regions, especially those inhabited by people, which is equally important [11]. The measurements performed by Orak et al. [12] in Marcellus Shale (USA) at the turn of 2015 and 2016 indicated significant increase of NOx during the fracturing phase and of C2H6 (ethane) and CH4 for the flowback phase. The values of the 25th percentile equaled 12.5 ppb, 24 ppb and 2.5 ppm for NOx, C2H6 and CH4, respectively, which were significantly higher than for other phases.
Shale gas is regarded as a low-carbon fuel, and its exploitation results in lower CO2, NOx, and sulfur oxides (SOx) emissions than conventional coal or oil combustion [12,13,14]. According to Song et al. [15] replacing coal with shale gas in China reduced the mean monthly concentration of PM2.5 during the time period 2013–2014 by 70%, compared with the same period in 2011–2012. Similar results were obtained by Zhang et al. [16]. Substitution of coal and biomass fuel by shale gas reduced PM emissions in 2017 in China by 123 kt/y. According to Wang et al. [17], the full environmental impact (mainly associated with water consumption and combustion of diesel fuel during the development of the well) of shale gas is 12.5% greater than natural gas. In contrast, replacing fossil fuels (especially coal and biomass) by shale gas could efficiently reduce emission of PM. Among the conventional/unconventional methods of generating electricity and related to the reduction of harmful effects of pollutant emissions on the environment, the investigations with thermal power plants are worth mentioning. Investigations of Abadan Combined Cycle Power Plant (CCPP) in Iran aimed at improving its performance and reducing emissions, indicated that introducing parabolic solar collectors and a multi-effect desalination cycle decreased CO2 emissions from 58.05 to 54.97 kg/MWh of electricity [18]. Further investigations of Abadan (CCPP) showed that the implementation of a multigeneration system caused environmental effects of the modeled system on air pollution production (NOx, CO, CO2) equal to 8.975 kg/s, which indicated a decrease of pollutant emissions equal to 9.8% [19]. Specific configurations of Abadan CCPP with a solar power tower [20] could reduce emissions into the atmosphere up to 345 63 tons/year. There are numerous objections to uncontrolled methane leakages from natural gas systems during exploration activities [13,21]. CH4 is the main component of shale gas and the second main greenhouse gas (GHG), and it is much more efficient in absorbing radiation than CO2 [21]. Many studies have focused on the contribution of natural gas production to CH4 emissions and GHG budget, although their results differ from each other. According to observations by Karion et al. [22], the total fugitive emission of CH4 from shale gas operations can be up to 17.3% of total gas production. Research by Shaw et al. [23] showed the results of air quality measurements performed near the UK’s first shale gas hydraulic fracturing area (Lancashire, England) at the turn of 2018 and 2019. Observations during flowback operations revealed elevated CH4 mixing ratios (30-min mean of up to 5 ppm) after artificial nitrogen lifting of a well. The estimate of mean CH4 flux using modeling approaches was approximately 16 g s−1.
In the vicinity of the borehole area, a significant emission of NOx and VOC contributes to surface O3 formation and results in high O3 episodes [24,25]. Schnell et al. [1] analyzed surface O3 concentrations in the rural area of the Upper Green River Basin, Wyoming (UGRBWY), during February 2008. Under conditions of a high-pressure system, freezing temperature, snow cover, low wind, and clear sky, the hourly average O3 concentration rose from 30 ppb (night) to above 140 ppb (afternoon). Measurements performed by Rappenglück et al. [26] in UGRBWY from January to March 2011 also indicated repeated episodes of the air quality standard (>75 ppb) being exceeded during the winter season. Large amounts of NOx and NMHC emitted during fossil fuel activities accumulate below the shallow stable (≈50 m.a.g.l.) boundary layer, thereby affecting surface O3 formation, which reaches maximum 1-h average values above 160 ppb. Cheadle et al. [27] reported that in the North Colorado Range Front (NCRF) in Colorado, on individual summer days with high O3 levels, emission related to oil and natural gas extraction activities could contribute up to 30 ppb to O3 formation locally. According to Benedict et al. [28] elevated VOC concentrations associated with greater oil and natural gas determine episodes of high O3 concentration in NCFR in Colorado. Emission connected with gas exploitation contributed up to 20 ppbv excess O3. Similar results were obtained by Lindaas et al. [29]. They found that amounts of long-lived VOCs, related to oil and natural gas exploitation, are so high that processes of their oxidation during the day produce sufficient peroxy radicals to significantly impact on O3 production in NCFR in Colorado during high O3 days in this region.
The success of shale gas in the United States at the beginning of the 21st century has become a phenomenon on a global scale, mainly due to advantages of geological and natural physical conditions (shallow shale gas formation with a high index of total organic carbon and relatively big thickness) [30]. Additionally, favorable government policies and regulation systems supported development of local shale infrastructure, advanced techniques and technology, especially hydraulic fracturing. In comparison with the United States, the natural conditions in Europe, in terms of geology, topography and water abundance, are more complicated [31]. Additionally, local hostility and lack of acceptance by the population, due to concerns about extreme environmental and health risks, has appeared in many places where attempts have been made to exploit shale gas [32]. Selected shale gas research approaches, with regards air pollutant measurements are listed in Table 1.
The shale rock formation located in the Wysin area belongs to the Baltic-Podlasie-Lublin Basin. This is the largest unconventional gas reservoir in Poland. According to the most reliable estimations, it contains 346–768 × 109 m3 of recoverable shale gas resources (on average four times more than the documented conventional gas resources in Poland) with a maximum of 1.92 × 1012 m3 [34]. The shale gas exploratory activities in Poland began a few years ago. The first horizontal drilling was performed in 2011 near the Łebień village in the Pomeranian voivodeship (about 90 km from Wysin). According to a series of measurements of potential air pollution indicators (SO2, NOx, CO, C6H6, and H2S) made in Łebień by the Regional Inspectorate for Environmental Protection in Gdańsk, no exceeding of the acceptable level of air pollution were observed [35].
This study covers two- years (from August 2015 to July 2017) of continuous observation of the main atmospheric air pollutants (CO, O3, PM10, NO, NO2, and NOx), as well as pollutants directly associated with shale gas exploration activities, such as drilling and hydraulic fracturing (CH4, NMHC, and CO2), from the Stary Wiec station located in the village of Wysin, about 1 km from the borehole area. The measurement campaign was performed within a framework of the SHEER (Shale gas exploration and exploitation induced risk) project realized under the Horizon 2020 program. The main goal of this work was to present the variation in air pollutant levels at different stages of the natural gas extraction activities and to indicate the potential impact of shale gas exploration activities on the air quality in relation to background levels. For this purpose, a number of statistical analyses were performed: indication of average pollution concentration in different time scales, estimation of high pollutant concentrations, determination of episodes with a particular emphasis on pollutants directly linked to shale gas exploration, and comparison of results obtained at the Stary Wiec station with those from neighboring stations.

2. Materials and Methods

An air pollution monitoring station was set up in the small village of Stary Wiec (54.1° N, 18.3° E, 170 m a.s.l.), the closest inhabited place to the boreholes. The distance between the wells area and the station was approximately 1100 m. (Figure 1).
Location of the station relative to the wells considering the prevailing wind directions in this region. Approximately 55% of the winds originated from the western directions during the measurement campaign. Moreover, only for 0.8% of the measurement time was the wind speed less than 1 m/s. In this situation, the possible emissions of air pollutants from the drilling site had a good chance of being recorded by the measuring equipment installed at the station. The closest vicinity of the station is dominated by small villages with a population of up to 600 inhabitants. The nearest small town (Skarszewy), with a population of more than 7000, is located 8 km east of the station. Local anthropogenic sources of air pollution were small and typical for rural conditions: several neighboring houses and a small local road located 50 m from the station. Residential buildings mainly included single family houses heated by coal- and wood-burning stoves and several multifamily houses connected to the heat distribution networks. The measuring station was located in the northwest part of the village. In the immediate vicinity of the station, there was a closed pig farm. Since the pig farm was nonfunctional, it had no effect on the level of pollution in the vicinity of the station. During the westerly winds (which prevail in this region), when there was an inflow of air from the direction of the borehole area, the combustion processes in these houses had no significant effect as they are on the leeward side of the station. The road was occasionally used, with traffic of a dozen cars per hour. Additional heavy traffic, related to activities at the wells area, occurred during the period of exploitation of the wells. In addition, a park with a small pond located in the vicinity of the station was a potential natural source of gases, such as CO2 and CH4. The surrounding landscape was dominated by farmlands and small forests. The nearest significant urban sources of air pollution were located at a distance of 25–50 km from the station (see Figure 1).
The station’s container was equipped with a standard set of analyzers measuring concentrations of the air pollutants: NO, NO2, NOx (sum of NO and NO2), CO, O3, and PM10. Gaseous pollutants were measured by the API analyzer, models 200, 300, and 400, while PM10 concentration was measured by a TEOM analyzer, model 1400. This standard air pollution monitoring measurement program was extended by measurements of the concentration of CO2 using a Thermo 410i analyzer and concentrations of CH4 and NMHC (sum of NMHCs) with a Horiba APHA370 analyzer. Technical specifications of the analyzers employed during the measuring campaign are listed in Table A1 in the Appendix A. Meteorological parameters were recorded by a collocated automatic weather station.
To ensure adequate data quality, i.e., total measurement uncertainty less than 15% in accordance with the requirements used in air quality monitoring networks in Poland, standard quality assurance procedures analogous to those routinely used in air pollution monitoring networks were applied during measurements, including regular calibrations with certified gas mixtures and a certified photometer (in the case of ozone), data verification and validation, and regular control of, and servicing of, analyzers and container systems. Instantaneous values of all pollutants’ concentrations were converted into 1-h averages, and these averages were used for further analysis.
The measuring campaign in Stary Wiec covered a wide range of air pollutants, both main gaseous pollutants measured as standard at state air quality monitoring stations (NO, NO2, NOx, CO, O3, PM10) and, additionally, those directly related with shale gas related activities (CO2, CH4, NMHC). The measurements were performed continuously for two years, from August 2015 to July 2017. This long period of uninterrupted measurements ensured detection of possible air pollution episodes related to all stages of the life cycle of wells: preparatory phase, drilling, hydrofracturing, and closing activities. Additionally, 24 months of monitoring covered two full seasonal cycles, which was important for identifying the natural local variability of each pollutant and possible local sources of pollution not related to the wells’ activities, thus allowing the determination of the background values of each analyzed pollutant. In this work we investigated the shale gas industrial development in a rural area in the North of Poland in terms of accompanying air pollutants and GHGs emissions. In this connection, the analysis of the potential impact of activities related to shale gas exploration and exploitation processes on ambient air quality, in the context of environmental protection and the local community, was performed. The comparative analysis of pollutant concentrations during an entire work cycle of the wells (as well as in the period before and after), taking into account meteorological conditions, was considered in terms of the following: (1) whether the National ambient Air Quality Criteria were exceeded (2) occurrence of peaks and deviations from the background values (3) comparison concentrations of all pollutants during different shale gas operation life cycles (drilling and hydrofracturing) relative to the reference period (in this article defined as a remaining period). In order to perform this, we first used the data on the levels of air pollution and meteorological conditions from several of the nearest (7 to 30 km) stations to the regional air pollution monitoring network to assess air quality in the region before and during the measurement campaign. The analysis of the measurement data (from 2012 to 2017), presented in Figure 2, indicated the lack of significant difference in ambient air pollutants between the baseline period (2012–2014) and the period of intensive activity in the shale gas area (2015–2017). To enable the estimation of the potential impact of shale gas exploration on air quality, an examination on the influence of wind direction on the measured mixing ratios of each pollutant was performed (Figure A1, Figure A2, Figure A3, Figure A4 and Figure A5 in Appendix A). Air pollution data were assigned to one of eight wind direction compartments (N, NE, E, SE, S, SW, W, NW) depending on the phase of activity at the area of the wells. Furthermore, an additional compartment for wind coming directly (±5°) from the borehole area was determined. Wind rose graphs were used to present the amounts of pollutants averaged for each sector of the wind and to compare these with the amounts estimated for the wind coming directly from the shale gas site. To determine whether or not the air quality standards for human health and plant protection were met the compliance of NO, NO2, CO, O3 and PM10 with the Air Quality Directive was evaluated. Finally, diurnal and seasonal analyses of CH4, NMHC and CO2 as the gases directly related to shale gas operations were carried out (Figure 3). The detailed analysis of the high temporal resolution data revealed the existence of episodes of elevated concentrations of CH4 and NMHC registered in the period 1–4 September 2016 (Figure 4). To determine the origin (artificial or biogenic) of high GHGs concentration calculations of CO2/CH4 ratio were performed (Figure 5). Additionally, statistical analysis of episodes recorded for the period 1–4 September 2016, based on correlation coefficient, p-value and CH4/NMHC ratio values, was performed to confirm whether or not the shale gas related emissions were a potential source of episodes of elevated concentrations of CH4.

3. Results

3.1. Analysis of Measurement Results from the Regional Network of Air Pollution Monitoring

Measurement results of pollutants, namely NO, NO2, CO, O3, and PM10, from the five nearest air pollution monitoring stations, namely, Gdańsk, Starogard Gdański, Kościerzyna, Tczew (reflecting urban background conditions), and Liniewko Kościerskie (reflecting rural background conditions), were used for the general assessment of air quality level in the region before starting any activities related to shale gas exploration and during the measurement campaign in Stary Wiec. Data from the period 2012–2017 were used for the analysis. Results of the available measurements (monthly means) for primary (NO) and secondary (O3) pollutants are shown in Figure 2.
The analysis of data revealed distinct natural seasonal and year-to-year variations in concentrations of all measured air pollutants in the region. This variation resulted especially from changes in meteorological conditions during the year and changes in emissions from local anthropogenic sources (e.g., traffic and heating systems of buildings in the winter). The period 2015–2017 did not differ from previous years in terms of concentrations of individual pollutants (see Table A2 in the Appendix A). Statistics of pollutant levels measured at Stary Wiec and at nearby rural and urban stations during the measurement campaign are presented in Table 2.
Air pollutant levels measured in Stary Wiec were representative of a clean, rural environment. Very low, close to zero, levels of NO, a pollutant with a lifetime of approximately 1 h, indicated the lack of significant engine exhaust type sources of air pollution in the vicinity of the station. Pollutants such as CO, NO2, and PM10 were the lowest, or one of the lowest, in the region during the measurement period. However, the levels of the secondary pollutant O3 were at a relatively high level, which is a known typical phenomenon for rural areas [36]. Higher O3 concentrations at rural stations than at urban stations are due to less efficient mechanisms for O3 destruction [37], the most important of which is the lack of efficient NO emission sources that can reduce the reactions of O3 destruction (NO+O3→NO2+O2). Most episodes of elevated pollutant concentrations registered at the Stary Wiec station during the 2015–2017 period resulted from regional scale phenomena, such as favorable meteorological conditions or transport of pollutants from regional sources, e.g., large cities of the Pomeranian Voivodeship.

3.2. Analysis of Air Pollutant Measurement Results from Stary Wiec Station

As described in Section 1, there were no significant sources of air pollution in the immediate vicinity of the station; hence, it can be assumed that the levels of pollutants in the area of the station were largely due to their transportation to the place of measurement from more distant source areas. Thus, pollutants emitted in the Wysin wells area could have been registered by the station’s equipment if they were transported to the sampling point when favorable wind conditions occurred.
To assess the possible impact of pollutant emissions from the wells area on the ambient air quality in the vicinity of the station, an analysis of the distribution of air pollutant concentrations recorded at the station, depending on wind direction, was performed. For this purpose, the average concentrations of individual air pollutants were calculated for eight wind direction compartments. In addition, the average concentrations of pollutants for wind blowing directly (±5°) from the wells were determined. The measurement time was divided into three periods in accordance with the actions taken at the boreholes: the drilling period (September–October 2015), the hydrofracturing period (June–July 2016), and the remaining period of measurements (August 2015–July 2017, excluding drilling and hydrofracturing periods).
Significance of differences between the average values of concentrations of pollutants from the NW (north-west) sector and for air originating directly from the well was confirmed by the F-test (to determine whether the analyzed samples were characterized by equal/unequal variance) and the Student’s t-test with two samples assuming equal or unequal variances (to determine whether the differences of mean values in two analyzed samples were statistically significant at the significance level of 0.05). Results are shown in Table 3.
Measurement results are shown in Figure A1, Figure A2, Figure A3, Figure A4 and Figure A5 in the Appendix A. These figures show the dependence of the concentrations of pollutants on the wind direction for a previously defined period. It can be observed that the cleanest air conditions occurred during advection of air from the western directions, while the highest levels of pollutants appeared while air was coming to the measuring point from the east or south. This situation roughly corresponded to the distribution of sources of air pollution around the measurement site. Cleaner air coming from the West was often of North or West Atlantic origin, while the air from the eastern sectors, mainly of continental origin, passed over relatively close significant sources of air pollution (several larger cities and the motorway located nearby) on its way to the measuring point.
During the low activity period in the wells area, the concentrations of air pollutants registered at the Stary Wiec station did not increase significantly when direct transport from the area of the wells occurred, in comparison to that when the air was coming from the entire NW sector. The relative differences in concentrations were small, in the order of a few percent (up to 5%), except for NO (−13%). For most air pollutants, nonsignificant changes in concentrations were recorded in the air coming directly from the well region as compared to their average values for the whole sector. Significantly lower values were only recorded for O3; however, the absolute differences remained very small (approximately 1 ppb).
A different situation occurred during the periods of increased activity at the wells, especially during hydrofracturing. Most of the pollutants occurred in lower concentrations during this period, because of natural seasonal variation. Hydrofracturing was performed in the summer, when the levels of air pollutants (except for O3) were usually at the lowest of the year. However, the levels of some of the pollutants in the air arriving at the measurement point directly from the wells area increased significantly according to the t-test results, compared to the average values registered at the same time for the air coming from the NW sector during the hydrofracturing period. A significant relative increase (108% and 18%) occurred for NO and CO, respectively. Increase in the levels of NO2 (12%), PM10 (7%), and NMHC (21%), although not statistically significant, were also detected. All the aforementioned pollutants were the components of engine exhausts. Therefore, an observed increase in air pollution levels could be the result of intense work of, for example, the pump systems. Moreover, smaller (by 11%) levels of O3 observed simultaneously with higher NO concentrations confirmed this hypothesis, as O3 levels are usually anti-correlated with NO levels, especially when close to the source of emission, because of chemical reactions occurring between them. Slight increases in the concentrations of CH4 and CO2 were also detected.
For the drilling period, the results were ambiguous. A statistically nonsignificant increase in NOx (by approximately 20%) occurred. However, for the other pollutants, the concentrations were usually lower than the average for the NW sector, up to 19% lower for PM10. In two cases, those relating to PM10 and CH4, the lower values were statistically significant.
The analysis of the above results showed that, despite a relatively high increase in the concentrations of some pollutants during the hydrofracturing period, their absolute concentrations remained at a very low level; for example, 1.2 ppb for NO or 5.0 ppb for NO2, during the hydrofracturing and drilling periods, respectively. Thus, even in the periods of intense activity at the boreholes, the air in the vicinity of the Stary Wiec station remained relatively clean. The analysis of compliance with the Air Quality Directive applicable to the pollutants (Table 4) also showed that the relevant air quality criteria were met for most pollutants.
By analyzing the air quality standards for human health protection, it can be stated that during the two-year measuring period, most noticeable were the relatively high values in ambient O3 concentrations. The information target value was exceeded twice in August 2015. Although the air quality limit values for O3 were frequently exceeded during summer months, it should be noted that exceedances registered in the period 2015–2017 originated from episodes of a larger, regional scale, which were not related to works at the area of the wells. Permissible levels of the remaining pollutants (NO2, PM10, and CO) for the protection of human health were not exceeded.

3.3. Methane, NMHCs, and CO2 Concentration Measurement Results

The monitoring of CH4, NMHC, and CO2 levels was included into Stary Wiec station’s measurement program to track the possible emissions of these gases accompanying the exploration, and subsequent exploitation, of shale gas wells in Stary Wiec. The analysis of the data collected during the two years of measurements showed the existence of well-known seasonal (summer minimum and winter maximum) and diurnal (minimum during the day and maximum at night) cycles of CO2 and CH4 related to the activity of natural sources and sinks, namely photosynthesis and respiration of the biosphere for CO2 and methanogenesis for CH4. CO2 is the basic substrate of the photosynthesis process occurring with the participation of light energy, most often solar energy. Photosynthesis processes are an important sink of CO2 in the natural ecosystem during the daytime. In contrast, the processes of respiration of plants and soil are significant natural sources of CO2 and contribute to its increase during the night (in the absence of photosynthesis). Methanogenesis is the process of biological formation of CH4 during decomposition of organic compounds (e.g., acetate, formate, CO2, hydrogen) by anaerobic microorganisms (mainly from the Archaea domain). Figure 3 shows the results of measurements of diurnal variations of CO2 and CH4 performed over the entire experimental period.
Both gas species had similar, distinct diurnal distribution, characterized by maximum values during late night–early morning (02:00–04:00 for CO2 and 03:00–05:00 for CH4), followed by a gradual decrease, reaching a minimum around 12:00–16:00, and then a gradual increase through the evening and night until the maximum peak the next day. Thus, it can be stated that the concentration values of both gases (CO2 and CH4) were cumulated and constantly increasing in the afternoon through evening and night until an early morning maximum on the next day. The maximum values of CO2 varied between 431.5 ppm and 433.6 ppm, while the lowest value was around 399 ppm. Daily concentration values of CO2 were characterized by a very high amplitude in the summer months, with the maximum reaching around 100 ppm. Hence, high values probably resulted from the presence of important natural sources and sinks in the vicinity of the station (pond and park with deciduous trees). During winter, the amplitude values were almost negligible, oscillating around several ppm. The maximum values of CH4 increased to around 1.962 ppm, while the minimum dropped to 1.916 ppm. Stronger amplitudes of diurnal variation were noted during summer, with the maximum values in July 2016, reaching 0.125 ppm.
The mean values at Stary Wiec for the recording period of August 2015–July 2017 were 1.936 ± 0.064 ppm for CH4, 414 ± 29 ppm for CO2, and 0.039 ± 0.034 ppmc for NMHC. These levels were equal to, or slightly higher than, those reported for the Poland remote mountain region [38,39].
Concentrations of all three analyzed components in the air coming directly from the wells were higher than the average values for the NW sector in the hydrofracturing period, ranging from 1% for CH4 up to 21% for NMHC. However, this increase was not statistically significant. For the remaining period, the corresponding differences were below 1% for CH4 and CO2 and below 5% for NMHC.
Variability in CH4 concentrations recorded at Stary Wiec was relatively small (standard deviation 0.06 ppm); hence, more attention was paid to the analysis of episodes of elevated CH4 concentrations, as potentially coming from the wells. The highest concentration of CH4 (1-h average), 2.81 ppm with more than 10 standard deviations from the average value, was registered on 1 September 2016, outside the drilling and hydrofracturing period. A more detailed analysis of the one-minute concentrations of CH4, CO2, and NMHC for this, and for several subsequent days, revealed a series of short CH4 and NMHC concentration peaks, up to 7.45 ppm for CH4 and up to 3.03 ppm for NMHC. CO2 concentrations remained at mean levels. The results of the measurements are shown in Figure 4.
The configuration and times of occurrence of the concentration peaks were similar for three consecutive days, with decreasing concentration values, indicating their artificial origin. The wind direction analysis for this period indicated the wells area as a possible methane source.
Analysis of CO2 and CH4 concentration ratio (ppm/ppm) is sometimes used as a tracer in investigating the origin of GHGs [40]. A low CO2/CH4 ratio is often interpreted as an indicator of significant anthropogenic impact [41]. This analysis was performed for the time of occurrence of the abovementioned CH4 and NMHC episodes. The graph of these ratios is shown in Figure 5.
Usually, its values remained relatively stable during the daytime, and the mean value for the entire measurement period was 219 ± 12, with an increase during nights, because of strong CO2 production through vegetation in a park near the station. During the duration of the episodes, the ratio systematically dropped down to 53, indicating a non-biogenic source of CH4, probably due to some technical procedures at the wells, which was associated with the release of certain amounts of CH4 directly into the atmosphere. Unfortunately, the site operator did not provide information about the type of activity being conducted at this time at the Wysin site.
Additionally, to investigate the predominance of the thermogenic origin of increased CH4 concentration during the episodes, an analysis, based on the relationship between CH4 and NMHC, was performed. Characteristics of the total number of episodes observed for the period 1–4 September 2016 are presented in Table 5.
The concentration values of CH4 (1-min average) qualified as episodes were greater than the average value by at least 5 deviations, reaching the maximum value of 7.45 ppm (mean 2.91 ± 0.92 ppm), whereas almost all the values of NMHC were greater than the average value by at least 3 standard deviations. The maximum value of NMHC was 3.03 ppm (mean 0.43 ± 0.49 ppm). A significant correlation was observed between CH4 and NMHC levels (p < 0.05) with correlation coefficient above 0.6, and a low value of the CH4/NMHC ratio was recorded in almost all episode cases (except episode numbers 8 and 9). Correlation coefficient values, in most cases, were above 0.9, thus indicating a strong association between species. The CH4/NMHC ratio values were significantly lower than the mean value for the entire experiment (90) and dropped down to approximately 5. These calculations indicated the possible influence of a thermogenic source of the CH4 episodes. The origin of episode numbers 8 and 9 is not clear. Correlation coefficient values were slightly smaller than the others (0.680 and 0.859, respectively), and p-values (0.064 and 0.062, respectively) were above the adopted significance (p < 0.05). The CH4/NMHC ratio values were the highest of all (15.46 and 23.01, respectively), but remained significantly lower than the average value for the entire measurement period. The potential sources of elevated CH4 levels in episodes 8 and 9 may indicate that the CH4/NMHC ratio included both thermogenic and biogenic contributions, but these inferences are not unequivocal.
For episode number 1 the methane plume possibly emitted from the well was modeled by the HYSPLIT model. The results are shown on Figure A6 in the Appendix A.

4. Conclusions

This work presents the first long-term (two year) measurement assessment of the air quality impact of activities performed at a shale gas site in Poland.
The development of a shale gas exploration and exploitation site in Wysin (over the period 2015–2017) consisted of a full cycle of technical works, including the following phases, each of which has potential for pollution emission to the atmosphere: pre-operational phase, well drilling, hydro-fracturing, closing the well and site restoration. Continuous monitoring of levels of air pollutants in the vicinity of shale gas wells supported the hypothesis that emission from shale gas extraction activities in the Wysin area has a negligible impact on ambient air quality at a distance of about 1100 m from the wells on their leeward side.
Concentrations of pollutants measured at the station strongly depend on the kind of activities performed at the borehole area. The highest, partly statistically significant increases in level of selected pollutants were noted during increased activity at the wells (during drilling and hydro-fracturing periods) in the conditions of direct transport of air from the area of the wells to the measurement point. The works of drilling and hydro-fracturing involved intensified vehicle activity, due to delivery of equipment to the site and its operation, and led to significant increases, especially in NOx emissions. Concentration of NO during intensive activity at the borehole area increased to an average value of 1.72 ppb (during drilling) and 1.21 ppb (during hydro-fracturing), compared to 0.83 ppb measured during the remaining period when no activities related to the borehole were performed.
Based on the observations of CH4, CO2 and NMHC, the gases directly related to shale gas exploration activities, the average concentrations registered for the period 2015–2017 were equal to 1.936 ± 0.064 ppm for CH4, 414 ± 29 ppm for CO2, and 0.039 ± 0.034 ppm for NMHC. During the measurement period, several episodes of elevated CH4 concentration values, with maximum one-hour average equal to 2.81 ppm (more than 10 standard deviations from the average value). occurred in September 2016. Detailed analysis of the one-minute concentrations of CH4, CO2 and NMHC revealed a series of short CH4 and NMHC concentration peaks, up to 7.45 ppm for CH4 and up to 3.03 ppm for NMHC. The time of appearance of these peaks was similar for three consecutive days, with decreasing concentration values, which could indicate their artificial origins. CH4/CO2 ratio values during the measurement period were characterized by relatively stable values ≈ 219 ± 12. Over 1–4 September CH4/CO2 ratio values dropped down to 53. Due to lack of information from the site’s operator, we cannot unambiguously state what kind of operations were performed during these days at the borehole area and what exactly caused the high emission of CH4 and NMHC pollutants. However, analysis of wind direction during these days indicated the borehole area as a potential CH4/NMHC source.
Results of the analysis of measurements taken during shale gas exploration activities indicated that significant increases of concentrations of air pollutants (expect of NOx, CH4 and NMHC, mentioned before) were not observed. Taking into account the Air Quality Criteria, applicable to individual pollutants, no exceedances of threshold values (with the exception of O3) were observed. Relatively low levels of concentrations of air pollutants measured during drilling, hydraulic fracturing and shale gas extraction at Stary Wiec could be determined, especially by the scale of operations across the time and space, which were relatively small. Additionally, exploration activities concerned only two wells. It can be assumed that the use of shale gas extraction technology, taking into account: risk-mitigation measures and technical quality control, could ensure environmentally safe exploration and exploitation of shale gas.

Author Contributions

Conceptualization, J.J. and I.P.; methodology, J.J.; formal analysis, I.P. and J.G.; writing—original draft preparation, I.P. and J.J.; writing—review and editing, J.J.; visualization, I.P.; corrections and consultations, A.P.; supervision, J.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by SHEER: “Shale Gas Exploration and Exploitation Induced Risks” project funded from Horizon 2020–R&I Framework Programme, call H2020-LCE-2014-1, Grant Agreement No. 640896 and partially supported by the Polish Ministry of Education and Science (Decision no. 13/E-41/SPUB/SP/2020) and the Chief Inspectorate of Environmental Protection, Poland under Grant number GIOŚ/19/2021/DMŚ/NFOŚ.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Air pollutant concentrations measured at the Stary Wiec station vs. wind direction
Figure A1. Average levels of NO and NO2 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. All concentration units are in parts per billion (ppb). The red bars indicate the airmass coming directly from the area of the wells.
Figure A1. Average levels of NO and NO2 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. All concentration units are in parts per billion (ppb). The red bars indicate the airmass coming directly from the area of the wells.
Atmosphere 13 01228 g0a1
Figure A2. Average levels of CO and PM10 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. Concentration units are in parts per billion (ppb), except for PM10 (μg/m3). The red bars indicate the airmass coming directly from the area of the wells.
Figure A2. Average levels of CO and PM10 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. Concentration units are in parts per billion (ppb), except for PM10 (μg/m3). The red bars indicate the airmass coming directly from the area of the wells.
Atmosphere 13 01228 g0a2aAtmosphere 13 01228 g0a2b
Figure A3. Average levels of O3 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. All concentration units are in parts per billion (ppb). The red bars indicate the airmass coming directly from the area of the wells.
Figure A3. Average levels of O3 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. All concentration units are in parts per billion (ppb). The red bars indicate the airmass coming directly from the area of the wells.
Atmosphere 13 01228 g0a3
Figure A4. Average levels of CO2 and CH4 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. Concentrations units for CO2 and CH4 are parts per million (ppm).
Figure A4. Average levels of CO2 and CH4 vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. Concentrations units for CO2 and CH4 are parts per million (ppm).
Atmosphere 13 01228 g0a4
Figure A5. Average levels of NMHC vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. Concentration units for NMHC are ppmc.
Figure A5. Average levels of NMHC vs. wind direction measured at the Stary Wiec station during drilling (DR), hydrofracturing (HF), and for the remaining period of measurements. Concentration units for NMHC are ppmc.
Atmosphere 13 01228 g0a5
Methane was the pollutant for which there was the most distinct influence of shale gas extraction activities on ambient air quality at the Stary Wiec station. From 1 to 3 September 2016, a series of short (registered on one-min average data) methane concentration peaks, with a maximum value equal to 7.45 ppm, were noticed. The way in which the repeatability of these peaks occurred indicated a probable artificial origin of the episode, whereas analysis of wind direction revealed the borehole area as the potential source. To confirm these assumptions, a 3D numerical simulation of dispersion of CH4 in the area surrounding a station was performed. The borehole area located to the north-west of the station was selected as an emission point. The results of modeling for 01.09.2016 7 pm are shown in Figure A6 A pollutant plume emitted at the well location extended from the NW and passed through the measurement site. The calculation of the pollutant emission rate on the basis of hourly mean pollutant concentration had height 4 m a.g.l. points to the value of 120 kg/h.
The numerical study used the Weather Research and Forecasting (WRF) Model [42] as a source of meteorological data and the HYSPLIT dispersion model [43] for the CH4 modeling.
The WRF Model is a mesoscale numerical weather prediction system, dedicated to atmospheric research and forecasting of applications. It contains two dynamical cores and serves as a data assimilation system and software architecture, enabling parallel calculations and system extensibility. The Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT), developed by NOAA’s Air Resources Laboratory, is one of the most widely used models for atmospheric trajectory and dispersion calculations. This model is widely used for both determining atmospheric trajectory and dispersion calculations.
The initial and boundary meteorology conditions for the WRF Model were used from the NCAR GFS [44] system.
  • WRF settings:
  • Domain solution
  • d1: 14 km ×14 km 51 vertical level 1 h
  • d2: 4.5 km × 4.5 km 51 vertical level 1 h
  • d3: 1.2 km × 1.2 km 51 vertical level 15 min
  • d4: 0.4 km × 0.4 km 51 vertical level 1 min
  • Simulation start: 01.09.2016 12 UTC end 02.09.2016 00 UTC
The Physics and microphysics schema configuration for WRF Model was performed based on the study by Werner et al. [45].
  • HYSPLIT settings:
  • Domain solution:
  • Setting for the domain was 0.4 km × 0.4 km. Vertically it was 40 levels between 0 m to 224 m above ground level.
Figure A6. Methane plume modeled by the HYSPLIT model for 1 September 2016, revealed the following values: CH4 = 2.81 ppm (hourly mean), NMHC = 0.2 ppm, estimated CH4 emission rate = 120 kg/h. Left and horizontal axis-geographical coordinates, right axis-additional methane concentration above the Stary Wiec background level, i.e., approximately 1.9 ppm. Red circle–well location, black circle-monitoring site.
Figure A6. Methane plume modeled by the HYSPLIT model for 1 September 2016, revealed the following values: CH4 = 2.81 ppm (hourly mean), NMHC = 0.2 ppm, estimated CH4 emission rate = 120 kg/h. Left and horizontal axis-geographical coordinates, right axis-additional methane concentration above the Stary Wiec background level, i.e., approximately 1.9 ppm. Red circle–well location, black circle-monitoring site.
Atmosphere 13 01228 g0a6
Table A1. Technical specifications of instruments employed during measuring campaign.
Table A1. Technical specifications of instruments employed during measuring campaign.
Measured PollutantNO/NO2/NOxCOO3PM10CO2CH4/NMHC
InstrumentAPI 200API 300API400TEOM 1400Thermo 410iHoriba APHA370
Measurement techniqueChemiluminescenceInfrared absorptionAbsorption of UV lightTapered Element Oscillating Microbalance (TEOM)NDIR technologySelective combustion+cross modulation+flame ionization detection
Zero drift<0.5 ppb/24 h<20 ppb/24 h<1.0 ppb/24 h ±1 ppm±0.1 ppmC/day
Lower detectable limit0.4 ppb<20 ppb<0.6 ppb<5 µg/m31 ppm0.02 ppmC (3σ)
Precision0.5%0,5%0,5%0.75%1%1%
Range0–50 ppb to 0–20 ppm0–100 ppb to 0–100 ppm100 ppb to 10 ppm0–5 mg/m3 0 ppmC–50 ppmC
Table A2. Statistics for air pollution data from air quality monitoring stations based on hourly means, collected from January 2012 to July 2015. Mean values for this period, together with standard deviations (1σ), are presented.
Table A2. Statistics for air pollution data from air quality monitoring stations based on hourly means, collected from January 2012 to July 2015. Mean values for this period, together with standard deviations (1σ), are presented.
LiniewkoKościerzynaTczewStarogardGdańsk
NO (ppb)0.86 ± 0.905.72 ± 14.334.30 ± 10.187.37 ± 11.562.93 ± 6.64
NO2 (ppb)3.38 ± 2.448.30 ± 6.007.30 ± 5.608.34 ± 5.977.67 ± 6.03
O3 (ppb)30.14 ± 11.8922.20 ± 12.81 23.97 ± 12.13
CO (ppb)255.36 ± 107.12437.45 ± 307.92309.61 ± 189.51386.57 ± 278.28266.34 ± 136.86
PM10 (µg/m3) 39.73 ± 36.8022.41 ± 17.3739.65 ± 39.3117.06 ± 12.81

References

  1. Schnell, R.C.; Oltmans, S.J.; Neely, R.R.; Endres, M.S.; Molenar, J.V.; White, A.B. Rapid photochemical production of ozone at high concentrations in a rural site during winter. Nat. Geosci. 2009, 2, 120–122. [Google Scholar] [CrossRef]
  2. Katzenstein, A.S.; Doezema, L.A.; Simpson, I.J.; Blake, D.R.; Rowland, F.S. Extensive regional atmospheric hydrocarbon pollution in the southwestern United States. Proc. Natl. Acad. Sci. USA 2003, 100, 11975–11979. [Google Scholar] [CrossRef] [PubMed]
  3. Skone, T.J. Life Cycle Greenhouse Gas Inventory of Natural Gas Extraction, Delivery and Electricity Production; Final Report; National Energy Technology Laboratory: Pittsburgh, PA, USA, 2011; Available online: http://www.canadiancleanpowercoalition.com/files/4713/2630/3388/LA7%20-%20NG-GHG-LCI.pdf (accessed on 27 March 2018).
  4. Sun, R.; Wang, Z. A comprehensive environmental impact assessment method for shale gas development. Nat. Gas Ind. B 2015, 2, 203–210. [Google Scholar] [CrossRef]
  5. Roy, A.A.; Adams, P.J.; Robinson, A.L. Air pollutant emissions from the development, production, and processing of Marcellus Shale natural gas. J. Air Waste Manag. Assoc. 2014, 64, 19–37. [Google Scholar] [CrossRef]
  6. U.S. Environmental Protection Agency (EPA). Nonroad Engines, Equipment and Vehicles: Locomotives. 2013. Available online: http://www.epa.gov/otaq/locomotives.htm (accessed on 22 March 2018).
  7. U.S. Environmental Protection Agency (EPA). Nonroad Engines, Equipment and Vehicles: Nonroad Diesel Engines. 2013. Available online: http://www.epa.gov/otaq/nonroad-diesel.htm (accessed on 22 March 2018).
  8. Purvis, R.M.; Lewis, A.C.; Hopkins, J.R.; Wilde, S.E.; Dunmore, R.E.; Allen, G.; Pitt, J.; Ward, R.S. Effects of ‘pre-fracking’ operations on ambient air quality at a shale gas exploration site in rural North Yorkshire, England. Sci. Total Environ. 2019, 673, 445–454. [Google Scholar] [CrossRef] [PubMed]
  9. Brandt, A.R.; Heath, G.A.; Kort, E.A.; O’Sullivan, F.; Pétron, G.; Jordaan, S.M.; Tans, P.; Wilcox, J.; Gopstein, A.M.; Arent, D.; et al. Methane Leaks from North American Natural Gas Systems. Science 2014, 343, 733–735. [Google Scholar] [CrossRef] [PubMed]
  10. Armendariz, A. Emissions from Natural Gas Production in the Barnett Shale Area and Opportunities for Cost-Effective Improvements. Ph.D. Thesis, Southern Methodist University, Dallas, TX, USA, 2009. [Google Scholar]
  11. Kemball-Cook, S.; Bar-Ilan, A.; Grant, J.; Parker, L.; Jung, J.; Santamaria, W.; Mathews, J.; Yarwood, G. Ozone impacts of Natural gas development in the Haynesville shale. Environ. Sci. Technol. 2010, 44, 9357–9363. [Google Scholar] [CrossRef] [PubMed]
  12. Orak, N.H.; Reeder, M.; Pekney, N.J. Identifying and quantifying source contributions of air quality contaminants during unconventional shale gas extraction. Atmos. Chem. Phys. 2021, 21, 4729–4739. [Google Scholar] [CrossRef]
  13. Burnham, A.; Clark, C. Examining the Impacts of Methane Leakage on Life-Cycle Greenhouse Gas Emissions of Shale and Conventional Natural Gas; Air & Waste Management Association: Pittsburgh, PA, USA, 2012; pp. 8–12. [Google Scholar]
  14. De Gouw, J.A.; Parrish, D.D.; Frost, G.J.; Trainer, M. Reduced emissions of CO2, NOx, and SO2 from US power plants owing to switch from coal to natural gas with combined cycle technology. Earth’s Future 2014, 2, 75–82. [Google Scholar] [CrossRef]
  15. Song, W.; Chang, Y.; Liu, X.; Li, K.; Gong, Y.; He, G.; Wang, X.; Christie, P.; Zheng, M.; Dore, A.J.; et al. A Multiyear Assesssment of Air Quality Benefits from China’s Emerging Shale Gas Revolution: Urumqi as a Case Study. Environ. Sci. Technol. 2015, 49, 2066–2072. [Google Scholar] [CrossRef] [PubMed]
  16. Zhang, Y.; Wang, H.; Han, Y.; Wang, D.; Zhu, G.; Lu, X. Shale gas development in China: Implications for indor and outdoor air quality and greenhouse gas emissions. Environ. Int. 2020, 141, 105727. [Google Scholar] [CrossRef]
  17. Wang, S.; Tang, X.; Wang, J.; Zhang, B.; Sun, W.; Höök, M. Environmental impacts from conventional and shale gas and oil development in China considering regional differences and well depth. Resour. Conserv. Recycl. 2021, 167, 105368. [Google Scholar] [CrossRef]
  18. Javadi, M.A.; Khalaji, M.; Ghasemials, R. Exergoeconomic and environment al analysis of a combined Power and water desalination plant with patrabolic solar collector. Desalination Water Treat. 2020, 193, 212–223. [Google Scholar] [CrossRef]
  19. Ghasemiasl, R.; Abhari, M.K.; Javadi, M.A.; Ghomashi, H. 4E investigating of a combined power plant and converting it to a multigeneration system to reduce environmental pollutant production and sustainable development. Energy Convers. Manag. 2021, 245, 114468. [Google Scholar] [CrossRef]
  20. Javadi, M.A.; Najafi, N.J.; Abhari, M.K.; Jabery, R.; Pourtaba, H. 4E analysis of three different configurations of a combined cycle power plant integrated with a solar power tower system. Sustain. Energy Technol. Assess. 2021, 48, 101599. [Google Scholar] [CrossRef]
  21. Meng, Q. The impacts of fracking on the environment: A total environmental study paradigm. Sci. Total Environ. 2017, 580, 953–957. [Google Scholar] [CrossRef]
  22. Karion, A.; Sweeney, C.; Pétron, G.; Frost, G.; Hardesty, R.M.; Kofler, J.; Miller, B.R.; Newberger, T.; Wolter, S.; Banta, R.; et al. Methane emissions estimate from airborne measurements over a western United States natural gas field. Geophys. Res. Lett. 2013, 40, 4393–4397. [Google Scholar] [CrossRef]
  23. Shaw, J.T.; Allen, G.; Pitt, J.; Shah, A.; Wilde, S.; Stamford, L.; Fan, Z.; Ricketts, H.; Williams, P.I.; Bateson, P.; et al. Methane flux from flowback operations at a shale gas site. J. Air Waste Manag. Assoc. 2020, 70, 1324–1339. [Google Scholar] [CrossRef] [PubMed]
  24. Ahmadov, R.; McKeen, S.; Trainer, M.; Banta, R.; Brewer, A.; Brown, S.; Edwards, P.M.; de Gouw, J.A.; Frost, G.J.; Gilman, J.; et al. Understanding high wintertime ozone pollution events in an oil—And natural gas—Producing region of the western US. Atmos. Chem. Phys. 2015, 15, 411–429. [Google Scholar] [CrossRef]
  25. Oltmans, S.; Schnell, R.; Johnson, B.; Pétron, G.; Mefford, T.; Neely III, R. Anatomy of wintertime ozone associated with oil and natural gas extraction activity in Wyoming and Utah. Elem. Sci. Anthr. 2014, 2, 24. [Google Scholar] [CrossRef]
  26. Rappenglück, B.; Ackermann, L.; Alvarez, S.; Golovko, J.; Buhr, M.; Field, R.A.; Soltis, J.; Montague, D.C.; Hauze, B.; Adamson, S.; et al. Strong wintertime ozone events in the Upper Green River basin Wyoming. Atmos. Chem. Phys. 2014, 14, 4909–4934. [Google Scholar] [CrossRef]
  27. Cheadle, L.C.; Oltmans, S.J.; Pétron, G.; Schnell, R.C.; Mattson, E.J.; Herndon, S.C.; Thompson, A.M.; Blake, D.R.; McClure-Begley, A. Surface ozone in the Colorado northern Front Range and the influence of oil and gas development during FRAPPE/DISCOVER-AQ in summer 2014. Elem. Sci. Anthr. 2017, 5, 61. [Google Scholar] [CrossRef]
  28. Benedict, K.B.; Zhou, Y.; Sive, B.C.; Prenni, A.J.; Gebhart, K.A.; Fisher, E.V.; Evanoski-Cole, A.; Sullivan, A.P.; Callahan, S.; Schichtel, B.A.; et al. Volatile organic compounds and ozone in Rocky Mountain National Park during FRAPPÉ. Atmos. Chem. Phys. 2019, 19, 499–521. [Google Scholar] [CrossRef]
  29. Lindaas, J.; Farmer, D.K.; Pollack, I.B.; Abeleira, A.; Flocke, F.; Fisher, E.V. Acyl Peroxy Nitrates Link Oil and Natural Gas Emissions to High Ozone Abundances in the Colorado Front Range During Summer 2015. J. Geophys. Res. Atmos. 2019, 124, 2336–2350. [Google Scholar] [CrossRef]
  30. Arthur, J.D.; Langhus, B.; Alleman, D. An overview of modern shale gas development in the United States. All Consult. 2008, 3, 14–17. [Google Scholar]
  31. Le, M.T. An assessment of the potential for the development of the shale gas industry in countries outside of North America. Heliyon 2018, 4, 2. [Google Scholar] [CrossRef]
  32. Patel, S. Europe Embraces Shale Gas. Power Magazine, 1 April 2013. Available online: https://www.powermag.com/europe-embraces-shale-gas/(accessed on 26 May 2022).
  33. Shah, A.; Ricketts, H.; Pitt, J.R.; Shaw, J.T.; Kabbabe, K.; Leen, J.B.; Allen, G. Unmanned aerial vehicle observations of cold venting from exploratory hydraulic fracturing in the United Kingdom. Environ. Res. Commun. 2020, 2, 021003. [Google Scholar] [CrossRef]
  34. Państwowy Instytut Badawczy (PIG). Ocena Zasobów Wydobywalnych Gazu Ziemnego i Ropy Naftowej w Formacjach Łupkowych Dolnego Paleozoiku w Polsce (Basen Bałtycko-Podlasko-Lubelski). Raport Pierwszy; 2012. Available online: https://www.pgi.gov.pl/docman-tree-all/aktualnosci-2012/zasoby-gazu/771-raport-pl/file.html (accessed on 23 March 2018).
  35. Polish Geological Institute. Environmental Aspects of Hydraulic Fracturing Treatment Performed on the Łebień LE-2H Well. Final Report; 2011. Available online: https://www.pgi.gov.pl/en/dokumenty-pig-pib-all/kopalnia-wiedzy/gaz-lupkowy/786-the-lebien-report/file.html (accessed on 10 September 2021).
  36. Pawlak, I.; Jarosławski, J. The influence of selected meteorological parameters on the concentration of surface ozone in the central region of Poland. Atmos. Ocean. 2015, 53, 126–139. [Google Scholar] [CrossRef]
  37. Saitanis, C.J. Background ozone monitoring and phytodetection in the greater rural area of Corinth–Greece. Chemosphere 2003, 51, 913–923. [Google Scholar] [CrossRef]
  38. Nęcki, J.M.; Gałkowski, M.; Chmura, Ł.; Gerbig, C.; Zimnoch, M.; Zięba, D.; Bartyzel, J.; Wołkowicz, W.; Różański, K. Regional Representativeness of CH4 and NO2 Mixing Ratio Measurements at High-Altitude Mountain Station Kasprowy Wierch, Southern Poland. Aerosol Air Qual. Res. 2016, 16, 568–580. [Google Scholar] [CrossRef]
  39. Różański, K.; Chmura, Ł.; Gałkowski, M.; Nęcki, J.; Zimnoch, M.; Bartyzel, J. Monitoring of greenhouse gases in the atmosphere—A Polish perspective. Pap. Glob. Chang. 2016, 23, 111–126. [Google Scholar] [CrossRef]
  40. Popa, M.E.; Vollmer, M.K.; Jordan, A.; Brand, W.A.; Pathirana, S.L.; Rothe, M.; Röckmann, T. Vehicle emissions of greenhouse gases and related tracers from a tunnel study: CO:CO2, N2O:CO2, CH4:CO2, O2:CO2 ratios, and the stable isotopes 13C and 18O in CO2 and CO. Atmos. Chem. Phys. 2013, 14, 2105–2123. [Google Scholar] [CrossRef]
  41. Stanley, E.H.; Casson, N.J.; Christel, S.T.; Crawford, J.T.; Loken, L.C.; Oliver, S.K. The ecology of methane in streams and rivers: Patterns, controls, and global significance. Ecol. Monogr. 2016, 86, 146–171. [Google Scholar] [CrossRef]
  42. Skamarock, W.C.; Klemp, J.B.; Dudhie, J.; Gill, D.O.; Liu, Z.; Berner, J.; Wang, W.; Powers, J.G.; Duda, M.G.; Barker, D.M.; et al. A Description of the Advanced Research WRF Model Version 4.1; National Center for Atmospheric Research: Boulder, CO, USA, 2019. [Google Scholar] [CrossRef]
  43. Stein, A.F.; Draxler, R.R.; Rolph, G.D.; Stunder, B.J.B.; Cohen, M.D.; Ngan, F. NOAA’s HYSPLIT atmospheric transport and dispersion modeling system. Bull. Amer. Meteor. Soc. 2015, 96, 2059–2077. [Google Scholar] [CrossRef]
  44. NCEP/National Weather Service/NOAA/U.S. Department of Commerce, NCEP GFS 0.25 Degree Global Forecast Auxiliary Grids Historical Archive; Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory: Boulder, CO, USA, 2015. [CrossRef]
  45. Werner, M.; Kryza, M.; Guzikowski, J. Can Data Assimilation of Surface PM2.5 and Satellite AOD Improve WRF-Chem Forecasting? A Case Study for Two Scenarios of Particulate Air Pollution Episodes in Poland. Remote Sens. 2019, 11, 2364. [Google Scholar] [CrossRef]
Figure 1. Map of Poland with marked Pomeranian region (red square) and locations of the air quality monitoring network sites and Stary Wiec station. (Available online: https://maps.google.com (accessed on 19 July 2022)).
Figure 1. Map of Poland with marked Pomeranian region (red square) and locations of the air quality monitoring network sites and Stary Wiec station. (Available online: https://maps.google.com (accessed on 19 July 2022)).
Atmosphere 13 01228 g001
Figure 2. Monthly means of nitrogen oxide (a) and surface O3 (b) concentrations at Stary Wiec and the surrounding stations for the period 2012–2017.
Figure 2. Monthly means of nitrogen oxide (a) and surface O3 (b) concentrations at Stary Wiec and the surrounding stations for the period 2012–2017.
Atmosphere 13 01228 g002
Figure 3. Diurnal variation (mean value together with 1 standard deviation) of CO2 (a) and CH4 (b) at the Stary Wiec station, August 2015–July 2017.
Figure 3. Diurnal variation (mean value together with 1 standard deviation) of CO2 (a) and CH4 (b) at the Stary Wiec station, August 2015–July 2017.
Atmosphere 13 01228 g003
Figure 4. Concentration of CH4 and NMHC over 1st (a), 2nd (b) and 3rd (c) September 2016, measured at the Stary Wiec station.
Figure 4. Concentration of CH4 and NMHC over 1st (a), 2nd (b) and 3rd (c) September 2016, measured at the Stary Wiec station.
Atmosphere 13 01228 g004
Figure 5. Values of the CO2/CH4 ratio over 1–4 September 2016, measured at the Stary Wiec station.
Figure 5. Values of the CO2/CH4 ratio over 1–4 September 2016, measured at the Stary Wiec station.
Atmosphere 13 01228 g005
Table 1. List of shale gas research approaches.
Table 1. List of shale gas research approaches.
Author, Date of PublicationCountryPeriod of MeasurementsMeasured PollutantsMeasurement MethodMain Goal of the Study
1Oltmans et al., 2014 [25]US2011O3, NOx, NMHC-measurements from tethered balloon
-observations from a tall tower
Study of the genesis of episodes of high O3 concentration values.
2Karion et al., 2013 [22]US2012CH4-measurements from instrumented single-engine turboprop aircraftUsing mass balance method to estimate CH4 average natural leak rate.
3Shah et al., 2020 [33]UK2018–2019CH4-measurements using unmanned aerial vehicle (UAV)Detection of CH4 leak rate during first exploratory fracking process.
4Purvis et al., 2019 [8]UK2016–2018NOx, O3, SO2, NMHC, PM-in situ monitoring stationObservation of air pollutants before, during and after preparations for hydrofracturing.
5This study
Jarosławski et al., 2022
Poland2015–2017NO, NO2, NOx, O3, CO, PM10, CO2, CH4, NMHC-in situ measuring stationComparative analysis of pollutants during different stages of well’s work cycle and checking meeting of National Air Quality Criteria
Table 2. Statistics for air pollution data from Stary Wiec and nearby rural and urban stations, based on hourly means, collected from August 2015 to July 2017. Mean values for this period, together with standard deviations (1σ), are presented.
Table 2. Statistics for air pollution data from Stary Wiec and nearby rural and urban stations, based on hourly means, collected from August 2015 to July 2017. Mean values for this period, together with standard deviations (1σ), are presented.
Stary WiecLiniewkoKościerzynaTczewStarogardGdańsk
NO (ppb)0.87 ± 1.591.03 ± 0.756.05 ± 13.284.73 ± 13.607.86 ± 13.193.72 ± 9.75
NO2 (ppb)4.20 ± 2.943.15 ± 2.358.57 ± 5.417.21 ± 6.557.76 ± 8.048.74 ± 6.79
O3 (ppb)27.37 ± 13.1826.68 ± 12.7022.01 ± 12.21 22.93 ± 12.04
CO (ppb)242.5 ± 159.5316.4 ± 96.6484.9 ± 400.0285.7 ± 168.4361.4 ± 297.6323.8 ± 146.6
PM10 (µg/m3)16.18 ± 13.65 31.02 ± 28.4520.35 ± 17.3142.13 ± 44.6417.60 ± 15.45
Table 3. Statistical significance according to the t-test results (at α = 0.05) of differences between concentrations of individual air pollutants for the air coming to the measurement point directly from the wells area and from the NW sector during drilling, hydrofracturing, and the remaining period. Note: Statistically significant cases (p-value < 0.05) are bolded. Units as shown in Table 2, except CH4 (ppm) and NMHC (ppmc).
Table 3. Statistical significance according to the t-test results (at α = 0.05) of differences between concentrations of individual air pollutants for the air coming to the measurement point directly from the wells area and from the NW sector during drilling, hydrofracturing, and the remaining period. Note: Statistically significant cases (p-value < 0.05) are bolded. Units as shown in Table 2, except CH4 (ppm) and NMHC (ppmc).
DrillingHydrofracturingRemaining Period
NWWell NWWell NWWell
MeanMeanp-ValueMeanMeanp-ValueMeanMeanp-Value
NO1,411.720.2170.581.210.0440.950.830.316
NO24.185.040.0652.743.060.1703.513.630.256
CO215.5187.70.085132.96157.030.005229.08218.050.114
PM1016.0213.020.04612.3313.180.59912.9113.070.728
O320.1720.890.64727.2924.240.28129.6628.420.019
CO2426.78420.680.265406.49423.240.135412.23410.570.181
CH41.991.950.0231.931.950.1971.931.930.521
NMHC0.060.060.850.020.020.1470.030.030.216
Table 4. Air quality criteria applicable to the pollutants measured at the Stary Wiec station (1) Information threshold, (2) Alarm threshold, (3) Accumulated O3 exposure indicator (AOT40), (4) Protection of vegetation. Normal font-standards for human health protection, italics-standards for plant protection.
Table 4. Air quality criteria applicable to the pollutants measured at the Stary Wiec station (1) Information threshold, (2) Alarm threshold, (3) Accumulated O3 exposure indicator (AOT40), (4) Protection of vegetation. Normal font-standards for human health protection, italics-standards for plant protection.
PollutantAveraging PeriodAcceptable Level [ppb]Acceptable Number of Occurrences during Calendar YearNumber of Exceedances in Stary WiecAverage or Maximum Concentration Value in Stary Wiec
NO21 h104.6 ppb180 (2015–2017)
NO21 year20.9 ppb 4.3 ppb (2016)
PM1024 h50 µg/m33511 (2015–2017)
PM101 year40 µg/m3 15.3 µg/m3 (2016)
CO8 h running mean8620 ppb max 1385 ppb
O38 h running mean60 ppb25 days8 (2015)
13 (2016)
2 (2017)
O31 h90 ppb (1) 2 (2015)
O31 h120 ppb (2) 0 (2015–2017)
O3May–July9000 ppb*h (3) 7696 ppb*h (2016)
2817 ppb*h (2017)
NOx1 year15.7 ppb (4) 5.0 ppb (2016)
Table 5. Characteristic of CH4 episode events observed for 1–4 September 2016.
Table 5. Characteristic of CH4 episode events observed for 1–4 September 2016.
Episode NumberDateMax CH4 Concentration [ppm]Max NMHC Concentration [ppm]Correlation Coefficient CH4-NMHCp-ValueCH4/NMHC Ratio
11 September 20167.453.030.9980.0006.95
22 September 20162.70.390.9800.0207.52
32 September 20162.810.450.9970.0037.01
42 September 20162.320.230.9950.00010.55
52 September 20165.361.620.9990.0004.50
63 September 20162.360.200.9060.00013.14
73 September 20162.750.360.9580.00010.67
83 September 20162.300.150.6800.06415.46
93 September 20162.350.110.8590.06223.01
103 September 20162.560.310.9970.00011.46
113 September 20163.420.690.9980.0007.80
123 September 20163.470.730.9990.0007.54
133 September 20162.740.390.9990.0008.11
143 September 20162.680.350.9920.0019.36
153 September 20163.350.660.9940.0006.42
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jarosławski, J.; Pawlak, I.; Guzikowski, J.; Pietruczuk, A. Impact of Shale Gas Exploration and Exploitation Activities on the Quality of Ambient Air—The Case Study of Wysin, Poland. Atmosphere 2022, 13, 1228. https://doi.org/10.3390/atmos13081228

AMA Style

Jarosławski J, Pawlak I, Guzikowski J, Pietruczuk A. Impact of Shale Gas Exploration and Exploitation Activities on the Quality of Ambient Air—The Case Study of Wysin, Poland. Atmosphere. 2022; 13(8):1228. https://doi.org/10.3390/atmos13081228

Chicago/Turabian Style

Jarosławski, Janusz, Izabela Pawlak, Jakub Guzikowski, and Aleksander Pietruczuk. 2022. "Impact of Shale Gas Exploration and Exploitation Activities on the Quality of Ambient Air—The Case Study of Wysin, Poland" Atmosphere 13, no. 8: 1228. https://doi.org/10.3390/atmos13081228

APA Style

Jarosławski, J., Pawlak, I., Guzikowski, J., & Pietruczuk, A. (2022). Impact of Shale Gas Exploration and Exploitation Activities on the Quality of Ambient Air—The Case Study of Wysin, Poland. Atmosphere, 13(8), 1228. https://doi.org/10.3390/atmos13081228

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