地球与空间科学 · 第 3 册 · 本科一年级

University Earth and Space Sciences — Year 1

书籍暂未推出简体中文版——当前显示英文版。

University Earth and Space Sciences — Year 1 · 本科一年级

6Field and Lab Methods I

A geologist mapping the Varden Hills carries little: a notebook, a pencil, a satellite receiver, a compass-clinometer, a hand lens and a few sample bags. With them she fixes where she stands, measures the attitude of the beds, names the rocks, and brings back the samples that the laboratory will cut, grind and look at through a microscope. Every number of the earlier chapters — an attitude on a map, a dip in a section — began as a reading of this kind, with its own error. This chapter describes how professionals make these observations, and how repeated readings are averaged and their scatter turned into an uncertainty. It ends with a season of readings on the east limb of the Varden anticline.

You already know

The attitude of a plane is written dip direction/dip, as 081/12081/12, or as a strike by the right-hand rule with a dip (Chapter 4). The magnetic declination is the angle between geographic north and the direction in which a compass needle points (physics, grade 11). Minerals are recognised in the hand by their lustre, cleavage, hardness and colour (Book 2, grade 11). The mean and the standard deviation of a set of numbers summarise their centre and their scatter (mathematics).

6.1 The notebook and the map

Definition 6.1 (Field station, field notebook)

A field station is a numbered place where a geologist makes and records observations: its position, the rocks seen, their attitude, and the samples taken. The field notebook is the bound record of the stations, in the order they were visited, with their numbers, sketches and measurements, written on the spot.

The notebook is the primary record: what is not written at the station is lost, and a map drawn in the evening is only as good as the notes behind it. Each station has a number that appears on the field map, in the notebook, and on every sample bag taken there. A sketch shows what a list of numbers cannot — how two units meet, which way younger beds lie, where a reading was taken on a folded surface — and records its scale and orientation. Facts and interpretations are kept apart: “limestone on shale, contact sharp” is an observation; “unconformity” is an inference, written as such.

In the field — A day of mapping

A mapping geologist plans a traverse across the strike, so as to meet every unit in turn, and stops wherever a contact, a change of rock or a measurable surface appears. At each station she writes the station number and the position, sketches what she sees, measures attitudes, and notes the samples. Positions come from a satellite receiver: an ordinary receiver, such as a smartphone’s, is typically accurate to within about 4.9 m4.9\,\mathrm{m} under open sky, worse near cliffs and trees, while survey receivers using two frequencies or correction services reach a few centimetres. At a map scale of 1:10 000, 4.9 m4.9\,\mathrm{m} is half a millimetre: an ordinary receiver is good enough for mapping, but not for measuring the slow motion of the ground (Chapter 12).

6.2 Measuring orientations

The attitude of a plane — a bed, a fault, a joint — is measured with a compass-clinometer: a magnetic compass to read azimuths, and a pendulum or level to read angles from the horizontal. Its edge is laid along the strike of the surface to read the strike, or the instrument is set on the surface along its steepest line to read the dip and, from the azimuth of that line, the dip direction.

Measuring an attitude (schematic). (a) The edge of the compass laid horizontally on the bed gives the strike; the dip direction is perpendicular to it, down the bed. (b) Set on its edge along the steepest line, the clinometer gives the dip .
Measuring an attitude (schematic). (a) The edge of the compass laid horizontally on the bed gives the strike; the dip direction is perpendicular to it, down the bed. (b) Set on its edge along the steepest line, the clinometer gives the dip δ\delta.

A line — a fold axis, a scratch on a fault surface, the long axis of pebbles — has an orientation of its own.

Definition 6.2 (Trend and plunge)

The trend of a line is the azimuth of its horizontal projection, measured toward the end of the line that points downward. Its plunge is the angle between the line and the horizontal, measured in the vertical plane that contains it. A line plunging 20∘20{}^{\circ} toward 015∘015{}^{\circ} is written 20→01520\to015.

A line of trend tt and plunge pp has the downward unit vector (cos⁡psin⁡t,cos⁡pcos⁡t,−sin⁡p)(\cos p\sin t, \cos p\cos t, -\sin p) in east, north and up components; a line lying in a plane of dip direction α\alpha and dip δ\delta plunges by the apparent dip of that plane in the vertical plane of azimuth tt, so that tan⁡p=tan⁡δcos⁡(α−t)\tan p = \tan\delta\cos(\alpha - t) (Theorem 5.4).

6.2.1 The magnetic declination

A compass needle points to magnetic north, not to geographic north; the angle between the two is the magnetic declination, counted positive when magnetic north lies east of geographic north. Maps and the attitudes of this book use geographic azimuths.

Method 6.3 (Correcting a compass azimuth)

  1. Find the magnetic declination DD for the place and the date, from a model of the Earth’s field (east positive, west negative).
  2. Add it to every azimuth read on the compass: geographic azimuth == magnetic azimuth +D+ D, taken modulo 360∘360{}^{\circ}.
  3. Leave dips and plunges alone: they are measured from the horizontal, with gravity, not with the needle.

