Despite efforts to make the transition transparent, some code in some calculators turned out to be less reliable on the new site and some changes made to make them reliable. QSRF was the main troublesome app, but is now stable and the update service has been improved to two hourly in the revision.

It is a quite complex web site, and a huge amount of work, some 10,000+ files, several databases, lots of custom PHP applications, and subtle but significant changes in PHP and apache environments.

Hopefully it has been worth the work (~150h) in terms of reliable affordable service for the future. The new hosting environment is certainly easier to work with and properly supports remote secured access… so that has to be good.

]]>With a change in hosting and the changes to adapt to supported and unsupported features, one of the positives is that the ‘short term’ fetching of data on which Quiet solar radio flux interpolations calculator is based is now happening two hourly at 5 minutes past the even hours.

If you wish, you could email SWPC.Webmaster@noaa.gov thanking them for their service and encouraging them to restore the rad.txt and 45day_rad.txt files to HTTP access quickly.

The data remains available on an FTP server, but my web server cannot access it for security reasons, and it would be unacceptably slow anyway.

I wrote to NOAA and they replied promptly:

You wrote to request restoring the 45day_rad.text and rad.txt to be accessible via HTTP. We are working to make a number of our products currently available only via FTP also available via our services at http://services.swpc.noaa.gov/text/. The files you’ve requested will be among those, but I do not yet have a schedule for the update.

As you can see, there is no time given for what ought to be a fairly trivial fix to restore availability under HTTP, albeit on a different URL.

In the mean time, I have implemented a temporary measure to take a snapshot of the FTP server’s file once per day (midday Sydney time) which should be up to date for the previous UTC day.

If you wish, you could email SWPC.Webmaster@noaa.gov thanking them for their service and encouraging them to restore the rad.txt and 45day_rad.txt files to HTTP access quickly.

The factor I is λ^2*Φ/(8*π*k) where λ is wavelength of the measurement, Φ is solar flux at **that** wavelength, k is Boltzmann constant.

The authors give a set of approximations for I in terms of solar flux at 10.7cm, based they say on a polynomial curve fit on experimental data. They do not give further information on their approximations.

Above is a scatter diagram of observations at the Palehua observatory during 2011 of solar flux at 2695MHz and 410MHz.

While these frequencies are not exactly those used in (Bertelsmeier & Magnin 1992), they are in almost the same ratio as 2800MHz and 432MHz (their formula 3), but displaced about 4% lower in frequency. They probably have very similar characteristics to 2800MHz vs 432MHz.

Examining the chart, there is obviously some relationship between the values, but the question are what is it, and more importantly, how strong is that relationship.

Fitting a simple linear relationship gives the formula SF410=0.1431*SF2694+19.936 as shown on the chart, but the more interesting statistic is the correlation coefficient R^2=0.1954. It is generally accepted that strong relationships are indicated by R^2>0.8, and that below that, any relationship is weak. The R^2 for this curve fit at 0.1954 is stunningly low and indicates to me that there will be large uncertainty in using the linear model.

Nevertheless, I have calculated the RMS error in using the model to predict the 341 data points in the original dataset, and it is 16.7SFU. This is an appalling error which in terms of calculated flux around 40SFU is an error of 3.8dB at 1σ confidence interval, even worse at higher confidence.

Let us compare the G/T calculated by (Bertelsmeier & Magnin 1992) for the most recent solar observations (27/04/14) given a Y factor of 6dB.

Above is the raw data from NOAA.

Above are the interpolated flux values at 432 and 1296MHz using Quiet sun radio flux interpolations.

Above the results of calculation from the same set of solar observations using the methods of (Bertelsmeier & Magnin 1992) (B&M) and (Duffy 2007),(Duffy 2009), (Duffy 2014), (Duffy 2014b) (D). There is around 4dB difference, quite substantial.

On the basis of this analysis, I would expect (Bertelsmeier & Magnin 1992) to be subject to similar uncertainty, and the uncertainty is too large to be used to verify small improvements or otherwise in a system.

Their approximations may account apparent errors in many analyses done by Hams, and may be embedded in calculation tools that users apparently do not understand and have not verified.

