Environmental Geosciences III: Hydrology and Hydrogeology, the Basins and Groundwater Provinces of India, Darcy's Law, Hydraulic Conductivity, Fluctuations, Tracers, Subsidence, Over-Use and Pollution of Groundwater, the Ghyben-Herzberg Relation, Resource Exploitation and Conservation, and Natural Hazards

This chapter completes Unit IV. It covers hydrology and hydrogeology, the types of aquifer and the basins and groundwater provinces of India; Darcy's law, its validity and the hydraulic conductivity of earth materials, with a discharge and a seepage velocity worked; groundwater fluctuations, tracers, land subsidence, the effects of over-use, groundwater quality and pollution, and the categories used to assess them in India; the Ghyben-Herzberg relation between fresh and saline water, worked; the exploration and exploitation of natural resources and their environmental costs, and the conservation of non-renewable resources; and the catastrophic geological hazards, floods, landslides, earthquakes, volcanism, avalanches, tsunamis and cloudbursts, with their prediction and mitigation.

1. Hydrology, hydrogeology and the groundwater provinces of India

Hydrology studies the occurrence and movement of water at and above the surface; its basic accounting is the water balance of a catchment, P = ET + R + ΔS, precipitation being divided among evapotranspiration, runoff and changes in storage. Hydrogeology studies water below the surface. Below the water table the ground is saturated; above it lies the unsaturated (vadose) zone. An aquifer is a saturated formation that stores and transmits useful quantities of water, such as sand, gravel, fractured rock or karst limestone; an aquitard transmits water slowly, an aquiclude (clay) stores but hardly transmits it, and an aquifuge (unfractured granite) neither stores nor transmits it.

Aquifer typeDescriptionBehaviour of a well
UnconfinedIts upper surface is the water table, open to the atmosphere through the unsaturated zone; recharged directly from aboveWater stands at the level of the water table; vulnerable to surface pollution
ConfinedBounded above and below by aquicludes or aquitards; water under pressure; recharged where the aquifer outcropsWater rises above the top of the aquifer to the piezometric (potentiometric) surface, and flows at the surface in an artesian well if that surface is above ground
PerchedA small saturated lens held above the main water table by a local impermeable layerLimited, seasonal supply

Not all the water in pores can be drained: specific yield is the fraction that drains by gravity and specific retention the fraction held back by surface tension, and together they equal porosity. Clay has a high porosity but a very low specific yield, which is why it makes a poor aquifer. Storativity is the volume released per unit area per unit fall in head, and transmissivity, T = Kb, the product of hydraulic conductivity and saturated thickness, measures how much water the whole aquifer thickness can transmit. India's groundwater occurs in a few broad provinces. The alluvial aquifers of the Indus-Ganga-Brahmaputra plains are thick, extensive and highly productive. The hard-rock terrain of the peninsula, granites, gneisses and the Deccan basalts, which underlies roughly two-thirds of the country, has little primary porosity, so water is held only in the weathered mantle and fractures, and yields are low and quickly exhausted. Coastal sedimentary aquifers are productive but threatened by salinity, and the desert and Himalayan regions have their own limited systems. The major river basins, the Ganga (the largest), the Indus and the Brahmaputra in the north and the Godavari, Krishna, Mahanadi, Cauvery, Narmada and Tapi in the peninsula, frame the surface-water resource.

2. Darcy's law, its validity and hydraulic conductivity

Henry Darcy, experimenting with sand filters for the fountains of Dijon in 1856, found that the discharge through a porous medium is proportional to the hydraulic gradient and the cross-sectional area: Q = K i A, where K is the hydraulic conductivity, i = Δh ÷ L the hydraulic gradient (the fall in head over the flow length) and A the area. The Darcy flux, q = Q ÷ A = K i, has units of velocity but is not the speed at which water moves, because water flows only through the pores; the average linear (seepage) velocity is v = K i ÷ nₑ, where nₑ is the effective porosity, and it is the seepage velocity that sets how fast a pollutant travels. For example, in an aquifer with K = 10 m/day, a head fall of 2 m over 1,000 m (i = 0.002) and a cross-section 50 m wide and 10 m thick (A = 500 m²), Q = 10 × 0.002 × 500 = 10 m³/day; with nₑ = 0.25, the seepage velocity is 10 × 0.002 ÷ 0.25 = 0.08 m/day, so a contaminant takes over 30 years to travel a kilometre.

