DIMENSIONLESS BANKFULL HYDRAULIC RELATIONS FOR EARTH AND TITAN

DIMENSIONLESS BANKFULL HYDRAULIC RELATIONS FOR EARTH AND TITAN European Space Agency Gary Parker Dept. of Civil & Environmental Engineering and Dept. of Geology University of Illinois 1

UNTIL RECENTLY TITAN WAS SHROUDED IN MYSTERY What we knew or could reasonably infer: 1. Larger than Mercury 2. Atmospheric pressure ~ 1. 5 Earth atmospheres near surface 3. ~ 95 K near surface 4. Atmosphere of nitrogen (mostly), methane, ethane 5. Crustal material of water/ice 6. Near triple point of methane/ethane: possibility of a. methane/ethane oceans b. methane/ethane precipitation as liquid/solid 7. Possibility of rivers of liquid methane carrying sediment of solid water ice! But a thick shroud of smog produced by the breakdown of methane under ultraviolet light prevented any surface visualization. 2

AND THEN JANUARY 14, 2005 ARRIVED! This and other images of Titan courtesy European Space Agency and NASA Cassini/Huygens Mission: very strong evidence for rivers of liquid methane carrying sediment of water ice I was glued to the internet! I had waited for years! 3

MARS VERSUS TITAN Mars shows evidence of ancient rivers of flowing water that carried sediment similar to that 4 of the Earth’s crust.

MARS VERSUS TITAN contd. But the era of flowing rivers was a long time ago, as evidenced by the fairly intense impact cratering of Mars, and may not has lasted very 5 long as compared to Earth.

MARS VERSUS TITAN contd. Tectonic ridges? Titan shows evidence of active tectonics, vulcanism, aeolian and fluvial reworking, and has very few impact craters: so its surface is likely active 6 in modern geological time!

MARS VERSUS TITAN contd. Volcano? Titan shows evidence of active tectonics, vulcanism, aeolian and fluvial reworking, and has very few impact craters: so its surface is likely active 7 in modern geological time!

MARS VERSUS TITAN contd. Aeolian dunes? Titan shows evidence of active tectonics, vulcanism, aeolian and fluvial reworking, and has very few impact craters: so its surface is likely active 8 in modern geological time!

MARS VERSUS TITAN contd. River drainage basin? Titan shows evidence of active tectonics, vulcanism, aeolian and fluvial reworking, and has very few impact craters: so its surface is likely active 9 in modern geological time!

MARS VERSUS TITAN contd. Impact crater Titan shows evidence of active tectonics, vulcanism, aeolian and fluvial reworking, and has very few impact craters: so its surface is likely active 10 in modern geological time!

ALLUVIAL GRAVEL-BED RIVERS ON TITAN? The evidence suggests that at least near where Huygens touched down, there is a plethora of alluvium in the gravel and sizes. The gravel presumably consists of water ice and appears to be fluvially rounded. 11

CAN OUR KNOWLEDGE OF ALLUVIAL GRAVEL-BED RIVERS ON EARTH HELP US MAKE INFERENCE ABOUT TITAN? 12

IF WE KNEW THE PHYSICS BEHIND RELATIONS FOR BANKFULL GEOMETRY HERE ON EARTH • Bankfull Depth Hbf ~ (Qbf)0. 4 Bbf ~ (Qbf)0. 5 • Bankfull Width • Bed Slope S ~ (Qbf)-0. 3 where Qbf = bankfull discharge we might be able to extend the relations to Titan. 13

WE BEGIN WITH EARTH The Parameters: Qbf Qb. T, bf Bbf Hbf S D s R = = = = = g = = bankfull discharge (m 3/s) volume bedload transport rate at bankfull discharge (m 3/s) bankfull width (m) bankfull depth (m) bed slope (1) surface geometric mean or median grain size (m) density of water (kg/m 3) density of sediment (kg/m 3) ( s/ ) – 1 = submerged specific gravity of sediment ~ 1. 65 (1) gravitational acceleration (m/s 2) kinematic viscosity of water (m 2/s) The forms sought: dimensionless versions of Why dimensionless? In order to allow scaling between Earth and Titan! 14

Meet my friends the DIMENSIONLESS PARAMETERS Particle Reynolds number Dimensionless bankfull discharge Dimensionless bankfull depth Dimensionless bankfull width Down-channel bed slope Dimensionless bankfull bedload transport rate Bankfull Shields number at threshold of motion 15

