4.2 Leveling, Earthworks, and Mass Haul Diagrams

Key Takeaways

  • Differential leveling establishes vertical control through height of instrument (HI = Elev + BS) and foresight reductions (Elev = HI - FS), with the mandatory mathematical check: ∑BS−∑FS=Last Elev−First Elev\sum \text{BS} - \sum \text{FS} = \text{Last Elev} - \text{First Elev}.

  • The combined effect of earth curvature and atmospheric refraction displaces optical sightlines upward, requiring a subtractive rod correction of hcr=0.0675K2h_{cr} = 0.0675 K^2 (with hcrh_{cr} in meters and sight distance KK in kilometers).

  • Three-level cross-sectional areas are evaluated via A=12[c(wL+wR)+b2(hL+hR)]A = \frac{1}{2}[c(w_L + w_R) + \frac{b}{2}(h_L + h_R)], and prismoidal volumes between successive stations are computed via Vp=Ve−CpV_p = V_e - C_p, where Cp=L12(c1−c2)(w1−w2)C_p = \frac{L}{12}(c_1 - c_2)(w_1 - w_2).

  • Borrow pit excavations are computed using the unit base method: V=Aunit4(∑h1+2∑h2+3∑h3+4∑h4)V = \frac{A_{\text{unit}}}{4}(\sum h_1 + 2\sum h_2 + 3\sum h_3 + 4\sum h_4), where corner heights are weighted by the number of adjoining grid squares.

  • On a mass haul diagram, ascending sections indicate excavation (cut) and descending sections indicate embankment (fill); the limit of economic haul is governed by LEH=FHD+CborrowCoverhaul\text{LEH} = \text{FHD} + \frac{C_{\text{borrow}}}{C_{\text{overhaul}}}.

Last updated: October 2026

4.2 Leveling, Earthworks, and Mass Haul Diagrams

Vertical control and volumetric computations constitute the core of highway geometric design, railway grading, site development, and drainage infrastructure. This section reviews differential leveling protocols, curvature and refraction corrections, three-level cross-section area formulas, prismoidal volume adjustments, borrow pit grading, and the economic optimization of mass haul diagrams.


1. Principles of Differential Leveling

Differential leveling determines elevation differences between points by measuring vertical distances on graduated leveling rods relative to a horizontal line of sight established by an automatic level, digital level, or tilting level.

Essential Definitions

  • Bench Mark (BM): A permanent or semi-permanent physical marker of precisely known elevation above a reference datum (e.g., Mean Sea Level).
  • Turning Point (TP): An intermediate, stable ground point used to advance the level line forward, serving successively as a foresight and a backsight.
  • Backsight (BS or Plus Sight): A rod reading taken on a point of known elevation to determine the Height of Instrument (HI).
  • Height of Instrument (HI): The elevation of the horizontal line of sight passing through the optical center of the telescope: HI=Elevation+BS\text{HI} = \text{Elevation} + \text{BS}
  • Foresight (FS or Minus Sight): A rod reading taken on a turning point or benchmark of unknown elevation to transfer the vertical datum: Elevationnew=HI−FS\text{Elevation}_{\text{new}} = \text{HI} - \text{FS}
  • Intermediate Foresight (IFS): Rod readings taken along cross-sections, building corners, or topography from a single instrument setup where no subsequent backsight is made.

Leveling Field Book Structure & Arithmetic Check

In professional surveying and board exam evaluations, leveling notes must satisfy a strict mathematical identity. The algebraic sum of all backsights minus the algebraic sum of all foresights must equal the final elevation minus the initial elevation: ∑BS−∑FS=ElevationFinal−ElevationInitial\sum \text{BS} - \sum \text{FS} = \text{Elevation}_{\text{Final}} - \text{Elevation}_{\text{Initial}}

