Skip to main content

Thiessen weights for semicircular catchment

Short answer

Perpendicular bisectors of stations at (0,0), (r/2,r/2), (-r/2,r/2) divide the semicircle. O's polygon is a triangle with area r²/4; total area πr²/2 gives weight 0.5/π. The rest splits symmetrically, giving each remaining gauge weight 0.5 - 0.25/π.

The catchment area of a basin may be approximated as a semicircle of radius r km with respect to the coordinate axis set up with its origin at the center of the circle and the x-axis coincident with the diameter; the area lies in the first and second quadrants and the position coordinates of the rain gauge stations are (0,0)(0,0), (r2,r2)\left(\frac{r}{2},\frac{r}{2}\right) and (−r2,r2)\left(\frac{-r}{2},\frac{r}{2}\right) km. Show that the Thiessen weights of the gauges are given by 0.5π\frac{0.5}{\pi}, (0.5−0.25/π)(0.5 - 0.25/\pi) and (0.5−0.25/π)(0.5 - 0.25/\pi) respectively.

Setup: Semicircle of radius rr in the upper half-plane, area A=πr2/2A = \pi r^2/2. Gauges: O(0,0)O(0,0), P1(r/2,r/2)P_1(r/2, r/2), P2(−r/2,r/2)P_2(-r/2, r/2).

        y
        |   P2 .  | . P1
        |   (-r/2,r/2) (r/2,r/2)
   -----+----O-----+---- x
       -r          r

Step 1: Perpendicular bisectors

  • Between P1P_1 and P2P_2: they are at the same height, so the bisector is x=0x = 0.
  • Between OO and P1P_1: points equidistant from both satisfy x2+y2=(x−r/2)2+(y−r/2)2x^2+y^2 = (x-r/2)^2 + (y-r/2)^2, which gives x+y=r/2x + y = r/2.
  • Between OO and P2P_2: similarly −x+y=r/2-x + y = r/2.

Step 2: Polygon of the gauge at O

The area closer to OO than to P1P_1 or P2P_2 is bounded by the x-axis, y=r/2−xy = r/2 - x and y=r/2+xy = r/2 + x. This is a triangle with vertices (−r/2,0)(-r/2, 0), (r/2,0)(r/2, 0) and (0,r/2)(0, r/2) (all inside the semicircle):

AO=12×r×r2=r24A_O = \tfrac{1}{2} \times r \times \frac{r}{2} = \frac{r^2}{4} WO=AOA=r2/4πr2/2=0.5πW_O = \frac{A_O}{A} = \frac{r^2/4}{\pi r^2/2} = \frac{0.5}{\pi}

Step 3: Remaining two gauges

The rest of the semicircle is divided by x=0x = 0 into two equal halves (by symmetry):

AP1=AP2=12(πr22−r24)A_{P_1} = A_{P_2} = \frac{1}{2}\left(\frac{\pi r^2}{2} - \frac{r^2}{4}\right) WP1=WP2=12(1−0.5π)=0.5−0.25πW_{P_1} = W_{P_2} = \frac{1}{2}\left(1 - \frac{0.5}{\pi}\right) = 0.5 - \frac{0.25}{\pi}

Numerically: WO=0.159W_O = 0.159, WP1=WP2=0.4204W_{P_1} = W_{P_2} = 0.4204; the sum =1.000= 1.000. Hence proved.

✓ Answer checked by Claude Sonnet 5.5

All Engineering Hydrology chapter 2 questions →