Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

`plumax` equation registry

Offline bookkeeping index of every labelled equation across the roadmap

UNEP
IMEO
MARS

plumax — equation registry

Developer aid. This file is a bookkeeping index of every labelled equation across the plumax roadmap notes. It exists so that authors can sanity-check label uniqueness, plan future cross-references, and audit equation numbering at a glance.

Scope: all {math} blocks under roadmap/. Convention: every block label is eq-<page>-<short-name>. <page> matches the file basename (tier1, tier2, …) or sub-page tag (va, vb, vc, vd, rtm, prereqs). Total: 77 labelled equations across 11 pages.


00 — Prerequisites — roadmap/00_prerequisites.md

LabelFormSection
eq-ak-operatorŷ = A (h^T x + (1−h^T) x_a)Averaging-kernel operator

I — Tier I (Gaussian family) — roadmap/01_tier1_gaussian.md

LabelFormSection
eq-tier1-gaussian-plumeSteady-state Gaussian plume with PBL image-source sumGaussian plume
eq-tier1-gaussian-puffSum of 3D Gaussian puffs advected by the windGaussian puff
eq-tier1-column-aky_model = A (∫c dz + c_bg) — column + AK forwardColumn + AK
eq-tier1-likelihoodHeteroscedastic Gaussian observation likelihoodLikelihood
eq-tier1-amortizedNPE / amortized predictor signatureAmortized inference
eq-tier1-pbl-wellmixedVertically well-mixed limit c → Q / (ū L √(2π) σ_y) exp(−y²/2σ_y²)Validation — PBL capping

II — Tier II (Lagrangian particle) — roadmap/02_tier2_lagrangian.md

LabelFormSection
eq-tier2-langevinMarkov-1 Langevin SDE on particle velocity + positionLangevin dynamics
eq-tier2-dt-boundAdaptive timestep bound from τL\tau_L, CFL, diffusionTime-stepping
eq-tier2-footprintBackward-mode source–receptor footprint integralFootprint definition
eq-tier2-forwardy = A col_z(F q) + c_bg + εForward observation operator
eq-tier2-likelihoodε ~ N(0, R_retr + R_repr)Likelihood
eq-tier2-posterior-meanLognormal posterior mean — Gaussian–Gaussian linearisedClosed form
eq-tier2-posterior-covLognormal posterior covariancesame
eq-tier2-amortizedPredictor signature (y, met) → p(log q(x))Amortized inference

III — Tier III (Eulerian FV) — roadmap/03_tier3_eulerian.md

LabelFormSection
eq-tier3-conservationρ\rho-weighted advection–diffusion–source–sinkConservation eq.
eq-tier3-forwardy_t = A_t col_z(c(S, c₀, t)) + c_bg + ε_tForward operator
eq-tier3-4dvar-cost3-term 4D-Var cost (source + IC + obs sum)4D-Var cost
eq-tier3-likelihoodPer-time heteroscedastic GaussianLikelihood
eq-tier3-adjointConservative adjoint transport equationAdjoint
eq-tier3-incrementalIncremental 4D-Var outer-iterate updateIncremental 4D-Var
eq-tier3-control-xformχ=B1/2(SSb)\boldsymbol{\chi} = \mathbf{B}^{-1/2}(S − S_b) control transformControl transform
eq-tier3-amortizedSequence predictor signatureAmortized inference

RTMRTM stack — roadmap/04_rtm_stack.md

LabelFormSection
eq-rtm-radianceSWIR Beer–Lambert radianceBeer–Lambert
eq-rtm-tau-totalTwo-way airmass-factor optical depthsame
eq-rtm-tauLayer-integrated τ\tau from HAPI cross-sectionsame
eq-rtm-thermal-irTIR variant (emission + path radiance)Thermal-IR addendum
eq-rtm-state-vectorJoint retrieval state vectorJoint state
eq-rtm-gauss-newtonIterative Gauss–Newton updateGauss–Newton
eq-rtm-gainGain matrix Gk\mathbf{G}_k + Jacobian definitionsame
eq-rtm-convergenceRodgers (2000) §5.7 convergence criterionsame
eq-rtm-posterior-covS\mathbf{S}^* posterior covarianceInfo content
eq-rtm-averaging-kernelA=GK\mathbf{A} = \mathbf{G}\mathbf{K}, DOFs = trsame
eq-rtm-info-contentShannon HH + posterior contraction ΔS\Delta\mathbf{S}same
eq-rtm-factorised-lutFactorised LUT (gas × surf + scatt) decompositionFactorised LUT
eq-rtm-jax-gradjax.grad ‖NeuralRTM(x) − y‖² end-to-end gradientEmulator-based inference
eq-rtm-amortizedDirect-retrieval predictor signatureAmortized inference