StationBS (m)HI (m)FS (m)Elevation (m)Remarks
BM-11.842101.842—100.000Known datum point
TP-12.115102.7211.236100.606Turning point 1
TP-20.954101.8101.865100.856Turning point 2
BM-2——2.41899.392Closing benchmark
SUM4.911—5.519—Check: 4.911−5.519=−0.608 m4.911 - 5.519 = -0.608\text{ m}

Check Verification: 99.392−100.000=−0.608 m99.392 - 100.000 = -0.608\text{ m}. The notes are mathematically closed.


2. Earth Curvature and Atmospheric Refraction

Over sight distances exceeding 100 m100\text{ m}, the earth's curvature and atmospheric refraction introduce systematic optical errors:

  1. Earth Curvature (hch_c): The earth's surface curves away from the tangent line of sight. The rod reading is displaced upward by: hc=K22R≈0.0785K2(m, with K in km)h_c = \frac{K^2}{2R} \approx 0.0785 K^2 \quad (\text{m, with } K \text{ in km})
  2. Atmospheric Refraction (hrh_r): Density gradients in the atmosphere bend optical rays downward toward denser air, causing the sightline to hit the rod lower than the geometric tangent line by approximately one-seventh of curvature: hr≈0.0110K2(m, with K in km)h_r \approx 0.0110 K^2 \quad (\text{m, with } K \text{ in km})
  3. Combined Correction (hcrh_{cr}): Combining both effects yields the standard net correction: hcr=hc−hr=0.0785K2−0.0110K2=0.0675K2(or 0.067K2)h_{cr} = h_c - h_r = 0.0785 K^2 - 0.0110 K^2 = 0.0675 K^2 \quad (\text{or } 0.067 K^2) Where KK is the horizontal sight distance in kilometers (K=D/1000K = D / 1000).

Application to Rod Readings

Because curvature dominates refraction, the optical line of sight strikes the rod too high. Therefore, the corrected rod reading is: Corrected Reading=Observed Reading−hcr\text{Corrected Reading} = \text{Observed Reading} - h_{cr} True Elevation=HI−(FS−hcr)=HI−FS+hcr\text{True Elevation} = \text{HI} - (\text{FS} - h_{cr}) = \text{HI} - \text{FS} + h_{cr}

Tip

In differential leveling loops, balancing backsight and foresight distances (∑KBS=∑KFS\sum K_{\text{BS}} = \sum K_{\text{FS}}) automatically cancels out curvature and refraction errors completely.


3. Three-Wire Leveling & Cross-Sectional Geometry

Three-Wire Leveling

Precision leveling employs reticles equipped with upper, middle, and lower stadia hairs:

  • Stadia Interval (ss): s=Upper Reading−Lower Readings = \text{Upper Reading} - \text{Lower Reading}.
  • Sight Distance (DD): D=100×sD = 100 \times s (using standard internal focusing instrument stadia constant Ks=100K_s = 100).
  • Rod Reading Check: To detect misreadings, ensure that Middle=Upper+Lower2\text{Middle} = \frac{\text{Upper} + \text{Lower}}{2}. The tolerance is typically ±1 mm\pm 1\text{ mm} to 2 mm2\text{ mm}.

Three-Level Road Cross-Sections

A standard three-level highway section consists of a roadbed of width bb, a centerline cut or fill cc, and slope stakes set where the side slopes intersect original ground:

  • Left slope stake: height hLh_L, horizontal offset wL=b2+sshLw_L = \frac{b}{2} + s_s h_L
  • Right slope stake: height hRh_R, horizontal offset wR=b2+sshRw_R = \frac{b}{2} + s_s h_R

Dividing the cross-section into four triangles (two with base cc and heights wL,wRw_L, w_R; two with base b/2b/2 and heights hL,hRh_L, h_R), the exact area formula is: A=12[c(wL+wR)+b2(hL+hR)]A = \frac{1}{2} \left[ c(w_L + w_R) + \frac{b}{2}(h_L + h_R) \right]


