Next Article in Journal
Hybrid-Scale Hierarchical Transformer for Remote Sensing Image Super-Resolution
Previous Article in Journal
Monitoring of Land Cover and Vegetation Changes in Juhugeng Coal Mining Area Based on Multi-Source Remote Sensing Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ground-Based Microwave Measurements of Mesospheric Ozone Variations over Moscow Region during the Solar Eclipses of 20 March 2015 and 25 October 2022

by
Sergey Rozanov
1,2,*,
Alexander Ignatyev
3 and
Alexey Zavgorodniy
2
1
P. N. Lebedev Physical Institute of the Russian Academy of Sciences, Leninskiy Prospekt 53, 119991 Moscow, Russia
2
Russian Metrological Institute of Technical Physics and Radio Engineering VNIIFTRI, Mendeleevo, Solnechnogorskiy District, 141570 Moscow, Russia
3
Federal State Unitary Enterprise “NII Kvant”, 4th Likhachevskiy Pereulok 15, 125438 Moscow, Russia
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(13), 3440; https://doi.org/10.3390/rs15133440
Submission received: 31 May 2023 / Revised: 30 June 2023 / Accepted: 4 July 2023 / Published: 7 July 2023
(This article belongs to the Section Atmospheric Remote Sensing)

Abstract

:
An increase in the ozone content in the mesosphere over the Moscow region during the solar eclipses of 20 March 2015 and 25 October 2022 was observed by means of a ground-based microwave radiometer operated at frequencies of the ozone spectral line of 142.175 GHz. Changes in ozone mixing ratio (OMR) at altitudes of 90 km and 65 km were estimated and compared with diurnal ozone variations measured on the dates closest to the events. It was found that the observed increase in the OMR at 90 km during the 20 March 2015 eclipse was almost two times greater than during the 25 October 2022 eclipse, although the maximum Sun’s obscurations of these eclipses were close to each other (0.625 and 0.646). Most likely, this difference can be explained by the difference in concentration of atomic hydrogen, which plays an important role in ozone destruction at altitudes of around 90 km and above.

Graphical Abstract

1. Introduction

Atmospheric ozone O3 is one of the most important minor gases in the atmosphere. Stratospheric ozone protects life on Earth from harmful UV-B solar radiation. The absorption of radiation by ozone determines the temperature profile of the stratosphere and influences its dynamics. Production and loss of ozone in the middle atmosphere in main features can be explained by oxygen chemistry [1]. In fact, ozone loss is controlled not only by photodissociation and reaction between O3 and atomic oxygen O. The catalytic reaction of ozone with radicals of odd hydrogen HOx (H, OH, and HO2), nitrogen NOx (N, NO, and NO2), and chlorine Clx (Cl and ClO) groups are very important as well [1,2]. In the mesosphere (altitudes A = 50–90 km), ozone destruction in reactions with the HOx radicals is usually dominant. Concentrations of the NOx and HOx components and rates of catalytic destruction of ozone by the species considerably increase under the ionization of the mesosphere in polar and subpolar latitudes due to energetic charged particles (electrons and protons) [3,4].
Solar eclipses are relatively short events, about 2 h at local ground-based observations. If the eclipse occurs at midday, both the solar zenith angle (SZA) and solar spectrum in the mesosphere change little during the eclipse. Thus, solar eclipses provide a unique possibility for experimentally studying the photochemistry of mesospheric ozone and testing photochemical models. The photochemical lifetime of ozone in the mesosphere is small, about 100 s [1,2], and ozone content follows the intensity of solar radiation with negligible delay. This case is much simpler for analysis than the diurnal ozone variations during sunrises and sunsets when significant changes in the optical depth of the atmosphere at SZA close to 90°, drastically changing the intensity and spectrum of solar radiation.
Up to now, only a few observations of mesospheric ozone during solar eclipses have been conducted using rockets [5], in near-IR [6], and at millimeter and submillimeter waves [7,8,9]. The most interesting and detailed results were reported in [9], where mesospheric ozone was observed during the solar eclipse of 15 January 2010 at frequencies of about 625 GHz from aboard the International Space Station. It was found in [9] that changes in the ozone mixing ratio (OMR) profiles at altitudes of 58–70 km in dependence on the Sun’s obscuration during the eclipse were in good agreement with predictions of the comprehensive work on ozone chemistry in the mesosphere and lower thermosphere (MLT region, altitudes A = 50–100 km) [2]. Other observations [5,6,7,8] were not as complete and precise as [9], and their data were not robustly analyzed. Atmospheric pressure and concentrations of ozone drop rapidly with increasing altitude. Thus, it is most convenient to use the OMR profile CO3 (A) expressed in ppm (parts per million, 10−6) units instead of the ozone concentration [O3] profiles expressed in cm−3.
At night in the MLT region, ozone is not destroyed by solar radiation, and its content increases greatly. The secondary maximum in the OMR profile appears near 90 km in addition to its main stratospheric maximum at 30–35 km [1,2,10]. Additionally, the tertiary OMR maximum may appear at night at altitudes of 60–70 km in polar and subpolar latitudes [11]. Ozone observations at the MLT region are carried out mainly from satellites in optical (UV, visible, IR), submillimeter, and millimeter (microwave) ranges [12,13,14,15,16,17]. The most effective ground-based technique is remote sensing at microwaves [18,19,20,21,22,23]. The advantages of ground-based microwave remote sensing of ozone and other minor gases of the middle atmosphere are related to its relatively low dependence on weather conditions, wide altitude range and possibility of both day- and night-time observations.
In this paper, we present and discuss the results of ground-based observations of mesospheric ozone over the Moscow region at a frequency of 142.2 GHz during the partial solar eclipses of 20 March 2015 and 25 October 2022. The main purposes of this work were to estimate changes in the OMR over the Moscow region during the eclipses and explain noticeable differences in the observational data obtained for the two events.

2. Microwave Spectrometer

The Microwave Ozone Spectrometer MOS-4 [24,25] used for observations during the solar eclipses operates at frequencies of 100,10–101,9 rotational ozone line centered at 142.175 GHz (wavelength of 2.1 mm). The instrument consists of a room-temperature heterodyne receiver based on a planar Schottky diode mixer and two spectrum analyzers. A simplified block diagram of the MOS-4 is presented in Figure 1, and its main parameters are listed in Table 1. Input Gaussian optics of the instrument provide a low loss of signal power. The chopper wheel switches the receiver input between the sky (or hot load 1 at the calibration stage, see lower) and reference load. The reference load includes cold (cooled by liquid nitrogen) and hot blackbody loads (hot load 2 in Figure 1) and a polarizing wire grid between them. The physical temperature of both hot loads is 5 K above room temperature. Thus, for polarization selected by the mixer waveguide brightness temperature of the reference load varies from TN = 89 K to T0 ≈ 300 K when rotating the grid from 0° to 90°. A single sideband (SSB) filter suppresses the image band of the mixer around 134.625 GHz. Both the SSB filter and diplexer are Martin–Puplett (polarizing Michelson) interferometers. Path length modulator (PLM) reduces distortions of the ozone spectra because of standing waves in the input optics. Frequencies of all local oscillators (LO) are stabilized by phase-locked loops (PLL) based on a high-quality 10 MHz quartz oscillator (not shown in Figure 1, see Table 1).
The first spectrum analyzer is a broadband acousto-optical spectrum analyzer (AOS) of about 0.5 GHz total band with a frequency resolution of 1 MHz. The AOS is used for the stratospheric ozone measurements. Additionally, the second analyzer is the SA-50 filter bank, providing high-frequency resolution from 0.1 to 0.6 MHz in the line center. The frequency characteristics of the SA-50 channels were sufficiently close to rectangular ones. The SA-50 is used for observations of mesospheric ozone. A part of the mixer output signal at the first intermediate frequency (IF1) is the square law detected through the 0.6 GHz frequency band, including the ozone line, to measure the average brightness temperature of the atmosphere in the band. Then, the measured value is used to estimate the tropospheric loss for further correction of the ozone spectra.
The MOS-4 is placed in a laboratory room at the P.N. Lebedev Physical Institute of the Russian Academy of Sciences (LPI) in Moscow (55°42′N, 37°34′E, 180 m above sea level) in front of a radio transparent window of two layers of thin (0.05 mm) polyethylene film (see Figure 1). A small reflection loss in the window of about 0.1 dB is taken into account when processing the ozone spectra. Local objects (the window frame, neighbouring houses and trees) limit the zenith angle (ZA) of the MOS-4 antenna to a range of 50–65°. Observations of the solar eclipses were carried out at ZA = 60°. The line of sight of the spectrometer was directed to the northwest with an azimuth of 315°.
An observation cycle of MOS-4 consists of a calibration stage and measurement stage(s). In the calibration stage, the input beam is closed by hot load 1 with a brightness temperature of T0 ≈ 300 K and the polarizing grid of the reference load is rotated to a position providing the minimum brightness temperature TN = 89 K (see Figure 1). Thus, at this stage, signals corresponding to brightness temperature “step” T0TN ≈ 211 K are recorded in the two spectrum analyzers and broadband detector. At the measurement stage(s), hot load 1 is removed from the input beam and the polarizing grid in the reference load is rotated to a position providing a brightness temperature close to the average brightness temperature of the sky within the AOS and broadband detector frequency bands. This tuning of the reference load minimizes distortions of the ozone line during the measurements. All signals recorded during the measurement stage are scaled by the temperature “step” recorded at the calibration stage and are referred to as the reference load brightness temperature. The result of the measurement is that the sum of the ozone line is reduced in the troposphere and the tropospheric background due to molecular oxygen and water vapour emission.
During observations of the solar eclipses, the calibration/measurement time was 200/200 s on 20 March 2015 and 80/80 s on 25 October 2022. The number of measurement stages was 6 on 20 March 2015 and 15 on 25 October 2022. Thus, the observation cycle duration was about 25 min in both cases.

3. Eclipses

All information about the solar eclipses of 20 March 2015 and 25 October 2022 in this section and hereafter is referred to as point P75 located on the MOS-4 line of sight at an altitude of 75 km above the Earth. The altitude was selected as the average one for the MLT region. The vertical projection of the point to the Earth’s surface has coordinates of 56°32′N, 36°05′E and is located 130 km to the northwest of the LPI. Local Moscow time (MT) will be used in the paper. The Moscow time is MT = UTC + 3 h.
The 20 March 2015 eclipse at point P75 started at 12:11 MT and ended at 14:26 MT. The maximum magnitude of 0.694 and obscuration of 0.625 were reached at 13:19 MT at a solar zenith angle (SZA) of 57.2°. (The magnitude and obscuration of the eclipse are fractions of the Sun’s diameter and area occulted by the Moon, respectively.) The solar noon on this day occurred at 12:43 MT, and the eclipse culmination was only 36 min after the noon. During the eclipse, the SZA varied from 56.8° to 60.4°. The ratio of the Moon-to-Sun diameter at the eclipse was equal to 1.049. The time dependencies of the 20 March 2015 eclipse magnitude and obscuration are presented in Figure 2. For convenience, time in the figure is expressed in percent of the total duration of the eclipse as being equal to 135 min.
The solar eclipse of 25 October 2022 started at point P75 at 12:21 MT, just after the solar noon, which was at 12:19 MT and ended at 14:49 MT. The maximum magnitude of 0.720 and obscuration of 0.646 were reached at 13:36 MT at the SZA of 70.5°. Thus, the eclipse culmination was 77 min after the noon. During the eclipse, the SZA varied from 68.7° to 75.3°. The ratio of the Moon-to-Sun diameters at the eclipse was equal to 0.992. The time dependencies of the eclipse magnitude and obscuration are also presented in Figure 2. The total duration of the eclipse was 148 min.
It is seen through a comparison of the data presented in Figure 2 that the maximum magnitudes and obscurations of the two eclipses were close to each other. The greatest difference between the events was in the SZA values. In October 2022, the SZA was 12–15° higher than in March 2015. At altitudes of 50–80 km, the mesosphere is not completely transparent for optical and UV solar radiation. Thus, the increase in SZA up to 70° and more could result in some decrease in the intensity of the solar radiation and changes in its spectrum.

4. Observational Data

Measurements of mesospheric ozone with the MOS-4 in March 2015 were carried out on all nights from 15 to 20 March. The night measurements started shortly after the sunset at the ground level, continued throughout the night, and finished 0.5–1 h after the sunrise at the ground level. It allowed us to record both the night- and daytime ozone spectra. On 20 March, the measurements were continued in the daytime without turning off the spectrometer from 10:20 MT to 16:30 MT to cover both the solar eclipse and approximately two-hour time intervals before and after the eclipse. Thirty-three units of ozone spectra with an integration time of 200 s were recorded by the SA-50 during the eclipse. In March 2015, it was impossible to use the AOS and SA-50 simultaneously, and priority was given to the SA-50. On all the days from 13 to 20 March 2015, the weather was consistently clear, sunny, and stable, providing favourable conditions for the observations. During the solar eclipse of 20 March 2015, ground level temperature was +9 °C. The tropospheric loss in the eclipse culmination was estimated to be 1.84 dB for the antenna ZA = 60°.
On the contrary, in October 2022, the weather was not stable, and 25 October was a rare clear sunny day among cloudy and gloomy ones. The ground level temperature during the solar eclipse was +4°…+5 °C, lower than on 20 March 2015. This resulted in a decreased tropospheric loss of 1.48 dB in the eclipse culmination. Thus, the loss difference between the two eclipses was 0.36 dB (9%). On 25 October 2022, the observations started at 10:50 MT and were performed continuously until 22:15 MT. Thus, ozone spectra were recorded in the daytime before, during and after the solar eclipse as well as in the evening and in the night-time. A reduced measurement time of 80 s was used when recording the unit spectra on this day, and 89 unit of spectra were measured during the eclipse.
The average daytime ozone spectrum for 25 October 2022 (out of the eclipse period) measured by the MOS-4 in a frequency band of about 0.5 GHz is shown in Figure 3. The spectrum is a fusion of the AOS and SA-50 spectra. It has a narrow central spike and broad wings. The shape of the ozone line corresponds to ozone emission from different layers of the stratosphere and the mesosphere, where pressure exponentially drops with increasing altitude. In this paper, our attention will be focused on changes in brightness temperature TB around the ozone line peak during the solar eclipses and transitions from day- to night-time. Computer simulations [21,22] have shown that the main contributions to the line from ozone at altitudes above 60 km fall into frequency offsets δf from the line centre within approximately ±2 MHz.
It should be noted here that both the spectrum in Figure 3 and all other spectra considered in this paper were corrected for tropospheric loss due to the absorption of microwaves in water vapour and molecular oxygen. The contribution from the water vapour and molecular oxygen to the experimental spectra and correspondent tropospheric loss were determined from the brightness temperature measured using the MOS-4 broadband detector. The contribution was subtracted from the spectra, and the residual (ozone line itself) was multiplied by the tropospheric loss. Thus, all the spectra in this paper are so-called “out-of-troposphere” ones.
Changes in the ozone line spike for the solar eclipses are presented in Figure 4 and Figure 5. In these figures, red lines are daytime ozone spectra from the SA-50 just before and after the eclipses, blue lines are the spectra measured during the central third (33–67%) of the eclipses time, and black lines are average night-time spectra for the closest nights. The ozone line shown in Figure 3 is symmetrical for frequency offsets up to several hundred megahertz. Thus, the spectra in Figure 4 and Figure 5 were obtained by folding their left and right wings to increase the signal-to-noise ratio by √2 times. The zigzag feature seen on the spectra in Figure 4 at frequency offsets of around 1.5 MHz appeared as a result of detuning two channels in the SA-50. This malfunction was eliminated after technical maintenance and tuning of the device in 2016.
The unit spectra of the SA-50 measured during the eclipses were divided into several groups to obtain partial spectra for different phases of the eclipses (see Section 5 for more details). It should be noted here that the noise levels of the partial spectra and accuracy of the OMR estimates from the spectra are inversely proportional to the root square of the total integration time of the spectra included in a group. Therefore, the number of parts into which the total eclipse time could be divided was limited by spectrometer noise.
R E = T B ( 0 ,   Eclipse ) T B ( 0 ,   Day ) T B ( 0 ,   Night ) T B ( 0 ,   Day )
It can be clearly seen through a comparison of Figure 4 and Figure 5 that the observed increase in the amplitude of the ozone line (the line amplitude is brightness temperature TB (0) in the line centre at δf = 0) during the 20 March 2015 eclipse was significantly greater than during the 25 October 2022 eclipse. In fact, the difference between blue and red lines at zero frequency offset TB (0, Eclipse)–TB (0, Day) was 1.2 K on 20 March 2015 and only 0.5 K on 25 October 2022. Since the changes in the ozone line amplitude resulted from changes in mesospheric ozone content and distribution, this means that during the 20 March 2015 eclipse, the increase in mesospheric ozone was significantly greater than for the 25 October 2022 eclipse. As a result, the ratio was equal to 0.41 for 20 March 2015 and only 0.13 for 25 October 2022. Here, TB (0, Night) is the amplitude of the average night-time ozone line. This discrepancy was surprising because the obscurations of the eclipses were close to each other (see Figure 2). This fact was the main subject of further discussion in this paper.
The similar ratio
R ( t ) = T B ( 0 ,   t ) T B ( 0 ,   Day ) T B ( 0 ,   Night ) T B ( 0 ,   Day )
was calculated for evenings of March 2015 and 25 October 2022 in dependence on time t to see temporal changes in the ozone line amplitude during transitions from day- to night-time. In Figure 6, the data are presented for the evenings of 17–19 March 2015 and 25 October 2022. Zero at the horizontal axis corresponds to sunset at an altitude of 100 km at the spectrometer line of sight, i.e., this point marks the disappearance of direct solar radiation at all altitudes lower than 100 km at the line. The vertical dashed line indicates the sunset at the ground level at the LPI. Numbers 0 and 1 at the vertical axis correspond to average day- and night-time amplitudes of the ozone line in accordance with the relation (2). The horizontal bars indicate time intervals of the observations, and the vertical bars denote errors calculated for the ratio R(t) based on the spectrometer noise corresponding to the intervals. It is seen from Figure 6 that in the three evenings preceding the solar eclipse of 20 March 2015, the amplitude of the ozone line increased faster than in the evening of 25 October 2022.