Most field compasses can be set to add DD mechanically; the setting, and the year it was made for, are written in the notebook.

Figures as of 2026

On 1 January 2026, the World Magnetic Model (version 2025) gives the magnetic declination as −0.47∘-0.47{}^{\circ} (west) at Edinburgh, changing by +0.20∘+0.20{}^{\circ} a year; −11.28∘-11.28{}^{\circ} at Reykjavik, changing by +0.30∘+0.30{}^{\circ} a year; +14.20∘+14.20{}^{\circ} (east) at Anchorage, changing by −0.28∘-0.28{}^{\circ} a year; and +23.23∘+23.23{}^{\circ} at Wellington. A geologist who forgot the correction at Wellington would misplace every azimuth by more than 23∘23{}^{\circ}; one who reused an Edinburgh setting from ten years before would carry an error of about 2∘2{}^{\circ}.

6.2.2 Repeated readings

A single reading of a bedding plane carries errors: the surface is never a perfect plane, the instrument is not perfectly set, and the reading is rounded to the degree. A geologist therefore measures the same surface several times, or several surfaces of the same limb, and averages. Averaging azimuths needs care: the readings 355∘355{}^{\circ} and 005∘005{}^{\circ} describe directions 10∘10{}^{\circ} apart, but their arithmetic mean, 180∘180{}^{\circ}, points the opposite way. The right mean works with vectors.

Definition 6.4 (Mean orientation (vector mean))

For a plane of dip direction α\alpha and dip δ\delta, let

n⃗=(sin⁡δsin⁡α, sin⁡δcos⁡α, cos⁡δ)\vec n = (\sin\delta\sin\alpha,\ \sin\delta\cos\alpha,\ \cos\delta)

be its upward unit normal in east, north and up components. The mean orientation (or vector mean) of NN readings of normals n⃗1,…,n⃗N\vec n_1, \dots, \vec n_N is the plane whose normal is m⃗=S⃗/∥S⃗∥\vec m = \vec S/\|\vec S\|, where S⃗=∑in⃗i\vec S = \sum_i \vec n_i is their resultant; Rˉ=∥S⃗∥/N\bar R = \|\vec S\|/N is the mean resultant length. The mean dip direction and dip are αˉ=atan2⁡(SE,SN)\bar\alpha = \operatorname{atan2}(S_E, S_N) and δˉ=arccos⁡(SU/∥S⃗∥)\bar\delta = \arccos(S_U/\|\vec S\|).

That n⃗\vec n is the normal follows from Chapter 4: the plane z=−tan⁡δ (xsin⁡α+ycos⁡α)z = -\tan\delta\,(x\sin\alpha + y\cos\alpha) has the upward normal (tan⁡δsin⁡α,tan⁡δcos⁡α,1)(\tan\delta\sin\alpha, \tan\delta\cos\alpha, 1), whose length is 1/cos⁡δ1/\cos\delta.

Proposition 6.5 (The vector mean is the closest direction)

Among all unit vectors u⃗\vec u, the mean normal m⃗\vec m maximises ∑in⃗i⋅u⃗=∑icos⁡θi(u⃗)\sum_i \vec n_i\cdot\vec u = \sum_i \cos\theta_i(\vec u), where θi(u⃗)\theta_i(\vec u) is the angle between n⃗i\vec n_i and u⃗\vec u; equivalently, it minimises ∑i∥n⃗i−u⃗∥2\sum_i \|\vec n_i - \vec u\|^2. Moreover

Rˉ=1N∑icos⁡θi,θi=θi(m⃗),\bar R = \frac1N\sum_i \cos\theta_i, \qquad \theta_i = \theta_i(\vec m),

so that 0<Rˉ≤10 < \bar R \le 1, with Rˉ=1\bar R = 1 only if all readings agree, and, when the scatter is small, the root-mean-square angle between the readings and their mean is

θrms≈2(1−Rˉ)(radians).\theta_{\text{rms}} \approx \sqrt{2(1 - \bar R)} \quad\text{(radians)}.

Proof. ∑in⃗i⋅u⃗=S⃗⋅u⃗≤∥S⃗∥\sum_i \vec n_i\cdot\vec u = \vec S\cdot\vec u \le \|\vec S\| by the Cauchy–Schwarz inequality, with equality exactly when u⃗=m⃗\vec u = \vec m. Since ∥n⃗i−u⃗∥2=2−2 n⃗i⋅u⃗\|\vec n_i - \vec u\|^2 = 2 - 2\,\vec n_i\cdot\vec u, the two problems are the same. With u⃗=m⃗\vec u = \vec m the bound is reached: ∥S⃗∥=∑icos⁡θi\|\vec S\| = \sum_i\cos\theta_i, which is the formula for Rˉ\bar R. Finally cos⁡θ≈1−θ2/2\cos\theta \approx 1 - \theta^2/2 for small θ\theta, so 1−Rˉ≈1N∑iθi2/2=θrms2/21 - \bar R \approx \frac1N\sum_i\theta_i^2/2 = \theta_{\text{rms}}^2/2. ∎

