More topic actionsEdit   Attach

Calculation of the distance between the earth and the Sun from measurements taken in occasion of a transit of Venus


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 M1 and M2 (see Figure 1) and at the same time t, Venus is projected on two different positions V1 and V2 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 V1 and V2 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.

Projeccions de Venus sobre el Sol
Figure 1. Observation of the transit of Venus from two different locations M1 and M2 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 V1 and V2 the observed centres of the Venus' projections as seen from points M1 and M2,respectively, the angles D1 and D2, the angular separations between each centre of Venus and the centre of the Sun as seen again from M1 and M2, respectively; that is to say, the parallax angles CM1V1 and CM2V2. Analogously, we can define the angles πs and πv as the angular separations between M1 and M2 as seen from the Sun and from Venus, respectively; that is, the angles M1CM2 and M1VM2.

Given that the four points M1, M2, C and V are not in the same plane (the most usual case is that M1 and M2 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 (x1,y1) and (x2,y2) 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 x1 and y1 by the factor Ã'(arcseconds)/Ã'1(mm) and each magnitude x2 and y2 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) = [ (x2 – x1)2+ (y2 – y1)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) = [ (x2 – x1)2+ (y2 – y1)2 ]1/2·Ã'(segons d'arc)/Ã'(mm) [1bis]

This is the supposition adopted from now on.

posicions de les projeccions de Venus
Figure 2. Positions of the Venus' projections over the solar disc.
Obtained from a page by P. Rocher (IMCCE)

Let's suppose that rV and rT 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 M1 and M2 over the plane perpendicular to OC is small if compared with the distances Earth-Sun and Earth-Venus, we can approximate:

πs = d/rT
πv = d/(rT -rV)

and from this, the following is deduced:

πv = πs rT/(rT -rV)
Δπ = πs (rT/(rT -rV)-1) = πs rv/(rT -rV)

and so,

πs = d/rT = Δπ (rT/rV - 1)

This last formula clearly express that, if we know the angular distance Δπ between the centres V1 and V2, and the ratio rT/rV 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 rT 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 rT/rV 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 rT/rV 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:

rp(t)=Rp (1 - ep cos Ep(t))

where Rp is the semimajor axis of the orbit, ep the eccentricity and Ep(t) the eccentric anomaly at the instant t. According to this:

rT/rv=[RT (1 - eT cos ET)] / [RV (1 - eV cos EV)]

The third law of Kepler gives a relation between the semimajor axes of the orbits and its correspounding periods Pp:

(RT / RV)3 = (PT / PV)2

in such a way that

rT/rv=(PT / PV)2/3 (1 - eT cos ET) / (1 - eV cos EV) [2]

Until this point, we have been able to determine πs and rT; 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 (RT) and the correspounding average parallax πo, which are related by means of the equatorial radius R according to:

πo ≈≈ R/RT

some other considerations must be done.

If the projection d of the distance between M1 and M2 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) / (rT /RT)] (R/RT) ≈ [(d/R) / (rT /RT)] πo

The ratio rT/RT can be deduced from the first law of Kepler as:

rT/RT = 1 - eT cos ET(t)

and so, we only have to calculate d/R (see Figure 3).

Projeccions en el pla normal a la direcció Terra-Sol Figure 3. Projection of the distance between M1 and M2 on the plane normal to the direction Earth-Sun.
Obtained from a page by P. Rocher (IMCCE)

Calculating the vectorial product between vectors M1M2 and OC, the value of sin θ is obtained, due to

M1M2 × OC = |M1M2| rT sin θ

The next expression is clear in figure 3:

d = |M1M2| cos (90 – θ) = |M1M2| sin θ

and then,

d = M1M2 × OC / rT

Now we just have to find M1M2 × 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=rT cos δ cos α
y=rT cos δ sin α
z=rT sin δ

Calculation of the vector M1M2

The position of each observer can be expressed as (see figure 4):

x=R cos φ cos (λ+TG)
y=R cos φ sin (λ+TG)
z=R sin φ

where φ and λ are the geographical coordinates (latitude and longitude) and TG=TG(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 – 2h in June.

Diagrama explicatiu de les coordenades equatorials 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 M1M2 can be easily found as:

X=x1 – x2
Y=y1 – y2
Z=z1 – z2

Practical example

One of the programs available on 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.

3. For more information

This site is powered by the TWiki collaboration platform Powered by Perl