Using satellite data to study the relationship between rainfall and diarrheal diseases in a Southwestern Amazon basin

O uso de dados de satélite para estudar a relação entre chuva e doenças diarreicas em uma bacia na Amazônia Sul-Ocidental

Paula Andrea Morelli Fonseca Sandra de Souza Hacon Vera Lúcia Reis Duarte Costa Irving Foster Brown About the authors

Abstract

The North region is the second region in Brazil with the highest incidence rate of diarrheal diseases in children under 5 years old. The aim of this study was to investigate the relationship between rainfall and water level during the rainy season principally with the incidence rate of this disease in a southwestern Amazon basin. Rainfall estimates and the water level were correlated and both of them were correlated with the diarrheal incidence rate. For the Alto Acre region, 2 to 3 days’ time-lag is the best interval to observe the impact of the rainfall in the water level (R = 0.35). In the Lower Acre region this time-lag increased (4 days) with a reduction in the correlation value was found. The correlation between rainfall and diarrheal disease was better in the Lower Acre region (Acrelândia, R = 0.7) and rainfall upstream of the city. Between water level and diarrheal disease, the best results were found for the Brasiléia gauging station (Brasiléia, R = 0.3; Epitaciolândia, R = 0.5). This study’s results may support planning and financial resources allocation to prioritize actions for local Civil Defense and health care services before, during and after the rainy season.

Infant diarrhea; Rainy season; Remote sensing; Environmental health surveillance; Brazilian Amazon

Resumo

A região Norte é a segunda no Brasil com a maior taxa de incidência de doenças diarreicas em crianças menores de 5 anos. O objetivo deste estudo foi investigar a relação entre chuva e nível do rio, principalmente durante a estação chuvosa, com a taxa de incidência da referida doença em uma bacia no sudoeste da Amazônia. Estimativas de chuva e nível do rio foram correlacionadas e ambos correlacionados com a taxa de incidência da diarreia. Para a região do Alto Acre, 2 a 3 dias de defasagem é o melhor intervalo para observar o impacto da chuva no nível do rio (R = 0.35). Na região do Baixo Acre essa defasagem aumentou (4 dias) com redução na correlação. A correlação entre chuva e doenças diarreicas foi melhor na região do Baixo Acre (Acrelândia, R = 0.7) e a chuva rio acima da cidade. Entre o nível do rio e as doenças diarreicas, os melhores resultados foram encontrados para a estação de Brasiléia (casos em Brasiléia, R = 0.3 e Epitaciolândia, R = 0.5). Os resultados deste estudo podem dar apoio ao planejamento e alocação de recursos financeiros para priorizar ações para Defesa Civil e serviços de saúde antes, durante e depois da estação chuvosa.

Diarreia infantil; Estação chuvosa; Sensoriamento remoto; Vigilância em saúde ambiental; Amazônia Brasileira

Introduction

Under normal climate conditions, diarrheal diseases and pneumonia are the major cause of mortality among children under five years, especially in poorer countries, which represent 29% of total deceased in an annual base worldwide11. World Health Organization (WHO). World Health Statistics 2013. [accessed 2014 Dec 1]. Available at: http://www.who.int/gho/publications/world_health_statistics/2013/en/
http://www.who.int/gho/publications/worl...
. However, these diseases have been widely reported after geophysical disasters, and hydrometeorological events such as floods22. Floret N, Viel JF, Mauny F, Hoen B, Piarroux R. Negligible risk for epidemics after geophysical disasters. Emerg Infect Dis 2006; 12(4):543.. According to statistics obtained from the Mortality Information System in Brazil33. Ministério da Saúde, Secretaria de Vigilância em Saúde MS/SVS/SIM. Sistema de Informações de Mortalidade (SIM). [accessed 2013 Dec 1]. Available at: http://tabnet.datasus.gov.br/cgi/deftohtm.exe?idb2012/c06.def
http://tabnet.datasus.gov.br/cgi/deftoht...
under the Informatics’ Department of the Health System (DATASUS) from 1996 to 2012, diarrhea and gastroenteritis arise as one of the main causes of children mortality in Brazil. In Northeast region, this disease is the 4th cause of death (6.2 cases per 1,000 children). Nevertheless it must be considered the higher number of children in this specific region. In the North and Central-West regions diarrhea is the 8th main cause of infant mortality with an incidence rate of 2.4 and 1.7 per 1,000 children respectively.

Several studies have been trying to establish a link between the occurrence of water borne infectious diseases (as diarrhea) and climatic and hydrological patterns44. Epstein PR. Climate Change and emerging infectious diseases. Microbes Infect 2001; 3(9):747-754.,55. Kovats RS, Bouma MJ, Hajat S, Worral E, Haines A. El Niño and health. Lancet 2003; 362(9394):1481-1489.. Some studies analyzed data related to climate extremes and flood episodes to discuss the impacts on health66. Alderman K, Turner LR, Tong S. Floods and human health: A systematic review. Environ Int 2012; 47:37-47.,77. Kozama S, Aizawa T, Watanabe T, Ranjan P, Gunawardhana L, Amano A. A quantitative risk assessment of waterborne infectious disease in the inundation area of a tropical monsoon region. Sustainable Sciences 2012; 7(1):45-54.. Some of the water diseases related to climatic extreme events are cholera, leptospirosis, hepatitis, bacillary dysentery and typhoid fever. These extremes episodes lead to a decrease in hygiene habits due to shortage of fresh water and the damage caused by the lack of sewage system. In addition, the high number of people accommodated in shelters during the flood period increases altogether the incidence of these diseases.