How do errors in the two angles move the normal? Differentiating n⃗\vec n,

∂n⃗∂δ=(cos⁡δsin⁡α, cos⁡δcos⁡α, −sin⁡δ),∂n⃗∂α=(sin⁡δcos⁡α, −sin⁡δsin⁡α, 0):\frac{\partial\vec n}{\partial\delta} = (\cos\delta\sin\alpha,\ \cos\delta\cos\alpha,\ -\sin\delta), \qquad \frac{\partial\vec n}{\partial\alpha} = (\sin\delta\cos\alpha,\ -\sin\delta\sin\alpha,\ 0):

the first has length 1, the second length sin⁡δ\sin\delta, and the two are perpendicular. Small errors Δδ\Delta\delta and Δα\Delta\alpha therefore turn the normal by the angle

Δθ≈Δδ2+sin⁡2δ Δα2.\Delta\theta \approx \sqrt{\Delta\delta^2 + \sin^2\delta\,\Delta\alpha^2}.

On a gently dipping bed an error in the dip direction hardly changes the plane: on a bed dipping 12∘12{}^{\circ}, a 3∘3{}^{\circ} error in dip direction turns the plane by only 3sin⁡12∘=0.6∘3\sin12{}^{\circ} = 0.6{}^{\circ}. The converse is the field geologist’s difficulty: the dip direction of a gentle bed is poorly fixed by its plane, and small undulations of the surface scatter it widely. The standard error of a mean of NN independent readings of standard deviation ss is s/Ns/\sqrt N (mathematics): to halve it, four times as many readings are needed.

Twenty-five simulated readings of the bedding of the Varden east limb, whose model attitude is 080.7/12.2, each with errors of 3 standard deviation on both angles, rounded to the degree. Solid line: the mean; dashed: one standard deviation either side (080.8 ± 2.9 and 12.4 ± 2.1). The vector mean of the 25 readings is 080.9/12.3.
Twenty-five simulated readings of the bedding of the Varden east limb, whose model attitude is 080.7/12.2080.7/12.2, each with errors of 3∘3{}^{\circ} standard deviation on both angles, rounded to the degree. Solid line: the mean; dashed: one standard deviation either side (080.8±2.9080.8 \pm 2.9 and 12.4±2.112.4 \pm 2.1). The vector mean of the 25 readings is 080.9/12.3080.9/12.3.

6.3 The hand lens and field identification

Most rocks are named in the field, before any laboratory work, from the minerals a hand lens shows and from the way the rock breaks. The geologist looks at a fresh surface, since weathering coats and discolours the rock, and notes the size of the grains (whether they can be seen with the naked eye, with the lens, or not at all), their shapes, whether they interlock like a mosaic or are cemented like sand, and any alignment. The minerals are recognised by the properties recalled at the start of the chapter: the glassy lustre and absence of cleavage of quartz, the two cleavages and dull whiteness of feldspar, the flakes of mica, the rhombic cleavage of calcite. Hardness is tried against a steel point: quartz is harder than steel and calcite much softer.

Professionals carry a small dropper bottle of dilute hydrochloric acid: a carbonate rock releases bubbles of carbon dioxide where the drop falls (chemistry). Calcite reacts briskly; dolomite reacts weakly unless powdered. The acid is handled with gloves and eye protection, and never by a reader of this book: the description is of professional practice.

6.4 Samples

Definition 6.6 (Oriented sample)

An oriented sample is a rock sample on which the geologist has marked, before removing it, the attitude of one of its surfaces (a strike line and a dip arrow, or the trend of a line), so that the sample can be put back in its field position in the laboratory.

A sample is useful only if it is representative of the unit it stands for — fresh, not weathered, not taken from a vein or a fallen block — and if its origin is certain. Each one receives the station number and a sample number, written on the rock and on the bag, and is listed in the notebook with its position, its unit and the reason it was taken. Samples for chemical or isotopic analysis, whose results may be published decades later, are tracked from the field to the laboratory by a written record of every hand they pass through — a chain of custody — so that no doubt remains about which rock gave which number. Oriented samples are needed whenever the direction of something inside the rock matters: an alignment of minerals, a fabric seen under the microscope, or the direction of the magnetisation the rock acquired when it formed.

6.5 Thin sections and the polarising microscope

Definition 6.7 (Thin section)

A thin section is a slice of rock, glued to a glass slide and ground to a standard thickness of 30 µm30\,\text{µ}\mathrm{m}, thin enough for most minerals to transmit light.

In the lab — Making a thin section

