Error propagation

Predicted traverse closure error against the observed misclosure

What misclosure should a five-course traverse produce, given the instrument that ran it? Propagate the per-course errors, predict the closure, and compare it with what the field actually returned.

Analyze· about 30 minutes by hand· 5 steps

Given

  • A five-course closed boundary traverse, observed with a total station:
  • A–B azimuth 26°59′20″, 423.12 ft
  • B–C azimuth 124°10′45″, 519.74 ft
  • C–D azimuth 199°03′50″, 431.68 ft
  • D–E azimuth 248°22′05″, 379.80 ft
  • E–A azimuth 344°32′50″, 480.37 ft
  • Each distance carries a standard error of σ_D = ±0.020 ft, and each course azimuth a standard error of σ_α = ±10″, treated as independent between courses.
  • The traverse as observed returned ΣLat = +0.0488 ft and ΣDep = −0.0467 ft.
  • Distances in US survey feet, perimeter 2234.71 ft.

Required

  • The predicted standard error of the latitude and departure of each course.
  • The predicted linear misclosure of the traverse and the precision it implies.
  • A judgement on whether the observed misclosure is consistent with that prediction, or points to a blunder.

Work it through yourself before reading on — the solution below shows every step, so there is no way to skim it without giving the answer away.

Worked solution

Per-course propagation

Each course contributes error to the closure through both its distance and its azimuth. The propagation into latitude and departure is the same as for a single course: partial derivative times standard error, combined in quadrature.

Convert the angular standard error to radians first: 10″ divided by 206,264.81 gives 4.84814 × 10⁻⁵ rad. Multiplied by a course length this is the across-course displacement, and it is the term that grows with the size of the traverse.

On these courses the two sources contribute comparably. The distance term is a flat 0.020 ft, and the angular term D·σ_α runs from 0.018 to 0.025 ft — so neither the distances nor the angles dominate, which is what makes the compass rule the right adjustment for this traverse.

σ_α = 10″ / 206,264.81 = 4.84814 × 10⁻⁵ rad σ_lat = √( (cos α · σ_D)² + (D sin α · σ_α)² ) σ_dep = √( (sin α · σ_D)² + (D cos α · σ_α)² )

Course by course

Work the two propagations for each of the five courses and square the results, because it is the variances that will be summed.

Note that σ_lat and σ_dep differ course by course according to the azimuth: on a course running nearly north the angular error lands almost entirely in the departure, and on a course running nearly east it lands almost entirely in the latitude. Across five courses spread around a closed figure those preferences largely even out, which is why the two summed variances come out close to one another.

Predicted per-course standard errors
CourseD (ft)D · σ_α (ft)σ_lat (ft)σ_dep (ft)
A–B423.120.020510.020110.02041
B–C519.740.025200.023680.02177
C–D431.680.020930.020100.02083
D–E379.800.018410.018640.01979
E–A480.370.023290.020250.02307
Σ of squares0.00212660.0022486

Predicted misclosure

The closure in latitude is the sum of five independent latitudes, so its variance is the sum of their variances; likewise for departure. The linear misclosure is the resultant of the two, so its variance is the sum of both.

That gives a predicted standard misclosure of 0.066 ft, and a predicted precision of about 1:33,800. This is what the instrument and the geometry can deliver — not a target, but an expectation.

Understand what the prediction is and is not. It is a root-mean-square figure for a two-dimensional quantity, so a substantial share of otherwise blameless traverses will close worse than it. The value at which a traverse becomes suspect is 1.73 times it — the 95 per cent point of the two-dimensional misclosure distribution — which is 0.114 ft here.

Σσ_lat² = 0.0021266 ft² Σσ_dep² = 0.0022486 ft² e_pred = √( 0.0021266 + 0.0022486 ) = √0.0043752 = 0.0661 ft predicted precision = 2234.71 / 0.0661 = 1 : 33,800 95 % two-dimensional bound = 1.73 × 0.0661 = 0.114 ft, i.e. 1 : 19,600

The whole sum collapses