5. Ozone Profiles

Collisional broadening of microwave spectral lines of ozone and other minor gases in the stratosphere and lower mesosphere is strongly dependent on pressure. This dependence is the key feature allowing for the retrieval of the altitude profiles of the species [18]. At the LPI, ozone profiles are retrieved from the broadened 142.175 GHz ozone spectral line using an iterative algorithm based on Tikhonov’s method [21,22]. However, at altitudes above 60 km, the dependence of the ozone line width on pressure (on altitude) weakens because of the growing role of the Doppler broadening. Above 75 km, the Doppler broadening becomes greater than the collisional one. This results in a restriction on the maximum altitude up to which the ozone profile can be retrieved by the algorithm without using a priori information about the true shape of the profile [21,22]. Thus, additional data about the ozone profile at altitudes above 60 km is necessary for the profile retrieval both at night and during the solar eclipses.
It was assumed in [21] that the predominant part of the secondary night OMR maximum at an altitude of 90 km C90CO3 (A = 90 km) was a Gaussian one with a 10 km full width at half maximum (FWHM). The amplitude of the Gaussian contribution must be slightly reduced compared to C90 to take into account the small daytime OMR value, typically of order 0.1 ppm. Similarly, the night OMR value at 65 km C65CO3 (A = 65 km) was assumed to be the sum of the daytime OMR at the altitude and the tertiary night maximum Gaussian contribution centered at 65 km with 15 km FWHM. Then, the a priori night OMR profile was built as the sum of the daytime ozone profile and the two night-time Gaussians. In [21], the profile was used as the initial one for the following iterative retrieval procedure. However, it was found that the procedure practically did not change the ozone profile around 90 km, and noticeable changes appeared only at altitudes lower than 75 km. In this work, as will be seen from the further discussion, the most valuable information about mesospheric ozone variations during the solar eclipses was obtained only for altitudes around 90 km. Thus, we did not perform the full retrieval procedure for the OMR profiles during the night-time and time of the eclipse and used the above-described ozone profiles consisting of the three components.
The daytime ozone profiles for 20 March 2015 and 25 October 2022 retrieved using the iterative Tikhonov algorithm [21,22] are presented in Figure 7 and Figure 8 through the use of green lines. It should be noted that above 75 km, the profiles are close to the initial smooth reference profiles used in the retrieval procedure. To estimate the OMR values at altitudes of 90 and 65 km C90 and C65 at night-time and during the solar eclipses, we used a method initially proposed in [21] and modified in [23]. An experimental ozone spectrum TBf) was compared with the model spectra TB mod (C90, C65, δf) from the “out-of-troposphere” database calculated for some typical daytime ozone profile and a broad range of the OMR maxima values within C90 = 0…20 ppm and C65 = 0…5 ppm at the standard zenith angle ZA = 60° and a temperature of 206 K at an altitude of 85 km. Brightness temperatures of the ozone line in the database were calculated for the SA-50 frequency offsets δfi = 0, 0.05, 0.15, …, 0.85 MHz, i = 1, 2, …, 10. In the method, the maxima C90 and C65 are determined as values minimizing the five-term quadratic form
Q = i = 1 5 ( Δ T B i Δ T B   mod   i ) 2 min
where
Δ T B i = T B ( δ f i ) T B ( δ f i + 1 )
and
Δ T B   mod   i = T B   mod ( C 90 ,   C 65 ,   δ f i ) T B   mod ( C 90 ,   C 65 ,   δ f i + 1 )
are differences in brightness temperature between the neighbouring points with frequency offsets of δfi and δfi+1 for the folded experimental spectrum and the model spectrum. The largest frequency offset used in the form Q is δf6 = 0.45 MHz. This means that Formula (3) takes into account all of the frequencies within the ±0.5 MHz band around the ozone line center. The computer simulations [23] showed that the rms errors of the C90 and C65 values estimated using Formula (3) were proportional to the noise σ of the experimental ozone spectra at least for small σ < 0.1 K. Dependences of the OMR errors on the σ were found to be
σ C 90 = ξ 90   σ ,
σ C 65 = ξ 65   σ
where ξ90 = 14 ppm/K and ξ65 = 2.4 ppm/K. A significant increase in σ above 0.1 K may result in positive shifts (overestimation) of the calculated values of the OMR maxima C90 and C65 because of the appearance of non-monotonic parts on the experimental spectra. Slight non-monotonicity is noticeable in the blue curves in Figure 4 (σ = 0.17 K) and Figure 5 (σ = 0.15 K). For the correct determination of the “out-of-troposphere” σ values for the spectra shown in Figure 3, Figure 4 and Figure 5, the spectrometer noise was multiplied by the same tropospheric loss as the ozone spectra themselves. Using the temperature of 206 K at an altitude of 85 km instead of its real values in the eclipse days estimated within 179–191 K (see Section 6 for more details) could result in additional errors of the OMR values C90 and C65 within 10–15% [21].
Using this method, the average night-time OMR values C90 and C65 were found. At an altitude of 90 km, the C90 was 6.3 ± 0.6 ppm in the night of 19–20 March 2015 and 9.65 ± 1.1 ppm in the night-time late on 25 October 2022. The increased rms error for 25 October 2022 was explained by the shorter observation time. The average OMR values at 90 km for all the nights from 15 to 20 March 2015 are presented in Figure 9. The rms errors of the average night-time ozone spectra were σ ≈ 0.04 K. In accordance with relation (6), the calculated rms errors of the C90 were σC90 ≈ 0.6 ppm. No strong changes in the OMR at 90 km during the St. Patrick’s severe two-step geomagnetic storm of 17–18 March 2015 [26,27] are seen in the figure. During the period of 15–20 March 2015, the average night OMR at 90 km smoothly varied between 7.3 and 6.1 ppm.
At an altitude of 65 km, the average night OMR C65 was 0.53 ± 0.1 ppm on the night of 19–20 March 2015 and 0.69 ± 0.2 ppm in the night-time late on 25 October 2022. On the other nights of 15–20 March 2015, the C65 varied within (0.53–0.74) ppm ± 0.1 ppm.
To study the changes in ozone content at 90 and 65 km during the eclipses, the sets of eclipse unit spectra were divided into groups in four ways and averaged in every group. The photochemical lifetime of ozone in the mesosphere is about 100 s [1,2], and this small delay in ozone response to changes in solar radiation flux was neglected. The first way (I) was to average all of the spectra for every eclipse. In the second way (II), the total eclipse time was divided into four equal parts, tE = (0–25), (25–50), (50–75), and (75–100). Hereafter the eclipse time tE is expressed in percent, as in Figure 2. It can be seen from Figure 2 that the eclipse curves are symmetrical with respect to the culmination point tE = 50%. Thus, the spectra falling into the first and fourth quarters contained identical information in terms of ozone, just like the spectra falling into the second and third quarters. For this reason, the spectra of these parts were united. The united groups were named (0, 25) and (25, 50) after the names of the first subgroups in every pair. Similarly, the groups in way (I) were named (0, 50) because, formally, it contained the spectra of two symmetrical parts, tE = (0–50) and (50–100). In the third (III) and fourth (IV) ways, the total eclipse time was divided into six and eight equal parts, tE = (0–16.7), …, (83.3–100) and tE = (0–12.5), …, (87.5–100), respectively. After combining the unit spectra falling into parts symmetrical with respect to the point tE = 50%, three and four groups were obtained for the ways. Similar to ways (I) and (II), these groups were named (0, 16.7), (16.7, 33.3), (33.3, 50) for the way (III) and (0, 12.5), (12.5, 25), (25, 37.5), (37.5, 50) for the way (IV). It should be noted that the total numbers of the unit spectra (33 for 20 March 2015 and 89 for 25 October 2022), in most cases, could not be divided into groups with quite equal numbers of spectra. Thus, the division was performed as accurately and reasonably as possible.
Then, the ozone mixing ratios C90 (tE1, tE2) and C65 (tE1, tE2) at altitudes of 90 and 65 km were found for all groups using the above-described method. Here, tE1 and tE2 are the time boundaries of the groups in percent. The OMR values at 90 km are presented in Figure 10, together with estimates of their rms errors σC90 calculated using the relation (6). Successive rows of panels in Figure 10 from top to bottom correspond to the division of the unit spectra by the ways from (I) to (IV). Data on the OMR C90 (tE1, tE2) during the solar eclipse of 20 March 2015 are shown in the left column (panels a–d), and the same data on the eclipse of 25 October 2022 are presented in the right column (panels e–h). It is seen that the OMR in the left column, as a whole, is greater than that in the right one. This directly corresponds to the differences in the observed mesospheric ozone responses to the eclipses (see Figure 4 and Figure 5).
Additionally, it is noticeable in Figure 10 that the rms errors σC90 are large and increase from the top to the bottom panels. The errors are inversely proportional to the root square of the number of the unit spectra in the groups. It is seen that on 25 October 2022, the errors were a little less than on 20 March 2015. This occurred because of two reasons. Firstly, the eclipse of 25 October 2022 lasted a little longer than the eclipse on 20 March 2015 (148 min versus 135 min), and secondly, the troposphere on 25 October 2022 was a little more transparent (loss of 1.48 dB, see Section 4) than on 20 March 2015 (loss of 1.84 dB). Typically, at an altitude of 90 km, the errors σC90 were comparable or even exceeded the OMR values C90 close to the eclipse culminations.
Unfortunately, those calculated by the groups’ OMR values C65 at an altitude of 65 km were not large and regular enough to reveal their dependences on the solar obscuration η. It was especially difficult since the expected dependence C65(η), in accordance with Formula (8), lower with exponent γ ≈ 0.5, was rather weak itself. For the solar eclipse on 20 March 2015, the OMR C65 (33.3, 66.7) was found to be 0.46 ± 0.42 ppm, but for the eclipse of 25 October 2022, it was estimated to be less than 0.2 ppm, with an rms error of 0.36 ppm. Thus, further in this paper, we discuss only ozone changes at an altitude of 90 km.

