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.
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.
| Course | D (ft) | D · σ_α (ft) | σ_lat (ft) | σ_dep (ft) |
|---|---|---|---|---|
| A–B | 423.12 | 0.02051 | 0.02011 | 0.02041 |
| B–C | 519.74 | 0.02520 | 0.02368 | 0.02177 |
| C–D | 431.68 | 0.02093 | 0.02010 | 0.02083 |
| D–E | 379.80 | 0.01841 | 0.01864 | 0.01979 |
| E–A | 480.37 | 0.02329 | 0.02025 | 0.02307 |
| Σ of squares | 0.0021266 | 0.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,600The 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 ✓ identicalCompare 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.
| Quantity | Predicted | Observed |
|---|---|---|
| Closure in latitude (ft) | ±0.0461 (1σ) | +0.0488 |
| Closure in departure (ft) | ±0.0474 (1σ) | −0.0467 |
| Linear misclosure (ft) | 0.0661 | 0.0676 |
| Relative precision | 1 : 33,800 | 1 : 33,074 |
| Observed / predicted | 1.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- Standard deviation, standard error and a rejection decision on repeated distancesTen EDM observations of one baseline, reduced to a residuals table, a standard deviation and a standard error of the mean — and one suspect observation tested against a criterion that can actually reject it.
- Error in a taped distance measured in n tape lengths, against a single EDM shotThe same 648 ft line measured two ways. Series propagation over seven tape settings against one EDM observation with a constant-plus-ppm specification, and what the comparison says about when taping is still adequate.
- Error propagated into a computed area from errors in two sidesA rectangular parcel measured on two sides, and the uncertainty those two measurements put into the acreage — with the fractional form of the same result as a check and the worst case as a bound.