4. Earthwork Volume Computations

Between consecutive cross-sections separated by horizontal distance LL, earthwork volume is calculated via prismoidal geometry.

Average End-Area Method (VeV_e)

The average end-area formula assumes the volume varies linearly between end sections: Ve=(A1+A22)LV_e = \left( \frac{A_1 + A_2}{2} \right) L While simple and accepted in standard DPWH contracts, this formula systematically overestimates volume whenever A1≠A2A_1 \neq A_2.

Prismoidal Formula (VpV_p)

The mathematically rigorous prismoidal formula accounts for parabolic warping: Vp=L6(A1+4Am+A2)V_p = \frac{L}{6} \left( A_1 + 4 A_m + A_2 \right) Where AmA_m is the middle area computed from linear dimensions averaged between Station 1 and Station 2: cm=c1+c22,wL,m=wL1+wL22,wR,m=wR1+wR22,hL,m=hL1+hL22,hR,m=hR1+hR22c_m = \frac{c_1 + c_2}{2}, \quad w_{L, m} = \frac{w_{L1} + w_{L2}}{2}, \quad w_{R, m} = \frac{w_{R1} + w_{R2}}{2}, \quad h_{L, m} = \frac{h_{L1} + h_{L2}}{2}, \quad h_{R, m} = \frac{h_{R1} + h_{R2}}{2}

Caution

Never compute AmA_m as A1+A22\frac{A_1 + A_2}{2}! Doing so reduces the prismoidal formula back to the end-area formula.

Prismoidal Correction (CpC_p)

For three-level sections, the difference between the average end-area volume and the exact prismoidal volume is the Prismoidal Correction (CpC_p): Cp=Ve−Vp=L12(c1−c2)(w1−w2)C_p = V_e - V_p = \frac{L}{12} (c_1 - c_2)(w_1 - w_2) Where w1=(wL1+wR1)w_1 = (w_{L1} + w_{R1}) and w2=(wL2+wR2)w_2 = (w_{L2} + w_{R2}). Then: Vp=Ve−CpV_p = V_e - C_p Because (c1−c2)(c_1 - c_2) and (w1−w2)(w_1 - w_2) almost always share the same sign, CpC_p is positive, confirming that Ve>VpV_e > V_p.


5. Borrow Pit Excavation (Unit Base / Grid Method)

For large foundation excavations, basements, and site leveling, the site is staked in a uniform grid of square or rectangular prisms of horizontal area AunitA_{\text{unit}}. Corner elevations are surveyed before and after excavation to determine net vertical cuts (hh):

V=Aunit4(∑h1+2∑h2+3∑h3+4∑h4)V = \frac{A_{\text{unit}}}{4} \left( \sum h_1 + 2 \sum h_2 + 3 \sum h_3 + 4 \sum h_4 \right) Where:

  • h1h_1 = corner cuts common to exactly one grid unit
  • h2h_2 = corner cuts common to two adjoining grid units
  • h3h_3 = corner cuts common to three adjoining grid units
  • h4h_4 = corner cuts common to four adjoining grid units

6. Mass Haul Diagrams & Economic Analysis

A Mass Haul Diagram (MHD) is a continuous graphical plot of stationing along the horizontal axis versus cumulative net earthwork volume along the vertical axis.

Fundamental Properties of the Mass Curve

  • Cut (Excavation): Slopes upward from left to right (cumulative volume increases).
  • Fill (Embankment): Slopes downward from left to right (cumulative volume decreases).
  • Peak (Summit): Represents the transition point from cut to fill (maximum net cut).
  • Trough (Valley): Represents the transition point from fill to cut (maximum net fill).
  • Zero Slope: Indicates zero earthwork cross-section (grade point).

Haul Terminology and Definitions

  • Haul: The work performed in transporting material: Haul=Volume (m3)×Distance (m or stations)\text{Haul} = \text{Volume } (\text{m}^3) \times \text{Distance } (\text{m or stations}).
  • Balance Line: A horizontal line drawn across the mass curve. Where the balance line intersects the curve, cut equals fill (Vcut=VfillV_{\text{cut}} = V_{\text{fill}}).
  • Free Haul Distance (FHD): The fixed distance (typically 50 m50\text{ m} to 100 m100\text{ m}, or 4 to 5 stations4\text{ to }5\text{ stations}) within which the contractor moves earth without additional transportation payment beyond the unit excavation price.
  • Overhaul (OH): The transport of excavated material beyond the Free Haul Distance. Measured in m3-stations\text{m}^3\text{-stations} or m3-meters\text{m}^3\text{-meters}: Overhaul Distance (OHD)=Distance between Cut/Fill Centroids−FHD\text{Overhaul Distance (OHD)} = \text{Distance between Cut/Fill Centroids} - \text{FHD}
  • Borrow: Material imported from external borrow pits when local excavation is insufficient to satisfy embankment requirements.
  • Waste: Surplus excavation dumped off-site when cut exceeds nearby fill and haul costs exceed disposal costs.

Limit of Economic Haul (LEH)

The Limit of Economic Haul (LEH) represents the maximum haul distance at which it is economically preferable to haul on-site cut rather than wasting the cut and purchasing borrow: Cost of Hauling=Cost of Borrowing\text{Cost of Hauling} = \text{Cost of Borrowing} Coverhaul×(LEH−FHD)=CborrowC_{\text{overhaul}} \times (\text{LEH} - \text{FHD}) = C_{\text{borrow}} LEH=FHD+CborrowCoverhaul\text{LEH} = \text{FHD} + \frac{C_{\text{borrow}}}{C_{\text{overhaul}}} Where:

  • CborrowC_{\text{borrow}} = unit cost of borrow (₱/m3\text{m}^3)
  • CoverhaulC_{\text{overhaul}} = unit cost of overhaul (₱/m3-station\text{m}^3\text{-station} or ₱/m3-m\text{m}^3\text{-m})

7. Comprehensive Worked Example

Problem: Three-Level Earthwork & Prismoidal Correction

Given: Two consecutive highway stations separated by L=60.00 mL = 60.00\text{ m} have base width b=10.00 mb = 10.00\text{ m}, side slopes 1.5:11.5:1 (H:V). The cross-section notes are:

  • Station 1+000: Center cut c1=2.00 mc_1 = 2.00\text{ m}, Left: hL1wL1=1.60 m7.40 m\frac{h_{L1}}{w_{L1}} = \frac{1.60\text{ m}}{7.40\text{ m}}, Right: hR1wR1=2.40 m8.60 m\frac{h_{R1}}{w_{R1}} = \frac{2.40\text{ m}}{8.60\text{ m}}
  • Station 1+060: Center cut c2=3.60 mc_2 = 3.60\text{ m}, Left: hL2wL2=2.80 m9.20 m\frac{h_{L2}}{w_{L2}} = \frac{2.80\text{ m}}{9.20\text{ m}}, Right: hR2wR2=4.20 m11.30 m\frac{h_{R2}}{w_{R2}} = \frac{4.20\text{ m}}{11.30\text{ m}}

Calculate:

  1. Cross-sectional areas A1A_1 and A2A_2.
  2. Volume by Average End-Area (VeV_e).
  3. Prismoidal correction CpC_p and exact prismoidal volume VpV_p.