DATA SETS FOR GRAVEL-BED RIVERS ON EARTH 1. Alberta streams, Canada 1 2. Britain streams (mostly Wales)2 3. Idaho streams, USA 3 4. Colorado River, USA (reach averages) Kellerhals, R. , Neill, C. R. and Bray, D. I. , 1972, Hydraulic and geomorphic characteristics of rivers in Alberta, River Engineering and Surface Hydrology Report, Research Council of Alberta, Canada, No. 72 -1. 2 Charlton, F. G. , Brown, P. M. and Benson, R. W. , 1978, The hydraulic geometry of some gravel rivers in Britain, Report INT 180, Hydraulics Research Station, Wallingford, England, 48 p. 3 Parker, G. , Toro-Escobar, C. M. , Ramey, M. and Beck S. , 2003, The effect of floodwater extraction on the morphology of mountain streams, Journal of Hydraulic Engineering, 129(11), 2003. 4 Pitlick, J. and Cress, R. , 2002, Downstream changes in the channel of a large gravel bed river, Water Resources Research 38(10), 1216, doi: 10. 1029/2001 WR 000898, 2002. 1 16

WHAT THE DATA SAY: WIDTH, DEPTH, SLOPE The four independent sets of data form a coherent set! 17

REGRESSION RELATIONS BASED ON THE DATA To a high degree of approximation, Remarkable, no? 18

WHAT DOES THIS MEAN? 19

WHAT THE DATA SAY: BANKFULL SHIELDS NUMBER 20

THE PHYSICS BEHIND IT ALL Assume the following relations. Manning-Strickler resistance relation Parker-Einstein bedload relation Relation for bankfull Shields number Channel form relation of type of Parker (1978) “Gravel yield” relation 21

THE RELATIONS OF THE PREVIOUS SLIDE YIELD PRECISELY THE OBSERVED DIMENSIONLESS RELATIONS! 22

GENERALIZATION FOR OTHER PLANETS/SATELLITES Manning-Strickler resistance relation Parker-Einstein bedload relation The presence of g and R allow us to go from Relation for bankfull Shields number Channel form relation of type of Parker (1978) to “Gravel yield” relation (volume to mass) 23

BACK-CALCULATED DIMENSIONALLY HOMOGENEOUS BANKFULL HYDRAULIC RELATIONS FOR ALLUVIAL GRAVEL RIVERS ON S ARBITRARY HEAVENLY BODIE The presence of g and R allow us to go from to 24

FROM Pressure E-atmo Temperature K Grav. accel. m/s 2 Fluid dens. kg/m 3 Sed. Dens. kg/m 3 ( s/ ) - 1 Kin. Viscosity m 2/s TO Parameter p T g Earth 1 ~ 293 9. 81 Titan 1. 5 ~ 95 1. 40 s R 1000 2650 1. 65 1. 00 x 10 -6 446 931 1. 09 4. 04 x 10 -7 25

CONSIDER A STREAM WITH THE SAME BANKFULL DISCHARGE Qbf AND CHARACTERISTIC GRAIN SIZE D HOW SHOULD TITAN COMPARE WITH EARTH? E = Earth, T = Titan From to 26

CONSIDER A STREAM WITH THE SAME BANKFULL DISCHARGE Qbf AND CHARACTERISTIC GRAIN SIZE D HOW SHOULD TITAN COMPARE WITH EARTH? E = Earth, T = Titan = 1. 48 x 0. 83 = 1. 23 = 1. 57 x 1. 56 = 2. 46 = 0. 72 x 0. 80 = 0. 57 27

SO FOR THE SAME BANKFULL DISCHARGE Qbf AND CHARACTERISTIC GRAIN SIZE D A gravel-bed river on might be 1. 23 x the bankfull depth, 2. 46 x the bankfull width and 0. 57 x the down-channel slope of a gravel-bed river on Could braiding be more common on Titan? 28

BUT WAIT A MINUTE! IS “GRAVEL” ON TITAN GRAVEL ON EARTH? For dynamic similarity in grain Reynolds number or or So the answer is “yes” to a reasonable approximation! 29

GRAIN REYNOLDS INVARIANCE Besides, the dynamics of sediment transport becomes approximately invariant to particle Reynolds number for or D >~ 8. 8 mm on Earth or D >~ 10. 6 mm on Titan based on the condition c*/ c, asymp* 0. 90 using 30

WHAT ABOUT AEOLIAN PROCESSES ON TITAN? Let Ua = wind velocity, a = atmospheric density, Cf = drag coefficient, s = sediment density, D = grain size. Scaling for mobility of grain size D: Atmospheric density Earth 293 K 1 E-atmo, a = 1. 21 kg/m 3 Titan (nitrogen) 95 K 1. 5 E-atmo, a = 5. 39 kg/m 3 Assuming Reynolds invariance (Cf constant), critical velocity Uac to blow around size D scales as: Much easier to blow sediment around on Titan! But much less solar heating to drive meteorology! 31

QUESTIONS OR COMMENTS? 32