rc_sea_spray

Attributes

A_w

h

U10

C

d_p

rho_p

k_a

Functions

emission_factor(dp_um)

Computes the emission factor (E_f) based on particle diameter (dp) in micrometers. After fitting a power law distribution to empirical data.

settling_velocity(d_p, rho_p[, rho_w, mu, g])

Compute the Stokes settling velocity for a particle in water.

aerosolization_flux(U10, C, d_p, rho_p, E_f_base)

Compute the aerosolization flux incorporating particle density.

aerosolization_rate_constant(U10, C, d_p, rho_p, A_w, h)

Computes the first-order aerosolization rate constant k_a.

Module Contents

rc_sea_spray.emission_factor(dp_um)

Computes the emission factor (E_f) based on particle diameter (dp) in micrometers. After fitting a power law distribution to empirical data.

Parameters:
  • E_f_size_dict – dict A dictionary containing the particle diameter (dp) in micrometers (dp_um) as the key and the corresponding emission factor (E_f) as the value.

  • dp_um – float or array-like Particle diameter in micrometers (µm)

Returns:

float or array-like

Emission factor in m³ m⁻² s⁻¹

Return type:

E_f

rc_sea_spray.settling_velocity(d_p, rho_p, rho_w=1000, mu=0.001002, g=9.81)

Compute the Stokes settling velocity for a particle in water.

Parameters: d_p : Particle diameter (m) rho_p : Particle density (kg/m³) rho_w : Water density (kg/m³) (default = 1000) mu : Dynamic viscosity of water (Pa·s) (default = 1.002e-3) g : Gravitational acceleration (m/s²) (default = 9.81)

Returns: v_s : Settling velocity (m/s)

rc_sea_spray.aerosolization_flux(U10, C, d_p, rho_p, E_f_base)

Compute the aerosolization flux incorporating particle density.

Parameters: U10 : Wind speed at 10 m height (m/s) C : MNP concentration in water (m⁻³) d_p : Particle diameter (m) rho_p : Particle density (kg/m³) E_f_base : Baseline emission factor (m³ m⁻² s⁻¹) from the paper

Returns: F : Aerosolization flux (m⁻² s⁻¹)

rc_sea_spray.aerosolization_rate_constant(U10, C, d_p, rho_p, A_w, h)

Computes the first-order aerosolization rate constant k_a.

Parameters: U10 : Wind speed at 10 m height (m/s) C : MNP concentration in water (m⁻³) d_p : Particle diameter (m) rho_p : Particle density (kg/m³) E_f_base : Baseline emission factor (m³ m⁻² s⁻¹) A_w : Surface area of the water compartment (m²) h : Depth of the water compartment (m)

Returns: k_a : Aerosolization rate constant (s⁻¹)

rc_sea_spray.A_w = 1000000.0
rc_sea_spray.h = 10
rc_sea_spray.U10 = 5
rc_sea_spray.C = 1000000.0
rc_sea_spray.d_p = 1e-06
rc_sea_spray.rho_p = 900
rc_sea_spray.k_a