Water Resources Engineering Interview Questions

The hydrological cycle describes continuous movement of water between Earth's surface and atmosphere. Main components: precipitation (rain, snow) reaches Earth's surface, some water becomes runoff flowing to streams and oceans, some infiltrates into soil becoming groundwater, plants absorb water releasing it through transpiration, and water evaporates from surfaces. Combined evaporation and transpiration is called evapotranspiration. This cycle maintains water balance and drives all water resources engineering.

A catchment area (watershed or drainage basin) is the land area that drains water to a common outlet point like a river, lake, or ocean. It is delineated by topographic divides or ridgelines where water flows in opposite directions. Catchment area directly affects runoff volume - larger catchments produce more runoff for the same rainfall. Understanding catchment boundaries is essential for water resource planning, flood estimation, and environmental assessment.

Dam types include: Gravity dams - concrete mass resisting water by weight, triangular cross-section. Arch dams - curved concrete, transfer load to abutments, suitable for narrow valleys. Buttress dams - concrete face supported by buttresses, economical for wide valleys. Embankment dams - earth or rockfill with impervious core or facing, most common worldwide. Selection depends on valley topography, foundation conditions, material availability, and purpose (storage, diversion, flood control).

Runoff is the portion of precipitation that flows over the land surface to streams and rivers instead of infiltrating into soil. Factors affecting runoff include: rainfall intensity and duration, soil type and moisture content, land slope, vegetation cover, land use (urbanization increases runoff), and catchment shape and size. Runoff coefficient is the ratio of runoff to rainfall, ranging from 0.1 for forests to 0.9 for paved areas.

Open channel flow is classified by time and space variation: Steady vs Unsteady (depth changes with time), Uniform vs Non-uniform (depth changes along channel). Combinations: steady uniform (constant depth everywhere, like long canal), steady non-uniform (gradually or rapidly varied flow, like spillway), and unsteady non-uniform (flood waves). Also classified by Froude number: subcritical (Fr<1, tranquil), critical (Fr=1), and supercritical (Fr>1, rapid).

Irrigation methods include: Surface irrigation - water flows over soil (flood, furrow, basin) - simplest but least efficient (40-60%). Sprinkler irrigation - water sprayed through nozzles simulating rain (60-75% efficiency). Drip/trickle irrigation - water delivered directly to plant roots through emitters (85-95% efficiency). Sub-surface irrigation - water applied below ground surface. Selection depends on crop type, water availability, topography, soil type, and economics.

Manning's equation calculates velocity in open channels: V = (1/n) * R^(2/3) * S^(1/2), where V = velocity (m/s), n = Manning's roughness coefficient, R = hydraulic radius (area/wetted perimeter), S = channel bed slope. Discharge Q = V * A. Manning's n depends on channel surface: 0.010 for smooth concrete, 0.025 for natural streams, 0.05 for heavily vegetated. This empirical equation is fundamental for channel design and flood analysis.

Unconfined (water table) aquifer has a free water surface open to atmosphere, bounded below by impermeable layer. Water level rises and falls with recharge. Confined aquifer is bounded above and below by impermeable layers, water under pressure. When tapped, water rises above the aquifer top; if pressure sufficient, artesian flow occurs. Unconfined aquifers are more susceptible to contamination but easier to recharge; confined aquifers provide more consistent water quality but recover slowly.

A spillway is a structure that safely passes flood water over, around, or through a dam to prevent overtopping, which could cause catastrophic failure of embankment dams. Types include: overflow (ogee) spillway built into dam, side channel spillway, chute spillway for earth dams, shaft/morning glory spillway for limited space. Spillway capacity must exceed maximum expected flood. Spillways include energy dissipators (stilling basins, flip buckets) to prevent erosion downstream.

Return period (recurrence interval) is the average time between occurrences of a hydrological event of specified magnitude. A 100-year flood has 1% probability of occurring in any given year, not that it occurs exactly every 100 years. Return period T = 1/P, where P is annual exceedance probability. Design return periods vary: 2-10 years for minor drainage, 25-50 years for highways, 100 years for major structures, and up to 10,000 years for large dams. Longer return periods mean rarer, more severe events.

