Next Article in Journal
Integrating Cost-Effective Measurements and CFD Modeling for Accurate Air Quality Assessment
Next Article in Special Issue
Comparative Study of CALPUFF and CFD Modeling of Toxic Gas Dispersion in Mountainous Environments
Previous Article in Journal
As Regular as Clockwork: Alexander von Humboldt, Robert de Lamanon and the Beginning of the Scientific Investigation of the Tidal Barometric Oscillation
Previous Article in Special Issue
Application of Effective Conversion Rates between NO and NO2 in a Standard Airport Dispersion Model System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized 3D Model of Crosswind Concentrations and Deposition in the Atmospheric Boundary Layer

Department of Mathematics, Faculty of Sciences, University Chouaib Doukkali, El Jadida 24100, Morocco
*
Author to whom correspondence should be addressed.
Atmosphere 2024, 15(9), 1054; https://doi.org/10.3390/atmos15091054
Submission received: 12 July 2024 / Revised: 18 August 2024 / Accepted: 28 August 2024 / Published: 31 August 2024

Abstract

:
In this paper, we introduce a comprehensive solution aimed at enhancing our understanding of three-dimensional atmospheric pollutant dispersion. This innovative solution involves the development of a generalized model that extends previous research and is applicable to all parameterization schemes of these equations, including wind speed profiles and turbulent diffusion coefficients, while incorporating the dry deposition criterion. Our methodology involves subdividing the atmospheric boundary layer into distinct sub-layers, which facilitates a detailed examination of pollutant dispersion dynamics. Extensive validation with data from the Hanford experiment has demonstrated the accuracy of this solution in simulating pollutant concentrations. The results demonstrate that there is a strong correlation between the projected and observed concentrations, underscoring the statistical reliability of our approach. This validation situates the statistical indices of our solution within an acceptable range, confirming its accuracy in predicting atmospheric pollutant dispersion. These findings thus establish our solution as a valid and effective method for studying complex environmental phenomena.

1. Introduction

Understanding the dispersion, transport, and fate of atmospheric pollutants necessitates a thorough analysis of deposition processes. These mechanisms are essential in determining plume behavior and concentration patterns within the atmospheric boundary layer ABL) [1,2]. However, Gaussian plume model estimates often fall short due to particle settling, deposition, or vertical gradients in diffusivity and wind speed [3]. The Earth’s surface acts as the primary reservoir for trace gases and particles, which transfer through wet and dry deposition pathways. Wet deposition is driven by precipitation, whereas dry deposition involves the adhesion or chemical reaction of trace gases and particles on surfaces like vegetation, soil, water bodies, and structures. Factors such as atmospheric turbulence intensity, the properties of gases and particles, and surface characteristics influence the efficiency of dry deposition [4]. In the absence of turbulence, even surfaces with high adsorptive capacity are ineffective in removing suspended materials, as the process is limited by molecular diffusion. Gravitational settling also plays a role in the dry deposition of particles, further complicating these dynamics [5].
The role of deposition processes in the dispersion of passive contaminants within the ABL is conventionally modeled using the advection–diffusion equation, as dictated by gradient transport theory, with deposition fluxes incorporated as boundary conditions at the surface [6]. To address this, a modified Gaussian model was formulated, initially presuming constant wind speed and homogeneous eddy diffusivity [7]. For enhanced applicability, eddy diffusivities were subsequently adjusted based on dispersion parameters relative to the downwind distance, and wind speed was modeled using a power-law relationship with vertical height [8]. However, this method introduces mathematical inconsistencies due to the employment of unrealistic parameterizations [9].
Given these challenges, there is a clear need for more flexible and accurate models that can better represent the complex dynamics of pollutant dispersion in the ABL. In this study, we address this need by introducing a generalized analytical solution that holistically integrates various wind speed profiles, turbulent diffusion coefficients, and dry deposition criteria. This model subdivides the ABL into distinct sub-layers, each parameterized with specific wind speed profiles and turbulent diffusivity coefficients. This segmented approach allows for a detailed analysis of vertical and horizontal variations in atmospheric properties, thereby providing a more precise understanding of the mechanisms of pollutants’ transport and deposition.
The key contributions of this work include the development of a generalized analytical model, the incorporation of closed-form solutions, and the meticulous validation of these models using empirical data. By addressing the intricacies of turbulent dispersion and dry deposition, our research offers a comprehensive framework that enhances the predictability and our understanding of atmospheric pollutant behavior. The results demonstrate that using 4 to 6 sub-layers offers an optimal balance between accuracy and computational efficiency, with strong correlations to observed data from the Hanford experiment.
To achieve the objectives of this study, this paper is structured as follows: Section 2 covers the modeling of the advection–diffusion equation. Section 3 explains the construction of the solution and the employed parameterizations. Section 4 showcases the results and provides a statistical performance analysis. Finally, Section 5 concludes the study.

2. Analytical Solutions to the Atmospheric Pollutant Equation

The advection–diffusion equation related to atmospheric pollution is fundamentally expressed as a principle of the conservation of suspended particles. This statement is mathematically formalized by the following expression [10]:
D c t o t D t = c t o t t + · ( U t o t c t o t ) = S
where
  • D D t represents the substantial derivative, capturing the rate of change in pollutant concentration from the perspective of a moving air parcel, embodying both the temporal and spatial variations due to airflow.
  • C t o t stands for the total concentration of pollutants, which includes both the average concentration and the transient fluctuations attributable to atmospheric turbulence.
  • U t o t denotes the total wind velocity vector represented by ( U , V , W ) T , where U indicates the east–west component, V the north–south component, and W the vertical component, encompassing both the average wind speed and its instantaneous fluctuations.
  • / t indicates the partial derivative with respect to time, denoting how the concentration of pollutants changes at a fixed point as time progresses.
  • is the gradient operator, employed to calculate the spatial gradient of the pollutant’s concentration.
  • · ( U t o t c t o t ) describes the divergence of the wind velocity and pollutant concentration product, representing the net movement of pollutants carried by the wind across the spatial domain.
  • S is the source term within the equation, quantifying the addition or removal of pollutants from the atmosphere via various processes such as industrial emissions, chemical reactions, or their deposition and uptake by terrestrial surfaces and vegetation.
