Assessing Peak ‐ to ‐ Mean Ratios of Odour Intensity in the Atmosphere Near Swine Operations

: Odour in the atmosphere is usually characterized by an intermittent time series of high peaks and periods of low (or zero) concentrations. The peak ‐ to ‐ mean ratio (PMR) is commonly used to estimate short ‐ term peaks from long ‐ term averages to assess the odour impact. The objective of this study was to quantify the peak ‐ to ‐ mean ratio of odour intensity (PMR_OI) in the atmosphere near swine operations. Fifteen human assessors (sniffers) were trained to use an 8 point odour intensity scale to measure odour intensity in the ambient air near two swine operations. In each measurement session, the sniffers were placed 0  (in the direction of wind), 30  , and 45  from the wind directions at 100, 500, and 1000 m from the swine operations to sniff odour in the air every 10 s for 30 minutes. The results showed that odour in the atmosphere was intermittent. The intermittency (% of time when odour was detected) increased with the averaging time and decreased with the distance from the odour source and the direction from the wind. The measured intermittency ranged from 13% to 85%. The PMR_OI increased with the averaging time, the distance from the source, and the direction from the wind. In the wind direction, the largest difference in PMR_OI between 1 and 30 minute averaging times was 68% (2.5 vs. 4.2), which occurred at 1000 m from the odour source under stability class B. The average PMR_OI increased from 1.5 at 100 m to 3.5 at 1000 m. Atmospheric stability had a noticeable effect on PMR_OI. At 1000 m, the 30 minute PMR_OI decreased from 4.2 at stability class B (unstable) to 2.4 at E (slightly stable). Q.Z.; data curation, X.Z. and Q.Z.; Q.Z.; Q.Z. and Q.Z. and


