BAAO 2023, Теория — Задача 1. Solar Analemma
British Astronomy and Astrophysics Olympiad 2022-2023 — Astronomy & Astrophysics Competition Paper — Monday 6th February 2023 · 6 февруари 2023 г. · 35 т.
Указания
This question paper must not be photographed or taken out of the exam room
Time: 3 hours ( minutes for Q1, minutes for Q2 and minutes for Q3).
Questions: All three questions should be attempted. Each question contains independent parts so that later parts can be attempted even if earlier parts are incomplete.
Solutions: Answers and calculations are to be written on loose paper. YOU WILL NEED GRAPH PAPER. Students should ensure their name and school is clearly written on the first answer sheet and that all pages are numbered. EACH QUESTION ANSWERED must be started on a new page. A standard formula booklet with standard physical constants may be used if desired.
Instructions: To accommodate students sitting the paper at different times, please do not discuss any aspect of the paper on the internet until 8 am Saturday 11th February.
Clarity: Solutions must be written legibly, in black pen, and working down the page. Scribble will not be marked and overall clarity is an important aspect of this exam.
Eligibility: The International Olympiad will be held in August 2023; all sixth form students are eligible to participate.
Calculators: Any standard calculator may be used, but calculators cannot be programmable and must not have symbolic algebra capability.
Training Dates and the IOAA (Chorzów, Poland, 10th - 20th August 2023)
The team will be selected from students taking this paper. The best students that are eligible to represent the UK at the IOAA will be invited to attend the Training Camp to be held in Oxford from Saturday 1st to Wednesday 5th April 2023. Astronomy material will be covered; problem solving skills and observational skills (telescope and naked eye observations) will be developed. At the Training Camp a Data Analysis exam along with a Round 3 theory paper will be sat. A team of five students (plus one reserve) will be selected for further training, including additional training camps in the summer.
Important Constants
| Constant | Symbol | Value |
|---|---|---|
| Speed of light | ||
| Earth’s rotation period | 1 day | 24 hours |
| Earth’s orbital period | 1 year | 365.25 days |
| parsec | pc | |
| Astronomical Unit | au | |
| Radius of the Sun | ||
| Radius of the Earth | ||
| Mass of the Sun | ||
| Mass of the Earth | ||
| Luminosity of the Sun | ||
| Stephan-Boltzmann constant | ||
| Gravitational constant | ||
| Boltzmann constant | ||
| Permittivity of free space | ||
| Permeability of free space | ||
| Planck’s constant | ||
| Elementary charge | ||
| Proton rest mass | ||
| Electron rest mass | ||
| Wien’s displacement law | ||
| Avagadro’s constant |
Basic calculus formulae:
| Chain rule | |
| Product rule | |
| Quotient rule | |
| Integration by parts | |
| Standard integral |
Important Formulae
You might find the diagram of an elliptical orbit below useful in solving some of the questions:
Elements of an elliptic orbit:
| semi-major axis | |
| semi-minor axis | |
| eccentricity | |
| focus | |
| periapsis distance (shortest distance from ) | |
| apoapsis distance (longest distance from ) |
Kepler’s Third Law: For an elliptical orbit, the square of the period, , of an object about the focus is proportional to the cube of the semi-major axis, (as defined above), such that
where is the total mass of the system (typically dominated by the central object) and is the universal gravitational constant.
Vis-Viva Equation: For an elliptical orbit, the speed of an object at a distance from the focus is related to the semi-major axis, , total mass of the system, , and universal gravitational constant, , (as defined above), such that
Magnitudes: The apparent magnitudes of two objects, and , are related to their apparent brightnesses, and , via the formula
The absolute magnitude of an object, , is the same as its apparent magnitude when viewed from 10 pc, hence the relationship between apparent and absolute magnitude and distance is
where is measured in parsecs.
Note for markers
BAAO 2022/23 Solutions and Marking Guidelines
Note for markers:
- Answers to two or three significant figures are generally acceptable. The solution may give more in order to make the calculation clear. Units should be present on final answers when appropriate.
- There are multiple ways to solve some of the questions; please accept all good solutions that arrive at the correct answer. Students getting the answer in a box will get all the marks available for that calculation / part of the question (as indicated in red), so long as there are no unphysical / nonsensical steps or assumptions made (students may not explicitly calculate the intermediate stages and should not be penalised for this so long as their argument is clear).
Условие
Plotting the position of the Sun in the sky at the same time every day, you get an interesting figure-of-eight shape known as an analemma (see Figure 1). For observers in the Northern hemisphere, you might expect to always see the Sun due South at midday, however on some days the Sun has already passed through that bearing and on others it needs a few more minutes before it gets there. This is due to two effects: the axial tilt of the Earth, and the fact the Earth’s orbit is not perfectly circular

