On spatially explicit models of cholera epidemics

Modelling techniques in cholera epidemiology: A systematic and critical review

Diverse modelling techniques in cholera epidemiology have been developed and used to (1) study its transmission dynamics, (2) predict and manage cholera outbreaks, and (3) assess the impact of various control and mitigation measures. In this study, we carry out a critical and systematic review of various approaches used for modelling the dynamics of cholera. Also, we discuss the strengths and weaknesses of each modelling approach. A systematic search of articles was conducted in Google Scholar, PubMed, Science Direct, and Taylor & Francis. Eligible studies were those concerned with the dynamics of cholera excluding studies focused on models for cholera transmission in animals, socio-economic factors, and genetic & molecular related studies. A total of 476 peer-reviewed articles met the inclusion criteria, with about 40% (32%) of the studies carried out in Asia (Africa). About 52%, 21%, and 9%, of the studies, were based on compartmental (e.g., SIRB), statistical (time series and regression), and spatial (spatiotemporal clustering) models, respectively, while the rest of the analysed studies used other modelling approaches such as network, machine learning and artificial intelligence, Bayesian, and agent-based approaches. Cholera modelling studies that incorporate vector/housefly transmission of the pathogen are scarce and a small portion of researchers (3.99%) considers the estimation of key epidemiological parameters. Vaccination only platform was utilized as a control measure in more than half (58%) of the studies. Research productivity in cholera epidemiological modelling studies have increased in recent years, but authors used diverse range of models. Future models should consider incorporating vector/housefly transmission of the pathogen and on the estimation of key epidemiological parameters for the transmission of cholera dynamics.

A two-phase fluid model for epidemic flow

We propose a new mathematical and computational modeling framework that incorporates fluid dynamics to study the spatial spread of infectious diseases. We model the susceptible and infected populations as two inviscid fluids which interact with each other. Their motion at the macroscopic level characterizes the progression and spread of the epidemic. To implement the two-phase flow model, we employ high-order numerical methods from computational fluid dynamics. We apply this model to simulate the COVID-19 outbreaks in the city of Wuhan in China and the state of Tennessee in the US. Our modeling and simulation framework allows us to conduct a detailed investigation into the complex spatiotemporal dynamics related to the transmission and spread of COVID-19.

Computation of the basic reproduction numbers for reaction-diffusion epidemic models

We consider a class of $ k $-dimensional reaction-diffusion epidemic models ($ k = 1, 2, \cdots $) that are developed from autonomous ODE systems. We present a computational approach for the calculation and analysis of their basic reproduction numbers. Particularly, we apply matrix theory to study the relationship between the basic reproduction numbers of the PDE models and those of their underlying ODE models. We show that the basic reproduction numbers are the same for these PDE models and their associated ODE models in several important scenarios. We additionally provide two numerical examples to verify our analytical results.

Mathematical Models for Cholera Dynamics-A Review

Cholera remains a significant public health burden in many countries and regions of the world, highlighting the need for a deeper understanding of the mechanisms associated with its transmission, spread, and control. Mathematical modeling offers a valuable research tool to investigate cholera dynamics and explore effective intervention strategies. In this article, we provide a review of the current state in the modeling studies of cholera. Starting from an introduction of basic cholera transmission models and their applications, we survey model extensions in several directions that include spatial and temporal heterogeneities, effects of disease control, impacts of human behavior, and multi-scale infection dynamics. We discuss some challenges and opportunities for future modeling efforts on cholera dynamics, and emphasize the importance of collaborations between different modeling groups and different disciplines in advancing this research area.

Spatial effects in parasite-induced marine diseases of immobile hosts

Emerging marine infectious diseases pose a substantial threat to marine ecosystems and the conservation of their biodiversity. Compartmental models of epidemic transmission in marine sessile organisms, available only recently, are based on non-spatial descriptions in which space is homogenized and parasite mobility is not explicitly accounted for. However, in realistic scenarios epidemic transmission is conditioned by the spatial distribution of hosts and the parasites' mobility patterns, calling for an explicit description of space. In this work, we develop a spatially explicit individual-based model to study disease transmission by waterborne parasites in sessile marine populations. We investigate the impact of spatial disease transmission through extensive numerical simulations and theoretical analysis. Specifically, the effects of parasite mobility into the epidemic threshold and the temporal progression of the epidemic are assessed. We show that larger values of pathogen mobility imply more severe epidemics, as the number of infections increases, and shorter timescales to extinction. An analytical expression for the basic reproduction number of the spatial model, R ~ 0 , is derived as a function of the non-spatial counterpart, R ₀, which characterizes a transition between a disease-free and a propagation phase, in which the disease propagates over a large fraction of the system.

Modeling epidemic flow with fluid dynamics

In this paper, a new mathematical model based on partial differential equations is proposed to study the spatial spread of infectious diseases. The model incorporates fluid dynamics theory and represents the epidemic spread as a fluid motion generated through the interaction between the susceptible and infected hosts. At the macroscopic level, the spread of the infection is modeled as an inviscid flow described by the Euler equation. Nontrivial numerical methods from computational fluid dynamics (CFD) are applied to investigate the model. In particular, a fifth-order weighted essentially non-oscillatory (WENO) scheme is employed for the spatial discretization. As an application, this mathematical and computational framework is used in a simulation study for the COVID-19 outbreak in Wuhan, China. The simulation results match the reported data for the cumulative cases with high accuracy and generate new insight into the complex spatial dynamics of COVID-19.

Cholera

Cholera was first described in the areas around the Bay of Bengal and spread globally, resulting in seven pandemics during the past two centuries. It is caused by toxigenic Vibrio cholerae O1 or O139 bacteria. Cholera is characterised by mild to potentially fatal acute watery diarrhoeal disease. Prompt rehydration therapy is the cornerstone of management. We present an overview of cholera and its pathogenesis, natural history, bacteriology, and epidemiology, while highlighting advances over the past 10 years in molecular epidemiology, immunology, and vaccine development and deployment. Since 2014, the Global Task Force on Cholera Control, a WHO coordinated network of partners, has been working with several countries to develop national cholera control strategies. The global roadmap for cholera control focuses on stopping transmission in cholera hotspots through vaccination and improved water, sanitation, and hygiene, with the aim to reduce cholera deaths by 90% and eliminate local transmission in at least 20 countries by 2030.

A reaction-advection-diffusion model of cholera epidemics with seasonality and human behavior change

Cholera is a water- and food-borne infectious disease caused by V. cholerae. To investigate multiple effects of human behavior change, seasonality and spatial heterogeneity on cholera spread, we propose a reaction-advection-diffusion model that incorporates human hosts and aquatic reservoir of V. cholerae. We derive the basic reproduction number [Formula: see text] for this system and then establish a threshold type result on its global dynamics in terms of [Formula: see text]. Further, we show that the bacterial loss at the downstream end of the river due to water flux can reduce the disease risk, and describe the asymptotic behavior of [Formula: see text] for small and large diffusion in a special case (where the diffusion rates of infected human and the pathogen are constant). We also study the transmission dynamics at the early stage of cholera outbreak numerically, and find that human behavior change may lower the infection level and delay the disease peak. Moreover, the relative rate of bacterial loss, together with convection rate, plays an important role in identifying the severely infected areas. Meanwhile spatial heterogeneity may dilute or amplify cholera infection, which in turn would increase the complexity of disease spread.

Epidemicity of cholera spread and the fate of infection control measures

The fate of ongoing infectious disease outbreaks is predicted through reproduction numbers, defining the long-term establishment of the infection, and epidemicity indices, tackling the reactivity of the infectious pool to new contagions. Prognostic metrics of unfolding outbreaks are of particular importance when designing adaptive emergency interventions facing real-time assimilation of epidemiological evidence. Our aim here is twofold. First, we propose a novel form of the epidemicity index for the characterization of cholera epidemics in spatial models of disease spread. Second, we examine in hindsight the survey of infections, treatments and containment measures carried out for the now extinct 2010–2019 Haiti cholera outbreak, to suggest that magnitude and timing of non-pharmaceutical and vaccination interventions imply epidemiological responses recapped by the evolution of epidemicity indices. Achieving negative epidemicity greatly accelerates fading of infections and thus proves a worthwhile target of containment measures. We also show that, in our model, effective reproduction numbers and epidemicity indices are explicitly related. Therefore, providing an upper bound to the effective reproduction number (significantly lower than the unit threshold) warrants negative epidemicity and, in turn, a rapidly fading outbreak preventing coalescence of sparse local sub-threshold flare-ups.

Structural identifiability analysis of age-structured PDE epidemic models

Computational and mathematical models rely heavily on estimated parameter values for model development. Identifiability analysis determines how well the parameters of a model can be estimated from experimental data. Identifiability analysis is crucial for interpreting and determining confidence in model parameter values and to provide biologically relevant predictions. Structural identifiability analysis, in which one assumes data to be noiseless and arbitrarily fine-grained, has been extensively studied in the context of ordinary differential equation (ODE) models, but has not yet been widely explored for age-structured partial differential equation (PDE) models. These models present additional difficulties due to increased number of variables and partial derivatives as well as the presence of boundary conditions. In this work, we establish a pipeline for structural identifiability analysis of age-structured PDE models using a differential algebra framework and derive identifiability results for specific age-structured models. We use epidemic models to demonstrate this framework because of their wide-spread use in many different diseases and for the corresponding parallel work previously done for ODEs. In our application of the identifiability analysis pipeline, we focus on a Susceptible-Exposed-Infected model for which we compare identifiability results for a PDE and corresponding ODE system and explore effects of age-dependent parameters on identifiability. We also show how practical identifiability analysis can be applied in this example.

Threshold dynamics of a nonlocal and delayed cholera model in a spatially heterogeneous environment

A nonlocal and delayed cholera model with two transmission mechanisms in a spatially heterogeneous environment is derived. We introduce two basic reproduction numbers, one is for the bacterium in the environment and the other is for the cholera disease in the host population. If the basic reproduction number for the cholera bacterium in the environment is strictly less than one and the basic reproduction number of infection is no more than one, we prove globally asymptotically stability of the infection-free steady state. Otherwise, the infection will persist and there exists at least one endemic steady state. For the special homogeneous case, the endemic steady state is actually unique and globally asymptotically stable. Under some conditions, the basic reproduction number of infection is strictly decreasing with respect to the diffusion coefficients of cholera bacteria and infectious hosts. When these conditions are violated, numerical simulation suggests that spatial diffusion may not only spread the infection from high-risk region to low-risk region, but also increase the infection level in high-risk region.

A spatial epidemic model with a moving boundary

We present a new mathematical model to investigate the spatial spread of an infectious disease. The model consists of a nonlinear PDE system with an unknown velocity field, defined on an epidemic domain that changes with time. The moving boundary of the domain represents the wavefront of the epidemic. We conduct an equilibrium analysis to the simplified models represented by ODE systems. We also perform a numerical study on the original PDE system for a range of scenarios, including one under a realistic epidemic setting.

Tracking the time course of reproduction number and lockdown's effect on human behaviour during SARS-CoV-2 epidemic: nonparametric estimation

Understanding the SARS-CoV-2 dynamics has been subject of intense research in the last months. In particular, accurate modeling of lockdown effects on human behaviour and epidemic evolution is a key issue in order e.g. to inform health-care decisions on emergency management. In this regard, the compartmental and spatial models so far proposed use parametric descriptions of the contact rate, often assuming a time-invariant effect of the lockdown. In this paper we show that these assumptions may lead to erroneous evaluations on the ongoing pandemic. Thus, we develop a new class of nonparametric compartmental models able to describe how the impact of the lockdown varies in time. Our estimation strategy does not require significant Bayes prior information and exploits regularization theory. Hospitalized data are mapped into an infinite-dimensional space, hence obtaining a function which takes into account also how social distancing measures and people's growing awareness of infection's risk evolves as time progresses. This also permits to reconstruct a continuous-time profile of SARS-CoV-2 reproduction number with a resolution never reached before in the literature. When applied to data collected in Lombardy, the most affected Italian region, our model illustrates how people behaviour changed during the restrictions and its importance to contain the epidemic. Results also indicate that, at the end of the lockdown, around [Formula: see text] of people in Lombardy and [Formula: see text] in Italy was affected by SARS-CoV-2, with the fatality rate being 1.14%. Then, we discuss how the situation evolved after the end of the lockdown showing that the reproduction number dangerously increased in the summer, due to holiday relax, reaching values larger than one on August 1, 2020. Finally, we also document how Italy faced the second wave of infection in the last part of 2020. Since several countries still observe a growing epidemic and others could be subject to other waves, the proposed reproduction number tracking methodology can be of great help to health care authorities to prevent SARS-CoV-2 diffusion or to assess the impact of lockdown restrictions on human behaviour to contain the spread.

