Dengue nowcasting in Brazil by combining official surveillance data and Google Trends information

Introduction

Dengue is a viral disease transmitted by mosquitoes of the Aedes genus such as Aedes aegypti and Aedes albopictus, which poses a significant global health threat, particularly in tropical and subtropical regions like Latin America, Southeast Asia, the Pacific Islands, parts of the Middle East, and Africa (1). The disease ranges in severity from mild symptoms such as high fever, severe headache, retro-orbital pain, joint and muscle pain, and rashes, to severe forms like dengue hemorrhagic fever. This severe condition can lead to critical complications such as bleeding, plasma leakage, and organ failure, significantly increasing the risk of mortality (2).

Dengue transmission occurs through mosquito bites, with peak activity early in the morning and before dusk. The disease is highly contagious and can rapidly spread among populations, particularly in densely populated urban areas where mosquitoes breed in stagnant water. Several factors contribute to the prevalence of dengue, including climate, urbanization, and socio-economic disparities (3). Additionally, global warming and increased international travel have expanded the geographic range of dengue (4).

Annually, there are approximately 400 million infections, with 100 million cases exhibiting clinical symptoms ranging from mild to severe (2; 5). In 2024, Brazil is facing a severe dengue outbreak, with 9.48 million suspected cases and 5.32 million confirmed cases as of August (6). This surge has made Brazil the most affected country in the Americas. The increased incidence is attributed to factors like early transmission seasons, climate change, and the presence of all four dengue serotypes.

Challenges in dengue prevention and control include the absence of specific antiviral treatments and the complexity of developing vaccines for the virus’s four serotypes (7). Moreover, dengue symptoms can resemble those of other diseases such as chikungunya and Zika leading to misdiagnosis and complicating effective control measures. Dengue prevention efforts focus on reducing mosquito populations and avoiding bites through community engagement, personal protection, and environmental measures.

The rapid transmission and extensive impact of dengue underscore the importance of research and the ability to predict its outbreaks. Surveillance systems play a crucial role for guiding strategies for prevention and control, but traditional surveillance systems often rely on delayed or incomplete data due to underreporting, healthcare infrastructure limitations, and time lags in laboratory testing. For example, in Brazil, a suspected dengue case is required by law to be reported to the government authorities. The notification is made by authorized health professionals in Notifiable Diseases Information System (SINAN) (8). The notification data is then gathered in Brazilian Ministry of Health and the most up-to-date information is used by the Epidemiological Surveillance teams to understand current status of dengue cases through the country. However, according to (9), less than 50% of dengue cases are reported within the first week, no more than 75% are reported within four weeks, and fewer than 90% are reported within nine weeks. This latency necessitates nowcasting methods to estimate occurred-but-not-yet-reported disease cases for real-time decision-making.

Several nowcasting techniques have been developed in different settings. For example, (10) employed reverse-time discrete hazard functions and maximum likelihood estimation to handle reporting delays and nowcast AIDS cases in Canada. A Bayesian hierarchical model was proposed by (11) to improve prediction and management of a Shiga toxin producing Escherichia coli that caused a major outbreak in Germany in 2011. The model combined a survival regression model for the delay distribution, and a quadratic spline for the epidemic curve, utilizing the generalized Dirichlet distribution for flexibility in handling uncertainty. (12) proposed a Bayesian hierarchical model that jointly estimates the expected number of deaths, and the reporting delay distribution to nowcast COVID-19 fatalities in Sweden. This model offers enhanced predictive performance and flexibility by incorporating leading indicators such as the number of reported cases and COVID-19 associated ICU admissions.

A nowcasting by Bayesian smoothing approach capable of producing nowcasts in multiple disease settings was developed by (13). This approach learns the reporting delay distribution and the time evolution of the epidemic to produce nowcasts in both stable and time-varying case reporting settings. The approach was tested on dengue in Puerto Rico and influenza-like illness (ILI) in the United States. (14) presented a framework for addressing reporting delays in malaria surveillance in Guyana. The method combines a data imputation model and network models to refine case estimates using historical data, neighboring region data, and precipitation levels. (15) introduced a Bayesian framework with sliding windows for dengue surveillance in Bangkok, Thailand, addressing reporting delays by accounting for spatial and temporal variations. A Bayesian hierarchical model for dengue nowcasting in Brazil was developed by (9). This approach uses a Negative Binomial distribution for the reported cases with mean explained by spatial, temporal, and delay information, offering a robust correction to the reported cases. This model is used by the InfoDengue system to nowcast dengue in Brazilian municipalities (16).