6. Discussion

M. Allen et al. [2] have derived formulas for the ozone concentration at different layers of the MLT region. For altitudes of 50–80 km, the formulas can be used only in daylight hours, and at 90 km and above—they can be used around the clock. K. Imai et al. [9] proposed a simplified form of the formula for ozone at altitudes between 50 and 70 km to present their data on the OMR profiles at the altitudes during the solar eclipse of 15 January 2010 in dependence on the Sun’s obscuration η
C O 3 ( A ,   η ) = C O 3 ( A ,   0 ) ( 1 η ) γ
In accordance with [2], at altitudes of 90 km and above, the concentration of ozone is equal to
[ O 3 ] = [ O ] [ O 2 ] ( k 12 [ O ] + k 13 [ O 2 ] + k 14 [ N 2 ] ) J 10 + J 11 + k 16 [ H ] + k 15 [ O ] .
Here, J10 and J11 are rates of ozone photodissociation due to the solar radiation with frequency ν
O 3 +   h ν         O 2 +   O ,
O 3 +   h ν         O 2 +   O ( D 1 ) ,
h is the Plank constant, [O], [O2], [N2], and [H] are concentrations of atomic and molecular oxygen, molecular nitrogen, and atomic hydrogen, respectively, and k12k16 are rate constants of the ozone producing and destructive chemical reactions with these species:
O   +   O 2 +   O         O 3 +   O ,
O   +   O 2 +   O 2         O 3 +   O 2 ,
O   +   O 2 +   N 2         O 3 +   N 2 ,
O   +   O 3         2 O 2 ,
H   +   O 3         OH   +   O 2 .
These constants and their dependencies on temperature T [2] are listed in Table 2.
The rate constants k13 and k14 are equal to each other (see Table 2), but at 90 km, the concentration of atomic oxygen is much less than that of molecular oxygen: [O] << [O2]. Thus, the first term in the parenthesis of the nominator in Formula (9) is small. Further, at altitudes of about 90 km, the contribution from the reaction (16) is 1–2 orders higher than one from the reaction (15). As a result, the next notations can be introduced:
S = [ O ] [ O 2 ] ( k 12 [ O ] + k 13 [ O 2 ] + k 14 [ N 2 ] ) [ O ] [ O 2 ] ( k 13 [ O 2 ] + k 14 [ N 2 ] )
D = k 16 [ H ] + k 15 [ O ] k 16 [ H ]
J = J 10 + J 11 ,
and one can write Formula (9) in daylight hours as follows:
[ O 3 ] ( Day ) = S J + D
Here, S is expressed in cm−3 s−1, while J and D are expressed as s−1. The lifetimes of atomic hydrogen and oxygen at 90 km are of about 0.5 h and 1 week, respectively. The lifetimes highly increase above 90 km [1,2]. Molecular oxygen and nitrogen are stable components. At night the radiative term J is equal to zero and
[ O 3 ] ( Night ) = S D .
Provided that chemical term D does not change noticeably from day to night, one can rewrite Equation (13) for daytime ozone concentration as follows:
[ O 3 ] ( Day ) = D [ O 3 ] ( Night ) J + D .
Dividing the left and right parts of relation (21) by the total concentration of molecules in the air results in a formula for daytime OMR values:
C O 3 ( A ,   Day ) = D C O 3 ( A ,   Night ) J + D
The large night-time OMR maximum C90 (Night) at an altitude of 90 km is well-defined from the microwave measurements (see Figure 9) in contrast to ill-defined small daytime OMR values C90 (Day) (typically C90 (Day) < 1.2 ppm [10]). Since the ratio C90 (Night)/C90 (Day) >> 1, the chemical term, D, should be small, D << J, as clearly seen from Formula (23). It follows from Formulas (17), (18) and (21) that the provided temperature and pressure are constant, and both the night-time ozone concentration and its mixing ratio at 90 km are proportional to the ratio of atomic-oxygen-to-atomic-hydrogen concentrations [O]/[H].
Similar to Equation (8), during the solar eclipses, the denominator in Formula (23) should be reduced because of the Sun’s obscuration:
C O 3 ( A ,   η ) = D C O 3 ( A ,   Night ) J ( 1 η ) + B + D
An additional positive term, B, was inserted into the denominator of Formula (24) to model the possible contribution of indirect (scattered and reflected) solar radiation to mesospheric ozone destruction. It was assumed that the term was small B << J and constant during the eclipse because the scattered and reflected radiation could come from the regions far outside of the Moon’s shadow. At night, the term disappeared together with solar radiation J. However, the transition from day to night for the indirect radiation could be completed later. Among the others, the term B in Formula (24) turned out to be a suitable fitting parameter when analyzing the ozone data at 90 km. For convenience, we divided both the nominator and denominator in the right part of relation (24) by J
C O 3 ( A ,   η ) = D J C O 3 ( A ,   Night ) 1 η + B J + D J
to reveal the dependences of the OMR during the eclipses on the dimensionless ratios BJ = B/J, DJ = D/J, the above-mentioned night-time OMR values at 90 km C90 = 6.3 ppm and C90 = 9.65 ppm were used in Formula (25) for the eclipses of 20 March 2015 and 25 October 2022.
The solar eclipses of 20 March 2015 and 25 October 2022 were separated by a time interval of about 7.5 years and fell into different stages of two adjacent 11-year solar cycles. The first eclipse occurred at the decline of the 24th solar cycle, and the second—at the beginning of the 25th one. However, the F10.7 indexes of solar activity [28] on the dates of eclipses were close to each other, 100.6 on 20 March 2015 and 103.6 on 25 October 2022. It was mentioned above that the solar zenith angles were different during the eclipses. The SZA varied from 57° to 60° on 20 March 2015 and from 68° to 75° on 25 October 2022. However, at altitudes of around 90 km, the atmosphere is sparse, and this difference could only slightly weaken the solar radiation on 25 October 2022, partly compensating for the difference in fluxes indicated by the F10.7 indexes. Thus, it was assumed that direct fluxes of solar radiation J in Formulas (20), (22), and (23) were equal during the two eclipses. On the contrary, it seems possible that the increased SZA on 25 October 2022 could result in some increase in indirect radiation flux B due to broadening the area where the reflection of solar radiation occurred at the sliding angles. However, any quantitative estimates of this effect are far from the scope of this paper.
To find dimensionless parameters BJ and DJ, which provide the best fitting of relation (25) to the experimental data C90 (tE1, tE2) presented in Figure 10, the parameters varied within sufficiently broad ranges of BJ = −0.37…1 and DJ = 0.025…0.3 to minimize the quadratic form of F(BJ, DJ)
F ( B J ,   D J ) = n = 1 N [ C 90 ( t E 1 n ,   t E 2 n ) 1 t E 2 n t E 1 n   t E 1 n t E 2 n D J C 90 ( Night ) 1 η ( t E ) + B J + D J     d   t E ] 2   min .
The number N of terms summed up in the right part of Formula (26) increases from N = 1 for way (I) to N = 4 for way (IV). Two sets of curves BJ (DJ) providing the fitting are presented in Figure 11. The blue lines refer to the eclipse of 20 March 2015 and the red lines to the eclipse of 25 October 2022. The letters near the curves correspond to the notations of panels in Figure 10. It can be seen that, despite the significant spread of the blue curves, they all move much to the right of the more compact set of red curves. This means that the magnitude of the chemical term DJ in Equation (25) was much larger during the first eclipse. To estimate this term, we separately calculated for every curve BJ (DJ) and averaged through the blue and red sets of the curves the DJ values for points where the curves crossed the levels BJ = 0, 0.1, and 0.2. These small levels of indirect solar radiation seemed to be more or less realistic ones. It was found that for the eclipse of 20 March 2015, DJ 20 was equal to 0.14, 0.16, and 0.18 for the three levels of BJ. Similarly, the same values of the DJ 25 for the eclipse on 25 October 2022 were 0.046, 0.055, and 0.064. As a result, the ratio of DJ 20/DJ 25 was found to be approximately constant for all three levels of BJ and equal to 2.9–3.0.
Therefore, it was concluded that the magnitude of indirect solar radiation B, even if it differed from zero and was different for the two solar eclipses, had little effect on ozone variations at 90 km during the eclipses. Thus, the main reason for the strongly different observed responses of ozone at 90 km during the two eclipses with approximately equal obscurations of the Sun should have large differences in terms of the values of the chemical terms D and DJ in Formulas (24) and (25).
The error of ratio DJ for every curve in Figure 11 was estimated through the use of computer simulations performed at BJ = 0. The rms errors σC90 of the OMR values C90 shown in Figure 10 were used in the simulations. As one might expect, these errors turned out to be large. They were lying within 0.21–0.26 for the blue curves a–d and within 0.11–0.13 for the red curves e–h, with some tendency in terms of decreased errors when there was an increase in N from 1 to 4. For the blue and red curves, the errors were correlated because the same sets of unit ozone spectra were used for the calculation of DJ at different N. However, it was too difficult to model the errors of the average DJ values. This would require varying the brightness temperatures in every channel of every unit spectrum for every eclipse. Thus, the errors of the average DJ values at BJ = 0 were estimated as falling in intervals (0.12, 0.24) on 20 March 2015 and (0.06, 0.12) on 25 October 2022, where the smaller numbers correspond to uncorrelated errors within the sets of blue and red curves, and the larger ones correspond to fully correlated errors. Despite the large errors, we believe that the data shown in Figure 11 quite clearly discriminated the values of the DJ for the two solar eclipses.
Three possible explanations for the three-times difference in the magnitudes of the chemical term D for the two solar eclipses were considered:
(1) The severe geomagnetic storm of 17–18 March 2015 [26,27] could have noticeably accelerated the chemical destruction of ozone in the upper mesosphere because of the increased HOx content (mainly atomic hydrogen H) and NOx (mainly nitric oxide NO) radicals resulting from the ionization of the upper mesosphere by charged particles of solar plasma [3,4,27]. Excessive atomic hydrogen should have disappeared after several hours [1,4], but the lifetime of mesospheric NO is about 2 days [1]. Thus, the increase in NO content could have been partially maintained by 20 March 2015 and ozone destruction in the reaction
NO   +   O 3         NO 2   +   O 2
could increase the chemical term D by the time of the solar eclipse on 20 March 2015. However, any strong decrease in night ozone was not seen at 90 km after 16 March 2015, as shown in Figure 9. The night-time OMR C90 (Night) decreased by 8%, from 6.9 ppm on the night of 16–17 March, just before the geomagnetic storm, to 6.35 ppm on the night of 17–18 March during the storm. Then, the decrease in ozone content slowed down to 6.1 ppm on the night of 18–19 March and was replaced by a slow increase up to 6.3 ppm on the night of 19–20 March. Therefore, the geomagnetic storm did not significantly increase the chemical term D at an altitude of 90 km on 20 March 2015. Possible changes in D resulted from the storm hardly exceeding approximately 10%. It should also be noted that the maximum destruction of ozone by NO molecules is located lower, at altitudes of about 80 km [3,4], and ozone at 90 km is less sensitive to possible increases in NO content;
(2) The rate constants k12k16 of chemical reactions (12)–(16) are temperature-dependent (see Table 2) [2]. The reaction (15) has a very strong dependence on temperature, but it plays only a minor role in ozone destruction at altitudes of about 90 km. The highest points of the temperature profiles available from the ERA5 reanalysis database [29] are located at altitudes of about 78 km. On 20 March 2015, the upper temperature points of the profile were 200.53 K at 78.09 km and 212.83 K at 72.43 km. Similarly, on 25 October 2022, the values were 192.71 K at 77.63 km and 201.21 K at 72.21 km. All the data referred to 12:00 MT. Linear interpolation of the ERA5 temperature profiles to the lower mesopause altitude of 85 km resulted in temperatures of about 183 K and 179 K for the two dates. Otherwise, spline interpolation of the ERA5 profiles using the average model temperature of 191 K at an altitude of 100 km [30] resulted in temperatures at 90 km of 191 K on 20 March 2015 and 179 K on 25 October 2022. According to Table 2, the increase in temperature by 12 K from 179 K to 191 K results in an increase in D by 9% and a decrease in S by 19%. For the smaller temperature increase of 4 K, the increase in D is only 3%, and the decrease in S is 6%. Obviously, the temperature dependence of the chemical term D is too small to explain the change by approximately three times in terms of magnitude;
(3) Then, it remains only to assume that the main reason for the large difference in the magnitudes of D was the significantly lower concentration of atomic hydrogen H at altitudes of about 90 km on 25 October 2022 in comparison with 20 March 2015. Taking into account that the above-discussed temperature effects were small, the hydrogen concentration should have decreased by almost three times. This assumption seems consistent with the results of modeling the distribution of atomic hydrogen in the mesosphere presented in [1]. In that book, at latitudes of about 60°, the calculated concentration of atomic hydrogen at 90 km in summer was considerably higher than in winter (approximately 108 cm−3 versus 3 × 107 cm−3). These seasonal variations depend on both eddy transport and photochemical processes [1,2]. Since 20 March 2015 was close to the vernal equinox, but 25 October 2022 was more than a month later than the autumnal equinox, we can assume that the concentration of atomic hydrogen in the second case could be lower.
Substituting the average values of DJ = 0.14 for 20 March 2015 and DJ = 0.046 for 25 October 2022 into Formula (25) at BJ = 0, we obtained the daytime OMR values of (out of the eclipse time) C90(Day) at an altitude of 90 km for the dates. They were at 0.77 and 0.42 ppm, respectively. Similarly, the calculated OMR values at the eclipse culminations C90max) were 1.71 ppm on 20 March 2015 and 1.11 ppm on 25 October 2022. The values satisfactorily matched the data in relation to C90(η) presented in Figure 10. It is interesting to note that the ratio of C90max)/C90(Day) was equal to 2.2 on 20 March 2015 and 2.6 on 25 October 2022. This means that relative changes in ozone at 90 km during the second eclipse were even greater, but they could not be detected with sufficient accuracy because of spectrometer noise.
Thus, the observed change in the ozone content (concentration or OMR) at an altitude of 90 km during a solar eclipse depends mainly on the night-time ozone content and on the concentration of atomic hydrogen at this altitude, determining the chemical term D in Formula (24). For partial solar eclipses with obscuration η < 1 − DJ2, the greater the concentration of atomic hydrogen, the greater the absolute increase in ozone content during the eclipse. (For values DJ < 0.2 found in this work the rule should be correct for all eclipses with η < 0.96). On the contrary, relative changes in the ozone content during any eclipse will increase with a decrease in atomic hydrogen concentration because of the decrease in daytime ozone content.
Additionally, the changes in the magnitude of the chemical term D allowed us to explain the differences in the transition curves from day to night presented in Figure 6. Let us use Formula (24) at B = 0 for the evening time as well. Then, the factor (1 − η) corresponds to the reduction in solar flux J during sunset. It can be seen from (24) that the smaller the D, the less sunlight should remain; therefore, the OMR C90 approaches its night value C90 (Night). Therefore, on the evening of 25 October 2022, at a low value of D, the transition to night conditions occurred later than on the evenings of 16–18 March 2015. An alternative explanation of the differences in the transition curves for March 2015 and 25 October 2022 is the assumption of a much larger amount of indirect radiation flux B on the evening of 25 October 2022. However, this requires unrealistically large values B > 0.2J. Therefore, such an explanation seems unlikely.

