16-Civ-B4 Engineering Hydrology · May 2018
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Paper format. National Examination, May 2018, 16-Civ-B4 Engineering Hydrology, three hours’ duration, closed book with one two-sided candidate-prepared aid sheet (8½″ × 11″) and one approved Casio or Sharp calculator whose model designation must be written on the first inside left-hand sheet of the work book. Seven problems are printed. Page-1 Note 4 states that any five (5) questions constitute a complete paper and that only the first five answers appearing in the work book will be marked; Note 5 weights each problem at twenty (20) points, so the examinable total is 5 × 20 = 100 points. All seven problems are solved here, because this set is a study resource rather than a timed sitting. Sub-part mark values below are the printed ones from the page-6 marking scheme, which on this sitting is internally consistent — every problem’s sub-parts sum to twenty, and the marks printed in the page margins agree with the scheme.
Reference texts. V. T. Chow, D. R. Maidment and L. W. Mays, Applied Hydrology (hydrologic cycle, unit hydrographs, conceptual models, routing, frequency analysis); L. W. Mays, Water Resources Engineering, 3rd ed. (urban and highway drainage, stormwater management, reservoir operation); W. Viessman and G. L. Lewis, Introduction to Hydrology, 5th ed. (precipitation and streamflow measurement, hydrologic modelling); C. W. Fetter, Applied Hydrogeology, 4th ed. (Darcy’s law and aquifer hydraulics); V. T. Chow, Open-Channel Hydraulics (1959) (flood-wave propagation and unsteady flow). Canadian practice references: Environment and Climate Change Canada Engineering Climate Datasets (short-duration rainfall and IDF curves) and the Water Survey of Canada HYDAT archive; the WMO Manual on Stream Gauging (WMO-No. 1044) and ISO 1100-2 (stage–discharge ratings); the Transportation Association of Canada Guide to Bridge Hydraulics and provincial highway drainage manuals; provincial stormwater management planning and design manuals (e.g. Ontario MOECC 2003, British Columbia Stormwater Planning Guidebook); and the Canadian Dam Association Dam Safety Guidelines (inflow design flood and reservoir routing).
Check — which numbers come from the paper and which are the solver’s. This sitting supplies numerical data in only three places: Problem 2(iii) (the 7 ha suburban development, its 30-minute time of concentration and the IDF relation), Problem 5(iii) (the printed IDF chart, from which an intensity must be read graphically), and Problem 6(iii) (the 50-year return period and the 10-year exposure). Those three answers are computed from the paper’s own data. Every other number below appears inside a short illustrative example whose inputs are stated in an explicit Given line as the solver’s own representative Canadian values; they exist to make a discussion answer concrete and checkable, and they are not exam data. A candidate who assumed different but reasonable values and carried them through consistently would receive the same marks.
Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.
Whatever method is used, the weights must be non-negative and must sum to unity, so that a uniform rainfall field returns its own depth. The two methods described below are the ones in routine Canadian practice; the figure shows both applied to the same three-gauge network.
Method 1 — Thiessen polygons. The gauges are joined by straight lines, the perpendicular bisectors of those lines are constructed, and the polygon formed around each gauge encloses every point in the basin that is nearer to that gauge than to any other. The weight of a gauge is its polygon area divided by the basin area, $\lambda_i = A_i/A$. Its first advantage is objectivity and reproducibility: the construction is purely geometric, two engineers obtain identical weights, and the whole procedure is automated in any GIS as a Voronoi tessellation. Its second advantage is that, because the weights depend only on the gauge positions and the basin boundary, they are computed once and reused for every storm and every month of record, which makes long-record processing cheap; and they properly handle a network of uneven spacing, including gauges outside the basin whose polygons intrude into it, which a simple arithmetic mean cannot do.
Method 2 — the isohyetal method. The observed depths are plotted at the gauge locations, lines of equal depth (isohyets) are drawn between them, and the mean depth is obtained as the area-weighted average of the mean depths of the bands between successive isohyets. Its first advantage is that the analyst may impose knowledge that the data alone do not contain: isohyets are drawn with reference to topography, to the known track and size of the storm cell, and to radar imagery, so orographic enhancement on a windward slope or a convective cell between two gauges can be represented. This makes it the most accurate method in mountainous terrain and for localised storms, which is the usual British Columbian case. Its second advantage is that it represents the rainfall as a continuous field with gradients rather than as uniform blocks, so it uses all the information in the record, including gauges outside the basin, and the resulting map is itself a useful product for checking the plausibility of the storm and for distributing rainfall to sub-basins.
Other admissible answers include the arithmetic mean ($\lambda_i = 1/N$), adequate only for a dense, uniform network in flat terrain; inverse-distance weighting, where $\lambda_i \propto 1/d_i^{\,2}$ normalised to unit sum; and kriging, in which the weights are chosen to minimise the estimation variance given a fitted semivariogram, which is the only method that also returns a map of the uncertainty of the estimate.
Given. Three gauges on an 80 km² basin recording 62, 48 and 85 mm in one storm, with Thiessen polygon areas of 40, 25 and 15 km² respectively. Find. The Thiessen mean areal depth, and its comparison with the arithmetic mean.
Difference 1 — what is actually measured, and what the hydrologic model needs. A point method reports the depth caught in a single orifice roughly 200 mm across, and treats that as representative of the surrounding tens or hundreds of square kilometres; the measurement is direct, continuous in time, and precise at the point. An areal method produces a single depth representing the average over the whole basin, either by combining several point records with weights or by integrating a spatially continuous field from weather radar or satellite. Only the areal quantity is dimensionally what a rainfall–runoff model needs, because runoff volume is a depth times an area; the point value enters only as one of its ingredients.
Difference 2 — the error structure, and the systematic smoothing of extremes. A point measurement carries instrument and exposure errors — wind under-catch, which can reach 5 to 15 per cent for rain and far more for snow, wetting and evaporation losses, splash, and gauge siting effects — and these are local, well understood and correctable gauge by gauge. An areal estimate inherits all of those and adds a sampling error that depends on the density and layout of the network relative to the size of the storm: a convective cell a few kilometres across may fall entirely between gauges and be missed, or fall on one gauge and be extrapolated over the whole basin. The consequence is a systematic one that must be allowed for in design: the areal mean depth is always less than the maximum point depth for the same storm, and the ratio of the two — the areal reduction factor — decreases as the area grows and as the duration shortens. Applying a point IDF value uniformly over a large watershed therefore over-estimates the design rainfall, and it is only for very small urban catchments that the point value may be used directly.
A third difference worth a line is temporal: a point gauge gives a continuous high-resolution intensity trace, which is what frequency analysis and IDF construction require, whereas areal methods are normally applied storm by storm or month by month and inherit whatever temporal resolution the weighting scheme can support.
Definition. Stream flow discharge is the volume of water passing a specified cross-section of a channel per unit of time, with dimensions $\mathrm{L}^3\mathrm{T}^{-1}$ and reported in Canadian practice in cubic metres per second. It is the integral of the velocity normal to the section over the flow area, $Q=\int_A v\,dA$, so that every method of measuring it is in essence a way of evaluating that integral, directly or by proxy.
Method 1 — the velocity–area method (current meter or acoustic Doppler profiler). The section is divided into twenty to twenty-five vertical panels, each carrying no more than about five per cent of the total flow. In each vertical the depth is sounded and the mean velocity is measured, either at six-tenths of the depth from the surface or as the average of readings at two-tenths and eight-tenths. The mid-section computation then sums the panel discharges, $Q=\sum b_i d_i \bar{v}_i$. It is a direct measurement, accurate to about five per cent for a careful gauging, and applicable over a very wide range of channel sizes from a wading measurement to a cableway or a boat-mounted acoustic profiler. Against that, it is slow (half an hour to several hours), labour-intensive, hazardous or impossible in a large flood, disturbed by ice cover, and it yields only an instantaneous discharge — one point on a hydrograph.
Method 2 — the stage–discharge rating curve. Because a continuous discharge record is what is really wanted, the standard arrangement is to record water level continuously with a float in a stilling well or a pressure transducer, and to convert it with a rating curve of the form $Q = a(h - h_0)^{\,b}$, where $h_0$ is the stage of zero flow at the control. The rating is established from many velocity–area gaugings spread over the range of stage. Its advantages are decisive for a network: it produces a continuous, telemeterable, low-cost record from a single cheap measurement, and it makes real time forecasting possible. Its limitations are equally clear: it depends on a stable control, so the rating must be re-verified after every significant flood or channel change and shifted if the control has moved; the highest flows are almost always extrapolated beyond the largest gauging, which is exactly where the design values come from; and during a steep flood wave the rating is not single-valued at all, since the water-surface slope is greater on the rising limb than on the falling limb, producing the familiar looped rating.
Given. A three-panel wading gauging with panel widths of 2.0 m each, depths of 1.20, 1.80 and 1.40 m, and mean panel velocities of 0.65, 0.95 and 0.72 m/s; at the same site a rating $Q = 12.5\,(h-0.30)^{2.15}$ with $h$ and the datum in metres. Find. The gauged discharge, and the discharge the rating gives at a stage of 2.10 m.
Other methods that would be credited are dilution (tracer) gauging, which suits steep turbulent mountain streams where a current meter cannot be used; a calibrated weir or flume for a small research catchment; and the indirect slope–area method using Manning’s equation and surveyed high-water marks, which is often the only way to estimate a peak that has already passed.