Solution:

  1. Area Computations: A1=12[c1(wL1+wR1)+b2(hL1+hR1)]A_1 = \frac{1}{2} \left[ c_1(w_{L1} + w_{R1}) + \frac{b}{2}(h_{L1} + h_{R1}) \right] w1=7.40+8.60=16.00 m,hL1+hR1=1.60+2.40=4.00 mw_1 = 7.40 + 8.60 = 16.00\text{ m}, \quad h_{L1} + h_{R1} = 1.60 + 2.40 = 4.00\text{ m} A1=12[2.00(16.00)+5.00(4.00)]=12[32.00+20.00]=26.00 m2A_1 = \frac{1}{2} [ 2.00(16.00) + 5.00(4.00) ] = \frac{1}{2} [ 32.00 + 20.00 ] = 26.00\text{ m}^2

    A2=12[c2(wL2+wR2)+b2(hL2+hR2)]A_2 = \frac{1}{2} \left[ c_2(w_{L2} + w_{R2}) + \frac{b}{2}(h_{L2} + h_{R2}) \right] w2=9.20+11.30=20.50 m,hL2+hR2=2.80+4.20=7.00 mw_2 = 9.20 + 11.30 = 20.50\text{ m}, \quad h_{L2} + h_{R2} = 2.80 + 4.20 = 7.00\text{ m} A2=12[3.60(20.50)+5.00(7.00)]=12[73.80+35.00]=54.40 m2A_2 = \frac{1}{2} [ 3.60(20.50) + 5.00(7.00) ] = \frac{1}{2} [ 73.80 + 35.00 ] = 54.40\text{ m}^2

  2. Average End-Area Volume (VeV_e): Ve=(A1+A22)L=(26.00+54.402)×60.00=40.20×60.00=2,412.00 m3V_e = \left( \frac{A_1 + A_2}{2} \right) L = \left( \frac{26.00 + 54.40}{2} \right) \times 60.00 = 40.20 \times 60.00 = 2,412.00\text{ m}^3

  3. Prismoidal Correction and Volume (VpV_p): Cp=L12(c1−c2)(w1−w2)C_p = \frac{L}{12}(c_1 - c_2)(w_1 - w_2) c1−c2=2.00−3.60=−1.60 mc_1 - c_2 = 2.00 - 3.60 = -1.60\text{ m} w1−w2=16.00−20.50=−4.50 mw_1 - w_2 = 16.00 - 20.50 = -4.50\text{ m} Cp=60.0012(−1.60)(−4.50)=5.00×7.20=36.00 m3C_p = \frac{60.00}{12} (-1.60)(-4.50) = 5.00 \times 7.20 = 36.00\text{ m}^3 Vp=Ve−Cp=2,412.00−36.00=2,376.00 m3V_p = V_e - C_p = 2,412.00 - 36.00 = 2,376.00\text{ m}^3

Loading diagram...
Mass Haul Diagram Profile and Haul Optimization Balance
Test Your Knowledge

A differential leveling circuit across a river valley records a backsight of 2.450 m on BM-A (Elevation 125.400 m) at a sight distance of 100 m. A foresight of 1.820 m is taken on TP-1 situated 1,200 m away on the opposite ridge. Incorporating the combined curvature and refraction correction (h_cr = 0.0675 K²), what is the corrected elevation of TP-1?

A

126.127 m

B

126.030 m

C

125.933 m

D

127.850 m

Test Your Knowledge

At Station 10+000, a three-level road cut has a center cut c = 2.40 m and roadbed width b = 10.0 m. The left slope stake reads h_L = 1.80 m at offset w_L = 7.70 m, and the right slope stake reads h_R = 3.20 m at offset w_R = 9.80 m. What is the cross-sectional area of the cut?

A

42.00 m²

B

33.50 m²

C

26.85 m²

D

46.00 m²

Test Your Knowledge

In a road grading project, the unit cost of borrow excavation is ₱180.00/m³, while the cost of overhaul is ₱15.00/m³-station (where 1 station = 20 meters). The contract specifies a Free Haul Distance (FHD) of 4 stations (80 meters). What is the Limit of Economic Haul (LEH)?

A

8 stations (160 m)

B

12 stations (240 m)

C

16 stations (320 m)

D

20 stations (400 m)

Sections you finish are checked off in the contents.