There is a simplification worth knowing, because it turns a twenty-line table into one line. Add σ_lat² and σ_dep² for a single course and the trigonometry cancels: the distance terms give cos²α + sin²α = 1 times σ_D², and the angular terms give sin²α + cos²α = 1 times (D σ_α)².

So the predicted misclosure of any traverse is the square root of the sum, over all courses, of σ_D² + (D σ_α)². The azimuths drop out entirely. Only the lengths and the two standard errors matter.

Use the long table when a per-course breakdown is wanted — to see which course is the weak one — and the collapsed form when only the total is needed. They must agree exactly, which makes each a check on the other.

σ_lat² + σ_dep² = σ_D² (cos²α + sin²α) + (D σ_α)² (sin²α + cos²α) = σ_D² + (D σ_α)² Σ σ_D² = 5 × 0.020² = 5 × 0.000400 = 0.002000 Σ (D σ_α)² = 0.02051² + 0.02520² + 0.02093² + 0.01841² + 0.02329² = 0.0023752 e_pred = √( 0.002000 + 0.0023752 ) = √0.0043752 = 0.0661 ft ✓ identical

Compare with the field

The traverse actually closed with ΣLat = +0.0488 ft and ΣDep = −0.0467 ft, giving an observed misclosure of 0.0676 ft and a precision of 1:33,074.

Observed over predicted is 1.02. The traverse closed almost exactly as well as the instrument and the geometry say it should — the strongest possible evidence that the misclosure is accumulated random error and contains no blunder. There is nothing to look for and the compass rule may be applied with confidence.

Read the comparison the other way to see its value. Had the same traverse closed at 0.30 ft — a precision of 1:7,400, which many specifications would still pass — that would be 4.5 times the predicted misclosure, which accumulated random error alone would essentially never produce. A specification would have accepted it; the prediction would have caught it. This is the discipline the comparison is for: a tolerance tells you whether a traverse is good enough, and only a prediction tells you whether it is behaving.

Predicted against observed
QuantityPredictedObserved
Closure in latitude (ft)±0.0461 (1σ)+0.0488
Closure in departure (ft)±0.0474 (1σ)−0.0467
Linear misclosure (ft)0.06610.0676
Relative precision1 : 33,8001 : 33,074
Observed / predicted1.02

Answer

  • Predicted per-course standard errors run 0.019 to 0.024 ft in latitude and 0.020 to 0.023 ft in departure, with Σσ_lat² = 0.0021266 ft² and Σσ_dep² = 0.0022486 ft².
  • Predicted linear misclosure = √(0.0021266 + 0.0022486) = 0.0661 ft, a predicted precision of about 1 : 33,800.
  • Observed misclosure = √(0.0488² + 0.0467²) = 0.0676 ft, a precision of 1 : 33,074.
  • Observed / predicted = 1.02. The traverse closed exactly as well as its instrumentation predicts, so the misclosure is accumulated random error and contains no blunder. Adjust by the compass rule.
  • The whole prediction collapses to e = √Σ(σ_D² + (D σ_α)²): the azimuths cancel out, and only the course lengths and the two standard errors matter.

Check

The collapsed one-line form and the five-row table both give 0.0661 ft, from arithmetic that shares no intermediate quantity except the course lengths. A slip in any single σ_lat or σ_dep would break the agreement.

Check the two summed variances against each other: Σσ_lat² = 0.0021266 and Σσ_dep² = 0.0022486, differing by only 5 per cent. For a closed figure whose courses run in all directions, the north–south and east–west variances should be comparable; a large disparity would indicate the traverse is a long thin loop rather than a compact one, and would be worth understanding before adjusting.

Check the plausibility of the observed misclosure statistically. In two dimensions the ratio of an observed misclosure to its predicted root-mean-square value exceeds 1.73 only 5 per cent of the time. The observed ratio is 1.02, so the closure sits in the middle of the expected distribution.

Check the angular contribution scale: 10″ over an average course of 447 ft displaces the far end by 447 × 10 / 206,264.81 = 0.0217 ft, which is the size of the D·σ_α column. If that column were an order of magnitude out, the radian conversion would be the place to look.

More error propagation

All in this category