Outer Planets and Ocean Worlds Environment

Beyond the asteroid belt, solar flux and communication rates fall as the inverse square of heliocentric distance, round-trip light time runs to hours, and the two accessible surfaces are a giant planet entry corridor and Titan, a 1.35 m/s2 body [1] under a dense cold nitrogen atmosphere lit at roughly a thousandth of Earth daylight. What follows sets out the available solar flux and communication geometry, the giant planet entry regime as measured by the only probe that has flown it, the Titan atmosphere and surface as measured by the only lander that has reached it, and the reference atmosphere models design work is run against.

Energy and communication geometry

Solar flux falls as the inverse square of heliocentric distance. At Titan the combined effect is larger than distance alone: about 100 times less flux than at Earth from the heliocentric distance, and about 10 times less again through the haze, giving a surface illumination roughly 1000 times below Earth daylight and comparable to about 1000 times full moonlight, concentrated in red and near infrared [1]. Photovoltaics are not a practical primary source for a Titan surface vehicle.

Communication cost scales with distance. Missions using a high gain antenna require empirically about 5 mJ per bit per astronomical unit to acquire and send science data to Earth [1]. The dependence is linear rather than inverse square because antenna aperture has historically grown with distance. The Huygens relay failure and its recovery from the ground illustrate the other consequence: at Saturn range, a probe transmitting through a relay has no margin to spare, and the Doppler Wind Experiment was recovered only because Earth-based radio telescopes could detect the probe directly [14].

Round trip light time is hours, so fault response cannot be supervised from Earth. A Titan surface vehicle experiences a solar day of 384 h, 16 Earth days, and from Titan the Earth is always within 6 degrees of the Sun, so any operation requiring a real-time link is confined to local daytime [1].

Titan

Property	Surface value	Source
Diameter	5150 km	[1]
Surface gravity	1.35 m/s2, 1/7 Earth	[1]
Distance from Saturn	1.2e6 km, 20 Saturn radii	[1]
Rotation and orbital period	15.945 days	[1]
Solar day	384 h	[1]
Surface pressure	1.47 bar	[1]
Surface temperature	94 K	[1]
Surface density	5.4 kg/m3, 4 times Earth sea level air	[1]
Composition	95 percent N2, 5 percent CH4, 0.1 percent H2, trace organics	[1]
Methane mole fraction near the surface	4.5 +/- 0.5 percent	[6]
Speed of sound	195 m/s against 340 m/s on Earth	[1]
Viscosity	6e-6 Pa s, about 3 times lower than Earth air	[1]
Obliquity to the Sun	26 degrees	[1]
Pressure scale height	about 40 km against about 10 km on Earth and Mars	[2]
Visibility near the surface	about 10 km	[1]

Four times the density in one seventh the gravity with three times lower viscosity is the whole aerial mobility argument. An airfoil of a given size and speed operates at several times the Reynolds number it would on Earth, so a 1 m rotor behaves like a much larger terrestrial one and a wind turbine blade section is an appropriate choice [1]. Flight power for a given mass is about 40 times lower than on Earth. The low speed of sound is the limit at the other end: tip Mach number constrains rotor size and speed. Multirotors in this regime remain susceptible to vortex ring state, which is entered when descent rate approaches the rotor wake speed and which is avoided by capping vertical descent rate or keeping forward speed up [12].

Winds and turbulence

Near-surface winds are predicted by global circulation models at 1 to 2 m/s maximum [1], and the Huygens measurements agree: below 2 m/s in the lowest kilometer and generally below 1 m/s [5]. Surface wind from parachute tracking was 0.3 +/- 0.1 to 0.4 +/- 0.1 m/s, with an upper limit of 0.25 m/s within 1 m of the ground derived from probe cooling [4].

Shear rather than mean wind is the design driver. The strongest large-scale shear Huygens encountered was about 5 m/s per km [5]. The engineering specification the probe was designed against gives a 95th percentile wind gradient by altitude:

Altitude	95th percentile gradient	Altitude	95th percentile gradient
below 10 km	3.1 m/s per km	80 km	7.3 m/s per km
20 km	4.6 m/s per km	90 km	6.2 m/s per km
30 km	5.6 m/s per km	100 km	5.6 m/s per km
40 km	6.3 m/s per km	120 km	5.1 m/s per km
50 km	7.1 m/s per km	140 km	4.8 m/s per km
60 km	9.2 m/s per km	160 km	4.5 m/s per km
70 km	9.1 m/s per km

All values from [5]. Measured shear in the lowest 2 km was 1 m/s per km, with a possible near-surface value around 10 m/s per km, and the shear regime changes at about 0.3 km and about 3 km. A specification of this shape, an envelope rather than a profile, is how both wind and density uncertainty are carried into Titan design analysis, because the measurement uncertainty is comparable to the real variability [2].