Introduction
Odour emission from swine operations has been a major cause of confrontation between producers and the neighbouring residents. One of the most critical issues in dealing with odour is to assess or predict odour impact near swine operations, and a commonly used assessment tool is dispersion modelling. When dispersion models are used to assess the downwind odour impact, the predicted downwind concentration represents an average concentration over a period of time (typically one hour as dictated by the meteorological data). However, large fluctuations and shortterm peaks of odour in the ambient air have been consistently reported by researchers, and they are considered a characteristic feature of the atmospheric dispersion process [1]. Short-term concentrations of odorous compounds are usually characterised by an intermittent time series of high peaks and periods of low (or zero) concentrations [2]. The instantaneous peaks may be two orders of magnitude higher than the average concentration, and a few seconds of odour exposure at high levels are enough to cause human annoyance [3]. In other words, the peaks that exceed the odour detection threshold can cause odour nuisance, even if the long-term average concentration is well below the odour detection. Therefore, the short-term odour peaks should be taken into consideration in assessing odour impact.
Intermittent bursts of odour in the atmosphere near swine operations are a fundamental characteristic of pollutant dispersion in the atmosphere. Based on the dispersion of a tracer gas, Mylme and Mason [4] showed that the time series of gas concentration in the atmosphere consisted of intermittent occurrences of non-zero concentrations interspersed by periods of essentially zero concentration. The intermittency (factor), defined as the proportion of time of non-zero concentration, is often used to characterize the intermittent occurrence of pollutants in the atmosphere [4].
The peak-to-mean ratio (PMR) has been commonly used to estimate the short-term concentrations of air pollutants from the long-term averages predicted by dispersion models [5][6][7][8]. However, there are no universally accepted values for the odour PMR; different PMR values are used in different jurisdictions for odour regulations [9]. For dispersion of pollutants in the atmosphere, a one-fifth power law with sampling time for sampling periods from 3 minutes to half an hour was suggested by Stewart et al. [10]. The power law was also reported by Cramer [11] to be applicable for sampling times from 3 seconds to 10 minutes. Gifford [1] reported that for a source and receptor located at the same level, the PMR could be expected to be in the range from 1 to about 5. However, when concentrations are observed at levels considerably off the release height, or off the centre line, the PMR could be much higher than those given by the power law. Singer [5] reported that the atmospheric stability and the type of terrain also affected the PMR. Hino [12] examined the relationships between ground-level pollutant concentrations and sampling times and concluded that the concentration was proportional to the sampling time to the -0.5 power for sampling periods between 10 minutes and 5-6 hours. A power function suggested by Smith [6] is commonly used in dispersion modelling to estimate the peak concentration (Cp) during a short time interval (tp) from the mean concentration (Cm) over a long time interval (tm): , (1) where ψ = peak-to-mean ratio of odour (pollutant) concentration (PMR_OC); Cp = peak concentration; Cm = mean concentration; tp = averaging time for peak concentration (short-term); tm = averaging time for mean concentration (long-term); q = exponent (empirical constant). The exponent q in the above equation takes different values for different stability classes of the atmosphere [13]. Schauberger et al. [14] further incorporated the downwind distance into the PMR calculation using a relationship developed by Mylne and Mason [4] as follows: where ψx = peak-to-mean ratio of concentration at distance x from source; T = travel time to distance x at the mean wind speed; tl = σ/; σ = variance of wind speed averaged over three wind components;  = rate of dissipation of turbulent energy.
The PMR for odour from livestock operations is influenced not only by atmosphere dispersion but also source characteristics (e.g., frequent changes in odour emission rate). Best et al. [2] pointed out that industrial and agricultural odour PMRs depended on source characteristics, downwind distance, and atmospheric stability. Based on the review of published literature, Schauberger et al. [8] listed the following variables that influence the peak-to-mean ratio for odour: atmospheric stability; intermittency; travel time or distance from the source; lateral distance from the axis of the wake; and geometry of the source (emission height and source configuration).
Although there is a general agreement that the peak values of odour from livestock operations may be much higher than the average over a certain time period, there is a need for field data to quantify the relationship between the peak and the average odour levels. The current methods of estimating the PMR for odour are mostly based on experiments using tracer gases (e.g., [4]). The tracer methods may provide high accuracy of measurement, but the tracer release may not accurately reflect odour emission from livestock operations (i.e., the emission characteristics of livestock operations may be different from the controlled releases of tracer gases). Furthermore, odour quantification goes beyond measuring the concentration of odorants (or tracers)-it involves human perception and sensitivity.
The objective of this study was to use human sniffers (odour assessors) to determine the peakto-mean ratios of odour intensity (PMR_OI) in the atmosphere near large-scale commercial swine operations. The effects of such variables as the averaging time, downwind distance, and atmospheric stability were investigated. The intermittency of odour in the field conditions was also assessed.
It should be noted that the PMR_OI determined in this study was based on odour intensity measurement and might not be applicable directly to odour concentrations because the relationship between the odour intensity and odour concentration is not linear. For example, the Stevens power law is commonly used to relate odour intensity to concentration as follows [15]: where I = odour intensity; C = odour concentration; n = characteristic exponent; k = proportionality constant. Using the Stevens power law, the peak-to-mean ratio for odour concentration (PMR_OC) may be estimated from PMR_OI as follows: where Ip = peak odour intensity; Im = mean odour intensity.
Although it is theoretically possible to transform PMR_OI to PMR_OC, the exponent n is not readily available for livestock odour.

Description of Study Sites
Field odour measurements were conducted on two commercial swine farms of 3,000 sow farrowing operations, located in southern Manitoba, Canada. The two farms were almost identical in layout, and the only difference between the two farms was that one had an extra quarantine room, which was normally empty and did not contribute to odour emission. Each farm had 17 production rooms that were mechanically ventilated with 84 wall-mounted exhaust fans. Manure on both farms was handled as liquid, which was stored in under-floor shallow gutters and then removed to outdoor earthen manure storage once every week from gestation/breeding rooms and once every three weeks from farrowing rooms. The topography around the two farms was flat cropland.

Odour Measurement Grid
To measure odour intensity in the air downwind from the farms, 15 trained human sniffers were hired in each measurement (sniffing) session and placed in pre-determined positions based on the wind direction at the times ( Figure 1). For each sniffing session, a weather station was set up first to determine wind direction. Based on the wind direction, 15 trained sniffers were placed in a three-row grid system downwind of the farm, with the first row (5 assessors) at 100 m from the farm, the second row at 500 m, and the third row at 1000 m ( Figure 1). These distances were similar to those used by Mylme and Mason [4] in their study of concentration fluctuations of a tracer gas in the atmosphere. Three sniffers at positions 1-3, 2-3, and 3-3 were directly downwind from the farm (odour source) at 100, 500, and 1000 m, respectively. Each sniffer carried a GPS unit to guide them to the predetermined grid point and then recorded their exact position based on the longitude and latitude readings from the GPS.