Since 2012, southwestern Amazon has been severely affected by extreme floods88. Espinoza JC, Ronchail J, Frappart F, Lavado W, Santini W, Guyot JL. The Major Floods in the Amazonas River and Tributaries (Western Amazon Basin) during the 1970-2012 Period: A Focus on the 2012 Flood. J Hydrometeorol 2013; 14:1000-1008.,99. Espinoza JC, Marengo JA, Ronchail J, Carpio JM, Flores LN, Guoyt JL. The extreme 2014 flood in the South-western Amazon basin: the role of tropical-subtropical South Atlantic SST gradient. Environ Res Lett 2014; 9:124007.. In 2014, Rio Branco’s mayor, the capital of Acre State, declared the emergency state alongside with Peruvian, Bolivian and other Brazilian cities. More than 2.000 houses were affected in the capital of Acre. Recently, during the 1st semester of 2015 another flood event affected this region and again exposed the Amazonian population to the impact of waterborne diseases1010. Portal de Notícias da Globo (G1). União reconhece calamidade em Rio Branco; 300 militares dão apoio. [accessed 2015 Mar 1]. Available at: http://g1.globo.com/ac/acre/noticia/2015/03/uniao-reconhece-calamidade-em-rio-branco-300-militares-serao-enviados.html
http://g1.globo.com/ac/acre/noticia/2015...
.

The number of meteorological stations and rain gauges is very limited in Amazonia considering the nature, scale, dynamic and microphysics involved in the rainfall events in the tropics, usually caused by local convection1111. Wallace JM, Hobbs PV. Cloud Microphysics. In: Wallace JM, Hobbs PV. Atmospheric science: an introductory survey. Burlington: Academic press; 2006. p. 209-263. (Vol. 92).. Satellite data provides a useful alternative to allow filling the absence of locally measured data1212. Santos e Silva CM, Lúcio OS, Spyrides MHC. Distribuição espacial da precipitação sobre o Rio Grande do Norte: estimativa via satélites e medidas por pluviômetros. Rev. Bras. Meteorol. 2012; 27(3):337-346.

13. Collischonn B, Allasia D, Collischonn W, Tucci CEM. Desempenho do satélite TRMM na estimativa de precipitação sobre a bacia do Paraguai Superior. Rev. Bras. Cartograf. 2007; 59(1):93-99.

14. Zulkafli Z, Buytaert W, Onof C, Manz B, Tarnavsky E, Lavado W, Guyot JL. A Comparative Performance Analysis of TRMM 3B42 (TMPA) Versions 6 and 7 for Hydrological Applications over Andean-Amazon River Basins. J Hydrometeorol. 2014; 15:581-592.
-1515. Kummerow C, Simpson J, Thiele O, Barnes W, Chang ATC, Stocker E, Adler RF, Hou A, Kakar R, Wentz F, Ashcroft P, Kozu T, Hong Y, Okamoto K, Iguchi T, Kuroiwa H, Haddad EIZ, Huffman G, Ferrier B, Olson WS, Zipser E, Smith EA, Wilheit TT, North G, Krishnamurti T, Nakamura K. The status of the Tropical Rainfall Measuring Mission (TRMM) after two years in orbit. J Appl. Meteorol. 2000; 39(12):1965-1982.. Therefore in Acre basin there is no previous work validating any satellite rainfall estimates. Thus, this work aims to a) evaluate the satellite rainfall estimates for the Acre Basin; b) analyze the correlation between these estimates and: 1) water level in three cities along river Acre, 2) Investigate the incidence rates of infant diarrheal diseases across all municipalities in the river Acre Basin.

Materials and methods

Study Area

The Acre State, which is located at the southwestern Amazon, is characterized by a seasonal rainfall regime monsoon77. Kozama S, Aizawa T, Watanabe T, Ranjan P, Gunawardhana L, Amano A. A quantitative risk assessment of waterborne infectious disease in the inundation area of a tropical monsoon region. Sustainable Sciences 2012; 7(1):45-54.,1616. Marengo JA, Liebmann B, Grimm AM, Misra V, Silva Dias PL, Cavalcanti IFA, Carvalho LMV, Berbery EH, Ambrizzi T, Vera CS, Saulo AC, Nogues-Paegle J, Zipser E, Seth A, Alves LM. Recent developments on the South American monsoon system. Int J Climatol 2012; 32(1):1-21., with an annual rainfall approx. 1,960 mm1717. Duarte AF. Aspectos da Climatologia do Acre, Brasil, com base no intervalo 1971-2000. Rev. Bras. Meteorol. 2006; 21(3b):308-317.. It is in the border with the Peruvian and Bolivian departments of Madre de Dios and Pando respectively, as shown in Figure 1 a. The Acre basin, which is in the eastern side of Acre State, comprises approximately 30,000km22. Floret N, Viel JF, Mauny F, Hoen B, Piarroux R. Negligible risk for epidemics after geophysical disasters. Emerg Infect Dis 2006; 12(4):543.,1818. Ministério Público do Acre (MPE/AC). Caracterização Socioambiental das Bacias Hidrográficas do Estado do Acre. Rio Branco: MPE/AC; 2009.. This basin is basically divided in Xapuri and Riozinho do Rola sub-basins. It covers 11 municipalities and they are distributed in two geopolitical regions Alto (Upper) and Baixo (Lower) Acre. Rio Branco is the Acre’s state capital localized downstream Riozinho do Rola River that has a crucial role on the flood process in this city. Approximately 65% of the ~770,000 Acre State population lives in these regions and 6% of them are children under 5 years old, half of this percentage lives in the east of Acre1919. Instituto Brasileiro de Geografia e Estatística (IBGE). [cited 2013 Dec 1]. Available at: http://www.ibge.gov.br/home/
http://www.ibge.gov.br/home/...
.