The implementation of these methods are complex and require extensive data on the historical weekly reported cases for the disease under consideration. In recent years, methods utilizing digital data, such as Google Trends indices and Twitter (now X) information on disease-related keywords, have shown significant effectiveness in understanding and predicting disease activity levels. These methods leverage real-time search query data to enhance the accuracy of traditional models, allowing for more timely and reliable public health responses.

For example, (17) utilized used search query logs and modeling techniques such as Elastic Net regularized regression and Gaussian Process regression to nowcast influenza-like illness in the USA. (18) used an ARIMA model augmented with Google Flu Trends data for nowcasting influenza outbreaks in the USA. They showed the incorporation of real-time search query data improves prediction accuracy compared to a model that uses only case data.

A Hidden Markov Model combining cases and Google Trends information for disease prediction was proposed by (19). The model incorporated an autoregressive component describing case counts and a linear covariate representing Google Trends. The model was applied to predict dengue in Brazil, Mexico, Thailand, Singapore, and Taiwan, as well as influenza-like illness in the USA (20).

In (21), authors utilized Baidu search query data, which is similar to Google Trends, to nowcast hand, foot, and mouth disease across China. They utilized a meta-learning framework to dynamically select among predictive models including Principal Component Analysis, LASSO, Ridge Regression, and ARIMA. They showed the inclusion of Baidu Index data enhances prediction accuracy by providing real-time public interest metrics correlated with hand, foot, and mouth cases.

(22) used Twitter data to monitor dengue in Brazil. First, they analyzed tweet sentiments to filter tweets indicative of actual cases, and found a high correlation of the number of dengue-related tweets with official dengue data. Then, they constructed a regression model for predicting the number of dengue cases using the proportion of dengue-related tweets, and used it to develop a monitoring system that generated weekly heatmaps of dengue across cities in Brazil.

We are not aware of any approach that utilizes Google Trends data in Brazil. In fact, the only existing dengue tracking system in Brazil is InfoDengue, which provides weekly case reports of dengue and other arboviruses across Brazilian municipalities, along with nowcasts generated by a Bayesian hierarchical model (16). Therefore, it is crucial to investigate alternative methods.

To this end, in this paper we assess the value of Google Trends to nowcast weekly dengue cases in Brazilian states by fitting models that integrate reported dengue cases in Brazilian states provided by InfoDengue https://info.dengue.mat.br/, and Google Trends indices of dengue-related keywords at the state level in Brazil which can be obtained from https://trends.google.com/trends/. Then, we compare nowcasts produced by models that utilize only reported dengue cases, only Google Trends data, and a combination of both, with the model provided by InfoDengue (16). As a baseline, we also use a naive approach where nowcasts are considered as the reported cases the previous week. Our aim is to evaluate the value of real-time Google Trends data in producing accurate nowcasts using simple models that do not rely on incomplete recent case data, and to determine whether these nowcasts can compete with more complex and time-consuming models.

