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 (∼65\sim 65 minutes for Q1, ∼70\sim 70 minutes for Q2 and ∼45\sim 45 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
ConstantSymbolValue
Speed of lightcc3.00×108 m s−13.00 \times 10^{8}\ \mathrm{m\,s^{-1}}
Earth’s rotation period1 day24 hours
Earth’s orbital period1 year365.25 days
parsecpc3.09×1016 m3.09 \times 10^{16}\ \mathrm{m}
Astronomical Unitau1.50×1011 m1.50 \times 10^{11}\ \mathrm{m}
Radius of the SunR⊙R_{\odot}6.96×108 m6.96 \times 10^{8}\ \mathrm{m}
Radius of the EarthR⊕R_{\oplus}6.37×106 m6.37 \times 10^{6}\ \mathrm{m}
Mass of the SunM⊙M_{\odot}1.99×1030 kg1.99 \times 10^{30}\ \mathrm{kg}
Mass of the EarthM⊕M_{\oplus}5.97×1024 kg5.97 \times 10^{24}\ \mathrm{kg}
Luminosity of the SunL⊙L_{\odot}3.85×1026 W3.85 \times 10^{26}\ \mathrm{W}
Stephan-Boltzmann constantσ\sigma5.67×10−8 J m−2 K−45.67 \times 10^{-8}\ \mathrm{J\,m^{-2}\,K^{-4}}
Gravitational constantGG6.67×10−11 m3 kg−1 s−26.67 \times 10^{-11}\ \mathrm{m^{3}\,kg^{-1}\,s^{-2}}
Boltzmann constantkBk_{\mathrm{B}}1.38×10−23 J K−11.38 \times 10^{-23}\ \mathrm{J\,K^{-1}}
Permittivity of free spaceε0\varepsilon_08.85×10−12 F m−18.85 \times 10^{-12}\ \mathrm{F\,m^{-1}}
Permeability of free spaceμ0\mu_04π×10−7 H m−14\pi \times 10^{-7}\ \mathrm{H\,m^{-1}}
Planck’s constanthh6.63×10−34 J s6.63 \times 10^{-34}\ \mathrm{J\,s}
Elementary chargeee1.60×10−19 C1.60 \times 10^{-19}\ \mathrm{C}
Proton rest massmpm_{\mathrm{p}}1.67×10−27 kg1.67 \times 10^{-27}\ \mathrm{kg}
Electron rest massmem_{\mathrm{e}}9.11×10−31 kg9.11 \times 10^{-31}\ \mathrm{kg}
Wien’s displacement lawλmax⁡T\lambda_{\max}T2.90×10−3 m K2.90 \times 10^{-3}\ \mathrm{m\,K}
Avagadro’s constantNAN_A6.02×1023 mol−16.02 \times 10^{23}\ \mathrm{mol^{-1}}

Basic calculus formulae:

Chain ruleddxf(g(x))=f′(g(x))g′(x)\frac{\mathrm{d}}{\mathrm{d}x}f(g(x)) = f'(g(x))g'(x)
Product ruleddx(uv)=dudxv+udvdx\frac{\mathrm{d}}{\mathrm{d}x}(uv) = \frac{\mathrm{d}u}{\mathrm{d}x}v + u\frac{\mathrm{d}v}{\mathrm{d}x}
Quotient ruleddx(uv)=dudxv−udvdxv2\frac{\mathrm{d}}{\mathrm{d}x}\left(\frac{u}{v}\right) = \frac{\frac{\mathrm{d}u}{\mathrm{d}x}v - u\frac{\mathrm{d}v}{\mathrm{d}x}}{v^2}
Integration by parts∫udvdx dx=uv−∫vdudx dx\int u\frac{\mathrm{d}v}{\mathrm{d}x}\,\mathrm{d}x = uv - \int v\frac{\mathrm{d}u}{\mathrm{d}x}\,\mathrm{d}x
Standard integral∫1x dx=ln⁡∣x∣+C\int \frac{1}{x}\,\mathrm{d}x = \ln\vert x\vert + C
Important Formulae

You might find the diagram of an elliptical orbit below useful in solving some of the questions:

Elements of an elliptic orbit:

a=OA (=OP)a = \mathbf{OA}\ (= \mathbf{OP})semi-major axis
b=OB (=OC)b = \mathbf{OB}\ (= \mathbf{OC})semi-minor axis
e=1−b2a2e = \sqrt{1 - \frac{b^2}{a^2}}eccentricity
F\mathbf{F}focus
PF=a(1−e)\mathbf{PF} = a(1-e)periapsis distance (shortest distance from F\mathbf{F})
AF=a(1+e)\mathbf{AF} = a(1+e)apoapsis distance (longest distance from F\mathbf{F})

Kepler’s Third Law: For an elliptical orbit, the square of the period, TT, of an object about the focus is proportional to the cube of the semi-major axis, aa (as defined above), such that

T2=4π2GMa3 ,T^2 = \frac{4\pi^2}{GM}a^3\,,

where MM is the total mass of the system (typically dominated by the central object) and GG is the universal gravitational constant.

Vis-Viva Equation: For an elliptical orbit, the speed vv of an object at a distance rr from the focus is related to the semi-major axis, aa, total mass of the system, MM, and universal gravitational constant, GG, (as defined above), such that

v2=GM(2r−1a) .v^2 = GM\left(\frac{2}{r} - \frac{1}{a}\right)\,.

Magnitudes: The apparent magnitudes of two objects, m1m_1 and m0m_0, are related to their apparent brightnesses, b1b_1 and b0b_0, via the formula

b1b0=10−0.4(m1−m0) .\frac{b_1}{b_0} = 10^{-0.4(m_1 - m_0)}\,.

The absolute magnitude of an object, M\mathcal{M}, is the same as its apparent magnitude when viewed from 10 pc, hence the relationship between apparent and absolute magnitude and distance is

m−M=5log⁡(d10),m - \mathcal{M} = 5\log\left(\frac{d}{10}\right),

where dd 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

Photograph of the Callanish standing stones on grass under a blue sky, with the low winter Sun glowing behind them; above the stones a vertical figure-of-eight of white dots (the Sun's noon positions) forms the analemma, with a small upper loop and a large lower loop. Watermark 'gmrphotographer' in the lower right corner.
Figure 1: The analemma above was composed from images taken every few days at noon near the village of Callanish in the Outer Hebrides in Scotland. In the foreground are the Callanish Stones and the main photo was taken on the winter solstice (when the maximum angle the Sun reaches above the horizon is the lowest of the year, so is at the bottom of the analemma). Credit: Giuseppe Petricca.

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, δ\delta, 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 α\alpha 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 αtilt\alpha_{\mathrm{tilt}} is the contribution due to the axial tilt and αecc\alpha_{\mathrm{ecc}} is the contribution due to the Earth’s orbital eccentricity, then

α=αtilt+αecc .\alpha = \alpha_{\mathrm{tilt}} + \alpha_{\mathrm{ecc}}\,.

If the angle of the axial tilt is ε\varepsilon and the eccentricity of the Earth’s orbit is ee, and we assume that both are small enough that the sinusoidal approximation of δ\delta, αtilt\alpha_{\mathrm{tilt}}, and αecc\alpha_{\mathrm{ecc}} apply, then we find the following boundary conditions:

  • δ\delta has a period of 1 year, an amplitude of ε\varepsilon, is maximum at the summer solstice (21st June) and minimum at the winter solstice (21st December)
  • αtilt\alpha_{\mathrm{tilt}} has a period of 0.5 years, an amplitude (in radians) of tan⁡2(ε/2)\tan^2(\varepsilon/2), 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
  • αecc\alpha_{\mathrm{ecc}} has a period of 1 year, an amplitude (in radians) of 2e2e, is zero at the perihelion (4th January) and the aphelion (6th July), and (using our sign convention) negative just after the perihelion

Given the nthn^{\mathrm{th}} day of the year, a value can be calculated for δ\delta and α\alpha, and these are the co-ordinates for the analemma (it is drawn by these parametric equations). For the Earth, ε=23.44∘\varepsilon = 23.44^{\circ} and e=0.0167e = 0.0167.

a. (i) Although α\alpha is really an angle in radians (where 2π2\pi radians =360∘= 360^{\circ}), 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 αtilt\alpha_{\mathrm{tilt}} and αecc\alpha_{\mathrm{ecc}} for the Earth into minutes. [3 т.]

a. (ii) Determine equations for δ\delta (in degrees) and αtilt\alpha_{\mathrm{tilt}} and αecc\alpha_{\mathrm{ecc}} (both in minutes) as a function of the day of the year, nn. Take n=1n = 1 to be 1st January and n=365n = 365 to be 31st December. [3 т.]

a. (iii) Sketch the analemma you would see if the Earth’s orbit was circular (i.e. e=0e = 0) but with the current value of ε\varepsilon, with α\alpha on the xx-axis (in minutes) and δ\delta on the yy-axis. [3 т.]

b. (i) Using your equations derived in the previous part, sketch on GRAPH PAPER on the same set of axes:

αecc\alpha_{\mathrm{ecc}} (vertical axis) in minutes against nn (horizontal axis). [1 т.]

b. (ii) αtilt\alpha_{\mathrm{tilt}} in minutes against nn. [1 т.]

b. (iii) Hence, carry out the superposition of those two waves to get α\alpha against nn. [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 α\alpha? [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: α≠0\alpha \neq 0 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.

Yellow closed curve on a grey background with black horizontal and vertical axes. The vertical axis has ticks labelled +15 degrees (top) and -15 degrees (bottom); the horizontal axis has ticks labelled -0h26m47s (left) and 0h18m56s (right). Yellow circles with month letters lie on the curve: going round from the top, F, J, D, N, O, S, A, J, J, M, A, M. The loop reaches the top near F, the far left near A, the far right near O and the bottom near A and S.
Figure 2: The analemma of the Sun at midday as seen by an observer on Earth 2.0. In this situation, α\alpha ranges from −26-26 mins 4747 secs to 1818 mins 5656 secs. The circled letters correspond to the same (unknown) day of each month (for example 5th Jan, 5th Feb, 5th March etc.). Credit: Bob Urschel.

d. (i) Considering data from Figure 2, determine for Earth 2.0:

The axial tilt, ε\varepsilon. [1 т.]

d. (ii) The month of the vernal equinox. [1 т.]

d. (iii) The eccentricity of the orbit, ee. [5 т.]

Отговори

Решение

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