Obukhov Length Research Articles

Scaling of vertical coherence and logarithmic energy profile for wall-attached eddies during sand and dust storms

High-frequency observation data, including all three components of instantaneous fluctuating velocity, temperature, as well as particulate matter 10 ( $PM_{10}$ ), collected from the unstable atmospheric surface layer at $z/L = -0.11$ and $-$ 0.12, $L$ being the Obukhov length, during sand and dust storms (SDS), were used to explore the scaling of vertical coherence and the logarithmic energy profile for wall-attached eddies. The present results demonstrate good agreement with the self-similar range of the wall-attached features for velocity and temperature components, as well as for $PM_{10}$ at lower heights ( $z<15$ m) during SDS. Following the idea depicted by Davenport (Q. J. R. Meteorol., vol. 372, 1961, pp. 194–211), an empirically derived transfer kernel comprises implicit filtering via a scale-dependent gain and phase, parametrically defined as $|H_L^2(f)|=\exp (c_1-c_2\delta /\lambda _x)$ , where $c_1$ and $c_2$ are parameters, $\delta$ is the boundary layer thickness and $\lambda _x$ is the streamwise wavelength. Linear coherence spectrum analysis is applied as a filter to separate the coherent and incoherent portions. After this separation procedure, the turbulence intensity decay for wall-attached eddies is described in a log–linear manner, which also identifies how the scaling parameter differs between the measured components. These findings present abundant features of wall-attached eddies during SDS which further are used to improve/enrich existing near-wall models.

Influence of atmospheric stability conditions on wind energy density in Ali Al-Gharbi region of Iraq

Atmospheric stability is considered as one of the most important factors affecting the increase or decrease in wind speed in the atmosphere and thus affects the wind energy density. This study aims to analyze the stability of the atmospheric conditions in southern part of Iraq, specifically in the Ali Al-Gharbi region, using one of methods to determine atmospheric stability called Monin-Obukhov length . Field data of horizontal and vertical wind speed and air temperature collected at three heights (10m, 30m, and 50m) in 2017 from tower installed in Ali Al-Gharbi region. The results show that the vertical wind speed increases with height up to 30 m and decreases at 50 m, while the horizontal wind speed increases as the height increases. The friction velocity , surface heat flux , momentum flux and shear stress were calculated and the stability conditions were determined. Results revealed that stable atmospheric condition was the most frequent with about 59% occasions occurring during the year followed by unstable (40%) and neutral (1%) conditions. The highest wind energy density was in stable conditions (0 ≤ L< 200), with percentage (58% - 57%) in spring season at both heights 10m and 50m, respectively, followed by unstable conditions (-200 < L ≤ 0), while the lowest wind energy density was in neutral conditions (-200 ≤ L ≥ 200) with percentage (1.3% - 0.8%) in autumn season.Atmospheric stability is considered as one of the most important factors affecting the increase or decrease in wind speed in the atmosphere and thus affects the wind energy density. This study aims to analyze the stability of the atmospheric conditions in southern part of Iraq, specifically in the Ali Al-Gharbi region, using one of methods to determine atmospheric stability called Monin-Obukhov length . Field data of horizontal and vertical wind speed and air temperature collected at three heights (10m, 30m, and 50m) in 2017 from tower installed in Ali Al-Gharbi region. The results show that the vertical wind speed increases with height up to 30 m and decreases at 50 m, while the horizontal wind speed increases as the height increases. The friction velocity , surface heat flux , momentum flux and shear stress were calculated and the stability conditions were determined. Results revealed that stable atmospheric condition was the most frequent with about 59% occasions occurring during the year followed by unstable (40%) and neutral (1%) conditions. The highest wind energy density was in stable conditions (0 ≤ L< 200), with percentage (58% - 57%) in spring season at both heights 10m and 50m, respectively, followed by unstable conditions (-200 < L ≤ 0), while the lowest wind energy density was in neutral conditions (-200 ≤ L ≥ 200) with percentage (1.3% - 0.8%) in autumn season.

The coupling of winds, ocean turbulence, and High Salinity Shelf Water in the Terra Nova Bay Polynya