Reynolds decomposition [11] is a pivotal technique used for analyzing turbulent flows. It separates the total wind velocity U t o t and pollutant concentration c t o t into their average components, U and c, and their fluctuating components, U and c , respectively. This separation is mathematically represented as
U t o t = U + U
and
c t o t = c + c
where the prime notation (’) indicates the fluctuating part of the quantity. By doing so, this facilitates an understanding of the complex interplay between the mean flow and the superimposed random fluctuations that characterize the turbulent diffusion of contaminants in the atmosphere.
By incorporating the decomposed expressions for c t o t and U t o t , the advection–diffusion Equation (1) is reformulated as
( c + c ) t + · ( U + U ) ( c + c ) = S
Applying the averaging operator across the equation, the mean of the fluctuations alone cancels out, as they are, by definition, zero:
c t ¯ = 0 , U ¯ = 0 , c ¯ = 0
Expanding the term · ( U c ) , we find that
· ( U c ) = · ( U c + U c + U c + U c )
After averaging, only the terms involving the averages and the term U c remain:
· ( U c ) ¯ = · ( U c ) + · U c ¯
By assembling the mean terms and subtracting the turbulent diffusion term · U c ¯ from the advection term · ( U c ) , we ultimately derive the advection–diffusion equation for turbulent flow:
c t + U · c = · U c ¯ + S
This equation now accurately represents the transport and spread of pollutants in a turbulent atmosphere, accounting for the mean motion of the air as well as the chaotic eddies that enhance diffusion.
In atmospheric modeling, the fluctuating wind component, denoted as U = ( u , v , w ) T , captures the erratic behavior of air movement. The gradient transport hypothesis, also known as K-theory, elucidates the mechanism by which turbulence drives materials against their concentration gradients. This is mathematically expressed as follows:
U c ¯ = K c
where the eddy diffusivity matrix K , defined as diag ( K x , K y , K z ) , incorporates the eddy diffusivity coefficients K x , K y , and K z for the x-, y-, and z-directions, respectively. This setup enriches our understanding of turbulent diffusion and its impact on the spatial distribution of pollutants.
Furthermore, the equation U c ¯ = u c ¯ , v c ¯ , w c ¯ T accurately models turbulent fluxes, underscoring the model’s effectiveness in simulating atmospheric pollutants’ transport and distribution, and highlights the pivotal role of turbulent diffusion as described by Fick’s laws:
u c ¯ = K x c x , v c ¯ = K y c y , w c ¯ = K z c z .
Integrating Equation (2) with the mass continuity principle yields advection–diffusion Equation (1), as rigorously documented by [10]:
t c + U · c = · ( K c ) + S
In the physics and engineering disciplines, the concept of point forces—actions highly localized in both space and time—is a recurrent theme.To mathematically model such phenomena, the Dirac delta distribution is utilized, enabling the precise identification of point sources [12,13]. For instance, a point source emitting pollutants at a continuous rate Q at a specific location ( x s , y s , z s ) is represented by the source term S ( x s , y s , z s ) , which is formulated as
S ( x s , y s , z s ) = Q · δ ( x x s ) δ ( y y s ) δ ( z z s ) ,
where δ denotes the Dirac delta function [14] and ⊗ symbolizes the tensor product. In practical scenarios, especially with industrial chimneys, a coordinate system is often used when the origin is at the base of the emission source. For a source at x s = 0 , y s = 0 , and z s = H s , the source term simplifies to S ( 0 , 0 , H s ) , emphasizing the source’s position at the origin of the horizontal plane and at height H s above the ground.
To accurately model the transport of emitted pollutants, the following boundary conditions are considered:
(i)
The apex of the mixing or inversion layer functions as an impervious boundary preventing the diffusion of contaminants. Consequently, a reflective flux boundary condition is postulated at the height H m i x of the inversion layer, articulated mathematically as
K z ( x , z ) c z = 0 , at z = H m i x .
(ii)
Addressing the scenario where a pollutant is partially absorbed by the ground surface, Calder [15] delineates an applicable boundary condition in flux terms. It is posited that the ground surface may facilitate the partial removal or absorption of the pollutant via deposition, characterized by a deposition velocity V g :
K z ( x , z ) c z = V g c , at z = z 0 .
where z 0 represents the surface roughness length.
(iii)
The lateral diffusive flux diminishes with increasing distance from the source, approaching zero as y extends towards the lateral boundary L y , formalized by
K y ( x , z ) c y 0 ; | y | L y ,
Within the context of this article, the key assumptions include the following:
  • The system reaches a state of equilibrium dispersion ( c / t = 0 );
  • Contributions from v ( c / y ) and w ( c / z ) are omitted, given that the x-axis aligns with the mean wind direction, rendering the w and v components of wind velocity minor;
  • Turbulent diffusion along the mean wind direction is disregarded in comparison to advective transport, implying that
    u c x x K x c x .
These assumptions lead to the steady-state advection–diffusion equation, which is as follows:
( PS ) : u ( z ) c x = y K y c y + z K z c z + Q δ ( x ) δ ( y ) δ ( z H s ) ; ( x , y , z ) ] 0 , L x [ × ] L y , L y [ × ] 0 , H m i x [ , Under   boundary   conditions : c ( x , y , z ) = 0 ; x = 0 , K y ( x , z ) c ( x , y , z ) y 0 ; | y | L y , K z ( x , z ) c ( x , y , z ) z = V g c ( x , y , z ) ; z = z 0 , K z ( x , z ) c ( x , y , z ) z = 0 ; z = H m i x .
H m i x represents the height of the ABL.
Theorem 1.
The problem ( PS ) is equivalent to the problem ( PS ) defined by
( P ) : u ( z ) c x = y K y c y + z K z c z ; ( x , y , z ) ] 0 , L x [ × ] L y , L y [ × ] 0 , H m i x [ , U n d e r   b o u n d a r y   c o n d i t i o n s : u ( z ) c ( x , y , z ) = Q · δ ( x ) δ ( y ) δ ( z H s ) ; x = 0 , K y ( x , z ) c ( x , y , z ) y 0 ; | y | L y , K z ( x , z ) c ( x , y , z ) z = V g c ( x , y , z ) ; z = z 0 , K z ( x , z ) c ( x , y , z ) z = 0 ; z = H m i x .
Proof of Theorem 1.
The steps for establishing equivalence are outlined as follows:
  • Consider c P ( x , y , z ) as the solution for x > 0 of the system ( P ) , and define c PS ( x , y , z ) as follows:
    c PS ( x , y , z ) = c P ( x , y , z ) for x > 0 , 0 for x 0 .
    We will establish that c PS satisfies the requirements of system ( PS ) .
  • To improve solution regularity, we apply mollifiers ρ n , as referenced in [16], over ] 0 , + [ , which are defined as follows:
    c PS , n ( x , y , z ) = c PS ( x , y , z ) ρ n ( x ) = c P , n ( x , y , z ) ; x > 0 , 0 ; otherwise , = c P , n ( x , y , z ) · H ( x )
    with c P , n ( x , y , z ) = c P ( x , y , z ) ρ n ( x ) and using the Heaviside function H ( x ) , which is 1 for x > 0 and 0 otherwise, where * denotes convolution.
  • We differentiate c PS , n with respect to x, leading to
    c PS , n x = c P , n x · H ( x ) + c P , n · H ( x ) = c P , n x · H ( x ) + c P , n ( 0 , y , z ) · δ ( x ) ,
  • Integrating the advection and diffusion terms, we find
    u ( z ) c PS , n x · ( K c PS , n ) = u ( z ) c P , n x · ( K c P , n ) + u ( z ) c P , n ( 0 , y , z ) · δ ( x ) = Q · δ ( x ) δ ( y ) δ ( z H s ) .
    This confirms that c P corresponds to the system ( P ) when n approaches infinity, thus completing the proof of the equivalence between the systems ( PS ) and ( P ) .
The problem ( P ) can be transformed into a series of specific advection–diffusion problems defined for each sub-layer of the atmospheric boundary layer, which is subdivided into H levels. For each level, the formulations of the turbulent diffusivity coefficients and wind profile adopt average values, as cited in [17]. Furthermore, two boundary conditions are applied at the ABL’s extremities, at z = z 0 and at z = H m i x , supplemented by continuity conditions for both the concentration and its flux at the interfaces between sub-layers. These continuity conditions, as shown in Figure 1, are expressed mathematically as follows:
  • The concentration at the end of a sub-layer must equal the concentration at the beginning of the subsequent sub-layer to ensure the continuity of concentration across sub-layer interfaces, as denoted by the equation
    c h ( z = z h ) = c h + 1 ( z = z h ) , for h = 1 , 2 , , H 1
  • The concentration flux, which is the product of diffusivity ( K z ) and the concentration gradient with respect to the vertical axis (z), must also be continuous across the interfaces between sub-layers. This continuity is captured by the equation
    K z , h c h z ( z = z h ) = K z , h + 1 c h + 1 z ( z = z h ) , for h = 1 , 2 , , H 1 .
We consider the atmospheric boundary layer’s intricate dynamics, where the eddy diffusivities adopt separable formulations that illuminate the interplay between horizontal and vertical mixing processes [18]. Specifically, these formulations are expressed as
K y , h ( x , z ) = ζ y ( x ) u h ,
K z , h ( x , z ) = ξ ( x ) φ z , h .
With this foundation, the equation ( P H ) outlines the advection–diffusion dynamics across each sub-layer of the atmospheric boundary layer, represented through a set of partial differential equations, which are expressed as follows:
( P H ) : u h c h x = ζ y ( x ) u h 2 c h y 2 + ξ ( x ) φ z h 2 c h z 2 ; ( x , y , z ) ] 0 , L x [ × ] L y , L y [ × ] z h 1 , z h [ , h { 1 , , H } Under the conditions : φ z 1 c 1 ( x , y , z 0 ) z = V g c 1 ( x , y , z 0 ) ; c h 1 ( x , y , z h 1 ) = c h ( x , y , z h 1 ) , h { 2 , , H } , φ z h 1 c h 1 ( x , y , z h 1 ) z = φ z h c h ( x , y , z h 1 ) z φ z H c H ( x , y , z H ) z = 0
Theorem 2.
( Solution to the P H problem)
The concentration within the H t h sub-layer, denoted by c H ( x , y , z ) , which is the solution to the problem P H , is given by the following expression:
c H ( x , y , z ) = Q 2 π 0 x ζ y ( s ) d s exp y 2 4 0 x ζ y ( s ) d s · [ n = 1 exp γ n 2 0 x ξ ( s ) d s P H , n ( z ) · { h = 2 H P h , n ( H s ) + 2 cos ( λ 1 , n ( H s z 0 ) ) + v g sin ( λ 1 , n ( H s z 0 ) ) u 1 1 + v g 2 ( z 1 z 0 ) } + 1 h = 1 H u h ( z h z h 1 ) ]
where
P h , n = 1 h = 1 H u h ( z h z h 1 ) 1 / 2 ; i f γ n = 0 α h , n cos ( λ h , n z ) + β h , n sin ( λ h , n z ) P h , n ; i f γ n 0
and
  • The coefficient α h , n is determined by the following relationship:
    α h , n = φ z h λ h , n φ z h 1 λ h 1 , n 2 φ z h λ h , n [ ( cos ( ( λ h , n + λ h 1 , n ) z h 1 ) + φ z h λ h , n + φ z h 1 λ h 1 , n φ z h λ h , n φ z h 1 λ h 1 , n cos ( ( λ h , n λ h 1 , n ) z h 1 ) ) α h 1 , n + ( sin ( ( λ h , n + λ h 1 , n ) z h 1 ) φ z h λ h , n + φ z h 1 λ h 1 , n φ z h λ h , n φ z h 1 λ h 1 , n sin ( ( λ h , n λ h 1 , n ) z h 1 ) ) β h 1 , n ] ;
  • The coefficient β h , n is determined by the following relationship:
    β h , n = φ z h λ h , n φ z h 1 λ h 1 , n 2 φ z h λ h , n [ ( sin ( ( λ h , n + λ h 1 , n ) z h 1 ) + φ z h λ h , n + φ z h 1 λ h 1 , n φ z h λ h , n φ z h 1 λ h 1 , n sin ( ( λ h , n λ h 1 , n ) z h 1 ) ) α h 1 , n ( cos ( ( λ h , n + λ h 1 , n ) z h 1 ) φ z h λ h , n + φ z h 1 λ h 1 , n φ z h λ h , n φ z h 1 λ h 1 , n cos ( ( λ h , n λ h 1 , n ) z h 1 ) ) β h 1 , n ] ;
  • The norm P h , n is defined as follows:
    P H , n 2 = h = 1 H z h 1 z h u h P h , n ( s ) 2 d s = h = 1 H u h 2 λ h , n [ sin ( λ h , n ( z h z h 1 ) ) ( ( α h , n 2 β h , n 2 ) cos ( λ h , n ( z h + z h 1 ) ) + 2 α h , n β h , n sin ( λ h , n ( z h + z h 1 ) ) ) + λ h , n ( α h , n 2 + β h , n 2 ) ( z h z h 1 ) ]
  • The parameter λ h , n is characterized by the equation
    λ h , n = γ n u h φ z h
    and as the associated eigenvalues the eigenvalues γ n , for n N , associated with this problem are real and discrete. Additionally, the corresponding eigenfunctions P h , n are mutually orthogonal.
Lemma 1.
(Transformed P H Problem Formulation)
The P H problem, when a Fourier transform is applied to the lateral coordinate y, is reformulated in terms of the transformed variable χ h ω , leading to the following system:
u h χ h ω x = ξ ( x ) φ z h 2 χ h ω z 2 ; ( x , z ) ] 0 , L x [ × ] z h 1 , z h [ , h { 1 , , H } U n d e r t h e c o n d i t i o n s : φ z 1 χ 1 ω ( x , z 0 ) z = V g χ 1 ω ( x , z 0 ) ; χ h 1 ω ( x , z h 1 ) = χ h ω ( x , z h 1 ) ; , h { 2 , , H } φ z h 1 χ h 1 ω ( x , z h 1 ) z = φ z h χ h ω ( x , z h 1 ) z ; φ z H χ H ω ( x , z H ) z = 0 .
where χ h ω ( x , z ) is defined as
χ h ω ( x , z ) = c ^ h ω ( x , z ) · exp ( 2 π ) 2 ω 2 0 x ζ y ( s ) d s
and c ^ h ω ( x , z ) denotes the Fourier transformation of c h with respect to y, which is given by
c ^ h ω ( x , z ) = + c h ( x , y , z ) e 2 i π ω y d y , h { 1 , , H }
Proof of Lemma 1.
In proving this lemma, we adhere to the following steps:
  • We begin by applying the Fourier transform with respect to the lateral coordinate y, as presented in Equation (16), to the P H problem. The resulting equation is
    u h c ^ h ω x + ( 2 π ) 2 ω 2 ζ y ( x ) c ^ h ω = ξ ( x ) φ z h 2 c ^ h ω z 2 , z z h 1 , z h .
  • We calculate the derivative of χ h ω ( x , z ) with respect to x, obtaining
    χ h ω x = c ^ h ω x + ( 2 π ) 2 ω 2 ζ y ( x ) c ^ h ω exp ( 2 π ) 2 ω 2 0 x ζ y ( s ) d s .
  • We multiply both sides of Equation (17) by the exponential term exp ( 2 π ) 2 ω 2 0 x ζ y ( s ) d s . Given that u h and φ z h remain constant within each sub-layer, this simplifies to
    u h χ h ω x = ξ ( x ) φ z h 2 χ h ω z 2 , z z h 1 , z h
  • We apply the same method to the boundary conditions, expressing them in terms of χ h ω to ensure consistency with the transformed system.
Having established the lemma, we now proceed to employ its results to demonstrate the theorem concerning the solution to the P H problem.
Proof of Theorem 2.
The construction of the solution for the P H problem unfolds through the following steps:
  • Applying Lemma 1, we transform the P H problem into a formulation that depends on χ h ω .
  • Representing χ h ω ( x , z ) as a series, we express it in a separated-variables form
    χ h ω ( x , z ) = n = 0 G h , n ω ( x ) P h , n ( z ) , h { 1 , , H } .
  • In verifying that the functions G h , n ω and P h , n satisfy their respective differential equations, we must therefore confirm that G h , n ω satisfies the equation
    d G h , n ω d x + γ n 2 ξ ( x ) G h , n ω = 0 , x ] 0 , L x [
    and P h , n ( z ) adheres to the system
    φ z h d 2 P h , n d z 2 + γ n 2 u h P h , n = 0 ; z ] z h 1 , z h [ , h { 1 , , H } U n d e r t h e c o n d i t i o n s : φ z 1 d P 1 , n ( z 0 ) d z = V g P 1 , n ( z 0 ) ; P h 1 , n ( z h 1 ) = P h , n ( z h 1 ) ; , h { 2 , , H } φ z h 1 d P h 1 , n ( z h 1 ) d z = φ z h d P h , n ( z h 1 ) d z ; φ z H d P H , n ( z H ) d z = 0 .
  • Determining the specific forms of G h , n and P h , n ensures that they are, respectively, expressed as
    G h , n ω ( x ) = μ n ( ω ) · exp γ n 2 0 x ξ ( s ) d s
    where μ n is an arbitrary function depending on ω , and
    P h , n ( z ) = α h , n cos λ h , n z + β h , n sin λ h , n z
    where λ h , n is expressed in Equation (14). This formulation must meet the system’s criteria, resulting in the determination of α h , n and β h , n , as specified in Equations (11) and (12), respectively.
  • The continuity condition between sub-layers necessitates a recursion of the eigenvalues λ h , n and the eigenfunctions P h , n for h { 1 , , H } . This process begins by determining the eigenvalues associated with the first layer, λ 1 , n , using the conditions related to P 1 , n . The resulting expressions are enumerated below:
    The function P 1 , n ( z ) is expressed as
    P 1 , n ( z ) = cos ( λ 1 , n ( z z 0 ) ) + v g sin ( λ 1 , n ( z z 0 ) )
    where v g = V g φ z 1 λ 1 , n .
    The norm of P 1 , n is given by
    P 1 , n 2 = u 1 ( z 1 z 0 ) ; when γ n = 0 u 1 1 + v g 2 z 1 z 0 2 ; when γ n 0
    The recursive generation of eigenvalues for the intermediate and final layers depends on the known values of the previous layers, in accordance with their respective conditions. This approach makes it possible to construct the associated eigenfunctions, succinctly maintaining continuity between the layers. For further details, please refer to Appendix A.1.
  • The Sturm–Liouville theorem [19] allows for the expanded representation of χ H ω as
    χ H ω ( x , z ) = n = 0 a n exp γ n 2 0 x ξ ( s ) d s P H , n ( z ) .
    where the eigenfunctions P h , n form a complete and mutually orthogonal set. Consequently, χ H ω is detailed as above. The coefficients a n are determined by leveraging the orthogonality of P h , n and calculated as follows:
    a n = = 1 H z 1 z u χ ω ( 0 , z ) P , n ( z ) d z = Q h = 1 H u h ( z h z h 1 ) , when γ n = 0 Q = 1 H P , n ( H s ) , when γ n 0
  • Applying the inverse Fourier transform to χ H ω ( x , z ) from Equation (16), and considering that the inverse transform of
    exp ( 2 π ) 2 ω 2 0 x ζ y ( s ) d s
    is given by
    1 2 π 0 x ζ y ( s ) d s exp y 2 4 0 x ζ y ( s ) d s ,
    leads to the final expression of the solution c H to the problem P H .
The crosswind-integrated concentration (CIC) model is a tool for analyzing the dispersion of pollutants in the atmosphere, accounting for the horizontal wind component crucial for pollutant transport [20,21]. It integrates the concentration across the lateral wind direction, y, yielding a two-dimensional concentration map along the wind direction (x) and crosswind direction (z) [4]. Mathematically, the CIC, c H y ( x , z ) , is defined by integrating the pollutant concentration c H ( x , y , z ) over y from to + [22]. This integration reflects the dispersal of pollutants due to the horizontal wind and assesses their impact in both the windward and lateral directions [23].
Corollary 1.
The crosswind-integrated concentration (CIC), c H y ( x , z ) , is equal to χ H ω ( x , z ) as follows:
c H y ( x , z ) = Q · [ n = 1 exp γ n 2 0 x ξ ( s ) d s P H , n ( z ) · { 2 cos ( λ 1 , n ( H s z 0 ) ) + v g sin ( λ 1 , n ( H s z 0 ) ) u 1 1 + v g 2 ( z 1 z 0 ) + h = 2 H P h , n ( H s ) } + 1 h = 1 H u h ( z h z h 1 ) ]
Proof. 
To confirm the corollary, it is sufficient to verify the Gaussian integral identity
+ 1 2 π 0 x ζ y ( s ) d s exp y 2 4 0 x ζ y ( s ) d s d y = 1 .
Corollary 2.
The solution c H for the P H problem encompasses and generalizes the previously established solutions under specific parameterizations:
1. 
The Single-Layer and Two-Dimensional Solution:
(a) 
For a single layer without discretization, where ζ y 1 and ξ 1 , our method yields the solution documented by Kumar [24].
(b) 
Transitioning to v g = 0 , ζ y 1 , ξ 1 , and z 0 = 0 aligns with the solution documented by Seinfeld [22].
2. 
The Multi-Layer and Three-Dimensional Solution:
(a) 
Setting v g = 0 retrieves the solution identified by Farhane et al. [18].
(b) 
Further simplification to v g = 0 and ξ 1 leads to the solution presented by Marie et al. [25].
Proof. 
For a detailed proof of this corollary, please refer to Appendix B. □
Corollary 3.
In the scenario of a single-layer framework, where z 0 approaches 0 and z 1 extends to H m i x , the crosswind-integrated concentration (CIC), c H y ( x , z ) , is defined as follows:
c 1 y ( x , z ) = Q u ¯ H m i x [ 1 + 2 1 + v g 2 n = 1 exp ( n π ) 2 φ z ¯ u ¯ H m i x 2 x { 1 + v g 2 cos n π H m i x ( z H s ) + 1 v g 2 cos n π H m i x ( z + H s ) + 2 sin n π H m i x ( z + H s ) } ]
where
  • u ¯ represents the average wind speed over the height of the mixing layer, calculated as
    u ¯ = 1 H m i x 0 H m i x u ( s ) d s ,
  • φ z ¯ denotes the average eddy diffusivity over the same height, expressed as
    φ z ¯ = 1 H m i x 0 H m i x φ z ( s ) d s .
Proof. 
To validate the crosswind-integrated concentration (CIC) for a scenario where H = 1 (a single-layer framework), we follow these steps:
  • Substitute Single-Layer Values into the CIC Equation—begin by inserting the conditions of a single-layer setup into the original CIC formula from Equation (27):
    c 1 y ( x , z ) = Q 1 u 1 ( z 1 z 0 ) + 2 n = 1 exp γ n 2 0 x ξ ( s ) d s cos λ 1 , n ( z z 0 ) + v g sin λ 1 , n ( z z 0 ) u 1 ( 1 + v g 2 ) ( z 1 z 0 ) + cos λ 1 , n ( H s z 0 ) + v g sin λ 1 , n ( H s z 0 ) u 1 ( 1 + v g 2 ) ( z 1 z 0 )
    where γ n 2 = ( n π ) 2 φ z ¯ u ¯ ( z 1 z 0 ) 2 .
  • Apply Limits as z 0 Approaches 0 and z 1 Extends to H m i x —evaluate the limits to reflect the shift to a single-layer framework that spans the entire mixing height H m i x :
    c 1 y ( x , z ) = Q u ¯ H m i x 1 + 2 1 + v g 2 n = 1 exp ( n π ) 2 φ ¯ z u ¯ H m i x 2 0 x ξ ( s ) d s cos n π z H m i x + v g sin n π z H m i x · cos n π z H m i x + v g sin n π z H m i x = Q u ¯ H m i x [ 1 + 2 1 + v g 2 n = 1 exp ( n π ) 2 φ z ¯ u ¯ H m i x 2 0 x ξ ( s ) d s { 1 + v g 2 cos n π H m i x ( z H s ) + 1 v g 2 cos n π H m i x ( z + H s ) + 2 sin n π H m i x ( z + H s ) } ]

3. Turbulent Parameters in Pollutant Dispersion Modeling

(a)
Parameterizing the Wind Speed Profile,  u ( z ) .
The wind speed profile u ( z ) is principally derived using surface layer theory enhanced by Monin–Obukhov (M-O) similarity theory [26]. While traditionally limited to the surface layer, extending the M-O scaling across the entire atmospheric boundary layer (ABL) involves maintaining a constant wind speed profile above this initial layer, ensuring the model’s extended applicability and continuity.
Gryning et al. [27] have refined this approach to create a comprehensive wind speed profile applicable over the entire ABL, which is particularly suited for homogeneous terrains. Their advanced model employs a local mixing length composed of three component scales, leading to a robust formulation that incorporates adjustments for atmospheric stability throughout the entire ABL:
u ( z ) = u k ln z z 0 + b z L M O 1 z 2 H m i x + z L M B L z H m i x z 2 L M B L ,
where b = 4.7 is an empirically determined constant, and L M O , the Monin–Obukhov length [26], plays a critical role in modeling stability within the boundary layer. The mid-boundary layer length, L M B L , crucial for stable conditions, is given by
L M B L = H m i x 2 ln u f c z 0 B ( μ ) 2 + A 2 ( μ ) 1 / 2 ln H m i x z 0 b H m i x 2 L M O 1 ,
In this formulation, f c represents the Coriolis parameter, k = 0.4 is the von Karman constant, and A 4.9 and B 1.9 are resistance law functions dependent on the stability parameter μ . Additionally, Gryning et al. [27] proposed an empirical formula for L M B L that enhances the model’s ability to adapt to various atmospheric stability conditions:
L M B L = u f c 2 ln u f c z 0 + 55 exp u / f c z 0 2 400 .
(b)
Parameterizing the Diffusivity Profiles
(b.1)
Vertical Diffusivity:  K z ( x , z )
In atmospheric dispersion models, the parameterization of eddy diffusivity, K z , significantly influences the model’s accuracy and performance. Traditionally, K z has been modeled as solely a function of the height above the ground, z. However, recent advancements have extended this parameterization to include the downwind distance x, recognizing its impact on near-source dispersion dynamics.
An innovative approach by Mooney and Wilson [28] introduced a composite model for K z , blending both spatial variables into a single expression:
K z ( x , z ) = φ z ( z ) 1 exp x L 1 ,
where φ z represents the vertical eddy diffusivity as a function of height and L 1 = u ( H s ) τ ( H s ) denotes the along-wind length scale. This scale delineates the transition across which diffusivity reaches its extensive field value. Here, τ signifies the Lagrangian time scale at the source height, reflecting the turbulence’s persistence in the atmosphere.
The parameterization of τ at the source height H s is given by
τ ( z ) = φ z ( z ) σ w 2 ,
with σ w representing the vertical turbulent intensity. In stable atmospheric conditions, the formulation for σ w is
σ w = 1.4 u 1 z H m i x 0.75 ,
which demonstrates how stability affects turbulence characteristics.
Mangia et al. [29] further refined φ z based on local similarity and statistical diffusion theories, resulting in the following relation:
φ z ( z ) = 0.3 u z ( 1 z / H m i x ) 1 + 3.7 z / Λ
where Λ represents the local Obukhov length and is expressed as
Λ = L M O 1 z H m i x ( 3 α 1 / 2 ) α 2 ,
utilizing the constants α 1 = 3 / 2 and α 2 = 1 , which are derived from the second-order closure model by Nieuwstadt [30], which accounts for stationary turbulence under a constant cooling rate.
(b.2)
Horizontal Diffusivity:  K y ( x , z )
In the realm of atmospheric dispersion modeling, lateral eddy diffusivity plays a crucial role in determining the spread of pollutants in the crosswind direction. This diffusivity is fundamentally expressed through a relationship that ties it to the wind speed profile and the variance in lateral dispersion:
K y ( x , z ) = 1 2 u ( z ) d σ y 2 ( x ) d x ,
where σ y represents the standard deviation of the pollutant’s distribution in the lateral or crosswind direction and is predominantly a function of the downwind distance x. This formulation underscores the dependency of lateral diffusivity on both the vertical wind speed profile and the rate of change in lateral dispersion as pollutants travel away from their source.
This parameterization, based on the foundational work cited in Huang’s theory [31], encapsulates the dynamic interactions between wind speed and atmospheric turbulence. By integrating the vertical wind profile u ( z ) with the derivative of the squared standard deviation σ y 2 with respect to x, the model effectively captures how lateral diffusion varies with both altitude and distance from the source.

4. Results: Data Analysis and Model Validation

The Hanford experiments provided a valuable dataset for the evaluation of our model of crosswind dispersion and deposition dynamics within the atmospheric boundary layer. These experiments, carried out in May and June of 1983, took place in a semi-arid region of southeastern Washington characterized by a largely flat terrain. A detailed documentation of these experiments is available in the work by Doran and Horst [32]. Measurements were taken from six dual-tracer releases positioned at distances of 1100, 200, 800, 1600, and 3200 m from the release point, under conditions ranging from moderately stable to near-neutral. Notably, the deposition velocity was only evaluated for the three farthest distances. There were six test runs in total, as summarized in the compiled Table 1. Each run had a release duration of 30 minutes, except for the fifth run, which was 22 minutes. The data collected, detailed in Doran et al. [33], were formatted as crosswind-integrated tracer concentrations.
Samplers were positioned at angular intervals of 8°, 4°, 4°, 2°, and 3° on concentric circles centered around the release point, with radii corresponding to 100, 200, 800, 1600, and 3200 m, respectively. Crosswind-integrated concentrations offer significant advantages over individual azimuthal values, as they provide a smoother dataset by integrating over multiple samplers, thereby reducing the impact of minor variations in sampler performance. This method also mitigates issues arising from non-Gaussian crosswind tracer distributions. The terrain had a roughness length of 3 cm. The experiment involved the simultaneous release of two tracers: zinc sulfide (ZnS) as the depositing tracer and sulfur hexafluoride (SF6) as the non-depositing tracer, both released from a height of 2 m. The release points for SF6 and ZnS were less than 1 meter apart. Table 1 provides the micrometeorological data, including the Monin–Obukhov length ( L M O ), friction velocity ( u ), and the height of the mixing layer ( H mix ). It also lists the crosswind-integrated tracer concentrations, adjusted for the release rate (Q) and the computed effective deposition rates ( V g ) using the surface depletion method described by Doran and Horst [32].
The ZnS tracer, being a polydisperse aerosol, presented challenges in achieving consistent-size distribution estimates. This inconsistency hinders the accurate prediction of ZnS deposition rates, as these depend on its mass distribution. The deposition velocity ( V g ) is critical for model evaluation, but estimates of V g in the literature vary with particle size, atmospheric turbulence, and surface characteristics. Combined with uncertainties in the tracer’s particle size distribution, this variability complicates the accurate prediction of deposition velocity. Since V g changes with particle size, both size distribution and deposition velocity vary with downwind distance. An effective deposition velocity ( V g ) for a given distance incorporates these local variations in V g from the source to that distance. For comprehensive details on calculating effective deposition rates and wind velocity, see Doran and Horst [32].

Results and Discussions

The model was assessed by comparing the crosswind-integrated concentrations of the pollutants ZnS and SF 6 at a height of 1.5 m above the ground, normalized by the emission rate Q. The comparison of the predicted and observed values is illustrated in Figure 2. The vertical eddy diffusivity was parameterized according to the scheme provided in Equation (30), while the wind profile followed the formulation described in Equation (29).
Figure 2 presents scatter plots of the predicted versus observed concentrations of the non-depositing material ( SF 6 ) (Figure 2a,b) and the depositing material ( ZnS ) (Figure 2c). Data points located within the outer two dotted lines indicate agreement within a factor of 2, while the central diagonal line represents perfect agreement. It is evident that all simulated concentration data for both the non-depositing material ( SF 6 ) and the depositing material ( ZnS ) fall within the dotted lines, indicating a good fit and suggesting that our model is more accurate when deposition processes are included. This contrast between ZnS and SF 6 underscores the significant impact of deposition on model accuracy.
Figure 2a,b analyze the non-depositing material ( SF 6 ) using a logarithmic scale and a linear scale, respectively. The logarithmic scale emphasizes lower concentration ranges, providing a more detailed insight into the effects of deposition, while the linear scale offers a straightforward comparison across the entire concentration range. Similarly, Figure 2c,d present the depositing material ( ZnS ) using both logarithmic and linear scales, demonstrating the model’s consistency and reliability in different deposition scenarios.
The results further demonstrate that increasing the discretization of the model improves the accuracy of its predictions. Higher levels of discretization provide a better average representation of the vertical profiles of pollutants, thereby enhancing the model’s capability to predict pollutant dispersion accurately. This improvement is particularly noticeable when comparing predicted and observed values, as the model with higher discretization shows a closer fit to the empirical data.
Overall, this analysis reveals that the model performs well in simulating pollutant dispersion, particularly for scenarios involving deposition. The findings emphasize the importance of using appropriate parameterizations and the benefits of higher discretization in improving model accuracy. This comprehensive evaluation underscores the robustness of the model and its potential for generating reliable predictions in atmospheric pollutant dispersion studies.
Table 2 presents the performance metrics of the models with and without deposition, evaluated using well-established statistical indices. These include the NMSE (Normalized Mean Square Error), MRSE (Mean Relative Square Error), COR (Correlation Coefficient), FB (Fractional Bias), FS (Fractional Standard Deviation), MG (Geometric Mean Bias), VG (Geometric Variance), and FAC2 (Factor of Two). These indices provide insights into model accuracy and reliability [34].
For models without deposition ( V g = 0 ), increasing the number of discretization levels H improves their performance. At H = 6 , the model shows a very low NMSE and MRSE, high COR, and minimal FB and FS, indicating strong correlation, minimal bias, and variability. When H = 4 , its performance slightly decreases but remains high, with a low NMSE and high COR. At H = 2 , its NMSE and MRSE increase, and its COR decreases slightly, but its performance is still acceptable. For models with deposition ( V g 0 ), similar trends are observed. Higher H values result in a better performance, with H = 6 showing a high COR and low bias. This analysis indicates that all parameterizations simulate observed concentrations well, with NMSE, FB, and FS values close to zero and COR and FAC2 values close to one. Models without deposition generally perform better, but those with deposition also show good performance, particularly with higher discretization levels, demonstrating the model’s robustness and accuracy in predicting pollutant dispersion under various conditions.
Table 2 also presents comparisons between the proposed model and previous models established by Marie et al. [25] and Laaouaoucha et al. [35] for V g = 0 , as well as those of Moreira et al. [36] and Kumar and Sharan [24] for V g 0 . The results indicate the excellent performance of our proposed solution, as evidenced by its favorable statistical indices across various conditions and discretization levels. This demonstrates the effectiveness of the proposed model in accurately predicting pollutant dispersion, confirming its superiority and reliability.
Figure 3 illustrates the variations in downwind ground-level concentrations for various dry deposition velocities from a source height of 25 m, as well as the influence of the number of discretizations (H) on these concentrations. Figure 3a, Figure 3b, and Figure 3c correspond to H = 2 , H = 4 , and H = 6 , respectively.
As H increases, the concentration profiles, both with and without deposition, show improved alignment with the observed data. This improvement is attributable to a more precise representation of wind speed and turbulent diffusivity profiles thanks to higher levels of discretization. A greater number of layers allows for the better capturing of the average characteristics of these profiles.
These profiles show that the highest concentrations occur for contaminants that are completely reflected by the ground. Introducing dry deposition reduces potential human exposure by decreasing airborne concentrations [37]. As dry deposition velocities increase, the location of the peak concentration moves closer to the source. This shift is more evident in the sub-figures as the depletion of pollutants becomes more significant. This behavior is consistent with the predictions obtained using the Gaussian plume model [37].
This analysis demonstrates that increasing the number of discretization layers (H) enhances the model’s performance in predicting pollutant dispersion. Higher levels of discretization result in a more accurate depiction of vertical profiles, thereby reinforcing the credibility of the model’s predictions. This finding underscores the importance of detailed vertical discretization in atmospheric modeling to achieve accurate and reliable results.
Figure 4 illustrates the numerical convergence of the proposed solution with the measured concentration as the number of eigenvalues increases across different values of V g . The sub-figures correspond to (Figure 4a) V g = 0 , (Figure 4b) V g = 0.01 , (Figure 4c) V g = 0.02 , and (Figure 4d) V g = 0.03 , each evaluated at three downwind distances ( x = 800 m, x = 1600 m, and x = 3200 m). This analysis demonstrates that the solution converges rapidly with an increasing number of eigenvalues.
The results confirm that this convergence behavior is consistent across different deposition velocities. The concentration values stabilize after a certain number of eigenvalues, indicating the robustness and efficiency of the proposed method. This quick convergence ensures that the model’s predictions are reliable and can be computed efficiently, even for varying environmental conditions and distances from the source.

5. Conclusions

In this study, we developed a comprehensive model of the dispersion of atmospheric pollutants within the atmospheric boundary layer, incorporating various parameterizations for wind speed profiles, turbulent diffusivity, and dry deposition processes. Our approach subdivides the atmospheric boundary layer into distinct sub-layers, enabling a detailed examination of pollutant dispersion dynamics. The analytical solutions obtained were validated against empirical data, demonstrating the model’s strong correlation and reliability. The results underscore the model’s robustness and accuracy, particularly with higher levels of discretization. This work not only extends existing models but also provides a framework for more precise predictions of pollutant behavior under different atmospheric conditions. Future research could focus on further refining these parameterizations and exploring their applicability in diverse environmental scenarios.

Author Contributions

M.F. was the lead author responsible for the writing, editing, and corrections. O.S. supervised the project. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ABLAtmospheric boundary layer
CICCrosswind-integrated concentration

Appendix A. Construction of the Solution Associated with the First Sub-Layer and an Eigenvalue Analysis of the Problem (uid56)

Appendix A.1. Construction of the Solution Associated with the First Sub-Layer of the Problem (uid56)

The solution to Problem (21) is articulated in Equation (23), where α h , n and β h , n are determined based on the specified conditions for h { 2 , , H } . The continuity conditions establish a recursive relationship between the eigenvalues γ n = λ h , n φ z h u h and the associated eigenfunctions P h , n for h { 1 , , H } . Initially, the eigenvalues for the first sub-layer (when h = 1 ) must be calculated using the associated condition for P 1 , n , which is
φ z 1 d P 1 , n ( z 0 ) d z = V g P 1 , n ( z 0 ) .
Considering the expression for P 1 , n is
P 1 , n ( z ) = α 1 , n cos λ 1 , n z + β 1 , n sin λ 1 , n z ,
substituting P 1 , n ’s expression into the boundary condition yields
α 1 , n sin ( λ 1 , n z 0 ) + v g cos ( λ 1 , n z 0 ) β 1 , n cos ( λ 1 , n z 0 ) v g sin ( λ 1 , n z 0 ) = 0 ,
where v g = V g φ z 1 λ 1 , n .
To satisfy this equation, we can select
α 1 , n = cos ( λ 1 , n z 0 ) v g sin ( λ 1 , n z 0 ) β 1 , n = sin ( λ 1 , n z 0 ) + v g cos ( λ 1 , n z 0 )
which, when substituted into the expression for P 1 , n , yield
P 1 , n ( z ) = cos ( λ 1 , n ( z z 0 ) ) + v g sin ( λ 1 , n ( z z 0 ) ) .

Appendix A.2. Eigenvalue Analysis of the Problem (uid56)

In the analysis of the first sub-layer, the eigenvalues, represented by γ n > 0 , are inferred from the norm of P 1 , n . Indeed,
P 1 , n 2 = z 0 z 1 u 1 P 1 , n ( s ) 2 d s = u 1 4 λ 1 , n [ 2 λ 1 , n 1 + v g 2 z 1 z 0 + 1 v g 2 sin 2 λ 1 , n z 1 z 0 v g cos 2 λ 1 , n z 1 z 0 1 ]
It is paramount that the norm does not depend on any specific function form. The norm, representing a measure of magnitude or distance, should remain constant across different functions to ensure uniformity and comparability. Therefore, the eigenvalues γ n are determined under conditions that guarantee this property, leading to
sin 2 λ 1 , n ( z z 0 ) = 0 , cos 2 λ 1 , n ( z z 0 ) 1 = 0 ,
which ultimately define the eigenvalues as
λ 1 , n = n π z 1 z 0 , n N ,
yielding a norm expression that is as follows:
P 1 , n 2 = u 1 ( 1 + v g 2 ) z 1 z 0 2 , n N .
For n = 0 , the norm is specifically evaluated by considering P 1 , 0 = 1 in Equation (A5).
Following the complete assessment of eigenvalues for the first sub-layer, our focus shifts to the intermediate and final layers. For the intermediate layers, by implementing the specific condition
P h 1 , n ( z h 1 ) = P h , n ( z h 1 ) ; , h { 2 , , H } φ z h 1 d P h 1 , n ( z h 1 ) d z = φ z h d P h , n ( z h 1 ) d z ,
we arrive at the equation governing the eigenvalues for these layers, which is presented as
φ z h λ h , n [ α h , n sin ( λ h , n z h 1 ) β h , n cos ( λ h , n z h 1 ) ] φ z h 1 λ h 1 , n [ α h 1 , n sin ( λ h 1 , n z h 1 ) β h 1 , n cos ( λ h 1 , n z h 1 ) ] = 0 .
Regarding the last layer, employing the relevant condition
φ z H d P H , n ( z H ) d z = 0
leads to the formulation of the equation that yields the eigenvalues for this ultimate layer:
α H , n sin ( λ H , n H m i x ) β H , n cos ( λ H , n H m i x ) = 0 .
To facilitate a clear understanding of the eigenvalue generation process for each sub-layer, the accompanying flowchart in Figure A1 concisely outlines the analytical steps involved, capturing the recursive determination of eigenvalues and the construction of the eigenfunctions as discussed.
Figure A1. Flowchart of the process for calculating eigenvalues and eigenfunctions in the subdivided media of the atmospheric boundary layer.
Figure A1. Flowchart of the process for calculating eigenvalues and eigenfunctions in the subdivided media of the atmospheric boundary layer.
Atmosphere 15 01054 g0a1

Appendix B. Derivation of Specific Solutions for cH

In the special case where v g = 0 , Equation (9) simplifies to yield the solution as determined by Farhane et al. [18], which is defined as follows:
c ( x , y , z ) = 1 2 π 0 x ζ y ( s ) d s exp y 2 4 0 x ζ y ( s ) d s n = 0 a n exp γ n 2 0 x ξ ( s ) d s P H , n ( z ) .
where the coefficients a n are as defined in Equation (26). Further analysis, with ζ y 1 , leads to a solution analogous to that provided by Marie et al. [25].
By transitioning to two dimensions and employing the CIC model for a single layer, we utilize the boundary condition at the top of the atmospheric boundary layer, specifically φ z H d P H ( x , z H ) d z = 0 . This condition derives the parameters α H , n = cos ( λ 1 , n H m i x ) and β H , n = sin ( λ 1 , n H m i x ) , leading to the expression P H , n = cos ( λ H , n ( z H m i x ) ) . The norm of P H , n is then given as
P H , n 2 = ( H m i x z 0 ) 2 1 + sin 2 λ H , n ( H m i x z 0 ) 2 λ H , n ( H m i x z 0 ) ,
matching the solution identified by Kumar [24].
For Seinfeld’s solution [22], we have
c 1 y ( x , z ) = Q u 1 ( z 1 z 0 ) 1 + 2 n = 1 exp ( n π ) 2 φ z 1 u 1 ( z 1 z 0 ) 2 x cos ( λ 1 , n ( z z 0 ) ) cos ( λ 1 , n ( H s z 0 ) ) = z 0 0 Q u 1 z 1 1 + 2 n = 1 exp ( n π ) 2 φ z 1 u 1 z 1 2 x cos ( λ 1 , n z ) cos ( λ 1 , n H s ) = z 1 H m i x Q u 1 H m i x 1 + 2 n = 1 exp ( n π ) 2 φ z 1 u 1 H m i x 2 x cos n π H m i x z cos n π H m i x H s

References

  1. Tirabassi, T. Operational advanced air pollution modeling. PAGEOPH 2003, 160, 5–16. [Google Scholar] [CrossRef]
  2. Viúdez-Moreiras, D. Editorial for the Special Issue “Atmospheric Dispersion and Chemistry Models: Advances and Applications”. Atmosphere 2023, 14, 1275. [Google Scholar] [CrossRef]
  3. Smith, R.B. A K-theory of dispersion, settling and deposition in the atmospheric surface layer. Boundary-Layer Meteorol. 2008, 129, 371–393. [Google Scholar] [CrossRef]
  4. Pasquill, F.; Smith, F.B. Atmospheric Diffusion, 3rd ed.; Ellis Horwood Ltd.: Hertfordshire, UK, 1983. [Google Scholar]
  5. Hewitt, C.N.; Jackson, A. Handbook of Atmospheric Science; Blackwell Publishing: Hoboken, NJ, USA, 2003; p. 633. [Google Scholar]
  6. Hanna, S.R. Applications in Air Pollution Modelling. In Atmospheric Turbulence and Air Pollution Modelling; Nieuwstadt, F.T.M., van Dop, H., Eds.; Reidel: Boston, MA, USA, 1982. [Google Scholar]
  7. Irwin, J.S. A theoretical validation of the wind profile power law exponent as a function of surface roughness and stability. Atmos. Environ. 1979, 13, 191–194. [Google Scholar] [CrossRef]
  8. Llewelyn, R.P. An analytical model for the transport, dispersion and elimination of air pollutants emitted from a point source. Atmos. Environ. 1983, 17, 249–256. [Google Scholar] [CrossRef]
  9. Sharan, M.; Singh, M.P.; Yadav, A.K. A mathematical model for the atmospheric dispersion in low winds with eddy diffusivities as linear functions of downwind distance. Atmos. Environ. 1996, 30, 1137–1145. [Google Scholar] [CrossRef]
  10. Zannetti, P. Air Pollution Modeling: Theories, Computational Methods and Available Software; Springer Science & Business Media: Berlin, Germany, 2013; Volume 6. [Google Scholar]
  11. Visscher, A. Air Dispersion Modeling: Foundations and Applications; Wiley–Blackwell: Hoboken, NJ, USA, 2013; Volume 11. [Google Scholar]
  12. Dirac, P.A.M. The Principles of Quantum Mechanics, 4th ed.; Oxford University Press: Oxford, UK, 1958. [Google Scholar]
  13. Lighthill, M.J. Introduction to Fourier Analysis and Generalized Functions; Cambridge University Press: Cambridge, UK, 1958. [Google Scholar]
  14. Kanwal, R.P. Generalized Functions: Theory and Applications, 3rd ed.; Birkhäuser: Basel, Switzerland, 2004. [Google Scholar]
  15. Calder, K. Atmospheric diffusion of particulate material, considered as a boundary value problem. J. Meteorol. 1961, 18, 413–415. [Google Scholar] [CrossRef]
  16. Hörmander, L. The Analysis of Linear Partial Differential Operators; Springer: Berlin, Germany, 1990. [Google Scholar]
  17. Farhane, M.; Souhar, O. A Novel Development in Three-Dimensional Analytical Solutions for Air Pollution Dispersion Modeling. Environ. Sci. 2023, 9. [Google Scholar] [CrossRef]
  18. Farhane, M.; Alehyane, O.; Souhar, O. Three-dimensional analytical solution of the advection-diffusion equation for air pollution dispersion. ANZIAM J. 2022, 64, 40–53. [Google Scholar]
  19. Boyce, W.; DiPrima, R.; Meade, D. Elementary Differential Equations and Boundary Value Problems; Wiley: Hoboken, NJ, USA, 2017. [Google Scholar]
  20. Arya, S.P. Air Pollution Meteorology and Dispersion; Oxford University Press: Oxford, UK, 1999. [Google Scholar]
  21. Csanady, G.T. Turbulent Diffusion in the Environment; D. Reidel Publishing Company: Dordrecht, The Netherlands, 1973. [Google Scholar]
  22. Seinfeld, J.H.; Pandis, S.N. Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  23. Stockie, J.M. The Mathematics of Atmospheric Dispersion Modeling. SIAM Rev. 2011, 53, 349–372. [Google Scholar] [CrossRef]
  24. Kumar, P.; Sharan, M. A Generalized Analytical Model for Crosswind-Integrated Concentrations with Ground-Level Deposition in the Atmospheric Boundary Layer. Environ. Model. Assess. 2014, 19, 487–501. [Google Scholar] [CrossRef]
  25. Marie, E.; Hubert, B.; Patrice, E.; Zarma, A.; Pierre, O. A Three-Dimensional Analytical Solution for the Study of Air Pollutant Dispersion in a Finite Layer. Boundary-Layer Meteorol. 2015, 155, 289–300. [Google Scholar] [CrossRef]
  26. Obukhov, A. Turbulence in an Atmosphere with a Non-uniform Temperature. Bound.-Layer Meteorol. 1971, 2, 7–29. [Google Scholar] [CrossRef]
  27. Gryning, S.E.; Batchvarova, E.; Brümmer, B.; Jørgensen, H.; Larsen, S. On the extension of the wind profile over homogeneous terrain beyond the surface boundary layer. Bound.-Layer Meteorol. 2007, 124, 251–268. [Google Scholar] [CrossRef]
  28. Mooney, C.J.; Wilson, J.D. Disagreements between gradient-diffusion and Lagrangian stochastic dispersion models, even for surface near the ground. Bound.-Layer Meteorol. 1993, 64, 291–296. [Google Scholar] [CrossRef]
  29. Mangia, C.; Moreira, D.M.; Schipa, I.; Degrazia, G.A.; Tirabassi, T.; Rizza, U. Evaluation of a New Eddy Diffusivity Parameterization from Turbulent Eulerian Spectra in Different Stability Conditions. Atmos. Environ. 2002, 36, 67–76. [Google Scholar] [CrossRef]
  30. Nieuwstadt, F.T.M.; van Ulden, A.P. A Numerical Study on the Vertical Dispersion of Passive Contaminants from a Continuous Source in the Atmospheric Boundary Layer. Atmos. Environ. 1978, 12, 2119–2124. [Google Scholar] [CrossRef]
  31. Huang, C. A Theory of Dispersion in Turbulent Shear Flow. Atmos. Environ. (1967) 1979, 13, 453–463. [Google Scholar] [CrossRef]
  32. Doran, J.C.; Horst, T.W. An evaluation of Gaussian plume-depletion models with dual-tracer field measurements. Atmos. Environ. 1985, 19, 939–951. [Google Scholar] [CrossRef]
  33. Doran, J.C.; Abbey, O.B.; Buck, J.W.; Glover, D.W.; Horst, T.W. Field Validation of Exposure Assessment Models; EPA/600/3-84/092A; Data Environmental Science Research Laboratory, U.S. Environmental Protection Agency: Research Triangle Park, NC, USA, 1984; p. 177.
  34. Kadyri, A.; Kandoussi, K.; Souhar, O. An approach on mathematical modeling of PV module with sensitivity analysis: A case study. J. Comput. Electr. 2022, 21, 1365–1372. [Google Scholar] [CrossRef]
  35. Laaouaoucha, D.; Farhane, M.; Essaouini, M.; Souhar, O. Analytical Model for the Two-Dimensional Advection-Diffusion Equation with the Logarithmic Wind Profile in Unstable Conditions. Int. J. Environ. Sci. Technol. 2022, 19, 6825–6832. [Google Scholar] [CrossRef]
  36. Moreira, D.M.; Tirabassi, T.; Vilhena, M.T.; Goulart, A.G. A Multi-layer Model for Pollutant Dispersion with Dry Deposition to the Ground. Atmos. Environ. 2010, 44, 1859–1865. [Google Scholar] [CrossRef]
  37. Ermak, D.L. An analytical model for air pollutant transport and deposition from a point source. Atmos. Environ. (1967) 1977, 11, 231–237. [Google Scholar] [CrossRef]
Figure 1. A representation of the vertical subdivision of the atmospheric boundary layer into H sub-layers and the continuity of concentration and flux at the interface of each sub-layer.
Figure 1. A representation of the vertical subdivision of the atmospheric boundary layer into H sub-layers and the continuity of concentration and flux at the interface of each sub-layer.
Atmosphere 15 01054 g001
Figure 2. Comparative analysis of predicted versus observed concentrations under different deposition conditions, illustrating the model’s performance across a range of environmental settings. (a) Scatter plot of predicted vs. observed concentrations without deposition using a logarithmic scale. (b) Scatter plot of predicted vs. observed concentrations without deposition using a linear scale. (c) Scatter plot of predicted vs. observed concentrations with deposition using a logarithmic scale. (d) Scatter plot of predicted vs. observed concentrations with deposition using a linear scale.
Figure 2. Comparative analysis of predicted versus observed concentrations under different deposition conditions, illustrating the model’s performance across a range of environmental settings. (a) Scatter plot of predicted vs. observed concentrations without deposition using a logarithmic scale. (b) Scatter plot of predicted vs. observed concentrations without deposition using a linear scale. (c) Scatter plot of predicted vs. observed concentrations with deposition using a logarithmic scale. (d) Scatter plot of predicted vs. observed concentrations with deposition using a linear scale.
Atmosphere 15 01054 g002
Figure 3. Discretization of the ABL into two, four, and six sub-layers, illustrating improved adherence to model parametrizations and a closer approximation to experimental measurements with increased discretization. (a) Discretization of the ABL into two sub-layers. (b) Discretization of the ABL into four sub-layers. (c) Discretization of the ABL into six sub-layers.
Figure 3. Discretization of the ABL into two, four, and six sub-layers, illustrating improved adherence to model parametrizations and a closer approximation to experimental measurements with increased discretization. (a) Discretization of the ABL into two sub-layers. (b) Discretization of the ABL into four sub-layers. (c) Discretization of the ABL into six sub-layers.
Atmosphere 15 01054 g003
Figure 4. Analysis of the concentration’s convergence in a discretized ABL model across different distances from the source at H s = 3 cm, demonstrating the impact of varying the deposition velocity V g . (a) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0 . (b) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0.01 . (c) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0.02 . (d) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0.03 .
Figure 4. Analysis of the concentration’s convergence in a discretized ABL model across different distances from the source at H s = 3 cm, demonstrating the impact of varying the deposition velocity V g . (a) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0 . (b) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0.01 . (c) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0.02 . (d) Convergence at x = 800 , x = 1600 , and x = 3200 m for V g = 0.03 .
Atmosphere 15 01054 g004
Table 1. Micrometeorological data and observed crosswind-integrated concentrations c y / Q ( 10 3 sm 2 ) from the Hanford diffusion experiment.
Table 1. Micrometeorological data and observed crosswind-integrated concentrations c y / Q ( 10 3 sm 2 ) from the Hanford diffusion experiment.
RunDate L MO u H mix Distance c SF 6 y / Q c ZnS y / Q V g ( ZnS )
(DD/MM/YY)( m )( ms 1 )( m )( m )( 10 3 sm 2 )( 10 3 sm 2 )( 10 2 ms 1 )
8003.732.244.21
118 May 19831650.4032516002.140.984.05
32001.300.573.65
80012.907.471.93
226 May 1983440.2613516009.083.251.80
32007.222.311.74
8005.913.063.14
305 June 1983770.2718216003.311.323.02
32001.790.662.84
80020.108.041.75
412 June 1983340.20104160013.104.261.62
32009.153.141.31
80010.505.251.56
524 June 1983590.2615716008.613.381.47
32006.642.921.14
80013.407.231.17
627 June 1983710.3018516006.152.521.15
32003.111.251.10
Table 2. Statistical measures achieved with models with and without deposition.
Table 2. Statistical measures achieved with models with and without deposition.
Model without Deposition ( V g = 0 )
ModelsNMSEMRSECORFBFSMGVGFAC2
( H = 6 ) 0.00591 0.00591 0.94745 0.014 0.00710 0.99044 1.0157 1
( H = 4 ) 0.00990 0.00989 0.93503 0.04653 0.00961 0.93302 1.0302 1
( H = 2 ) 0.01142 0.01141 0.93413 0.05485 0.00969 0.91911 1.036 1
Ref. [25] ( H = 4 ) 0.08 - 0.86 0.02 0.05 -- 1.0
Ref. [25] ( H = 2 ) 0.1 - 0.82 0.1 0.04 -- 0.92
Ref. [35] ( H = 1 ) 0.069 -- 0.009 0.051 0.996 1.055 1.009
Model with Deposition ( V g 0 )
ModelsNMSEMRSECORFBFSMGVGFAC2
( H = 6 ) 0.01746 0.01746 0.94619 0.00732 0.00877 1.0221 1.0388 1
( H = 4 ) 0.02003 0.02003 0.92379 0.02882 0.01897 0.98231 1.0411 1
( H = 2 ) 0.05188 0.05179 0.89382 0.08679 0.01440 0.91068 1.0962 0.94444
Ref. [36] ( H = 1 ) 0.09 - 0.95 0.22 0.08 --1
Ref. [24] ( H = 1 ) 0.088 - 0.96 0.196 0.32 --1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Farhane, M.; Souhar, O. Generalized 3D Model of Crosswind Concentrations and Deposition in the Atmospheric Boundary Layer. Atmosphere 2024, 15, 1054. https://doi.org/10.3390/atmos15091054

AMA Style

Farhane M, Souhar O. Generalized 3D Model of Crosswind Concentrations and Deposition in the Atmospheric Boundary Layer. Atmosphere. 2024; 15(9):1054. https://doi.org/10.3390/atmos15091054

Chicago/Turabian Style

Farhane, Mehdi, and Otmane Souhar. 2024. "Generalized 3D Model of Crosswind Concentrations and Deposition in the Atmospheric Boundary Layer" Atmosphere 15, no. 9: 1054. https://doi.org/10.3390/atmos15091054

APA Style

Farhane, M., & Souhar, O. (2024). Generalized 3D Model of Crosswind Concentrations and Deposition in the Atmospheric Boundary Layer. Atmosphere, 15(9), 1054. https://doi.org/10.3390/atmos15091054

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