Geomapping generalized eigenvalue frequency distributions for predicting prolific Aedes albopictus and Culex quinquefasciatus habitats based on spatiotemporal field-sampled count data.

Marked spatiotemporal variabilities in mosquito infection of arboviruses require adaptive strategies for determining optimal field-sampling timeframes, pool screening, and data analyses. In particular, the error distribution and aggregation patterns of adult arboviral mosquitoes can vary significantly by species, which can statistically bias analyses of spatiotemporal-sampled predictor variables generating misinterpretation of prolific habitat surveillance locations. Currently, there is a lack of reliable and consistent measures of risk exposure based on field-sampled georeferenced explanatory covariates which can compromise quantitative predictions generated from arboviral mosquito surveillance models for implementing larval control strategies targeting productive habitats. In this research we used spatial statistics and QuickBird visible and near-infra-red data for determining trapping sites that were related to Culex quinquefasciatus and Aedes albopictus species abundance and distribution in Birmingham, Alabama. Initially, a Land Use Land Cover (LULC) model was constructed from multiple spatiotemporal-sampled georeferenced predictors and the QuickBird data. A Poisson regression model with a non-homogenous, gamma-distributed mean then decomposed the data into positive and negative spatial filter eigenvectors. An autoregressive process in the error term then was used to derive the sample distribution of the Moran's I statistic for determining latent autocorrelation components in the model. Spatial filter algorithms established means, variances, distributional functions, and pairwise correlations for the predictor variables. In doing so, the eigenfunction spatial filter quantified the residual autocorrelation error in the mean response term of the model as a linear combination of various distinct Cx. quinquefasciatus and Ae. albopictus habitat map patterns. The analyses revealed 18-27% redundant information in the data. Prolific habitats of Cx. quinquefasciatus and Ae. albopictus can be accurately spatially targeted based on georeferenced field-sampled count data using QuickBird data, LULC explanatory covariates, robust negative binomial regression estimates and space-time eigenfunctions.

A heteroskedastic error covariance matrix estimator using a first-order conditional autoregressive Markov simulation for deriving asympotical efficient estimates from ecological sampled Anopheles arabiensis aquatic habitat covariates.

BACKGROUND: Autoregressive regression coefficients for Anopheles arabiensis aquatic habitat models are usually assessed using global error techniques and are reported as error covariance matrices. A global statistic, however, will summarize error estimates from multiple habitat locations. This makes it difficult to identify where there are clusters of An. arabiensis aquatic habitats of acceptable prediction. It is therefore useful to conduct some form of spatial error analysis to detect clusters of An. arabiensis aquatic habitats based on uncertainty residuals from individual sampled habitats. In this research, a method of error estimation for spatial simulation models was demonstrated using autocorrelation indices and eigenfunction spatial filters to distinguish among the effects of parameter uncertainty on a stochastic simulation of ecological sampled Anopheles aquatic habitat covariates. A test for diagnostic checking error residuals in an An. arabiensis aquatic habitat model may enable intervention efforts targeting productive habitats clusters, based on larval/pupal productivity, by using the asymptotic distribution of parameter estimates from a residual autocovariance matrix. The models considered in this research extends a normal regression analysis previously considered in the literature. METHODS: Field and remote-sampled data were collected during July 2006 to December 2007 in Karima rice-village complex in Mwea, Kenya. SAS 9.1.4 was used to explore univariate statistics, correlations, distributions, and to generate global autocorrelation statistics from the ecological sampled datasets. A local autocorrelation index was also generated using spatial covariance parameters (i.e., Moran's Indices) in a SAS/GIS database. The Moran's statistic was decomposed into orthogonal and uncorrelated synthetic map pattern components using a Poisson model with a gamma-distributed mean (i.e. negative binomial regression). The eigenfunction values from the spatial configuration matrices were then used to define expectations for prior distributions using a Markov chain Monte Carlo (MCMC) algorithm. A set of posterior means were defined in WinBUGS 1.4.3. After the model had converged, samples from the conditional distributions were used to summarize the posterior distribution of the parameters. Thereafter, a spatial residual trend analyses was used to evaluate variance uncertainty propagation in the model using an autocovariance error matrix. RESULTS: By specifying coefficient estimates in a Bayesian framework, the covariate number of tillers was found to be a significant predictor, positively associated with An. arabiensis aquatic habitats. The spatial filter models accounted for approximately 19% redundant locational information in the ecological sampled An. arabiensis aquatic habitat data. In the residual error estimation model there was significant positive autocorrelation (i.e., clustering of habitats in geographic space) based on log-transformed larval/pupal data and the sampled covariate depth of habitat. CONCLUSION: An autocorrelation error covariance matrix and a spatial filter analyses can prioritize mosquito control strategies by providing a computationally attractive and feasible description of variance uncertainty estimates for correctly identifying clusters of prolific An. arabiensis aquatic habitats based on larval/pupal productivity.

Describing Anopheles arabiensis aquatic habitats in two riceland agro-ecosystems in Mwea, Kenya using a negative binomial regression model with a non-homogenous mean.

This research illustrates a geostatistical approach for modeling the spatial distribution patterns of Anopheles arabiensis Patton (Patton) aquatic habitats in two riceland environments. QuickBird 0.61 m data, encompassing the visible bands and the near-infra-red (NIR) band, were selected to synthesize images of An. arabiensis aquatic habitats. These bands and field sampled data were used to determine ecological parameters associated with riceland larval habitat development. SAS was used to calculate univariate statistics, correlations and Poisson regression models. Global autocorrelation statistics were generated in ArcGISfrom georeferenced Anopheles aquatic habitats in the study sites. The geographic distribution of Anopheles gambiae s.l. aquatic habitats in the study sites exhibited weak positive autocorrelation; similar numbers of log-larval count habitats tend to clustered in space. Individual rice land habitat data were further evaluated in terms of their covariations with spatial autocorrelation, by regressing them on candidate spatial filter eigenvectors. Each eigenvector generated from a geographically weighted matrix, for both study sites, revealed a distinctive spatial pattern. The spatial autocorrelation components suggest the presence of roughly 14-30% redundant information in the aquatic habitat larval count samples. Synthetic map pattern variables furnish a method of capturing spatial dependency effects in the mean response term in regression analyses of rice land An. arabiensis aquatic habitat data.

Hydrological modeling of geophysical parameters of arboviral and protozoan disease vectors in Internally Displaced People camps in Gulu, Uganda.

BACKGROUND: The aim of this study was to determine if remotely sensed data and Digital Elevation Model (DEM) can test relationships between Culex quinquefasciatus and Anopheles gambiae s.l. larval habitats and environmental parameters within Internally Displaced People (IDP) campgrounds in Gulu, Uganda. A total of 65 georeferenced aquatic habitats in various IDP camps were studied to compare the larval abundance of Cx. quinquefasciatus and An. gambiae s.l. The aquatic habitat dataset were overlaid onto Land Use Land Cover (LULC) maps retrieved from Landsat imagery with 150 m x 150 m grid cells stratified by levels of drainage. The LULC change was estimated over a period of 14 years. Poisson regression analyses and Moran's I statistics were used to model relationships between larval abundance and environmental predictors. Individual larval habitat data were further evaluated in terms of their covariations with spatial autocorrelation by regressing them on candidate spatial filter eigenvectors. Multispectral QuickBird imagery classification and DEM-based GIS methods were generated to evaluate stream flow direction and accumulation for identification of immature Cx. quinquefasciatus and An. gambiae s.l. and abundance. RESULTS: The main LULC change in urban Gulu IDP camps was non-urban to urban, which included about 71.5 % of the land cover. The regression models indicate that counts of An. gambiae s.l. larvae were associated with shade while Cx. quinquefasciatus were associated with floating vegetation. Moran's I and the General G statistics for mosquito density by species and instars, identified significant clusters of high densities of Anopheles; larvae, however, Culex are not consistently clustered. A stepwise negative binomial regression decomposed the immature An. gambiae s.l. data into empirical orthogonal bases. The data suggest the presence of roughly 11% to 28 % redundant information in the larval count samples. The DEM suggest a positive correlation for Culex (0.24) while for Anopheles there was a negative correlation (-0.23) for a local model distance to stream. CONCLUSION: These data demonstrate that optical remote sensing; geostatistics and DEMs can be used to identify parameters associated with Culex and Anopheles aquatic habitats.

Remote and field level quantification of vegetation covariates for malaria mapping in three rice agro-village complexes in Central Kenya.

BACKGROUND: We examined algorithms for malaria mapping using the impact of reflectance calibration uncertainties on the accuracies of three vegetation indices (VI)'s derived from QuickBird data in three rice agro-village complexes Mwea, Kenya. We also generated inferential statistics from field sampled vegetation covariates for identifying riceland Anopheles arabiensis during the crop season. All aquatic habitats in the study sites were stratified based on levels of rice stages; flooded, land preparation, post-transplanting, tillering, flowering/maturation and post-harvest/fallow. A set of uncertainty propagation equations were designed to model the propagation of calibration uncertainties using the red channel (band 3: 0.63 to 0.69 microm) and the near infra-red (NIR) channel (band 4: 0.76 to 0.90 microm) to generate the Normalized Difference Vegetation Index (NDVI) and the Soil Adjusted Vegetation Index (SAVI). The Atmospheric Resistant Vegetation Index (ARVI) was also evaluated incorporating the QuickBird blue band (Band 1: 0.45 to 0.52 microm) to normalize atmospheric effects. In order to determine local clustering of riceland habitats Gi*(d) statistics were generated from the ground-based and remotely-sensed ecological databases. Additionally, all riceland habitats were visually examined using the spectral reflectance of vegetation land cover for identification of highly productive riceland Anopheles oviposition sites. RESULTS: The resultant VI uncertainties did not vary from surface reflectance or atmospheric conditions. Logistic regression analyses of all field sampled covariates revealed emergent vegetation was negatively associated with mosquito larvae at the three study sites. In addition, floating vegetation (-ve) was significantly associated with immature mosquitoes in Rurumi and Kiuria (-ve); while, turbidity was also important in Kiuria. All spatial models exhibit positive autocorrelation; similar numbers of log-counts tend to cluster in geographic space. The spectral reflectance from riceland habitats, examined using the remote and field stratification, revealed post-transplanting and tillering rice stages were most frequently associated with high larval abundance and distribution. CONCLUSION: NDVI, SAVI and ARVI generated from QuickBird data and field sampled vegetation covariates modeled cannot identify highly productive riceland An. arabiensis aquatic habitats. However, combining spectral reflectance of riceland habitats from QuickBird and field sampled data can develop and implement an Integrated Vector Management (IVM) program based on larval productivity.