|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Recently, we found an association between antecedent drought, coincident wetting, and transmission of SLE virus in Indian River County Florida.2 We used a dynamic hydrology model3 to hindcast mean area water table depth (WTD) in Indian River County, and compared this model simulation to sentinel chicken seroconversion data. Seroconversion of sentinel chickens, as measured by serum assay for hemagglutination inhibition (HI) antibodies to SLE virus, has been strongly correlated with the clinical disease in humans.4 Using logistic regression, we found the probability of sentinel chicken seroconversion, i.e., transmission of SLE virus, to be strongly associated with low WTDs 1117 weeks prior and higher WTDs 02 weeks prior.2
A mechanism for this empirical relationship was suggested by mosquito collection data obtained in Indian River County from 1986 to 1991. Culex nigripalpus Theobald is the demonstrated enzootic and epidemic vector of SLE virus in south Florida.57 Collections of Cx. nigripalpus were made in densely vegetated "hammock" (an island of dense native vegetation) habitats used by this species for daytime resting.8 During the driest conditions (low modeled WTD) preceding heavy SLE virus transmission, Cx. nigripalpus collections dramatically increased.2 Rather than indicate an increase of mosquito abundance, these data suggest that drought restricts Cx. nigripalpus activity to the more humid hammock habitats. Extreme drought periods in south Florida tend to occur during the spring, a time when nesting wild birds also make use of the hammocks. Thus, drought drives the mosquitoes and birds into contact with one another. This forced interaction of vector mosquitoes and susceptible avian amplification hosts provides an ideal environment for the rapid epizootic amplification of SLE virus. In addition, confinement of blood fed and gravid Cx. nigripalpus females to the hammock habitats for extended periods forces more infected females to complete the extrinsic incubation of acquired arboviruses during a single gonotrophic cycle.1,8 Subsequently, when the drought ends and water resources increase, infected mosquitoes and birds disperse and carry the virus from the hammocks. Gravid female mosquitoes oviposit, re-feed, and if infective, they initiate the early transmission phase of the Florida SLE cycle.
Additional supporting evidence was found from wild bird serologic data obtained during the 1990 SLE epidemic. We showed that the presence of HI antibodies to SLE virus in wild birds was significantly associated with modeled antecedent drought, coincident wetting, and the emergence of the Florida SLE virus vector, Cx. nigripalpus, five weeks prior.9 Our findings suggest that three factors conspired to create the 1990 epidemic: 1) a large population of susceptible wild birds; 2) severe springtime drought, which facilitated amplification of the SLE virus among the Cx. nigripalpus and a portion of the wild bird population; 3) continued rainfall and wetting of the land surface in the summer and early fall, which sustained a large, active Cx. nigripalpus population. The continued biting and reproductive activity of Cx. nigripalpus maintained levels of epizootic transmission throughout the summer and early fall in Indian River County. These constant high levels of SLE virus amplification facilitated spillover transmission to humans throughout much of southcentral Florida.
In this study, we shift our analysis to human cases of SLE reported throughout southcentral Florida. We hypothesize that the spatial-temporal organization of human SLE cases in Florida is tied to hydrologic variability. Specifically, we propose that the epidemic transmission with Indian River County as the epicenter in 1990, the focal transmission with Lee and Collier Counties as the focus in 1993, and the sporadic transmission in 1997, correspond to the spatial-temporal patterns of drought and wetting during these years. We expect that differences in the number and location of human cases of SLE among years reflect differences in land surface wetness conditions across southcentral Florida. We further hypothesize that foci of drought and wetting should be demonstrable within or near Indian River County during 1990, at the Lee/Collier County border in 1993, and scattered among several counties during 1997. To support our model of SLE virus amplification and transmission, these hydrologic epicenters should be manifest as spring droughts followed by wetting concurrent or prior to that of surrounding areas.
Finally, we hypothesize that antecedent drought is a necessary, although not sufficient, condition for SLE virus amplification and transmission in Florida. Thus, all counties recording human SLE cases should include or be near localities that experienced drought prior to the onset of human infections. To test these hypotheses we acquired 19881998 meteorologic data from available meteorologic stations in all Florida counties lying between 26.5°N and 29°N. These data were used to force a dynamic hydrology model and retrospectively simulate surface wetness conditions in the area of each station.
| METHODS AND MATERIALS |
|
|
|---|
The TBH model combines a soil column, which simulates the vertical movement of water and heat within the soil and between the soil surface, vegetation and the atmosphere, with the TOPMODEL approach,10 which incorporates the statistics of topography to track the horizontal movement of shallow groundwater from the uplands to the lowlands. The TOPMODEL formulations permit calculation of both the saturated fraction within the watershed (partial contributing area), and the groundwater flow that supports this area, from knowledge of the mean WTD and a probability density function for soil moisture deficit derived from topographic statistics. This approach to modeling the land surface has been validated at several catchments, ranging in scale from the Red Arkansas Basin (570,000 km2)11 to the Black Rock Forest catchment (1.34 km2).12
Hourly meteorologic data. Seven hourly meteorologic variables, precipitation, temperature, surface pressure, wind speed, mixing ratio, incident solar radiation, and downwelling longwave radiation, are needed to force the TBH model. Because hourly data is of limited availability we use daily precipitation and temperature data and a resampling procedure to generate the hourly forcing data sets. This resampling procedure is performed to accommodate the TBH model interface; in so doing, daily variability is maintained in its entirety. See Shaman and others13 for a more complete explanation of this procedure.
The hourly meteorologic data used for the resampling procedure were assembled from National Climate Data Center (NCDC) archives for Vero Beach, Florida. Gaps in the hourly Vero Beach record were filled with hourly data from NCDC archives for Melbourne and West Palm Beach, Florida. In addition, hourly meteorologic data were assembled from the NCDC Solar and Meteorological Surface Observation Network (SAMSON) dataset for Daytona Beach, Miami, Tampa, and West Palm Beach. For the Vero Beach hourly record, solar radiation data were provided by the Northeast Regional Climate Center (NRCC) from analysis of the NCDC data using the NRCC solar energy model.14 Solar radiation data are included in the SAMSON dataset.
Daily meteorological data.
Daily precipitation and temperature data were assembled from NCDC archives for all stations with near continuous records spanning 19881998 for the 28 Florida counties lying partly or wholly between 26.5°N and 29°N. Small gaps in the daily records were filled with data from adjacent stations. A total of 42 records of 19881998 daily data were assembled. Table 1
provides a listing of the daily record station sites. Using the closest hourly meteorologic record (Daytona Beach, Miami, Tampa, Vero Beach, or West Palm Beach) and the resampling procedure,13 we then created hourly forcing data sets for each of the 42 station sites.
|
At the Vero Beach 4W site, the TBH model was run from 1949 through 2001, and model output provided a daily time series of mean WTD and river runoff for this site area. Model validation was performed using groundwater well measurements and surface (canal) water levels, provided by the St. Johns Water Management District and USGS sources. Partitioning of runoff and evapotranspiration matched bulk estimates derived from USGS sources.
The TBH model simulations. Model simulations at all 42 station sites were performed using the calibrations established at the Vero Beach 4W site. Most of southcentral Florida has similar flat topography, similar tropical vegetation, and is subject to channelization and water control. These commonalities permit this regional validation and application of the TBH model to these other areas, many of which lack local validation data.
Mean area WTD, as simulated by the TBH model, is not representative of a measurable quantity (it is the mean for the area) but rather serves as an index of relative wetness for the area through time. Use of the single, regional model calibration permits intercomparison of mean area WTD among the site simulations.
Human case data.
Summaries of mosquito-borne arboviral human case data are compiled, analyzed, and reported weekly, monthly, and annually by the Florida Department of Health (Tallahassee, FL). Table 2
shows the monthly human cases of SLE as reported by county for 19901998.
|
Bivariate logistic regression analysis was used to associate the probability of countywide human SLE incidence or prevalence with lag combinations of half-monthly modeled WTD at each station site within the county. All counties were analyzed in aggregate. Whole model goodness-of-fit was measured by log-likelihood ratio and the pseudo r2 (uncertainty) coefficient. Individual parameter estimates were made using a maximum likelihood procedure. Wald chi-square tests were used to determine whether these estimates were significantly different from zero.
Additional analyses were performed to determine whether drought or drought and wetting were necessarily associated with human SLE cases. Between 1990 and 1998 and for the counties represented in this study, there were 71 instances in which a county reported one or more monthly cases of human SLE, i.e., human SLE incidence. We hypothesized that antecedent drought should have occurred within or near all of these counties. We therefore calculated the number of these 71 occurrences for which antecedent drought fell below a given mean monthly WTD (e.g., below 1.35 meters, representative of deep drought) somewhere within the county. The antecedence lag length used for this analysis was the best fit of the logistic regression analysis. Many counties have more than one station record, and for this analysis we only required that at least one station meet the antecedent drought criterion. A similar calculation was also performed for both antecedent drought and near-coincident wetting (e.g., a WTD below 1.4 meters, and a near coincident wetting above 1.2 meters).
Bootstrap confidence intervals were then estimated using a Monte Carlo procedure. A total of 10,000 combinations of the 71 county-months were randomly sampled among the years 19901998 (1988 and 1989 were excluded to allow for TBH model spin-up), and a distribution of the totals meeting each defined criterion of drought (or drought and wetness) was constructed. For example, because Indian River County recorded human SLE cases in July 1990, a random selection of antecedent WTDs among the years 19901998 (only for Indian River County and July) would be made, and whether the drought criterion was met within the county for that randomly selected year would be determined. This random selection would then be repeated for the remaining 70 county-months of human SLE incidence, the total recorded, and the process repeated 10,000 times. The significance of the actual number of antecedent drought instances (or antecedent drought and near-coincident wetting) occurring with human SLE incidence was then assessed based on this distribution of 10,000. The null hypothesis is that the number of counties meeting the drought (or drought and wetness) criterion is no greater than that due to chance.
| RESULTS |
|
|
|---|
![]() |
where P(SLE) is the probability of human SLE incidence, WTD4 is the half-monthly mean WTD four months prior, and WTD0.5 is the half-monthly mean WTD one-half month prior. This pattern of drought and wetting associated with human SLE incidence is similar to the pattern of drought and wetting found to be associated with SLE virus transmission to sentinel chickens in Indian River County.2 Using the best-fit model time lags found for the above analysis, we repeated our logistic regression analysis using human SLE prevalence. Drought four months prior and wetting one-half month prior were significantly associated with all levels of human SLE prevalence.
Figure 1a
shows a scatter plot of mean WTD four months prior versus mean WTD one-half month prior for all station sites for 19901998. The human SLE incidence associated with each plot value is indicated in color: blue if no human SLE cases occurred within the county, and red if one or more human SLE cases occurred within the county. Stations for which there was human SLE incidence within the county are clustered at the driest conditions. As a comparison, Figure 1b
shows the same plot but with an equal number of station sites chosen randomly and colored red.
|
|
Our previous analyses2 compared human SLE incidence with drought using fixed time lags. However, our working hypothesis, based on the biology of the system, is that drought during spring brings nesting wild birds and Cx. nigripalpus into increased contact facilitating amplification.2 We therefore repeated our analysis, but instead of using drought four months prior to the onset of human SLE cases, we restricted the analysis of drought to the preceding May (Table 4
). Our findings are significant and indicate that the onset of human SLE cases was predicated on drought during the previous May. Specifically, a drought of defined magnitude < 1.35 meters WTD occurred in May prior to the onset of human SLE cases in 59 of 71 possible instances (83%).
|
To illustrate the possible effects of both these factors, we present plots of land surface wetness for the three study years in which human SLE cases were reported in south Florida: the 1990 epidemic; the 1993 focal event; and the 1997 sporadic event. Figure 2
shows the May 1990 mean modeled WTD for the 42 station sites, the June 1630, 1990 mean modeled WTD for the 42 station sites, and the July occurrence of human SLE cases in Florida by county. Indian River County was the only county reporting human SLE cases. Pervasive drought is evident in May throughout southcentral Florida; most of the sites show a mean WTD of 1.3 meters or lower. The Vero Beach station in eastern Indian River County shows a May drought, but with no visible wetting in late June. Western Indian River County has no meteorologic station, yet just over the border at Fort Drum in Okeechobee County, a strong May drought and late June wetting is apparent. This wetting (denoted by the arrow), by virtue of its location, may have facilitated mosquito dispersal into Indian River County. It also suggests that the network of meteorologic stations may be too sparse for comprehensive modeling of land surface wetness (i.e., wetting may have occurred, but was undetected, in western Indian River County).
|
|
|
|
We also mapped the sporadic SLE transmission event in 1997. Figure 6
shows the May 1997 mean modeled WTD for the 42 stations, the August 1631, 1997 mean modeled WTD for the stations, and the September occurrence of human SLE cases by county. Human cases were reported in September by 4 of the 28 counties within the study area. May drought was pervasive in the northern half of the study area, and in patches along the southern extreme. By the second half of August, most stations reported a modeled WTD above 1.2 meters. May drought and wetting one-half month prior to the onset of SLE cases occurred in all four counties where human SLE cases were reported in 1997.
|
| DISCUSSION |
|
|
|---|
We have shown that human SLE cases in southcentral Florida from 1990 to 1998 were significantly associated with drought 36 months prior and wetting 01 months prior to onset. This range of lags reflects the slow variability of modeled WTD, particularly the slow drop of the WTD as the model soil column dries during drought conditions. The best-fit model of human SLE incidence was shown to be drought four months prior and wetting one-half month prior to case onset, and these lags were also significant predictors of the prevalence of human SLE cases. These logistic regression models illustrate that the spatial-temporal distribution of human SLE cases is indeed associated with spatial-temporal patterns of drought and wetting in south-central Florida. These empirical models using TBH model simulations may also enable development of skillful forecasts of human SLE epidemics in southcentral Florida.
Using maps of modeled WTD we were able to identify the spatial-temporal evolution of land surface wetness in south-central Florida from 1990 through 1998. In 1990, there was pervasive drought throughout the region. Within the network of meteorologic stations used for this study, we were able to detect initial surface wetting in late June just west of Indian River County. The first human SLE case in 1990 was reported in eastern Indian River County. Wetting may also have occurred in the western region of this county where there was no meteorologic record. In 1993, a focus of drought was identified in Lee and Hendry Counties. This drought may have also affected northern Collier County where again, there was no meteorologic station record. In 1997, May drought occurred in many areas of southcentral Florida, and sporadic human SLE cases were reported in counties affected by drought and subsequent wetting. These maps (Figures 2
6
) demonstrate the resolution with which modeled land surface wetness can be used to detect loci of SLE virus amplification and transmission. A denser network of meteorologic stations would provide a higher resolution picture of evolving land surface wetness conditions.
Differences among the years suggest that the spatial extent of May drought delimits where amplification occurs and the type of transmission that will follow (i.e., epidemic, focal, or sporadic). Thus, the spatial structure of the May 1993 drought appears to provide an explanation of why human SLE cases occurred focally in Lee and Collier Counties. Similarly, the May 1990 drought and subsequent wetting throughout south-central Florida separates this year from 1993 and indicates that the 1990 SLE epidemic was broadly hydrologically preconditioned and driven.
In our study, we also show that human SLE cases are less strongly associated with near-coincident wetting than with antecedent (May) drought. This finding is consistent with previous results, in which near-coincident wetting was more weakly associated than antecedent drought with the transmission of SLE virus to sentinel chickens.2 Within our working hypothesis of the Florida epidemic SLE transmission cycle, drought induces amplification of SLE virus in select habitats. In addition, these drought conditions allow more Cx. nigripalpus females to become infective during a single gonotrophic cycle.1 Once land surface wetting occurs these infectious Cx. nigripalpus can range more broadly and thus are more likely to encounter humans. However, rainfall may also permit dispersal of infected Cx. nigripalpus.1,8 While a single pulse of rain may not be sufficient to increase the availability of oviposition sites, it can, and does, raise humidity levels allowing host-seeking females to disperse farther. The daily flight behavior of Cx. nigripalpus is significantly correlated with increased relative humidity one hour after sunset.15 Local rainfall increases humidity and allows Cx. nigripalpus females, infected and uninfected alike, to increase their daily flight time and flight range.16 Thus, the behavioral response of Cx. nigripalpus to local rainfall events might explain the weaker association found between near-coincident land surface wetting and SLE virus transmission.
Our final hypothesis is that drought should always precede the onset of human SLE cases. We found that drought occurred four months prior to the onset of human cases in 73% of the counties reporting SLE cases. By focusing exclusively upon May, when wild birds are nesting and when we hypothesize the bulk of SLE amplification occurs, we showed that drought occurred in 83% of the counties reporting human SLE cases. However, if drought-induced amplification of SLE virus is a necessary precondition for human SLE cases, pre-transmission drought should always be present in counties reporting human SLE cases.
This shortcoming may be a consequence of a mismatch between the spatial variability of subtropical precipitation and the incomplete geographic coverage presently provided by the existing network of meteorologic stations. To first order, hydrologic conditions in subtropical Florida should be localized at the spatial scales at which storms organize. With the exception of irrigation, precipitation is the sole source of water to near surface soils, and it is these wet soils that support the surface pooling that forms the preferred breeding sites of Cx. nigripalpus. Because much of the rainfall in Florida is convective, land surface wetness can vary over short spatial scales. Examination of TBH model-simulated land surface wetness maps for 1990, 1993, and 1997 in southcentral Florida (Figures 2
6
) suggests that not all of the spatial variability of wetness is being captured by the present network of meteorologic stations. Where station sites are close, modeled land surface wetness can vary considerably over short distances, on the order tens of kilometers. However, over much of the study area, the network of meteorologic stations is too sparse to comprehensively detect such smaller scale variability. These gaps in our simulation of the land surface wetness seem to explain why May drought is not found within all counties subsequently reporting human SLE cases. Mosquito, wild bird, and human movement in and out of a county may also have introduced noise to our analysis.
While it has not been clearly established to be a necessary precondition for human SLE incidence, overall, the analysis does demonstrate that May drought is significantly associated with subsequent human SLE incidence. This finding supports our working hypothesis that during drought mosquitoes and nesting wild birds are brought into close contact, facilitating epizootic amplification and generating the mosquito infection rates necessary to support high levels of SLE virus transmission. Further testing of this working hypothesis requires more evidence from the field including sampling of mosquito and wild bird abundance and infection rates within southcentral Florida hammock habitats. Using these data we will be able to document arboviral amplification during periods of drought, and more precisely determine how mosquitoes respond to rainfall and land surface wetting. We will also be able to explore the ability of the TBH model to predict the reproductive response of mosquitoes to short periods of wetting.
Future predictions of human SLE incidence may also need to consider virologic factors, such as increased herd immunity to the SLE virus from heavy exposure in preceding years. The effects of a variety of behaviors on human SLE incidence might also be explored; these could include occupational and avocational habits, such as whether SLE is more prevalent in agricultural workers and outdoorsmen.
Our findings support three hypotheses: 1) the spatial-temporal variability of antecedent drought and near-coincident wetting is associated the spatial-temporal variability of human SLE cases; 2) epicenters of human SLE are identifiable as foci of drought and wetting; and 3) drought is a measurable precondition associated with subsequent human SLE incidence.
Large-scale epidemics are of primary concern to public health workers. The findings of this study suggest that pervasive May drought and summer wetting throughout southcentral Florida enables such an epidemic event. Such large-scale hydrologic events are easily monitored using records from the current network of meteorologic stations to model land surface wetness. A denser network of stations may be necessary to monitor and predict smaller-scale focal and sporadic transmission events in southcentral Florida.
Received July 3, 2003. Accepted for publication January 1, 2004.
Acknowledgment: We thank Carina Blackmore, Caroline Collins and their staff at the Florida Department of Health in Tallahassee for their work tracking and reporting human arbovirus cases in Florida. This report is Florida Agricultural Experiment Station Journal Series R-101091.
Financial support: This work was supported by National Aeronautics and Space Administration Earth System Science Fellowship NGT5-50323 and the National Oceanic and Atmospheric Administration Postdoctoral Program in Climate and Global Change, administered by the University Corporation for Atmospheric Research.
Authors addresses: Jeffrey Shaman, Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138, E-mail: jshaman{at}fas.harvard.edu. Jonathan F. Day, Florida Medical Entomology Laboratory, Institute of Food and Agricultural Sciences, University of Florida, Vero Beach, FL 32962, E-mail: JFDA{at}mail.ifas.ufl.edu. Marc Stieglitz, Lamont-Doherty Earth Observatory, Columbia University, Palisades, NY 10964, E-mail: marc{at}ldeo.columbia.edu.
| REFERENCES |
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |