Dry Deposition

Model

This is an implementation of a box model used to calculate changes in gas species concentration due to dry deposition.

Running the model

Here's an example of how concentration of different species, such as SO₂, O₃, NO₂, NO, H₂O₂ and CH₂O change due to dry deposition.

We can create an instance of the model in the following manner:

using AtmosphericDeposition
using ModelingToolkit
using DifferentialEquations
using EarthSciMLBase
using Unitful

@parameters t [unit = u"s", description="Time"]
model = DrydepositionG(t)

Before running any simulations with the model we need to convert it into a system of differential equations.

sys = structural_simplify(model)
tspan = (0.0, 3600*24) 
u0 = [2.0,10.0,5,5,2.34,0.15] # initial concentrations of SO₂, O₃, NO₂, NO, H₂O₂, CH₂O
sol = solve(ODEProblem(sys, u0, tspan, []),AutoTsit5(Rosenbrock23()), saveat=10.0) # default parameters

which we can plot as

using Plots
plot(sol, xlabel="Time (second)", ylabel="concentration (ppb)", legend=:outerright)
Example block output

Parameters

The parameters in the model are:

parameters(sys) # [z, z₀, u_star, L, ρA, G, T, θ]

where z is the top of the surface layer [m], z₀ is the roughness length [m], u_star is friction velocity [m/s], and L is Monin-Obukhov length [m], ρA is air density [kg/m3], T is surface air temperature [K], G is solar irradiation [W m-2], Θ is the slope of the local terrain [radians].

Let's run some simulation with different value for parameter z.

@unpack O3 = sys

p1 = [50,0.04,0.44,0,1.2,300,298,0]
p2 = [10,0.04,0.44,0,1.2,300,298,0]
sol1 = solve(ODEProblem(sys, u0, tspan, p1),AutoTsit5(Rosenbrock23()), saveat=10.0)
sol2 = solve(ODEProblem(sys, u0, tspan, p2),AutoTsit5(Rosenbrock23()), saveat=10.0)

plot([sol1[O3],sol2[O3]], label = ["z=50m" "z=10m"], title = "Change of O3 concentration due to dry deposition", xlabel="Time (second)", ylabel="concentration (ppb)")
Example block output