Selection and Training of Field Human Odour Sniffers
Quantifying (measuring) odour in the ambient air is challenging. We followed the principle of the field inspection methods in the European Standard [16,17]-using trained human sniffers (assessors). Over 40 students were recruited from the University of Manitoba as potential human sniffers. The recruits were screened based on their sensitivity and ability to determine odour intensities of n-butanol solutions (in water). The 8 point Odour Intensity Reference Scale in the ASTM (American Society for Testing and Materials) standard (Table 1) [18] was adopted to screen and train human sniffers. According to the concentrations in Table 1, eight samples of n-butanol solutions were prepared in 45 mL glass bottles with Teflon-coated lids. In the screening test, each participant was asked to sniff the eight samples in a random order and place them in the order from the weakest to the strongest odour levels. The inversion (error) value was then calculated, and those who scored 0 or 1 were selected for further training. Of the 40 recruits, 22 passed the screening test and were selected for further training. Table 1. ASTM standard 8 point n-butanol odour intensity reference scale [18].

N-butanol in water (ppm) 0 120 240 480 960 1940 3880 7750 15500
Intensity level 0 1 2 3 4 5 6 7 8 A series of training sessions were conducted to train the selected recruits, and each training session consisted of three phases. In phase 1, every sniffer was provided with a set of eight n-butanol samples at the concentrations described in Table 1 and asked to sniff the samples one by one from #1 (0 ppm) to #8 (15500 ppm) repeatedly until they felt confident that they had "memorized" the concentration scale. To minimize the effect of any odour in the ambient air, the sniffers wore a carbon filtered mask for 10-20 s before each sniffing. In phase 2, each sniffer was given 3 to 6 coded samples of certain concentrations (but unknown to the sniffer) and asked to evaluate these samples one at a time and assign a scale (0 to 8) to each sample. Those who correctly rated all samples were asked to check with the standard solution bottles of n-butanol to reinforce their ratings. Sniffers who incorrectly rated a sample(s) had to sniff both the standard and the coded samples to "feel" the difference until she/he could correctly rate all samples. In phase 3, three samples of swine manure odour simulant [19] were presented to each sniffer for assessment one at a time. Each sniffer assigned a scale (0-8) to each of the three samples of simulated swine manure odour, and the process was repeated until a group consensus was reached for each of the three samples. The results were collected after each training session to evaluate the performance of sniffers over the entire training period. Each sniffer was allowed to be one level off for the wrong identification; otherwise further training had to be conducted (repeated).

Field Odour Intensity Measurement Procedure
At the beginning of each odour measurement session (before leaving for field), the sniffers "checked" their noses by sniffing the eight standard n-butanol solutions and recorded their ratings on a labelled magnitude scale sheet [20], which was kept to track the history of each sniffer's performance in each field session. After reaching the field measurement positions as illustrated in Figure 1, the sniffers waited for instructions from a central coordinator through a two-way radio to start a sniffing session. Following a similar procedure to that outlined by Mannebeck et al. [21], the sniffers inhaled the ambient air to assess (sniff) odour every 10 seconds and repeated sniffing for 10 minutes. Three 10 minute sessions were conducted consecutively with a 10 minute break between each session. The timing of each sniffing was synchronized by a central coordinator, i.e., through the two radios, the coordinator broadcast the instruction to the sniffers every 10 s to start and stop sniffing. To prevent the nose from being "saturated" by odour in the air, every sniffer wore a carbon filtered mask and only took off the mask briefly every 10 s to sniff odour. Each sniffer recorded her/his sniffing result (an intensity value of between 0 and 8) every 10 s and a brief description of odour characteristics on a data recording sheet.
In total, 51 odour sniffing sessions were carried out on the two farms between June and October, when swine odour was considered to be a problem for the local residents. Two or three sessions were normally conducted in one day to cover different times of the day, specifically, from early morning to noon or from afternoon to early evening. The selected measurement dates and times covered a range of atmospheric stability conditions (Table 2). However, the stability class was not known until the collected weather data was analysed after each field session. Although some sessions were conducted in the early morning (before sunrise) with the intention of capturing stable atmospheric conditions, no stability class F was measured.

