How Earth's Hadley Cells Work: An Interactive Atmospheric Circulation Laboratory
MERIDIAN is an original, exploratory three-dimensional model of Earth's tropical atmospheric circulation. Move the season, strengthen the imposed temperature contrast, or change the chosen cell width, then observe how simulated meridional flow, the rising branch near the intertropical convergence zone (ITCZ), and solar energy gradients respond. The aim is not to reproduce today's weather. It is to connect the physics of solar geometry and large-scale atmospheric overturning to visually interpretable patterns.
1. Why does the tropical atmosphere circulate?
The tropics receive more solar energy, on average over a year, than the polar regions. Solar radiation warms the surface and the lower atmosphere. Warm, moisture-bearing tropical air can rise in regions where the atmospheric column is unstable and surface winds converge. As the rising air cools, water vapor can condense and form deep clouds. Aloft, air can move toward the subtropics. In an idealized Hadley circulation, the air eventually sinks, warms by compression, and returns toward the equator near the surface.
This circulation transfers energy and moisture, shapes trade winds and tropical rain belts, and contributes to the relatively dry climate often found in subtropical regions. Real-world Hadley cells are not perfectly symmetric, continuous tubes. They vary with longitude, season, land–sea contrasts, topography, weather systems, and the distribution of clouds and oceans. MERIDIAN intentionally suppresses those complications to make the large-scale mechanism visible.
2. Reading the 3D planet
The central ocean globe has a latitude–longitude reference grid. Amber flow trajectories mark the idealized Northern Hemisphere Hadley cell; cyan trajectories mark the Southern Hemisphere cell. A bright equatorial ring denotes the illustrative latitude of concentrated tropical ascent. Pale cloud markers identify convective regions and higher-latitude cloud belts. The flow ribbons rise outward from the globe and sweep meridionally to show a vertical cross-section repeated around multiple longitudes.
The simulated troposphere is assigned a reference depth of 16 km. Earth's radius is 6,371 km, but the atmospheric height in the visualization is deliberately exaggerated by about 53 times so students can actually see the loops. A displayed path height is therefore not a literal proportion of Earth's radius. Model positions represent a conceptual spherical diagram, and the background ocean texture is decorative: it is not a geographic land map, temperature observation, or satellite product.
3. Try the season control
The solar declination is the latitude directly under the Sun at local solar noon. As Earth orbits the Sun, axial tilt causes the declination to migrate between approximately 23.44° north and 23.44° south. The simulator uses the convenient approximation
declination = 23.44° × sin[2π(day_of_year − 80) / 365.2422].
At the March equinox, the declination is close to the equator. Near the June solstice, it approaches 23.44° N. Near the December solstice, it approaches 23.44° S. Adjust the day of year and watch the daily mean top-of-atmosphere (TOA) solar flux change by latitude. The latitude of the ITCZ also shifts in this model, but by a smaller amount than the declination. This smaller migration is a chosen teaching approximation rather than a universal prediction; the real ITCZ responds to ocean heat storage, regional monsoons, cross-equatorial energy transport, and many other influences.
4. Calculating daily solar energy
The energy curve uses the standard spherical-astronomy equation for daily mean TOA insolation. With latitude φ, solar declination δ, solar constant S₀ and the sunset hour angle H₀:
H₀ = arccos[clamp(−tan(φ) tan(δ), −1, 1)]
Q = S₀ / π × [H₀ sin(φ) sin(δ) + cos(φ) cos(δ) sin(H₀)].
This lab sets S₀ = 1361 W/m², a rounded modern nominal value. At an equinox, the equatorial value is approximately S₀/π ≈ 433 W/m². This is daily average incoming solar power at the top of the atmosphere, not the temperature of the surface, the midday solar maximum, or the amount absorbed by the land or ocean. Reflection by clouds, atmospheric absorption, day-to-day weather and planetary heat transport are not included in that number. The formula naturally handles polar night and midnight Sun through the sunset-hour-angle clamp.
5. Flow in the cell: a conservation-based sketch
MERIDIAN uses a two-dimensional latitude–altitude streamfunction for each of the two hemispheric overturning cells, then projects many copies around a spherical globe. Define x as normalized distance from the illustrative ITCZ to the selected poleward cell edge, z as normalized altitude in the 16 km troposphere, V₀ as a reference meridional velocity and L as cell half-width in metres. A simplified streamfunction is
ψ(x,z) = −(V₀ H / π) sin(πx) sin(πz).
The associated meridional and vertical velocities satisfy the two-dimensional incompressibility condition by construction:
v = ∂ψ/∂height = −V₀ sin(πx) cos(πz);
w = −∂ψ/∂distance = V₀ (H/L) cos(πx) sin(πz).
The lower branch returns equatorward, air rises near the ITCZ, the upper branch returns poleward, and the circulation descends near the subtropical edge. This is a kinematic teaching model. It does not resolve moisture thermodynamics, density stratification, eddy momentum exchange, turbulent convection, radiation balance, mass-weighted pressure coordinates, or the full fluid equations of atmospheric motion. The displayed peak vertical velocity is a scale implied by the chosen V₀, troposphere depth and cell width; it must not be mistaken for a measured convective updraft speed.
6. Understanding trade winds and upper-level winds
Air moving toward the equator near the surface is deflected by Earth's rotation. In a rotating frame, the Coriolis effect turns Northern Hemisphere equatorward flows toward the west, helping create northeasterly trade winds; Southern Hemisphere equatorward flows also acquire a westward component, creating southeasterly trades. Poleward transport aloft tends to develop an eastward component, and the subtropical jet is associated with that momentum evolution and atmospheric temperature gradients.
A useful limiting calculation is axial angular-momentum conservation. If frictionless air originates at latitude φ₀ and subsequently moves to latitude φ, the idealized zonal wind relative to Earth is
u = Ωa [cos²(φ₀)/cos(φ) − cos(φ)],
where Ω = 7.2921159 × 10⁻⁵ s⁻¹ and a = 6,371,000 m. This idealization becomes unrealistic at sufficiently high latitude because eddies, pressure gradients and mixing strongly affect actual winds. The model shows an intentionally capped illustrative wind scale in its diagnostic calculation rather than claiming to solve a realistic jet-stream forecast.
7. Why the air sinks near the subtropics
Persistent rising air in tropical convective regions must be accompanied by compensating motion elsewhere if mass is to be conserved. In the textbook Hadley-cell picture, air transported poleward in the upper troposphere can descend into broad subtropical high-pressure regions, where it warms adiabatically and tends to suppress cloud formation. This does not mean that all subtropical latitudes are deserts, or that there is always sinking motion at exactly 30 degrees. Ocean currents, mountains, monsoons, and seasonal weather change local outcomes considerably.
8. Suggested hands-on experiments
Start with day 80, a 30° half-width and 1.00× differential heating. Compare the daily solar flux at 25° N and 25° S: near the equinox they should be similar. Then select day 172. The subsolar latitude moves north; the curve shows stronger summertime insolation in the Northern Hemisphere. Observe the approximate ITCZ migrate north in the 3D schematic. Repeat for day 355, and compare the reversal.
Next, increase differential heating from 0.50× to 1.50× while holding the day and width fixed. The flow speed parameter scales as the square root of this factor in the simple model, which changes the vertical velocity scale. Finally, widen the cell from 18° to 38°. The same chosen meridional flow has to span a longer horizontal distance: the vertical velocity scale V₀H/L becomes smaller. This illustrates an important dimensional relationship, not a validated prediction that a warmer Earth must produce narrower or wider Hadley cells.
Use the play/pause, drag-to-orbit, zoom, flow, cloud and grid controls to isolate visual structures. The screenshot function lets educators capture frames. The lower graphs present seasonal solar energy latitude profiles and a conceptual meridional cross-section, which makes comparison easier than a 3D perspective alone.
9. What the experiment cannot tell you
The simulator has no observational assimilation, grid-resolved weather, land–sea temperature contrasts, clouds that feed back radiatively, interactive ocean dynamics, water-vapor transport, or real-time climate information. It is deliberately axisymmetric in its core equations. It should not be used to predict drought, storms, regional rainfall, future warming or local wind speeds. The moving particles are animated representations of flow trajectories and are not tracked air parcels with resolved density, moisture or turbulence.
10. Educational outcomes
After using this lab, learners should be able to explain the rising and sinking branches of a Hadley cell; identify the location and approximate meaning of the ITCZ; distinguish solar declination from surface temperature; identify trade-wind directions in both hemispheres; explain why Earth's rotation influences upper-level winds; and interpret the effect of changing cell width in a continuity-based illustrative model. More advanced learners can reimplement the insolation and streamfunction equations, test mass-continuity properties, and compare qualitative outcomes against an observational reanalysis dataset.
Frequently asked questions
1. What is a Hadley cell?
A Hadley cell is a large-scale tropical atmospheric overturning circulation with rising air in low latitudes, poleward transport aloft, sinking motion in subtropical regions and an equatorward low-level return flow.
2. Why are there two Hadley cells?
The classic schematic has a Northern and a Southern Hemisphere circulation, one on each side of the equatorial rain belt. In reality, their sizes, strengths and cross-equatorial connections change with season and location.
3. What does ITCZ mean?
The Intertropical Convergence Zone is a broad zone where low-level tropical winds converge, encouraging rising air, clouds and precipitation. Its latitude is highly variable around the globe.
4. How does Earth’s tilt change the simulation?
Axial tilt shifts the latitude receiving the strongest midday sunlight throughout the year, changing the daily mean energy profile. MERIDIAN also shifts its illustrative ITCZ proxy with solar declination.
5. Are the displayed wind speeds measured on Earth?
No. They are values from a transparent, simplified kinematic model. They are not station measurements, forecasts or outputs of a full general circulation model.
6. Why does the air sink around 30 degrees latitude?
In a simplified model, air rising in the tropics must eventually descend elsewhere to conserve mass. Broader atmospheric dynamics often place descending branches in subtropical regions, though they vary considerably.
7. Does stronger heating always make a real Hadley cell stronger?
No general conclusion follows from that control. Its square-root speed scaling is an explicit teaching choice, not a universal climate sensitivity law.
8. Why is the atmosphere drawn so thick?
A literal 16 km troposphere would be too thin to see on a globe representing a 6,371 km radius Earth. The visual altitude is exaggerated by about 53 times and is labeled accordingly.
9. Can I run the experience without an internet connection?
The viewer includes a native WebGL2 fallback which can run without a JavaScript CDN when opened through a local HTTP server. The enhanced Three.js engine imports a pinned library from a public CDN and therefore requires that dependency to be reachable.
10. Is this a map of observed climate or cloud cover?
No. Globe coloring, particles and cloud clusters are original illustrative artwork. None are georeferenced measurements or satellite-derived records.
Model metadata: Earth radius 6,371 km; idealized tropopause 16 km; solar constant 1361 W/m²; rotation rate 7.2921159 × 10⁻⁵ rad/s; original random generation seed 654105. Authorship: Original procedural educational visualization for MySimulator.uk, 2026-11-04.