In the preparation laboratory, a technician cuts a small block from the sample with a diamond saw, grinds one face perfectly flat, and glues it to a glass slide with an epoxy resin whose refractive index, 1.54 to 1.55, is almost that of quartz. Most of the block is then sawn off, and the slice is ground down on a rotating lap with ever finer abrasive until it is 30 µm30\,\text{µ}\mathrm{m} thick. The technician checks the thickness under the microscope from the colours shown by quartz grains (Chapter 8), and finishes with a thin cover glass or a polished surface. Porous or crumbly rocks are first soaked in resin so that they hold together.

The thin section is studied with a microscope fitted with two polarising filters. Linear polarisation is the subject of physics: here it is enough that a polarising filter transmits only light vibrating in one direction, and that two filters set at right angles transmit nothing.

Definition 6.8 (Polarising microscope, plane-polarised light, crossed polars)

A polarising microscope is a transmitted light microscope with a polarising filter (the polariser) below a stage that rotates through 360∘360{}^{\circ} against a graduated scale, and a second polarising filter (the analyser) between the objective and the eyepiece, which can be moved in or out of the light path, its vibration direction at right angles to the polariser’s. With the analyser out, the specimen is seen in plane-polarised light; with the analyser in, it is seen between crossed polars.

The light path of a polarising microscope (schematic, from bottom to top). The polariser and the analyser vibrate at right angles (double arrows); with nothing between them, crossed polars give a black field.
The light path of a polarising microscope (schematic, from bottom to top). The polariser and the analyser vibrate at right angles (double arrows); with nothing between them, crossed polars give a black field.

Between crossed polars, the glass of the slide, the epoxy and any isotropic mineral (one whose optical properties are the same in every direction, as in cubic crystals) stay black whatever the rotation of the stage. Most minerals are not isotropic: they change the polarisation of the light that crosses them, and appear in interference colours that change as the stage turns and vanish four times in a full turn. Why they do so, and how the colours and the angles identify a mineral, is the subject of Chapter 8. Minerals in thin section are usually examined at magnifications from about 16×16\times to 400×400\times.

A granite in thin section between crossed polars. Biotite is brown, muscovite shows bright pink and green interference colours, plagioclase is striped by its twinning, and quartz is white to grey. Photograph Darla Sondrol, CC0 (Wikimedia Commons).
A granite in thin section between crossed polars. Biotite is brown, muscovite shows bright pink and green interference colours, plagioclase is striped by its twinning, and quartz is white to grey. Photograph Darla Sondrol, CC0 (Wikimedia Commons).

Safety

Field and laboratory work are made safe by professional rules, not by luck. In the field, geologists work in pairs or leave a written plan of their route and return time; they carry weather forecasts, maps and first aid, keep away from the foot and the edge of cliffs and quarry faces, and wear a helmet and eye protection wherever rock can fall or splinter; work on active sites follows the operator’s rules. In the laboratory, saws and laps run behind guards with water to keep dust down, resins and acids are handled with gloves and goggles under extraction, and every chemical has its safety data sheet. None of these activities is for a reader of this book alone: the chapter describes professional practice.

6.6 Exercises

Exercise 6.1 ★

A satellite receiver gives positions within about 4.9 m4.9\,\mathrm{m}. How large is that on a map at 1:10 000? At 1:25 000?

Solution

Solution of Exercise 6.1.

At 1:10 000, 4.9 m4.9\,\mathrm{m} becomes 4.9/10 000=0.49 mm4.9/10\,000 = 0.49\,\mathrm{mm}; at 1:25 000, 4.9/25 000=0.20 mm4.9/25\,000 = 0.20\,\mathrm{mm}, about the width of a pencil line.

Exercise 6.2 ★

A compass with no declination setting reads the dip direction 052∘052{}^{\circ} on a bed. Give the geographic dip direction if the reading was made at Reykjavik on 1 January 2026, then at Wellington the same day (dated box).

Solution

Solution of Exercise 6.2.

At Reykjavik D=−11.28∘D = -11.28{}^{\circ}: 52−11.28=040.7∘52 - 11.28 = 040.7{}^{\circ}. At Wellington D=+23.23∘D = +23.23{}^{\circ}: 52+23.23=075.2∘52 + 23.23 = 075.2{}^{\circ}. The same compass reading means directions 34.5∘34.5{}^{\circ} apart.

Exercise 6.3 ★

A fold axis plunges 20∘20{}^{\circ} toward 015∘015{}^{\circ}. Write its orientation as in the text and give its downward unit vector in east, north and up components.

Solution

Solution of Exercise 6.3.

20→01520\to015. With t=15∘t = 15{}^{\circ} and p=20∘p = 20{}^{\circ}, the downward unit vector is

(cos⁡20∘sin⁡15∘, cos⁡20∘cos⁡15∘, −sin⁡20∘)=(0.2432, 0.9077, −0.3420).(\cos20{}^{\circ}\sin15{}^{\circ},\ \cos20{}^{\circ}\cos15{}^{\circ},\ -\sin20{}^{\circ}) = (0.2432,\ 0.9077,\ -0.3420).