Darcy's law is valid only for laminar flow, in which viscous forces dominate: in practice where the Reynolds number based on grain size is below about 1 to 10, which covers almost all natural groundwater flow through sands and silts. It breaks down in turbulent flow through coarse gravels, large fractures and karst conduits near wells, and it assumes a saturated, rigid medium; in very fine clays at low gradients the flow may also depart from strict proportionality. Hydraulic conductivity depends on both the medium and the fluid (its density and viscosity): typical values are 10⁻³ to 10⁻¹ m/s for clean gravel, 10⁻⁵ to 10⁻³ m/s for clean sand, 10⁻⁸ to 10⁻⁶ m/s for silt and below 10⁻⁹ m/s for unfractured clay, a range of more than ten orders of magnitude. It is measured in the laboratory with constant-head and falling-head permeameters and in the field by pumping tests, analysed by the Thiem and Theis (1935) methods, and by slug tests.

⚠️ Darcy flux is not the velocity of the water
q = K i is the discharge per unit of total cross-section, including the solid grains. Because water passes only through the connected pores, it actually moves faster, at K i ÷ nₑ. Using q as the travel speed of a contaminant plume underestimates its speed by a factor of 1 ÷ nₑ, typically three to five times.

3. Fluctuations, tracers, subsidence, over-use, quality and pollution of groundwater, and the Ghyben-Herzberg relation

Groundwater levels fluctuate seasonally, rising after monsoon recharge and falling through the dry season as pumping and natural discharge continue, so levels are monitored before and after the monsoon; a long-term decline in pre-monsoon levels signals extraction beyond recharge. Tracers reveal the origin, age and paths of groundwater. Natural tracers include tritium and carbon-14 for age, the stable isotopes of oxygen and hydrogen (δ¹⁸O and δ²H) for the source and altitude of recharge, chlorofluorocarbons and chloride; artificial tracers include fluorescent dyes such as fluorescein and rhodamine and salts such as bromide. When heads fall in confined aquifers interbedded with clay, the clays compact and the land subsides, often irreversibly; Mexico City and parts of the Central Valley of California have sunk by several metres.

India is the world's largest user of groundwater, which supplies most of its irrigation and rural drinking water. Satellite gravity measurements from the GRACE mission showed rapid depletion across north-western India, in Punjab, Haryana, Rajasthan and Delhi. The effects of over-use are falling water tables, dry shallow wells, rising pumping costs, subsidence, the intrusion of saline water in coastal aquifers and, in places, the mobilisation of fluoride and arsenic. The Central Ground Water Board assesses each assessment unit by its stage of groundwater extraction, annual extraction as a percentage of the annual extractable resource: up to 70% is safe, 70 to 90% semi-critical, 90 to 100% critical and above 100% over-exploited. Groundwater quality problems are both geogenic, fluoride in much of the peninsula and Rajasthan, arsenic in the Ganga-Brahmaputra plains, iron in the north-east and salinity in arid and coastal areas, and anthropogenic, nitrate from fertilisers and sewage, pathogens from on-site sanitation, and leachate from landfills and industry.

In a coastal or island aquifer, fresh groundwater floats on denser sea water. The Ghyben-Herzberg relation, derived independently by Badon Ghyben (1888) and Herzberg (1901), balances the two columns hydrostatically: z = [ρf ÷ (ρs − ρf)] h, where h is the height of the water table above sea level, z the depth of the fresh-salt interface below sea level, and ρf and ρs the densities of fresh (1.000 g/cm³) and sea water (1.025 g/cm³). Then z = 1.000 ÷ 0.025 × h = 40h. If the water table stands 1.5 m above sea level, fresh water extends about 60 m below sea level. The corollary is the danger: lowering the water table by 1 m by pumping raises the interface by about 40 m, drawing salt water into wells, which is why coastal aquifers must be pumped with great care.

4. Natural resource exploration and exploitation, and the conservation of non-renewable resources

Mineral and fuel resources are found by geological mapping, geochemical sampling of soils, stream sediments and plants, geophysical surveys (gravity, magnetic, electrical and seismic) and remote sensing, and then proved by drilling. Their exploitation carries environmental costs at every stage. Open-cast mining, the commonest method for coal, iron ore and bauxite in India, removes forest and topsoil, creates overburden dumps and voids, and generates dust; underground mining can cause subsidence and fires, as in the Jharia coalfield. Acid mine drainage forms when sulphide minerals such as pyrite are exposed to air and water: pyrite oxidises to sulphuric acid and dissolved iron, with bacteria such as Acidithiobacillus ferrooxidans accelerating the reaction, and the acid leaches heavy metals into streams. Mining also lowers water tables, silts rivers and displaces communities; oil and gas extraction adds spills, gas flaring and produced water.

Historically, resource use grew with each industrial revolution, from charcoal and wind to coal, then oil and gas, and the worry of exhaustion is old: M. King Hubbert predicted in 1956 that United States oil production would peak around 1970, which it did for conventional oil, and Hotelling (1931) showed that the scarcity rent of an exhaustible resource should rise at the rate of interest. Conservation of non-renewable resources relies on using less (efficiency and demand management), using longer (durable design and repair), using again (reuse, recycling and recovery of metals from scrap and waste), substituting renewable or more abundant materials, and planning mine closure and reclamation so that land and water are restored.