The Terra Nova Bay (TNB) Polynya in the Western Ross Sea of Antarctica is a major producer of High Salinity Shelf Water (HSSW), a precursor to Antarctica Bottom Water (AABW). Processes occurring in and around the polynya can therefore effect change in the lower limb of overturning circulation in this region. Here, we use data from a densely-instrumented upper-ocean mooring, deployed for 1 year in a region of active HSSW formation within TNB, to examine the coupling of surface brine rejection and vertical mixing to katabatic wind forcing. We find a high correlation between salinity and winds during the wintertime HSSW production season at the mooring site, with a lag-response of 20 h in near-surface (∼47 m) salinity to winds measured at the nearby Automatic Weather Station (AWS) Manuela. Salinity and temperature measurements show a fully destratified water column by June, with a lag-response of near-seabed (∼360 m) salinity to near-surface salinity of just 5 h. Measurements of turbulent kinetic energy (TKE) dissipation rate (ɛ) from moored pulse-coherent acoustic Doppler current profilers (ADCPs) show general agreement with classic boundary layer scaling (BLS), and calculations of a vertical mixing timescale using the Obukhov length scale average to ∼2.5 h during austral winter, consistent with the 5-h lag time observed in the salinity data. Comparisons to data from concurrent mooring deployments along the southern boundary of TNB, as well as to previously published assessments of model simulations and data from Climatic Long-term Interaction for the Mass-balance in Antarctica (CLIMA) moorings, allow us to explore spatial variability in the coupling of winds and salinity across TNB and to speculate on possible HSSW circulation pathways.

Assessment of Wind over Complex Terrain Considering the Effects of Topography, Atmospheric Stability and Turbine Wakes

This study proposes a microscale flow model to estimate mean wind speed, fluctuating wind speed and wind direction over complex terrain considering the effects of topography, atmospheric stability, and turbine wakes. Firstly, the effect of topography is considered using Computational Fluid Dynamics (CFD). Next, a mesoscale model is presented to account for the effect of atmospheric stability. The effect of turbine wakes on the mean and fluctuating wind speeds are then represented by an advanced wake model. The model is validated using the measurement data of a wind farm located in the North of Japan. The measured wind data by Lidar at a reference height are horizontally extrapolated to a nearby met mast hub height and validated by a cup anemometer. Moreover, a novel averaging method is proposed to calculate a directional equivalent Monin–Obukhov length scale to account for the effect of atmospheric stability. Finally, the measured wind data at the reference height are vertically extrapolated and validated at the lidar location. The predicted mean and fluctuating wind speeds show good agreement with the measurements.

A numerical analysis of wind farm wake characteristics in the southern part of the North Sea

The rapid expansion of wind farms in the North Sea requires a better understanding of the wind turbines’ wake effects. In this study, the classification of wake patterns under different atmospheric boundary layer stability conditions is investigated. For this purpose, the Weather Research and Forecast model is utilized to calculate wind farm wake effects over the southern North Sea for a year-long period. The atmospheric stability is characterized by the value of Monin-Obukhov length, and seven different classes are considered. The results have shown that the predominance of the stability condition depends on seasonality. In autumn and winter months very unstable conditions prevail, while in spring and summer periods near-neutral and stable events occur more frequently. The different atmospheric stability conditions have distinct effects on the averaged wind-speed deficits. Specifically, in near-neutral and stable stratification, wakes propagate further downwind of the wind turbines affecting neighboring wind farms, while in the case of unstable conditions, these effects are weaker.

Wind Shear Model Considering Atmospheric Stability to Improve Accuracy of Wind Resource Assessment

An accurate wind shear model is an important prerequisite in extrapolating the wind resource from lower heights to the increasing hub height of wind turbines. Based on the 1-year dataset (collected in 2014) consisting of 15-minute intervals collected at heights of 2, 10, 50, 100, and 150 m on an anemometer tower in northern China, the present study focuses on the time-varying relationship between the wind shear coefficient (WSC) and atmospheric stability and proposes a wind shear model considering atmospheric stability. Through the relationship between Monin–Obukhov (M-O) length and gradient Richardson number, the M-O length is directly calculated by wind data, and the WSC is calculated by combining the Panofsky and Dutton (PD) models, which enhances the engineering practicability of the model. Then, the performance of the model is quantified and compared with two alternative methods: the use of annual average WSC and the use of stability change WSC extrapolation. The analysis demonstrates that the proposed model outperforms the other approaches in terms of normal root mean square error (NRMSE) and normal bias (NB). More specifically, this method reduces the NRMSE and NB by 24–29% and 76–95%, respectively. Meanwhile, it reaches the highest extrapolation accuracy under unstable and stable atmospheric conditions. The results are verified using the Weibull distribution.

Small-scale variation of atmospheric dynamics applying chaos theory, case study

Characterization and knowledge of the variability of atmospheric dynamics on a small scale in the city of Riobamba, Ecuador, are achieved through the chaos theory. Meteorological data is taken every hour during four years, including variables such as wind speed, wind direction, incident radiation, temperature, and humidity, from the ESPOCH, SAN JUAN, and QUIMIAG weather stations in the canton of Riobamba. The van Ulden and Hostlang models are used to calculate the Obukhov length, surface heat fluxes, and latent heat flux. The chaos theory is applied to study the variation of atmospheric microdynamics. The Lyapunov coefficients, Kolmogorov-Sinai entropy, and Kaplan-Yorke fractal dimension are determined. Before analysis, noise reduction is necessary due to the lack of correlation, especially in the Obukhov length. This research follows a longitudinal design and employs quantitative and explanatory methods based on data analysis, statistical-mathematical techniques, and inductive-deductive approaches. The results indicate a highly variable system, reflected in a high number of Lyapunov coefficients, fractional dimensions, and entropy variations. The microdynamic parameters exhibit hyperchaotic behavior, as indicated by the presence of more than one positive Lyapunov coefficient. The variables also demonstrate a fractional fractal dimension, highlighting the irregularity in the geometric representation of the system.

Addressing Methane Venting: Strategies and Implications on the Environment

Industrial processes, including oil and gas production, landfill operations, and agriculture, emit methane, a potent greenhouse gas. This article examined methane emissions in Rivers State's Soku, Agbada, and Oyigbo oil and gas extraction communities. The net methane emissions(mg/l) for facilities A, B, and C are 0.90, 0.28, and 1.03, respectively. The corresponding temperature rise for the host communities over the same period were 1.87 oC, 0.37 oC, and 1.16 oC, respectively. All three sites have near-neutral unstable methane dispersion using Monin-Obukhov length. This article examined new methane venting mitigation technologies and their effects on bioenergy and environmental engineering. The study analyzed existing methane venting trends in several Rivers State oil production sites and their effects on local temperature and the environment, emphasizing the necessity for adequate mitigation. This paper also examined new methane capture, use, and sequestration technologies. The article also examines methane mitigation and bioenergy industry synergies, stressing co-benefit strategies that improve sustainability and economic feasibility. Methane mitigation strategies in Rivers State increase energy efficiency, carbon footprint, and air quality, according to the paper. This study helps reduce methane emissions by providing possible alternatives and their effects on bioenergy development and environmental stewardship. Interdisciplinarity and innovation are encouraged to accelerate sustainable and resilient development.

Three dimensions advection-diffusion equation in unstable condition

The advection-Diffusion Equation in three dimensions was solved using the Separation variable, Substitution and Laplace transform methods, taking into account that eddy diffusivity along the horizontal, lateral coordinates and wind velocity are constant respectively, while the vertical eddy diffusivity is a function of the vertical height "z", Monin- Obukhov length "L" friction velocity at different emission rate and operation of iodine (I-131) in an unstable condition. We took the measured data from Inshas, Cairo, Egypt and compared the concentration of the calculated results with the results obtained from the solution equation at heights of 0.7, 27, and 43 meters. We found that all models are within a factor of two with the observed data from statistical evaluations in unstable conditions, we found that all models fell into a factor of two with the observed data. Whereas with NMSE and FB, the current models correlate better with the observed data.

Symmetry Analysis of Mean Velocity Distribution in Stratified Atmospheric Surface Layers

The mean velocity distributions of unstably and stably stratified atmospheric surface layers (ASLs) are investigated here using the symmetry approach. Symmetry groups for the mean momentum and the Reynolds stress equations of ASL are searched under random dilation transformations, which, with different leading order balances in different flow regions, lead to a set of specific scalings for the characteristic length ℓ13 (defined by Reynolds shear stress and mean shear). In particular, symmetry analysis shows that in the shear-dominated region, ℓ13 scales linearly with the surface height z, which corresponds to the classical log law of mean velocity. In the buoyancy-dominated region, ℓ13/L∼z/L4/3 for unstably stratified ASL and ℓ13/L∼const for stably stratified ASL, where L is the Obukhov length. The specific formula of the celebrated Monin–Obukhov similarity function is obtained, and hence an algebraic model of mean velocity profiles in ASL is derived, showing good agreement with the datum from the QingTu Lake observation array (QLOA) in China.

Obukhov Length Estimation From Spaceborne Radars

AbstractTwo air‐sea interaction quantification methods are employed on synthetic aperture radar (SAR) scenes containing atmospheric‐turbulence signatures. Quantification performance is assessed on Obukhov length L, an atmospheric surface‐layer stability metric. The first method correlates spectral energy at specific turbulence‐spectrum wavelengths directly to L. Improved results are obtained from the second method, which relies on a machine‐learning algorithm trained on a wider array of SAR‐derived parameters. When applied on scenes containing convective signatures, the second method is able to predict approximately 80% of observed variance with respect to validation. Estimated wind speed provides the bulk of predictive power while parameters related to the kilometer‐scale distribution of spectral energy contribute to a significant reduction in prediction errors, enabling the methodology to be applied on a scene‐by‐scene basis. Differences between these physically based estimates and parameterized numerical models may guide the latter's improvement.

Stratification and temporal evolution of mixing regimes in diurnally heated river flows

Direct numerical simulations of stratified open channel flows subject to a varying surface heat flux are performed. The influence of the diurnal heating time on the spatial and temporal variation of mixing in the flow and the characteristics of the mean flow state are examined. The control parameters are the bulk stability parameter lambda_{B}, defined through the ratio of the channel height delta and a bulk Obukhov length scale mathscr{L}_{B}, and the diurnal time scale hat{t}, defined as the ratio of the heating time to an eddy turnover time. The Prandtl number Pr and Reynolds number Re_{tau} have values of 1 and 400. Simulations are performed over hat{t} = 1 to 24 and lambda_{B} = 0.6 to 26. Two key flow features are used to classify the flow regimes observed, namely the laminar layer depth (LLD) and stratified layer depth (SLD) where the LLD is defined as the depth from the free surface when the buoyancy Reynolds number Re_{B} approx 7 and the SLD is the depth from the free surface when the turbulent Froude number Fr approx 1. This study attempts to characterise how these length scales vary across the diel cycle. The LLD is a viscous length scale and a regime map of a viscous parameter, the bulk Obhukov Reynolds number Re_mathscr{L}, and hat{t} is presented to classify the LLD behaviour. A regime map of lambda _{B} and hat{t} is presented to classify the behaviour of the SLD. Three classifications for each layer depth behaviour within a diel cycle form the basis of the regime maps for this paper: a neutral flow where the LLD or SLD does not exist (denoted by NL and NS), a stratified flow where the LLD or SLD are diurnally varying (denoted as DL and DS) and a persistent layer of the LLD or SLD (denoted as PL and PS). The transition between the NL to DL is hat{t} propto Re_{mathscr{L}}^{4.5}, DL to PL is hat{t} propto Re_{mathscr{L}}^{- 0.5}, NS to DS is hat{t} propto lambda _{B}^{0} and DS to PS is hat{t} propto lambda _{B}^{1}. The regime maps may be used as a predictive tool to determine when suppressed mixing regimes occur in rivers. At each flow depth, the flow sweeps though a range of mixing states across the diel cycle. The local mixing efficiency are briefly assessed and found to scale well with the instantaneous Fr number according to the regimes proposed by Garanaik and Venayagamoorthy (J. Fluid Mech., vol. 867, 2019, pp. 323-333).

Evidence of Mixed Scaling for Mean Profile Similarity in the Stable Atmospheric Surface Layer

Abstract A new mixed scaling parameter Z = z/(Lh)1/2 is proposed for similarity in the stable atmospheric surface layer, where z is the height, L is the Obukhov length, and h is the boundary layer depth. In comparison with the parameter ζ = z/L from Monin–Obukhov similarity theory (MOST), the new parameter Z leads to improved mean profile similarity for wind speed and air temperature in large-eddy simulations. It also yields the same linear similarity relation for CASES-99 field measurements, including in the strongly stable (but still turbulent) regime where large deviations from MOST are observed. Results further suggest that similarity for turbulent energy dissipation rate depends on both Z and ζ. The proposed mixed scaling of Z and relevance of h can be explained by physical arguments related to the limit of z-less stratification that is reached asymptotically above the surface layer. The presented evidence and fitted similarity relations are promising, but the results and arguments are limited to a small sample of idealized stationary stable boundary layers. Corroboration is needed from independent datasets and analyses, including for complex and transient conditions not tested here.

On the Retrieval of Surface-Layer Parameters from Lidar Wind-Profile Measurements

