Saira Karimi1, Naeem Akhtar1, Olodo Katiana2, Sajjad Haider3, Burhan Ahmed4,5, Mujeeb Khan1, Muhammad Umer1
1Department of Pathology, Rawalpindi Medical University, Rawalpindi; 2Université Cheikh Anta Diop de Dakar (UCAD), Senegal; 3Pakistan Meteorological Department, Islamabad; 4Flood Forecasting Division, Pakistan Meteorological Department, Lahore, 5Infection Disease Unit, Holy Family Hospital, Rawalpindi
Objective: Climate warming and vector born infectious diseases are a growing concern under Pakistan's recent climate warming scenario. This retrospective study considers a three-year (2019-2023) comprehensive analysis of dengue virus cases reported at Holy Family Hospital and anticipate the area’s future dengue dynamics through climate modeling.
Methodology: Patient demographic features, age, location, gender, lab profiles for Dengue Fever (DF), Dengue Shock Syndrome (DSS), and Dengue Hemorrhagic Fever (DHF), were collected and analyzed to understand the prevalence and hotspots of dengue virus.
Results: The study revealed that DF was more frequent among individuals above 50 age group, emphasizing the age-dependent nature of dengue vulnerability. Furthermore, our findings highlighted a gender-based vulnerability, indicating that males were more prone to DF (73.4%). Based on these findings, we predicted the impact of climate change (Temperature) on dengue transmission suitable days (DTSD). The proposed predictive model incorporates the baseline (2019-2023) and future (2025–2035, 2041–2070, and 2071–2099) periods under Representative Concentration Pathway (RCP4.5 and RCP8.5) scenarios. CMIP5 models, downscaled and bias-corrected with the quantile delta mapping technique, were employed to project the potential spatiotemporal shifts and DSTD due to climate change.
Conclusion: Our predictive analysis contributes to proactive public health measures by anticipating and preparing for the evolving dynamics of dengue in the context of a changing climate.
Dengue fever is caused by an arbovirus belonging to flaviviruses mainly transmitted via mosquito.1 It is one of the most rapidly spreading diseases worldwide and is becoming a public health concern affecting more than 100 million humans yearly. Dengue also causes 20 to 25,000 deaths, mainly affecting children, and exists in more than 100 countries. At the same time, the mortality rate is 1-2%. These are enveloped single-stranded RNA viruses more prevalent in tropical and subtropical areas due to optimum climatic conditions. Infection with DENV results in varying degrees of pathological conditions, ranging from mild asymptomatic dengue fever (DF) to severe dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS) which may turn fatal.1 In some cases, DF may accompany bleeding complications such as gingival bleeding, epistaxis, gastrointestinal bleeding, hematuria, and menorrhagia (in the case of women). Dengue hemorrhagic fever (DHF) is characterized by symptoms of DF along with thrombocytopenia, hemorrhagic manifestations, and plasma leakage.2 The epidemiological triangle of dengue consists of host, pathogen, and mosquito vectors (including Ae. aegypti and Ae. Albopictus) along with their interactions in the environment.3 Climate influences the propagation of dengue as the virus has to complete part of its development in the vectors responsible for disease transmission that transmit the disease. The major vectors are mosquitoes such as Aedes egypti or albopictus whose lifecycle is dependent on ambient temperature and rainfall.4 The climate is constantly changing which increases the chances to expand the geographical distribution of several mosquito-borne diseases 5. The evidence from around the areas with frequent dengue outbreaks indicate a close association between increasing incidence of dengue with rainfall, relative humidity, and temperature.6 Several studies have indicated an increase in the transmission potential of dengue fever under climate change based on the association between meteorological factors and disease transmission. Meteorological factors have an intricate influence on disease transmissions, and patterns of disease epidemics, such as epidemic size and peak time, may be subject to change under global warming.7 Climate factors such as precipitation also significantly affect the life cycle of dengue mosquitoes thus, these factors need to be included in mathematical models of dengue transmission. Furthermore. Similarly, rainfall has been shown to influence larval density and larval population size of vector of dengue.4,6,8,9 A substantial increase in the occurrence of dengue has been shown to increase globally in the last two decades (Dengue Worldwide Overview, 2022). Since 2023 an increase in dengue transmission occurred due to urbanization, the distribution pattern of vectors, and changing climate patterns. Almost five million cases and more than 5000 dengue-related deaths were reported in more than 80 countries/territories and WHO regions such as Africa, Americas, Southeast Asia, Wester, Pacific, and Eastern Mediterranean Regions globally (Dengue – Pakistan WHO (2022) including vulnerable countries such as Afghanistan, Pakistan, Sudan, Somalia, and Yemen.10 Pakistan (n= 20 072), Saudi Arabia, and Oman have reported the highest number of confirmed cases. (Dengue - Pakistan, WHO,.2023; Xu et al., 2017). Asia represents about 70% of the world’s dengue disease burden11 and has significant negative financial and sociological consequences with the potential to considerably impede the development of economies, politics, and society. In the past decade, the prevalence of dengue fever has surged globally (Khan et al., 2022.) At present, dengue is widespread in over 129 countries and cases are rising constantly over time. It has become imperative to delve into policies that can devise the dengue manifestation, detection, and future impacts at the local level.

Our research area, consisting of Rawalpindi and Islamabad, is a major metropolitan region situated in the northern part of Pakistan, specifically in the Punjab province. The area's geography is varied, with a mix of plains and the foothills of the Himalayas. The terrain of Islamabad is mostly flat, particularly in its central areas, but it gradually becomes hillier towards the Margalla Hills. Despite being near the Soan River, the cities are not situated directly on its banks. A search strategy was formulated to comprise diagnoses conducted at Holy Family Hospital during the years 2019 to 2023, with a focus on positive cases of DF, DHF, and DSS. The inclusion criteria encompassed studies that reported on the prevalence of dengue fever, employed NS1 antigen, IgG, and IgM detection methods, and offered sufficient demographic information. Data extraction was executed systematically to capture pertinent information from the selected studies. Variables of interest included the number of dengue cases, gender distribution, age distribution, district-wise distribution, and diagnostic techniques employed. The data extraction was performed independently, and any discrepancies were reconciled through consensus. The extracted data were subjected to a robust statistical analysis to unravel trends and patterns in the prevalence of dengue fever. Categorization of data was performed based on gender, age, and district, with further exploration of the abundance of cases. Descriptive statistics, chi-square testing, and logistic regression analysis using R programming were applied. Gender-specific prevalence rates were examined to identify any notable variations. Pie charts were employed to illustrate the overall distribution of dengue cases, emphasizing the proportion of each infection type. Mapping tools were employed to visualize the data of age, gender, and frequency of disease at different locations of Rawalpindi to estimate the total disease burden and its socio-economic details employed on in our study focusing on identifying suitable transmission days for dengue fever, we leveraged temperature data from both historical records and future projections obtained from the RCA4 CORDEX South Asia dataset. This dataset enabled us to examine temperature patterns over past periods as well as project potential future scenarios, providing essential insights into the environmental conditions conducive to dengue transmission. By utilizing this comprehensive temperature data, which covers historical periods and future projections, we aimed to better understand the relationship between temperature fluctuations and the suitability of conditions for dengue transmission. This analysis allowed us to identify periods with optimal temperature conditions for the proliferation of dengue vectors, aiding in the assessment of dengue risk and the development of targeted intervention strategies. In our research, we employed a combination of statistical bias correction and downscaling techniques to enhance the accuracy and resolution of climate data used in our study. Statistical bias correction involves adjusting climate model outputs to reduce systematic errors or biases relative to observations. This process ensures that the simulated climate variables align more closely with historical data, improving the reliability of future projections. Additionally, we utilized downscaling methods to refine the spatial and temporal resolution of climate model outputs. Downscaling involves generating finer-scale climate information from coarser-resolution global climate models or reanalysis datasets. By downscaling the data, we were able to capture local-scale variations in climate variables, which are essential for assessing the impacts of climate change on regional and local scales. The combined approach of statistical bias correction and downscaling allowed us to produce high-quality climate data suitable for our specific research objectives. These refined climate datasets provided a more accurate representation of past climate conditions and enabled robust projections of future climate scenarios, facilitating a comprehensive analysis of the potential impacts of climate change on dengue transmission dynamics. To empirically adjust variable biases originating from regional climate model simulations we used the widely used quantile mapping approach (Gudmundsson et al., 2012). We identified shape distributions of different quantile mapping methods. Afterwards, investigation of the most appropriate distribution of the quantile mapping based on the least biases was executed (Burhan and Rasul, 2018). The selected statistical transformation strived to acquire a function that mapped a modeled parameter Pm such that its post-processed distribution equaled the distribution of the observed parameter Po. Adhering to Piani et al. (2010), we deployed a generic transformation
expressed as in Eqn. (1) and (2). In the following mathematical expressions for the quantile mapping approach, the parameters Po and Pm indicate observed and modeled climatic variables, respectively
Po=c(Pm) (1)
Po= Fo-1 (Fm (Pm)) (2)
where Fm is the Cumulative Distribution Function (CDF) of Pm and Fo−1 is the inverse CDF (or quantile function) linking to Po.
Disease burden and distribution of Dengue fever (DF), Dengue Hemorrhagic Fever (DHF), and dengue shock syndrome (DSS) in different areas of Rawalpindi shows that in central pindi the abundance of dengue fever was more as compared to the other parts of the district (Table I, Figure 2).


Table I: Month wise distribution of dengue-positive cases in holy Family hospital from 2019-2023 |
|||||||||||||
Year |
Jan |
Feb |
Mar |
Apr |
May |
Jun |
Jul |
Aug |
Sep |
Oct |
Nov |
Dec |
Total |
2019 |
0 |
0 |
0 |
0 |
1 |
2 |
27 |
684 |
4686 |
5581 |
961 |
0 |
1192 |
2022 |
0 |
0 |
2 |
0 |
0 |
5 |
9 |
420 |
1912 |
2106 |
565 |
20 |
5039 |
2023 |
0 |
0 |
1 |
1 |
3 |
9 |
25 |
148 |
185 |
- |
- |
- |
372 |





According to multiple studies, the temperature in Islamabad-Rawalpindi is likely to increase significantly under the high emissions scenario RCP8.5 in comparison to the moderate emissions pathway RCP4.5.18 Studies have found that Islamabad and its surroundings are projected to experience a notable warming trend under RCP8.5, which represents a continued high greenhouse gas emission trajectory.19 It was studied that climate projections for Pakistan under different emission scenarios and highlighted the substantial temperature increases expected by the end of the century under RCP8.5.19,20 Similarly, Islamabad's climate projections emphasized the heightened temperature rise under RCP8.5 compared to RCP4.5.17,21-23 These studies underscore the significance of accounting for emission scenarios when evaluating future climate impacts on Islamabad. The high emissions trajectory of RCP8.5 indicates a more severe warming trend, which poses greater challenges for the city in terms of heat stress, water availability, and ecosystem resilience. In contrast, RCP4.5 represents a pathway with moderated emissions due to mitigation efforts, resulting in a comparatively lower rate of temperature increase.24 Although this still leads to warming and associated impacts, the magnitude and severity are expected to be less pronounced compared to RCP8.5.25 The chart depicts the annual temperature cycles in Islamabad-Rawalpindi for historical and future periods, assuming scenario RCP8.5. The pattern of each line is similar, with temperatures increasing from January, reaching a peak in June, and then decreasing towards December. However, the temperatures in future periods appear to be higher than in the past, indicating an overall increase in mean temperatures. This graph is a typical representation of climate change projections, with RCP8.5 being just one of the possible scenarios for future climate. This particular scenario assumes the highest greenhouse gas emissions compared to other pathways. Mean temperatures will likely increase even more under RCP8.5 than under RCP4.5, exacerbating the effects of climate change in the studied area. All of the future projections show higher temperature ranges than historical observations.9,17,19-21
Climate projections for Islamabad indicate that the shoulder seasons, which include spring and autumn, are likely to experience higher mean temperature changes in comparison to the summer season.26 Temperature trends across different seasons were analyzed by Amnuaylojaroen et al., 2022; Hausfather & Peters, 2020; Sharif et al., 2024, their finding suggested that while summer temperatures are expected to increase, the rate of temperature rise during the shoulder seasons, particularly in spring and autumn, is projected to be even higher. Ali et al. (2020) also found significant temperature changes anticipated during the shoulder seasons, with potentially greater impacts on local ecosystems, agriculture, and human health.26
Research studies indicate that the climate change effects on the transmission of dengue fever in Islamabad-Rawalpindi are alarming.29 The projected periods under both RCP8.5 and RCP4.5 scenarios are expected to significantly increase the season favorable for dengue transmission. Analysis of the impact of various emission scenarios on dengue suitability in Islamabad and found that the season conducive to dengue transmission is projected to increase substantially in the future.16,17,30 Similarly, when examined the impact of climate change on dengue transmission dynamics it was concluded that regardless of the emissions trajectory, the period characterized by suitable temperature conditions for dengue transmission is expected to expand significantly. This expansion of the dengue transmission period poses significant challenges for dengue control and public health efforts.30
Recent research suggests that the incidence of DSTD may rise during the transitioning months, while decreasing in the summer.17 When projected to analyze the link between climate change and dengue fever it was indicate that the number of days with temperatures suitable for dengue transmission is expected to increase notably in spring and autumn.5,16 Similarly, analyzed a significant decrease in the frequency of dengue-suitable temperature days during the summer, but a rise in spring and autumn. Recent analyses have revealed that days with suitable temperature for dengue fever in Islamabad have historically followed a unimodal pattern, meaning there was a single peak during a specific season. However, future projections indicate a shift towards a bimodal pattern, which means that there will be increased occurrences of the disease during both spring and autumn. This shift is caused by the historical unimodal pattern converting to bimodal, which means that both transitioning seasons are expected to exhibit higher occurrences of the disease. This shift can create an alarming situation, leading to a surge in patients at the Holy Family Hospital.
This study emphasizes the importance of considering seasonal variations and the potential impacts of climate change on dengue transmission dynamics. The shift towards a bimodal pattern indicates a significant change in the timing and distribution of dengue-suitable temperature days, which has implications for dengue control strategies and public health interventions in Islamabad-Rawalpindi. Adapting surveillance and vector control measures to account for the increased occurrences during spring and autumn will be crucial for effectively managing dengue transmission in the future. Understanding these shifting patterns can inform policymakers and public health authorities in developing targeted strategies to mitigate the impacts of climate change on dengue incidence in Islamabad and ensure the health and well-being of the population.

An Official Publication of
Islamabad Medical & Dental College
Volume 13 Issue 2
Saira Karimi
Email:
saira.karimi@rmur.edu.pk
Cite this article. Karimi S, Akhtar N, Katiana O, Haider S, Ahmed B, Khan M, Umer M Decoding Dengue: A Comprehensive Analysis of Cases at Holy Family Hospital (2019–2023) and Anticipating Pakistan's Future Dengue Dynamics under Climate Change. J Islamabad Med Dental Coll. 2024; 13(2): 345-354. DOI: https://doi.org/10.35787/jimdc.v13i2.1262