Meteorological Condition Monitoring
A portable weather station (WatchDog Model 550, Spectrum Technologies, Inc., Plainfield, IL) was used to record the on-site meteorological conditions, including solar radiation, temperature, relative humidity, and wind speed and direction. The data was recorded every minute. Based on the wind direction at the time, an open area on the farm was selected to place the weather station at a height of 2 m above the ground level.

Data Analysis
The odour intensity measurement method used in this study was based on magnitude estimation scaling. Therefore, the measured intensity data could be analysed statistically by using parametric methods, such as computation of means, variance analysis (ANOVA), and regression [22]. Furthermore, the mean of intensity values was calculated as the arithmetic mean [23][24][25]. Minitab (Minitab 19, Minitab Inc., State College, PA) was used to perform the statistical analyses.
A sample of field measurement data for a 5 minute period is presented in Figure 2 to illustrate the calculations of PMR_OI and intermittency (). This sample data had a total of 31 intensity values recorded by a sniffer, with a maximum of 6.0 and an arithmetic mean of 1.6 (Im). In the context of odour dispersion modelling, the peak value is estimated from a long-term average value by using empirical equations, such as Equation (1). This estimated peak value is commonly used to define the acceptable odour level. However, researchers are still debating the definition of the peak value [8]. In the current study, the peak intensity value (Ip) was defined as the 90th percentile of measured odour intensity values during the time period of interest [8], which was calculated to be Ip = 4.0 in this example. Therefore, the PMR_OI in this example was calculated to be PMR_OI = Ip/Im = 4.0/1.6 = 2.5. In addition, 14 zero values (no odour) were reported, based on which the intermittency was calculated to be  = (1 − 14/31) = 0.548 (or 54.8%). The conditioned mean odour intensity, defined as the mean of the non-zero readings [4], was calculated to be 2.9 and the 90th percentile 5.0. Therefore, the corresponding conditioned peak-to-mean ratio (CPMR_OI) was determined to be 5.0/2.9 = 1.7. To improve the data quality, retrospective screening was performed on the datasets in the following three steps: (1) the mean and standard deviation of the data at each of the 15 measurement positions (1-1 to 3-5) were calculated for the 30 minute session; (2) any data points that were outside three standard deviations were removed; (3) the means and standard deviations were calculated again for the screened dataset, and the process was repeated until no change in the mean was observed.

Preliminary Test
Visual cues could modulate olfactory perceptions [26]. A preliminary test was conducted to check if seeing the odour source (farm) might bias the odour rating by human sniffers. Two odour sniffing sessions were conducted, one upwind and one downwind of a farm. It was expected that the sniffers would report odour only when they were downwind from the farm if they were not biased. Of the 2745 data points reported by 15 sniffers, 93% were at level 0 (no odour), 5% at level 1, and about 1% each at levels 2 and 3, when the sniffers were upwind from the odour source ( Figure 3). The sniffers reported much higher odour intensity levels when they were downwind from the farm. Specifically, the frequencies of reporting intensity levels of 4, 5, and 6 were 14%, 17%, and 17%, respectively ( Figure 3). The average reported odour intensity was close to zero (0.1) for upwind and 2.9 for downwind. The above two sets of data confirmed that seeing the odour source did not bias the odour assessments by human sniffers.

Intermittency
One of the important characteristics of odour dispersion in the atmosphere is "intermittency", i.e., it occurs in bursts with low or no odour between bursts. It should be noted that the intermittency of odour in the atmosphere around swine operations is dependent on both the atmospheric conditions and the variability of odour emission (e.g., changes in ventilation in swine barns). This intermittent characteristic of odour in the atmosphere was clearly shown by the measured data ( Figure 4). It was observed that more zeros (lower intermittency) were reported as the distance from the odour source increased. The calculated intermittency values based on the data from all 51 field measurement sessions are presented in Figure 5 for different distances and average times along the wind direction.  Two-way ANOVA indicated that both the averaging time and the downwind distance had a significant (P<0.05) effect on the intermittency. The intermittency generally decreased with the distance from the odour source, and the largest decreases were observed for the averaging time of 1 minute. Specifically, the intermittency was 75%, 42%, and 17% at 100, 500, and 1000 m, respectively, for a 1 minute averaging time. This trend agrees with the observation reported by Mylne and Mason [4]. In comparison, the changes in intermittency with distance were less pronounced for a 30 minute averaging time (85%, 87%, and 57% at 100, 500, and 1000 m, respectively).
The intermittency increased with the averaging time, and the increases were more pronounced at longer distances. At 100 m, the average intermittency value increased from 75% to 80%, 83%, and 85% as the averaging time increased from 1 to 5, 10, and 30 minutes, respectively. In comparison, the increases in intermittency at 1000 m were more pronounced: 17%, 33%, 50%, and 57% for the four averaging times, respectively, or an increase of 3.4 times when the averaging time increased from 1 to 30 minutes.
The intermittency decreased with the angle from the wind direction, i.e., the highest intermittency occurred along the wind direction (0, Table 3). The differences in intermittency between the 0 and 45 angles were larger at longer distances from the odour source. The intermittency at 45 was on average 8% lower than that at 0 at 100 m, 17% at 500 m, and 30% at 1000 m. ANOVA indicated that the difference in intermittency between 0 and 45 was significant (P<0.05) at 500 and 1000 m, but not significant at 100 m. This was because odour was relatively strong and less intermittent at locations close to the source.