7. Conclusions

In this work, we observed microwaves and studied variations in mesospheric ozone O3 over the Moscow region during two partial solar eclipses on 20 March 2015 and 25 October 2022. Both eclipses occurred at midday, with clear weather favourable for observations. The maximum obscuration of the Sun during the eclipses was 0.625 on 20 March 2015 and 0.646 on 25 October 2022. The brightness temperature of the 100,10–101,9 rotational ozone spectral line centered at 142.175 GHz was measured with a frequency resolution of 0.1 MHz near the line center by the ground-based low-noise spectrometer MOS-4 located at the LPI in Moscow. The unit ozone spectra with durations of 200 or 80 s were recorded before, during, and after the eclipses. The spectra were corrected for the tropospheric loss estimated from the measured average brightness temperature of the sky in a broad band of 0.6 GHz. Then, the spectra recorded within the eclipse time were divided into groups with different durations and averaged by the groups to reveal the dependence of the ozone mixing ratio (OMR) in the mesosphere on the obscuration of the Sun at altitudes of 90 and 65 km, where the secondary and tertiary maxima of the OMR profile are located. Special techniques were used to calculate the OMR values of C90 and C65 at these altitudes. It was surprising that despite the close values of the maximum eclipse obscurations, the increases in both brightness temperatures of the ozone line and in the OMR values at 90 km were approximately two times greater on 20 March 2015 than on 25 October 2022.
Unfortunately, the calculated OMR values of C65 at an altitude of 65 km were not large and regular enough to reveal their dependence on solar obscuration. For the solar eclipse on 20 March 2015, the C65 value for the central third of the eclipse time was found to be 0.46 ± 0.42 ppm, but for the eclipse of 25 October 2022, it was estimated to be less than 0.2 ppm, with an rms error of 0.36 ppm.
To understand the origin of the large difference in results at 90 km, the formulas for ozone content in the mesosphere derived by M. Allen et al. [2] were used. Their formula for ozone concentration at altitudes of 90 km and above was modified for use with the OMR values during solar eclipses. Fitting the modified formula to the calculated C90 values resulted in the conclusion that on 25 October 2022 at 90 km, the relative role of the chemical destruction of ozone (compared with its photodissociation by solar radiation) was approximately three times less than on 20 March 2015. Thus, at approximately equal solar fluxes, we had to conclude that absolute rates of the chemical ozone destruction, mainly by atomic hydrogen, differed by three times for the two eclipses.
Night-time measurements of the mesospheric ozone from 15 to 20 March 2015 showed the possibility that the increase in chemical ozone destruction resulted from the severe geomagnetic storm of 17–18 March 2015 preceding the solar eclipse on 20 March 2015, which did not exceed approximately 10% at an altitude of 90 km. Similarly, an estimated decrease in temperature at altitudes of around 90 km within 4–12 K on 25 October 2022 compared to 20 March 2015 could not decrease the rate constant of the reaction of ozone with atomic hydrogen by more than 9%. Consequently, the most probable reason for the three-fold difference in chemical ozone destruction rates was a decrease by almost three times in atomic hydrogen concentration on 25 October 2022 in comparison to 20 March 2015.
Using the modified formula for the OMR at 90 km, we calculated the daytime OMR values C90(Day) (out of the eclipse time) of 0.77 and 0.42 ppm, respectively, for 20 March 2015 and 25 October 2022. At the eclipse culminations, the calculated OMR values C90max) were 1.71 ppm on 20 March 2015 and 1.11 ppm on 25 October 2022. Thus, the ratio C90max)/C90(Day) was equal to 2.2 on 20 March 2015 and 2.6 on 25 October 2022. This means that relative changes in ozone at 90 km during the eclipse of 25 October 2022 were even greater, but they could not be detected with sufficient accuracy because of spectrometer noise.
Thus, the observed changes in the ozone content (concentration or OMR) at an altitude of 90 km during a solar eclipse depend mainly on the night-time ozone content and on the concentration of atomic hydrogen at this altitude. For partial solar eclipses with obscuration η < 1 − DJ2, the greater the concentration, the greater the absolute increase in the ozone content during the eclipse. Otherwise, relative changes in the ozone content during any eclipse will increase with the decrease in the atomic hydrogen concentration. Additionally, the decrease in the rate of chemical ozone destruction results in a later transition from day to night in the mesospheric ozone profile.
The research confirmed that solar eclipses provide unique possibilities for experimentally studying the photochemistry of mesospheric ozone and testing photochemical models. Interesting events, such as sufficiently deep partial and, especially, total solar eclipses, are rare in one place. In Moscow, the next eclipse comparable in the Sun’s obscuration to the eclipses discussed in this paper will not occur until 2039. Probably, microwave spectrometers operating onboard satellites have more chances to observe and investigate the changes in mesospheric ozone profiles during solar eclipses.