Duty of water is the area of land that can be irrigated with a unit discharge of water. Expressed in hectares per cumec (ha/m3/s). Higher duty means water is used more efficiently. Factors affecting duty: climate (evaporation losses), soil type (percolation), crop type (water requirement), conveyance efficiency, and irrigation method. Delta of water is total depth of water required by crop during growing season. Relationship: Duty = 8.64 * crop period (days) / Delta (cm).

Specific energy (E) is the total energy per unit weight of water measured from channel bed: E = y + V^2/(2g) = y + Q^2/(2gA^2), where y = flow depth, V = velocity. For given discharge, specific energy varies with depth, reaching minimum at critical depth. Below critical depth, flow is supercritical (high velocity, shallow); above, flow is subcritical (low velocity, deep). Specific energy concepts are fundamental for analyzing channel transitions, contractions, and hydraulic jumps.

Rain gauge types include: Non-recording (standard) gauge - simple cylinder measuring accumulated rainfall, read manually daily. Recording gauges provide continuous record: weighing type (weight of collected water recorded), tipping bucket (tips after fixed amount, counts pulses), siphon/float type (float rises with water level). Modern automatic weather stations use electronic gauges with telemetry. Radar and satellite provide spatial precipitation estimates. Gauge placement away from obstructions is critical for accurate measurement.

A hydraulic jump is a rapid transition from supercritical (high velocity, shallow) to subcritical (low velocity, deep) flow, accompanied by significant energy loss and turbulence. It occurs when downstream conditions force subcritical flow but upstream conditions create supercritical flow. Common locations: at base of spillways, below sluice gates, in stilling basins. The jump dissipates excess kinetic energy preventing downstream erosion. Jump height ratio depends on upstream Froude number.

Reservoir storage zones from bottom up: Dead storage - below lowest outlet, stores sediment, cannot be used. Active/live storage - between dead storage level and normal pool level, for water supply, irrigation, or power generation. Flood control storage - between normal pool and maximum pool levels, reserved for temporary flood detention. Total storage includes all zones. Reservoir operation rules balance these uses. Sedimentation gradually fills dead storage, eventually affecting active storage if not managed.

The Rational Method estimates peak discharge for small catchments (<200 ha): Q = CIA/360, where Q = peak discharge (m3/s), C = runoff coefficient (0.1-0.95), I = rainfall intensity (mm/hr) for duration equal to time of concentration, A = catchment area (ha). Time of concentration is time for water to travel from most distant point to outlet. Method assumes: uniform rainfall over area, peak occurs when entire catchment contributes, and constant C during storm. Limited to small catchments where assumptions are reasonable.

Gradually varied flow (GVF) has gradual depth changes along channel with constant discharge. The GVF equation: dy/dx = (S0 - Sf)/(1 - Fr^2), where S0 = bed slope, Sf = friction slope, Fr = Froude number. Flow profile types depend on bed slope (mild, steep, critical, adverse, horizontal) and depth relative to normal and critical depths. Computation uses standard step method (known sections) or direct step method (known depths). Profiles predict water surface through channel transitions, over weirs, and approaching controls.

Darcy's law describes groundwater flow through porous media: Q = -KA(dh/dl), where Q = discharge, K = hydraulic conductivity, A = cross-sectional area, dh/dl = hydraulic gradient. Negative sign indicates flow from high to low head. Darcy velocity v = Q/A is apparent velocity; actual pore velocity = v/porosity. Law valid for laminar flow (Reynolds number < 1 based on grain diameter). Applications include well yield estimation, seepage under dams, and contaminant transport calculations.

A unit hydrograph (UH) shows the direct runoff response of a catchment to one unit (usually 1 cm) of excess rainfall occurring uniformly over a specified duration. Assumptions: linear system (proportional response), time invariance, excess rainfall uniform in space and time. To get runoff hydrograph: apply UH to each rainfall increment and sum (superposition). Synthetic UHs (SCS, Snyder) can be developed from catchment characteristics when data unavailable. UH is fundamental tool for flood hydrograph prediction.

