Modulation of Polar Mesospheric Summer Echoes (PMSE) during HF Heating
- 1Department of Physics and Technology, UiT Arctic University of Norway, Tromsø, Norway
- 2Physics Faculty, TU Dortmund University, Germany
- 3Institude of Geophysics, University of Tehran, Iran
- 4Swedish Institute of Space Physics, IRF, Kiruna, Sweden
- 5EISCAT Scientific Association, Kiruna, Sweden
- 6EISCAT Scientific Association, Ramfjord, Norway
- 1Department of Physics and Technology, UiT Arctic University of Norway, Tromsø, Norway
- 2Physics Faculty, TU Dortmund University, Germany
- 3Institude of Geophysics, University of Tehran, Iran
- 4Swedish Institute of Space Physics, IRF, Kiruna, Sweden
- 5EISCAT Scientific Association, Kiruna, Sweden
- 6EISCAT Scientific Association, Ramfjord, Norway
Abstract. The formation of Polar Mesospheric Summer Echoes (PMSE) is linked to the presence of charged dust/ice particles in the mesosphere and the PMSE modulation is a possible way to observe the effect of dust charging. We investigate the modulation of PMSE by HF radio waves based on measurements carried out with the EISCAT Heating facility and the EISCAT VHF radar in August 2018 and August 2020 toward the end of the PMSE season. The measurements were made during the night with reduced solar illumination. The EISCAT Heating was operated in subsequent 48s on and 168s off intervals. In our observations, the PMSE modulation by the HF heating disappears during ionospheric conditions of energetic particle precipitation. We observe more than half of the cycles being influenced by the heating with a reduced PMSE power when the heater is on and a similar amount of cycles showing an increase in power when the heater is turned off. In less than half of the observed cycles we see an overshoot and it seems to be largely influenced by the overall PMSE power; with smaller or nonexistent overshoot when the PMSE power is high. We observe instances of very large overshoots, where background PMSE power seems to be reduced. During periods when we observe modulation, they often vary strongly from one cycle to the next; they are highly variable on spatial scales smaller km and time scales of minutes that are shorter than the scales assumed for the variation of dust parameters. Averaged curves over several heating cycles are similar to the overshoot curves predicted by theory and observed previously. Some individual curves however, show a stronger overshoot than observed in previous studies. A possible explanation for this difference can lie in the dust charging conditions that are different during the night or other conditions might be at play. We observe two possible instances of sporadic E-layers that are influenced by heating but do not show overshoots, as is to be expected.
Tinna Gunnarsdottir et al.
Status: closed
-
RC1: 'Comment on angeo-2022-16', Anonymous Referee #1, 20 Jun 2022
This paper provides an analysis of PMSE modulation by HF heating. It presents data from several days under varying background ionospheric conditions and analyses the response of the different PMSE layers to heating, including a comparison with simulation to explain the size of the overshoot.
There is some nice, thorough analysis and definitely merit in this work. However, there are also some major issues with the submission, mostly in the way that the paper is presented. The comments I provide below are given in the spirit of trying to improve the manuscript.
The main subject:
The first major issue is that it is not at all clear what the main point of the paper is. The title is very generic and gives no clue about what the new/noteworthy and exciting science is. Many papers have been written on HF modulation of PMSE (as the authors show in their introduction), why and what does this one add to the current state of knowledge?
I think that the main takeaway is that this is the first time that analysis has been performed on modulated PMSE under relatively low solar illumination, at late local times and at the end of the traditional PMSE season. The reference to this in the abstract (and how it may affect the observations) is very vague (line 15): “Some individual curves however, show a stronger overshoot than observed in previous studies. A possible explanation for this difference can lie in the dust charging conditions that are different during the night or other conditions might be at play.” Essentially you offer a suggestion and then undermine it by saying it might be something else, but without saying what the something else might be.
I highly recommend changing the title to better reflect the main aim/results of the paper (e.g., “Artificial modulation of PMSE during low solar illumination”). I recommend highlighting the unique aspects of the work in the abstract, introduction and discussion. Make it clear why this is different to previous work, or what it does to build on and support it.
In connection to this, the authors often digress into discussions of sporadic E (it is even mentioned in the abstract). This paper is about PMSE and I recommend removing the discussion and analysis of sporadic E since it is a distraction to the main topic (see line 173, lines 232 onwards). You highlight the E layers as an interval of interest in figure 3c and discuss it as Area 3 in Observations 3 (line 195). Please remove this, and simply state that “the layers above 90 km after 22:45 UT are not PMSE, they are examples of Sporadic E” perhaps with a reference. I understand the temptation to discuss interesting features, but it only detracts from the main focus; if there is something noteworthy about this interval, I recommend writing a separate article on it.
The style:
The manuscript reads often reads more like a student report rather than a journal article, a clear example of which is the inclusion of a table of contents. The submission includes a sizeable appendix, which is effectively the same size as the main body of the article. I leave it to the editor to determine whether the appendix is appropriate.
The submission requires significant proof reading for some spelling but mostly for grammar, with attention on the appropriate use of articles (‘the’, ‘a’). I do recognise that English is not the easiest language when it comes to this issue and that it can be far from intuitive, but I note that the paper has experienced authors and at least one who has English as a first language and so I am sure they can be helpful.
Additional Comments
Line 61 you give a range of temperatures (130 – 150 K) and then say, ‘and lower’. Why not say “150 K and lower”.
Line 69: You state that “The charges and sizes of the dust particles are largely unknown”. Is this really true or is it more nuanced than that? Lots of work has been done by Havnes (for example - Havnes, O. and Kassa, M.: On the sizes and observable effects of dust particles in polar mesospheric winter echoes, J. Geophys. Res, 114, D09209, doi:10.1029/2008JD011276, 2009.) on examining dust sizes, others have also done a lot (e.g., Rapp, Scales, etc)
Line 71: ‘shown’ instead of ‘suggested’
Line 130 onwards: you describe the observations taking place under dusk and night conditions. This gives the initial impression of some of the observations taking place in darkness, but are they really? You go onto discuss the reduction in solar illumination of an order of magnitude compared to midday in mid-Summer, but I am not sure if a reader would have a feel for what this really means. To put it in context, how does the range of illumination vary through the year (lowest to highest) – where on the range do these observations sit? Did you consider including a sentence on the solar zenith angle (just the average if it doesn’t vary much) at the times of the observations? This is a measure that I suspect the readers may be more familiar with and are more likely to intuitively understand the difference in degrees between mid-summer and these observations.
Given that the level of illumination is potentially the major selling point of the paper, it would be wise to make sure the reader understood how unique this might be in terms of past observations.
Line 138: Is it possible to give an estimate of the change in the photo-emission current? Or indicate that is something that you will consider in the analysis. Saying it ‘should translate…’ without providing more detail seems very vague at this stage.
Line 139 – 144: this section seems to describe observations that were mostly irrelevant. Except for the later comment about lack of NLC limiting likelihood of dust size I would have recommended removing this. This could be shortened to: “weather conditions were poor for measuring NLC, though cameras were operating at Kiruna and Nikkaluokta (~200km south of Tromso) to allow triangulation. During 15/16 August, faint NLC were visible from Kiruna, close to the horizon, but westward of the EISCAT site (closer to Andoya)”
I do have a question about this: how confident are you that the lack of NLC above Tromso was not due to low clouds in the line of site of the camera? This is quite important since your dismissal of larger particle sizes later on relies on this.
Figure 2 and lines 145 onwards: these measurements have been selected as they are closest to the radar measurements at the time of the PMSE observations. This is unacceptably vague, how far away were they? 1 degree of latitude, 2, 10? It is important to give the right context for the reader to judge.
Lines 153-160: I thought it might be useful to include an example of some of the grammatical issues that occur throughout the paper, along with an indication of some alternate wording and removal of over-precise measurements.
“The radar observations were made in the zenith direction with the EISCAT VHF (224 MHz) antennas radar near Tromsø (69.59°N, 19.23°E). The VHF and Heating parameters are given in Table 1. Details on the radar experiments were provided by the EISCAT documentation (Tjulin, 2022*). EISCAT Heating (Rietveld et al., 1993, 2016) was operated with a vertical beam at 5.423 MHz with a nominal 80 kW per transmitter, which corresponds to Effective Radiated Power (ERP) between 500 and 580 MW. and X-mode polarization was used with HF pumping in a sequence of 48 s on and 168 s off. The beam width of the Heating facility (7°) is much larger than that of the VHF radar (~4°) and given that the vertical winds and velocity fluctuations of the PMSE are only a few m/s and horizontal winds possibly upto a few 10s m/s (Strelnikova and Rapp, 2011), the radar is always measuring PMSE influenced by the heating.”
*include the documentation as a reference, but make sure you pick the relevant edition for the year of the experiments.
Line 167: remove “we consider the back-scattered power in units…” since you have effectively said this on line 163.
Line 185 onwards: I am curious about the early part of this observation, it is very hard to see because of the scales, but even the weak precipitation from 2000 – 2110 UT seem to shows signs of modulation at higher altitudes. I know that precipitation can show distinct natural modulation that leads to regular Ne patterns, but I wonder if this is a problem with the analysis not fully taking into account that the system power is fluctuating as the heater is turned on and off.
I would also be wary about suggesting the PMSE in (1) ‘seems that it was trigged by the heating’, beyond the coincidence of timing is there a reason for this speculation?
Power modulation might also be visible in Observation 3 around 21 UT. Is it possible that this might also be responsible for the modulation of the sporadic E later in that interval, there are hints the striations continue to higher altitudes? I recommend talking to your co-authors Reitveld and Haggstrom to double check whether this might be the case.
Line 220 onwards: if these figures are important enough to discuss in the main text, why are they not included in the main body of the paper?
Line 223: why was this value selected? You discuss the threshold is perhaps too high, in which case why did you not try the analysis with a different threshold and publish that?
Line 262 and elsewhere: when you refer to particular data that you have presented, you should refer to that figure (i.e. fig3c)
Lines 388 onwards: I found the description a bit confusing. When you say that the ‘simulations assume an initial plasma temp…’ do you mean that you set that as the initial temperature based on the IRI value for that time? Or did you mean that the simulations are hard wired with that temperature which thankfully happened to match what IRI would expect? I assume not the latter but that is sort of how it reads. Also, how confident are you that the IRI temperature is itself reasonable. Did you look at the EISCAT temperatures for heater-off times around the PMSE region to get an estimate for comparison?
How did you determine that 3nm dust particles were the best fit – did you judge it by eye, or did you fit the model to the data by altering the dust size until some measure of best fit was achieved?
Figure 16. Why is the X-axis scale different between the data and the simulation?
Line 472 onwards: is there a way to explore how much of an effect absorption of the heater beam might play? The way this is written you seem to be suggesting that you have an answer for the large overshoot in interval A, but that other effects may be important and perhaps even fully explain the large overshoot, thus negating your initial statement. It is really good to consider other factors and suggest alternative explanations, however it is better if you can provide an indication of how confident you are in your explanation and why. As it stands it sort of reads as if you are saying: “analysis suggests that X can explain the overshoot effect, but it could also be explained by Y instead”. Stating uncertainties and explaining the limits of analysis are of course very good, but at the same time you can make it clearer how confident you are in the analysis.
“we know that Y potentially contributes to the effect at an unquantified level; however our analysis shows that N% of the effect is explained by X”
I was surprised that there is no reference to:
Havnes, O., Pinedo, H., La Hoz, C., Senior, A., Hartquist, T. W., Rietveld, M. T., and Kosch, M. J.: A comparison of overshoot modelling with observations of polar mesospheric summer echoes at radar frequencies of 56 and 224 MHz, Ann. Geophys., 33, 737–747, https://doi.org/10.5194/angeo-33-737-2015, 2015.
who also found varying overshoot for different conditions (albeit with different radar frequencies as well).
Summary
There is some good analysis and interesting data in this submission. Some work is required to bring the level of presentation up to standard and in doing so highlight the important aspects of the work. The title needs rethinking, and the key results need to be better highlighted in the abstract, introduction and discussion to make them stand out more. I think its important to check that the guisdap analysis has properly taken the system porwer changes into account in each instance.
I am confident that the authors can address these issues and produce a good paper.
-
AC1: 'Reply on RC1', Tinna Gunnarsdottir, 19 Sep 2022
Dear Anonymous Referee #1
We the authors thank you for the time you have taken to read and review our preprint as well as your detailed comments and suggestions. We agree that the paper and title needs to be better focused on the main point of the paper. In addition, the overall language/grammar will be revised. We will remove the sporadic E discussion since it is not relevant to the paper as well as the table of contents and shorten the appendix where possible.
Comments:
Line 61: Yes that is better.
Line 69: The sizes and possibly also the charge states are possibly different between summer and winter echoes but yes you are right, this sentence makes it sound like there is no previous work done when that is not the case. We will change this to better reflect previous work and current understanding on the subject.
Line 130: Yes thank you for highlighting this point, we will make sure that the modified version will include a more detailed description of the solar illumination and the UV flux for these sets of observations. The solar zenith angles are all in the range 88-97 degrees with the main part of the observations larger than 90 degrees. This would mean that the sun is below the horizon but only by a few degrees.
Line 138: Yes we can consider this for the analysis and give relevant estimates.
Line 139-144: We can keep this on the shorter side. Yes you are right, we cannot conclude that there were no NLC present for the observations with bad weather. For the 15-16 of august 2018 observation however, there were no NLC visible above the EISCAT site.
Figure 2 and lines 145+: For Fig.2a, the horizontal distance between these profiles and Tromsø is 490 km. For Fig. 2b, the horizontal distance is 293 km. The geographical coordinates are in the figures but we forgot to specify this. We can include this information in the text for the reader.
Lines 153-160: Thank you for your example, we will review the text as a whole and make the necessary changes. In addition, add in the reference of the EISCAT documentation.
Line 167: Yes we can remove this.
Line 185+: We will investigate if it is possible that the system power could have influenced the analysis for this observation. As the modulation is visible only in this area of low electron precipitation and not at others we assumed that this is due to the heater modulating the higher electron density in this altitude range.
Our suggestion here was that if there would not have been heating there most likely would not have been PMSE visible in area 1 and that the overshoot caused by the heating make the faint PMSE visible.
There were some operational problems in the early part of observation 3 that should be fixed before start of area 1. We will investigate if there was any lingering problems that caused the modulation of the sporadic E and could have influenced area 1 and 2.
Line 220+: Yes we can add them to the main body.
Line 223: We use this threshold since other authors have used it to indicate the presence of PMSE and thus we can better compare our results to theirs. In addition, by lowering the threshold we would introduce many instances of data were there is no PMSE present.
Line 262 and elsewhere: Yes we can include this.
Lines 388 onwards: The simulations are set at 150 K and the heating temperature calculated as a multiple of that initial temperature. By using 150 K the simulations can be compared to previous publications that have simulated the overshoot since many use 150 K as an initial temperature. This value was also quite close to both the IRI model values as well as the satellite temperature measurements from figure 2. We can look into using EISCAT data for estimating the temperature around the PMSE layers.
The simulation was run for several dust sizes and 3 nm deemed the best match – judged by eye.
Figure 16: This is just to show that the simulation goes back to equilibrium when you let it run longer. This is also what we expect from the simulations. For interval A the simulation suggest that this takes longer than what the observations suggests.
Line 472+: Yes you are right, we will review our statements here and make it more clear how confident we are in the analysis and determine if absorption can have a determinable effect.
Thank you for this the suggestion to the Havnes et al (2015) paper; this is a very interesting paper with relevant aspects. We will make sure to include it in our discussion.
-
AC1: 'Reply on RC1', Tinna Gunnarsdottir, 19 Sep 2022
-
RC2: 'Comment on angeo-2022-16', Stephan C. Buchert, 21 Jun 2022
The authors display various aspects of observed PMSE under articial HF heating and stretching over ~4-6 hours each of several nights and a thorough analysis of the data. A simulation complements this analysis. Generally I also think that the main points do not become immediately obvious to the reader, owing to the large amount of presented material and the sometimes lengthy and too pondering language. Appendices show material that the authors consider as details and perhaps less important, which I think is a good division. Nevertheless, some "streamlining" of the opus could improve it.
Specifically,
Figure 3c and line 198: "Area 1" starts at 21:30 UT, but the PMSE seems be present already at the start of the measurements at 20:30 UT, partially with precipitation. What is the reason that this period was not included in "Area 1" or added as separate "Area"? The period would be perhaps be long enough to make a small difference for some of the statistical analysis.
Line 454: "... λirreg is the wavelength of the irregularity and is given by 2 ∗ λR (the radar wavelength) ..." According to the Bragg condition (https://en.wikipedia.org/wiki/Bragg%27s_law) the irregularity wavelength would be half of radar wavelength , i. e. "given by λR/2", just 67 cm for the EISCAT VHF.
Is this just a textual mistake, or are the columns "τdff at 0 sec" and "τdff at 48 sec" of Table 3, I think, too large by a factor 4? According to my rough calculation the latter is the case.
Which values of Te/Ti and Zd*nd0/ne0 in equation (1) were used to produce the numbers is Table 3?
Te/Ti: Were the temperatures obtained from the GUISDAP analysis? In the standard configuration GUISDAP would always give values Te=Ti≈150 K. However, the HF heating increases Te dramatically, which is the main cause of the suppression of PMSE at "heater on" and the overshoot following "heater off". This increase of Te probably cannot be estimated from the radar data. Chen and Scales (2005) had assumed that Te/Te0 (=Te/Ti) would reach 10. But the values in Table 3, ratio of columns "τdff at 0 sec" and "τdff at 48 sec", seem to be inconsistent with such large heating. A plot of the for the simulation assumed Te/Te0 over time might be helpful.
Similarly, which model of dust charging was simulated/assumed in the simulations? The numbers in Table 3 suggest that Zd*nd0/ne0 did not change much between 0 and 48 sec, but the overshoot would depend strongly on the amount of dust charging?
These minor issues perhaps don't affect the main conclusions of the paper, but should of course be checked and explained. Perhaps assuming a larger Te/Te0 and stronger dust charging could explain the large observed overshoots?
I tried to verify some of these details by looking at the code, but the link given on line 529 does not work for me.
-
AC2: 'Reply on RC2', Tinna Gunnarsdottir, 19 Sep 2022
We the authors thank you Stephan C. Buchert for your comments and help in making our manuscript better. We will review the paper and try to make it more straightforward and to the point.
Figure 3c and line 198: Unfortunately due to operation problems with the heater this part of the observation could not be included in the analysis.
Line 454: Yes you are right, it is λirreg = λR/2. An unfortunate mistake that has caused the calculated diffusion timescales to be much too large.
No the temperature ratio is taken from the simulation and not the radar data. Values for Te/Ti at 0 seconds is assumed to be 1 before the heater is turned on, for values at 48 seconds we use the ratios given in figure 16 (simulation on right). Te/Ti = 8 for interval A and Te/Ti = 5 for intervals B and C. Values for Zd come from figure 17 (left) taken at 0 seconds and 48 seconds and nd0/ne is taken from given values in figure 16(right).
For the estimation of tdff in table 3 only Zd is assumed to change from 0 sec to 48 sec not the ratio nd0/ne0 where we have used the value for Zd from figure 17. The ratio for nd0/ne0 we have used the values from figure 16(right) and assumed that it is the same for both 0 sec and 48 sec.
We have decided to remove table 3 from the manuscript, as its main purpose was to examine the neutral density effect of the diffusion timescale however the diffusion time scale should be calculated using numerical results and not estimated in a few instances as we have done here. We will rather focus on the figures that show the ratio of the two timescales connected to the simulation and only mention involvement of neutral density.
We apologize for the link not working we will make sure to fix this.
Thank you for your discussion and comments.
-
AC2: 'Reply on RC2', Tinna Gunnarsdottir, 19 Sep 2022
Status: closed
-
RC1: 'Comment on angeo-2022-16', Anonymous Referee #1, 20 Jun 2022
This paper provides an analysis of PMSE modulation by HF heating. It presents data from several days under varying background ionospheric conditions and analyses the response of the different PMSE layers to heating, including a comparison with simulation to explain the size of the overshoot.
There is some nice, thorough analysis and definitely merit in this work. However, there are also some major issues with the submission, mostly in the way that the paper is presented. The comments I provide below are given in the spirit of trying to improve the manuscript.
The main subject:
The first major issue is that it is not at all clear what the main point of the paper is. The title is very generic and gives no clue about what the new/noteworthy and exciting science is. Many papers have been written on HF modulation of PMSE (as the authors show in their introduction), why and what does this one add to the current state of knowledge?
I think that the main takeaway is that this is the first time that analysis has been performed on modulated PMSE under relatively low solar illumination, at late local times and at the end of the traditional PMSE season. The reference to this in the abstract (and how it may affect the observations) is very vague (line 15): “Some individual curves however, show a stronger overshoot than observed in previous studies. A possible explanation for this difference can lie in the dust charging conditions that are different during the night or other conditions might be at play.” Essentially you offer a suggestion and then undermine it by saying it might be something else, but without saying what the something else might be.
I highly recommend changing the title to better reflect the main aim/results of the paper (e.g., “Artificial modulation of PMSE during low solar illumination”). I recommend highlighting the unique aspects of the work in the abstract, introduction and discussion. Make it clear why this is different to previous work, or what it does to build on and support it.
In connection to this, the authors often digress into discussions of sporadic E (it is even mentioned in the abstract). This paper is about PMSE and I recommend removing the discussion and analysis of sporadic E since it is a distraction to the main topic (see line 173, lines 232 onwards). You highlight the E layers as an interval of interest in figure 3c and discuss it as Area 3 in Observations 3 (line 195). Please remove this, and simply state that “the layers above 90 km after 22:45 UT are not PMSE, they are examples of Sporadic E” perhaps with a reference. I understand the temptation to discuss interesting features, but it only detracts from the main focus; if there is something noteworthy about this interval, I recommend writing a separate article on it.
The style:
The manuscript reads often reads more like a student report rather than a journal article, a clear example of which is the inclusion of a table of contents. The submission includes a sizeable appendix, which is effectively the same size as the main body of the article. I leave it to the editor to determine whether the appendix is appropriate.
The submission requires significant proof reading for some spelling but mostly for grammar, with attention on the appropriate use of articles (‘the’, ‘a’). I do recognise that English is not the easiest language when it comes to this issue and that it can be far from intuitive, but I note that the paper has experienced authors and at least one who has English as a first language and so I am sure they can be helpful.
Additional Comments
Line 61 you give a range of temperatures (130 – 150 K) and then say, ‘and lower’. Why not say “150 K and lower”.
Line 69: You state that “The charges and sizes of the dust particles are largely unknown”. Is this really true or is it more nuanced than that? Lots of work has been done by Havnes (for example - Havnes, O. and Kassa, M.: On the sizes and observable effects of dust particles in polar mesospheric winter echoes, J. Geophys. Res, 114, D09209, doi:10.1029/2008JD011276, 2009.) on examining dust sizes, others have also done a lot (e.g., Rapp, Scales, etc)
Line 71: ‘shown’ instead of ‘suggested’
Line 130 onwards: you describe the observations taking place under dusk and night conditions. This gives the initial impression of some of the observations taking place in darkness, but are they really? You go onto discuss the reduction in solar illumination of an order of magnitude compared to midday in mid-Summer, but I am not sure if a reader would have a feel for what this really means. To put it in context, how does the range of illumination vary through the year (lowest to highest) – where on the range do these observations sit? Did you consider including a sentence on the solar zenith angle (just the average if it doesn’t vary much) at the times of the observations? This is a measure that I suspect the readers may be more familiar with and are more likely to intuitively understand the difference in degrees between mid-summer and these observations.
Given that the level of illumination is potentially the major selling point of the paper, it would be wise to make sure the reader understood how unique this might be in terms of past observations.
Line 138: Is it possible to give an estimate of the change in the photo-emission current? Or indicate that is something that you will consider in the analysis. Saying it ‘should translate…’ without providing more detail seems very vague at this stage.
Line 139 – 144: this section seems to describe observations that were mostly irrelevant. Except for the later comment about lack of NLC limiting likelihood of dust size I would have recommended removing this. This could be shortened to: “weather conditions were poor for measuring NLC, though cameras were operating at Kiruna and Nikkaluokta (~200km south of Tromso) to allow triangulation. During 15/16 August, faint NLC were visible from Kiruna, close to the horizon, but westward of the EISCAT site (closer to Andoya)”
I do have a question about this: how confident are you that the lack of NLC above Tromso was not due to low clouds in the line of site of the camera? This is quite important since your dismissal of larger particle sizes later on relies on this.
Figure 2 and lines 145 onwards: these measurements have been selected as they are closest to the radar measurements at the time of the PMSE observations. This is unacceptably vague, how far away were they? 1 degree of latitude, 2, 10? It is important to give the right context for the reader to judge.
Lines 153-160: I thought it might be useful to include an example of some of the grammatical issues that occur throughout the paper, along with an indication of some alternate wording and removal of over-precise measurements.
“The radar observations were made in the zenith direction with the EISCAT VHF (224 MHz) antennas radar near Tromsø (69.59°N, 19.23°E). The VHF and Heating parameters are given in Table 1. Details on the radar experiments were provided by the EISCAT documentation (Tjulin, 2022*). EISCAT Heating (Rietveld et al., 1993, 2016) was operated with a vertical beam at 5.423 MHz with a nominal 80 kW per transmitter, which corresponds to Effective Radiated Power (ERP) between 500 and 580 MW. and X-mode polarization was used with HF pumping in a sequence of 48 s on and 168 s off. The beam width of the Heating facility (7°) is much larger than that of the VHF radar (~4°) and given that the vertical winds and velocity fluctuations of the PMSE are only a few m/s and horizontal winds possibly upto a few 10s m/s (Strelnikova and Rapp, 2011), the radar is always measuring PMSE influenced by the heating.”
*include the documentation as a reference, but make sure you pick the relevant edition for the year of the experiments.
Line 167: remove “we consider the back-scattered power in units…” since you have effectively said this on line 163.
Line 185 onwards: I am curious about the early part of this observation, it is very hard to see because of the scales, but even the weak precipitation from 2000 – 2110 UT seem to shows signs of modulation at higher altitudes. I know that precipitation can show distinct natural modulation that leads to regular Ne patterns, but I wonder if this is a problem with the analysis not fully taking into account that the system power is fluctuating as the heater is turned on and off.
I would also be wary about suggesting the PMSE in (1) ‘seems that it was trigged by the heating’, beyond the coincidence of timing is there a reason for this speculation?
Power modulation might also be visible in Observation 3 around 21 UT. Is it possible that this might also be responsible for the modulation of the sporadic E later in that interval, there are hints the striations continue to higher altitudes? I recommend talking to your co-authors Reitveld and Haggstrom to double check whether this might be the case.
Line 220 onwards: if these figures are important enough to discuss in the main text, why are they not included in the main body of the paper?
Line 223: why was this value selected? You discuss the threshold is perhaps too high, in which case why did you not try the analysis with a different threshold and publish that?
Line 262 and elsewhere: when you refer to particular data that you have presented, you should refer to that figure (i.e. fig3c)
Lines 388 onwards: I found the description a bit confusing. When you say that the ‘simulations assume an initial plasma temp…’ do you mean that you set that as the initial temperature based on the IRI value for that time? Or did you mean that the simulations are hard wired with that temperature which thankfully happened to match what IRI would expect? I assume not the latter but that is sort of how it reads. Also, how confident are you that the IRI temperature is itself reasonable. Did you look at the EISCAT temperatures for heater-off times around the PMSE region to get an estimate for comparison?
How did you determine that 3nm dust particles were the best fit – did you judge it by eye, or did you fit the model to the data by altering the dust size until some measure of best fit was achieved?
Figure 16. Why is the X-axis scale different between the data and the simulation?
Line 472 onwards: is there a way to explore how much of an effect absorption of the heater beam might play? The way this is written you seem to be suggesting that you have an answer for the large overshoot in interval A, but that other effects may be important and perhaps even fully explain the large overshoot, thus negating your initial statement. It is really good to consider other factors and suggest alternative explanations, however it is better if you can provide an indication of how confident you are in your explanation and why. As it stands it sort of reads as if you are saying: “analysis suggests that X can explain the overshoot effect, but it could also be explained by Y instead”. Stating uncertainties and explaining the limits of analysis are of course very good, but at the same time you can make it clearer how confident you are in the analysis.
“we know that Y potentially contributes to the effect at an unquantified level; however our analysis shows that N% of the effect is explained by X”
I was surprised that there is no reference to:
Havnes, O., Pinedo, H., La Hoz, C., Senior, A., Hartquist, T. W., Rietveld, M. T., and Kosch, M. J.: A comparison of overshoot modelling with observations of polar mesospheric summer echoes at radar frequencies of 56 and 224 MHz, Ann. Geophys., 33, 737–747, https://doi.org/10.5194/angeo-33-737-2015, 2015.
who also found varying overshoot for different conditions (albeit with different radar frequencies as well).
Summary
There is some good analysis and interesting data in this submission. Some work is required to bring the level of presentation up to standard and in doing so highlight the important aspects of the work. The title needs rethinking, and the key results need to be better highlighted in the abstract, introduction and discussion to make them stand out more. I think its important to check that the guisdap analysis has properly taken the system porwer changes into account in each instance.
I am confident that the authors can address these issues and produce a good paper.
-
AC1: 'Reply on RC1', Tinna Gunnarsdottir, 19 Sep 2022
Dear Anonymous Referee #1
We the authors thank you for the time you have taken to read and review our preprint as well as your detailed comments and suggestions. We agree that the paper and title needs to be better focused on the main point of the paper. In addition, the overall language/grammar will be revised. We will remove the sporadic E discussion since it is not relevant to the paper as well as the table of contents and shorten the appendix where possible.
Comments:
Line 61: Yes that is better.
Line 69: The sizes and possibly also the charge states are possibly different between summer and winter echoes but yes you are right, this sentence makes it sound like there is no previous work done when that is not the case. We will change this to better reflect previous work and current understanding on the subject.
Line 130: Yes thank you for highlighting this point, we will make sure that the modified version will include a more detailed description of the solar illumination and the UV flux for these sets of observations. The solar zenith angles are all in the range 88-97 degrees with the main part of the observations larger than 90 degrees. This would mean that the sun is below the horizon but only by a few degrees.
Line 138: Yes we can consider this for the analysis and give relevant estimates.
Line 139-144: We can keep this on the shorter side. Yes you are right, we cannot conclude that there were no NLC present for the observations with bad weather. For the 15-16 of august 2018 observation however, there were no NLC visible above the EISCAT site.
Figure 2 and lines 145+: For Fig.2a, the horizontal distance between these profiles and Tromsø is 490 km. For Fig. 2b, the horizontal distance is 293 km. The geographical coordinates are in the figures but we forgot to specify this. We can include this information in the text for the reader.
Lines 153-160: Thank you for your example, we will review the text as a whole and make the necessary changes. In addition, add in the reference of the EISCAT documentation.
Line 167: Yes we can remove this.
Line 185+: We will investigate if it is possible that the system power could have influenced the analysis for this observation. As the modulation is visible only in this area of low electron precipitation and not at others we assumed that this is due to the heater modulating the higher electron density in this altitude range.
Our suggestion here was that if there would not have been heating there most likely would not have been PMSE visible in area 1 and that the overshoot caused by the heating make the faint PMSE visible.
There were some operational problems in the early part of observation 3 that should be fixed before start of area 1. We will investigate if there was any lingering problems that caused the modulation of the sporadic E and could have influenced area 1 and 2.
Line 220+: Yes we can add them to the main body.
Line 223: We use this threshold since other authors have used it to indicate the presence of PMSE and thus we can better compare our results to theirs. In addition, by lowering the threshold we would introduce many instances of data were there is no PMSE present.
Line 262 and elsewhere: Yes we can include this.
Lines 388 onwards: The simulations are set at 150 K and the heating temperature calculated as a multiple of that initial temperature. By using 150 K the simulations can be compared to previous publications that have simulated the overshoot since many use 150 K as an initial temperature. This value was also quite close to both the IRI model values as well as the satellite temperature measurements from figure 2. We can look into using EISCAT data for estimating the temperature around the PMSE layers.
The simulation was run for several dust sizes and 3 nm deemed the best match – judged by eye.
Figure 16: This is just to show that the simulation goes back to equilibrium when you let it run longer. This is also what we expect from the simulations. For interval A the simulation suggest that this takes longer than what the observations suggests.
Line 472+: Yes you are right, we will review our statements here and make it more clear how confident we are in the analysis and determine if absorption can have a determinable effect.
Thank you for this the suggestion to the Havnes et al (2015) paper; this is a very interesting paper with relevant aspects. We will make sure to include it in our discussion.
-
AC1: 'Reply on RC1', Tinna Gunnarsdottir, 19 Sep 2022
-
RC2: 'Comment on angeo-2022-16', Stephan C. Buchert, 21 Jun 2022
The authors display various aspects of observed PMSE under articial HF heating and stretching over ~4-6 hours each of several nights and a thorough analysis of the data. A simulation complements this analysis. Generally I also think that the main points do not become immediately obvious to the reader, owing to the large amount of presented material and the sometimes lengthy and too pondering language. Appendices show material that the authors consider as details and perhaps less important, which I think is a good division. Nevertheless, some "streamlining" of the opus could improve it.
Specifically,
Figure 3c and line 198: "Area 1" starts at 21:30 UT, but the PMSE seems be present already at the start of the measurements at 20:30 UT, partially with precipitation. What is the reason that this period was not included in "Area 1" or added as separate "Area"? The period would be perhaps be long enough to make a small difference for some of the statistical analysis.
Line 454: "... λirreg is the wavelength of the irregularity and is given by 2 ∗ λR (the radar wavelength) ..." According to the Bragg condition (https://en.wikipedia.org/wiki/Bragg%27s_law) the irregularity wavelength would be half of radar wavelength , i. e. "given by λR/2", just 67 cm for the EISCAT VHF.
Is this just a textual mistake, or are the columns "τdff at 0 sec" and "τdff at 48 sec" of Table 3, I think, too large by a factor 4? According to my rough calculation the latter is the case.
Which values of Te/Ti and Zd*nd0/ne0 in equation (1) were used to produce the numbers is Table 3?
Te/Ti: Were the temperatures obtained from the GUISDAP analysis? In the standard configuration GUISDAP would always give values Te=Ti≈150 K. However, the HF heating increases Te dramatically, which is the main cause of the suppression of PMSE at "heater on" and the overshoot following "heater off". This increase of Te probably cannot be estimated from the radar data. Chen and Scales (2005) had assumed that Te/Te0 (=Te/Ti) would reach 10. But the values in Table 3, ratio of columns "τdff at 0 sec" and "τdff at 48 sec", seem to be inconsistent with such large heating. A plot of the for the simulation assumed Te/Te0 over time might be helpful.
Similarly, which model of dust charging was simulated/assumed in the simulations? The numbers in Table 3 suggest that Zd*nd0/ne0 did not change much between 0 and 48 sec, but the overshoot would depend strongly on the amount of dust charging?
These minor issues perhaps don't affect the main conclusions of the paper, but should of course be checked and explained. Perhaps assuming a larger Te/Te0 and stronger dust charging could explain the large observed overshoots?
I tried to verify some of these details by looking at the code, but the link given on line 529 does not work for me.
-
AC2: 'Reply on RC2', Tinna Gunnarsdottir, 19 Sep 2022
We the authors thank you Stephan C. Buchert for your comments and help in making our manuscript better. We will review the paper and try to make it more straightforward and to the point.
Figure 3c and line 198: Unfortunately due to operation problems with the heater this part of the observation could not be included in the analysis.
Line 454: Yes you are right, it is λirreg = λR/2. An unfortunate mistake that has caused the calculated diffusion timescales to be much too large.
No the temperature ratio is taken from the simulation and not the radar data. Values for Te/Ti at 0 seconds is assumed to be 1 before the heater is turned on, for values at 48 seconds we use the ratios given in figure 16 (simulation on right). Te/Ti = 8 for interval A and Te/Ti = 5 for intervals B and C. Values for Zd come from figure 17 (left) taken at 0 seconds and 48 seconds and nd0/ne is taken from given values in figure 16(right).
For the estimation of tdff in table 3 only Zd is assumed to change from 0 sec to 48 sec not the ratio nd0/ne0 where we have used the value for Zd from figure 17. The ratio for nd0/ne0 we have used the values from figure 16(right) and assumed that it is the same for both 0 sec and 48 sec.
We have decided to remove table 3 from the manuscript, as its main purpose was to examine the neutral density effect of the diffusion timescale however the diffusion time scale should be calculated using numerical results and not estimated in a few instances as we have done here. We will rather focus on the figures that show the ratio of the two timescales connected to the simulation and only mention involvement of neutral density.
We apologize for the link not working we will make sure to fix this.
Thank you for your discussion and comments.
-
AC2: 'Reply on RC2', Tinna Gunnarsdottir, 19 Sep 2022
Tinna Gunnarsdottir et al.
Tinna Gunnarsdottir et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
553 | 79 | 17 | 649 | 7 | 6 |
- HTML: 553
- PDF: 79
- XML: 17
- Total: 649
- BibTeX: 7
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1