Figure 1
A - Municipalities within the Acre River Basin using the code number available in the shape file provided by UCEGEO/AC. The municipalities are part of 2 Acre Regions: Alto Acre (Upper Basin):16 – Assis Brasil, 17 – Brasiléia, 18 – Epitaciolândia, 19 – Xapuri; and Baixo Acre (Lower Basin):4 – Rio Branco, 20 - Capixaba, 8 – Bujari, 5 – Porto Acre, 2 – Senador Guiomard, 21 – Plácido de Castro, 22 – Acrelândia. Black dots represent the geographical locations where the water level series were obtained. Riozinho do Rola and Xapuri sub-basins are delimited in order to highlight the proportion and importance each one has when compared to the total extension of Acre Basin. B - Accumulated rainfall for the NOV-APR semester from 2006 to 2013 split between observed, estimated, and spatially between Upper and Lower Acre; C - Average accumulated rainfall in Upper and Lower Acre regions considering the geographical locations of the rain gauges. Values from NOV 2006 to APR 2013; D - Daily mean of the water level for each municipality, the values were calculated for specific year intervals, according to what was available for each location; E - Incidence rate for infant diarrhea diseases in Acre River basin municipalities from Jan/2000 to Dec/2012.

The river water level data sources

Water level daily data series were obtained from the National Water Agency2020. Agência Nacional de Águas (ANA). HidroWeb: Sistema de Informações Hidrológicas. [cited 2013 Dec 1]. Available at: http://hidroweb.ana.gov.br/
http://hidroweb.ana.gov.br/...
for three gauging station: in Brasiléia and Xapuri, both in the Xapuri sub basin, and data series available to 1998-2006 and 2001-2012 period, respectively, and in Rio Branco, Riozinho do Rola sub basin, with data available for 1998-2012 period.

Difficulties regarding maintenance of the equipment (station or rain gauges) and logistic problems are some of the reason why the data from some of the stations are not updated.

Rainfall observed data source

This study used daily records of the 12 rain gauges from AcreBioClima Project2121. Grupo de Estudos e Serviços Ambientais (AcreBioClima). [cited 2013 Dec 1]. Available at: http://acrebioclima.net/
http://acrebioclima.net/...
for the 2006-2013 period (Table 1).The selection was based on the availability of data during the period of November to April. This criteria was applied because 83% of the total annual rainfall is recorded during this semester1717. Duarte AF. Aspectos da Climatologia do Acre, Brasil, com base no intervalo 1971-2000. Rev. Bras. Meteorol. 2006; 21(3b):308-317..

These rain gauges were split into Upper and Lower Acre (six for each) with a balanced distribution of the available stations in order to show the series behavior.

Table 1
Rain gauges with data available for the NOV-APR semester to 2006-2013.

Rainfall estimated data source

The Tropical Rainfall Measuring Mission2222. Tropical Rainfall Measuring Mission Database (TRMM). [cited 2013 Dec 1]. Available at: ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/Derived_Products/3B42_V7/Daily/
ftp://disc2.nascom.nasa.gov/data/TRMM/Gr...
covers the latitudes from 50oN to 50oS and provides information on a 3-hour time through the 3B42RT product. In this work it was used a derived version of the 3B42 data provided in a daily basis available from 1998-2015-Feb1515. Kummerow C, Simpson J, Thiele O, Barnes W, Chang ATC, Stocker E, Adler RF, Hou A, Kakar R, Wentz F, Ashcroft P, Kozu T, Hong Y, Okamoto K, Iguchi T, Kuroiwa H, Haddad EIZ, Huffman G, Ferrier B, Olson WS, Zipser E, Smith EA, Wilheit TT, North G, Krishnamurti T, Nakamura K. The status of the Tropical Rainfall Measuring Mission (TRMM) after two years in orbit. J Appl. Meteorol. 2000; 39(12):1965-1982.,2323. Huffman GJ, Adler R, Bolvin D, Nelkin E. The TRMM Multi-Satellite Precipitation Analysis (TMPA). In: Gebremichael M, Hossain F, editors. Satellite Rainfall Applications for Surface Hydrology. New York: Springer; 2010. p. 3-22.. TRMM’s spatial resolution is 0.25o x 0.25o degrees and it corresponds to grid cells of 625 km2 approximately.

Infant Diarrheal Diseases Incidence Cases

The statistic of infant diarrheal diseases from the Acre State Health Secretary was obtained for the 2000-2012 period on a monthly basis. See PULSE-Brazil website2424. Impact of Climate Extremes on Ecosystem and Human Health in Brazil (PULSE-Brazil). [cited 2013 dec 1]. Available at: http://www.pulse-brasil.org/tool/
http://www.pulse-brasil.org/tool/...
in order to access this data series.

Our focus was children under 5 years old and the data were extracted from the Acute Diarrhea Daily Monitoring System also known by the acronym SIVEP-DDA. The dataset was used to calculate the incidence rate for each municipality. According to World Health Organization (WHO) diarrheal diseases is classified in the International Classification of Diseases (ICD) under code 10, in Chapter I, and it covers the A00-A09 subtypes, which range from cholera to gastroenteritis and colitis of infectious and unspecified origin. In other words, it includes the three main causes of the disease virus, protozoa and bacteria2525. World Health Organization. International Statistical Classification of Diseases and Related Health Problems 10th Revision (ICD-10)-2015-WHO. [cited 2015 Aug 1]. Available at: http://apps.who.int/classifications/icd10/browse/2015/en#/I
http://apps.who.int/classifications/icd1...
.

Satellite Rainfall Validation

Rainfall observed data series were used to evaluate and validate the satellite data for this geographic location to which TRMM data had never been tested before.

Rainfall estimates from TRMM were extracted using the AcreBioClima geographic coordinates for the stations available in the Acre basin and both datasets were compared against each other through a contingency table2626. Wilks DS. Nonprobabilistic Forecasts of Discrete Predictands. In: Wilks DS. Statistical Methods in the Atmospheric Sciences. 2nd ed. Burlington: Academic Press; 2006. p. 260-275.,2727. Warner TT. Verification Model. In: Warner TT. Numerical weather and climate prediction. Cambridge: Cambridge University Press; 2011. p. 294-320., as can be found in Wilks (2011, p. 261)2626. Wilks DS. Nonprobabilistic Forecasts of Discrete Predictands. In: Wilks DS. Statistical Methods in the Atmospheric Sciences. 2nd ed. Burlington: Academic Press; 2006. p. 260-275.. This verification procedure is essential to validate the quality and coherence of satellite-retrieved climate variables in the context of a specific spatial and seasonal/temporal scale of analysis. Finally, validated estimates are more comparable to other studies and products.

