Demographic trends
For the period 2001-2016, a total of 18,115 RRV notifications were analysed across the SEQ region. During this period, the mean annual RRV notification rate across all of SEQ was 41 cases/100,000 population. Mean annual rates were higher in females versus males, at 44/100,000 population compared to 34/100,000, respectively, though this difference was not significant (Mann-Whitney U=84, p=0.102). For both genders, the highest incidence rates occurred in the 40-44 and 45-49 age categories (69 and 66/100,000, respectively) and the lowest in the two age categories <10 years (1 and 3/100,000, respectively) (Fig. 2). The overall trend showed a gradual increase in incidence from birth up to age 29 years (mean incidence of 17/100,000 across these age groups), peaking between ages 30-64 years (mean 59/100,000), and dropping again ≥65 years (mean 27/100,000). This trend was consistent across all years of this study. Statistical comparisons within and between these three broad age groups indicated that rates did not differ significantly between adults from 30-64 years (Kruskal-Wallis H=8, p=0.239), but that this group’s rates were significantly higher than those aged ≤ 29 years (Mann-Whitney U=1175, p<0.0001), and ≥ 65 years (U=1589, p<0.0001).
Fig. 2. Mean annual RRV incidence in South East Queensland by age-group, 2001-2016.
Age-adjusted mean annual incidence rates are grouped into 18 age groups of SEQ, between 0 and 95 years. Mean rates per group are indicated by X, and the median by the horizontal line across each box. Whiskers indicate the estimated minimum (lower whisker) and maximum (upper whisker) incidence values per age group (equal to the 1st and 3rd quartile -/+ 1.5 x the inter-quartile range, respectively). Outlying points indicate extreme values reported during large outbreak years, e.g. during the outbreaks of 2014 and 2015.
Temporal trends
RRV disease occured in SEQ throughout all months of the year; however, a distinct seasonal pattern in the timing of annual notification peaks was shared amongst all LGAs. The peak annual notification period was typically between February and May, and this was consistent across years with higher and lower notifications (Fig. 3). These months also showed the highest variability in case numbers (with a monthly average of 186 cases across SEQ, SD 246). Conversely, notifications during the winter and spring months between June-November were generally low and relatively stable across the region (monthly average of 46 cases, SD 35). The relative magnitude of RRV outbreaks across the region were variable between years, with larger and smaller outbreaks occurring in intermittent years (the monthly temporal trend is shown in Additional file 1). Long-term trends indicated that outbreaks generally occurred synchronously across the region, rather than initially occurring in one LGA and then spreading to another. This was also the case during largest recorded outbreak of 2015, which began earlier than usual (in late 2014) and peaked in February-March 2015, with timing consistent across all SEQ LGAs.
Fig. 3. Mean monthly trend of RRV case notifications in South East Queensland, 2001-2016.
Mean monthly case numbers across the 16-year period are shown, with extreme case report values from large outbreak years indicated as outliers. Mean rates per group are indicated by X, and median by the horizontal line crossing each box. Whiskers indicate the estimated maximum and minimum case numbers per month across the 16 years (equal to the 1st and 3rd quartile -/+ 1.5 x the inter-quartile range, respectively). *Total cases reported in February 2015 (1,357) were the highest ever reported in a single month.
Spatial trends
Spatiotemporal patterns of RRV disease varied between the 9 LGAs, with the highest case numbers in the most populated LGA, Brisbane City (5,352 total cases), and the lowest in the sparsely populated Scenic Rim Region (393 total cases) (Table 1). Annual case numbers for each LGA can be found in Additional file 2. Mean annual incidence rates across the 16-year period varied from 30/100,000 in Gold Coast City up to 130/100,000 in Noosa Shire (Table 1), with rates in Noosa Shire being significantly higher than all other LGAs (Kruskal Wallis H=50, p<0.0001; Mann-Whitney U=58, p=0.007 for pairwise comparison between Noosa Shire and Scenic Rim LGAs). Rates were also routinely high in the Sunshine Coast Region and sporadically high in the Scenic Rim Region. These 3 areas (Noosa Shire, Scenic Rim Region and Sunshine Coast Region) had the highest incidence rates overall, and had more rural characteristics (lower population density and proportion of urban SSCs) compared with the 4 largest cities: Brisbane, Gold Coast, Ipswich and Logan (Table 1). These 4 higher-density cities had the lowest rates of all LGAs.
Table 1. Summary characteristics of each Local Government Area (LGA) of South East Queensland during the study period, 2001-2016.
LGA name*
|
Average population (2001-2016)
|
Population density
(per km2)
|
Total cases
(2001-2016)
|
Proportion of urban SSCs
|
Mean annual incidence rate
(per 100,000)
|
Brisbane City
|
1,023,663
|
762
|
5,352
|
94.2%
|
33
|
Gold Coast City
|
475,169
|
356
|
2,252
|
74.1%
|
30
|
Ipswich City
|
157,963
|
146
|
1,044
|
58.0%
|
41
|
Logan City
|
271,082
|
283
|
1,456
|
81.7%
|
34
|
Moreton Bay Region
|
357,532
|
175
|
3,097
|
58.5%
|
54
|
Noosa Shire
|
41,008
|
47
|
850
|
29.2%
|
130
|
Redland City
|
134,328
|
250
|
991
|
63.6%
|
46
|
Scenic Rim Region
|
35,443
|
8
|
393
|
2.8%
|
69
|
Sunshine Coast Region
|
256,982
|
114
|
2,680
|
50.4%
|
65
|
Total
|
2,753,170
|
188
|
18,115
|
61.0%
|
41
|
* See Fig. 1 for location of each LGA within the study area.
|
Fig. 4. Mean annual RRV incidence of State Suburb Codes (SSCs) in South East Queensland, 2001-2016.
a) Mean annual incidence rate, and b) smoothed mean annual incidence rate for SSCs of South East Queensland over the 16-year study period. Seventeen unpopulated SSCs are indicated by the white (zero case) locations in b).
Fig. 5. Mean annual RRV incidence of rural and urban State Suburb Codes (SSCs), 2001-2016.
a) raw incidence rates, and b) smoothed incidence rates for the 757 SSCs of South East Queensland, shown according to the 9 Local Government Areas (LGAs) they are located within, and all SSCs combined. Rural and urban SSCs are indicated by blue and red circles, respectively. Two outlying SSCs (in Ipswich City and Noosa Shire) had incidence values beyond the scale shown here, but are listed in Additional file 4. See Fig. 1 for location of each LGA within the study area.
Hot spot analyses
Using raw mean annual incidence rates, a total of 127/757 (17%) SSCs were identified as hot spots and 352/757 (47%) as cold spots in any individual year. These increased to 178 and 458 SSCs, respectively, using smoothed rates. A comparison of raw versus smooth analyses for annual hot and cold spots across individual years was shown in Additional file 5. In either raw or smoothed analyses, 86 hot spots and 272 cold spots were persistent (present in ≥2 years). There were 14 SSCs that were both hot and cold spots in ≥2 years that were excluded, leaving 72 hot and 258 cold spots. Of these, 45 hot spots and 154 cold spots were identified as persistent in both raw and smooth analyses (shown in Additional file 6). Persistent hot spots were similar to those identified using mean annual rates (n=56 mean hot spots), while mean cold spots differed (n= 47 mean cold spots). Although hot spots were geographically dispersed across all LGAs, the SSCs with the most persistent hot spots were in the Sunshine Coast Region and Noosa Shire LGAs. The same hot spots were rarely detected in consecutive years, although some were detected in multiple years, up to 7/16 years (Additional file 7). Conversely, cold spots tended to persist more in the same SSCs across several years, particularly in the Scenic Rim Region where there were very low populations and cases (Additional file 6).
A visualisation of hot spots relative to rural and urban areas of SEQ is shown in Fig. 6a. Hot spots tended to be most focused around the edges of where major urban and rural areas intersect. Of the 72 persistent hot spots detected in either raw or smoothed analyses, 35 were located in urban and 37 in rural SSCs (for the 45 hot spots shared between both analyses, 19 were urban and 26 rural). There also appeared to be diverse land use types within or adjacent to hot spot SSCs (Fig. 6b). All hot spots contained some degree of urban infrastructure, and many were located in close proximity to either dryland agriculture and plantations or to major water bodies. While some more inland hot spots were surrounded largely by conservation and natural environments, these environments appeared to be most often identified as cold spots (Fig. 6, Additional file 6). Very few hot spots were located centrally within major urban areas, rather than on the edge. Only one persistent hot spot SSC was identified within Brisbane City LGA (Chelmer, identified in 2/16 years) in the raw rate analysis only. The largest cities including Brisbane, Gold Coast and Logan LGAs were more commonly dominated by cold spots rather than hot spots (Additional files 5 and 6).
Fig. 6. Persistent high incidence hot spots in South East Queensland, 2001-2016.
Persistent hot spots identified for State Suburb Codes (SSCs) are shown relative to a) urban and rural areas of South East Queensland; and b) different land use types of South East Queensland. Points indicate the centroid of each hot spot SSC. Hot spots detected using raw incidence rates only (yellow points; n=27 SSCs) and those detected in both raw and smoothed incidence analyses (black points; n=45 SSCs) are indicated. Note: in b) ‘Intensive uses’ refers to residential areas and urban infrastructure. See methods for further description.