Calculation of the distance between the earth and the Sun from measurements taken in occasion of a transit of Venus
Back
By
Dr. Carme Jordi i Nebot 
Translation by Víctor Gómez Jiménez
Astronomy and Meteorology Department of the University of Barcelona
June 2nd, 2004
1. Introduction
The main goal of this document is to describe the grounds on which the calculation of the distance Earth-Sun from the observation of a transit of Venus is founded. To describe the rigorous calculation is not intended, because it involves a level of mathematical knowledge that the most students and general public have not still acquired. Who is interested in this rigorous treatment can check the bibliography at the end of the document.
The determination of the distance Earth-Sun is based on the perspective effect by which, from two different locations, Venus is projected over two distinct points of the solar disc. Therefore, observations from distant places on Earth must be combined: the perspective effect is the more important the further this two localions are, and so a more accurate measurement of the astronomical unit can be derived.
Observations must be complemented by the laws of Kepler which describe the orbits of planets around the Sun and were deduced by Johannes Kepler from very numerous observations of the planetary motion. The law of Universal Gravitation, formulated by sir Isaac Newton, applied to the case of two bodies moving around a common centre of mass explains the three empirical laws of Kepler.
In this document, a simplified method to calculate the distance between the Earth and the Sun is described. It needs simultaneous observations of the transit from two distant locations on the surface of Earth, in order to measure the distance between the centre of Venus over the solar disc as seen from each point. Another method to carry out this calculation is founded on the comparison of the time Venus spends on crossing the solar disc from side to side (the duration of the transit) when observed from different places. The mathematical approach of both methods is similar, but the last one is quite more demanding from an observational and mathematical point of view.
2. Using simultaneous observations from different places
This method is founded on knowing the distances between the centre of Venus, projected over the solar disc, in a given instant, as seen from two different locations. Therefore, to carry out simultaneous observations from two different places on Earth (the more distant the better) is imperatively needed to use this simplified method, and the easiest way to achieve it is to take a photograph of the Sun at the same moment from both locations.
From two different places M
1 and M
2 (see Figure 1) and at the same time t, Venus is projected on two different positions V
1 and V
2 over the solar disc due to the perspective effect. Here is happening the same than when we put a finger very near our eyes and look at it closing alternatively one eye or the other. Then, we see the finger projected over the background in a different position with each eye. In the case of the transit of Venus, the two observers are equivalent to the eyes, Venus to the finger and the Sun to the background. Depending on the place from where we observe, Venus is projected over a point or another of the solar disc.
The measurement of the distance between V
1 and V
2 in terms of the solar radius makes possible to determine the average distance Earth-Sun, but not in an easy way, in spite of the simplifications we adopt here.

Figure 1. Observation of the transit of Venus from two different locations M
1 and M
2 at the same time.
Obtained from a page by P. Rocher
(IMCCE)
Let's present the mathematical formulation.
2.1 Determination of the distance in the moment of the observation
Let us suppose that O is the centre of the Earth, C the centre of the Sun and V
1 and V
2 the observed centres of the Venus' projections as seen from points M
1 and M
2,respectively, the angles D
1 and D
2, the angular separations between each centre of Venus and the centre of the Sun as seen again from M
1 and M
2, respectively; that is to say, the parallax angles CM
1V
1 and CM
2V
2. Analogously, we can define the angles π
s and π
v as the angular separations between M
1 and M
2 as seen from the Sun and from Venus, respectively; that is, the angles M
1CM
2 and M
1VM
2.
Given that the four points M
1, M
2, C and V are not in the same plane (the most usual case is that M
1 and M
2 are not on the same meridian and the Earth, Venus and the Sun not perfectly aligned), the geometry of the problem becomes a bit complicated, in such a way that, as can be seen in figure 2, the distance Δπ between the two centres of Venus is the only observable magnitude. It corresponds to Δπ = π
v – π
s and makes possible to calculate the distance to the Sun.
The carrying out of the measurement of Δπ from the two photographs can be done by finding the position of the centre of Venus with respect to a common point of reference on the solar disc (a spot, for instance), and then comparing the difference with the total size of the disc. Measurements on the photographs are taken in length units (
mm for instance) and we must transform them to angles. So, we need to know the angular size of the Sun.
Let's call (x
1,y
1) and (x
2,y
2) to the spaces, in
mm, between the centre of Venus and the reference spot in the horizontal and vertical directions for each photograph. Then, the spaces in arcseconds are obtained by multiplying each magnitude x
1 and y
1 by the factor Ã'(arcseconds)/Ã'
1(
mm) and each magnitude x
2 and y
2 by Ã'(arcseconds)/Ã'
2(
mm) where Ã'(arcseconds) and Ã'(
mm) are the diameter of the Sun expressed in arcseconds and
mm, respectively. Ã'
1(
mm) and Ã'
2(
mm) have the same value if both photographs share the scale.
The distance between the centres of Venus in the two photographs is:
Δπ(segons d'arc) = [ (x
2 – x
1)
2+ (y
2 – y
1)
2 ]
1/2 [1]
If the photographs are taken with telescopes that give the same scale and the solar disc is situated exactly in the same place of the photographs, then a corner of the photograph can be taken as a reference point. Besides, in this case Ã'
1(
mm) = Ã'
2(
mm) and the previous equation can be written as:
Δπ(segons d'arc) = [ (x
2 – x
1)
2+ (y
2 – y
1)
2 ]
1/2·Ã'(segons d'arc)/Ã'(
mm) [1bis]
This is the supposition adopted from now on.

