3.1. Backscatter Coefficient Results over the Decade
The monthly mean backscatter coefficients were calculated over the decade for each of the four BSS. For brevity, only the LO results are shown for the monthly maps of the backscatter coefficient over the decade. However, the data reduction and detailed analysis are presented for all four sounders.
Figure 4 and
Figure 5 show maps of the monthly mean backscatter coefficients for the LO sounder from January 2011 through to December 2020. Higher backscatter coefficients are more commonly seen in the middle of the year (from around June to August). Conversely, little variation is noted from year-to-year, indicating that while there is considerable seasonal variation, there is little variation over the 11 year solar cycle. The lack of variation over the solar cycle is expected as the solar cycle is a driver of the ionospheric propagation conditions and has little effect on terrestrial weather which is the main driver of the sea state. It is the variation of the sea state which affects the backscatter coefficient over the ocean.
We can make an initial assessment of the difference between land and sea backscatter by examining the monthly standard deviation of the backscatter coefficients for each range azimuth cell calculated over the decade.
Figure 6 shows the decade mean of the monthly standard deviation in the backscatter coefficients for each of the sounders. The standard deviation for LO was generally greater over the ocean (especially for the Arafura Sea (around 8°S, 136°E) and the Pacific Ocean (around 5°N, 135°E) where the mean standard deviation is as large as 9 dB). This suggests the sea backscatter coefficients change throughout each month more than the land backscatter coefficients, which is expected from changes in the sea conditions due to weather (as discussed in our earlier paper [
15]). A range dependence in the standard deviation can be seen, which is attributed to limitations of the backscatter ionogram synthesis. Larger errors occur in the propagation model for longer distances. Consequently, the difference between the real and synthetic backscatter ionograms is potentially larger for longer distances (see our earlier paper [
2]) which will lead to a larger random error in the calculation of the backscatter coefficient. This results in a larger standard deviation in the results at longer ranges.
The decade mean of the monthly mean backscatter coefficients for each of the sounders is shown in
Figure 7. Various topographic features are visible in these maps of the backscatter coefficient. Of particular note is the strong backscatter from New Guinea (at around 5°S, 140°E) which is due to a large mountain range which runs centrally along the length of the island. This backscatter is stronger in the LAE sounder than LO which is due to the LAE sounder observing the mountain range orthogonally. In general, the lowest backscatter coefficients occurred over desert regions (such as at 20°S, 130°E and 23°S, 140°E), while the highest land backscatter coefficients occurred where the terrain was rough with mountainous features or highly vegetated (such as at 5°S, 140°E). We do not discuss this further here as land backscatter coefficients were investigated in detail in our earlier paper [
2]. An aspect dependence can also be seen in the backscatter coefficients, where multiple sounders observe the same location from different directions (for example at 16°S, 130°E and at 5°S, 140°E).
Figure 8 shows the mean number of data points contributing to each range-azimuth cell each month. The maximum number of data points available occurred at ranges of 1500–3000 km as expected, as this is a common one hop range due to the ionospheric propagation support. This is important as areas with more data points have better statistics for the calculation of the surface backscatter coefficient, while the results for areas where few data points are available must be treated with caution as discussed later.
3.2. Seasonal Variation in the Surface Backscatter Coefficients
To investigate the seasonal variation in the surface backscatter coefficients, the changes in the monthly mean backscatter coefficients for each location over time need to be examined. Before we can do this for the entire data set we first investigate the monthly backscatter coefficients at a single location: beam 4 of the LO sounder at a range of 3000 km (3°S, 127°E). This location was chosen as it exhibits a large standard deviation in the monthly mean backscatter coefficient over the decade of data. The entire time series of the monthly mean backscatter coefficient for this location is plotted in
Figure 9 (left) and as a superposed epoch against the month of year in
Figure 9 (right). The backscatter coefficient followed similar trends for all of the years. A clear trend of higher backscatter coefficients in the middle of the year can be seen. A Fourier transform was used to find the peak frequencies in this time series and the resultant amplitude spectrum is shown in
Figure 10. No window was applied and missing data was filled with the mean value of the nearest data either side, weighted by their temporal offset to the missing data points.
A clear peak in frequency at one cycle per year was found. A weaker biannual peak was also seen. This was not surprising as typically ocean wave heights follow a sinusoidal annual cycle with larger significant wave heights in winter due to seasonal changes in high latitude storm patterns [
25]. Colosi, Villas Bôas and Gille [
25] investigated the annual and biannual cycles in the significant wave height and wind speed globally. They showed there was a strong annual cycle in the wave height and wind speed in this location, and a weaker biannual cycle. It was suggested that the South Asian monsoon’s biannual occurrence may play a role in the biannual cycles in this region. The phase of the FFT indicates that the annual cycle peaks in July and the biannual cycle peaks in May and November.
The frequency and amplitude of the primary peak was calculated in the manner described above for each range azimuth cell of the LO data and plotted in
Figure 11. Locations were not considered for this analysis if the monthly mean data availability was less than 50%. In the left panel, data with low statistical significance (see below) are not plotted. The primary-peak amplitude over the Australian landmass was, as expected, much less than over the sea. For most locations (considering both the land and sea) the peak frequency was one cycle per year. A peak frequency of two cycles per year was seen over Northern Australia; however, the amplitude of this peak was low, which shows there is little seasonal variation in this area.
Figure 12 shows a histogram of the primary peak frequencies for the LO range-azimuth cells. This shows that a peak frequency of one cycle per year was the most common.
To determine the significance of the spectral peaks we employed a statistical technique described by [
26]. Gaussian noise was generated with the same mean, standard deviation and number of data points as the monthly mean backscatter coefficients data for each range-azimuth cell. The fast Fourier transform was applied to this noise, and this was repeated for 1000 sets of noise. An amplitude threshold was then set at a level which is exceeded by less than 1% of the noise data. This amplitude threshold was then used to determine which peaks from the backscatter coefficient data were significant. Peaks in the Fourier amplitude spectrum were considered significant if the peak was greater than the amplitude threshold. The largest peaks in the Fourier amplitude spectrum for majority of the range-azimuth cells were considered significant.
This analysis shows that while the amplitude of the biannual peaks from the Northern Australian landmass is low, these peaks are statistically significant. We currently do not have an explanation for this variation as the rainfall and consequently changes to ground conductivity in this region has a strong annual variation peaking in late January. It is possible that seasonal variations in the ground properties at the backscatter sounder transmit and receive sites may affect the antenna patterns and thus the power transmitted/received, especially for low elevation rays. This could contribute to weak seasonal variations in the power of the backscatter ionograms that has not been accounted for by the model. Consequently, this could affect the backscatter coefficients, and may introduce the observed weak biannual cycle over land if the wet seasons at the transmit and receive sites and the observed locations are out of phase. However, the rainfall at Longreach, Queensland, is low and peaks in early February with a minimum in August, i.e., in phase with the observed locations and so we reject this hypothesis as an explanation of the observed weak biannual variation.
The frequency and amplitude of the second largest peak of the backscatter amplitude spectrum was also calculated in the manner described above for each range-azimuth cell of the LO data and plotted in
Figure 13. Again, only the peaks considered statistically significant were included in the left most panel of
Figure 13. The secondary-peak amplitude for all locations was low, especially in comparison with the amplitude of the primary peaks over the ocean. This indicates that most of the temporal variation in the backscatter coefficients occurs annually. A weak biannual variation is seen in the sea backscatter coefficients in certain regions.
Figure 14 shows a histogram of the secondary-peak frequencies for the LO range-azimuth cells. Many of the secondary peaks in the Fourier amplitude spectrum were not considered significant. The frequencies of the statistically significant secondary peaks were annual and biannual. Some of the spectral peaks at lower frequencies were also statistically significant. However, their amplitude was low and further it should be noted that these were generally for locations where a high percentage of data points were missing, and so these results should be treated with caution.
Results from a Fourier analysis of data from each of the sounders are shown in
Figure 15. The Fourier analysis for AS was done on a shorter nine-year data set, as no data was available before 2012. As with the Longreach sounder, the most common peak frequency in the Fourier spectra of the other sounders was also 1 year
−1. This peak was most prominent for the ocean around Indonesia (at 4°S, 125°E) and New Guinea (at 1°N, 145°E), which indicates that the largest seasonal variations in the sea surface backscatter coefficient occurs in these locations. Hemer et al. [
27] show that for the Australian region the largest inter-annual variation in the significant wave height (normalized by the mean significant wave height) occurs around the location of these most prominent peaks: in the Gulf of Carpentaria (approximately 13°S, 139°E), the Arafura Sea (approximately 9°S, 136°E) and off the Queensland coast (approximately 15°S, 148°E). A yearly cycle in the surface backscatter coefficients over the land is also seen for many locations, although for this case the variation is weak as previously discussed. The biannual variation in the backscatter coefficient noted previously in Northern Australia (around 15°S, 135°E) in the LO data is also seen in the LAE results. Again, while statistically significant, this variation is weak as indicated by the peak frequency amplitude plots. The rainfall at the LAE BSS sites is low with little seasonal variation which peaks in November and so, as with the LO results, we cannot explain this biannual variation as being due to an instrumental effect caused by ground conductivity changes at the BSS transmit and receive sites.
Over the Australian land mass the amplitude of the peak frequency was low, indicating the seasonal change in the surface backscatter coefficients was small. Thus, for most applications a constant mean value for land locations will be sufficient to describe the backscatter coefficient. However, to describe the mean sea backscatter coefficient over time, monthly values must be used, noting this would not capture the day to day variations in the sea surface backscatter coefficients.
The phases associated with the annual and biannual variations in the backscatter coefficients are shown in
Figure 16. In general the annual cycle (
Figure 16 left) peaks around June to August for most ocean locations, which is expected [
25]. Over the land the annual cycle tends to peak earlier in the year around March to May. This is around the end of the wet season in Northern Australia, when it is expected the ground would be at its wettest and hence have maximum conductivity.
The phase associated with the biannual cycle (
Figure 16 right) is generally around 4 to 5 months, such that the backscatter coefficients peak around May/June and November/December each year. However, near the coast of Australia in the Indian Ocean (15°S, 117°E) around to the Timor Sea (11°S, 139°E) and also in the Gulf of Carpentaria (14°S, 139°E) the peak in the backscatter coefficient associated with the biannual cycle occurs around January and July each year.
3.3. Monthly Maps of the Surface Backscatter Coefficients
The analysis in the previous sections indicates that the variability in the sea backscatter coefficient is seasonal with little year-to-year variation and little variation at all for the land backscatter. Thus, we can construct monthly maps of the backscatter coefficient by combining the entire decade of data. We do this by calculating the decade mean of the monthly backscatter coefficients at each location.
The mean monthly backscatter coefficient maps for LO are shown in
Figure 17. Changes in the backscatter coefficient throughout the year are clearly seen. The sea backscatter coefficients around Indonesia (at 4°S, 125°E) and New Guinea (at 5°S, 140°E) are significantly larger in the middle of the year, from around June to August, than at other times of the year. The backscatter coefficients over the Northern Australia landmass remained relatively constant throughout the year, with the lowest backscatter coefficients appearing over the Tanami desert (located around 20°S, 130°E).
Figure 18,
Figure 19 and
Figure 20 show the mean monthly backscatter coefficients for LAE, LAW and AS. Similar trends to the LO results were found for these sounders, with larger sea backscatter coefficient generally seen in the middle of the year. In many areas the backscatter coefficients for the LAE sounder are significantly lower than the other sounders. This can be attributed to the field of view of LAE containing much of central Australia, which is predominantly desert and so has a lower backscatter coefficient. This is discussed in detail in our earlier paper [
2].
The LAE observations also display an interesting and somewhat surprising feature with a strong seasonal variation in the backscatter coefficient from the island of New Guinea. The backscatter coefficient in central New Guinea is ~7dB stronger in August than in January. The single-sided amplitude spectrum from the FFT on the backscatter coefficient data from this location shows strong annual and biannual peaks (
Figure 21). This behavior is in contrast to other land areas which show no seasonal variation in the backscatter coefficient. The strong backscatter from the inland region of this island is, as discussed earlier, due to the large mountain range that runs along the length of the island. However, we have no explanation for the seasonal variation of the backscatter from this region. Indeed, we would expect any variation in the backscatter coefficient to be due to changes in the ground conductivity with a maximum during the rainy season. However, the rainy season in New Guinea occurs from December through to April when the backscatter coefficient is typically lower.
The standard deviation between the months was also calculated. We have not included the plots here. However, the standard deviation of the land backscatter coefficients is less than 2 dB for most locations, while the standard deviation in the sea backscatter coefficients is larger, rising to 5–6 dB in some locations. This is in line with the Fourier analysis results discussed in the previous section and shows the sea backscatter varies considerably more throughout the year than the land backscatter. This is not surprising given the sea state is governed by the meteorological conditions. Hence, while the use of a single temporally invariant backscatter coefficient for land locations may be sufficient for many applications and locations, a seasonally varying backscatter coefficient is required for sea backscatter.