We revisit two recent methodologies based on Monin–Obukhov Similarity Theory (MOST), the 2D method and Hybrid-Wind (HW), which are aimed at estimation of the Obukhov length, friction velocity and kinematic heat flux within the surface layer. Both methods use wind-speed profile measurements only and their comparative performance requires assessment. Synthetic and observational data are used for their quantitative assessment. We also present a procedure to generate synthetic noise-corrupted wind profiles based on estimation of the probability density functions for MOST-related variables (e.g., friction velocity) and the statistics of the noise-corrupting perturbational amplitude found during an 82-day IJmuiden observational campaign. In the observational part of the study, 2D and HW parameter retrievals from floating Doppler wind lidar measurements are compared against those from a reference mast. Overall, the 2D algorithm outperformed the HW in the estimation of all the three parameters above. For instance, when assessing the friction-velocity retrieval performance with reference to sonic anemometers, determination coefficients of ρ2D2=0.77 and ρHW2=0.33 were found under unstable atmospheric stability conditions, and ρ2D2=0.81 and ρHW2=0.07 under stable conditions, which suggests the 2D algorithm as a prominent method for estimating the above-mentioned surface-layer parameters.

Influence of Atmospheric Stability on Wind Turbine Energy Production: A Case Study of the Coastal Region of Yucatan

Wind energy production mainly depends on atmospheric conditions. The atmospheric stability can be described through different parameters, such as wind shear, turbulence intensity, bulk Richardson number, and the Monin–Obukhov length. Although they are frequently used in micrometeorology and the wind industry, there is no standard comparison method. This study describes the atmospheric stability of a coastal region of Yucatan, Mexico, using these four parameters. They are calculated using six-month data from a meteorological mast and a marine buoy to determine atmospheric stability conditions and compare their results. The unstable atmospheric condition was predominant at the site, with an 80% occurrence during the measurement period, followed by 12% in neutral and 6% in stable conditions. Wind speed estimations were performed for each atmospheric stability scenario, and the variation in the energy produced was derived for each case. Unstable atmospheric conditions deliver up to 8% more power than stable conditions, while neutral conditions deliver up to 9% more energy than stable conditions. Therefore, considering a neutral state may lead to a considerably biased energy production estimation. Finally, an example calculation indicates that atmospheric stability is a crucial parameter in estimating wind energy production more accurately.

Intermittency and critical mixing in internally heated stratified channel flow

Through direct numerical simulations we investigate the effects of spatiotemporal intermittency as a result of stable stratification in surface heated stratified open channel flow. By adapting the density inversion criterion method of Portwood et al. [J. Fluid Mech., vol. 807, 2016, R2] for our flow, we demonstrate that the flow may be robustly separated into regions of active turbulence for which $Re_B \gtrsim {O}(10)$ and surrounding quiescent fluid, where $Re_B$ is the buoyancy Reynolds number. The intermittency in the flow spontaneously manifests as a deformed horizontal interface between the upper quiescent and lower turbulent flow, characterised by vigorous mixing from ‘overturning’ shear instabilities. The resulting vertical intermittency profile is accurately predicted by a local Monin–Obukhov length normalised by viscous wall units $\varLambda ^+$ such that the flow displays intermittency within the parameter range of $2.5 \lesssim \varLambda ^+ \lesssim 260$ . By considering conditional averages of the ‘turbulent’ and ‘quiescent’ flow separately, we find the ‘turbulent’ flow within this region to be described by constant critical gradient Richardson and turbulent Froude numbers of $Ri_{g,c} \approx 0.2$ and $Fr_c \approx 0.3$ . We find that the turbulent flow continues to display a $\varGamma \sim Fr^{-1}$ relationship when $Fr < Fr_c$ , whereas the quiescent flow shows no correlation between $\varGamma$ and $Fr$ , where $\varGamma$ is the flux coefficient. Hence, we demonstrate directly that for our flow, the emergence of an asymptotic ‘saturated’ $\varGamma$ regime in the limit of a low ‘global’ $Fr$ occurs due to intermittency and increasing contributions to measurements of $\varGamma$ from the quiescent flow.

Some Further Aspects of Stable Boundary-Layer Simulation in a Stratified-Flow Wind Tunnel

