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
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.
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
| Hazard | Causes and character | Prediction and mitigation |
|---|---|---|
| Floods | Riverine floods from prolonged monsoon rain, flash floods from intense storms, coastal floods from storm surges; aggravated by deforestation, encroachment on floodplains and silting | Forecasting by the Central Water Commission, floodplain zoning, embankments, reservoirs, wetland conservation and flood-proofing |
| Landslides | Steep slopes, weak or fractured rock, saturation by heavy rain, earthquakes, and road cutting or deforestation that removes support | Landslide hazard zonation (the Geological Survey of India is the nodal agency), drainage, retaining walls, rock bolting and afforestation |
| Earthquakes | Sudden slip on faults, mostly at plate boundaries; the Himalaya is highly active because of the India-Eurasia collision | Cannot 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 |
| Volcanism | Lava 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 1991 | Monitoring of seismicity, ground swelling and gas emissions allows warning and evacuation; India's only active volcano is on Barren Island in the Andaman Sea |
| Avalanches | Rapid slides of snow on steep slopes after heavy snowfall, warming or disturbance | Forecasting for Himalayan passes and roads, snow fences, deflecting structures and afforestation of starting zones |
| Tsunamis | Long-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 Sumatra | The Indian Tsunami Early Warning Centre at INCOIS, Hyderabad (from 2007), coastal zoning, mangrove belts and evacuation drills |
| Cloudbursts | Extremely 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 slopes | Very 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.
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.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.Darcy's law ceases to be valid for groundwater flow when
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.Which statements about a confined aquifer are correct? Select all that apply.
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.Under the Central Ground Water Board's assessment, an assessment unit is classed as over-exploited when its stage of groundwater extraction is
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.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.
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.India's only active volcano is located on
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.Which of the following are natural (environmental) tracers used in groundwater studies? Select all that apply.
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.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
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.Which earth material typically has the lowest hydraulic conductivity?
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.