Coupled Boundary-Layer Science

This page describes the physical equations, scaling formulations, and validation studies of Shardian Atmos's boundary-layer scheme (EML-SR).

1. Monin-Obukhov Similarity Theory (MOST)

In the surface layer (lowest 10% of the planetary boundary layer), turbulent fluxes of momentum and heat are parameterized using Monin-Obukhov similarity relations:

\[u_* = \frac{\kappa u}{\ln\left(\frac{z}{z_{0m}}\right) - \psi_m\left(\frac{z}{L}\right)}\]

\[\theta_* = \frac{\kappa \Delta\theta}{\ln\left(\frac{z}{z_{0h}}\right) - \psi_h\left(\frac{z}{L}\right)}\]

Where: * \(u_*\) is the friction velocity (momentum scale). * \(\theta_*\) is the temperature scale (heat flux scale). * \(\kappa \approx 0.4\) is the von Kármán constant. * \(L\) is the Obukhov length, representing the height where shear production equals buoyant production of turbulent kinetic energy. * \(\psi_m\) and \(\psi_h\) are the stability correction functions for momentum and heat.

2. The Shardian EML-SR Coupling Model

Traditional models use empirical relationships (such as Dyer-Businger or Beljaars) to define \(\psi_m\) and \(\psi_h\). These formulations break down in stable nocturnal regimes (\(z/L > 0\)), leading to excessive heat loss or runaway cooling.

Shardian Atmos replaces empirical correction terms with a non-linear symbolic model discovered by symbolic regression (EML-SR). The scheme adaptively calculates the surface exchange coefficients (\(C_m\) and \(C_h\)):

\[C_m = \left(\frac{u_*}{u}\right)^2, \quad C_h = \frac{u_* \theta_*}{u \Delta\theta}\]

Adaptive Canopy Clamping

In grid cells with a high vegetation cover fraction (\(\text{VEGFRA} > 0.7\)), momentum extraction is regulated using a canopy-clamping factor to prevent artificial velocity deceleration near the surface:

\[z_{0m,\text{clamped}} = \max\left(z_{0m,\text{min}}, \; z_{0m} \cdot \exp\left(-\beta_{\text{canopy}} \cdot \text{VEGFRA}\right)\right)\]

This bounds the shear stress, providing realistic wind speeds at turbine hub heights (80m to 120m).

3. Validation over Global Datasets

Shardian Atmos was validated using five global meteorological datasets representing diverse boundary layer regimes:

graph TD
    BOMEX[BOMEX: Marine Boundary Layer] --> A[EML-SR Validation]
    BUBBLE[BUBBLE: High-Roughness Urban Layer] --> A
    COASTAL[COASTAL: Land-Water Transition] --> A
    GABLS1[GABLS1: Highly Stable Polar Boundary Layer] --> A
    WANGARA[WANGARA: Semi-Arid Diurnal Boundary Layer] --> A

Key Scientific Findings:

GABLS1 (Stable Boundary Layer): Traditional schemes decouple the boundary layer, producing a flat surface temperature profile. EML-SR maintains a weak but continuous turbulent mixing rate, reducing nocturnal cold bias by 51.3%.
BUBBLE (Urban Canopy): By adjusting the roughness length scales over highly uneven topography, Shardian Atmos resolves the urban heat island effect with a 36.5% lower RMSE compared to MOST.
WANGARA (Diurnal Cycle): Resolves the transition from unstable daytime convective boundary layers to stable nighttime layers without numerical oscillations.