Are these realistic daisies?

No. The model is deliberately minimal: no latitude, no clouds, no soil chemistry, and no evolution — only the qualitative coupling between life and reflectivity.

Why can T(L) differ from the bare-rock curve?

Because ⟨α⟩ depends on b and w, while b and w depend on local temperatures shifted by their albedos. That closed loop can buffer T against changes in L compared with a constant albedo.

Daisyworld (Climate Regulation)

Daisyworld couples a zero-dimensional Stefan–Boltzmann balance σT⁴ = (1−⟨α⟩) L S₀/4 to two daisy species whose cover fractions evolve on available bare ground. Each species has its own albedo and a simple parabolic growth rate versus local patch temperature T_i = T + q(⟨α⟩−α_i), so black daisies are warmer than the mean and white daisies cooler. Death provides turnover. Sweeping luminosity L reveals how the inhabited planet can track a narrower surface-temperature band than the same planet with fixed bare-soil albedo alone — the pedagogical “regulation” effect.

Who it's for: Introductory biogeoscience, astrobiology, or nonlinear dynamics modules.

Key terms

Daisyworld
Gaia hypothesis
Albedo feedback
Homeostasis
Luminosity

A minimal homeworld lesson: life can couple to albedo and thermostatically soften climate sensitivity to stellar output — still a cartoon, not a biosphere model.

Live graphs

Gray: equilibrium temperature with fixed bare albedo only. Green: after settling daisies at each L (continuation upward in L). Regulation narrows the T(L) curve vs the lifeless planet.

Planet & luminosity

Solar luminosity factor L1 ×

Albedo α (black daisies)0.25

Albedo α (white daisies)0.75

Albedo α (bare ground)0.5

Global balance σT⁴ = (1−⟨α⟩)L S₀/4 with ⟨α⟩ = b α_b + w α_w + (1−b−w) α_g. Growth uses local patch temperatures T_i = T + q(⟨α⟩−α_i): black daisies warm relative to the mean, white daisies cool relative to it; bare area 1−b−w is colonized.

Growth & death

Black growth β_b1.1

White growth β_w1.1

Death rate γ0.35

Integration speed1.5 ×

Comfort windows (K)

Black T low265 K

Black T high285 K

White T low278 K

White T high302 K

Show T(L) panel

Measured values

Surface T-39.11 °C

Black patch T-34.1 °C

White patch T-44.1 °C

Mean ⟨α⟩0.500

Black cover b8.0 %

White cover w8.0 %

Bare ground84.0 %

T if only bare soil-39.11 °C

How it works

Watson–Lovelock-style Daisyworld: black and white daisies modify planetary albedo and buffer global temperature as luminosity changes.

Frequently asked questions

Are these realistic daisies?: No. The model is deliberately minimal: no latitude, no clouds, no soil chemistry, and no evolution — only the qualitative coupling between life and reflectivity.
Why can T(L) differ from the bare-rock curve?: Because ⟨α⟩ depends on b and w, while b and w depend on local temperatures shifted by their albedos. That closed loop can buffer T against changes in L compared with a constant albedo.

Other simulators in this category — or see all 32.

View category →

NewSchool

ENSO Delayed Oscillator (Toy)

Scalar delay differential equation ẋ = a·x(t−τ) − b·x − c·x³ + seasonal forcing: explore El Niño– / La Niña–like swings vs damping in a Suarez–Schopf–style cartoon.

Launch Simulator

NewUniversity / research

Rossby Waves on a Beta-Plane

Linear planetary wave ψ = A cos(kx+ly−ωt) with ω = −βk/(k²+l²+Rd⁻²): streamfunction bands and geostrophic arrows show westward phase propagation.

Launch Simulator

NewUniversity / research

Geostrophic Balance & Thermal Wind

Pressure-gradient force balanced by Coriolis: isobars set geostrophic wind, while meridional temperature gradients produce vertical shear between 1000 and 300 hPa.

Launch Simulator

NewUniversity / research

Groundwater Contaminant Plume

Advection-dispersion pulse in groundwater with longitudinal/transverse spreading, retardation factor R, and a monitoring-well breakthrough curve.

Launch Simulator

NewUniversity / research

Infinite Slope Stability

Limit-equilibrium factor of safety for a shallow planar slide: slope angle, cohesion, friction, depth, and rainfall-driven pore pressure set the landslide threshold.

Launch Simulator

NewUniversity / research

Earthquake Aftershocks: Omori + Gutenberg-Richter

Modified Omori aftershock decay n(t)=K/(t+c)^p combined with Gutenberg-Richter magnitude-frequency curves, b-value, and a synthetic catalog.

Launch Simulator

Daisyworld (Climate Regulation)

Live graphs

How it works

Frequently asked questions

More from Biophysics, Fluids & Geoscience

ENSO Delayed Oscillator (Toy)

Rossby Waves on a Beta-Plane

Geostrophic Balance & Thermal Wind

Groundwater Contaminant Plume

Infinite Slope Stability

Earthquake Aftershocks: Omori + Gutenberg-Richter

Daisyworld (Climate Regulation)

Live graphs

How it works

Frequently asked questions