Mean winds through the descent fell only slightly, from 4.9 m/s at 4 km to 4.5 m/s at the surface [5], against a design envelope that had to accommodate an order of magnitude more shear at altitude.

The boundary layer is shallow. The planetary boundary layer top sits at 300 m, the diurnal layer is about 1 km deep and the seasonal layer a few times thicker [5]. A rotorcraft operating below 1 km therefore spends its whole flight inside the convective layer, which is where the profiling flight profiles are defined [1]. Turbulent fluctuations in the lowest 4 km are about 0.2 m/s: the standard deviation of Doppler wind differences over 3 s intervals is 0.21 m/s, and the tilt sensor standard deviation of 2.5 degrees implies a characteristic fluctuation near 0.1 m/s [5]. Turbulent dissipation rate at about 500 m is 16 cm2/s3 and the convective velocity scale is 0.15 m/s. A first-order autoregression with coefficient 0.96 reproduces the observed wind time series.

Surface mechanics, measured once

Huygens is the only measurement of a Titan surface. It arrived under parachute at 4.5 m/s vertical with about 1 m/s horizontal [4]; the Surface Science Package records a vertical impact velocity of 4.60 m/s [6].

Quantity	Value	Source
Penetration depth	12 cm	[4]
Inferred hole radius	53 cm	[4]
Bounce duration	0.4 s	[4]
Post-bounce speed	0.8 +/- 0.1 m/s horizontal, 0.56 m/s vertical	[4]
Slide distance and duration	30 to 40 cm over about 1.2 s	[4]
Coefficient of friction	0.4 +/- 0.1	[4]
Bearing strength	about 50 N/cm2 (500 kPa)	[4]
Penetrometer trench	1.5 cm wide, 2 to 3 cm deep	[4]
Fluffy surface layer	7 mm thick	[4]
Obstacle encountered	a 1 to 2 cm pebble	[4]
Probe pitch at rest	-3.1 +/- 0.5 degrees	[4]
Inferred local surface slope	2.3 +/- 0.4 degrees	[4]
Inferred surface radius of curvature	6.5 m	[4]

The probe wobbled with a period of about 2.0 s, against a free wobble period of 2.2 s, for five complete cycles with a damping time constant of 1.5 s, and motion remained detectable for 10 s after impact, on material described as fine-grained, damp and relatively soft, behaving like damp sand [4]. A bearing strength of 500 kPa over a fluffy 7 mm layer is a surface a landed vehicle can stand on but not one that resists a point load without sinking through the top centimeter.

Titan entry

Huygens remains the only Titan entry and is therefore the validation case for every model that follows: entry interface was at 1270 km altitude, three parachutes were deployed in sequence, pilot, main and stabilizer, and descent to the surface for a Dragonfly-class vehicle takes about 2.5 hours [7]. The current practice is to re-fly the Huygens trajectory in a six degree of freedom multi-body simulation against the atmosphere model a new mission intends to use, and to reconcile the differences [7], [8].

Giant planet entry

Only one probe has entered a giant planet atmosphere. The Galileo probe entered Jupiter on 7 December 1995 [9].

Quantity	Value	Source
Entry velocity	47.4 km/s	[9]
Entry flight path angle	84 +/- 1 degrees	[9]
Angle of attack	10 degrees	[9]
Spin rate	10.5 rpm	[9]
Peak heating rate	about 30 kW/cm2	[9]
Total heat load	about 300 kJ/cm2	[9]
Peak deceleration	250 g	[9]
Velocity at 450 km altitude	47.406 km/s	[9]
Velocity at 50 km altitude	0.833 km/s	[9]
Time between those two points	111.02 s	[9]

The thermal protection system is the part of that record with no substitute. The forebody nose cap was chopped-molded carbon phenolic, the frustum tape-wrapped carbon phenolic, and the aft heat shield phenolic nylon [9]. Forebody thickness ran from 14.6 cm at the centreline down to 5.1 cm at the front of the frustum, over a 22.2 cm nose radius and a 44.86 degree cone half angle. Measured recession was 4.45 +/- 0.25 cm at the nose, fell to a minimum of 2.74 +/- 0.15 cm at mid-frustum, then rose again toward the base, with a maximum recession rate of 0.4 cm/s [9]. Probe mass fell from 337 +/- 4.0 kg at entry to 248.1 +/- 8.7 kg, of which 79.0 +/- 4.0 kg was forebody ablation, 1.5 +/- 0.5 kg pyrolysis and 8.4 +/- 4.2 kg aft shield loss. The mid-frustum minimum contradicted the design prediction, which is why the experiment exists.