Weir discharge formulas: Rectangular weir: Q = Cd * (2/3) * sqrt(2g) * L * H^1.5, where Cd = discharge coefficient (~0.62), L = crest length, H = head above crest. V-notch (triangular): Q = Cd * (8/15) * sqrt(2g) * tan(theta/2) * H^2.5, commonly 90-degree notch (theta = 90). Broad-crested weir: Q = Cd * L * sqrt(g) * H^1.5. Actual discharge depends on approach velocity, nappe aeration, and end contractions. Weirs serve as flow measurement devices and control structures.

Crop water requirement (ETc) = reference evapotranspiration (ETo) * crop coefficient (Kc). ETo calculated from weather data using Penman-Monteith equation (FAO-56 method) based on temperature, humidity, wind, and solar radiation. Kc varies through crop growth stages: initial, development, mid-season, late season. Net irrigation requirement = ETc - effective rainfall - groundwater contribution. Gross requirement accounts for application efficiency and conveyance losses. Irrigation scheduling determines when and how much to irrigate.

Gravity dam stability checks: (1) Overturning - factor of safety (resisting moment/overturning moment) > 1.5 for normal, > 1.2 for extreme loads, (2) Sliding - friction resistance plus shear key resistance > driving forces, FOS > 1.5, (3) Stress - maximum compression < allowable (0.3f'c), no tension at heel (or limited tension with drain effectiveness), (4) Base pressure within foundation capacity. Loading conditions include: normal (dead load, water, uplift), unusual (ice, silt), and extreme (earthquake, maximum flood). Uplift reduced by drainage systems.

Theis equation for drawdown in confined aquifer under unsteady pumping: s = (Q/4piT) * W(u), where s = drawdown, Q = pumping rate, T = transmissivity, W(u) = well function, u = r^2S/(4Tt), r = distance from well, S = storativity, t = time since pumping started. Well function values from tables or software. Equation describes cone of depression expanding over time. Assumptions: homogeneous, isotropic, infinite aquifer, fully penetrating well, constant pumping. Used for aquifer testing and impact assessment.

Flood routing tracks flood wave movement through systems. Reservoir routing: Modified Puls (storage indication) method uses continuity and storage-outflow relationship iteratively. Channel routing: Muskingum method uses storage = K[xI + (1-x)O], where K = storage constant, x = weighting factor (0-0.5). Muskingum-Cunge provides physical basis for parameters. Kinematic wave routing for overland flow. Hydraulic (dynamic) routing solves full Saint-Venant equations for detailed analysis. Choice depends on required accuracy and data availability.

Irrigation canal design balances: (1) Non-silting velocity - sufficient to transport sediment (~0.5-0.8 m/s), (2) Non-scouring velocity - prevent erosion (<1.0 m/s for earth canals), (3) Minimum section for required discharge, (4) Stable side slopes (1.5:1 to 2:1 for earth). Kennedy's theory relates velocity to sediment transport; Lacey's regime theory provides complete channel dimensions. Lined canals allow higher velocities (2-3 m/s), reduce seepage losses. Design includes freeboard, berm, and bank dimensions.

Flood frequency analysis estimates flood magnitudes for various return periods: (1) Collect annual maximum series from stream gauge records, (2) Rank data and assign plotting positions (Weibull, Gringorten), (3) Fit probability distribution (Log-Pearson Type III, Gumbel, Log-Normal) using method of moments or maximum likelihood, (4) Use fitted distribution to estimate flood quantiles for design return periods. Regional analysis combines nearby stations to improve estimates. Results subject to uncertainty from limited record length and non-stationarity (climate change).

Seepage control in earth dams through: (1) Core - central impervious zone (clay) to minimize flow through dam, (2) Cutoff - extending impervious layer into foundation (trench, sheet pile, grout curtain), (3) Upstream blanket - extends impervious zone to increase seepage path, (4) Horizontal drain/filter - collects seepage water safely without piping, (5) Relief wells - reduce uplift pressure in foundation. Filters designed with gradation criteria preventing particle migration (D15 filter/D85 base < 4-5). Exit gradients must be below critical to prevent piping failure.

Detention basin design: (1) Determine inflow hydrograph for design storm using rainfall-runoff methods, (2) Set allowable peak outflow (often pre-development rate), (3) Route inflow through basin using storage-indication method to find required storage, (4) Size basin for required volume with adequate freeboard, (5) Design outlet structure (orifice, weir combinations) to achieve target outflow rate, (6) Include emergency spillway for larger storms. Multi-stage outlets can control multiple frequencies. Extended detention for water quality provides settling time before release.

Energy dissipation methods: (1) Stilling basin - hydraulic jump confined in basin with floor blocks and end sill to force jump (USBR Type I-IV designs), (2) Flip bucket/ski jump - throws jet away from dam, energy dissipated in plunge pool, (3) Roller bucket - submerged bucket creates reverse roller, (4) Stepped spillway - energy dissipated over multiple steps, (5) Impact-type - flow impacts floor or walls. Selection depends on tailwater conditions, space availability, downstream erodibility, and discharge magnitude. Proper design prevents erosion that could undermine dam.

Conjunctive use is coordinated management of surface water and groundwater to maximize water availability and reliability. Surface water (variable but large) used during wet periods; groundwater (stable but limited yield) supplements during dry periods. Methods include: aquifer storage and recovery (ASR), recharge basins, streambed infiltration, and managed aquifer recharge. Benefits include drought resilience, peak demand management, reduced infrastructure costs, and water quality blending. Requires understanding of aquifer properties and surface-groundwater interactions.

Depth-Area-Duration (DAD) curves show how average rainfall depth decreases as averaging area increases for storms of given duration. Point rainfall is always maximum; areal average is less due to storm spatial variability. DAD analysis uses isohyetal maps of recorded storms. Index values express areal depth as percentage of point depth. Used for: design storms for large catchments, PMP development, spillway design flood estimation. Longer duration storms have smaller DAD reduction because spatial variation averages out over time.

Critical depth occurs when specific energy is minimum for given discharge, or Fr = 1. For rectangular section: yc = (Q^2/(gB^2))^(1/3) = (q^2/g)^(1/3), where q = unit discharge. For triangular: yc = (2Q^2/(gm^2))^(1/5), m = side slope. For general sections, solve: Q^2*T/(g*A^3) = 1, where T = top width. Critical slope maintains critical depth as normal depth. Critical depth is control point for flow calculations, occurring at free overfall, channel transitions, and minimum energy state.

Drip irrigation design considers: (1) Emitter selection - flow rate (2-8 L/hr), pressure compensation for uniform distribution, (2) Lateral layout - spacing based on wetting pattern overlap and crop spacing, (3) Hydraulic design - lateral and manifold sizing to maintain pressure uniformity (emission uniformity > 85%), (4) Filtration - sand, disc, or screen filters to prevent emitter clogging, (5) Fertigation capability - injecting fertilizers through system, (6) Automation - controllers, sensors for scheduling. Operating pressure typically 50-200 kPa. System efficiency 85-95%.

Reservoir sedimentation estimation: (1) Sediment yield from catchment using USLE, sediment rating curves, or regional data, (2) Trap efficiency using Brune or Churchill curves based on capacity-inflow ratio, (3) Sediment distribution using empirical area-reduction method. Management strategies: (1) Watershed treatment to reduce erosion, (2) Sediment exclusion structures, (3) Sluicing through low-level outlets during floods, (4) Dredging (expensive), (5) Density current venting, (6) Flushing. Dead storage accounts for sedimentation; periodic bathymetric surveys monitor accumulation.

Head loss in pipes has friction and minor components. Darcy-Weisbach: hf = f*(L/D)*(V^2/2g), where f = friction factor from Moody diagram (function of Reynolds number and relative roughness). Hazen-Williams: V = 0.849*C*R^0.63*S^0.54, where C = roughness coefficient (150 for plastic, 100 for old iron). Minor losses (valves, bends, fittings): hm = K*(V^2/2g). Energy equation with losses used for pipe system analysis. Hardy-Cross method iteratively solves pipe network problems.

PMP is theoretically maximum precipitation for a given area and duration with no assigned probability. Estimation methods: (1) Meteorological approach - moisture maximization (adjust observed storms to maximum dew point) and transposition (move storms to study location with adjustments), (2) Statistical extrapolation - extending frequency analysis to extreme events (Hershfield), (3) Generalized estimates from hydrometeorological reports (WMO, Bureau of Meteorology). PMP used for spillway design flood (PMF) for high-hazard dams. Climate change may require reassessment of historical PMP values.

Dam break analysis involves: (1) Breach formation - geometry (depth, width, shape) and timing (regression equations or physical modeling), (2) Peak discharge estimation - empirical (Q = f(height, storage)) or physical (Q = Cd*B*(2g*H)^0.5), (3) Downstream routing - solve Saint-Venant equations (1D: HEC-RAS, MIKE 11; 2D: HEC-RAS 2D, MIKE 21) including floodplain storage and bridge effects, (4) Inundation mapping - depths, velocities, arrival times. Results inform emergency action plans (EAP), evacuation zones, and dam classification. Sensitivity analysis for breach parameters essential given uncertainty.

Sediment transport modes: suspended load (fine particles carried by turbulence), bed load (coarse particles rolling/sliding along bed), and wash load (very fine, supply limited). Analysis methods: (1) Bed load - empirical equations (Meyer-Peter-Muller, Einstein, Parker), (2) Suspended load - integration of concentration profile with velocity profile, (3) Total load - Ackers-White, Engelund-Hansen, Yang equations. Modeling requires sediment size distribution, flow hydraulics, and transport capacity. Mobile bed models (HEC-RAS sediment, MIKE 21C) simulate morphological changes over time.

Saint-Venant equations (shallow water equations) describe unsteady open channel flow: Continuity: dA/dt + dQ/dx = ql (lateral inflow), Momentum: dQ/dt + d(Q^2/A)/dx + gA*dh/dx + gA*Sf = 0. Terms: local acceleration, convective acceleration, pressure gradient, friction slope. Solution methods: finite difference (implicit schemes stable for larger timesteps), finite volume, or characteristics. Applications: flood wave propagation, dam break, tidal hydraulics, storm surge. Models: HEC-RAS unsteady, MIKE 11, ISIS. 2D equations add lateral dimension for floodplain modeling.

Distributed models represent spatial variation in catchment properties using grid cells or sub-catchments, unlike lumped models treating catchment as single unit. Types: physically-based (solve Richards equation, energy balance - MIKE SHE, GSSHA), conceptual distributed (HEC-HMS gridded, SWAT for water quality). Each cell has soil, land use, and topography parameters. Input includes spatial rainfall (radar, satellite). Advantages: capture spatial patterns, model ungauged basins, assess land use change impacts. Require extensive data and calibration. Appropriate for large, heterogeneous catchments.

Spillway Design Flood (SDF) determination: (1) Dam hazard classification based on downstream consequences (high, significant, low), (2) Select design standard - high hazard typically uses PMF (Probable Maximum Flood = runoff from PMP), lower hazard may use percentage of PMF or frequency-based flood (1:10,000), (3) Develop design hydrograph using PMP temporal distribution and unit hydrograph or distributed model, (4) Route through reservoir to get peak outflow and maximum water level, (5) Add freeboard for wind/wave setup. Incremental damage analysis may justify reduced standards where no incremental damage from dam failure.

Groundwater modeling process: (1) Conceptualization - define aquifer boundaries, layers, flow system, (2) Grid design - finite difference (MODFLOW) or finite element, refined where detail needed, (3) Boundary conditions - specified head, flux, or mixed at model edges and internal features, (4) Parameter assignment - K, storage, recharge for each cell/element, (5) Calibration - adjust parameters to match observed heads and flows using inverse methods (PEST, UCODE), (6) Sensitivity analysis - identify influential parameters, (7) Validation - test on independent dataset, (8) Prediction scenarios. Uncertainty quantification essential for decision-making.

Flood risk assessment combines hazard and vulnerability: (1) Hazard analysis - flood frequency/magnitude, inundation mapping for multiple scenarios (10, 50, 100, 500-year), depths, velocities, duration, (2) Exposure analysis - inventorying assets in flood zones (buildings, infrastructure, people), (3) Vulnerability assessment - depth-damage functions relating flood depth to damage percentage by building type, (4) Risk quantification - expected annual damage (EAD) = integral of damage-probability curve. Risk mapping guides land use planning, insurance, emergency response. Climate change scenarios assess future risk evolution.

Water hammer is pressure surge from rapid flow velocity change in pipelines. Causes: pump trip, valve closure, demand changes. Joukowsky equation: delta_P = rho*a*delta_V, where a = wave speed (typically 900-1200 m/s). Analysis methods: (1) Arithmetic - hand calculation for simple systems, (2) Method of characteristics - solving hyperbolic PDEs numerically (standard approach), (3) Commercial software (HAMMER, Surge, WANDA). Protection devices: surge tanks, air chambers, pressure relief valves, slow-closing valves, flywheels on pumps. Design must prevent negative pressures (column separation) and overpressure.

Reservoir optimization balances competing objectives (water supply reliability, flood control, hydropower, environment). Methods: (1) Simulation - evaluate operating rules through historical or synthetic inflows, (2) Linear programming - optimize releases subject to constraints for simplified systems, (3) Dynamic programming - find optimal policy considering state (storage) and decisions (release) over time, (4) Stochastic DP - incorporate inflow uncertainty, (5) Evolutionary algorithms - handle complex, nonlinear multi-objective problems. Rule curves define storage zones and release policies. Real-time operation combines forecasts with optimization. Trade-off analysis using Pareto fronts.

Continuous models simulate entire water balance over extended periods, tracking soil moisture states between events. Components: (1) Interception - canopy storage, (2) Infiltration - Green-Ampt, Philip, or conceptual (deficit-constant), (3) Soil moisture accounting - multiple layers, ET extraction, (4) Groundwater - linear reservoirs, (5) Routing - unit hydrograph or kinematic wave. Models: Sacramento (NWS), HBV (Scandinavia), GR4J (France). Calibration requires multi-year records including droughts and floods. Used for: water resource yield, low flow assessment, climate impact studies. More demanding than event-based but more realistic.

Environmental flows maintain river ecosystem health. Assessment methods: (1) Hydrological - percentage of natural flow (Tennant method), flow duration analysis, (2) Hydraulic - maintain habitat (wetted perimeter, depth criteria), (3) Habitat simulation - PHABSIM links flow to habitat suitability indices for target species, (4) Holistic - expert panels considering all ecosystem components (DRIFT, ELOHA). Key flow components: low flows for water quality, high flows for channel maintenance and floodplain connection, seasonal variation for life cycles. E-flows compete with consumptive uses, requiring integrated water management.

Saltwater intrusion occurs when pumping or sea level rise allows denser seawater to migrate into freshwater aquifer. Analysis: Ghyben-Herzberg approximation (interface depth = 40 * head above sea level for static conditions), sharp interface models, or variable-density flow models (SEAWAT, SUTRA) for transient analysis. Management: (1) Pumping reduction or redistribution, (2) Injection barriers creating freshwater ridge, (3) Extraction barriers pumping saltwater, (4) Recharge enhancement, (5) Monitoring wells with electrical conductivity sensors. Climate change increases vulnerability through sea level rise and reduced recharge.

Computational hydraulics solves partial differential equations numerically. Schemes: (1) Finite difference - approximate derivatives on structured grid (explicit: simple but stability-limited, implicit: stable but requires matrix solution), (2) Finite volume - conservative, handles discontinuities well (Godunov-type schemes for shock capturing), (3) Finite element - flexible geometry, higher-order accuracy. Time stepping: Crank-Nicolson (2nd order implicit), Preissmann (4-point implicit for rivers). Stability requires CFL condition for explicit schemes. Parallelization and GPU computing enable large-scale 2D/3D simulations. Model selection balances accuracy, computational cost, and application needs.

IWRM coordinates development and management of water, land, and related resources to maximize economic and social welfare equitably without compromising ecosystem sustainability. Implementation: (1) Enabling environment - policies, legislation, financing, (2) Institutional framework - basin organizations, stakeholder participation, (3) Management instruments - demand management, economic tools (pricing, markets), regulatory mechanisms, (4) Information systems - monitoring, modeling, decision support, (5) Capacity building. Challenges: cross-sectoral coordination, transboundary issues, climate uncertainty, balancing competing demands. UN SDG 6 promotes IWRM. Success requires political will and stakeholder engagement.