Interstellar Horizons

Atmospheric sky color models, compared across the day — Preetham, Nishita, Hosek-Wilkie, and friends

Atmosphere & Sky
Filters

Planetary Presets

Planet Physics
Rayleigh Only Nishita with Mie scattering disabled — pure molecular scattering baseline
Rayleigh Scattering Only
Iλ= Isun ∫0D βR(λ) ρR(h) PR(θ) e-τR ds
Scattering Coefficient
βR∝ 1λ4
No Mie Term
βM=0
Preetham (1999) Analytical — Perez sky luminance/chromaticity model
Perez Sky Function
F(θ,γ) = ( 1+A⁢ eBcosθ ) ( 1+ C⁢eDγ +E⁢cos2γ )
Sky Radiance
Lλ = Lz ⋅ F(θ,γ) F(0,θs)
Zenith Luminance
Yz= (4.0453T-4.971) tanχ -0.2155T+2.4192 χ= ( 49-T120 ) (π-2θs)
Nishita (1993) Ray-marching — single-scattering Rayleigh + Mie
Single-Scattering Integral
Iλ = Isun ∫0D [ βR(λ) ρR(h) PR(θ) + βM ρM(h) PM(θ) ] e -τ(x→,s→) -τ(x→,v→) ds
Rayleigh Phase
PR(θ) = 316π (1+cos2θ)
Henyey-Greenstein Mie Phase
PM(θ) = 3(1-g2) (1+cos2θ) 8π(2+g2) (1+g2-2gcosθ) 32
Density Profiles
ρR(h) = e-h/HR ρM(h) = e-h/HM
Hosek-Wilkie (2012) Enhanced analytical — improved horizon, ground albedo feedback
9-Parameter Sky Radiance
F(θ,γ) = ( 1+Ae Bcosθ+0.01 ) ( C +DeEγ +Fcos2γ +Gχ(H,γ) +Icosθ )
Multiple Scattering
Ims ≈ Iss ⋅ kms (1+ A2 )
Ground Albedo Feedback
Iground = A ⋅ Esun ⋅ sinαs ⋅ e -τ(g→s) ⋅ e -τ(g→v)
Nishita + Ozone Adds Chappuis band absorption — deeper blues during twilight
Total Optical Depth
τ= τR+ τM+ τO3
Ozone Density
ρO3(h) = e -2 |h-25km| 15km
Chappuis Absorption
τO3 = ∫ βO3(λ) ρO3(h) ds
CIE Standard Clear Sky (Type 12) Empirical luminance model — ISO 15469 architectural standard
Sky Luminance
L(θ,γ) Lz = φ(θ)⋅f(γ) φ(0)⋅f(θs)
Gradation
φ(θ)= 1+a eb/cosθ
Indicatrix
f(γ)= 1+c (edγ -edπ/2) +ecos2γ
Parameter Matrix Rayleigh Preetham Nishita Hosek-Wilkie Nishita+O₃ CIE Clear
Geometry
Sun position αs ✓ ✓ ✓ ✓ ✓ ✓
Viewing elevation θ ray direction Perez F(θ,γ)F(0,θs) ray direction ray direction ray direction φ(θ) gradation
Observer altitude h0 ray origin — ray origin ray origin ray origin —
Atmospheric Scattering
Rayleigh coefficient βR sole scatter term implicit in Perez fit β(λ) per bin β(λ) per bin β(λ) per bin —
Mie coefficient βM βM=0 implicit in Perez fit βM∝T βM∝T βM∝T —
Rayleigh scale height HR 8 km — 8 km 8 km 8 km —
Mie scale height HM — — 1.2 km 1.2 km 1.2 km —
Turbidity T scales ρR drives all Perez coeffs scales βM scales βM scales βM affects CCT
Phase Functions
Rayleigh phase PR 316π(1+cos2θ) baked into Perez E 316π(1+cos2θ) 316π(1+cos2θ) 316π(1+cos2θ) —
Mie phase PM — baked into Perez C,D Henyey-Greenstein Cornette-Shanks Henyey-Greenstein —
Asymmetry parameter g — — 0.76 (fixed) 0.76 (fixed) 0.76 (fixed) —
Absorption & Ozone
Ozone absorption βO3 — — — — Chappuis band, σ(λ) —
Ozone density ρO3 — — — — e-2|h-25|15 km —
Ozone column strength SO3 — — — — slider (0–3) —
Surface & Multi-Scatter
Ground albedo A — — — Ig=A⋅E☉e-τ — —
Multiple scattering Ims — implicit in empirical fit — order-2, kms=0.075 — implicit in empirical fit
Color & Tone Mapping
Color space 16-bin spectral → XYZ → RGB CIE Yxy → RGB 16-bin spectral → XYZ → RGB 16-bin spectral → XYZ → RGB 16-bin spectral → XYZ → RGB CCT → Yxy → RGB
Tone mapping k Reinhard, k=1.0 Reinhard, k=0.04 Reinhard, k=1.0 Reinhard, k=1.0 Reinhard, k=1.0 Reinhard, k=0.15
Twilight handling indirect glow fade to -6° indirect glow indirect glow indirect glow fade to -6°
CIE-Specific
Sky type coefficients a–e — — — — — Type 12 (clear, polluted)
Correlated color temp Tcc — — — — — 3500–15000 K
Preetham-Specific
Perez coefficients A–E — 3 sets (Y,x,y)×T — — — —
Zenith chromaticity xz,yz — polynomial in T,θs — — — —
Notes: All models show the sky color at the selected viewing elevation angle. Below the horizon the sun is shown as dark. Hosek-Wilkie uses an analytical approximation of the published model (full dataset is ~100KB of coefficients). Colors are tone-mapped to sRGB and may clip in extreme conditions.