Sedimentology: Texture and Statistics, Fluid Flow and Structures, Sedimentary Rocks, Facies and Cycles, Provenance and Basins
2. Fluid flow, entrainment, settling and bedforms
A fluid has density ρ and dynamic viscosity μ (Pa·s); kinematic viscosity is ν = μ/ρ. Water at about 20 °C has ρ ≈ 1000 kg/m³ and μ ≈ 1.0 × 10⁻³ Pa·s. Water and air are Newtonian; debris flows and muds are non-Newtonian (Bingham: they have a yield strength). The Reynolds number Re = ρUh/μ = Uh/ν compares inertial and viscous forces; in open channels (h = flow depth) flow is laminar below about 500 and turbulent above about 2000, and natural flows are almost always turbulent (U = 1 m/s, h = 1 m gives Re ≈ 10⁶). The Froude number Fr = U/√(gh) compares flow speed with the speed √(gh) of a shallow-water wave: Fr < 1 is subcritical (waves can travel upstream), Fr > 1 supercritical, and the change from supercritical to subcritical is a hydraulic jump.
For steady uniform flow the bed shear stress is τ₀ = ρgh·S (S = slope) and the shear velocity is u* = √(τ₀/ρ). Grains move when the Shields parameter θ = τ₀/[(ρs − ρ)gD] exceeds a critical value θc, about 0.03–0.06 for sand and gravel in turbulent flow. The Hjulström diagram shows the entrainment velocity lowest for sand of roughly 0.1–1 mm; finer silt and clay need higher velocities because cohesion binds them. Grains travel by traction, saltation or suspension, and stay suspended when turbulent eddy velocities, of order u*, exceed the settling velocity.
Stokes' law gives the settling velocity of a small smooth sphere in laminar flow: vs = (ρs − ρf) g d² / (18 μ). It holds for particle Reynolds number Re_p = ρf vs d/μ below about 1 (quartz in water finer than about 0.1 mm). Because vs varies as d², doubling the diameter quadruples the velocity. A 2 µm clay particle needs about three days to fall 1 m, so clay deposits only in still water or after flocculation (clumping in saline water). Coarser grains follow a drag law in which vs rises only as √d, and angular or platy grains settle slower than spheres.
| Bedform | Flow regime | Structure preserved |
|---|---|---|
| Current ripples | Lower | Small cross-lamination, in sand finer than about 0.6 mm |
| Dunes | Lower | Planar or trough cross-bedding; foresets dip downcurrent |
| Upper plane bed | Upper | Parallel lamination with parting lineation |
| Antidunes | Upper, Fr near or above 1 | Low-angle cross-stratification, rarely preserved |
| Wave ripples, hummocky bedforms | Oscillatory or combined flow | Symmetrical ripples; hummocky cross-stratification in storm beds |
At a fixed depth, with sand of about 0.3 mm, raising the velocity gives no motion, then ripples, dunes, upper plane bed and finally antidunes.
3. Sedimentary structures and sediment gravity flows
| Structure | Process | What it tells |
|---|---|---|
| Graded bedding | Normal grading from a waning flow; inverse grading in grain flows | Way-up and flow type |
| Cross-bedding | Migration of ripples and dunes | Palaeocurrent; bipolar (herringbone) means tides |
| Flaser, wavy, lenticular bedding | Alternating sand ripples and mud settling | Tidal flats; flaser is sandy, lenticular is muddy |
| Flute casts, groove casts | Scours in mud filled by sand; flute steep end faces upcurrent | Flow direction (grooves give axis only) |
| Mud cracks | Desiccation | Subaerial exposure |
| Load casts, convolute lamination | Unequal loading, dewatering of wet sediment | Rapid deposition, soft-sediment deformation |
Sediment gravity flows are driven by their own weight. Debris flows are cohesive, matrix-supported and poorly sorted, with floating clasts, ungraded or inversely graded; grain flows are supported by grain collisions on steep slopes; turbidity currents are turbulence-supported suspensions that deposit graded turbidites. Debris flows are Bingham fluids, turbidity currents behave as Newtonian density currents.
The Bouma sequence has five divisions. Ta is massive or graded sand; Tb parallel lamination (upper-regime plane bed); Tc ripple and convolute lamination (lower regime); Td upper parallel lamination of fine sand and silt; Te mud (pelite and hemipelagite). It records a decelerating flow, so proximal beds are thick, coarse and Ta-rich, while distal beds are thin, fine and begin at Tc or Td, with the lower divisions missing. Flute and groove marks usually occur on the base. Bed thickness is the sum of the divisions.
4. Petrology of common sedimentary rocks
Gravel rocks are conglomerate (rounded clasts) and breccia (angular clasts); a clast-supported rock is an orthoconglomerate and a matrix-supported one a paraconglomerate, and diamictite is a poorly sorted mixture, glacial or from mass flow. Sandstones are classified from the framework grains quartz (Q), feldspar (F) and lithic fragments (L), recalculated to 100, and the matrix (grains finer than about 0.03 mm): in Dott's scheme arenites have under 15% matrix, wackes 15–75%, and more than 75% is mudrock.
| Rock | Framework | Meaning |
|---|---|---|
| Quartz arenite | Quartz above 95% | Very mature, recycled or long worked, as on stable cratons |
| Subarkose, sublitharenite | Quartz 75–95%; feldspar over lithics, or lithics over feldspar | Moderately mature |
| Arkose | Quartz below 75%, feldspar over lithics | First-cycle granite or gneiss source, rapid erosion and burial |
| Litharenite | Quartz below 75%, lithics over feldspar | Orogenic source of volcanic, sedimentary or metamorphic fragments |
| Greywacke (wacke) | Matrix over 15%, poorly sorted | Immature; commonly turbidity-current deposits |
Mudrocks (silt and clay) are shale if fissile and mudstone if massive; clay minerals (kaolinite, illite, smectite, chlorite) record weathering and diagenesis, and black shales hold organic matter preserved in oxygen-poor bottom water. Carbonates are named by Dunham's textures: mudstone (under 10% grains, mud-supported), wackestone (over 10% grains, mud-supported), packstone (grain-supported with mud), grainstone (grain-supported, no mud, so high energy) and boundstone (organically bound, as in reefs). Folk names the allochems (ooids, peloids, bioclasts, intraclasts) and whether they sit in micrite (lime mud) or sparite (cement). Aragonite and high-Mg calcite are metastable and invert to low-Mg calcite, and dolomitisation replaces Ca by Mg.
Chemical and biochemical rocks include evaporites (gypsum, anhydrite, halite) from evaporated brine in arid basins; chert (microcrystalline silica, often from radiolarian or diatom ooze); banded iron formation (alternating iron-oxide and silica bands, mainly Precambrian); phosphorite; and coal, from buried plant matter (peat, lignite, bituminous coal, anthracite).
5. Facies, environments and cyclicity
A facies is a rock body distinguished by lithology, structures and fossils, and a facies association is the set of facies that occur together in one environment. Walther's law: facies that lie side by side in a continuous environment also follow each other vertically, provided there is no break such as an unconformity. A vertical profile is therefore read as a lateral model.
| Environment | Typical vertical profile | Diagnostic features |
|---|---|---|
| Meandering river | Fining upward: channel lag, cross-bedded point-bar sand, floodplain mud | Erosive base, lateral accretion, palaeosols |
| Braided river | Coarse sand and gravel, little mud, stacked channel fills | Weak fining, sheet-like bodies |
| Delta | Coarsening upward: prodelta mud, delta-front sand, mouth bar, delta plain | Progradation, Gilbert-type foresets |
| Tidal flat | Fining upward from sand flat to mud flat | Flaser, wavy, lenticular bedding; herringbone cross-bedding |
| Turbidite fan | Thickening-upward lobes; thinning-upward channel fills | Bouma beds, sole marks |
| Aeolian dune | Large-scale cross-bedding in thick sets | Well-sorted, rounded, frosted sand |
| Glacial | Diamictite, laminated mud with dropstones | Striated clasts, varves |
Cyclicity. Autocyclic cycles come from processes inside the system (channel avulsion, delta-lobe switching, tidal-flat progradation) and are local and irregular. Allocyclic cycles are driven from outside (tectonic pulses, climate, eustatic sea level) and correlate across basins; cyclothems (Pennsylvanian coal measures) alternate marine and non-marine strata with sea level. Milankovitch orbital cycles pace allocyclic cycles: precession about 19–23 kyr, obliquity about 41 kyr, eccentricity about 100 kyr (and 400 kyr). The test is bundling: about five precession cycles per eccentricity cycle (5:1). The time per cycle is the interval divided by the number of cycles.
6. Provenance and basin analysis
Provenance asks what and where the source was. The tools are framework petrography on Dickinson's QFL triangles (continental block: quartz-rich; magmatic arc: feldspar and volcanic lithics; recycled orogen: quartz with sedimentary and metamorphic lithics and chert), heavy minerals (zircon, tourmaline and rutile are stable; garnet, staurolite and kyanite come from metamorphic sources; chromite from ultramafic rocks), element ratios such as Th/Sc and Zr/Sc, and palaeocurrents. The chemical index of alteration, CIA = 100 × Al₂O₃/(Al₂O₃ + CaO* + Na₂O + K₂O) in moles (CaO* is silicate-bound CaO), is about 50 for fresh upper crust and rises towards 100 as feldspar alters to kaolinite. Detrital zircon U–Pb ages give a spectrum of source ages, and the youngest cluster sets a maximum depositional age.
| Basin | Mechanism | Fill and Indian example |
|---|---|---|
| Rift | Fault-bounded extension, half-grabens | Alluvial, lacustrine; Gondwana coal basins of India |
| Passive margin | Thermal subsidence after rifting | Wedge thickening seawards |
| Foreland | Flexure under thrust-sheet load | Asymmetric, deepest at the orogen; Siwalik molasse |
| Pull-apart | Releasing bend on a strike-slip fault | Narrow, deep, rapidly filled |
Subsidence arithmetic. Sedimentation rate is thickness divided by time (1 m/Myr = 1 mm/kyr), and rates measured over shorter intervals look higher because gaps count (the Sadler effect). Sediment shrinks on burial as porosity falls, so thickness must be decompacted. By Airy isostasy, for a basin filled to sea level with sediment of mean density ρs, the tectonic subsidence is Y = S(ρm − ρs)/(ρm − ρw), where S is the sediment thickness and ρm and ρw the mantle and water densities. Backstripping removes the sediment load to recover the tectonic part, and after rifting the thermal subsidence decays roughly exponentially with time (McKenzie model).
Key takeaways
- φ = −log₂ d (mm), larger φ is finer. Folk and Ward: Mz = (φ16 + φ50 + φ84)/3 and σI = (φ84 − φ16)/4 + (φ95 − φ5)/6.6, with σI below 0.35 very well sorted. Positive skewness in φ means excess fines. Palaeocurrents need the vector mean, never the arithmetic mean.
- Re = ρUh/μ separates laminar from turbulent flow; Fr = U/√(gh) separates subcritical from supercritical. Grains move when the Shields parameter exceeds about 0.03–0.06. Stokes: vs = (ρs − ρf)gd²/(18μ), valid for Re_p below about 1, and vs ∝ d².
- Raising velocity gives ripples, dunes, upper plane bed, antidunes. Bouma Ta to Te records a waning turbidity current; a matrix-supported, ungraded bed is a debris flow. Flutes show direction, grooves only the axis, herringbone cross-bedding means tides.
- Sandstones: quartz arenite above 95% quartz, arkose feldspar-rich, litharenite lithic-rich, wacke over 15% matrix. Dunham: mudstone, wackestone, packstone, grainstone, boundstone. Evaporites, chert, banded iron formation and coal are chemical or biochemical rocks.
- Walther's law links lateral and vertical facies if there is no break. Rivers fine upward, deltas and shorelines coarsen upward. Autocyclic cycles are local, allocyclic ones regional; Milankovitch periods are about 19–23, 41 and 100 kyr, with five precession cycles per eccentricity cycle.
- Provenance uses QFL, heavy minerals, geochemistry, CIA and detrital zircon ages, whose youngest cluster gives a maximum depositional age. Rift, passive-margin, foreland and pull-apart basins differ in mechanism. Airy subsidence for a basin filled to sea level is Y = S(ρm − ρs)/(ρm − ρw).
Practice questions (18)
Attempt each one before opening the answer. Every explanation names the tempting wrong option as well as the right one, because that is where marks are lost.
The median grain size of a sand is φ50 = 2.5. What is its median diameter and its Udden–Wentworth class?
Show answer
Answer: B — 0.18 mm, fine sand, since d = 2^(−2.5) mm lies between 0.125 and 0.25 mm
φ = −log₂ d, so d = 2^(−φ) = 2^(−2.5) = 0.177 mm, which lies between 0.125 and 0.25 mm, the fine-sand class. Very fine sand is 0.0625–0.125 mm, so the option that gives the right diameter with the finer class is the tempting slip. Using d = 2^φ inverts the scale and gives a pebble, and 2^(−1.5) would be φ = 1.5.Which of the following statements are correct?
Show answer
Answer: A — Larger φ values correspond to smaller grains; B — A positive Folk and Ward skewness indicates a tail of fine grains
φ increases as grain size decreases, so positive skewness in φ units means excess fines. σI below 0.35 is very well sorted, so 0.3 is the opposite of poor sorting. Azimuths are circular: 350° and 10° have a vector mean of 0° (north), and the arithmetic value of 180° points the opposite way.From the cumulative curve of a sand, φ5 = 0.5, φ16 = 1.2, φ50 = 2.1, φ84 = 3.3 and φ95 = 4.1. Use the Folk and Ward inclusive graphic standard deviation σI = (φ84 − φ16)/4 + (φ95 − φ5)/6.6 to find the sorting value. Give the answer to two decimal places.
Numerical answer — type the value.
Show answer
Answer: 1.07
(φ84 − φ16)/4 = (3.3 − 1.2)/4 = 2.1/4 = 0.525 and (φ95 − φ5)/6.6 = (4.1 − 0.5)/6.6 = 3.6/6.6 = 0.5455. Sum σI = 1.0705, which rounds to 1.07, so the sand is poorly sorted (1.00–2.00). The median φ50 is not used in σI.Cross-bed azimuths, already corrected for tilt, measure 350°, 10°, 30° and 50°. Use the vector mean θ̄ = atan2(Σ sin θ, Σ cos θ) to find the mean palaeocurrent direction, in degrees to the nearest whole degree, within the range 0 to 360.
Numerical answer — type the value.
Show answer
Answer: 20
Σ sin θ = −0.1736 + 0.1736 + 0.5000 + 0.7660 = 1.2660 and Σ cos θ = 0.9848 + 0.9848 + 0.8660 + 0.6428 = 3.4784. θ̄ = atan(1.2660/3.4784) = 20.0°. The readings are symmetric about 20° (−30°, −10°, +10°, +30°). The arithmetic mean, (350 + 10 + 30 + 50)/4 = 110°, is wrong because azimuths are circular.A river 2.0 m deep flows at a mean speed of 4.0 m/s. Take g = 9.8 m/s² and ν = 1.0 × 10⁻⁶ m²/s, and use Re = Uh/ν and Fr = U/√(gh). Which description is correct?
Show answer
Answer: A — Subcritical (Fr ≈ 0.90) and turbulent (Re ≈ 8 × 10⁶)
Fr = 4.0/√(9.8 × 2.0) = 4.0/4.43 = 0.90, below 1, so the flow is subcritical. Re = 4.0 × 2.0/10⁻⁶ = 8 × 10⁶, far above 2000, so it is turbulent. A fast river is not necessarily supercritical, which is the tempting error: the criterion compares speed with the wave speed √(gh) = 4.43 m/s. Fr ≈ 1.8 would need a depth of only 0.5 m.Use Stokes' law, vs = (ρs − ρf) g d²/(18 μ). A quartz grain with ρs = 2650 kg/m³ and diameter 0.05 mm settles in water with ρf = 1000 kg/m³ and μ = 1.0 × 10⁻³ Pa·s. Take g = 9.8 m/s². What is its settling velocity in mm/s? Give the answer to two decimal places.
Numerical answer — type the value.
Show answer
Answer: 2.25
d = 0.05 mm = 5 × 10⁻⁵ m, so d² = 2.5 × 10⁻⁹ m². vs = (1650)(9.8)(2.5 × 10⁻⁹)/(18 × 10⁻³) = 4.0425 × 10⁻⁵/0.018 = 2.246 × 10⁻³ m/s = 2.25 mm/s. The particle Reynolds number is about 0.11, so Stokes' law is valid.In the Stokes regime the settling velocity of a sphere, with density and viscosity fixed, is proportional to
Show answer
Answer: C — the square of the diameter, d²
Stokes' law is vs = (ρs − ρf) g d²/(18 μ), so vs ∝ d²: doubling the diameter quadruples the velocity. The √d dependence belongs to the turbulent drag regime of coarser grains, and d³ is the dependence of the grain's weight, not of its velocity.In a flume with sand of about 0.3 mm at constant depth, the flow velocity is raised steadily from below the threshold of motion until the Froude number exceeds 1. In which order do the bedforms appear?
Show answer
Answer: D — ripples, dunes, upper plane bed, antidunes
With increasing velocity and stress, fine to medium sand passes through the lower flow regime (ripples, then dunes) to the upper regime (plane bed, then antidunes at Fr near 1). The other orders put upper-regime forms before lower-regime ones, or reverse ripples and dunes.A 25 cm thick sandstone bed has an erosional base with flute casts. From the base upward it shows parallel lamination, ripple cross-lamination, thin parallel-laminated silt and then mud, with no massive graded division at the base. Which interpretation is best?
Show answer
Answer: B — A turbidity-current deposit with Ta missing, from a decelerating and more distal flow
Flute casts with a sequence of Tb (parallel lamination), Tc (ripples), Td (laminated silt) and Te (mud) is a Bouma sequence lacking Ta, as expected for a waning, more distal turbidity current. A debris flow is matrix-supported and ungraded without internal lamination, which is the tempting alternative because of the thickness and erosional base. Wave deposits do not carry flute casts or a full upward-fining turbidite succession, and a grain flow shows no such lamination series.A complete turbidite bed is 120 cm thick. Its Ta division is 45 cm thick, Tb is 25 cm, Tc is 20 cm and Td is 10 cm, and the rest of the bed is Te mud. A section contains 18 such beds. What is the total thickness of Te mud in the 18 beds, in metres? Give the answer to one decimal place.
Numerical answer — type the value.
Show answer
Answer: 3.6
Te = 120 − (45 + 25 + 20 + 10) = 120 − 100 = 20 cm = 0.20 m per bed. For 18 beds the total is 18 × 0.20 = 3.6 m.A sandstone has framework grains of quartz 88%, feldspar 7% and lithic fragments 5%, and 4% of the rock is clay matrix. In the Pettijohn–Dott scheme it is a
Show answer
Answer: C — subarkose, because it is an arenite with 75–95% quartz and more feldspar than lithics
Matrix of 4% is below 15%, so the rock is an arenite. Quartz of 88% lies between 75 and 95%, and feldspar (7%) exceeds lithics (5%), giving a subarkose. Quartz arenite needs more than 95% quartz and is the tempting answer for a quartz-rich rock. Arkose needs quartz below 75%, and a wacke needs more than 15% matrix.A limestone consists of ooids that touch one another, with lime mud filling part of the space between them. In Dunham's classification it is a
Show answer
Answer: A — packstone, because the grains support each other and lime mud is present between them
Grains in mutual contact make the rock grain-supported, and the lime mud between them separates a packstone from a grainstone, which has no mud. A wackestone is mud-supported, with grains floating in the mud, which does not fit touching ooids. Boundstone requires organic binding, as in reefs, not simply cement.Which of the following statements about facies and cycles are correct?
Show answer
Answer: A — Walther's law links lateral facies to vertical order only where no break, such as an unconformity, intervenes; B — A fining-upward sequence of channel lag, cross-bedded sand and mud is typical of a meandering-river point bar
Walther's law holds only for conformable successions, and a point bar fines upward from lag to floodplain mud. A coarsening-upward profile is the signature of progradation, as in deltas and shorelines. Cycles driven by orbital forcing are allocyclic, the external kind, and are the ones that correlate between basins; autocyclic cycles are local.In a carbonate succession, thin cycles that each represent about 20 kyr of precession occur in bundles of five. A bundle most likely records
Show answer
Answer: D — about 100 kyr of short eccentricity, giving the 5:1 bundling ratio
Five cycles of 20 kyr make 100 kyr, the period of short eccentricity, and the 5:1 bundling of precession cycles within eccentricity cycles is the diagnostic Milankovitch signature. A 400 kyr cycle would need about twenty precession cycles per bundle, and obliquity has a period of about 41 kyr, not 100.The youngest concordant detrital zircon U–Pb ages in a sandstone cluster at 520 Ma, with no younger grains. What does this tell us?
Show answer
Answer: B — The sandstone cannot be older than about 520 Ma, so 520 Ma is a maximum depositional age
Detrital zircons crystallised in the source rocks and were eroded and deposited later, so the rock cannot be older than its youngest grain: 520 Ma is a maximum age, and deposition may be much younger. Treating the age as the time of deposition is the tempting error. It is not a minimum age, and zircon records igneous as well as metamorphic crystallisation.A basin is filled to sea level with 3000 m of sediment of mean density 2400 kg/m³. Use Airy isostasy, Y = S(ρm − ρs)/(ρm − ρw), where S is the sediment thickness, ρm = 3300 kg/m³ is the mantle density and ρw = 1000 kg/m³ is the water density. What is the tectonic subsidence Y in metres, to the nearest whole metre?
Numerical answer — type the value.
Show answer
Answer: 1174
Y = 3000 × (3300 − 2400)/(3300 − 1000) = 3000 × 900/2300 = 1173.9 m, which rounds to 1174 m. Only about 39% of the sediment thickness is tectonic subsidence; the rest is accommodated by the isostatic response to the sediment load.Which of the following statements about sedimentary basins are correct?
Show answer
Answer: A — A foreland basin is an asymmetric wedge, deepest next to the thrust belt, where the lithosphere flexes under the thrust load; B — A passive-margin basin subsides mainly by thermal contraction after an earlier stage of fault-controlled stretching
Flexure under the thrust load makes the foreland basin deepest next to the orogen, and a passive margin subsides thermally after rifting. Pull-apart basins open at releasing bends, where extension thins the crust, not at restraining bends, which build uplift.A river 1.5 m deep flows over a bed of 2 mm quartz gravel with a water-surface slope of 0.0002. Take ρ = 1000 kg/m³, ρs = 2650 kg/m³, g = 9.8 m/s², τ₀ = ρgh·S, a Shields parameter θ = τ₀/[(ρs − ρ)gD] and a critical value θc = 0.045. Will the gravel move?
Show answer
Answer: C — Yes, because τ₀ ≈ 2.9 Pa gives θ ≈ 0.09, above the critical value
τ₀ = 1000 × 9.8 × 1.5 × 0.0002 = 2.94 Pa. θ = 2.94/(1650 × 9.8 × 0.002) = 2.94/32.34 = 0.091, which is greater than 0.045, so the gravel moves. A subcritical Froude number does not prevent entrainment, and the entrainment criterion is the Shields threshold, not simply turbulence. The option that gives θ ≈ 0.03 uses the same inputs but misplaces the arithmetic.