The atmosphere the probe was designed against was not the atmosphere it met. Nominal entry assumed 89 percent H2 and 11 percent He for a molecular weight of 2.22; the measured helium mole fraction was 13.6 +/- 0.3 percent for a molecular weight of 2.35 [9]. The measured He/H ratio of 7.85e-2 +/- 0.16e-2 remains the anchor for every giant planet entry analysis that follows [10].

Jupiter profile at the probe site	Pressure	Temperature	Source
450 km	7.21e-8 bar	527.8 K	[9]
250 km	1.60e-5 bar	163.2 K	[9]
100 km	7.00e-3 bar	160.6 K	[9]
50 km	7.51e-2 bar	112.4 K	[9]

Equilibrium cloud decks sit at 0.7 bar for ammonia ice, 2.2 bar for ammonium hydrosulfide and 5 bar for water ice and droplets [11]. The probe descended into a 5 micrometer hot spot and found water still increasing at 22 bar, with nitrogen and sulfur about 3 times solar at about 10 bar. Juno microwave radiometry of the equatorial zone later gave water at 2.5e3 ppm, 2.7 +/- 1.7 times protosolar, and ammonia at 351 ppm, 2.76 +/- 0.46 times solar [11]. A single entry site does not characterize the planet.

Elemental ratios measured by the probe are the reference set for any future entry payload: He/H 7.85e-2 +/- 0.16e-2, C/H 1.19e-3 +/- 0.29e-3, N/H 3.32e-4 +/- 1.27e-4, O/H 2.45e-4 +/- 0.80e-4, S/H 4.45e-5 +/- 1.05e-5 and Ar/H 9.10e-6 +/- 1.80e-6 [10]. For Saturn only He/H at 6.75e-2 +/- 1.25e-2 and C/H at 2.65e-3 +/- 0.10e-3 are available, from remote sensing.

Reference atmosphere models

Design work runs against the Global Reference Atmospheric Model suite. Each body has its own model, sharing a common software core with the others while keeping body-specific atmosphere data, and each returns mean density, temperature and pressure along a user-defined path selected by time, latitude and longitude [3]. These are engineering models, not forecasts. Their value is dispersion: Monte Carlo perturbation profiles for guidance, navigation and control development and for thermal systems analysis [2].

Titan-GRAM and Neptune-GRAM represent both real variability and the large measurement uncertainty at these bodies inside a single minimum-average-maximum density envelope, with one input parameter selecting where in that envelope a run falls, from -1 to +1 [2]. Titan profiles come from Yelle et al. 1997 and Neptune from Cruikshank 1995.

Source: [2]. Public domain (NASA).

Model	Data source	Known gaps	Source
Uranus-GRAM	NASA Ames Uranus Atmospheric Model, built on Voyager 2 radio science, IRIS and UVS data from the 24 January 1986 flyby	The winds on Uranus are unknown, so the model carries no wind data	[18]
Jupiter-GRAM	Galileo probe Atmospheric Structure Instrument data, after Seiff et al.	No wind model, no perturbation model and no constituent gas model	[3]
Titan-GRAM	Yelle et al. 1997	Envelope only	[2], [15]
Neptune-GRAM	Cruikshank 1995	Envelope only	[2], [16]

The suite has been rearchitected into a common software core, with ephemeris calculations moved to the NAIF SPICE toolkit, built against CSPICE version N0066 [3]. The rearchitecture is a code change rather than a new atmosphere: Neptune-GRAM was first released in September 2004, and the 2020 C++ version carries the same atmosphere data as the 2004 Fortran version, with the changes confined to input parameter names, SPICE ephemeris and output file formats [16]. The Titan-GRAM user guide is the reference for its inputs and outputs [15].

Radiation

The design environment away from a planetary magnetosphere is the interplanetary galactic cosmic ray field, taken from the Badhwar-O’Neill model with solar modulation derived from time-delayed sunspot number [17]. Inside the Jovian magnetosphere the trapped electron environment dominates instead, and no open engineering model of it was retrievable for this entry; the Jovian radiation environment is therefore not quantified here.

Ocean worlds

Europa and Enceladus have no landed measurement of any kind. Reaching a subsurface ocean requires penetrating kilometers of ice inside the Jovian trapped radiation environment, under planetary protection requirements that no viable terrestrial organism be carried. No such mission has been approved and no engineering environment specification comparable to Titan-GRAM exists for either body. Titan is the only outer solar system surface for which in-situ mechanical, thermal and wind data have been returned, and that data set is a single descent and a single landing [4], [5], [13].