- Duffy, O. 2007. Measuring system G/T ratio using Sun noise. VK1OD.net (offline).
- ———. 2009. Quiet sun radio flux interpolations. https://owenduffy.net/calc/qsrf/index.htm.
- ———. 2014. Measuring G/T. https://owenduffy.net/blog/?p=1490.
- ———. 2014b. Calculate G/T from Sun noise Y measurement. https://owenduffy.net/calc/sf2gt.htm.
- Bertelsmeier, R and Magnin, P. 1992. Performance evaluation for EME systems In DUBUS 3/92. http://dpmc.unige.ch/dubus/9203-3.pdf(accessed 29/04/14).
- ITU-R. 2000. Recommendation ITU-R S.733-2 (2000)Determination of the G/T ratio for earth stations operating in the fixed-satellite service .

The first system he describes is 4 x ZB7013 Yagis with claimed gain of 22.65dB, and measured Sun/ColdSky of 4.5dB at SF=68SFU.

We can calculate the G/T from these observations, G/T=-7.15dB/K.

If we assume the gain figure is correct, then the system noise temperature is 10^(G/10-G/T/10)=10^(22.65/10–7.15/10)=955K. This is an awful system noise figure.

The next system uses a single EF7015 Yagi with claimed gain of 17dB, and measured Sun/ColdSky of 4.0dB at SF=70SFU. (I note he states SFI=70 which questions whether he is giving solar flux at 432MHz… but lets hope it was a typo).

We can calculate the G/T from these observations, G/T=-8.08dB/K. It is not much worse than the higher gain 4 x ZB7013 system.

If we assume the gain figure is correct, then the system noise temperature is 10^(G/10-G/T/10)=10^(17.0/10–8.08/10)=322K. Though only a third of the previous case, this is still a relatively high system noise figure… perhaps double of what you might like.

The next system uses 4x EF7015 Yagi with claimed gain of 22.94dB, and he suggests Sun/ColdSky of 10.0dB at SF=70SFU. (I note he again states SFI=70 which questions whether he is giving solar flux at 432MHz… but again lets hope it was a typo).

We can calculate the G/T from these observations, G/T=-0.33dB/K. It is quite a leap from the 1 x EF7015 system, and it is not clear that he actually measured 10dB. Lets assume that he did, that he is not just making wild estimates (and a 8dB improvement would be a wild estimate).

If we assume the gain figure is correct, then the system noise temperature is 10^(G/10-G/T/10)=10^(22.94/10–0.33/10)=212K. System noise is getting closer to what might be expected.

With G/T around 0dB/K, it is at the low end of 432MHz EME station performance, G/T up to nearly 10dB/K would represent the state of the art from 4 x Yagi systems (though more than twice as long).

It is a common practice amongst Hams to use 10.7cm flux measurements as the flux density at other frequencies. This is wrong, and if the solar flux values quoted in (Haefner 2010) are 10.7cm flux measurements, then the above results will be wrong as flux at 432 might be more like 30-50% of the 10.7cm flux. He has not given dates for the experiments, so I cannot check archives. The simple fact is that there is no reliable accurate way to preduct solar flux at 432MHz from flux at 10.7cm, examination of historical data will demonstrate that fact.

To use Sunspot number (SSI, T index) give even worse error.

Sometimes the uncertainty can be resolved by calculating a solution assuming that solar flux was given for 10.7cm and testing its feasibility. If that is done for the 1 x EF7015 scenario, and we assume solar flux was around 30SFU, then system noise temperature would be around 140K which is quite believable and it does not help to resolve the ambiguity about the author’s meaning of solar flux values.

It seems Hams boast Sun/ColdSky Y factors without understanding that they have no stand alone value, but they are very useful if the solar flux at the frequency of observations is known and the G/T figure calculated.

If I was a betting person, I would bet on Ham behaviour and that all the solar flux observations were in fact for a frequency of 2800MHz, more than 6 times that of the Sun Y factor measurements.

- Duffy, O. 2007. Measuring system G/T ratio using Sun noise. VK1OD.net (offline).
- ———. 2009. Quiet sun radio flux interpolations. https://owenduffy.net/calc/qsrf/index.htm.
- ———. 2014. Measuring G/T. https://owenduffy.net/blog/?p=1490.
- Haefner, A. 2010. 432MHz EME with a small Antenna In DUBUS. http://www.g0ksc.co.uk/DJ3JJ.pdf(accessed 27/04/14).
- ITU-R. 2000. Recommendation ITU-R S.733-2 (2000)Determination of the G/T ratio for earth stations operating in the fixed-satellite service .

]]>

The report shows solar flux at 1296MHz to be 106SFU.

On returning the model form, the solar flux at 10.7cm has been updated to 127SFU in agreement with the previous form. There is no sign that the solar flux at 1296MHz has been copied back.

We have a model at 1296MHz:

- after fetching current solar flux, it is shown as 127 @10.7cm;
- antenna gain is 30dB (upper right);
- Tsys=60K (red text at middle left);
- Y factor for Sun/ColdSky is shown as 12.31dB.

From these we can calculate G/T=gain-10*log(Tsys)=12.22dB/K.

We can also calculate the solar flux at 1296MHz that would give rise to a 12.22dB Y measurement for a system with G/T=12.22dB/K (ITU-R S.733-2). The solar flux is 63SFU, about half of what is shown on the form (admittedly labelled as @10.7cm), and nothing like its reported 106SFU at 1296MHz.

There appears to be some secret magic used here, (McArthur 2008):

This has been a significant break through as previous calculations were only accurate over a small sfu range. The ability to perform this calculation with complete accuracy has been the â€œHoly Grailâ€ I had been trying (unsuccessfully) to achieve for many years.

The correct Y factor for 106SFU @ 1296MHz for a system with G/T=12.22dB/K is 14.49dB, not 12.13dB as shown in EMRCalc v9.09.

I have not raised this with the EMRCalc’s author, why would I bother when he says (McArthur 2008):

The figures calculated by the software are â€œfirst principlesâ€ and if you do not get close to those predicted then it is your system and not the software that is the problem!

- ———. 2007. Measuring system G/T ratio using Sun noise. VK1OD.net (offline).
- ———. 2009. Quiet sun radio flux interpolations. https://owenduffy.net/calc/qsrf/index.htm.
- ———. 2014. Measuring G/T. https://owenduffy.net/blog/?p=1490.
- McArthur, D. 2008. Sun noise and measurements. http://www.vk3um.com/Documents/SunNoise_Measurements.pdf (accessed 26/04/2014).
- ITU-R. 2000. Recommendation ITU-R S.733-2 (2000)Determination of the G/T ratio for earth stations operating in the fixed-satellite service .

]]>

Above is his graphical summary of the measurements.

This article concentrates on just the measurements of K5SO station and the underlying model.

It is worth noting that the observations of all three stations fall within less than 25% of the range shown in the graph… so most of the graph is an extrapolation, and over a relatively huge range.

K5SO gives a model for the behaviour, it is SN(dB)=10*log(((0.72*(SF-64))+47)/A).

Note that he defines SF to be the solar flux at 2800MHz whereas the noise measurements are made at 1296MHz.

This expression can be simplified to SN(dB)=10*log((0.72*SF+0.92)/A), or as a simple ratio SN=(0.72*SF+0.92)/A =0.72/A*SF+0.92/A …(eqn 1).

For the purpose of this article, I will refer to K5SO’s quantity SN as Y, the ratio of noise pointing at the Sun to that pointing to cold sky.

The power received when pointing to the Sun includes both noise due to the Sun itself (lets call it Th) plus other noise, external and internal to the system (lets call that Tc).

So, Y=(Th+Tc)/Tc=Th/Tc+1 …(eqn 2).

Note the similarity of form of eqn 2 and eqn 1. That suggests that K5SO’s factor of 0.92/A is in fact 1.

Above is a comparison of the two formulas. A value of A=0.55 has been used to calibrate to K5SO’s curve.

So, we now know that Th/Tc=0.72/0.55*SF=1.31*SF. The factor 1.31 accounts for station antenna gain, and Sun noise at 2800MHz relative to system noise (Tc) at 1296MHz. That implies that there is a constant relationship between received Sun noise at 1296MHz and that at 2800MHz… but is there?

Above is a comparison of SF at 1415MHz compared to 2695MHz as measured at Learmonth observatory over 45 days in 2014. Clearly there is not a fixed relationship between them (they are produced in different parts of the chromosphere), and there is likely to be a somewhat similar variation in solar flux at 1296MHz compared to 2800MHz if one was to measure it.

Using 2800MHz SF just adds statistical noise. Whilst solar flux is not measured in these observatories at 1296MHz, they measure a number of frequencies and It should be better to use a sensible interpolation of those measurements.

Above is a plot of the distribution of the ratio of solar flux at 1415MHz to 2695MHz for the same data. This variation would give rise to 3σ uncertainty of 0.5dB, a little worse probably for the actual case of 1296MHz from 2800MHz, a significant amount which can easily be reduced subtantially by a better projection of 1296MHz from the observatory data.