It is demonstrated that the vertical profile of gradient Richardson number, Ri, can be shaped by control of the working-section inlet temperature profile. In previous work (Hancock and Hayden in Boundary-Layer Meteorol 168:20–57, 2018; 175:93–112, 2020; 180:5–26, 2021) the inlet temperature profile had been specified but without control of the profile of Ri in the developed-flow region of the working section. Control of the inlet temperature profile is provided by 15 inlet heaters (spread uniformly across the height of the working section), allowing control of the temperature gradient over the bulk of the boundary layer, and the overall temperature level above that of the surface. The bulk Richardson number for the 11 cases covers the range 0.01–0.17 (there is no overlying inversion). In the upper approx 2/3 of the boundary layer the Reynolds stresses and turbulent heat flux are controlled by the gradient in mean temperature, while in the lower approx 1/3 they are controlled both by this gradient and by the level above the surface temperature. In three examples, Ri is approximately constant at 0.07, 0.10 and 0.13 across the bulk of the layer. The previous observation of horizontally homogenous behaviour in the temperature profiles in the top approx 2/3 of the boundary layer but not in the lower approx 1/3 is repeated here, except when, tentatively, Ri does not exceed 0.05 over the bulk of the boundary layer. Favourable validation comparisons are made against two sets of local scaling systems and field data over the full depth of the boundary layer, over the range 0.006 le Ri le 0.3, or, in terms of height and local Obukhov length, 0.005 le z/L le 1.

Study of Some Stability Parameters in the Atmosphere of Oil Al-Dura Refinery, Southeast Baghdad

Wind and temperature measurements at an oil refinery site located southeast of Baghdad city at two levels, 15 and 30 m, are presented. Three schemes are used to determine different stability classifications: Monin-Obukhov length, gradient, and bulk Richardson numbers. Meanwhile, vertical changes in air temperature and wind shear are also computed. There were lapse rate and inversion cases during the nights and days while favorable wind shear was dominant. The variation of stability in each scheme is large, covering the entire range of stability for a given class. The results of stability schemes are compared to each other. The results show that the schemes based on gradient and bulk Richardson numbers reasonably compare them.

Turbulent Inflow Generation for Large-Eddy Simulation of Winds around Complex Terrain

Accurate turbulent inflow conditions are needed to broaden the application of the large-eddy simulation technique to predict winds around arbitrarily complex terrain. We investigate the concept of buoyancy perturbations with colored noise to trigger turbulence in upstream flows approaching complex terrain regions. Random perturbations are imposed on the source term in the pseudo-temperature transport equation. These perturbations are effective within three-dimensional boxes and scaled using a bulk Richardson number defined for each box. We apply the turbulent inflow generation technique to predict winds around the Askervein and Bolund Hills under neutrally stratified conditions. We find that a common value for the bulk Richardson number works well for a variety of flow problems. Additionally, we show that the height of the perturbation box plays an important role in the accuracy of the predictions around complex terrain. We consistently obtained good results for both simulation cases when the perturbation box height was made a fraction of the Obukhov length scale.

Effect Dynamic Stability of Atmospheric Boundary Layer on Plume Downward Flux Emitted from Daura Refinery Stacks

The vertical motion of air that enhances or restricts the atmospheric turbulence through atmospheric boundary layer is known as "atmospheric stability", in which any movement of atmospheric components such as water vapor, aerosols, etc. are affected by the atmospheric stability. The aim of this research is to test the effect of atmospheric boundary layer stability on the amount of downward aersols flux at 10um (PM10) emission from stacks of Daura refinery, and estimated of deposition dust aerosols PM10 amount in area around the Daura refinery. In this study, hourly atmospheric stability based on similar theory of Monin-Obukhov length is calculated from archived data of the European Center for Medium Range Weather Forecast, and deposition velocity (Vd) for PM10 that is emitted from stacks Daura Refinery is calculated using the stability parameter (L). PM10 concentration is estimated according to the Gaussian model, which is used along with deposition velocity at this partical size in order to produce downward sedimentation flux (Fp) at distances 1000, 5000 and 10000 m from stack point sources emission. Results show that areas located to the south and southeast of the refinery receive large amounts of deposited flux values through stable weather conditions, where the accumulated PM10 amounts during one month have recorded 1.5 million μg /m2.s in January at a distance of 1000 m from refinery center stacks, while this amount reaches 532 million μg/m2.s during July due to the high emission rates resulting from burning fuel oil during July. The percentages of PM10 sedimentation decreased with the distance from the refinery to 1712 and 322839 μg /m2.s at a distance 10 km from the refinery in January and July, respectively. According to this method, the accumulated amount of PM10 in square meters can be estimated at any time, if atmospheric stability conditions and the domain of wind direction are known.