From the contingency table it was possible to calculate the following indexes: Proportion Correct (PC), Critical Success Index (CSI), Bias, False Alarm Ratio (FAR), Hit rate (H) and Probability of False Detection (POFD). The formulas for each of them are described in Equation 1-6.

Equation 1: Proportion correct (PC)

Where,

a = rainfall events both observed and estimated

d = rainfall events neither observed nor estimated

n = number of days analyzed (n), or sample size.

Equation 2: Critical success index (CSI)

Where,

b = total rainfall events observed and not estimated

c = rainfall events only observed.

Equation 3: Ratio of the estimated and precipitated rainfall events.

Equation 4: False Alarm Ratio

Equation 5: Hit rate

Equation 6: Probability of false detection.

Rainfall estimates and water level were correlated assuming that the daily estimates would directly influence the water level in the three gauging station chosen and a shorter or longer time lag would be better to represent that influence. For the correlation analysis between rainfall and diseases, it was considered that the monthly accumulated rainfall in each pixel within a municipality have an influence in the incidence rate of diarrheal diseases in this municipality. For the correlation between river water level and diarrhea diseases, it was chosen three parameters combination: monthly water level mean versus diarrhea incidence, monthly maximum water level versus diarrhea incidence and monthly minimum water level versus diarrhea incidence. It was maintained all cities independent of its location in the basin (whether if it was up or downstream from where the water level was measured).

Results

Testing TRMM dataset

The results obtained through the analysis of contingency table have PC = 73% and CSI = 67%, and the higher the result the more accurate are the estimates. The Critical Success Index value means that the total hits for the rainfall estimates are higher than all the times it missed or generated a false alarm (estimated a rainfall event that was not captured by the station). For bias it was obtained 1.07 which indicates that there is a balance and the satellite estimates are not sub/overestimating the number of events. It also accounts for estimation accuracy. To evaluate reliability and resolution it was calculated FAR, 22%, and H, 83%. Last index, POFD resulted in 49%.

Figure 1, graphic b, presents the rainfall data observed and estimated divided by region, Upper and Lower Acre regions. Median values for both data sources and for both regions vary from approx. 200 to 260 mm monthly. Maximum values floats around 400 mm (estimated) and 500 mm (observed). The minimum values present differences lower than 100 mm comparing the two sources / two regions and they float around 100 mm for the upper range limit.

In Figure 1 c the monthly estimated values basically follow the observed ones which implies in a good representation of the seasonal behavior. Correlation values obtained between rain gauges’ series (observed) and TRMM values were around 0.4.

Figure 1 d introduces the water level exploratory analysis showing the observed daily average during the year, which was built using the time series available. Characteristics as the location of the city relative to the main rivers and the tributaries were taken into consideration to analyze these results.

When a time-lag is applied it is observed a spatial coherence between the rainfall estimates and the water level at Brasiléia gauging station. The higher correlation values obtained was for 2 and 3 days lag, as shown in Figure 2, maps A and B (R ≈ 0.35). These values are the pixels in light grey on the left side of the white dot that indicates the station. The three pixels in the left bottom correspond to Assis Brasil municipality that is less than 100 km from Brasiléia. This confirms that the water transit time is higher the further the rainfall is upstream.

Figure 2
Best correlation values between rainfall estimates and water level taken in Brasiléia, Xapuri and Rio Branco. Check the white dots to identify where the river measurements where obtained.

Xapuri gauging station presents the strongest signals during the 3 and 4 day lag as shown in Figure 2, maps C and D (R ≈ 0.35) and it represents the critical time-lags for the water level in Xapuri to respond to rainfall events in the upstream basin.

Rio Branco gauging station shows the strongest signals for 4 and 5 day lag, as shown in Figure 2, maps E and F (R ≈ 0.31). The pixels located upstream the measurement point have stronger correlation coefficients as the time-lag increases. Short time-lags have low correlation values indicating that rainfall relates with water level in this gauging station more frequently after a minimum time-lag of 3 days.

From Figure 2 e, a pixel in the northwest part of the basin highlights an area of higher influence on the water levels (light grey). This particular pixel corresponds to Espalha River, a small tributary of the Riozinho do Rola River. This pixel maintains higher values until the final time-lag is applied even though its neighbor’s pixels presents lower values as the lag increases. It is not clear why this pixel maintains a higher correlation value despite of the decrease presented by the neighbor pixels. However, the fact that it is next to northern limit of the basin an area that probably has higher altitude values would explain why this area has a larger contribution to the water level variation.

Correlation: rainfall estimates versus diarrheal diseases cases

Box plots for the incidence rate of infant diarrheal diseases cases are shown in Figure 1 e for the eleven municipalities. As the ratio numbers are low (incidence rate) for most of the cities it is difficult to analyze the data. The median values for these cities are around 0.5 (less than one case per 1,000 children). Maximum extreme values are for Assis Brasil (> 4) and Acrelândia (> 3).

Figure 3 b and k show overlapping between the rainfall pixel and the cities, Epitaciolândia and Plácido de Castro, suggesting that the rainfall in that period had an influence in the diarrhea cases in these cities. However, the correlation is not as higher than Senador Guiomard and Acrelândia cities where the pixels with the highest values are located upstream the city.

Figure 3
Correlation results obtained from rainfall estimates for all the Acre river basin and the infant diarrhea diseases cases for all municipalities in the Acre basin (A-G) for the Dec-Jan-Feb.

Correlation: water level data versus diarrheal diseases cases

The results presented in this section comprise the correlation values obtained from the water level data and the diarrheal diseases cases (Figure 4 a-c).

Figure 4
Correlation between Acre River water level for Brasiléia (a), Xapuri (b) and Rio Branco (c) and the infant diarrhea in children in cities that comprises Upper and Lower Acre regions. Correlation values where obtained using maximum, medium and minimum monthly water level for the DJF trimester of each dataset.

Diarrheal diseases recorded in the neighboring cities, Epitaciolândia and Brasiléia, show good correlation values with the Brasiléia gauging station (0.7 and 0.5 respectively, Figure 4 a). However, Xapuri and Rio Branco gauging station show low correlation values with the diseases case registered in the respectively cities.

Discussion

Rainfall estimates validation

Rainfall estimates obtained through remote sensing methods have been used due to its enormous advances in areas with difficult access and logistic (such as Amazonia).

The rainfall estimates validation performed for TRMM using the contingency table found a bias next to a unit, and it is less than what was found in Santos e Silva et al.2828. Santos e Silva CM, Lúcio PS, Spyrides MHC. Distribuição espacial da precipitação sobre o Rio Grande do Norte: estimativa via satélites e medidas por pluviômetros. Rev Bras. Meteorol. 2012; 27(3):337-346..

For the FAR index, the value obtained is reasonable when compared to the 25% approx. obtained in Santos e Silva et al.2828. Santos e Silva CM, Lúcio PS, Spyrides MHC. Distribuição espacial da precipitação sobre o Rio Grande do Norte: estimativa via satélites e medidas por pluviômetros. Rev Bras. Meteorol. 2012; 27(3):337-346. but it was higher than 0.15 found in Scheel et al.2929. Scheel MLM, Rohrer M, Huggel C, Santos Villar D, Silvestre E, Huffman GG. Evaluation of TRMM Multi-satellite Precipitation Analysis (TMPA) performance in the Central Andes region and its dependency on spatial and temporal resolution. Hydrol. Earth Syst. Sci. 2011; 15:2649-2663. study which performed a similar validation process using 3B42-V6 in the Peruvian and Bolivian territory during the rainy season (DJF).

The value obtained for H means that the estimation achieves 4/5 rainfall events and it is a very good proxy in qualitative terms being higher than the ~77% found in Scheel et al.2929. Scheel MLM, Rohrer M, Huggel C, Santos Villar D, Silvestre E, Huffman GG. Evaluation of TRMM Multi-satellite Precipitation Analysis (TMPA) performance in the Central Andes region and its dependency on spatial and temporal resolution. Hydrol. Earth Syst. Sci. 2011; 15:2649-2663.. Moffitt et al.3030. Moffitt CB, Hossain F, Adler RF, Yilmaz KK, Pierce HF. Validation of a TRMM-based global Flood Detection System in Bangladesh. Int J Appl Earth Obs Geoinf 2011; 13(2):165-177. analyzing H, which authors named Probability of Detection index (POD), for the previous version of the TRMM product 3B42V6 (3h) found values of approximately 57% for rain gauges in the Ganges, Brahmaputra and Mehgna basins in Bangladesh. Using these values as a threshold it is assumed that the value obtained in this study is a reasonable match. Authors also found that there are more hits during the rainy season (monsoon) then on the dry season (no-monsoon). Due to the data scarcity, it was not possible to test this hypothesis in our study. The value obtained for POFD index means that there is a 50% percentage of a false event actually being a non-occurrence of rainfall. This value is more than twice of what was found in Moffitt et al.3030. Moffitt CB, Hossain F, Adler RF, Yilmaz KK, Pierce HF. Validation of a TRMM-based global Flood Detection System in Bangladesh. Int J Appl Earth Obs Geoinf 2011; 13(2):165-177.. Authors attributed it to an inherited uncertainty of the satellite estimates. Improvements in the newest version of 3B42 (3h) algorithm are analyzed in the Zulkafli et al.1414. Zulkafli Z, Buytaert W, Onof C, Manz B, Tarnavsky E, Lavado W, Guyot JL. A Comparative Performance Analysis of TRMM 3B42 (TMPA) Versions 6 and 7 for Hydrological Applications over Andean-Amazon River Basins. J Hydrometeorol. 2014; 15:581-592. study as the reduction of the negative bias especially during the rainy season and probably due to the inclusion of more rain gauges in the validation process. It was reassuring for our purposes the accuracy measures being above 60%, the fact that it does not overestimate the number of events and a small false alarm ratio (less than 30%).

This specific TRMM product used in our analysis aggregates the best of all TRMM set of sensors and even data from other satellites although it is still a proxy from the real data (indirect measurement) and it has limitations. The convection processes in this region has characteristics as the weak horizontal fluxes and the fact that most of the momentum is vertical and happens during one-hour time-space3131. Fisch G, Vendrame IF, Hanoaka PCDM. Variabilidade espacial da chuva durante o experimento LBA/TRMM 1999 na Amazônia. Acta Amazon 2007; 37(4):583-590. which are not properly accounted by the sensors and it may lead to over/underestimates. However, TRMM estimates provide a good dataset to analyze rainfall events in this area and its main advantages are due to their availability, spatial resolution, and the series length (15 yr. long). The estimated rainfall pattern in the Upper and Lower Acre is similar to the observed one, although the correlation value between both datasets is low (~0.4). Our hypothesis to the low correlation is the fact that the datasets were compared to each other directly, instead of using mean or accumulated values per municipality. No statistical technique or filters were applied in order to smooth possible failures in these series as the ones listed in Collischonn et al.1313. Collischonn B, Allasia D, Collischonn W, Tucci CEM. Desempenho do satélite TRMM na estimativa de precipitação sobre a bacia do Paraguai Superior. Rev. Bras. Cartograf. 2007; 59(1):93-99..

Although the correlation values are statistically low, we classified them as relevant or high when compared with the results presented in each map or graphic.

Rainfall Estimates versus Water Level

The results for the correlation between rainfall estimates and water level presented a good correspondence despite the low values for the three cities analyzed. The seasonal analysis (not shown) did not improve the results shown here. A hypothesis that might explain this is that longer time intervals would smooth other factors that might be influencing the river seasonal behavior.

There is no correlation between instantaneous values for the rainfall estimates analysis whether it is with water level or diseases. However the results improve when it is applied 1 to 4 days time-lag, for the rainfall and water level analysis. After 5 day time-lag the correlation values tends to be very low.

Comparing the three gauging stations in the first analysis performed, rainfall estimates and water level, the best correlation values was found in Xapuri for 3 day time-lag (R = 0.35).

The results of the Brasiléia gauging station might be influenced by the exclusion of the basin area that does not belong to the Brazilian territory.

The rainfall near the station influences the water level for a few hours when it does so, because there are cases when these amounts of rain join the river downstream or even infiltrates the soil. The effects of flash floods with rapid response on the water level are not registered in the official records, because it usually lasts less than 24 hours.

Diarrheal diseases analysis

It was found an acceptable correlation between the hydrologic variables (rainfall and water level) and diarrheal diseases. It can be due to the short time series available, data gaps and due to the temporal interval in which the analysis were performed (e.g. diarrheal cases are monthly grouped). The result might improve if the diseases statistics were available in a higher frequency (e.g. weekly) which was also suggest by Curriero et al.3232. Curriero FC, Patz JA, Rose JB, Lele S. The association between extreme precipitation outbreaks in the United States, 1948-1994. Am J Public Health Res 2001; 91(8):1194-1199..

This study was not able to show the relationship during the rainy period between rainfall and water level and the diarrheal diseases. However, for Epitaciolândia e Plácido de Castro the correlation values obtained for rainfall and diarrheal disease suggest that the rainfall around the city could influence the disease series. Carlton et al.3333. Carlton EJ, Eisenberg JN, Goldstick J, Cevallos W, Trostle J, Levy K. Heavy rainfall events and diarrhea incidence: the role of social and environmental factors. Am J Epidemiol 2013; 179(3):344-352. also suggested that the highest incidence of diarrheal diseases occurs during the rainy season and that the most affected cities are the ones with low sanitation coverage. Torres et al.3434. Torres RMC, Bittencourt SA, Oliveira RM De, Siqueira ASP, Sabroza PC, Toledo LM. Uso de indicadores de nível local para análise espacial da morbidade por diarreia e sua relação com as condições de vida. Cien Saude Colet 2013; 18(5):1441-1450. assert that the increase in population in urban areas leads to an increase in the pressure and demands for the sanitation system. And although the improvements in the basic Health Attention helped to diminish the number of deaths caused by diarrheal diseases, the lack of adequate sanitation systems maintained the number of cases. According to ACRE3535. Acre. Caracterização Socioambiental das Bacias Hidrográficas do Estado do Acre. Rio Branco: Ministério Público do Acre; 2009., Rio Branco as the state largest city provides a median infrastructure and health services in the state, and concentrates most of the state’s population. A study performed by Oliveira et al.3636. Oliveira AF, Leite IDC, Valente JG. Global burden of diarrheal disease attributable to the water supply and sanitation system in the State of Minas Gerais, Brazil: 2005. Cien Saude Colet 2015; 20(4):1027-1036. for Minas Gerais State, in the Southeast Brazilian region, showed that the increase in the coverage of the sanitation system combined with the actual coverage of the water system could reduce in almost 5% the fraction of diarrheal disease in children less than 5 years attributed to these two infrastructure characteristics.

There were no coherent values between rainfall and Brasiléia diarrheal diseases cases however high correlation values were found between the water level data in Brasiléia and diarrheal diseases cases in both Epitaciolândia (0.7) and Brasiléia (0.5). It can be due to the fact that the analysis with rainfall does not include the Bolivian part of the basin and on the other hand the water level data indirect includes the runoff provided by the excluded area or due to a different response time between the rainfall event and a possible impact on infant’s health. Another possibility to explain these findings may be the distance from the place of residence and health services. However, we could not get data in time for analysis.

Conclusion

This study showed that the rainfall satellite estimates is a good indicator of the rainfall distribution in the Acre basin. It allows the use of this data series to support local public policies and actions of various groups such as civil defense, health professionals to better plan and prioritize the financial and human resources available, and mitigate the impact of floods to the local population. Moreover, these estimate also allowed an analysis of the influence of rainfall upstream in the water level in three gauging stations in the Acre basin.

It is a first approach to aggregate satellite information on water level for children ´s diarrheal diseases analyses in the Amazon basin. It is necessary to deepen statistical analysis, which could allow a better understanding of these processes. Knowing the response time between a rainfall event and its impact on the river and consequences to the population’s health are important issues when it comes to planning response and early warning in case of floods or even extreme rainfall events. Recent studies using climate projections indicate an increase in extreme rainfall events in Amazon3737. Valverde MC, Marengo JA. Extreme Rainfall Indices in the Hydrographic Basins of Brazil. Open Journal of Modern Hydrology 2014; 4(1):10-26.,3838. Fonseca PAM, Veiga JP, Correia FS, Brito A, Queiroz M, Lyra A, Chan SC. Projecting Extreme Changes in Summer Rainfall in South America by the Middle of the 21st Century. Atmos Clim Sci 2014; 4:743-756. and on the discharge of the Amazonian rivers3939. Guimberteau M, Ronchail J, Espinoza JC, Lengaigne M, Sultan B, Polcher J, Drapeau G, Guyot JL, Ducharne A, Ciais P. Future changes in precipitation and impacts on extreme streamflow over Amazonian sub-basins. Environ Res Lett 2013; 8(1):014035.. This type of study can contribute to the projections of health impacts associated with the flooding of the rivers of Amazonia, and consequently help government entities in the strategies of early actions to the local people. Such analysis can improve the preparedness and help the government to adapt its strategies to help the population. The location, strength and time-lag of the spatial correlation calculated here indicate that rainfall is not the only environmental variable explaining the behavior of water level or the occurrence of infant diarrheal diseases cases.

Finally, considering the possible relationship between diarrheal diseases and rainfall behaviour and the fact that other states in the Brazil North region have worse diarrheal diseases ratios there is a scope to replicate this methodology in other contexts within the Brazilian Amazon.

Aknowledgments

The authors thank to the financial support from the National Council for Scientific and Technological Development (CNPq), FINEP project Rede CLIMA 2; Also to the Acre Government and PULSE-Brasil project NERC. We also would like to thank ENSP/FIOCRUZ (National School of Public Health / Oswaldo Cruz Foundation). Special thanks to Dr. Alejandro Duarte Fonseca (AcreBioClima Project) for maintaining such an important project for many years and thus providing important contribution for the society. To Dr. Elisa Armijos from PPG-CLIAMB/UEA-INPA for the revision support. To Diego Viana from the Acre State Health Secretary for the technical support.

References

  • 1
    World Health Organization (WHO). World Health Statistics 2013. [accessed 2014 Dec 1]. Available at: http://www.who.int/gho/publications/world_health_statistics/2013/en/
    » http://www.who.int/gho/publications/world_health_statistics/2013/en/
  • 2
    Floret N, Viel JF, Mauny F, Hoen B, Piarroux R. Negligible risk for epidemics after geophysical disasters. Emerg Infect Dis 2006; 12(4):543.
  • 3
    Ministério da Saúde, Secretaria de Vigilância em Saúde MS/SVS/SIM. Sistema de Informações de Mortalidade (SIM). [accessed 2013 Dec 1]. Available at: http://tabnet.datasus.gov.br/cgi/deftohtm.exe?idb2012/c06.def
    » http://tabnet.datasus.gov.br/cgi/deftohtm.exe?idb2012/c06.def
  • 4
    Epstein PR. Climate Change and emerging infectious diseases. Microbes Infect 2001; 3(9):747-754.
  • 5
    Kovats RS, Bouma MJ, Hajat S, Worral E, Haines A. El Niño and health. Lancet 2003; 362(9394):1481-1489.
  • 6
    Alderman K, Turner LR, Tong S. Floods and human health: A systematic review. Environ Int 2012; 47:37-47.
  • 7
    Kozama S, Aizawa T, Watanabe T, Ranjan P, Gunawardhana L, Amano A. A quantitative risk assessment of waterborne infectious disease in the inundation area of a tropical monsoon region. Sustainable Sciences 2012; 7(1):45-54.
  • 8
    Espinoza JC, Ronchail J, Frappart F, Lavado W, Santini W, Guyot JL. The Major Floods in the Amazonas River and Tributaries (Western Amazon Basin) during the 1970-2012 Period: A Focus on the 2012 Flood. J Hydrometeorol 2013; 14:1000-1008.
  • 9
    Espinoza JC, Marengo JA, Ronchail J, Carpio JM, Flores LN, Guoyt JL. The extreme 2014 flood in the South-western Amazon basin: the role of tropical-subtropical South Atlantic SST gradient. Environ Res Lett 2014; 9:124007.
  • 10
    Portal de Notícias da Globo (G1). União reconhece calamidade em Rio Branco; 300 militares dão apoio. [accessed 2015 Mar 1]. Available at: http://g1.globo.com/ac/acre/noticia/2015/03/uniao-reconhece-calamidade-em-rio-branco-300-militares-serao-enviados.html
    » http://g1.globo.com/ac/acre/noticia/2015/03/uniao-reconhece-calamidade-em-rio-branco-300-militares-serao-enviados.html
  • 11
    Wallace JM, Hobbs PV. Cloud Microphysics. In: Wallace JM, Hobbs PV. Atmospheric science: an introductory survey. Burlington: Academic press; 2006. p. 209-263. (Vol. 92).
  • 12
    Santos e Silva CM, Lúcio OS, Spyrides MHC. Distribuição espacial da precipitação sobre o Rio Grande do Norte: estimativa via satélites e medidas por pluviômetros. Rev. Bras. Meteorol 2012; 27(3):337-346.
  • 13
    Collischonn B, Allasia D, Collischonn W, Tucci CEM. Desempenho do satélite TRMM na estimativa de precipitação sobre a bacia do Paraguai Superior. Rev. Bras. Cartograf 2007; 59(1):93-99.
  • 14
    Zulkafli Z, Buytaert W, Onof C, Manz B, Tarnavsky E, Lavado W, Guyot JL. A Comparative Performance Analysis of TRMM 3B42 (TMPA) Versions 6 and 7 for Hydrological Applications over Andean-Amazon River Basins. J Hydrometeorol 2014; 15:581-592.
  • 15
    Kummerow C, Simpson J, Thiele O, Barnes W, Chang ATC, Stocker E, Adler RF, Hou A, Kakar R, Wentz F, Ashcroft P, Kozu T, Hong Y, Okamoto K, Iguchi T, Kuroiwa H, Haddad EIZ, Huffman G, Ferrier B, Olson WS, Zipser E, Smith EA, Wilheit TT, North G, Krishnamurti T, Nakamura K. The status of the Tropical Rainfall Measuring Mission (TRMM) after two years in orbit. J Appl. Meteorol 2000; 39(12):1965-1982.
  • 16
    Marengo JA, Liebmann B, Grimm AM, Misra V, Silva Dias PL, Cavalcanti IFA, Carvalho LMV, Berbery EH, Ambrizzi T, Vera CS, Saulo AC, Nogues-Paegle J, Zipser E, Seth A, Alves LM. Recent developments on the South American monsoon system. Int J Climatol 2012; 32(1):1-21.
  • 17
    Duarte AF. Aspectos da Climatologia do Acre, Brasil, com base no intervalo 1971-2000. Rev. Bras. Meteorol 2006; 21(3b):308-317.
  • 18
    Ministério Público do Acre (MPE/AC). Caracterização Socioambiental das Bacias Hidrográficas do Estado do Acre. Rio Branco: MPE/AC; 2009.
  • 19
    Instituto Brasileiro de Geografia e Estatística (IBGE). [cited 2013 Dec 1]. Available at: http://www.ibge.gov.br/home/
    » http://www.ibge.gov.br/home/
  • 20
    Agência Nacional de Águas (ANA). HidroWeb: Sistema de Informações Hidrológicas. [cited 2013 Dec 1]. Available at: http://hidroweb.ana.gov.br/
    » http://hidroweb.ana.gov.br/
  • 21
    Grupo de Estudos e Serviços Ambientais (AcreBioClima). [cited 2013 Dec 1]. Available at: http://acrebioclima.net/
    » http://acrebioclima.net/
  • 22
    Tropical Rainfall Measuring Mission Database (TRMM). [cited 2013 Dec 1]. Available at: ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/Derived_Products/3B42_V7/Daily/
    » ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/Derived_Products/3B42_V7/Daily/
  • 23
    Huffman GJ, Adler R, Bolvin D, Nelkin E. The TRMM Multi-Satellite Precipitation Analysis (TMPA). In: Gebremichael M, Hossain F, editors. Satellite Rainfall Applications for Surface Hydrology New York: Springer; 2010. p. 3-22.
  • 24
    Impact of Climate Extremes on Ecosystem and Human Health in Brazil (PULSE-Brazil). [cited 2013 dec 1]. Available at: http://www.pulse-brasil.org/tool/
    » http://www.pulse-brasil.org/tool/
  • 25
    World Health Organization. International Statistical Classification of Diseases and Related Health Problems 10th Revision (ICD-10)-2015-WHO. [cited 2015 Aug 1]. Available at: http://apps.who.int/classifications/icd10/browse/2015/en#/I
    » http://apps.who.int/classifications/icd10/browse/2015/en#/I
  • 26
    Wilks DS. Nonprobabilistic Forecasts of Discrete Predictands. In: Wilks DS. Statistical Methods in the Atmospheric Sciences 2nd ed. Burlington: Academic Press; 2006. p. 260-275.
  • 27
    Warner TT. Verification Model. In: Warner TT. Numerical weather and climate prediction. Cambridge: Cambridge University Press; 2011. p. 294-320.
  • 28
    Santos e Silva CM, Lúcio PS, Spyrides MHC. Distribuição espacial da precipitação sobre o Rio Grande do Norte: estimativa via satélites e medidas por pluviômetros. Rev Bras. Meteorol 2012; 27(3):337-346.
  • 29
    Scheel MLM, Rohrer M, Huggel C, Santos Villar D, Silvestre E, Huffman GG. Evaluation of TRMM Multi-satellite Precipitation Analysis (TMPA) performance in the Central Andes region and its dependency on spatial and temporal resolution. Hydrol. Earth Syst. Sci 2011; 15:2649-2663.
  • 30
    Moffitt CB, Hossain F, Adler RF, Yilmaz KK, Pierce HF. Validation of a TRMM-based global Flood Detection System in Bangladesh. Int J Appl Earth Obs Geoinf 2011; 13(2):165-177.
  • 31
    Fisch G, Vendrame IF, Hanoaka PCDM. Variabilidade espacial da chuva durante o experimento LBA/TRMM 1999 na Amazônia. Acta Amazon 2007; 37(4):583-590.
  • 32
    Curriero FC, Patz JA, Rose JB, Lele S. The association between extreme precipitation outbreaks in the United States, 1948-1994. Am J Public Health Res 2001; 91(8):1194-1199.
  • 33
    Carlton EJ, Eisenberg JN, Goldstick J, Cevallos W, Trostle J, Levy K. Heavy rainfall events and diarrhea incidence: the role of social and environmental factors. Am J Epidemiol 2013; 179(3):344-352.
  • 34
    Torres RMC, Bittencourt SA, Oliveira RM De, Siqueira ASP, Sabroza PC, Toledo LM. Uso de indicadores de nível local para análise espacial da morbidade por diarreia e sua relação com as condições de vida. Cien Saude Colet 2013; 18(5):1441-1450.
  • 35
    Acre. Caracterização Socioambiental das Bacias Hidrográficas do Estado do Acre Rio Branco: Ministério Público do Acre; 2009.
  • 36
    Oliveira AF, Leite IDC, Valente JG. Global burden of diarrheal disease attributable to the water supply and sanitation system in the State of Minas Gerais, Brazil: 2005. Cien Saude Colet 2015; 20(4):1027-1036.
  • 37
    Valverde MC, Marengo JA. Extreme Rainfall Indices in the Hydrographic Basins of Brazil. Open Journal of Modern Hydrology 2014; 4(1):10-26.
  • 38
    Fonseca PAM, Veiga JP, Correia FS, Brito A, Queiroz M, Lyra A, Chan SC. Projecting Extreme Changes in Summer Rainfall in South America by the Middle of the 21st Century. Atmos Clim Sci 2014; 4:743-756.
  • 39
    Guimberteau M, Ronchail J, Espinoza JC, Lengaigne M, Sultan B, Polcher J, Drapeau G, Guyot JL, Ducharne A, Ciais P. Future changes in precipitation and impacts on extreme streamflow over Amazonian sub-basins. Environ Res Lett 2013; 8(1):014035.

Publication Dates

  • Publication in this collection
    Mar 2016

History

  • Received
    15 Aug 2015
  • Reviewed
    04 Dec 2015
  • Accepted
    14 Dec 2015
ABRASCO - Associação Brasileira de Saúde Coletiva Rio de Janeiro - RJ - Brazil
E-mail: revscol@fiocruz.br