Author Contributions

Methodology, S.R. and A.I.; observations, S.R. and A.Z.; auxiliary data, A.I.; software, A.I.; data curation, S.R., A.I. and A.Z.; writing draft paper, S.R.; editing the paper, S.R., A.I. and A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brassuer, G.; Solomon, S. Aeronomy of the Middle Atmosphere, 3rd ed.; Springer: Dordrecht, Netherlands, 2005. [Google Scholar]
  2. Allen, M.; Lunine, J.I.; Yung, Y.K. The vertical distribution of ozone in the mesosphere and lower thermosphere. J. Geophys. Res. 1984, 89, 4841–4872. [Google Scholar] [CrossRef] [Green Version]
  3. Ozaki, M.; Shiokawa, K.; Kataoka, R.; Mlynczak, M.; Paxton, L.; Connors, M.; Yagitani, S.; Hashimoto, S.; Otsuka, Y.; Nakahira, S.; et al. Localized mesospheric ozone destruction corresponding to isolated proton aurora coming from Earth’s radiation belt. Sci. Rep. 2022, 12, 16300. [Google Scholar] [CrossRef] [PubMed]
  4. Grankin, D.V.; Mironova, I.A.; Rozanov, E.V. Polar Winter Mesospheric Ozone Depletion during Energetic Electron Precipitation. In Proceedings of the 27th International Symposium on Atmospheric and Ocean Optics, Atmospheric Physics, Moscow, Russia, 5–9 July 2021. [Google Scholar] [CrossRef]
  5. Randhawa, J.S. Mesospheric ozone measurements during a solar eclipse. J. Geophys. Res. 1968, 73, 493–495. [Google Scholar] [CrossRef]
  6. Agashe, V.V.; Rathi, S.M. Changes in the concentration of mesospheric ozone during the total solar eclipse. Planet. Space Sci. 1982, 30, 507–513. [Google Scholar] [CrossRef]
  7. Connor, B.J.; Siskind, D.E.; Tsou, J.J.; Parrish, A.; Remsberg, E.E. Ground-based microwave observations of ozone in the upper stratosphere and mesosphere. J. Geophys. Res. 1994, 99, 16757–16770. [Google Scholar] [CrossRef]
  8. Kulikov, Y.Y.; Krasil’nikov, A.A.; Demkin, V.M.; Ryskin, V.G. Variations in the concentration of mesospheric ozone during the total solar eclipse of March 29, 2006, from microwave radiometric data. Izv. Atmos. Oceanic Phys. 2008, 44, 486–490. [Google Scholar] [CrossRef]
  9. Imai, K.; Imamura, T.; Takahashi, K.; Akiyoshi, H.; Yamashita, Y.; Suzuki, M.; Ebisawa, K.; Shiotani, M. SMILES observations of mesospheric ozone during the solar eclipse. Geophys. Res. Lett. 2015, 42, 3576–3582. [Google Scholar] [CrossRef] [Green Version]
  10. Smith, A.K.; Marsh, D.R. Processes that account for the ozone maximum at the mesopause. J. Geophys. Res. 2005, 110, D23305. [Google Scholar] [CrossRef] [Green Version]
  11. Marsh, D.; Smith, A.; Brasseur, G.; Kaufmann, M.; Grossmann, K. The existence of a tertiary ozone maximum in the high-latitude middle mesosphere. Geophys. Res. Lett. 2001, 28, 4531–4534. [Google Scholar] [CrossRef] [Green Version]
  12. Bruhl, C.; Drayson, S.R.; Russell, J.M.; Crutzen, P.J.; Mcinerney, J.M.; Purcell, P.N.; Claude, H.; Gernandt, H.; Mcgee, T.J.; Mcdermid, I.S.; et al. Halogen Occultation Experiment ozone channel validation. J. Geophys. Res. 1996, 101, 10217–10240. [Google Scholar] [CrossRef]
  13. Kostsov, V.S.; Timofeev, Y.M. Mesospheric ozone from the CRISTA-1 satellite experimental data: 2. Spatial distributions and diurnal variations. Izvestiya. Atmos. and Oceanic Phys. 2005, 41, 191–202. [Google Scholar]
  14. Kaufmann, M.; Gusev, O.A.; Grossmann, K.U.; Martin-Torres, J.; Marsh, D.R.; Kutepov, A.A. Satellite observations of daytime and nighttime ozone in the mesosphere and lower thermosphere. J. Geophys. Res. 2003, 108, 4272–4285. [Google Scholar] [CrossRef]
  15. KKyrölä, E.; Tamminen, J.; Leppelmeier, G.W.; Sofieva, V.; Hassinen, S.; Seppälä, A.; Verronen, P.T.; Bertaux, J.L.; Hauchecorne, A.; Dalaudier, F.; et al. Nighttime ozone profiles in the stratosphere and mesosphere by the Global Ozone Monitoring by Occultation of Stars on Envisat. J. Geophys. Res. 2006, 111, D24306. [Google Scholar] [CrossRef] [Green Version]
  16. Waters, J.W.; Froidevaux, L.; Harwood, R.S.; Jarnot, R.F.; Pickett, H.M.; Read, W.G.; Siegel, P.H.; Cofield, R.E.; Filipiak, M.J.; Flower, D.A.; et al. The Earth observing system microwave limb sounder (EOS MLS) on the Aura satellite. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1075–1092. [Google Scholar] [CrossRef]
  17. Imai, K.; Manago, N.; Mitsuda, C.; Naito, Y.; Nishimoto, E.; Sakazaki, T.; Fujiwara, M.; Froidevaux, L.; von Clarmann, T.; Stiller, G.P.; et al. Validation of ozone data from the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES). J. Geophys. Res. Atmos. 2013, 118, 5750–5769. [Google Scholar] [CrossRef] [Green Version]
  18. Janssen, M.A. (Ed.) Atmospheric Remote Sensing by Microwave Radiometry; John and Wiley & Sons: New York, NY, USA, 1993. [Google Scholar]
  19. Zommerfelds, W.C.; Kunzi, K.F.; Summers, M.E.; Bevilacqua, R.M.; Strobel, D.F. Diurnal variations of mesospheric ozone obtained by ground-based microwave radiometry. J. Geophys. Res. 1989, 94, 12819–12832. [Google Scholar] [CrossRef]
  20. Kawabata, K.; Ogawa, H.; Yonekura, Y.; Suzuki, H.; Suzuki, M.; Iwasaka, Y.; Shibata, T.; Sakai, T. Ground-based radiometry of stratomesospheric ozone employing a superconductive receiver. J. Geophys. Res. 1997, 102, 1371–1377. [Google Scholar] [CrossRef]
  21. Ignatyev, A.N. Radiometry of the Atmospheric Ozone and Chlorine Monoxide. PhD Thesis, Moscow State University, Moscow, Russia, 2006. [Google Scholar]
  22. Solomonov, S.V.; Gaikovich, K.P.; Kropotkina, E.P.; Rozanov, S.B.; Lukin, A.N.; Ignatyev, A.N. Remote sensing of atmospheric ozone at millimeter waves. Radiophysics. Quantum. Electron. 2011, 54, 102–110. [Google Scholar] [CrossRef]
  23. Rozanov, S.B.; Zavgorodniy, A.S.; Ignatyev, A.N. Microwave measurements of variations in night mesospheric ozone over Moscow. In Proceedings of the SPIE Remote Sensing 2019, Strasbourg, France, 9–12 September 2019. [Google Scholar] [CrossRef]
  24. Rozanov, S.B.; Shtanyuk, A.M.; Bol’shakov, O.S.; Zavgorodniy, A.S.; Kropotkina, E.P.; Kruglov, S.K.; Logvinenko, S.V.; Lukin, A.N.; Nikiforov, P.L.; Saenko, I.I.; et al. Transportable millimeter-wave spectrometer for the atmospheric ozone monitoring. In Proceedings of the 2013 International Kharkov Symposium on Physics and Engineering of Microwaves, Millimeter and Submillimeter Waves, Kharkov, Ukraine, 23–28 June 2013; pp. 483–485. Available online: http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=6622107 (accessed on 27 June 2023).
  25. Rozanov, S.B.; Zavgorodniy, A.S.; Lukin, A.N. Microwave Spectrometer for Ground-Based Measurements of the Stratospheric and Mesospheric Ozone. In Proceedings of the 2019 European Microwave Conference in Central Europe (EuMCE), Prague, Czech Republic, 13–15 May 2019; pp. 614–617. Available online: https://ieeexplore.ieee.org/document/8874699 (accessed on 27 June 2023).
  26. Kamide, Y.; Kusano, K. No major solar flares but the largest geomagnetic storm in the present solar cycle. Space Weather. 2015, 13, 365–367. [Google Scholar] [CrossRef]
  27. Wu, C.-C.; Liou, K.; Lepping, R.P.; Hutting, L.; Plunkett, S.; Howard, R.A.; Socker, D. The first super geomagnetic storm of solar cycle 24: “The St. Patrick’s day event (17 March 2015)”. Earth Planets Space 2016, 68, 151. [Google Scholar] [CrossRef] [Green Version]
  28. Index of /stp/data/solar.act/flux10.7/daily. Available online: http://www.wdcb.ru/stp/data/solar.act/flux10.7/daily/ (accessed on 30 May 2023).
  29. ECMWF Reanalysis v5 (ERA5). Available online: https://www.ecmwf.int/en/forecasts/dataset/ecmwf-reanalysis-v5 (accessed on 30 May 2023).
  30. COSPAR International Reference Atmosphere (CIRA-86). Available online: https://www.sparc-climate.org/data-centre/data-access/reference-climatology/cira-86/ (accessed on 30 May 2023).