Stochastic models of infectious diseases in a periodic environment with application to cholera epidemics

Seasonal variation affects the dynamics of many infectious diseases including influenza, cholera and malaria. The time when infectious individuals are first introduced into a population is crucial in predicting whether a major disease outbreak occurs. In this investigation, we apply a time-nonhomogeneous stochastic process for a cholera epidemic with seasonal periodicity and a multitype branching process approximation to obtain an analytical estimate for the probability of an outbreak. In particular, an analytic estimate of the probability of disease extinction is shown to satisfy a system of ordinary differential equations which follows from the backward Kolmogorov differential equation. An explicit expression for the mean (resp. variance) of the first extinction time given an extinction occurs is derived based on the analytic estimate for the extinction probability. Our results indicate that the probability of a disease outbreak, and mean and standard derivation of the first time to disease extinction are periodic in time and depend on the time when the infectious individuals or free-living pathogens are introduced. Numerical simulations are then carried out to validate the analytical predictions using two examples of the general cholera model. At the end, the developed theoretical results are extended to more general models of infectious diseases.

Predicting Oil Production Sites for Planning Road Infrastructure: Trip Generation Using SIR Epidemic Model

Drilling activity produces a significant amount of road traffic through unpaved and paved local roads. Because oil production is an important contributor to the local economy in the state of North Dakota, the state and local transportation agencies make efforts to support local energy logistics through the expansion and good repair and maintenance of transportation infrastructure. As part of this effort, it is important to build new roads and bridges, maintain existing road pavement and non-marked road surface conditions, and improve bridge and other transportation infrastructure. Therefore, the purpose of this study is to review previous oil location prediction models and propose a novel geospatial model to predict drilling locations which have a significant impact on local roads, to verify and provide a better prediction model. Then, this study proposes a SIR (susceptible–infected–recovered) epidemic model to predict oil drilling locations which are traffic generators. The simulation has been done on the historical data from 1980 to 2015. The study found that the best fit parameters of β (contact rate) and μ (recovery rate) were estimated by using a dataset of historical oil wells. The study found that the SIR epidemic model can be applied to predict the locations of oil wells. The proposed model can be used to predict other drilling locations and can assist with traffic, road conditions, and other related issues, which is a much needed predictive model that is key in transportation planning and pavement design and maintenance.

COVID-19 epidemic control using short-term lockdowns for collective gain

While many efforts are currently devoted to vaccines development and administration, social distancing measures, including severe restrictions such as lockdowns, remain fundamental tools to contain the spread of COVID-19. A crucial point for any government is to understand, on the basis of the epidemic curve, the right temporal instant to set up a lockdown and then to remove it. Different strategies are being adopted with distinct shades of intensity. USA and Europe tend to introduce restrictions of considerable temporal length. They vary in time: a severe lockdown may be reached and then gradually relaxed. An interesting alternative is the Australian model where short and sharp responses have repeatedly tackled the virus and allowed people a return to near normalcy. After a few positive cases are detected, a lockdown is immediately set. In this paper we show that the Australian model can be generalized and given a rigorous mathematical analysis, casting strategies of the type short-term pain for collective gain in the context of sliding-mode control, an important branch of nonlinear control theory. This allows us to gain important insights regarding how to implement short-term lockdowns, obtaining a better understanding of their merits and possible limitations. Effects of vaccines administration in improving the control law's effectiveness are also illustrated. Our model predicts the duration of the severe lockdown to be set to maintain e.g. the number of people in intensive care under a certain threshold. After tuning our strategy exploiting data collected in Italy, it turns out that COVID-19 epidemic could be e.g. controlled by alternating one or two weeks of complete lockdown with one or two months of freedom, respectively. Control strategies of this kind, where the lockdown's duration is well circumscribed, could be important also to alleviate coronavirus impact on economy.

Modeling and analysis of different scenarios for the spread of COVID-19 by using the modified multi-agent systems - Evidence from the selected countries

Currently, there is a global pandemic of COVID-19. To assess its prevalence, it is necessary to have adequate models that allow real-time modeling of the impact of various quarantine measures by the state. The SIR model, which is implemented using a multi-agent system based on mobile cellular automata, was improved. The paper suggests ways to improve the rules of the interaction and behavior of agents. Methods of comparing the parameters of the SIR model with real geographical, social and medical indicators have been developed. That allows the modeling of the spatial distribution of COVID-19 as a single location and as the whole country consisting of individual regions that interact with each other by transport, taking into account factors such as public transport, supermarkets, schools, universities, gyms, churches, parks. The developed model also allows us to assess the impact of quarantine, restrictions on transport connections between regions, to take into account such factors as the incubation period, the mask regime, maintaining a safe distance between people, and so on. A number of experiments were conducted in the work, which made it possible to assess both the impact of individual measures to stop the pandemic and their comprehensive application. A method of comparing computer-time and dynamic parameters of the model with real data is proposed, which allowed assessing the effectiveness of the government in stopping the pandemic in the Chernivtsi region, Ukraine. A simulation of the pandemic spread in countries such as Slovakia, Turkey and Serbia was also conducted. The calculations showed the high-accuracy matching of the forecast model with real data.

A hybrid multi-scale model of COVID-19 transmission dynamics to assess the potential of non-pharmaceutical interventions

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel coronavirus that emerged in Wuhan, China in December 2019. It has caused a global outbreak which represents a major threat to global health. Public health resorted to non-pharmaceutical interventions such as social distancing and lockdown to slow down the spread of the pandemic. However, the effect of each of these measures remains hard to quantify. We design a multi-scale model that simulates the transmission dynamics of COVID-19. We describe the motion of individual agents using a social force model. Each agent can be either susceptible, infected, quarantined, immunized or deceased. The model considers both mechanisms of direct and indirect transmission. We parameterize the model to reproduce the early dynamics of disease spread in Italy. We show that panic situations increase the risk of infection transmission in crowds despite social distancing measures. Next, we reveal that pre-symptomatic transmission accelerates the onset of the exponential growth of cases. After that, we demonstrate that the persistence of SARS-CoV-2 on hard surfaces determines the number of cases reached during the peak of the epidemic. Then, we show that the restricted movement of the individuals flattens the epidemic curve. Finally, model predictions suggest that measures stricter than social distancing and lockdown were used to control the epidemic in Wuhan, China.

Basic Reproduction Numbers for a Class of Reaction-Diffusion Epidemic Models

We study the basic reproduction numbers for a class of reaction-diffusion epidemic models that are developed from autonomous ODE systems. We present a general numerical framework to compute such basic reproduction numbers; meanwhile, the numerical formulation provides useful insight into their characterizations. Using matrix analysis, we show that the basic reproduction numbers are the same for these PDE models and their associated ODE models in several important cases that include, among others, a single infected compartment, constant diffusion rates, uniform diffusion patterns among the infected compartments, and partial diffusion in the system.

Spatially explicit models for exploring COVID-19 lockdown strategies

This article describes two spatially explicit models created to allow experimentation with different societal responses to the COVID-19 pandemic. We outline the work to date on modeling spatially explicit infective diseases and show that there are gaps that remain important to fill. We demonstrate how geographical regions, rather than a single, national approach, are likely to lead to better outcomes for the population. We provide a full account of how our models function, and how they can be used to explore many different aspects of contagion, including: experimenting with different lockdown measures, with connectivity between places, with the tracing of disease clusters, and the use of improved contact tracing and isolation. We provide comprehensive results showing the use of these models in given scenarios, and conclude that explicitly regionalized models for mitigation provide significant advantages over a "one-size-fits-all" approach. We have made our models, and their data, publicly available for others to use in their own locales, with the hope of providing the tools needed for geographers to have a voice during this difficult time.

Modeling and analyzing cholera transmission dynamics with vaccination age

A new mathematical model is formulated to investigate the transmission dynamics of cholera under vaccination, with a focus on the impact of vaccination age. The basic reproduction number is derived and proved to be a sharp control threshold determining whether or not the infection is persistent. We conduct a rigorous analysis on the local and global stability properties of the equilibria in system. Meanwhile, we compare the results to those of the simplified model based on ordinary differential equations where the effects of vaccination age are not incorporated. Numerical simulation results verify our analytical findings.

Generation and application of river network analogues for use in ecology and evolution

Several key processes in freshwater ecology are governed by the connectivity inherent to dendritic river networks. These have extensively been analyzed from a geomorphological and hydrological viewpoint, yet structures classically used in ecological modeling have been poorly representative of the structure of real river basins, often failing to capture well-known scaling features of natural rivers. Pioneering work identified optimal channel networks (OCNs) as spanning trees reproducing all scaling features characteristic of natural stream networks worldwide. While OCNs have been used to create landscapes for studies on metapopulations, biodiversity, and epidemiology, their generation has not been generally accessible.Given the increasing interest in dendritic riverine networks by ecologists and evolutionary biologists, we here present a method to generate OCNs and, to facilitate its application, we provide the R-package OCNet. Owing to the stochastic process generating OCNs, multiple network replicas spanning the same surface can be built; this allows performing computational experiments whose results are irrespective of the particular shape of a single river network. The OCN construct also enables the generation of elevational gradients derived from the optimal network configuration, which can constitute three-dimensional landscapes for spatial studies in both terrestrial and freshwater realms. Moreover, the package provides functions that aggregate OCNs into an arbitrary number of nodes, calculate several descriptors of river networks, and draw relevant network features.We describe the main functionalities of the package and its integration with other R-packages commonly used in spatial ecology. Moreover, we exemplify the generation of OCNs and discuss an application to a metapopulation model for an invasive riverine species.In conclusion, OCNet provides a powerful tool to generate realistic river network analogues for various applications. It thereby allows the design of spatially realistic studies in increasingly impacted ecosystems and enhances our knowledge on spatial processes in freshwater ecology in general.

Microbiome and Functional Analysis of a Traditional Food Process: Isolation of a Novel Species (Vibrio hibernica) With Industrial Potential

Traditional food preservation processes are vital for the food industry. They not only preserve a high-quality protein and nutrient source but can also provide important value-added organoleptic properties. The Wiltshire process is a traditional food curing method applied to meat, and special recognition is given to the maintenance of a live rich microflora within the curing brine. We have previously analyzed a curing brine from this traditional meat process and characterized a unique microbial core signature. The characteristic microbial community is actively maintained and includes the genera, Marinilactibacillus, Carnobacterium, Leuconostoc, and Vibrio. The bacteria present are vital for Wiltshire curing compliance. However, the exact function of this microflora is largely unknown. A microbiome profiling of three curing brines was conducted and investigated for functional traits by the robust bioinformatic tool, Tax4Fun. The key objective was to uncover putative metabolic functions associated with the live brine and to identify changes over time. The functional bioinformatic analysis revealed metabolic enrichments over time, with many of the pathways identified as being involved in organoleptic development. The core bacteria present in the brine are Lactic Acid Bacteria (LAB), with the exception of the Vibrio genus. LAB are known for their positive contribution to food processing, however, little work has been conducted on the use of Vibrio species for beneficial processes. The Vibrio genome was sequenced by Illumina MiSeq technologies and annotated in RAST. A phylogenetic reconstruction was completed using both the 16S rRNA gene and housekeeping genes, gapA, ftsZ, mreB, topA, gyrB, pyrH, recA, and rpoA. The isolated Vibrio species was defined as a unique novel species, named Vibrio hibernica strain B1.19. Metabolic profiling revealed that the bacterium has a unique substrate scope in comparison to other closely related Vibrio species tested. The possible function and industrial potential of the strain was investigated using carbohydrate metabolizing profiling under food processing relevant conditions. Vibrio hibernica is capable of metabolizing a unique carbohydrate profile at low temperatures. This characteristic provides new application options for use in the industrial food sector, as well as highlighting the key role of this bacterium in the Wiltshire curing process.

Direct transmission via households informs models of disease and intervention dynamics in cholera

While several basic properties of cholera outbreaks are common to most settings-the pathophysiology of the disease, the waterborne nature of transmission, and others-recent findings suggest that transmission within households may play a larger role in cholera outbreaks than previously appreciated. Important features of cholera outbreaks have long been effectively modeled with mathematical and computational approaches, but little is known about how variation in direct transmission via households may influence epidemic dynamics. In this study, we construct a mathematical model of cholera that incorporates transmission within and between households. We observe that variation in the magnitude of household transmission changes multiple features of disease dynamics, including the severity and duration of outbreaks. Strikingly, we observe that household transmission influences the effectiveness of possible public health interventions (e.g. water treatment, antibiotics, vaccines). We find that vaccine interventions are more effective than water treatment or antibiotic administration when direct household transmission is present. Summarizing, we position these results within the landscape of existing models of cholera, and speculate on its implications for epidemiology and public health.

THE COMPUTATION OF REPRODUCTION NUMBERS FOR THE ENVIRONMENT-HOST-ENVIRONMENT CHOLERA TRANSMISSION DYNAMICS

This study presents a new model for the environment-host-environment transmission dynamics of V. cholerae in a community with an interconnected aquatic pond–river water network. For the case when the human host is the sole target of anti-cholera control and the volume of water in the pond is maximum, the disease-free equilibrium of the model is shown to be globally asymptotically stable whenever a certain epidemiological threshold, known as the basic reproduction number [Formula: see text], is less than unity. The epidemiological implication of this result is that cholera can be eliminated from the community if the control strategies implemented can bring (and maintain) [Formula: see text] to a value less than unity. Four scenarios, that represent different interpretations of the role of the V. cholerea pathogen within the environment, were studied. The corresponding basic reproduction numbers were shown to exhibit the same threshold property with respect to the value unity (i.e., if one is less (equal, greater) than unity, then the three others are also less (equal, greater) than unity. Further, it was shown that for the case where anti-cholera control is focused on the human host population, the associated type reproduction number of the model (corresponding to each of the four transmission scenarios considered) is unique. The implication of this result is that the estimate of the effort needed for disease elimination (i.e., the required herd immunity threshold) is unique, regardless of which of the four transmission scenarios is considered. However, when any of the other two bacterial population types in the aquatic environment (i.e., bacterial in the pond or river) is the focus of the control efforts, this study shows that the associated type reproduction number is not unique. Extensive numerical simulations of the model, using a realistic set of parameters from the published literature, show that the community-wide implementation of a strategy that focus on improved water quality, sanitation, and hygiene (known as WASH-only strategy), using the current estimated coverage of 50% and efficacy of 60%, is unable to lead to the elimination of the disease. Such elimination is attainable if the coverage and efficacy are increased (e.g., to 80% and 90%, respectively). Further, elimination can be achieved using a strategy that focuses on oral rehydration therapy and the use of antibiotics to treat the infected humans (i.e., treatment-only strategy) for moderate effectiveness and coverage levels. The combined hybrid WASH-treatment strategy provides far better population-level impact vis a vis disease elimination. This study ranks the three interventions in the following order of population-level effectiveness: combined WASH-treatment, followed by treatment-only and then WASH-only strategy.

Conditions for transient epidemics of waterborne disease in spatially explicit systems

Waterborne diseases are a diverse family of infections transmitted through ingestion of-or contact with-water infested with pathogens. Outbreaks of waterborne infections often show well-defined spatial signatures that are typically linked to local eco-epidemiological conditions, water-mediated pathogen transport and human mobility. In this work, we apply a spatially explicit network model describing the transmission cycle of waterborne pathogens to determine invasion conditions in metacommunities endowed with a realistic spatial structure. Specifically, we aim to define conditions under which pathogens can temporarily colonize a set of human communities, thus triggering a transient epidemic outbreak. To that end, we apply generalized reactivity analysis, a recently developed methodological framework for the study of transient dynamics in ecological systems subject to external perturbations. The study of pathogen invasion is complemented by the detection of the spatial signatures associated with the perturbations to a disease-free system that are expected to be amplified the most over different time scales. Understanding the drivers of waterborne disease dynamics over time scales that are relevant to epidemic and/or endemic transmission is a crucial, cross-disciplinary challenge, as large portions of the developing world still struggle to cope with the burden of these infections.

Rainfall as a driver of epidemic cholera: Comparative model assessments of the effect of intra-seasonal precipitation events

The correlation between cholera epidemics and climatic drivers, in particular seasonal tropical rainfall, has been studied in a variety of contexts owing to its documented relevance. Several mechanistic models of cholera transmission have included rainfall as a driver by focusing on two possible transmission pathways: either by increasing exposure to contaminated water (e.g. due to worsening sanitary conditions during water excess), or water contamination by freshly excreted bacteria (e.g. due to washout of open-air defecation sites or overflows). Our study assesses the explanatory power of these different modeling structures by formal model comparison using deterministic and stochastic models of the type susceptible-infected-recovered-bacteria (SIRB). The incorporation of rainfall effects is generalized using a nonlinear function that can increase or decrease the relative importance of the large precipitation events. Our modelling framework is tested against the daily epidemiological data collected during the 2015 cholera outbreak within the urban context of Juba, South Sudan. This epidemic is characterized by a particular intra-seasonal double peak on the incidence in apparent relation with particularly strong rainfall events. Our results show that rainfall-based models in both their deterministic and stochastic formulations outperform models that do not account for rainfall. In fact, classical SIRB models are not able to reproduce the second epidemiological peak, thus suggesting that it was rainfall-driven. Moreover we found stronger support across model types for rainfall acting on increased exposure rather than on exacerbated water contamination. Although these results are context-specific, they stress the importance of a systematic and comprehensive appraisal of transmission pathways and their environmental forcings when embarking in the modelling of epidemic cholera.

A Game-Theoretic Model of Cholera with Optimal Personal Protection Strategies

Cholera is an acute gastro-intestinal infection that affects millions of people throughout the world each year, primarily but not exclusively in developing countries. Because of its public health ramifications, considerable mathematical attention has been paid to the disease. Here we consider one neglected aspect of combating cholera: personal participation in anti-cholera interventions. We construct a game-theoretic model of cholera in which individuals choose whether to participate in either vaccination or clean water consumption programs under assumed costs. We find that relying upon individual compliance significantly lowers the incidence of the disease as long as the cost of intervention is sufficiently low, but does not eliminate it. The relative costs of the measures determined whether a population preferentially adopts a single preventative measure or employs the measure with the strongest early adoption.

Near real-time forecasting for cholera decision making in Haiti after Hurricane Matthew

Computational models of cholera transmission can provide objective insights into the course of an ongoing epidemic and aid decision making on allocation of health care resources. However, models are typically designed, calibrated and interpreted post-hoc. Here, we report the efforts of a team from academia, field research and humanitarian organizations to model in near real-time the Haitian cholera outbreak after Hurricane Matthew in October 2016, to assess risk and to quantitatively estimate the efficacy of a then ongoing vaccination campaign. A rainfall-driven, spatially-explicit meta-community model of cholera transmission was coupled to a data assimilation scheme for computing short-term projections of the epidemic in near real-time. The model was used to forecast cholera incidence for the months after the passage of the hurricane (October-December 2016) and to predict the impact of a planned oral cholera vaccination campaign. Our first projection, from October 29 to December 31, predicted the highest incidence in the departments of Grande Anse and Sud, accounting for about 45% of the total cases in Haiti. The projection included a second peak in cholera incidence in early December largely driven by heavy rainfall forecasts, confirming the urgency for rapid intervention. A second projection (from November 12 to December 31) used updated rainfall forecasts to estimate that 835 cases would be averted by vaccinations in Grande Anse (90% Prediction Interval [PI] 476-1284) and 995 in Sud (90% PI 508-2043). The experience gained by this modeling effort shows that state-of-the-art computational modeling and data-assimilation methods can produce informative near real-time projections of cholera incidence. Collaboration among modelers and field epidemiologists is indispensable to gain fast access to field data and to translate model results into operational recommendations for emergency management during an outbreak. Future efforts should thus draw together multi-disciplinary teams to ensure model outputs are appropriately based, interpreted and communicated.

River networks as ecological corridors: A coherent ecohydrological perspective

This paper draws together several lines of argument to suggest that an ecohydrological framework, i.e. laboratory, field and theoretical approaches focused on hydrologic controls on biota, has contributed substantially to our understanding of the function of river networks as ecological corridors. Such function proves relevant to: the spatial ecology of species; population dynamics and biological invasions; the spread of waterborne disease. As examples, we describe metacommunity predictions of fish diversity patterns in the Mississippi-Missouri basin, geomorphic controls imposed by the fluvial landscape on elevational gradients of species' richness, the zebra mussel invasion of the same Mississippi-Missouri river system, and the spread of proliferative kidney disease in salmonid fish. We conclude that spatial descriptions of ecological processes in the fluvial landscape, constrained by their specific hydrologic and ecological dynamics and by the ecosystem matrix for interactions, i.e. the directional dispersal embedded in fluvial and host/pathogen mobility networks, have already produced a remarkably broad range of significant results. Notable scientific and practical perspectives are thus open, in the authors' view, to future developments in ecohydrologic research.

Epidemiological link of a major cholera outbreak in Greater Accra region of Ghana, 2014

BACKGROUND

Cholera remains an important public health challenge globally. Several pandemics have occurred in different parts of the world and have been epidemiologically linked by different researchers to illustrate how the cases were spread and how they were related to index cases. Even though the risk factors associated with the 2014 cholera outbreak were investigated extensively, the link between index cases and the source of infection was not investigated to help break the transmission process. This study sought to show how the index cases from various districts of the Greater Accra Region may have been linked.

METHODS

We carried out a descriptive cross sectional study to investigate the epidemiological link of the 2014 cholera outbreak in the Greater Accra region of Ghana. An extensive review of all district records on cholera cases in the Greater Accra region was carried out. Index cases were identified with the help of line lists. Univariate analyses were expressed as frequency distributions, percentages, mean ± Standard Deviation, and rates (attack rates, case-fatality rates etc.) as appropriate. Maps were drawn using Arc GIS and Epi info software to describe the pattern of transmission.

RESULTS

Up to 20,199 cholera cases were recorded. Sixty percent of the cases were between 20 and 40 years and about 58% (11,694) of the total cases were males. Almost 50% of the cases occurred in the Accra Metro district. Two-thirds of the index cases ate food prepared outside their home and had visited the Accra Metropolis.

CONCLUSIONS

The 2014 cholera outbreak can be described as a propagated source outbreak linked to the Accra Metropolis. The link between index cases and the source of infection, if investigated earlier could have helped break the transmission process. Such investigations also inform decision-making about the appropriate interventions to be instituted to prevent subsequent outbreaks.

Climate-driven endemic cholera is modulated by human mobility in a megacity

Although a differential sensitivity of cholera dynamics to climate variability has been reported in the spatially heterogeneous megacity of Dhaka, Bangladesh, the specific patterns of spread of the resulting risk within the city remain unclear. We build on an established probabilistic spatial model to investigate the importance and role of human mobility in modulating spatial cholera transmission. Mobility fluxes were inferred using a straightforward and generalizable methodology that relies on mapping population density based on a high resolution urban footprint product, and a parameter-free human mobility model. In accordance with previous findings, we highlight the higher sensitivity to the El Niño Southern Oscillation (ENSO) in the highly populated urban center than in the more rural periphery. More significantly, our results show that cholera risk is largely transmitted from the climate-sensitive core to the periphery of the city, with implications for the planning of control efforts. In addition, including human mobility improves the outbreak prediction performance of the model with an 11 month lead. The interplay between climatic and human mobility factors in cholera transmission is discussed from the perspective of the rapid growth of megacities across the developing world.

Modeling Key Drivers of Cholera Transmission Dynamics Provides New Perspectives for Parasitology

Hydroclimatological and anthropogenic factors are key drivers of waterborne disease transmission. Information on human settlements and host mobility on waterways along which pathogens and hosts disperse, and relevant hydroclimatological processes, can be acquired remotely and included in spatially explicit mathematical models of disease transmission. In the case of epidemic cholera, such models allowed the description of complex disease patterns and provided insight into the course of ongoing epidemics. The inclusion of spatial information in models of disease transmission can aid in emergency management and the assessment of alternative interventions. Here, we review the study of drivers of transmission via spatially explicit approaches and argue that, because many parasitic waterborne diseases share the same drivers as cholera, similar principles may apply.

Dose-response relationships for environmentally mediated infectious disease transmission models

Environmentally mediated infectious disease transmission models provide a mechanistic approach to examining environmental interventions for outbreaks, such as water treatment or surface decontamination. The shift from the classical SIR framework to one incorporating the environment requires codifying the relationship between exposure to environmental pathogens and infection, i.e. the dose-response relationship. Much of the work characterizing the functional forms of dose-response relationships has used statistical fit to experimental data. However, there has been little research examining the consequences of the choice of functional form in the context of transmission dynamics. To this end, we identify four properties of dose-response functions that should be considered when selecting a functional form: low-dose linearity, scalability, concavity, and whether it is a single-hit model. We find that i) middle- and high-dose data do not constrain the low-dose response, and different dose-response forms that are equally plausible given the data can lead to significant differences in simulated outbreak dynamics; ii) the choice of how to aggregate continuous exposure into discrete doses can impact the modeled force of infection; iii) low-dose linear, concave functions allow the basic reproduction number to control global dynamics; and iv) identifiability analysis offers a way to manage multiple sources of uncertainty and leverage environmental monitoring to make inference about infectivity. By applying an environmentally mediated infectious disease model to the 1993 Milwaukee Cryptosporidium outbreak, we demonstrate that environmental monitoring allows for inference regarding the infectivity of the pathogen and thus improves our ability to identify outbreak characteristics such as pathogen strain.

Global stability and uniform persistence of the reaction-convection-diffusion cholera epidemic model

We study the global stability issue of the reaction-convection-diffusion cholera epidemic PDE model and show that the basic reproduction number serves as a threshold parameter that predicts whether cholera will persist or become globally extinct. Specifically, when the basic reproduction number is beneath one, we show that the disease-free-equilibrium is globally attractive. On the other hand, when the basic reproduction number exceeds one, if the infectious hosts or the concentration of bacteria in the contaminated water are not initially identically zero, we prove the uniform persistence result and that there exists at least one positive steady state.

Model distinguishability and inference robustness in mechanisms of cholera transmission and loss of immunity

Mathematical models of cholera and waterborne disease vary widely in their structures, in terms of transmission pathways, loss of immunity, and a range of other features. These differences can affect model dynamics, with different models potentially yielding different predictions and parameter estimates from the same data. Given the increasing use of mathematical models to inform public health decision-making, it is important to assess model distinguishability (whether models can be distinguished based on fit to data) and inference robustness (whether inferences from the model are robust to realistic variations in model structure). In this paper, we examined the effects of uncertainty in model structure in the context of epidemic cholera, testing a range of models with differences in transmission and loss of immunity structure, based on known features of cholera epidemiology. We fit these models to simulated epidemic and long-term data, as well as data from the 2006 Angola epidemic. We evaluated model distinguishability based on fit to data, and whether the parameter values, model behavior, and forecasting ability can accurately be inferred from incidence data. In general, all models were able to successfully fit to all data sets, both real and simulated, regardless of whether the model generating the simulated data matched the fitted model. However, in the long-term data, the best model fits were achieved when the loss of immunity structures matched those of the model that simulated the data. Two parameters, one representing person-to-person transmission and the other representing the reporting rate, were accurately estimated across all models, while the remaining parameters showed broad variation across the different models and data sets. The basic reproduction number (R₀) was often poorly estimated even using the correct model, due to practical unidentifiability issues in the waterborne transmission pathway which were consistent across all models. Forecasting efforts using noisy data were not successful early in the outbreaks, but once the epidemic peak had been achieved, most models were able to capture the downward incidence trajectory with similar accuracy. Forecasting from noise-free data was generally successful for all outbreak stages using any model. Our results suggest that we are unlikely to be able to infer mechanistic details from epidemic case data alone, underscoring the need for broader data collection, such as immunity/serology status, pathogen dose response curves, and environmental pathogen data. Nonetheless, with sufficient data, conclusions from forecasting and some parameter estimates were robust to variations in the model structure, and comparative modeling can help to determine how realistic variations in model structure may affect the conclusions drawn from models and data.

A REACTION–DIFFUSION MODEL FOR THE CONTROL OF CHOLERA EPIDEMIC

Understanding the spatio-temporal dynamics of cholera outbreaks may help in devising more effective control procedures. In this paper, we have considered a reaction–diffusion system for biological control of cholera epidemic. Firstly, we have focused on temporal evolution of cholera in a region and its control using lytic bacteriophage in the aquatic reservoirs. Then, we have explored the effect of spatial dispersion of populations on the disease dynamics. We have observed the onset of sustained oscillations via Hopf-bifurcation for the endemic state of temporal system. This onset of fluctuations in populations depends upon the phage adsorption rate. But in the spatially extended setting, all the populations stabilize i.e., the spatio-temporal distribution of all the populations becomes uniform. Some numerical computations have been done to verify analytical results.

Modelling Population Dynamics in Realistic Landscapes with Linear Elements: A Mechanistic-Statistical Reaction-Diffusion Approach

We propose and develop a general approach based on reaction-diffusion equations for modelling a species dynamics in a realistic two-dimensional (2D) landscape crossed by linear one-dimensional (1D) corridors, such as roads, hedgerows or rivers. Our approach is based on a hybrid "2D/1D model", i.e, a system of 2D and 1D reaction-diffusion equations with homogeneous coefficients, in which each equation describes the population dynamics in a given 2D or 1D element of the landscape. Using the example of the range expansion of the tiger mosquito Aedes albopictus in France and its main highways as 1D corridors, we show that the model can be fitted to realistic observation data. We develop a mechanistic-statistical approach, based on the coupling between a model of population dynamics and a probabilistic model of the observation process. This allows us to bridge the gap between the data (3 levels of infestation, at the scale of a French department) and the output of the model (population densities at each point of the landscape), and to estimate the model parameter values using a maximum-likelihood approach. Using classical model comparison criteria, we obtain a better fit and a better predictive power with the 2D/1D model than with a standard homogeneous reaction-diffusion model. This shows the potential importance of taking into account the effect of the corridors (highways in the present case) on species dynamics. With regard to the particular case of A. albopictus, the conclusion that highways played an important role in species range expansion in mainland France is consistent with recent findings from the literature.

Disease and disaster: Optimal deployment of epidemic control facilities in a spatially heterogeneous population with changing behaviour

Epidemics of water-borne infections often follow natural disasters and extreme weather events that disrupt water management processes. The impact of such epidemics may be reduced by deployment of transmission control facilities such as clinics or decontamination plants. Here we use a relatively simple mathematical model to examine how demographic and environmental heterogeneities, population behaviour, and behavioural change in response to the provision of facilities, combine to determine the optimal configurations of limited numbers of facilities to reduce epidemic size, and endemic prevalence. We show that, if the presence of control facilities does not affect behaviour, a good general rule for responsive deployment to minimise epidemic size is to place them in exactly the locations where they will directly benefit the most people. However, if infected people change their behaviour to seek out treatment then the deployment of facilities offering treatment can lead to complex effects that are difficult to foresee. So careful mathematical analysis is the only way to get a handle on the optimal deployment. Behavioural changes in response to control facilities can also lead to critical facility numbers at which there is a radical change in the optimal configuration. So sequential improvement of a control strategy by adding facilities to an existing optimal configuration does not always produce another optimal configuration. We also show that the pre-emptive deployment of control facilities has conflicting effects. The configurations that minimise endemic prevalence are very different to those that minimise epidemic size. So cost-benefit analysis of strategies to manage endemic prevalence must factor in the frequency of extreme weather events and natural disasters.

The impact of spatial arrangements on epidemic disease dynamics and intervention strategies

The role of spatial arrangements on the spread and management strategies of a cholera epidemic is investigated. We consider the effect of human and pathogen movement on optimal vaccination strategies. A metapopulation model is used, incorporating a susceptible-infected-recovered system of differential equations coupled with an equation modelling the concentration of Vibrio cholerae in an aquatic reservoir. The model compared spatial arrangements and varying scenarios to draw conclusions on how to effectively manage outbreaks. The work is motivated by the 2010 cholera outbreak in Haiti. Results give guidance for vaccination strategies in response to an outbreak.

A Theoretical Analysis of the Geography of Schistosomiasis in Burkina Faso Highlights the Roles of Human Mobility and Water Resources Development in Disease Transmission

We study the geography of schistosomiasis across Burkina Faso by means of a spatially explicit model of water-based disease dynamics. The model quantitatively addresses the geographic stratification of disease burden in a novel framework by explicitly accounting for drivers and controls of the disease, including spatial information on the distributions of population and infrastructure, jointly with a general description of human mobility and climatic/ecological drivers. Spatial patterns of disease are analysed by the extraction and the mapping of suitable eigenvectors of the Jacobian matrix subsuming the stability of the disease-free equilibrium. The relevance of the work lies in the novel mapping of disease burden, a byproduct of the parametrization induced by regional upscaling, by model-guided field validations and in the predictive scenarios allowed by exploiting the range of possible parameters and processes. Human mobility is found to be a primary control at regional scales both for pathogen invasion success and the overall distribution of disease burden. The effects of water resources development highlighted by systematic reviews are accounted for by the average distances of human settlements from water bodies that are habitats for the parasite’s intermediate host. Our results confirm the empirical findings about the role of water resources development on disease spread into regions previously nearly disease-free also by inspection of empirical prevalence patterns. We conclude that while the model still needs refinements based on field and epidemiological evidence, the proposed framework provides a powerful tool for large-scale public health planning and schistosomiasis management.

Dynamical models of schistosomiasis infections, even spatially explicit ones, have so far only addressed spatial scales encompassing at best a few villages and the disease transmission impacts of related short-range human mobility. Here, we build from existing models of disease dynamics and spread, including a proxy of the ecology of the intermediate host of the parasite, and from generalized reproduction numbers of SIR-type systems developed for epidemics of waterborne disease, to set up large-scale projections of spatial patterns of the disease at whole country level. We ground our study in Burkina Faso in sub-Saharan Africa, and its model of social and economic development including the infrastructure built to exploit water resources, especially irrigation schemes, which have been empirically linked to enhanced disease burden. We make extensive use of remotely sensed and field data, and capitalize on ecohydrological insight. We suggest that reliable nationwide patterns of disease burden can be projected in relation to the key roles of human mobility and water resources development subsuming exposure, and claim that the case at hand provides an insightful example towards the integration of development and environmental thinking not confined to ad-hoc indicators of human development.

Influence of human behavior on cholera dynamics

This paper is devoted to studying the impact of human behavior on cholera infection. We start with a cholera ordinary differential equation (ODE) model that incorporates human behavior via modeling disease prevalence dependent contact rates for direct and indirect transmissions and infectious host shedding. Local and global dynamics of the model are analyzed with respect to the basic reproduction number. We then extend the ODE model to a reaction-convection-diffusion partial differential equation (PDE) model that accounts for the movement of both human hosts and bacteria. Particularly, we investigate the cholera spreading speed by analyzing the traveling wave solutions of the PDE model, and disease threshold dynamics by numerically evaluating the basic reproduction number of the PDE model. Our results show that human behavior can reduce (a) the endemic and epidemic levels, (b) cholera spreading speeds and (c) the risk of infection (characterized by the basic reproduction number).

Influence of Environmental Factors on Vibrio spp. in Coastal Ecosystems

Various studies have examined the relationships between vibrios and the environmental conditions surrounding them. However, very few reviews have compiled these studies into cohesive points. This may be due to the fact that these studies examine different environmental parameters, use different sampling, detection, and enumeration methodologies, and occur in diverse geographic locations. The current article is one approach to compile these studies into a cohesive work that assesses the importance of environmental determinants on the abundance of vibrios in coastal ecosystems.

On the predictive ability of mechanistic models for the Haitian cholera epidemic

Predictive models of epidemic cholera need to resolve at suitable aggregation levels spatial data pertaining to local communities, epidemiological records, hydrologic drivers, waterways, patterns of human mobility and proxies of exposure rates. We address the above issue in a formal model comparison framework and provide a quantitative assessment of the explanatory and predictive abilities of various model settings with different spatial aggregation levels and coupling mechanisms. Reference is made to records of the recent Haiti cholera epidemics. Our intensive computations and objective model comparisons show that spatially explicit models accounting for spatial connections have better explanatory power than spatially disconnected ones for short-to-intermediate calibration windows, while parsimonious, spatially disconnected models perform better with long training sets. On average, spatially connected models show better predictive ability than disconnected ones. We suggest limits and validity of the various approaches and discuss the pathway towards the development of case-specific predictive tools in the context of emergency management.

Detection of Vibrio cholerae O1 and O139 in environmental waters of rural Bangladesh: a flow-cytometry-based field trial

Presence of Vibrio cholerae serogroups O1 and O139 in the waters of the rural area of Matlab, Bangladesh, was investigated with quantitative measurements performed with a portable flow cytometer. The relevance of this work relates to the testing of a field-adapted measurement protocol that might prove useful for cholera epidemic surveillance and for validation of mathematical models. Water samples were collected from different water bodies that constitute the hydrological system of the region, a well-known endemic area for cholera. Water was retrieved from ponds, river waters, and irrigation canals during an inter-epidemic time period. Each sample was filtered and analysed with a flow cytometer for a fast determination of V. cholerae cells contained in those environments. More specifically, samples were treated with O1- and O139-specific antibodies, which allowed precise flow-cytometry-based concentration measurements. Both serogroups were present in the environmental waters with a consistent dominance of V. cholerae O1. These results extend earlier studies where V. cholerae O1 and O139 were mostly detected during times of cholera epidemics using standard culturing techniques. Furthermore, our results confirm that an important fraction of the ponds' host populations of V. cholerae are able to self-sustain even when cholera cases are scarce. Those contaminated ponds may constitute a natural reservoir for cholera endemicity in the Matlab region. Correlations of V. cholerae concentrations with environmental factors and the spatial distribution of V. cholerae populations are also discussed.