Additionally, recognizing the importance of timely and accurate data for dengue surveillance, we developed the Dengue Tracker website (https://diseasesurveillance.github.io/dengue-tracker/). This site is updated with weekly nowcasts for each Brazilian state, and presents information through interactive maps and time trend plots to inform decision-makers and the public about dengue activity levels in real-time.

The rest of the paper is organized as follows. First, we detail the study region, data sources, nowcasting methods, and performance assessment measures. Next, results are presented through explanatory tables and figures, highlighting error and uncertainty measures across states. Then, we describe the Dengue Tracker website developed to support real-time decision-making. Finally, the paper concludes with a discussion of the findings and future work.

Materials and methods

Study region

Brazil, the largest country in South America and the fifth-largest in the world, spans a vast geographical area characterized by diverse climates and ecosystems (23). Its territory stretches from the equatorial Amazon basin in the north to the temperate regions in the south. This wide latitudinal range encompasses various climatic zones, including tropical, subtropical, and temperate, which significantly influence the epidemiology of vector-borne diseases like dengue fever. Brazil is divided into 26 states and one federal district, each with unique geographical and socio-economic characteristics that impact dengue transmission and control (Figure 1).

• Download figure
• Open in new tab

Fig 1.

Map of South America with Brazil highlighted (left) and map of 26 states and the federal district of Brazil (right).

Several factors significantly contribute to the proliferation of dengue outbreaks in Brazil (24). Climate change, characterized by increased temperatures, has diminished geographical barriers to dengue transmission, particularly in southern Brazil, by reducing the seasonal cold periods that typically inhibit mosquito propagation. In the Amazon region, climatic changes have made previously protected areas more susceptible to dengue outbreaks.

Rapid urbanization has led to high population densities in large cities, fostering environments conducive to mosquito breeding due to inadequate infrastructure, such as insufficient piped water and waste management systems (25). This results in standing water containers and uncollected garbage, serving as breeding grounds for mosquitoes. Furthermore, high connectivity within Brazil’s urban network exacerbates the risk of disease spread. Major regional hub cities with frequent transportation and logistical connections facilitate faster virus transmission, while enhanced connectivity between urban and rural areas contributes to the spread of the virus from cities to the countryside (24).

Figure 2 illustrates the seasonal pattern of monthly dengue incidence rates in Brazilian states from January 2010 to July 2024, aggregated at the state level, with data sourced from InfoDengue (26). The figure reveals a seasonal pattern, with dengue outbreaks typically occurring from January to May. In certain regions, such as Acre, Rondônia, Mato Grosso, and Goiás in southwestern Brazil, outbreaks may even commence in November or December. The spread of the disease is generally minimal during the winter months. Moreover, it is observed that dengue outbreaks tend to be more severe every three to four years. This periodicity could be influenced by El Niño, which significantly impacts weather patterns in certain areas of Brazil (27). Notably, the states of Acre, Espírito Santo, Goiás, Mato Grosso do Sul, Paraná, Rio Grande do Norte, São Paulo, and Tocantins experience a more substantial impact from dengue, as indicated by consistently higher incidence rates.

• Download figure
• Open in new tab

Fig 2.

Dengue incidence rate (cases per 100k people) on a log₁₀ scale in Brazilian states from January 2010 to July 2024.

Google Trends data

Google Trends (https://trends.google.com/trends/) is a tool that provides anonymized and aggregated insights into global search behaviors, ensuring user privacy by anonymizing individual searches and consolidating them into high-level trends. This data is valuable for identifying emerging topics and seasonal trends. The Google Trends index for a specific keyword at a given time ranges from 0 to 100, calculated by dividing the number of searches for that keyword by the total number of searches in a specific region and timeframe, enabling fair comparisons between search terms, locations, and periods. This standardization eliminates disparities due to population size or total search volume. However, it may also create misleading representations, as a term could appear less popular in larger regions due to a dilution effect. Moreover, aggregation might obscure specific nuances or emerging trends within smaller subgroups, potentially limiting the depth of analysis (29).

Here, we utilize Google Trends data to understand dengue search behavior patterns that could complement official dengue data. To select the keywords for the Google Trend indices to include in the models, we calculated the Pearson correlation between country-level aggregated dengue cases and the Google Trends index for various dengue-related keywords identified in (30). We used data from January 1, 2013, to December 31, 2023, with correlations computed using monthly resolution data, as Google Trends data cannot be obtained at a finer temporal resolution for periods exceeding five years and two months.

Figure 3 illustrates the correlation between the number of dengue cases and the Google Trends index corresponding to each of the keywords considered. We observe the highest correlation of dengue cases and Google Trends index corresponding to keywords “sintomas dengue” “dengue” and “sintomas de dengue” with correlations of 0.93, 0.90, and 0.89, respectively. The intercorrelations among these keywords were also significant, with the highest correlation being 0.97 between “sintomas dengue” and “dengue.”

• Download figure
• Open in new tab

Fig 3.

Correlation among dengue cases and Google Trends indices for several dengue-related keywords in Brazil.

In addition, to assess the current situation, we used Google Keyword Planner https://ads.google.com/intl/en/home/tools/keyword-planner/ to find keywords related to “dengue” in Brazil from January to April 2024. By providing the keyword dengue, the tool showed high search interest in “dengue” and “sintomas de dengue.” Consequently, we decided to incorporate Google Trends indices for the keywords “sintomas dengue” and “dengue” into our models.

Results

This section presents the error and uncertainty measures obtained for each method across all states. Results are not shown for the Espírito Santo as this state stopped reporting dengue cases to the federal governement since epidemiology week 16 of 2024 (April 14), resulting in missing case counts from that date onwards. Furthermore, the data for epidemiology week 24 of 2024 (June 9) was not uploaded by InfoDengue, so this week’s comparison was skipped in our analysis. For comparisons requiring this week’s data as the “true value”, the data from epidemiology week 25 of 2024 (June 16) was used.

Tables 1 to 4 present the error measures obtained for each state and nowcasting approach, where red indicates the best performance (smallest error), and blue represents the second-best performance (second smallest error) model. Figure 4 displays boxplots of the differences between the nowcasts and the true number of cases across states. In this figure, models are sorted by the mean of the differences being closest to zero, corresponding to the MAE. Maps with the best nowcasting models across states according to each error metric are shown in Figure 5.

• Download figure
• Open in new tab

Fig 4.

Boxplots of differences between the nowcasts and actual cases in Brazilian states. Models are sorted in ascending order of their mean absolute difference.

• Download figure
• Open in new tab

Fig 5.

Best nowcasting models in different states according to different metrics.

Table 1 presents the RMSE values for prediction errors of each model. Overall, the GT model performed the best among the five approaches, achieving the most accurate predictions in 14 out of 26 states evaluated. The InfoDengue and DCGT models showed comparable performance, excelling in 10 and 2 states, respectively. Conversely, the naive and DC approaches performed the worst. In several states, such as Acre, Pará, Paraíba, Piauí, Rio Grande do Sul, and Sergipe, the GT model significantly outperformed the second-best models (the error of the second-best performing model was more than 1.5 times the error of the best-performing model). The smallest errors were observed in Roraima, all below 100, whereas São Paulo exhibited the largest errors, ranging from approximately 50,000 to 100,000, reflecting the differences in population sizes.

For MAE (Table 2), the overall results were similar to those observed for RMSE. The GT model outperformed the other approaches in 18 out of 26 states, followed by InfoDengue and DCGT, which were better in 7 and 1 states, respectively. The naive approach and the DC model both had the poorest performance, with no states showing an advantage. For Acre, Ceará, Maranhão, Pará, Paraíba, Paraná, and Rio Grande do Sul, the difference between the best and second-best performances was significant.

RMSPE (Table 3) and MAPE (Table 4) values indicate that the top two models remain GT and InfoDengue, as they consistently produce the best nowcasts across all 26 states. For RMSPE, GT leads in 17 states, and InfoDengue in 11 states (both models produce the same RMSPE in Mato Grosso and Tocantins). Moreover, in many states, there is a notable difference between the performance of the best and second-best models. For example, in Alagoas, Amapá, Mato Grosso do Sul, and Rio Grande do Sul, the error of the second-best performing model was more than 2 times the error of the best-performing model. The results are similar for MAPE, with the GT model maintaining its lead in 16 states, followed by InfoDengue in 11 states (they are both the best in Acre).

Table 5 displays the 95% coverage probabilities, and Table 6 the average width of 95% uncertainty intervals obtained for each state and nowcasting approach. In these tables, results for InfoDengue are not included, as nowcasts for this method are provided by municipality level, while our analysis is conducted at the state level. We found that both the DC and GT models exhibited high coverage rates, achieving the highest coverage in 14 different states. In contrast, the DCGT model only achieved the highest coverage in 6 states. While the coverage rates for DC and GT are comparable, the average width of the uncertainty intervals in Table 6 shows that the DC model has a significantly larger average width. In some regions, such as the Distrito Federal and Maranhão, the prediction width for DC can be as much as twice that of GT.

Additionally, there is considerable variation in 95% coverage probabilities across states. For example, in Acre, Amazonas, and Tocantins, all models have high coverage rates, approaching or reaching 1. However, in states like Bahia and Rio Grande do Sul, the coverage rates for all models are only around 20% to 30%. The variation also exists in average width. This variation may be influenced by differences in population size, data quality, or the nature of dengue transmission across states, leading to significant discrepancies in results.

Figure 6 shows the weekly nowcasting results for each method. The green lines represent the reported cases for a given epidemiological week after 10 weeks. We consider these values as the true number of cases to benchmark the models’ performance. Notably, the suspected cases lines are consistently the lowest, serving as a practical “lower bound”, with various models employed to adjust this lower limit. In general, the red lines representing nowcasts obtained with Google Trends closely align with the true values, indicating that the GT model performs particularly well. In contrast, the purple line for InfoDengue nowcasts frequently appear above the true values, suggesting a tendency to overestimate potential cases. Despite this, the InfoDengue nowcasts successfully capture the trends indicated by the true values. Finally, the DCGT and GT nowcasts are generally less likely to accurately capture the trends of the true values.

• Download figure
• Open in new tab

Fig 6.

Predictions of the models used in the study from week 10 to week 23 of 2024 (March 3 to June 2, 2024). The ”True number of cases” line represents the reported cases after 10 weeks, which are used to evaluate the performance of the models.

Discussion

Reporting delays in surveillance data make traditional surveillance systems ineffective for planning and control. In this paper, we compared the usefulness of integrating Google Trends information for dengue nowcasting in Brazil with other analytical approaches that rely solely on reported data. Specifically, we evaluated the error and uncertainty produced by approaches that used only Google Trends information (GT), only dengue case data (DC), and a combination of both (DCGT). In addition, we compared these results with the nowcasting algorithm provided in InfoDengue and a naive approach, where the number of cases in a given week was nowcasted as the number of cases reported in the previous week.

Our study demonstrates the effectiveness of combining Google Trends information with reported case data for nowcasting dengue in Brazil. We show that using the reported number of cases as the nowcast for the following week, as done in the naive approach, is insufficient for real-time monitoring since it significantly underestimates the actual number of dengue cases. Overall, the GT model demonstrates the lowest error across most states, outperforming other models in all metrics, with InfoDengue ranking second. In contrast, neither DCGT nor DC provides significant improvements in nowcasting accuracy. The results, as illustrated through boxplots and accuracy metrics, confirm that incorporating Google Trends data effectively enhances predictive performance. Regarding uncertainty, we lack uncertainty intervals for InfoDengue since the nowcasts are provided at the municipality level. Among the remaining models, GT generally has the narrowest intervals while achieving the highest coverage rate. In contrast, DCGT generates significantly narrower uncertainty intervals compared to DC but slightly lowers the coverage rate. Overall, GT improves predictive accuracy and reduces uncertainty.

In our study, we excluded Espírito Santo from model comparisons because the state ceased reporting dengue cases. However, it is possible to fit an appropriate GT model using historical data and continue nowcasting with Google Trends data, even after reporting has stopped — this approach was actually implemented on the Dengue Tracker website. Another noteworthy case is Rio de Janeiro, where changes in the notification system’s infrastructure and workflow were introduced during the epidemic to achieve a faster response. This shift was not captured by InfoDengue, leading to nowcasts that significantly overestimated the actual number of cases from March to April, whereas GT’s nowcasts were much closer to the true figures.

One limitation of our study is that Google Trends data is inherently biased since not all individuals use Google to search for dengue-related information. This bias may result in underrepresentation of specific populations or regions, potentially affecting the accuracy of our models. Further research is needed to understand how different population groups use Google to search for dengue information, how to select dengue-related keywords that accurately reflect disease transmission, and to develop models that integrate Google Trends and dengue case data in the most effective way.

In this study, we employed statistical models that excluded the most recent weeks of incomplete information to generate dengue nowcasts. Although this approach allowed us to demonstrate the superior performance of approaches using Google Trends data compared to models relying only on reported cases, further work could be done to develop models that utilize incomplete data to further improve predictive accuracy.

In addition, the models used in this study are limited in their ability to detect sudden changes in dengue incidence, since they heavily rely on historical data. This limitation hinders their effectiveness in identifying abrupt outbreaks or sharp increases in cases. Future research will explore more flexible approaches to improve responsiveness.

Additionally, we will investigate the incorporation of variables like climate and socio-economic factors, known to influence dengue transmission, into future models. Moreover, we intend to develop spatial models that allow us to obtain nowcasts at finer geographical resolutions, such as microgregion or municipality levels (39). This enhancement will provide localized nowcasts, enabling more precise public health interventions.

Effective and timely communication of dengue activity levels is crucial for planning and response efforts in public health. To address this need, we also developed Dengue Tracker (https://diseasesurveillance.github.io/dengue-tracker/index.html), a system designed to aggregate, analyze, and visualize dengue case data. This system provides weekly nowcasts at the state level in Brazil, aiding decision-makers and the public in understanding current risk levels. The Dengue Tracker website is updated weekly and features interactive maps and time series plots that dynamically present the latest dengue information across Brazil. By integrating our nowcasting models that use Google Trends information into this platform, the website delivers real-time alerts and trend analyses for better disease prevention and control.

In conclusion, our study presents a promising approach to improving dengue surveillance in Brazil, emphasizing the potential of integrating digital data with traditional epidemiological models. This integration significantly enhances situational awareness for public health authorities and the general public alike, facilitating responses to changes in dengue activity levels, ultimately reducing the impact of dengue and improving the health and well-being of the population.