The vertical co-ordinate of a point in the analemma is entirely determined by the Earth’s axial tilt. This is known as the solar declination, , and varies sinusoidally throughout the year. The horizontal co-ordinate of a point in the analemma is determined by a combination of the Earth’s axial tilt and the eccentricity of the Earth’s orbit. Both of these individually vary sinusoidally, but the superposition of the two is no longer sinusoidal.
We will define as the angle between due South and the Sun at local midday as seen from Oxford, where a positive value means the Sun has already passed through due South (so is on the right of the figure above) whilst a negative value means the Sun has yet to pass through due South. If is the contribution due to the axial tilt and is the contribution due to the Earth’s orbital eccentricity, then
If the angle of the axial tilt is and the eccentricity of the Earth’s orbit is , and we assume that both are small enough that the sinusoidal approximation of , , and apply, then we find the following boundary conditions:
- has a period of 1 year, an amplitude of , is maximum at the summer solstice (21st June) and minimum at the winter solstice (21st December)
- has a period of 0.5 years, an amplitude (in radians) of , is zero at the solstices and the equinoxes (vernal equinox = 21st March, autumnal equinox = 21st September), and (using our sign convention) positive just after the vernal equinox
- has a period of 1 year, an amplitude (in radians) of , is zero at the perihelion (4th January) and the aphelion (6th July), and (using our sign convention) negative just after the perihelion
Given the day of the year, a value can be calculated for and , and these are the co-ordinates for the analemma (it is drawn by these parametric equations). For the Earth, and .
a. (i) Although is really an angle in radians (where radians ), it is normally more useful to convert it into time units (essentially the time since the Sun was due South, or the time until the Sun reaches due South). Taking the mean solar day to be exactly 24 hours:
Convert the amplitude of and for the Earth into minutes. [3 т.]
a. (ii) Determine equations for (in degrees) and and (both in minutes) as a function of the day of the year, . Take to be 1st January and to be 31st December. [3 т.]
a. (iii) Sketch the analemma you would see if the Earth’s orbit was circular (i.e. ) but with the current value of , with on the -axis (in minutes) and on the -axis. [3 т.]
b. (i) Using your equations derived in the previous part, sketch on GRAPH PAPER on the same set of axes:
(vertical axis) in minutes against (horizontal axis). [1 т.]
b. (ii) in minutes against . [1 т.]
b. (iii) Hence, carry out the superposition of those two waves to get against . [Note: a clear indication of the general shape will do - no need to spend very long on making this perfect] [2 т.]
c. (i) Using your graph to guide you to the relevant point of the year, but using your precise algebraic expressions (rather than reading off the graph):
What are the dates and values of the maximum and minimum values of ? [5 т.]
c. (ii) What are the dates when the Sun is due South at midday? [5 т.]
c. (iii) What are the dates when the Sun passes between the large and small loops in the figure-of-eight (i.e. the crossover point)? [Note: at that point, but is small] [5 т.]
Consider an alternative version of Earth, known as Earth 2.0. On this planet, the year is unchanged and the perihelion and aphelion are at the same time, but it has a different axial tilt, a different orbital eccentricity, and a different month for the vernal equinox (although it is still on the 21st day of that month). The analemma as viewed from Earth 2.0 is show in Figure 2 below.

d. (i) Considering data from Figure 2, determine for Earth 2.0:
The axial tilt, . [1 т.]
d. (ii) The month of the vernal equinox. [1 т.]
d. (iii) The eccentricity of the orbit, . [5 т.]
Отговори
a. (i) αtilt amplitude 9.86 minutes (= 9 m 51.8 s); αecc amplitude 7.65 minutes (= 7 m 39.3 s)
a. (ii)
c. (i) Minimum: 13th Feb (n = 44), α = −14.31 mins; maximum: 31st Oct (n = 304), α = +16.57 mins
c. (ii) 16th Apr, 13th Jun, 31st Aug and 25th Dec (n = 106, 164, 243, 359)
c. (iii) 14th Apr (n = 104) and 29th Aug (n = 241)
d. (i) 15 °
d. (ii) November
d. (iii) 0.0498
Using the other extremum, or the shorter method described in the official solution, gives e = 0.0499.
Решение
Q1 – Solar Analemma [35 marks]
a. (i)
To convert from radians to minutes of time need to multiply by a factor of:
[1]
Putting the numbers into the given formula:
[1]
[1] [3]
[If they only calculate both in radians, award 1 mark. If they confuse the needed conversion with arcminutes (giving = 148’ and = 115’), lose first mark but give ecf. Conversion done correctly to just one means max 2 marks. Allow ecf for later parts.]
a. (ii)
Considering the number of days between the winter solstice (21st Dec) and 1st Jan, which is a minimum, then
[1] [1]
Since the winter solstice is a zero point, and is negative after it, with double the frequency of :
[1] [1]
Since perihelion on 4th Jan is a zero point and is negative after it, with the same frequency as
[1] [1]
[Allow any equivalent expressions in sin or cos, or using degrees inside the trigonometric functions. In these equations we are assuming midday on 1st Jan has n = 1 (so n = 0 & n = 365 at midday on 31st Dec) – do not penalise students for choosing a different zero point (e.g. midnight on 1st Jan) and hence getting slightly different additive terms to their n (allow ±1)]
a. (iii)
Correct shape, with rotational and reflectional symmetry and passing through (0,0) [1]
axis labelled and extrema labelled as -9.86 mins and 9.86 mins (allow ecf from a)i)) [1]
axis labelled and extrema labelled as -23.44° and 23.44° [1]

[3]
b.

Part i) is the orange line [1] [1]
Part ii) is the blue line [1] [1]
Part iii) is the grey line
Has two maxima / minima (one bigger than the other) in ≈ the right place [1]
Largest maxima is larger amplitude than the lowest minima [1] [2]
[Allow ecf from earlier parts. Lose 0.5 marks for missing axes labels]
c. (i)
We can use the differential of to find turning points, and then we know from the graph we want the ones near n ≈ 50 and n ≈ 300
[1]
Looking where d/dn = 0 (using any suitable method e.g. binary search) we find turning points in at
[0.5 marks for each] [2]
We want the first and last values:
and [1]
and [1] [5]
[Do not penalise non-integer n. Allow n and dates consistent with their counting system (so may be different by ±1). For the final two marks it is 0.5 marks for a correct date and 0.5 marks for a correct ]
[Our equation for alpha is known as the Equation of Time, and positive values mean that the Sun is ahead of the clock (so to the right of due south in the photo). Despite our approximations for the two components of , we are relatively close to the real extrema: -14.25 mins on 11th Feb and +16.42 mins on 3rd Nov]
c. (ii)
These correspond to when = 0 [1]
Using any suitable method (e.g. binary search)
[0.5 marks for each] [2]
These correspond to the following dates:
and [0.5 marks for each] [2] [5]
[Do not penalise non-integer n. Allow n and dates consistent with their counting system (so may be different by ±1)]
[The real values are 15th Apr, 13th Jun, 1st Sep, and 25th Dec, so our approximations have done very well]
c. (iii)
We are told that is small and can see from the photo that it is not near the solstices, so must be near the zero points in far from the solstices i.e. 16th Apr and 31st Aug [1]
On 16th Apr is increasing, and on 31st Aug is decreasing. The know that at the crossover is equal on both dates (due to the equation for being double valued in the range of n) and between their two values of therefore the relevant dates must be earlier [1]
Using a suitable method (e.g. starting a binary search with these points as initial values)
[2]
so date and n = 241 so date [1] [5]
[A wide variety of numerical methods might be employed here, although binary search close to the roots of is probably the quickest. Allow ecf considerations as in earlier steps. Give 0.5 marks for each correct date]
[The real dates are 12th Apr and 29th Aug, so we have again done rather well]
d. (i)
The extrema in will be so just reading off the graph [1] [1]
d. (ii)
At the vernal equinox, = 0 and is increasing, so it must be [1] [1]
[Allow 0.5 marks for May, if it is clear it was chosen due to requiring = 0]
d. (iii)
With the information in the previous two parts, we can create an equation for
The component due to eccentricity has the same general form
[0.5]
Using the earlier equations for
[0.5]
Since the vernal equinox is now 21st November (30 days earlier than the winter solstice used to be), our sinusoidal equation can be modified to
[Accept equivalent forms] [0.5]
Putting in data for one extremum (e.g. -26 mins 47 secs = -26.783 mins) and knowing it is a turning point (so ) gives us two equations:
①
② [0.5]
Combining ① and ② to eliminate e, gives the equation
[1]
Using an appropriate numerical method to solve it for , the only valid values in the allowed region are
or [1]
Recognising that the relevant extremum is around April, = 96 is the correct one to choose.
Putting this into either ① or ② leads to [1] [5]
[Using the other extremum gets or and so recognising the latter is the correct one to use gives . The method outlined here works in general, however in this specific case an alternative (much shorter) method will work; has minima at n = 97 and 279, which almost exactly coincide with the turning points of at n = 95 (min) and 278 (max). Consequently, by considering how the graph of , and would look, it is clear that and . Taking the average of the modulus of the extrema in , and assuming that is 2e gives . Students that use this approach to take advantage of the greatly simplifying phase relationship between and receive full credit (same first 1.5 marks, 1 mark for turning points, 0.5 marks for recognising the consequences of the tuning points being close together, 1 mark for setting up equations, 1 mark for final answer). Those whose working implies they did not recognise this was a special case (and hence suggest this method would always work) get 1 mark only for a correct final answer]
The relevant graph of superposition is below, emphasising why in this case the simple method would work due to the alignment of the minima of with the max / min of .

Questions proposed by
END OF PAPER
Questions proposed by:
Dr Alex Calverley (Surbiton High School)
Josh Brown (University of Cambridge)
Оригинал в Архива: BAAO_2023.pdf · официални решения: BAAO_2023_solutions.pdf