Figure 1. Simplified block diagram of the MOS-4. Abbreviations are explained in the text. Central frequencies of the ozone line after frequency conversions to three intermediate frequencies (IF) are given together with the first LO frequency.
Figure 1. Simplified block diagram of the MOS-4. Abbreviations are explained in the text. Central frequencies of the ozone line after frequency conversions to three intermediate frequencies (IF) are given together with the first LO frequency.
Remotesensing 15 03440 g001
Figure 2. Time dependencies of magnitude (red) and obscuration (blue) of the solar eclipses of 20 March 2015 (dashed lines) and 25 October 2022 (solid lines) at the P75 point.
Figure 2. Time dependencies of magnitude (red) and obscuration (blue) of the solar eclipses of 20 March 2015 (dashed lines) and 25 October 2022 (solid lines) at the P75 point.
Remotesensing 15 03440 g002
Figure 3. The average daytime “out-of-troposphere” ozone spectrum measured using MOS-4 on 25 October 2022 (out of the solar eclipse period).
Figure 3. The average daytime “out-of-troposphere” ozone spectrum measured using MOS-4 on 25 October 2022 (out of the solar eclipse period).
Remotesensing 15 03440 g003
Figure 4. The SA-50 ozone spectra measured in daytime of 20 March 2015 before and after the solar eclipse (red), during the eclipse (33–67% of the eclipse time, blue), and in the night of 19–20 March 2015 (black).
Figure 4. The SA-50 ozone spectra measured in daytime of 20 March 2015 before and after the solar eclipse (red), during the eclipse (33–67% of the eclipse time, blue), and in the night of 19–20 March 2015 (black).
Remotesensing 15 03440 g004
Figure 5. The SA-50 ozone spectra were measured in the daytime of 25 October 2022 before and after the solar eclipse (red), during the eclipse (33–67% of the eclipse time, blue), and late on 25 October 2022 (black).
Figure 5. The SA-50 ozone spectra were measured in the daytime of 25 October 2022 before and after the solar eclipse (red), during the eclipse (33–67% of the eclipse time, blue), and late on 25 October 2022 (black).
Remotesensing 15 03440 g005
Figure 6. Relative changes in the amplitude of the ozone line measured in the evenings of 17 March 2015 (red), 18 March 2015 (green), 19 March 2015 (blue), and 25 October 2022 (black).
Figure 6. Relative changes in the amplitude of the ozone line measured in the evenings of 17 March 2015 (red), 18 March 2015 (green), 19 March 2015 (blue), and 25 October 2022 (black).
Remotesensing 15 03440 g006
Figure 7. The OMR profile for the night of 19–20 March 2015 and contributions to it from the daytime ozone profile for 20 March and the two Gaussian maxima at 90 and 65 km determined from the night measurements on 19–20 March.
Figure 7. The OMR profile for the night of 19–20 March 2015 and contributions to it from the daytime ozone profile for 20 March and the two Gaussian maxima at 90 and 65 km determined from the night measurements on 19–20 March.
Remotesensing 15 03440 g007
Figure 8. The OMR profile for night-time late on 25 October 2022 and contributions to it from the retrieved daytime ozone profile and the two Gaussian maxima at 90 and 65 km determined from the measurements after the sunset.
Figure 8. The OMR profile for night-time late on 25 October 2022 and contributions to it from the retrieved daytime ozone profile and the two Gaussian maxima at 90 and 65 km determined from the measurements after the sunset.
Remotesensing 15 03440 g008
Figure 9. The OMR C90 for the nights of 15–20 March 2015. Zero at the time scale corresponds to 00:00 MT 15 March. The dates are written above the timescale as 15.03, 16.03, …, 20.03. Night-time is highlighted with the hatching. Periods of the severe geomagnetic storm of 18 March and the solar eclipse of 20 March are marked with the hatching of the other style. The vertical bars indicate rms errors of the night OMR data.
Figure 9. The OMR C90 for the nights of 15–20 March 2015. Zero at the time scale corresponds to 00:00 MT 15 March. The dates are written above the timescale as 15.03, 16.03, …, 20.03. Night-time is highlighted with the hatching. Periods of the severe geomagnetic storm of 18 March and the solar eclipse of 20 March are marked with the hatching of the other style. The vertical bars indicate rms errors of the night OMR data.
Remotesensing 15 03440 g009
Figure 10. Calculated values of the OMR C90 (tE1, tE2) during the solar eclipses of 20 March 2015 (left panels (ad)) and 25 October 2022 (right panels (eh)). Rms errors σC90 of the calculated OMR values are written at every panel: (a) the eclipse of 20 March 2015, way (I); (b) the eclipse of 20 March 2015, way (II); (c) the eclipse of 20 March 2015, way (III); (d) the eclipse of 20 March 2015, way (IV); (e) the eclipse of 25 October 2015, way (I); (f) the eclipse of 25 October 2015, way (II); (g) the eclipse of 25 October 2015, way (III); (h) the eclipse of 25 October 2015, way (IV).
Figure 10. Calculated values of the OMR C90 (tE1, tE2) during the solar eclipses of 20 March 2015 (left panels (ad)) and 25 October 2022 (right panels (eh)). Rms errors σC90 of the calculated OMR values are written at every panel: (a) the eclipse of 20 March 2015, way (I); (b) the eclipse of 20 March 2015, way (II); (c) the eclipse of 20 March 2015, way (III); (d) the eclipse of 20 March 2015, way (IV); (e) the eclipse of 25 October 2015, way (I); (f) the eclipse of 25 October 2015, way (II); (g) the eclipse of 25 October 2015, way (III); (h) the eclipse of 25 October 2015, way (IV).
Remotesensing 15 03440 g010
Figure 11. Parameters BJ and DJ minimize the quadratic form (26). Blue lines refer to the solar eclipse of 20 March 2015, and red lines to the eclipse of 25 October 2022. The curves are marked with letters a–h, corresponding to the data shown in Figure 10 at panels notated with the same letters. The horizontal dashed line indicates BJ = 0.1 level.
Figure 11. Parameters BJ and DJ minimize the quadratic form (26). Blue lines refer to the solar eclipse of 20 March 2015, and red lines to the eclipse of 25 October 2022. The curves are marked with letters a–h, corresponding to the data shown in Figure 10 at panels notated with the same letters. The horizontal dashed line indicates BJ = 0.1 level.
Remotesensing 15 03440 g011
Table 1. Main parameters of the MOS-4.
Table 1. Main parameters of the MOS-4.
ParameterValue
Antenna pattern:
Beam width1.5° at −3 dB level
Side lobesLess than −30 dB
Antenna zenith angles38–100° nominal
50–65° available
Operating temperature295 K
Chopper wheel modulation frequency30 Hz
SSB noise temperature1100 K
SSB filter suppression for the image band 25–34 dB
PLL quarts oscillator stability±2.0 × 10−10 per day
±3 × 10−8 per year
±1.5 × 10−7 per decade
Acousto-optical spectrum analyzer (AOS):
Total band0.5 GHz
Number of channels1000
Frequency resolution1 MHz
Filter-bank spectrum analyzer (SA-50):
Total band10.8 MHz
Number of channels50
Frequency resolution0.1, 0.2, 0.4, 0.6 MHz
Table 2. Rate constants of the chemical reactions (12)–(16) [2].
Table 2. Rate constants of the chemical reactions (12)–(16) [2].
ReactionRate ConstantDimension
(12)2.15 × 10−34 e345/Tcm6 s−1
(13)2.15 × 10−34 e345/Tcm6 s−1
(14)8.82 × 10−35 e575/Tcm6 s−1
(15)1.5 × 10−11 e−2218/Tcm3 s−1
(16)1.4 × 10−10 e−270/Tcm3 s−1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rozanov, S.; Ignatyev, A.; Zavgorodniy, A. Ground-Based Microwave Measurements of Mesospheric Ozone Variations over Moscow Region during the Solar Eclipses of 20 March 2015 and 25 October 2022. Remote Sens. 2023, 15, 3440. https://doi.org/10.3390/rs15133440

AMA Style

Rozanov S, Ignatyev A, Zavgorodniy A. Ground-Based Microwave Measurements of Mesospheric Ozone Variations over Moscow Region during the Solar Eclipses of 20 March 2015 and 25 October 2022. Remote Sensing. 2023; 15(13):3440. https://doi.org/10.3390/rs15133440

Chicago/Turabian Style

Rozanov, Sergey, Alexander Ignatyev, and Alexey Zavgorodniy. 2023. "Ground-Based Microwave Measurements of Mesospheric Ozone Variations over Moscow Region during the Solar Eclipses of 20 March 2015 and 25 October 2022" Remote Sensing 15, no. 13: 3440. https://doi.org/10.3390/rs15133440

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop