Nghiên cứu được thực hiện với mục tiêu xác định mối liên quan giữa khí hậu và sốt xuất huyết (SXH), chúng tôi khảo sát các mối quan hệ giữa biến số khí hậu và số bệnh nhân nhập viện ở Cần Thơ từ 2001 đến 2011. Mời các bạn cùng tham khảo đề tài qua bài viết này.
Y Học TP. Hồ Chí Minh * Tập 18 * Phụ bản của Số 6 * 2014 Nghiên cứu Y học MỐI LIÊN QUAN GIỮA SỐT XUẤT HUYẾT VÀ BIẾN ĐỔI KHÍ HẬU TẠI CẦN THƠ 2001‐2011 Nguyễn Phương Toại*, Đặng văn Chính**, Amy Vittor***, Nguyễn Ngọc Huy****, Hồ Hữu Tính** TĨM TẮT Đặt vấn đề: Nhằm nỗ lực xác định mối liên quan giữa khí hậu và sốt xuất huyết (SXH), chúng tơi khảo sát các mối quan hệ giữa biến số khí hậu và số bệnh nhân nhập viện ở Cần Thơ từ 2001 đến 2011. Phương pháp: Số liệu nhập viện hàng tháng và các chỉ số trung bình nhiệt độ, lượng mưa, độ ẩm đuợc thu thập 2001‐2011. Sử dụng phân tích hồi quy và mơ hình Box‐Jenkins theo chuỗi thời gian với cấu tố tự hồi quy ‐ sai phân và trung bình trượt theo mùa (SARIMA) để xác định mối tương quan giữa số ca mắc và yếu tố thời tiết. Kết quả: Chỉ số mắc mới trong thời gian nghiên cứu có xu hướng tăng theo thời gian từ 2001 ‐ 2011. Số ca mắc mới SXH có mối quan hệ ý nghĩa với độ ẩm và lượng mưa. Tuy nhiên, sau khi hiệu chỉnh theo yếu tố mùa, mơ hình chỉ sử dụng yếu tố độ ẩm với độ trễ của chuỗi ca bệnh sau một tháng là phù hợp nhất để dự báo dịch. Kết luận và kiến nghị: Xu hướng dịch SXH vẫn còn tăng từ 2001 – 2011. Độ ẩm và và các chuỗi ca bệnh với độ trễ một tháng là các yếu tố có liên quan thống kê với số ca mắc mới trong thời gian nghiên cứu và là chỉ số quan trọng để dự đốn dịch. Từ khóa: Sốt xuất huyết, chỉ số mắc mới, lượng mưa, độ ẩm, độ trễ, mơ hình SARIMA. ABSTRACT ASSOCIATION BETWEEN DENGUE HOSPITALIZATION AND CLIMATE IN CAN THO, VIETNAM, 2001 – 2011 Nguyen Phuong Toai, Dang Van Chinh, Amy Vittor, Nguyen Ngoc Huy, Ho Huu Tinh * Y Hoc TP. Ho Chi Minh * Vol. 18 ‐ Supplement of No 6‐ 2014: 443 – 450 Background: In an effort to better characterize the relationship between climate and dengue, we examine here the associations between weather variables and dengue hospitalizations in Can Tho between 2001 and 2011. Methods: Monthly data on hospitalized dengue cases and means of temperature, rainfall and humidity were recorded from 2001 to 2011. We used the Box‐Jenkins approach to Seasonal Autoregressive Integrative Moving Average (SARIMA) modeling of time series to assess the association between dengue cases and these weather factors. Result: Trend of DF incidence was increasing between 2001‐ 2011. There was significant relationship between humidity, rainfall and the incidence of dengue case. After adjusting for seasonality, the incidence of dengue fever was significantly associated with humidity with a lag of one month. Conclusions and suggestions: The trend of DF incidence was increasing. Humidity and rainfall was asociated to incidence.These data suggest that humidity (with a one month lag) is an important determinant of dengue hospitalizations. Keywords: Dengue fever, incidence, rainfall, humidity, lag, SARIMA model. ĐẶT VẤN ĐỀ Sốt xuất huyết (SXH) là bệnh truyền nhiễm * do virus Dengue gây ra và lây truyền bởi một số loài muỗi thuộc giống Aedes, chủ yếu là A. Trường Cao đẳng Y tế Cần Thơ **Viện Y tế công cộng Tp.HCM ***Viện đại học Pennsylvania, Philadelphia PA 19104, USA ****Nghiên Cứu Chuyển Đổi Môi Trường và Xã Hội (ISET) Tác giả liên lạc: TS.Nguyễn phương Toại ĐT: 0939766866 Chun Đề Y Tế Cơng Cộng Email: phuongtoai@yahoo.com 443 Nghiên cứu Y học Y Học TP. Hồ Chí Minh * Tập 18 * Phụ bản của Số 6 * 2014 aegypti. Ước tính thế giới hàng năm có khoảng 50‐100 triệu người mắc, trong đó có 500.000 ca nặng với 22.000 ca tử vong. Mặc dù số ca nặng khơng nhiều, nhưng SXH có biến chứng nặng thường rơi vào nhóm trẻ nhỏ nên có thể gây tử vong. Hiện nay vẫn chưa có vắc‐xin phòng bệnh, nên việc loại trừ muỗi và phòng muỗi đốt chính là cách phòng bệnh chủ yếu. lượng việc chứng minh sự liên quan giữa BĐKH và SXH là rất khó khăn. Nhất là ở Việt Nam, do hiện nay cả nước chưa có hệ thống theo dõi diễn biến dịch bệnh với mối liên quan đến các yếu tố thời tiết. Các kế hoạch khống chế dịch vẫn chưa đánh giá tác động của BĐKH, chưa có các chương trình thích ứng để phòng ngừa và ứng phó với dịch bệnh SXH trong bối cảnh BĐKH(2). SXH ở Việt Nam có khuynh hướng gia tăng những năm gần đây. Các số liệu thống kê quốc gia gần đây cho thấy, năm 2008 có 96.451 ca mắc, năm 2009 tăng lên 105.370 ca và năm 2010 số ca được ghi nhận trên cả nước là 128.831. Do đó, SXH ở Việt Nam được xếp trong nhóm 10 bệnh truyền nhiễm có nguy cơ tử vong cao nhất(15). Để xác định giả thiết “các hiện tượng dịch SXH xảy ra ở Cần Thơ có bị tác động của các yếu tố khí hậu bất thường hay khơng?” chúng tơi đã tiến hành khảo sát sự liên quan giữa các yếu tố thời tiết với dịch SXH. Qua đó giúp ngành y tế hoạch định các chiến lược y tế thích hợp, chủ động hơn trong cơng tác phòng chống dịch, tìm ra biện pháp dự báo thích hợp nhằm ứng phó dịch bệnh SXH trong điều kiện BĐKH. Tại Cần Thơ, theo số liệu của Trung tâm Y tế Dự phòng (TTYTDP) TP Cần Thơ, năm 2008 có 1.235 ca mắc SXH, năm 2009 số ca mắc tăng lên 1.357. Sốt xuất huyết ở Cần Thơ chủ yếu xảy ra đối với trẻ em dưới 16 tuổi, phân bố hầu hết trên địa bàn các xã phường, nhất là các khu dân cư nghèo, điều kiện sống còn nhiều khó khăn, mơi trường ln trong trạng thái ơ nhiễm(13). Xét về khuynh hướng chung, thì tổng số ca mắc hàng năm vẫn chưa được khống chế giảm một cách vững chắc(11). Trước đây, SXH ở Cần Thơ thường chỉ tập trung xuất hiện vào mùa mưa, sau đó giảm dần và tăng cao trở lại vào đầu mùa mưa năm sau. Tuy nhiên, những năm gần đây dịch có khuynh hướng diễn ra ở cả hai mùa mưa nắng và có khuynh hướng gia tăng về mức đg bắt đầu gia tăng trong tháng 4 và cao nhất ở tháng 10 sau đó có xu hướng bắt đầu giảm và lập đi lập lại như thế trong thời gian nghiên cứu. 448 Xác định mơ hình ước lượng SARIMA Bảng 2: Ước lượng các tham số của mơ hình SARIMA Biến số độc lập Độ ẩm Mơ hình Tham số Hệ SARIMA số 2,290 AR(1), MAR(1, 12) AR(1) 0,567 MAR(1, 12) 0,421 Sai số Giá trị chuẩn p 1,12 0,042 0,096 0,108 0,005 0,005 Chun Đề Y Tế Cơng Cộng Y Học TP. Hồ Chí Minh * Tập 18 * Phụ bản của Số 6 * 2014 Nghiên cứu Y học Sau khi điều chỉnh theo mùa, kết quả mơ hình phù hợp nhất với chuỗi thời gian là SARIMA, AR(1); MAR(1,12), cho thấy độ ẩm, là yếu tố có liên quan để dự báo dịch SXH với chuỗi tự tưong quan độ trễ 1 tháng. Chỉ số AR(1) là các yếu tố tiên đốn có độ thích hợp với chuỗi tự tuơng quan số ca bệnh mắc mới sốt xuất huyết với độ trễ là 1 tháng MAR(1,2) thể hiện chuỗi dữ liệu đuợc chỉnh theo mùa của 12 tháng. Test chẩn đồn độ phù hợp của mơ hình Qua sử dụng White‐Noise Test cho các mơ hình dự báo thường dùng, cho thấy có tính tương thích cao của mơ hình dự báo ARIMA và SARIMA đối với bố dữ liệu trong nghiên cứu. Đồ thị phần dư của mơ hình dự báo dựa vào độ ẩm khơng còn tự tương quan (nhiễu trắng). 0.00 Cumulative periodogram for phandu 0.20 0.40 0.60 0.80 1.00 Cumulative Periodogram White-Noise Test 0.00 0.10 Bartlett's (B) statistic = 0.20 0.30 Frequency 0.40 0.50 0.72 Prob > B = 0.6745 Biểu đồ 7: White‐Noise Test tính tương thích cao của mơ hình dự báo SARIMA đổi bất thường khí hậu địa phương đặc biệt là Phân bố phần dư của mơ hình dự báo lượng mưa(7), do đó cần được lưu ý khi triển SARIMA, cho thấy mơ hình xây dựng là nhiễu khai hoạt động phòng chống dịch hàng năm trắng, khơng có sự tương quan trong dữ liệu được thực hiện kịp thời và đúng thời điểm hơn. chuỗi thời gian đang nghiên cứu. Nghĩa là mơ hình thích hợp với chuỗi số liệu đuợc phân tích. Qua bộ số liệu nghiên cứu, dịch SXH ở Cần Thơ trong 11 năm, chưa cho thấy dịch diễn ra có BÀN LUẬN tính chu kỳ năm. Trong đó, năm 2002 có số ca Về tổng qt, chưa phát hiện xu hướng của thấp nhất và năm 2004 có số ca cao nhất, chiếm biến đổi khí hậu đối với nhiệt độ, lượng mưa, và gần 1/5 tổng số ca của 11 năm quan sát. Về phát độ ẩm trong suốt 11 năm từ 2001 đến 2011, có triển dịch SXH cho thấy một xu hướng tăng thể do chuỗi số liệu dùng trong phân tích khơng hàng năm có ý nghĩa thống kê. Đặc biệt là các đủ dài để phát hiện ra các xu hướng của BĐKH. năm gần đây, sau 2004 đều duy trì số ca ở mức Tuy nhiên, qua phân tích hồi quy tuyến tính về cao, số ca nặng và tử vong có tăng, là điều cần khuynh hướng riêng theo từng yếu tố thời tiết được quan tâm. của từng tháng trong 11 năm nhận thấy có sự Về yếu tố độ trễ, kết quả phân tích tự tương thay đổi có ý nghĩa thống kê ở lượng mưa trung quan diễn biến số ca mắc mới cho thấy có sự bình của các tháng 10 giảm, nhiệt độ tối đa trung tương quan lớn nhất ở độ trễ 1 tháng.Yếu tố độ bình các tháng 8 tăng (số liệu khơng trình bày). trễ cho thấy các ca bệnh bộc phát ra có thể là hậu Điều này phù hợp với các ghi nhận có sự thay Chun Đề Y Tế Cơng Cộng 449 Nghiên cứu Y học Y Học TP. Hồ Chí Minh * Tập 18 * Phụ bản của Số 6 * 2014 quả của giai đoạn tiềm tàng về phát triển cơn trùng, mầm bệnh trước đó. Tác dụng trễ của khí hậu trên tỷ lệ mắc SXH có thể giải thích bởi những yếu tố khí hậu khơng ảnh hưởng trực tiếp trên số mắc mới nhưng chỉ ảnh hưởng gián tiếp qua tác động của nó trên q trình dịch. Kết quả này về độ trễ khá phù hợp với nhóm nghiên cứu New Caledonia phân tích giữa khí hậu, dịch SXH và A. aegypti từ 1971‐2010 nhằm xây dựng mơ hình dự báo dịch, cũng cho thấy độ trễ so với đỉnh dịch là 1‐2 tháng. Tương quan có ý nghĩa giữa dịch với nhiệt độ, lượng mưa, và độ ẩm(5). Điều này phù hợp với nghiên cứu của Hà Nội phân tích chỉ số mắc mới ở Hà Nội, Việt Nam, từ 1998 – 2009, để xem có sự tăng theo thời gian với các chỉ số thời tiết và tìm thấy độ trễ 1‐2 tháng; ở độ ẩm và tốc độ gió (với độ trễ 4‐5 tháng)(4). SXH ở Cần Thơ thể hiện mang tính mùa, dịch tăng vào mùa mưa và giảm vào tháng nắng. Phân tích hồi quy tuyến tính đơn biến cho thấy, khi lượng mưa tăng thì số ca mắc mới SXH có xu hướng tăng theo(10). Tuy nhiên trong phân tích đa biến SARIMA mối quan hệ này khơng có ý nghĩa thống kê. Về phương pháp dự báo thì phân tích phù hợp với cách dự báo của Brunkard và cộng sự, sử dụng mơ hình tự tương quan phân tích mối tương quan giữa chỉ số mắc mới tuần và nhiệt độ, lượng mưa và hiện tượng El Nino, từ 1995 – 2005. Đã tìm thấy tương quan dương giữa chỉ số mắc mới SXH với các yếu tố thời tiết này(3). Trong phân tích mơ hình dự báo SARIMA thì yếu tố lượng mưa khơng cho dự báo có ý nghĩa trên số ca SXH. Kết quả nầy phù hợp với nghiên cứu của Banu cộng sự, thực hiện thống kê các nghiên cứu SXH và thời tiết ở Châu Á, cho thấy có sự liên quan giữa SXH và lượng mưa nhưng không bền vững, mặc dù lượng mưa tăng quá lại có tương quan âm với dịch ở Thailand, Indonesia, Taiwan(1). Trong khi đó, yếu tố độ ẩm khơng có liên quan tuyến tính với dịch, nhưng khi phân tích các yếu tố thời tiết trong mơ hình SARIMA, cho 450 thấy độ ẩm trung bình có thể làm dự báo tốt hơn cho số ca mắc mới. Vì độ ẩm tương đối liên quan chặt chẽ với nhiệt độ và lượng mưa, mà trị số độ ẩm luôn hiện diện ở mức cao trong năm dù là trong mùa khơ, nên có ảnh hưởng trên số ca mắc mới. Kết quả này cũng phù hợp với nghiên cứu của Hales và cộng sựl(6), của Myriam Gharbi, và cộng sự(10) và nhất là phù hợp với kết quả của nhóm nghiên cứu New Caledonia(5). Mơ hình SARIMA được sử dụng tốt hơn các mơ hình tuyến tính thơng thường vì nó khắc phục tính tương quan trong chuỗi số liệu theo thời gian và tính mùa. KẾT LUẬN Mặc dù có thể hiện sự biến thiên của các yếu tố khí hậu, nhưng chuỗi 11 năm số liệu khí hậu chưa thể hiện rõ khuynh hướng của các yếu tố biến đổi khí hậu. Dịch SXH ở Cần Thơ vẫn còn xu huớng tăng trong thời gian 11 năm nghiên cứu với ý nghĩa thống kê. Các ca mắc mới có mối liên quan, độ ẩm và độ trễ của một tháng của ca bệnh sốt xuất huyết. Trong mùa nắng, mặc dù nhiệt độ tăng cao, nhưng do độ ẩm cao nên mật độ muỗi vẫn duy trì và có khả năng gây dịch. Độ ẩmvà ca bệnh với độ trễ 1 tháng là yếu tố thích hợp nhất được dùng tiên đốn sự gia tăng của ca bệnh sốt xuất huyết ở Cần Thơ. Đề xuất trong phòng chống dịch bệnh Biện pháp can thiệp phải thực hiện quanh năm, nhưng tập trung hơn diệt lăng quăng từ đầu mùa mưa và khi các gia bệnh bắt đầu gia tăng.Vì độ ẩm, lượng mưa và nhiệt độ ảnh hưởng với nhau, do đó, khi có một trong các yếu tố lượng mưa tăng hay độ ẩm tăng thì càng phải chú ý cơng tác diệt lăng quăng trong tháng trước đó, nhằm phát huy tác dụng làm giảm mật độ cơn trùng cho các tháng sau. Sử dụng mơ hình SARIMA trong nghiên cứu có thể góp phần cho giám sát và cho việc dự đốn trong phòng chống SXH với yếu tố độ ẩm, theo thời gian. Chun Đề Y Tế Cơng Cộng Y Học TP. Hồ Chí Minh * Tập 18 * Phụ bản của Số 6 * 2014 Nghiên cứu Y học TÀI LIỆU THAM KHẢO Banu S, Hu W, Hurst C, Tong S (2011) Dengue transmission in the Asia‐Pacific region: impact of climate change and socio‐ economic factors. Trop Med Int Health. 16(5):598‐607. Bộ Y tế (2007) Chương trình tổng thể dự án chiến lược quốc gia về y tế dự phòng và kế hoạch tổng thể 5 năm về phát triển hệ thống chăm sóc sức khỏeViệt Nam. Tr. 78‐94. Brunkard JMCE, Rothenberg SJ (2008) Assessing the roles of temperature, precipitation, and ENSO in dengue re‐ emergence on the Texas‐Mexico border region. Salud Publica Mex. 50:227‐234. Cuong HQ, Hien NT, Duong TN, Phong TV, Cam NN, Farrar J, Nam VS, Thai KT, Horby P (2011) Quantifying the emergence of dengue in Hanoi,Vietnam: 1998 – 2009. PLoS NTD. 5(9):e1322. Descloux E, Mangeas M, Menkes CE, Lengaigne M, Leroy A, Tehei T, Guillaumot L, Teurlai M, Gourinat AC, Benzler J, Pfannstiel A, Grangeon JP, Degallier N, Lamballerie XD (2012) Climate‐based models for understanding and forecasting dengue epidemics. PLoS NTD, 6(2):e1470. Hales SW, de Wet N, Maindonald J, Woodward A. (2002) Potential effect of population and climate changes on global distribution of dengue fever: an empirical model. Lancet. 360:830‐834. Intergovernmental Panel of Climate Change (2007) Climate Change: Impacts, Adaptation and Vulnerability. Cambridge: Cambridge University Press. Pp. 123‐234. Johansson MA, Dominici F, Glass GE,(2009) Local and global effects of climate on dengue transmission in Puerto Rico. PLoS NTD. 3(2):382 LuzPM,MendesBV,CodeỗoCT,StruchinerCJ,GalvaniAP. (2008) Time series analysis of dengue incidence in Rio de Janeiro, Brazil. Am J Trop Med Hyg. 79(6):933‐939. 10 Myriam GPQ, Joël G, Sylvie C, Guy LR, Laurent G, Laurence M (2011) Time series analysis of dengue incidence in Guadeloupe, French West Indies: Forecasting models using climate variables as predictors. BMC Infectious Diseases. 11(166). 11 Sở Y tế Cần Thơ (2011) Số liệu báo cáo công tác Y tế 5 năm 2006 – 2010. Thành phố Cần Thơ. Tr. 3‐4. 12 Thai KT, Anders KL (2011) The role of climate variability and change in the transmission dynamics and geographic distribution of dengue. ExpBiol Med (Maywood). 236(8):944‐ 954. 13 Trung tâm Y tế dự phòng Cần Thơ (2009) Báo cáo hoạt động YTDP năm 2008 – Kế hoạch YTDP năm 2009. Thành phố Cần Thơ. Tr. 12‐13. 14 Viện Pasteur thành phố Hồ Chí Minh (2012) Chẩn đốn sốt Dengue/sốt xuất huyết Dengue. Chun mục Y tế cộng đồng. Viện Pasteur thành phố Hồ Chí Minh. Tr. 2‐6. 15 Viện Vệ sinh dịch tễ trung ương (2010) Báo cáo số liệu thống kê sốt xuất huyết dengue. Hà Nội. Tr. 5‐8. Ngày nhận bài báo: 9/5/2014 Ngày phản biện nhận xét bài báo: 15/6/2014 Ngày bài báo được đăng: 14/11/2014 Chuyên Đề Y Tế Công Cộng 451 ... dịch bệnh SXH trong điều kiện BĐKH. Tại Cần Thơ, theo số liệu của Trung tâm Y tế Dự phòng (TTYTDP) TP Cần Thơ, năm 2008 có 1.235 ca mắc SXH, năm 2009 số ca mắc tăng lên 1.357. Sốt xuất huyết ở Cần Thơ chủ yếu xảy ra ... lượng việc chứng minh sự liên quan giữa BĐKH và SXH là rất khó khăn. Nhất là ở Việt Nam, do hiện nay cả nước chưa có hệ thống theo dõi diễn biến dịch bệnh với mối liên quan đến các yếu tố thời tiết. Các kế hoạch khống chế dịch vẫn chưa ...Brunkard và cộng sự, sử dụng mơ hình tự tương quan phân tích mối tương quan giữa chỉ số mắc mới tuần và nhiệt độ, lượng mưa và hiện tượng El Nino, từ 1995 – 2005. Đã tìm thấy tương quan dương giữa chỉ số mắc mới SXH với các yếu tố