K5SO gets a Y factor of about 21dB when 2800MHz solar flux is 100SFU. If we were to assume that at that time, solar flux at 1296MHz was say 85SFU, we can calculate G/T=20dB/K (making a beamwidth correction for a 39dB dish). Assuming antenna gain to be 39dB, that indicates a system noise temperature (Tc) of 78K which is credible.

W2UHI’s gain is probably around 36dB, G/T is probably around 16dB/K and assuming antenna gain of 36dB implies a system noise temperature (Tc) of 102K which is again believable.

Similarly, VK7MO’s gain is probably around 28dB, G/T is probably around 9dB/K and assuming antenna gain of 28dB implies a system noise temperature (Tc) of 113K which is again believable.

So, above is a graph similar to K5SO for a family of G/T values. The graph is constructed using the technicque and formulas given in Measuring G/T. (One must apply a beamwidth correction factors for narrow beamwidth antennas, and in the above chart that is calculated based on G/T and an assumed system noise temperature Ts.)

The similarity of the above chart and K5SO’s shows that conventional theory underlies K5SO’s experimental observations.

- Duffy, O. 2006. Effective use of a Low Noise Amplifier on VHF/UHF. VK1OD.net (offline).
- ———. 2007. Measuring system G/T ratio using Sun noise. VK1OD.net (offline).
- ———. 2009. Quiet sun radio flux interpolations. https://owenduffy.net/calc/qsrf/index.htm.
- ———. 2014. Measuring G/T. https://owenduffy.net/blog/?p=1490.
- ITU-R. 2000. Recommendation ITU-R S.733-2 (2000) Determination of the G/T ratio for earth stations operating in the fixed-satellite service .
- K5SO 2007. Sun noise. http://www.k5so.com/Using_sun_noise.html (accessed 24/04/14).

A metric that may be used to express the performance of an entire receive system is the ratio of antenna gain to total equivalent noise temperature, usually expressed in deciBels as dB/K. G/T is widely used in design and specification of satellite communications systems.

**G/T=AntennaGain/TotalNoiseTemperature 1/K**

**Example:** if AntennaGain=50 and TotalNoiseTemperature=120K, then G/T=50/120=0.416 1/K or -3.8 dB/K**.**

The utility of G/T is that receive S/N changes dB for dB with G/T, in fact you can calculate S/N knowing G/T, wavelength, bandwidth and the field strength of the signal (Duffy 2007).

**S/N=S*λ ^{2}/(4*π)*G/T/(k_{b}*B)** where:

S is power flux density;

λ is wavelength;

k

B is receiver quivalent noise bandwidth

Usage in this article is consistent with the industry standard meaning of G/T given at (ITU-R. 2000) (as opposed to the meaning used by some Hams who have appropriated the term for their own purpose).

The expression above can be rearranged to make G/T the subject, and if S/B is replaced with the known power flux density of a celestial object, measurement of S/N or Y factor (hot to cold sky) allows calculation of G/T.

Note that published figures for power flux density for celestial object are usually for a randomly polarised signal and they must be halved to determine the power received by a linear or circularly polarised antenna.

ITU-R S.733-2 gives the rearranged / revised formula.

ITU-R S.733-2 gives details of corrections that may be needed.

The Sun is a relatively high intensity source, and though its flux density is not sufficiently predictable, it is measured at several frequencies at several solar observatories each day at local noon and those measurements can be applied to Y measurements taken at about the same time.

The calculator embedded in Quiet sun radio flux interpolations can be used to fetch recent data for any of the observatories and calculate cubic spline interpolations to intermediate frequencies (eg Ham bands). It is a mistake to use the published 10.7cm flux figures, or simple extrapolations for this purpose, see (Duffy 2009).

A technique then is to find the recent solar flux values for the frequency of interest and predict the Y factor from modelled G/T to ensure that the measurements are practical (ie that expected Y can be measured with good accuracy, it isn’t too small or too large).

Measure Y several times around noon and find the highest value of that can be consistently obtained as the best capture of quiet solar conditions, and record all details of the experiment.

Later, use Quiet sun radio flux interpolations to get the measured noon solar flux at the observatory nearest to you and record it.

Then calculate the G/T figure from all the information you have… and share it. Very few hams share reliable G/T measurement data for their station, it is helpful to others in understanding the state of the art.