Figure 2. Positions of the Venus' projections over the solar disc.
Obtained from a page by P. Rocher
(IMCCE)
Let's suppose that r
V and r
T are the distances between the centre of the Sun and the centres of Venus and the Earth, respectively, at the moment t of the observation. As the projection d of the distance between M
1 and M
2 over the plane perpendicular to OC is small if compared with the distances Earth-Sun and Earth-Venus, we can approximate:
π
s = d/r
T π
v = d/(r
T -r
V)
and from this, the following is deduced:
π
v = π
s r
T/(r
T -r
V)
Δπ = π
s (r
T/(r
T -r
V)-1) = π
s r
v/(r
T -r
V)
and so,
π
s = d/r
T = Δπ (r
T/r
V - 1)
This last formula clearly express that, if we know the angular distance Δπ between the centres V
1 and V
2, and the ratio r
T/r
V between the distances Earth-Sun and Venus-Sun, the parallax π
s can be deduced; and that if we know the projected distance d between the two locations where observations are carried out, the distance r
T can be calculated. (In these expressions, the values of π
v, π
s and Δπ are given in radians. To transform them to arcseconds to make them compatible with equation [1], multiplying by number 64800/π is enough).
Δπ is the observable magnitude, d can be determined as is explained beneath and then, the only magnitude we have to know to solve the problem is the ratio r
T/r
V between the distances Earth-Sun and Venus-Sun.
The orbits of the Earth and Venus around the Sun are slightly elliptical and so, the ratio r
T/r
V is not constant in time. To know this ratio at the exact time t, we need to use the first law of Kepler, which states that the Sun is situated on one of the ellipse's focuses. Then, the distance to the planet is given by:
r
p(t)=R
p (1 - e
p cos E
p(t))
where R
p is the semimajor axis of the orbit, e
p the eccentricity and E
p(t) the eccentric anomaly at the instant t. According to this:
r
T/r
v=[R
T (1 - e
T cos E
T)] / [R
V (1 - e
V cos E
V)]
The third law of Kepler gives a relation between the semimajor axes of the orbits and its correspounding periods P
p:
(R
T / R
V)
3 = (P
T / P
V)
2
in such a way that
r
T/r
v=(P
T / P
V)
2/3 (1 - e
T cos E
T) / (1 - e
V cos E
V) [2]
Until this point, we have been able to determine π
s and r
T; the parallax and the distance Earth-Sun at the instant t of the observation.
2.2 Determination of the average distance
To determine the average distance Earth-Sun (R
T) and the correspounding average parallax π
o, which are related by means of the equatorial radius R according to:
π
o ≈≈ R/R
T
some other considerations must be done.
If the projection d of the distance between M
1 and M
2 on the plane normal to the direction Earth-Sun is expressed in terms of the equatorial terrestrial radius, and the distance Earth-Sun in terms of the average distance, we have:
π
s = [(d/R) / (r
T /R
T)] (R/R
T) ≈ [(d/R) / (r
T /R
T)] π
o
The ratio r
T/R
T can be deduced from the first law of Kepler as:
r
T/R
T = 1 - e
T cos E
T(t)
and so, we only have to calculate d/R (see Figure 3).


Figure 3. Projection of the distance between M
1 and M
2 on the plane normal to the direction Earth-Sun.
Obtained from a page by P. Rocher
(IMCCE)
Calculating the vectorial product between vectors M
1M
2 and OC, the value of sin θ is obtained, due to
M
1M
2 × OC = |M
1M
2| r
T sin θ
The next expression is clear in figure 3:
d = |M
1M
2| cos (90 – θ) = |M
1M
2| sin θ
and then,
d = M
1M
2 × OC / r
T
Now we just have to find M
1M
2 × OC.
Calculation of the vector OC
This vector can be expressed in terms of the equatorial coordinates of the Sun (α,δ) at the moment of the observation as:
x=r
T cos δ cos α
y=r
T cos δ sin α
z=r
T sin δ
Calculation of the vector M1M2
The position of each observer can be expressed as (see figure 4):
x=R cos φ cos (λ+T
G)
y=R cos φ sin (λ+T
G)
z=R sin φ
where φ and λ are the geographical coordinates (latitude and longitude) and T
G=T
G(0) + 1.00273791 t. The time t of the observation must be expressed in the scale of universal time (
UT). For the most countries in Europe
UT = official time – 2
h in June.


Figure 4. Positions of a celestial object (for instance, the Sun) and of an observer on Earth in equatorial coordinates.
Obtained from a page by P. Rocher
(IMCCE)
The coordinates of the vector M
1M
2 can be easily found as:
X=x
1 – x
2 Y=y
1 – y
2 Z=z
1 – z
2
Practical example
One of the programs available on
https://gaia.am.ub.es/Twiki/bin/view/ServiAstro/AgendaAstronomica2004 is founded on this formulation and can be used as a practical example. The second program is based on the comparison of the duration of the transit as observed from two different locations.