Exercise 6.4 ★

Three readings of a vertical joint give the azimuths of its normal 350∘350{}^{\circ}, 358∘358{}^{\circ} and 006∘006{}^{\circ}. Compute their arithmetic mean, then the azimuth of the sum of the three horizontal unit vectors. Which mean is right?

Solution

Solution of Exercise 6.4.

The arithmetic mean is (350+358+6)/3=238∘(350 + 358 + 6)/3 = 238{}^{\circ}, nearly opposite to all three readings. The horizontal unit vectors (sin⁡a,cos⁡a)(\sin a, \cos a) add up to (−0.1736−0.0349+0.1045, 0.9848+0.9994+0.9945)=(−0.1040,2.9787)(-0.1736 - 0.0349 + 0.1045,\ 0.9848 + 0.9994 + 0.9945) = (-0.1040, 2.9787), whose azimuth is atan2⁡(−0.1040,2.9787)=−2.0∘\operatorname{atan2}(-0.1040, 2.9787) = -2.0{}^{\circ}, that is 358.0∘358.0{}^{\circ}: the right mean, in the middle of the readings.

Exercise 6.5 ★★

Five readings of the dip of one bed are 34∘34{}^{\circ}, 37∘37{}^{\circ}, 35∘35{}^{\circ}, 33∘33{}^{\circ} and 36∘36{}^{\circ}. Compute their mean, their standard deviation and the standard error of the mean.

Solution

Solution of Exercise 6.5.

Mean 175/5=35.0∘175/5 = 35.0{}^{\circ}. Deviations −1,2,0,−2,1-1, 2, 0, -2, 1; sum of squares 10; s=10/4=1.58∘s = \sqrt{10/4} = 1.58{}^{\circ}. Standard error 1.58/5=0.71∘1.58/\sqrt5 = 0.71{}^{\circ}.

Exercise 6.6 ★★

A geologist averages one reading on each limb of a symmetric fold: 090/30090/30 and 270/30270/30. Compute the vector mean and Rˉ\bar R. What has gone wrong?

Solution

Solution of Exercise 6.6.

The normals are (sin⁡30∘,0,cos⁡30∘)=(0.5,0,0.866)(\sin30{}^{\circ}, 0, \cos30{}^{\circ}) = (0.5, 0, 0.866) and (−0.5,0,0.866)(-0.5, 0, 0.866); S⃗=(0,0,1.732)\vec S = (0, 0, 1.732) is vertical: the mean is a horizontal plane, and Rˉ=1.732/2=0.866\bar R = 1.732/2 = 0.866, each reading 30∘30{}^{\circ} from the mean. The two readings do not sample one surface with errors: they belong to two limbs of a fold, two different populations, and their average describes neither. Each limb is averaged on its own; together they define the fold axis, the line common to both planes, here horizontal and north–south.

Exercise 6.7 ★★

Scratches on a bedding surface of attitude 080/12080/12 have the trend 120∘120{}^{\circ}. What is their plunge? What would it be for the trend 170∘170{}^{\circ}?

Solution

Solution of Exercise 6.7.

tan⁡p=tan⁡12∘cos⁡(80∘−120∘)=0.2126×0.7660=0.1628\tan p = \tan12{}^{\circ}\cos(80{}^{\circ} - 120{}^{\circ}) = 0.2126\times0.7660 = 0.1628: p=9.2∘p = 9.2{}^{\circ}, the scratches plunging toward 120∘120{}^{\circ}. For the trend 170∘170{}^{\circ}, cos⁡(−90∘)=0\cos(-90{}^{\circ}) = 0: the line is horizontal, along the strike.

Exercise 6.8 ★★

At Edinburgh the magnetic declination is −0.47∘-0.47{}^{\circ} on 1 January 2026 and changes by +0.20∘+0.20{}^{\circ} a year. Assuming the change stays steady, when does it reach zero, and what is it on 1 January 2030?

Solution

Solution of Exercise 6.8.

0.47/0.20=2.350.47/0.20 = 2.35 years after 1 January 2026, that is in spring 2028. On 1 January 2030: −0.47+4×0.20=+0.33∘-0.47 + 4\times0.20 = +0.33{}^{\circ}, east of geographic north.

Exercise 6.9 ★★

Three readings of one bed are 080/10080/10, 080/12080/12 and 080/14080/14. Show that their vector mean is 080/12080/12, compute Rˉ\bar R and the estimate 2(1−Rˉ)\sqrt{2(1 - \bar R)} of θrms\theta_{\text{rms}}, and compare with the exact root-mean-square angle.

Solution

Solution of Exercise 6.9.

The three normals lie in the vertical plane of azimuth 080∘080{}^{\circ}, at 10∘10{}^{\circ}, 12∘12{}^{\circ} and 14∘14{}^{\circ} from the vertical; the outer two are symmetric about the middle one, so their sum lies along it: the mean is 080/12080/12, with θ=2∘,0,2∘\theta = 2{}^{\circ}, 0, 2{}^{\circ}. Then Rˉ=(1+2cos⁡2∘)/3=(1+1.998782)/3=0.999594\bar R = (1 + 2\cos2{}^{\circ})/3 = (1 + 1.998782)/3 = 0.999594, and 2(1−Rˉ)=0.000812=0.028 50 rad=1.633∘\sqrt{2(1 - \bar R)} = \sqrt{0.000812} = 0.028\,50\,\mathrm{rad} = 1.633{}^{\circ}; exactly, θrms=(4+0+4)/3=1.633∘\theta_{\text{rms}} = \sqrt{(4 + 0 + 4)/3} = 1.633{}^{\circ}.

Exercise 6.10 ★★★

A thin section cuts a spherical grain of radius rr at a distance hh from its centre, with hh equally likely anywhere between 0 and rr. Show that the mean radius of the circle seen in the section is πr/4\pi r/4. Quartz grains 0.50 mm0.50\,\mathrm{mm} across appear on average how large in thin section?

Solution

Solution of Exercise 6.10.

The section is a circle of radius r2−h2\sqrt{r^2 - h^2}. Its mean is 1r∫0rr2−h2 dh=1r⋅πr24=πr4\frac1r\int_0^r\sqrt{r^2 - h^2}\,\mathrm dh = \frac1r\cdot\frac{\pi r^2}{4} = \frac{\pi r}{4}, since the integral is the area of a quarter disc. Grains 0.50 mm0.50\,\mathrm{mm} across appear on average 0.50×π/4=0.39 mm0.50\times\pi/4 = 0.39\,\mathrm{mm} across: a thin section underestimates grain sizes.

Exercise 6.11 ★★★

Show that if every dip direction of a set of readings is increased by the same angle DD, the dip direction of their vector mean increases by DD and its dip does not change. Why does this matter for the declination correction?

Solution

Solution of Exercise 6.11.

The rotation by DD about the vertical, (x,y,z)↦(xcos⁡D+ysin⁡D, −xsin⁡D+ycos⁡D, z)(x, y, z) \mapsto (x\cos D + y\sin D,\ -x\sin D + y\cos D,\ z), sends (sin⁡δsin⁡α,sin⁡δcos⁡α,cos⁡δ)(\sin\delta\sin\alpha, \sin\delta\cos\alpha, \cos\delta) to (sin⁡δsin⁡(α+D),sin⁡δcos⁡(α+D),cos⁡δ)(\sin\delta\sin(\alpha + D), \sin\delta\cos(\alpha + D), \cos\delta): it adds DD to every dip direction. Being linear, it sends the resultant S⃗\vec S to the resultant of the rotated normals; it keeps lengths and the vertical component, so the mean dip arccos⁡(SU/∥S⃗∥)\arccos(S_U/\|\vec S\|) is unchanged, and it turns the horizontal part of S⃗\vec S by DD, so the mean dip direction increases by DD. Correcting each reading for the declination, or only their mean, gives the same result.

Exercise 6.12 ★★★

On a bed dipping 12∘12{}^{\circ}, the normal of a single reading scatters with a standard deviation of 2∘2{}^{\circ} in each of two perpendicular directions: along the line of steepest dip, which changes the dip, and across it, which changes the dip direction. How many readings are needed for a standard error of 0.5∘0.5{}^{\circ} on the mean dip? Of 1∘1{}^{\circ} on the mean dip direction?

Solution

Solution of Exercise 6.12.

The scatter along the line of steepest dip is a scatter of the dip: the standard error is 2/N≤0.52/\sqrt N \le 0.5 for N≥16N \ge 16. Across it, a turn of the normal by 2∘2{}^{\circ} changes the dip direction by 2/sin⁡12∘=9.62∘2/\sin12{}^{\circ} = 9.62{}^{\circ}; the standard error of the mean dip direction is 9.62/N≤19.62/\sqrt N \le 1 for N≥92.5N \ge 92.5, that is 93 readings. A gentle dip is quickly known; its direction is not.

6.7 Problem: A field season on the Varden east limb

Problem 6.1

Weekend problem — averaging twelve compass readings, measuring their scatter, correcting them, and comparing them with the boreholes

A geologist measured the bedding of the east limb of the Varden anticline at twelve stations, with a compass set for the declination, rounding to the degree. The table gives each reading (dip direction/dip) and the components of its upward unit normal n⃗=(nE,nN,nU)\vec n = (n_E, n_N, n_U), rounded to four decimals.

readingnEn_EnNn_NnUn_UreadingnEn_EnNn_NnUn_U
1076/11076/110.18510.04620.98167083/11083/110.18940.02330.9816
2081/13081/130.22220.03520.97448083/12083/120.20640.02530.9781
3078/16078/160.26960.05730.96139085/12085/120.20710.01810.9781
4081/11081/110.18850.02980.981610078/11078/110.18660.03970.9816
5076/10076/100.16850.04200.984811081/13081/13
6083/14083/140.24010.02950.970312082/11082/11

Part I — The mean attitude.

  1. Compute the arithmetic means of the twelve dip directions and of the twelve dips.
  2. Give an example of two readings for which the arithmetic mean of the dip directions is absurd, and say why it is safe here.
  3. Show that n⃗=(sin⁡δsin⁡α,sin⁡δcos⁡α,cos⁡δ)\vec n = (\sin\delta\sin\alpha, \sin\delta\cos\alpha, \cos\delta) is a unit vector perpendicular to the dip line and to the strike line of the plane α/δ\alpha/\delta.
  4. Complete the table for readings 11 and 12.
  5. Compute the resultant S⃗\vec S of the twelve normals.
  6. Compute ∥S⃗∥\|\vec S\| and the mean resultant length Rˉ\bar R.
  7. Deduce the vector mean: its dip direction and its dip.

Part II — The scatter.

  1. Compute the sample standard deviations of the dip directions and of the dips.
  2. Using Proposition 6.5, estimate the root-mean-square angle θrms\theta_{\text{rms}} between the readings and their mean, in degrees.
  3. Predict θrms\theta_{\text{rms}} from the two standard deviations of question 8 and the error formula of the chapter, and compare.
  4. Which of the two standard deviations contributes most to θrms\theta_{\text{rms}}, and why?
  5. Compute the standard errors of the mean dip and of the mean dip direction.
  6. Is the scatter of the dip directions, 2.9∘2.9{}^{\circ}, a sign of careless measurement? Explain.

Part III — The declination.

  1. A second geologist forgot to set the declination. Her twelve readings have the vector mean 083.1/12.1083.1/12.1 in magnetic azimuths. Suppose the magnetic declination in the Varden Hills that year is 2.5∘2.5{}^{\circ} west. What is her mean in geographic azimuths?
  2. Why is her mean dip right without correction?
  3. Could she have corrected each reading before averaging, rather than the mean? Does it change the answer (Exercise 6.11)?
  4. Ten years later, a student reuses her declination setting while the declination has changed by 0.2∘0.2{}^{\circ} a year. How large is the error, compared with the standard error of question 12?
  5. Give the strike of the mean attitude of the first geologist by the right-hand rule.

Part IV — Compass and boreholes.

  1. From three boreholes, Chapter 4 found the attitude 080.6/12.2080.6/12.2 for the base of the limestone. Compare it with the vector mean of question 7, angle by angle.
  2. Compute the angle between the two normals.
  3. The standard error of the mean direction is about θrms/N\theta_{\text{rms}}/\sqrt N. Is the difference of question 20 significant?
  4. How many readings would halve that standard error?
  5. Name two errors that no number of readings would reduce.
  6. Give the mean attitude of the east limb from the twelve readings, in the form dip direction/dip rounded to the degree.
Solution

Solution of Problem 6.1.

1. Dip directions: 967/12=80.6∘967/12 = 80.6{}^{\circ}; dips: 145/12=12.1∘145/12 = 12.1{}^{\circ}.

2. 355∘355{}^{\circ} and 005∘005{}^{\circ} average to 180∘180{}^{\circ}, the opposite direction. Here all the dip directions lie between 076∘076{}^{\circ} and 085∘085{}^{\circ}, far from the jump at 000∘000{}^{\circ}, so the arithmetic mean works.

3. ∥n⃗∥2=sin⁡2δ(sin⁡2α+cos⁡2α)+cos⁡2δ=1\|\vec n\|^2 = \sin^2\delta(\sin^2\alpha + \cos^2\alpha) + \cos^2\delta = 1. The downward dip line has the direction (cos⁡δsin⁡α,cos⁡δcos⁡α,−sin⁡δ)(\cos\delta\sin\alpha, \cos\delta\cos\alpha, -\sin\delta), and its dot product with n⃗\vec n is sin⁡δcos⁡δ−cos⁡δsin⁡δ=0\sin\delta\cos\delta - \cos\delta\sin\delta = 0; the strike line is horizontal and perpendicular to the dip direction, (cos⁡α,−sin⁡α,0)(\cos\alpha, -\sin\alpha, 0), and its dot product with n⃗\vec n is sin⁡δsin⁡αcos⁡α−sin⁡δcos⁡αsin⁡α=0\sin\delta\sin\alpha\cos\alpha - \sin\delta\cos\alpha\sin\alpha = 0.

4. Reading 11 is the same as reading 2: (0.2222,0.0352,0.9744)(0.2222, 0.0352, 0.9744). Reading 12:

(sin⁡11∘sin⁡82∘, sin⁡11∘cos⁡82∘, cos⁡11∘)=(0.1890, 0.0266, 0.9816).(\sin11{}^{\circ}\sin82{}^{\circ},\ \sin11{}^{\circ}\cos82{}^{\circ},\ \cos11{}^{\circ}) = (0.1890,\ 0.0266,\ 0.9816).

5. S⃗=(2.4747, 0.4082, 11.7294)\vec S = (2.4747,\ 0.4082,\ 11.7294).

6. ∥S⃗∥=6.1241+0.1666+137.5788=11.99457\|\vec S\| = \sqrt{6.1241 + 0.1666 + 137.5788} = 11.99457, and Rˉ=11.99457/12=0.999547\bar R = 11.99457/12 = 0.999547.

7. αˉ=atan2⁡(2.4747,0.4082)=80.6∘\bar\alpha = \operatorname{atan2}(2.4747, 0.4082) = 80.6{}^{\circ} and δˉ=arccos⁡(11.7294/11.99457)=arccos⁡0.977893=12.1∘\bar\delta = \arccos(11.7294/11.99457) = \arccos 0.977893 = 12.1{}^{\circ}: the vector mean is 080.6/12.1080.6/12.1, the same as the arithmetic means here.

8. Dip directions: s=2.94∘s = 2.94{}^{\circ}; dips: s=1.68∘s = 1.68{}^{\circ}.

9. θrms≈2×0.000453=0.0301 rad=1.72∘\theta_{\text{rms}} \approx \sqrt{2\times0.000453} = 0.0301\,\mathrm{rad} = 1.72{}^{\circ}.

10. 1.682+(sin⁡12.1∘×2.94)2=2.82+0.38=1.79∘\sqrt{1.68^2 + (\sin12.1{}^{\circ}\times2.94)^2} = \sqrt{2.82 + 0.38} = 1.79{}^{\circ}, close to 1.72∘1.72{}^{\circ}: the readings behave as the error model says, within the effects of rounding to the degree and of a small sample.

11. The dip: it contributes 2.82 of the 3.20 (in square degrees). The dip direction counts only through sin⁡δ=0.21\sin\delta = 0.21, because on a gentle bed a change of dip direction turns the plane very little.

12. Mean dip: 1.68/12=0.48∘1.68/\sqrt{12} = 0.48{}^{\circ}; mean dip direction: 2.94/12=0.85∘2.94/\sqrt{12} = 0.85{}^{\circ}.

13. No: 2.9∘2.9{}^{\circ} of dip direction on a bed dipping 12∘12{}^{\circ} is only 2.9sin⁡12.1∘=0.6∘2.9\sin12.1{}^{\circ} = 0.6{}^{\circ} of turn of the plane, less than the 1.7∘1.7{}^{\circ} of the dips. Any slight undulation of a gentle bedding surface scatters its dip direction widely.

14. D=−2.5∘D = -2.5{}^{\circ}: 83.1−2.5=080.6∘83.1 - 2.5 = 080.6{}^{\circ}; her mean is 080.6/12.1080.6/12.1, the same as the first geologist’s.

15. The dip is read with the clinometer, against the vertical given by gravity; the needle plays no part.

16. Yes: a rotation about the vertical moves every reading and the mean by the same angle (Exercise 6.11); the answer is the same.

17. 10×0.2=2∘10\times0.2 = 2{}^{\circ} of error on every dip direction, more than twice the standard error of 0.85∘0.85{}^{\circ}: a systematic error that the averaging cannot detect.

18. 80.6−90=350.6∘80.6 - 90 = 350.6{}^{\circ}; the attitude is strike 351∘351{}^{\circ}, dip 12∘12{}^{\circ} to the east.

19. Dip directions 080.6080.6 and 080.6080.6: no difference at this precision; dips 12.112.1 and 12.212.2: 0.1∘0.1{}^{\circ}.

20. The normals (sin⁡12.2∘sin⁡80.6∘,sin⁡12.2∘cos⁡80.6∘,cos⁡12.2∘)(\sin12.2{}^{\circ}\sin80.6{}^{\circ}, \sin12.2{}^{\circ}\cos80.6{}^{\circ}, \cos12.2{}^{\circ}) and S⃗/∥S⃗∥\vec S/\|\vec S\| make an angle of 0.13∘0.13{}^{\circ}.

21. 1.72/12=0.50∘1.72/\sqrt{12} = 0.50{}^{\circ}; the difference, 0.13∘0.13{}^{\circ}, is a quarter of it: not significant. The compass and the boreholes agree.

22. Four times as many: 48 readings.

23. For example: a wrong or outdated declination setting; a compass deflected by iron nearby (a vehicle, tools, magnetic rocks); a clinometer out of adjustment; readings taken only on the best exposed surfaces, which may not represent the limb. Averaging reduces random errors, not systematic ones.

24. The mean attitude of the east limb is 081/12: dip direction 081∘081{}^{\circ}, dip 12∘12{}^{\circ}.

本章定义的术语

在术语表中查看全部 957 个术语