PMR_OI versus CPMR_OI
Two definitions were used to calculate peak-to-mean ratios of odour intensity: the commonly used "true" PMR_OI that included all data points, and the less commonly used conditioned peak-tomean ratio (CPMR_OI) based on non-zero odour readings only. While this paper was focused on the PMR_OI, the CPMR_OI was used for comparison to illustrate the effect of intermittency on the calculation of the peak-to-mean ratio.
The measured PMR_OI varied markedly within a measurement session (with time and location) and from session to session, with values ranging from 1 to 9 (Figure 6a). A wider spread of PMR_OI was observed at longer distance from the source. For example, the PMR_OI was mostly between 1 and 2 at 100 m from the source, while the PMR_OI ranged from 2 to 9 at 1000 m. In comparison, the highest CPMR_OI was 6, and most CPMR_OI values were between 1 and 2 for all three distances (Figure 6b). This clearly indicated that zero odour readings between bursts (i.e., intermittency) had a drastic effect on the PMR_OI because the zero readings lowered the average odour intensity, thus resulting in a higher PMR_OI. The difference between the PMR_OI and CPMR_OI was statistically significant (P<0.05) at 500 and 1000 m, but not at 100 m (Figure 7). For the 30 minute averaging time, the PMR_OI was 1.8 and 2.3 times higher than the CPMR_OI for 500 and 1000 m, respectively.
(a) (b) Figure 6. Frequency distribution of the peak-to-mean ratio of odour intensity along the wind direction at 100, 500, and 1000 m from the odour source. (a) PMR_OI (true peak-to-mean ratio of odour intensity) and (b) CPMR_OI (conditioned peak-to-mean ratio).

Figure 7.
Comparison of measured PMR_OI (peak-to-mean ratio of odour intensity) and CPMR_OI (conditioned peak-to-mean ratio of odour intensity) (30 minute averaging time) along the wind direction at 100, 500, and 1000 m from odour source (bar: standard deviation).

Variation of PMR_OI with Distance to the Odour Source
It was observed that the PMR_OI increased with the distance from the odour source, while the CPMR_OI changed little with the distance (Figure 7). For example, the PMR_OI for a 30 minute averaging time increased from 1.5 to 3.1 and 3.6 as the distance increased from 100 to 500 and 1000 m. ANOVA indicated that the effect of distance from the odour source on PMR_OI was statistically significant (P<0.05). Comparing Figure 5 and Figure 7 indicated that the intermittency had a dominant effect on the PMR_OI. Specifically, decreasing intermittency at longer distances could result in a lower mean odour intensity and higher PMR_OI. To support this argument, the following set of simplified equations was proposed to estimate the PMR_OI: where PMR_OI * = peak-to-mean ratio of odour intensity as a function of intermittency; Ii = a reading of odour intensity at a time point; N = total number of readings within a time period. Based on the definition of PMR_OI * , all zero Ii values are included in the calculation. If the zeros are excluded, the conditioned mean Im,c is obtained. Given that the difference between the "true" odour intensity mean Im and the conditioned mean Im,c is due to the intermittency, it is reasonable to assume that the mean was proportional to the intermittency and the conditioned mean, that is, Im =Im,c. The PMR_OI * is now calculated as where Im,c = conditioned mean intensity (the mean of non-zero readings);  = intermittency;  = a proportional constant.
Given that both Ip and Im,c are independent of intermittency, it was further assumed that they both follow a similar exponential decay pattern in the atmosphere as the distance to the source increases (due to dilution): where = peak odour intensity (constant) at x=0; , = conditioned mean odour intensity (constant) at x=0; x = distance from the odour source; a = constant. Substituting Equation (7) into (6) yields: The above equation indicates that PMR_OI * increases as the intermittency decreases. The intermittency  decreases with distance x (Table 3), and therefore, PMR_OI * increases with distance. It should be noted that at longer distances, the actual odour levels were much lower even though the PMR_OI was higher. Taking the measurement session of August 5 as an example (the measurements of other dates had similar patterns), the PMR_OI was 1.8, 3.1, and 6.1 for 100, 500, and 1000 m, respectively, while the measured peak odour intensity was 6.0, 2.0, and 1.0 at the three distances, respectively ( Figure 8).

Effect of Wind Direction on PMR_OI
In the off-wind directions, more zero odour readings were generally recorded, and consequently, higher PMR_OIs were observed ( Figure 9). Two-way ANOVA, considering the effects of the angle from the wind direction and the distance from the odour source, indicated that the angle effect on the PMR_OI was not significant (P>0.05) at 500 m but was significant at 100 and 1000 m. The largest difference in PMR_OI between the 0 and 45 angles was about 17% (3.6 vs. 4.2) at 1000 m. Figure 9. Comparison of measured PMR_OIs (30 minute averaging time) between 0 (along the wind direction) and 45 from wind direction at 100, 500, and 1000 m from the odour source.

Effect of Atmospheric Stability on PMR_OI
Of the 51 measurement sessions, there were only two sessions at atmospheric stability class A, which did not generate sufficient data to include stability class A in the statistical analysis. Therefore, only the data for atmospheric stability classes B, C, D, and E were analysed. Two-way ANOVA indicated that the effect of atmospheric stability on the PRM_OI was not significant (P>0.05) at 100 and 500 m and was significant at 1000 m (Figure 10). At 1000 m, the more turbulent (unstable) the atmosphere was, the higher the PMR_OI. Specifically, the PMR_OI decreased from 4.2 to 3.2, 3.1, and 2.4 as the stability changed from B (unstable) to, C, D, and E (slightly stable).

Effect of Averaging Time on PMR_OI
The measured PMR_OI increased with the averaging time, and the largest increase was observed at 1000 m under stability class B: from 2.5 to 4.2 when the averaging time increased from 1 to 30 minutes. The increase of the PMR_OI with the averaging time followed a non-linear pattern ( Figure  11). Therefore, the same power function as defined in Equation 1, i.e., PMR_OI = (tm/tp) q , was used to predict the PMR_OI from the time ratio (tm/tp). To evaluate the exponent q in the equation, the PMR_OI was plotted against (tm/tp) in the log-log scale and the exponent q was determined as the slope of the regression line ( Table 4). The exponent q decreased from 0.33 for stability class B (unstable) to 0.19 for E (slightly stable). The measured exponent values were lower than those reported in the literature for unstable atmosphere conditions. Smith [13] gave the following q-values for three stability classes: 0.65 (stability class B), 0.52 (C), and 0.35 (D). Duffee et al. [27] reported values of q being 0.50 for stability class A or B, 0.33 for C, 0.20 for D, and 0.17 for E and F.

Conclusions
Odour in the atmosphere near (within 1000 m) swine operations was intermittent. The intermittency, defined as the proportion of time when odour was detected, generally decreased with the distance from the odour source, and the decreases were more pronounced for the shorter averaging time of 1 minute. The intermittency also increased with the averaging time, in particular at longer distances.
The peak-to-mean ratio of odour intensity (PMR_OI) in the atmosphere increased with the distance from the odour source. Along the wind direction, the 30 minute PMR_OI increased from 1.5 at 100 m to 3.6 at 1000 m. The higher PMR_OI at longer distance from the odour source was primarily due to the lower intermittency at longer distances.
The effect of atmospheric stability on PMR_OI was not apparent within 500 m from the odour source. At 1000 m, the PMR_OI decreased as the atmospheric stability changed from unstable (B) to slightly stable (E) conditions.
The PMR_OI increased with the averaging time in a non-linear pattern. The increase of the PMR_OI with the peak-to-mean time ratio could be reasonably predicted by the power law, with the exponent decreased from 0.33 to 0.19 as the atmospheric stability class changed from stability class B to E.