IV — Tier IV (Coupled E2E) — roadmap/05_tier4_coupled.md

LabelFormSection
eq-tier4-per-inst-forwardPer-instrument coupled forwardPer-instrument forward
eq-tier4-per-inst-likelihoodPer-instrument noise covariance decompositionsame
eq-tier4-state-vectorFull coupled state vector enumerationState vector
eq-tier4-fused-operatorJoint multi-instrument observation operatorMulti-instrument fusion
eq-tier4-q-ouOU process for Q(t)Q(t)Q(t)Q(t)
eq-tier4-gradientEnd-to-end multi-instrument gradientGradient
eq-tier4-cost3-term coupled 4D-Var cost (per-inst obs + prior + Q(t)Q(t) kernel)Cost function
eq-tier4-coupled-emulatorCoupled emulator gϕg_\phi signatureCoupled emulator
eq-tier4-amortizedOperational predictor signatureAmortized inference

V — Tier V index — roadmap/06_tier5_population.md

LabelFormSection
eq-tier5-loglikTMTPP three-term log-likelihoodLog-likelihood
eq-tier5-markPer-event mark integral PdLifdQ\int P_d L_i f \, \mathrm{d}QMark contribution
eq-tier5-importance-weightsImportance-weighted MC estimator f/πper-eventf / \pi_\text{per-event}same

V.A — Instantaneous emission — roadmap/06a_instantaneous.md

LabelFormSection
eq-va-full-posteriorPer-event posterior Liπper-event\propto L_i \cdot \pi_\text{per-event}Full-physics posterior
eq-va-wind-rescaleWind-source consistency rescalingCatalog Q
eq-va-mark-integralMark-likelihood integral formMark likelihood
eq-va-importance-mcImportance-weighted MC estimatorsame
eq-va-thinned-rateIntegrated thinned-rate termDetection-floor / non-detection
eq-va-multisourceMulti-source product over within-overpass sourcesMulti-source

V.B — TMTPProadmap/06b_point_process.md

LabelFormSection
eq-vb-lambda-constantHomogeneous PoissonTemporal
eq-vb-lambda-diurnalDiurnal sinusoidal intensitysame
eq-vb-lambda-stepStep intensity (valve fail)same
eq-vb-lambda-decayExponential decay (blowdown)same
eq-vb-lambda-hawkesHawkes / self-exciting kernelsame
eq-vb-lambda-lgcpLog-Gaussian Cox processsame
eq-vb-pod-hillHill function Pd(Q)=1/(1+(Q50/Q)k)P_d(Q) = 1/(1 + (Q_{50}/Q)^k)POD
eq-vb-pod-hier-priorHierarchical prior on (Q50,k)(Q_{50}, k)POD calibration
eq-vb-loglikCanonical TMTPP log-likelihoodLikelihood — canonical
eq-vb-mark-iwImportance-weighted MC for the mark integralPractical evaluation
eq-vb-loglik-pointPoint-regime simplificationPoint regime

V.C — Persistency — roadmap/06c_persistency.md

LabelFormSection
eq-vc-wait-homogeneousE[Δt]=1/λ0\mathbb{E}[\Delta t] = 1/\lambda_0Wait time
eq-vc-wait-inhomogeneousInhomogeneous-Poisson wait integralsame
eq-vc-occurrence-homogeneous1eλ0Δt1 - e^{-\lambda_0 \Delta t}Occurrence
eq-vc-occurrence-inhomogeneous1eλdt1 - e^{-\int \lambda \, \mathrm{d}t}same
eq-vc-hawkes-bumpConditional intensity given prior detectionConditional intensity
eq-vc-cumulativeE[N(0,T)]=0Tλdt\mathbb{E}[N(0,T)] = \int_0^T \lambda \, \mathrm{d}tCumulative count

V.D — Total emission — roadmap/06d_total_emission.md

LabelFormSection
eq-vd-mtotalMtotal(T)=E[Ntrue]E[Q]M_\text{total}(T) = \mathbb{E}[N_\text{true}] \cdot \mathbb{E}[Q]Corrected estimator
eq-vd-mnaiveBiased naive estimator (sum over detected)same
eq-vd-ndetectedE[Ndetected]\mathbb{E}[N_\text{detected}] formulasame
eq-vd-mtotal-posteriorPosterior samples for MtotalM_\text{total}Posterior
eq-vd-pod-unionUnion POD across KK satellitesMulti-satellite fusion

Maintenance

When you add or remove a {math} block:

  1. Add / remove the corresponding row in this file.

  2. Re-run the count grep -rE "^:label: eq-" projects/plume_simulation/notes/ | wc -l and update the total at the top.

  3. If a label changes, grep for {eq} references across the roadmap (grep -rE "\{eq\}\`" projects/plume_simulation/notes/`) and rename them too.