LA8LF gives information on the performance of his 1296MHz station with 5m parabolic dish including that the receiving system measures 16.5 dB of Sun Noise with a flux of 71.

Let us estimate the antenna gain assuming aperture efficiency of 60% to be 34.4dB. This will give rise to a very small beamwidth correction factor, 0.04dB.

Assuming that the solar flux reported is correct (ie for 1296MHz and concurrent with his measurements), and that the 16.5dB Y factor is accurate, we can calculate G/T to be 16.0dB/K.

Cross checking, it G/T is 16.0dB/K and Gain is 34.4dB, then T=69.2K which is a believable figure for such a system pointing high in the quiet sky.

A G/T model was constructed of the station.

From G/T and current solar flux, Y was expected to be around 10dB, accurate measurement was quite feasible.

(Duffy 2007) reported measurements of the 2.4m dish on 2304MHz. Gain was 32.6dB and Sun noise rise (Y) was measured at 6.28dB when flux at 2304MHz was 61SFU.

Calculated G/T was 10.4dB/K which reconciled well with the modelled G/T confirming station performance.

- Duffy, O. 2006. Effective use of a Low Noise Amplifier on VHF/UHF. VK1OD.net (offline).
- ———. 2007. Measuring system G/T ratio using Sun noise. VK1OD.net (offline).
- ———. 2009. Quiet sun radio flux interpolations. https://owenduffy.net/calc/qsrf/index.htm.
- ITU-R. 2000. Recommendation ITU-R S.733-2 (2000) Determination of the G/T ratio for earth stations operating in the fixed-satellite service .
- Thorrud, Anders (LA8LF). LA8LF EME Station Info. http://www.nitehawk.com/rasmit/la8lf.html (accessed 24/04/14).

Above is a sample report.

The report above was fetched from NOAA, and contains Learmonth data for the same day, 22nd April.

The immediate issue is that the figures are different, though represented as being from the same source.

The difference is not huge (0-17%), but one is left wondering why two source of the same information differ at all.

The first figure shows some calculated values for the Ham bands.

Above is a chart of EMRCalc’s Learmonth figures, and the interpolated / extrapolated values for the Ham bands. The Ham band values appear to be linearly interpolated between the observation data (a linear spline), the values at 900 and 1296MHz show that fairly clearly.

The problem with linear spline is that there is a discontinuity at each observation point, the slope of the interpolated characteristic changes instantly other side of the observation and this is most unlikely to happen in nature… ie the technique may not be a good estimator. The EMRCalc estimate for 900MHz is about 10% lower than a cubic spline estimate.

The estimates below 245MHz are clearly extrapolation in that they are outside the observation data. They appear to be on some kind of curve.

There appears to be no hints in the documentation in the methods used for interpolation and extrapolation.

The CS Interp line is constructed from a cubic spline interpolation of the observation points. It has the advantage that there is not a slope discontinuity at the observation points, it is a smoother curve and for that reason, a more believable explanation of the observations.

Using a cubic spline interpolation on the NOAA Learmonth dataset gives a solar flux figure 9% higher for 1296MHz than in the EMRCalc report for the same time.

Though there are several solar observatories around the world reported in the NOAA report, it seems EMRCalc uses only one. It may be that if your interest is in solar flux in another part of the world, noon observations from an observatory nearer your own location might be more accurate, eg for European locations, San Vito data might be better.

Above is a table of cubic spline interpolations for noon NOAA data from San Vito, same day but seven hours later. A data pair of [0,0] is prepended to the observation data to allow cubic spline interpolation below 245MHz, though in the strict sense, it is extrapolation (denoted by the asterisk).

Above is the NOAA data for Learmonth using the same interpolation methods. The values obtained are different to EMRCalc.

A follow up on the reconciliation issue of solar flux stats.

Above is the solar flux statistics in EMRCalc immediately after an update. The datestamp turns out to be that of the time of update, not that of the underlying data.

Above is the IFLUX report from Learmonth’s website at about the same time. The data reconciles with EMRCalc, but it is from the two days prior.

Above is from NOAA at about the same time, and it shows data from one day ago, in 6 hour it will have todays data.

So, the issue seems that EMRCalc’s source is stale, and datestamp misrepresented in its report.

If you make measurements at a certain time, you need to collect flux measurements that clearly and correctly identify the time so that you can use accurate figures.

- Duffy, O. 2009. Quiet sun radio flux interpolations. https://owenduffy.net/calc/qsrf/index.htm.