5. Catastrophic geological hazards: prediction and mitigation

HazardCauses and characterPrediction and mitigation
FloodsRiverine floods from prolonged monsoon rain, flash floods from intense storms, coastal floods from storm surges; aggravated by deforestation, encroachment on floodplains and siltingForecasting by the Central Water Commission, floodplain zoning, embankments, reservoirs, wetland conservation and flood-proofing
LandslidesSteep slopes, weak or fractured rock, saturation by heavy rain, earthquakes, and road cutting or deforestation that removes supportLandslide hazard zonation (the Geological Survey of India is the nodal agency), drainage, retaining walls, rock bolting and afforestation
EarthquakesSudden slip on faults, mostly at plate boundaries; the Himalaya is highly active because of the India-Eurasia collisionCannot be predicted in time and place; hazard mapping (India's seismic zones II to V, V the highest), earthquake-resistant building codes and early warning
VolcanismLava flows, pyroclastic flows, ash falls, mudflows (lahars) and gases; large eruptions inject sulphate aerosols that cool the climate, as Pinatubo did by about 0.5 °C after 1991Monitoring of seismicity, ground swelling and gas emissions allows warning and evacuation; India's only active volcano is on Barren Island in the Andaman Sea
AvalanchesRapid slides of snow on steep slopes after heavy snowfall, warming or disturbanceForecasting for Himalayan passes and roads, snow fences, deflecting structures and afforestation of starting zones
TsunamisLong-wavelength sea waves generated by undersea earthquakes, landslides or eruptions; the Indian Ocean tsunami of 26 December 2004 followed a magnitude 9.1 earthquake off SumatraThe Indian Tsunami Early Warning Centre at INCOIS, Hyderabad (from 2007), coastal zoning, mangrove belts and evacuation drills
CloudburstsExtremely intense, localised rainfall, of the order of 100 mm or more in an hour over a small area, typically where moist monsoon air is forced up steep mountain slopesVery hard to forecast; mitigation through hazard zonation, keeping settlements off stream beds and alluvial fans, and Doppler radar nowcasting

Earthquake size is expressed as magnitude, a measure of the energy released, and intensity, a measure of the effects at a place. Richter's local magnitude (1935) is logarithmic: each whole unit is a tenfold increase in the amplitude of ground motion and about a 32-fold increase in energy, so a magnitude 7 earthquake has 100 times the amplitude and about 1,000 times the energy of a magnitude 5; seismologists now use the moment magnitude for large events. Intensity is graded on the Modified Mercalli scale from I to XII. Mitigation in India is organised under the Disaster Management Act 2005, which created the National Disaster Management Authority, and follows the global Sendai Framework for Disaster Risk Reduction 2015–2030, which shifts the emphasis from responding to disasters to reducing risk before they happen.

Key takeaways

  • Aquifers are unconfined, confined or perched; specific yield plus specific retention equals porosity; India's alluvial aquifers are productive while the peninsular hard rocks store water only in weathered and fractured zones.
  • Darcy's law, Q = K i A, holds for laminar flow; the seepage velocity K i ÷ nₑ, not the Darcy flux, sets how fast pollutants travel; K ranges over ten orders of magnitude from gravel to clay.
  • Over-use lowers water tables, causes subsidence and saline intrusion and mobilises fluoride and arsenic; the stage of extraction classes units as safe (up to 70%), semi-critical, critical or over-exploited (above 100%).
  • Ghyben-Herzberg gives z ≈ 40h: fresh water extends about 40 times as far below sea level as the water table stands above it, so a 1 m fall in the water table raises the salt interface by about 40 m.
  • Mining brings land degradation and acid mine drainage; floods, landslides, earthquakes, volcanism, avalanches, tsunamis and cloudbursts are mitigated by zonation, codes and early warning, under the DM Act 2005 and the Sendai Framework.

Practice questions (10)

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.

  1. An aquifer has a hydraulic conductivity of 10 m/day and a hydraulic gradient of 0.002, and the flow crosses a section 50 m wide and 10 m thick. What is the discharge through it in m³/day?

    Numerical answer — type the value.

    Show answer

    Answer: 10

    By Darcy's law Q = K i A = 10 × 0.002 × (50 × 10) = 10 × 0.002 × 500 = 10 m³/day. Forgetting to multiply width by thickness, or using the gradient as 2 instead of 2 ÷ 1,000, are the usual slips.
  2. On a coastal island the water table stands 1.5 m above mean sea level. Taking fresh water as 1.000 g/cm³ and sea water as 1.025 g/cm³, how far below sea level, in metres, does the fresh-salt interface lie according to the Ghyben-Herzberg relation?

    Numerical answer — type the value.

    Show answer

    Answer: 60

    z = [ρf ÷ (ρs − ρf)] h = [1.000 ÷ 0.025] × 1.5 = 40 × 1.5 = 60 m. The total thickness of the fresh-water lens is h + z = 61.5 m; answering 61.5 gives the lens thickness rather than the depth below sea level that was asked.
  3. Darcy's law ceases to be valid for groundwater flow when

    1. the aquifer is made of fine sand
    2. the flow becomes turbulent, as in karst conduits, large fractures and coarse gravel near pumping wells
    3. the aquifer is confined rather than unconfined
    4. the hydraulic gradient is small
    Show answer

    Answer: B — the flow becomes turbulent, as in karst conduits, large fractures and coarse gravel near pumping wells

    Darcy's law assumes laminar flow dominated by viscous forces, so it fails when the Reynolds number rises and flow becomes turbulent, in wide openings at high velocity. It applies well to sands, small gradients and both confined and unconfined aquifers.
  4. Which statements about a confined aquifer are correct? Select all that apply.

    1. It is bounded above and below by low-permeability layers
    2. The water level in a well tapping it rises above the top of the aquifer
    3. Its upper surface is the water table
    4. It is recharged directly by rainfall over its whole area
    Show answer

    Answer: A — It is bounded above and below by low-permeability layers; B — The water level in a well tapping it rises above the top of the aquifer

    A confined aquifer lies between aquicludes or aquitards and holds water under pressure, so wells rise to the piezometric surface and may flow as artesian wells. It is recharged only where it outcrops, and a water table belongs to an unconfined aquifer.
  5. Under the Central Ground Water Board's assessment, an assessment unit is classed as over-exploited when its stage of groundwater extraction is

    1. between 90% and 100%
    2. between 70% and 90%
    3. above 100% of the annual extractable resource
    4. up to 70%
    Show answer

    Answer: C — above 100% of the annual extractable resource

    Over-exploited means annual extraction exceeds the annual extractable recharge, so storage is being mined. Up to 70% is safe, 70 to 90% semi-critical and 90 to 100% critical.
  6. Assertion (A): Heavy pumping of a coastal aquifer can draw sea water into wells. Reason (R): By the Ghyben-Herzberg relation, each metre that the water table falls lets the fresh-salt interface rise by about 40 m. Choose the correct option.

    1. Both A and R are true, and R is the correct explanation of A
    2. Both A and R are true, but R is not the correct explanation of A
    3. A is true, but R is false
    4. A is false, but R is true
    Show answer

    Answer: A — Both A and R are true, and R is the correct explanation of A

    Both are true and R explains A. Because z ≈ 40h, a small drawdown of the water table produces a forty-fold rise of the interface beneath the well, which is why saline intrusion appears suddenly in over-pumped coastal aquifers.
  7. India's only active volcano is located on

    1. Barren Island in the Andaman Sea
    2. Lakshadweep
    3. Narcondam Island
    4. the Deccan Plateau
    Show answer

    Answer: A — Barren Island in the Andaman Sea

    Barren Island has erupted repeatedly in recent decades. Narcondam is a volcano with no historical eruption, the Deccan Traps are ancient flood basalts about 66 million years old, and Lakshadweep is formed of coral atolls.
  8. Which of the following are natural (environmental) tracers used in groundwater studies? Select all that apply.

    1. Fluorescein dye injected into a well
    2. Chlorofluorocarbons dissolved from the atmosphere
    3. The stable isotope ratio δ¹⁸O
    4. Tritium
    Show answer

    Answer: B — Chlorofluorocarbons dissolved from the atmosphere; C — The stable isotope ratio δ¹⁸O; D — Tritium

    Tritium and CFCs date young groundwater and δ¹⁸O identifies the source and altitude of recharge, all entering groundwater from the environment. Fluorescein is an artificial tracer deliberately injected to trace flow paths, for example in karst.
  9. Compared with an earthquake of magnitude 5 on the Richter scale, one of magnitude 7 has ground-motion amplitude and released energy greater by factors of roughly

    1. 1,000 and 100
    2. 2 and 4
    3. 100 and 1,000
    4. 20 and 64
    Show answer

    Answer: C — 100 and 1,000

    Each unit of magnitude is a tenfold increase in amplitude and about a 32-fold increase in energy, so two units give 10² = 100 times the amplitude and about 32² ≈ 1,000 times the energy. The magnitude scale is logarithmic, not linear, which rules out 2 and 4.
  10. Which earth material typically has the lowest hydraulic conductivity?

    1. Clean gravel
    2. Karstic limestone
    3. Unfractured clay
    4. Clean sand
    Show answer

    Answer: C — Unfractured clay

    Clay's pores are tiny and poorly connected, giving hydraulic conductivity below about 10⁻⁹ m/s despite its high porosity, which is why clay is used to line landfills. Gravel, sand and karst limestone transmit water many orders of magnitude faster.