Introduction

Livestock production is a cornerstone of rural livelihoods across sub-Saharan Africa, contributing to economic stability, socio-cultural cohesion, and food security (Herrero et al., 2013). However, the sector is increasingly threatened by infectious diseases and climate change impacts – ranging from conditions intensifying vector proliferation to those exacerbating water and pasture scarcity (Stanimirova et al., 2019; Thornton et al., 2009). In this context, transhumant and nomadic practices across diverse settings serve as critical resilience strategies (Aryal et al., 2018; Motta et al., 2018), allowing pastoral communities to respond to seasonal resource challenges. In the West African Sahel, these transhumant movements typically lead to resource-rich regions along major rivers where abundance of water and high-quality pasture support large populations of livestock (Bassett & Turner, 2007; Freudenberger & Freudenberger, 1993). Despite their importance, these movements are poorly understood due to their largely unregulated nature and limited documentation (Belkhiria et al., 2019; Erdaw, 2023). The association between livestock movements and the spread of infectious diseases has long been recognised (Fevre et al., 2006; Kim et al., 2021; Prentice et al., 2017), allowing introduction of pathogens to local susceptible populations. The spread of infectious diseases is further exacerbated by the very limited veterinary surveillance and disease prevention services available to nomadic and transhumant pastoralists (Belkhiria et al., 2019; Schelling et al., 2008), who often inhabit remote and inaccessible regions.

The epidemiology of Rift Valley fever (RVF), a WHO priority zoonotic disease caused by the Rift Valley fever virus (RVFV) [Family: Phenuiviridae], exemplifies the interplay between ecological conditions, vector activity, and livestock mobility (Durand et al., 2020; Walsh et al., 2017). RVFV is primarily transmitted by mosquitoes, with competent species identified across the Aedes, Culex, and Anopheles genera (Linthicum et al., 2016; Lumley et al., 2017). In addition to vector transmission, both human and ruminant infections can also occur through direct contact with infectious tissues, and fluids such as blood and milk during outbreaks (Gerken et al., 2022). Outbreaks in livestock are characterised by devastating abortions and neonatal mortalities (Ksiazek et al., 1989; Munyua et al., 2010). Human cases range from asymptomatic or mild flu-like symptoms to a severe and sometimes fatal syndrome, with complications including haemorrhagic fever or encephalitis (Hoogstraal et al., 1979; Madani et al., 2003). Since its discovery in Kenya in 1931 (Daubney et al., 1931), RVF has caused repeated outbreaks across sub-Saharan Africa, with notable expansions into Egypt (Hoogstraal et al., 1979), the Arabian Peninsula (Balkhy & Memish, 2003) and Comoros Archipelago (Sissoko et al., 2009), underscoring its global threat potential.

In The Gambia, situated in the Sahelian zone of West Africa, the livestock sector contributes approximately 30% of the agricultural GDP and 10% of the national GDP (Rich et al., 2020). The majority of livestock management practices are traditional, with 98% of livestock-owning households engaged in either pastoral (entirely dependent on livestock production) or semi-pastoral systems (integrating crop husbandry and livestock rearing) (Department of Livestock Services, 2016, pers. comms.). The floodplains of the Gambia river serve as a perennial resource hub, providing water and pasture that attract transhumant herds migrating from the Sudano-Sahelian zone during the prolonged harsh dry season (November – June). Additionally, riverine ecosystems offers ideal breeding sites for RVFV-competent mosquito species in the Culex and Mansonia genera (Diallo et al., 2000; Saluzzo et al., 1984; Snow, 1983), supporting year-round mosquito proliferation. The Gambia river floodplains thus present an environment conducive to the spread of both vector-borne and non-vector-borne livestock diseases. RVF cases affecting both humans and livestock were first reported in The Gambia in 2002 (FAO, 2003), however, surveillance remains limited, leaving critical gaps in understanding of its epidemiology. Our recent serological investigation revealed ongoing RVFV transmission within ruminant livestock, and highlighted the importance of transhumant movements to the Gambia river as risk factors for livestock seropositivity (Jarra et al., 2025). These findings underscore the need to explore how livestock movement and ecological variability interact to shape RVF dynamics in The Gambia.

In a disease system where transmission is influenced by the movement of hosts, the incorporation of movement processes between habitats is essential for accurately modelling its dynamics (Merkle et al., 2017). Recent studies have integrated livestock movement patterns into the classical compartmental model framework to enhance our understanding of RVF dynamics. For instance, simulations have explored how nomadic ruminants, which occasionally or seasonally migrate from pond to pond according to the availability of water, sustain endemic RVFV circulation in the Ferlo region of Senegal (Durand et al., 2020). Similarly, Tennant et al. (2021) demonstrated that localised environmental conditions and inter-island livestock movements could sustain RVFV circulation in the Comoros Archipelago, without reintroduction from mainland Africa. Notably, this study used vegetation index as a proxy parameter for mosquito-borne transmission in a data-constraint context and achieved predictions consistent with empirical observations (Tennant et al., 2021). While these studies provide valuable insights, their focus on specific ecological and epidemiological contexts limits their applicability to The Gambia, which presents distinct ecological and livestock movement characteristics.

This study aimed to develop and apply an SIR model to simulate RVF dynamics in the cattle population of The Gambia, with a focus on incorporating seasonal transhumant movement as a key mechanism in explaining observed RVF epidemiology. We used household survey data to reconstruct network connections between cattle herds at water and grazing locations, and estimated seasonally dynamic transmission matrices for distinct ecoclimatic regions based on livestock movements. These transmission matrices were integrated into the model to estimate basic reproduction numbers and simulate temporal trends and infection magnitudes that approximate observed RVF dynamics. Using estimates of the force of infection in specific ecoclimatic regions, we predicted age-dependent seroprevalence, and estimated seropositivity waning rates consistent with observed data. Finally, we evaluated the potential impact of a transhumant movement ban on the modelled dynamics. By applying both deterministic and stochastic versions of the same model within an Approximate Bayesian framework, this study sought to provide insight into the intersection of livestock movement and ecological conditions that supports RVF dynamics in The Gambia.

Methods

RVFV seroprevalence and household survey data

A cross-sectional sero-epidemiological survey conducted in 2022 collected 3,602 serum samples from ruminant livestock (cattle, sheep and goats), linked to questionnaire data from 202 livestock-owning households in 52 villages (Jarra et al., 2025). Details of the study location and data collection procedures are provided in the Supplementary Text. Due to the very limited number of small ruminants (sheep and goats) involved in transhumant movement in our household survey data, we focus this study on the cattle population. The estimated RVFV-specific IgG seroprevalence in cattle was 36.8% (95% Confidence Interval [CI]: 31.6 – 40.0), with seroprevalence ranging from 0% to 69.8% (95% CI: 47.1 – 92.4) in cattle-owning villages (Jarra et al., 2025).

Epidemiological setting and cattle population structure

Following Sumaye et al. (2019), we define two regions with distinct ecological and climatic characteristics (hereafter referred to as ecoclimatic regions or eco-regions) in The Gambia: the Gambia river ecoclimatic region with its year-round suitability for RVFV transmission contrasts with the Sahelian ecoclimatic region, which experiences seasonal rainfall and hot arid dry seasons. Our data revealed that 30% of the cattle-owning households participate in transhumant movements from the Sahelian ecoclimatic region to the Gambia river during the preceding dry season. Resident cattle herds are assumed to remain at their homestead villages year-round. This observation led to the classification of cattle into three structured subpopulations; M: resident cattle in Sahelian villages, present year-round in this eco-region without access to the Gambia river. L: resident cattle in villages within the Gambia river ecoclimatic region, with year-round access to its floodplains. T: transhumant cattle in Sahelian villages which seasonally move to the Gambia river eco-region during the dry season (November - June) and return at the start of the wet season (July - October).

The observed RVFV IgG seroprevalence in each of the three populations was estimated as M: 30% (95% CI: 27 – 34), L: 36% (95% CI: 28 – 45) and T: 45% (95% CI: 41 – 48). A 2016 national livestock census of The Gambia estimated the total cattle population at 292,837 (Department of Livestock Services, 2016, pers. comms.). The livestock census includes complementary data on livestock management practices and grazing patterns. These data were used to disaggregate the total cattle population across the country into Sahelian resident (M), transhumant (T), and river resident (L) subpopulations, providing the basis for the subpopulation estimates in Table 1.

Preliminary parameter values for the SIR model.

Movement network analysis

The household herd was selected as the mapping unit to capture both daily local movements within the vicinity of the homestead villages and seasonal transhumant movements. Specifically, herds in the same village share common water and grazing locations and are also considered to share these locations with herds in neighbouring villages, creating homogenous mixing of cattle within neighbouring villages. The transhumant movements create epidemiological links between geographically distant herds that would otherwise be isolated. Therefore, to investigate how seasonal transhumant movements to risk areas with suitable eco-climatic conditions for RVFV-competent mosquito species might affect RVF dynamics in The Gambia, we first constructed an undirected network that represents the movement events within and across the three structured subpopulations.

The nodes of the network represented herds, and the edges connected to shared grazing and water locations. This allowed us to investigate important influencing movements in the network. At the node level, summary statistics were calculated for relevant normalised centrality measures such as degree, betweenness, and eigenvector centrality to quantify the importance of each node in the network (see Supplementary Table S1). The geographic distances of the edges were calculated to measure the total distance of travel between connected herds. This was performed by calculating the mean Euclidean distance between herds that shared resource areas. Finally, we used Spearman’s rank correlations to examine the relationship between network measures, the geographic distance of connected herds, and mean herd-specific RVFV seropositivity.

General overview of the deterministic model

To better understand and simulate the dynamics of RVFV transmission among the cattle population in The Gambia, we developed a deterministic susceptible-infectious-recovered (SIR) model that incorporated data from the 2022 survey and accounted for transhumant cattle movements across the two different ecoclimatic regions 𝑖 (𝑖 ∈ {Gambia river eco-region, Sahelian eco-region}) The model operated in continuous time with daily observations and incorporated matrices reflecting cattle movements and region-specific transmission parameters, under the assumption of random mixing among varyingly combined cattle subpopulations within each eco-region. The model was required to track the three cattle subpopulations (M, L, and T) as they fuse and fission over the wet and dry seasons in the different eco-regions. Figure 2 represents the three distinct processes in the model: seasonal cattle movements (arrival and departure of transhumant cattle herds), the population dynamics (birth and death of cattle), and the RVFV dynamics (infection and transmission) in each eco-region.

Spatial range of transhumant cattle movements in The Gambia identified from the household survey.

Movements originated from study villages (represented as red squares) and extended to destination villages, which are either other study villages or villages not selected for this study (represented as black triangles). This study was conducted during the dry season, and all transhumant herds were sampled at their destination villages. Distinct directional movements between homestead villages and their respective destinations are illustrated with blue arrows. The Gambia river is depicted as the meandering white line running the length of the country. The map was generated using the igraph package in R.

Epidemiological model of RVFV transmission and infection dynamics among the cattle subpopulations in The Gambia.

A. Schematic representation of ecoclimatic region and seasonal combinations that influences RVFV transmission between the Sahelian areas and the Gambia river eco-region. B. The transmission matrices that determined possible RVFV transmissions between the subpopulations during each season. C. Transmission framework of the within eco-region RVFV transmission, parameters are defined within the main text.

Distribution of degree, normalised betweenness and eigenvector centrality values based on cross-sectional household survey data on livestock movement.

Only a few of households showed high centrality values, highlighting their importance in the network.

The simulation assumed that RVFV was introduced as a single event into a fully susceptible population. Entry and exit of cattle through trade or gifting outside the system were not considered. Each structured cattle subpopulation 𝑗 (𝑗 ∈ {M, L, T}) was divided into number of susceptible (𝑆𝑗), infectious (𝐼𝑗), and recovered (𝑅𝑗) compartments. The processes driving population dynamics were per-capita birth (𝑏) and natural death (𝜇). Susceptible cattle became exposed to infection at a transmission rate 𝛽𝑖,𝑠𝑒𝑎𝑠𝑜𝑛 that depends on the eco-region (𝑖) and the season (Figure 2B). Infectious cattle recovered from RVF at rate 𝛿 or died of the disease at mortality rate 𝛾. The model accounted for temporal variations in population sizes due to transhumant movements, and explicitly tracks seasonal population changes in each region:

  • Sahelian eco-region (𝒔): During the wet season, both resident (M) and transhumant (T) cattle are combined (𝑁𝑠,𝑤𝑒𝑡(𝑡) = 𝑇 + 𝑀). Only resident cattle (M) remain during the dry season (𝑁𝑠,𝑑𝑟𝑦(𝑡) = 𝑀).

  • Gambia river eco-region (𝒓): During the wet season, only the resident cattle (L) are present (𝑁𝑟,𝑤𝑒𝑡(𝑡) = 𝐿), while transhumant cattle (T) join during the dry season (𝑁𝑟,𝑑𝑟𝑦(𝑡) = 𝐿 + 𝑇).

We modelled the RVFV transmission rate (𝛽𝑖,𝑠𝑒𝑎𝑠𝑜𝑛) to account for both mosquito-borne and direct transmission – contact with highly infectious abortion (calving) products. Cattle of all age groups were assumed to participate in the seasonal movement between the two eco-regions, a practice consistently observed among transhumant households over the years. Only movement of the transhumant subpopulation (T), serving as the primary epidemiological link between the resident cattle in the Sahelian (M) and those in the Gambia river eco-region (L), was considered. This necessitated changing the RVFV transmission matrices (Figure 2B) between wet and dry seasons to accurately reflect the effects of seasonal movement on RVF transmission dynamics.

Estimation of model parameters

We estimated basic reproduction numbers (𝑅0,𝑖) – the number of secondary infections caused by one infectious cattle in an entirely susceptible population in each ecoclimatic region – based on a standard expression adapted from Anderson and May (1991). This estimation was based on the proportion of the total cattle population in each eco-region that is susceptible (𝑠𝑖), assuming homogenous mixing of the cattle subpopulations. Preliminary estimates of seasonal transmission parameters 𝛽𝑟,𝑤𝑒𝑡 and 𝛽𝑟,𝑑𝑟𝑦 were then calculated as the remaining unknown that can be solved for from a standard expression for 𝑅0,𝑖 derived from the SIR model (see Supplementary Text). To align the wet season transmission in the Sahelian eco-region (𝛽𝑠,𝑤𝑒𝑡) with potential year-round transmission observed in the Gambia river eco-region, a scaling factor 𝜓 was applied, thus our preliminary estimate of 𝛽𝑠,𝑤𝑒𝑡 = βs,dry is considered to be negligible during the longer dry season (Chevalier et al., 2004), and thus set to zero.

Preliminary estimates of four key parameters (𝑏, 𝜇, 𝛿 and 𝛾) were obtained from existing literature (Table 1). The scaling factor 𝜓 was initially set to 2, due to the absence of an established value in the literature. A final re-estimation of all eight parameters was conducted by fitting the SIR model to the observed seroprevalence using a simulation-based Approximate Bayesian Computation (ABC) framework. Informative priors for each parameter were defined as triangular distributions, with modes based on the preliminary estimates with upper and lower bounds set to ±70% of these estimates. The rejection method of the ABC algorithm was implemented in the abc package in R [version 2.2.1] (Csillery et al., 2012). The priors were used to generate 500,000 randomly parameterised models, and a snapshot RVFV seroprevalence was predicted for each model (summary statistics) at the end of each simulation. These predicted seroprevalence values were then compared to the observed RVFV seroprevalence for each cattle subpopulation. Sets of proposed parameters were accepted for the 0.2% best fitting models (as determined by the Euclidean distance between predicted and observed seroprevalence values for each subpopulation). The final parameter values were taken as the means of the posterior distributions and used for subsequent model analysis and simulations.

The RVFV dynamics were mathematically formalised using ordinary differential equations (ODEs). A full description of the equations is provided in the Supplementary Text. An epidemic was initiated at t = 0 by introducing one infectious cattle into each previously RVFV-free subpopulation. The simulation was run for 20 years (representing the approximate period from the initial RVFV outbreak in The Gambia in 2002 to our serological survey in 2022). The model incorporated distinct seasonal values of the transmission parameters (𝛽𝑖,𝑠𝑒𝑎𝑠𝑜𝑛) (Figure 2B). The dry season spanned 245 days (days 1 – 245), while the wet season lasted 120 days (days 246 – 365) of each year. The model iteratively checked the time variable daily to apply the appropriate seasonal parameters. The ODEs are implemented and solved in R statistical software [v4.2.2] (R Core Team, 2022) using the package deSolve (Soetaert et al., 2010). All plots are generated using the ggplot2 (Wickham, 2009).

Estimating the per capita rate of decay of RVFV seropositivity

The force of infection [FOI] (𝜆𝑖,𝑠𝑒𝑎𝑠𝑜𝑛), representing the per capita rate at which susceptible cattle become infected with RVFV in each ecoclimatic region during wet and dry seasons, was estimated and a modified version of the model (where age of cattle was substituted in for time) was used to predict age-seroprevalence profiles for each subpopulation in each eco-region on the condition of quasi-equilibrium (see Supplementary Text). Given that cattle typically remain infectious with RVFV for about one week, the variation in 𝜆𝑖,𝑠𝑒𝑎𝑠𝑜𝑛 was assessed at weekly intervals over a 10-year period, which corresponds to the assumed lifespan of cattle in this study.

Cattle that recover from RVFV infection were assigned to an immune seropositive compartment (𝑃𝑗), and the age-seroprevalence curves predicted were compared with observed data. To adjust for inconsistencies between the predicted and observed curves, the predicted age-seroprevalence in each subpopulation was fitted to the observed seroprevalence by including a decay parameter (𝜋) which described the rate of waning of RVFV seropositivity in recovered cattle over time (Figure 2C). The proportion of immune cattle in which seropositivity decayed over time transitioned to an immune seronegative compartment (𝐷𝑗) (see Supplementary Text).

The age-seroprevalence model was initialised with a fully susceptible population of newly born calves, aging over a 10-year lifespan, and subjected to seasonal variations in the 𝜆𝑖,𝑠𝑒𝑎𝑠𝑜𝑛 observed at conditions of quasi-equilibrium. A preliminary decay rate of 0.005 per immune cattle per week was set, and a final value was estimated by fitting the predicted seroprevalence to the observed data using ABC. Triangular priors were assigned to 𝜋, with limits extending ±70% from the preliminary rate. The posterior mean of 𝜋 was used to calculate the half-life of RVFV seropositivity decay in cattle in The Gambia.

Estimation of basic and seasonal reproduction numbers for the entire system

The basic reproduction number (𝑅0) for the fully susceptible cattle population in the entire system during the wet and dry season was calculated by numerically solving for the leading eigenvalues of a matrix Q, using the next-generation matrix method described by Diekmann et al. (2010) (Supplementary Text). Additionally, to account for dynamic susceptibility and transmissibility within the cattle population, a seasonal reproduction number (𝑅𝑠𝑡) was computed at weekly intervals over the course of the 10-year simulation. At each weekly time step, the model output provided the number of cattle in each partially susceptible subpopulation, which is used to determine 𝑅𝑠𝑡.

Stochastic modelling and RVFV extinction risk

To complement the deterministic solution, a parallel stochastic simulation was conducted to capture demographic stochasticity arising from the inherent probabilistic dynamics of RVFV transmission. The model employed the tau-leap method (Gillespie, 2007), in combination with the Gillespie algorithm to improve computational efficiency. Within our model, the expected rate for each event type (e.g., population dynamics, transmission dynamics, and infection dynamics) was calculated based on the current state of the system, and the number of events occurring within a time step (𝜏 = 0.1 simulated day) was then assumed to be Poisson-distributed. The states of all subpopulations were iteratively updated according to the number of sampled events.

We approximated the extinction risk of RVFV in two contexts: (i) within the transhumant subpopulation (T) returning to the Sahelian eco-region at the start of each wet season (stochastic fade-out/local extinction), and (ii) across all three subpopulations (system-wide extinction). For this purpose, 1,000 independent stochastic realisations of the model were conducted, to ensure statistical robustness. RVFV extinction risks were quantified by recording the frequencies at which the number of infectious cattle declined to zero within the transhumant subpopulation (Supplementary Text) and across the entire system. It is important to note that the impact of environmental stochasticity was not considered in this analysis.

The impact of a transhumant movement ban on RVF dynamics

To evaluate how cattle movement to the Gambia river eco-region during the dry season influences RVF dynamics, we simulated a ban on transhumant movement. In this scenario, we assumed the subpopulation configuration is fixed in the wet season structure for each eco-climatic region for both the wet and dry season. The impact of this ban was assessed over the 20-year simulation period and compared to baseline RVF dynamics where cattle movement is unrestricted.

Elasticity Analyses

We evaluated the sensitivity of the predicted RVFV seroprevalence in each of the three subpopulations to variations in the eight posterior parameter values: 𝑏, 𝜇, 𝛾, 𝛿, 𝜓, 𝛽𝑠,𝑤𝑒𝑡, 𝛽𝑟,𝑤𝑒𝑡, 𝛽𝑟,𝑑𝑟𝑦. Each parameter was independently sampled 1,000 times from a uniform distribution, constrained within ±20% of its mean posterior value, while the remaining seven parameters are held at their respective mean posterior values. For each modified parameter set, a new predicted seroprevalence was estimated for each subpopulation. The percentage change in predicted seroprevalence was calculated as the difference between the predicted seroprevalence under the original posterior means and the predicted seroprevalence obtained with the modified parameter value. Likewise, the percentage change in parameter value was calculated as the deviation of the randomly sampled parameter value from its mean posterior value. The elasticity was defined as the percentage change in the seroprevalence (the response variable) per unit percentage change in the parameter value (the explanatory variable) and estimated using a general linear model (GLM). Additionally, a Loess smoothing function was applied to evaluate the adequacy of the linear approximation and to identify potential nonlinear trends.

Results

Household herd connectivity

The network analysis revealed variations in between-herd connectivity through shared water and grazing locations (Supplementary Figure S2). Centrality measures were disproportionately influenced by a small subset of transhumant herds with relatively high centrality values, highlighting their prominent role within the network (Figure 3). However, these connections show no significant correlation with herd-specific RVFV seropositivity (Supplementary Figure S4).

Time-series of the proportion of infectious cattle in the three subpopulations across the two eco-regions.

Main: The simulated RVF infection dynamics from the deterministic model (solid lines) together with twenty realisations of the stochastic model (dashed lines). Insert: A magnified view highlighting the seasonal dynamics of RVFV transmission focused on the last two years of the 20-year deterministic simulation in the cattle subpopulations, highlighting finer seasonal variations (wet season = beige; dry season = cyan). The proportion of infectious cattle peaked at the latter part of the wet season, but infections quickly disappeared in the dry season in the Sahelian eco-region. Once transhumant herds begin to arrive in the Gambia river region, infections are predicted to rise. I = infectious cattle in each subpopulation. A square root transformation was applied to the y-axis for visualisation purposes using coord_trans(y = "sqrt"), while the data remained in its original scale.

Transhumant movements during the dry season covered considerable distances, with herds directed towards the Gambia river from within the Sahelian eco-region. The pairwise geographic distances between connected herds from different study villages ranged from 3.9 to 81.4 km (median 13.5 km). A weak negative correlation (Spearman’s R = -0.073, p < 0.05) was observed between the geographic distance of connected herds and their mean RVFV seroprevalence (Supplementary Figure S4), suggesting that herds in closer proximity share common risk factors contributing to similar RVFV seropositivity.

Estimated model parameters

The posterior distributions of the eight estimated parameters, derived through the Approximate Bayesian framework, are shown in Supplementary Figure S5. The mean posterior values and their 95% credible intervals (CrI) are summarised in Table 2. These posterior means are used in all subsequent simulations.

Mean posterior values and 95% CrI of the estimated parameters.

Seasonal and regional RVF dynamics

We simulated the timing and magnitude of RVFV infections across three structured cattle subpopulations, which seasonally occupy two eco-regions with distinct ecological characteristics. At the end of the 20-year simulations, the model remains in a transient phase without fully attaining a quasi-equilibrium state (Figure 3). The results exhibited an oscillatory pattern, with recurrent spikes in the proportion of infectious cattle aligning with the wet season. In the Sahelian eco-region, infections peaked at approximately 5% of the cattle population between the 64th and 111th days of the wet season (September – October) (Figure 3 insert), followed by a secondary phase characterised by sharp decline of infection during the dry season, accompanied by gradual decrease in herd immunity.

While both ecoclimatic regions showed similar temporal RVFV infection dynamics during the wet season, simulations indicated sustained low-level transmission in the Gambia river eco-region throughout the dry season (Figure 3 insert). This continued virus activity suggests a potential “source-sink” effect, as a small fraction (approximately 1%) of the transhumant subpopulation (T) returned to the Sahelian eco-region at the start of the wet season while infectious.

At the end of the deterministic simulation, the maximum number of infectious cattle in the M, L, and T subpopulations during the wet season were approximately 4,562, 2,467, and 443, respectively. During the dry season, the minimum number of infectious cattle remaining in the L and T subpopulations within the Gambia river eco-region were approximately 37 and 4, respectively. By simulating seasonal transmission, and incorporating cattle movement matrices, our deterministic model provided a close approximation to the RVFV infection fluctuations demonstrated by the stochastic model (Figure 3).

The stochastic simulation indicated that RVFV infection undergoes periodic demographic fade-outs (local extinctions) in the transhumant (T) subpopulation, occurring in 73.8% of realisations (738 of 1,000) during the dry season in the Gambia river eco-region. Each of these realisations predicted between one and five distinct local extinction events (Supplementary Figure S6). However, RVFV infection persisted within the broader river eco-region through transmission in the river resident (L) subpopulation, facilitating recurrent re-infection of the T subpopulation in subsequent years. Overall, the stochastic simulation recorded 1,363 local extinction events with an average transmission duration of 13.9 years, which approximated to 9.8% local extinction probability in the T subpopulation. Accordingly, the annual re-introduction of RVFV into the Sahelian eco-region at the start of the wet season, via the T subpopulation returning from the Gambia river eco-region, was predicted to occur at a statistically high rate. Finally, in 40.2% of realisations (402 of 1,000) RVFV underwent complete extinction across the entire system before the conclusion of the 20-year simulation.

Our results showed that, in most realisations, the transhumant subpopulation experienced local extinctions within the first 1,000 days, coinciding with the initial transient phase of the RVF dynamics – big outbreaks within a fully susceptible subpopulation. As the system gradually converged on a quasi-equilibrium state, the probability of RVFV persistence in the T subpopulation stabilised (Figure 4A). In contrast, RVFV exhibited a higher initial persistence probability within the entire system, with most system-wide extinctions occurring beyond the 5,000th day of the simulation as infection dynamics neared a quasi-equilibrium state (Figure 4B).

Extinction rate of RVFV over time (red line with 95% CI – grey ribbon), in the transhumant subpopulation (A) and the entire cattle population within the system (B) based on 1,000 stochastic realisations, illustrating differences in the timing of local and system-wide extinctions, respectively.

Most of the local extinction occurred shortly after RVFV introduction into a fully susceptible population in the T subpopulation. Note: the timing of the local extinction in the T subpopulation depicted here was the first extinction event within this subpopulation, and re-infection occurs whence the subpopulation returned to the river.

Estimated 𝑹𝒔𝒕 and 𝝀𝒊,𝒔𝒆𝒂𝒔𝒐𝒏

The 𝑅0 computed for the entire system assuming a fully susceptible cattle population, was 3.15 during the wet season and 1.33 during the dry season. The seasonality of transmission was further illustrated by variations in the 𝑅𝑠𝑡 value (Figure 5), which revealed peaks (𝑅𝑠𝑡 > 1) during the wet season across the entire system, aligning with the amplification of infection as epidemiological conditions are assumed to be favourable for transmission. In the latter part of the dry season, 𝑅𝑠𝑡 > 1 values were also observed as the susceptible population recovered, supporting the localised persistence of transmission in the Gambia river eco-region (Figure 5). Table 3 summarises 𝑅𝑠𝑡 and 𝜆𝑖,𝑠𝑒𝑎𝑠𝑜𝑛 values across both eco-regions and seasons, providing insights into the RVFV dynamics across temporal and spatial scales.

The full 10-year simulation of the weekly 𝑅𝑠𝑡 (green), and force of infection, 𝜆𝑟𝑒𝑔𝑖𝑜𝑛,𝑠𝑒𝑎𝑠𝑜𝑛, in each eco-region (blue = Sahelian, red = Gambia river eco-region).

Shaded areas correspond to the seasons (wet season = beige; dry season = cyan) in the last two years of the simulation. A square root transformation was applied to the y-axis for visualization purposes.

Summary of estimated seasonal reproduction number (𝑅𝑠𝑡) and regio-specific force of infection (𝜆𝑖)

Predicted age-seroprevalence curve

Here, we fitted a model using the weekly 𝜆𝑖,𝑠𝑒𝑎𝑠𝑜𝑛 values estimated for each ecoclimatic region at conditions of quasi-equilibrium to explore the observed trend in age-seroprevalence. Age was used as a proxy for cumulative exposure duration, consistent with prior methodologies (Salje et al., 2016). A cohort of newly born calves was subjected to the specific 𝜆𝑖,𝑠𝑒𝑎𝑠𝑜𝑛 associated with the eco-region they occupy across their lifespan to estimate the cumulative exposure to RVFV in the absence of waning seropositivity. In the absence of any waning of seropositivity, the model substantially overestimated seroprevalence, predicting values of approximately M: 63%, L: 64% and T: 84% after 10-years (Figure 6A). These predicted seroprevalence values were then re-estimated accounting for potential waning of seropositivity in recovered cattle, enabling reconstruction of the age-stratified seroprevalence profile observed in The Gambia (Figure 6B). The posterior distribution of the rate of decay of RVFV seropositivity (𝜋) yielded a mean value of 0.00217 per week (95% CrI: 0.00216 – 0.00219). This decay rate corresponded to a seropositivity half-life of approximately 319 weeks (∼6 years). These results suggested that the observed seroprevalence was influenced by gradual loss of seropositivity over time in immune cattle.

Predicted RVFV seroprevalence in each population subject to FOI during the wet and dry seasons.

A: Cumulative growth of RVFV immunity over a 10-year period, representing the hypothetical lifespan of cattle, with no consideration for decay of RVFV seropositivity (𝜋 = 0). B: Proportion of immune cattle after introducing a seropositivity decay parameter (𝜋 = 0.0035), aligning the predicted seroprevalence with observed data. Observed seroprevalence across different age classes in the three structured cattle populations in The Gambia is shown as dots, with 95% confidence intervals (CI) as error bars. P: Cattle exposed to RVFV infection and recovered, assumed to be seropositive.

Impact of a transhumant movement ban on RVF dynamics

Under the ban, periodic large-scale RVF outbreaks during the wet season were predicted in the Sahelian eco-region, separated by inter-epidemic periods (IEPs) of approximately two to three years (Figure 7). While outbreaks became less frequent, the decline in herd immunity during these IEPs significantly increased the size of outbreaks, with infection peaks in both the Sahelian resident (M) and transhumant (T) cattle subpopulations reaching approximately 50%. However, RVFV seroprevalence in the M and T subpopulations in the Sahelian eco-region declined by 30% compared to the scenario allowing seasonal transhumant movement in the deterministic solutions. In contrast, the Gambia river eco-region exhibited a sustained transmission pattern similar to the dynamics observed when movement was permitted, though the proportion of infectious cattle increased by approximately 4% during the wet season.

Simulated RVF infection dynamic over a 20-year period within each structured cattle subpopulation, assuming a ban on seasonal transhumant movement.

In the Sahelian eco-region, periodic large outbreaks were separated by long inter-epidemic periods, allowing the susceptible subpopulations to regrow. In contrast, RVFV infection remained sustained and low in amplitude in the Gambia river eco-region – comparable to when movement was allowed. I = infectious cattle in each subpopulation. A square root transformation was applied to the y-axis for visualization purposes.

At the end of the 20-year simulation using the deterministic model, the maximum number of infectious cattle during an outbreak in the wet season was approximately 137,762 and 23,135 in the M and T subpopulations in the Sahelian eco-region, respectively. In the Gambia river eco-region, the wet-season infection count among the L subpopulation was approximately 3,883. These findings suggested the transhumant cattle movements played a critical role in modulating RVFV transmission dynamics, maintaining herd immunity, and mitigating large-scale outbreaks across eco-regions.

Elasticity Analysis

A 1% change in the mean posterior value of each of the eight posterior parameters resulted in a corresponding percentage change in predicted seroprevalence of each subpopulation (Table 4). The elasticity analysis indicated that the per-capita birth (𝑏) and natural death (𝜇) rates exerted the strongest influence on predicted seroprevalence (Table 4). In contrast, parameters related to transmission (𝛽𝑠,𝑤𝑒𝑡, 𝛽𝑟,𝑤𝑒𝑡, 𝛽𝑟,𝑑𝑟𝑦) and infection dynamics (𝜏, 𝛾 and 𝛿) exhibited relatively lower elasticities. Additionally, comparisons between GLM and Loess fits suggested that the relationship between the percentage change in the mean posterior parameter values and the corresponding percentage change in predicted seroprevalence was well-approximated as linear (Supplementary Figure S7).

The elasticity for all eight parameters in the GLM.

The parameters with the largest coefficient are in shown in bold.

Discussion

Previous modelling studies have significantly advanced our understanding of RVF epidemiology, particularly regarding the disease dynamics in ruminant hosts and mosquito vectors. These models have informed the development of predictive tools, surveillance frameworks, control interventions, as well as the evaluation of their impact (Anyamba et al., 2006; Gachohi et al., 2016; Gaff et al., 2007; Tennant et al., 2021). Despite these advancements, the persistence of RVFV in the Sahelian zone, where unfavourable dry season conditions make permanent RVFV circulation unlikely, remains poorly understood (Durand et al., 2020; Lumley et al., 2017). Increasing evidence suggests that livestock movements play a critical role in influencing RVF epidemiology in this zone (Belkhiria et al., 2019; Jarra et al., 2025). In this context, risk exposure to RVFV is primarily driven by movements into regions with high environmental suitability for RVFV-competent vectors in search of water and pasture, rather than transmission through direct livestock-to-livestock contact.

To investigate this hypothesis, we developed deterministic and stochastic models, parameterised using seroprevalence data collected in 2022 in The Gambia, to stimulate RVF dynamics. Our study focuses on the transhumant movement of cattle from the Sahelian ecoclimatic region of the country to the floodplains of the Gambia river as a potential mechanism accounting for the observed RVF dynamics among three cattle subpopulations. Key findings from our study include: (i) seasonal infection patterns peaking at approximately 3% and 1% of cattle in the wet and dry seasons, respectively, in the Gambia river eco-region, whereas infections peaked at approximately 5% in the wet season and disappeared in the dry season in the Sahelian region; (ii) local extinction and re-infection dynamics in the transhumant (T) subpopulation, with a high risk (90.2%) of re-introducing RVFV into their homestead villages in the Sahelian eco-region predicted by our stochastic model; (iii) the 𝑅0for the entire system was estimated as 3.15 during the wet season and 1.33 during the dry season, while the mean seasonal reproduction number (𝑅𝑠𝑡) in both seasons was >1; (iv) the age-seroprevalence data in combination with the predicted FOI in each eco-region implied a gradual decay in seropositivity in recovered cattle over time; and (v) a hypothetical movement ban had minimal impact on reducing infections over the long-term, as it led to less frequent but significantly more severe outbreaks in the Sahelian eco-region.

RVFV inter-epidemic transmission dynamics

Our modelling approach incorporates transhumant cattle movement, specific eco-regional transmission rates and ecological assumptions. Based on predictions that mosquito-borne transmission is responsible for 84% of RVFV infections in cattle (Nicolas et al., 2014), we pre-suppose the existence of an active mosquito population in the Gambia river eco-region (Snow, 1983) and show that RVFV infections can persist year-round therein. Despite the relatively low proportion of infectious cattle across the three structured subpopulations, RVFV transmission could be sustained in an endemic state over an extended period. Our simulations predicted there is 90.2% chance the transhumant subpopulation returning to their homestead villages at the start of the wet included a small proportion of infectious cattle (approximately 1% of the subpopulation). At the start and end of transhumance, GPS-tracked cattle herds in Cameroon exhibit active movement behaviour at an average daylight speed of 3 – 4 km/hour (Motta et al., 2018), a pattern corroborated by survey data from transhumant households in The Gambia. Given that RVFV incubation and viraemia last approximately seven to eight days, virus re-introduction to the Sahelian eco-region through the returning transhumant herds is plausible, posing a recurrent transmission risk. However, herd immunity among the Sahelian resident cattle population in our simulation prevents annual large-scale outbreaks.

The possibility of RVFV re-introduction between the two eco-regions is consistent with explanations given for the dissemination of the virus during outbreaks within and across countries through movement of potentially viraemic but asymptomatic ruminant livestock (Abdo-Salem et al., 2011; Shoemaker et al., 2002). In addition, similar patterns have been observed in northern Senegal (Cecilia et al., 2020; Durand et al., 2020) and Tanzania (Sumaye et al., 2013), where movement of infectious livestock between regions with favourable ecological conditions is thought to maintain RVFV endemicity. Our simulation results of the timing and magnitude of infections are supported by evidence showing RVFV amplification during the wet season, with peak circulation occurring at the latter part of the season in Senegal (Durand et al., 2020). The 𝑅𝑠𝑡 reflected this, with mean 𝑅𝑠𝑡 values greater than one in the first half of the wet season (9 of 17 weeks) in the entire system, and in over two thirds of the dry season (23 of 35 weeks) in the Gambia river eco-region. Earlier studies conducted during an inter-epidemic period by Metras et al., (2017) (maximum 𝑅𝑠𝑡 = 2.19) and Tennant et al., (2021) (maximum 𝑅𝑠𝑡 = 2.77) in Mayotte estimated seasonal reproduction numbers broadly consistent with those estimated by our study.

Age-seroprevalence curve and seropositivity decay

Acquired immunity following natural infection with many communicable diseases is often imperfect, with immunity potentially waning over time, offering partial protection or failing entirely (Le et al., 2021). Although acquired immunity to RVFV is often considered long-lasting (Swanepoel & Coetzer., 2004; Wilson, 1994), there is lack of empirical data regarding the nature and longevity of IgG seropositivity following natural infection in ruminant livestock. In a long-term cohort study, Wright et al. (2020) observed that RVFV-specific IgG and neutralising antibodies remained detectable in adult humans up to eleven years post-infection, although antibody titers progressively declined. This finding has important implications for the interpretation of serological data: while antibodies may be present, the ability to detect them (i.e. seropositivity) may decline over time. Serological diagnostics reliant on detecting specific antibody thresholds may classify immune individuals as seronegative, potentially underestimating the true extent of population immunity. This implication has been observed in other viral infections. For instance, in a 13-year longitudinal study of 34 healthcare workers infected with SARS-CoV, the 50th percentile exhibited anti-nucleocapsid IgG titers below the assay detection limit by year nine, and 75th percentile fell below the limit by year thirteen (Guo et al., 2020). These observations reinforce the importance of considering seropositivity decay when interpreting long-term antibody persistence.

Our sero-epidemiological study in The Gambia revealed an age-seroprevalence pattern in cattle that is characterised by high IgG seropositivity during early calfhood, that quickly plateaued or declined with advancing age (Jarra et al., 2025). To reconcile the discrepancies between this observed age-seroprevalence pattern and the predictions based on our FOI estimates from each eco-region, we incorporated a seropositivity decay parameter into the model. Specifically, RVFV-exposed cattle in which seropositivity decayed over time transitioned into an immune but seronegative compartment. This adjustment is supported by evidence that cellular immunity (T cell-mediated responses) can confer durable protection even in the presence of substantially low RVFV-specific antibody levels following both natural RVFV infection in humans (Wright et al., 2020) and wildtype RVFV challenge in vaccinated mice (Doyle et al., 2022).

Although earlier RVF seroprevalence-based modelling studies have not explicitly incorporated or contextualized the kinetics of IgG decay, Le et al. (2021) provides a valuable theoretical framework of imperfect infection-derived immunity by outlining seven distinct models. These models span between the two classical extremes: the "perfect" immunity assumed in SEIR frameworks (recovered individuals are permanently protected) and the “none at all” immunity characteristic of SEIS models (individuals become susceptible again after recovery). Each of these immune decay models produced unique implications for the overall disease dynamics within a population (Le et al., 2021). In the case of RVF, the lack of longitudinal data on the duration and durability of natural antibody responses in ruminant livestock represents a potential limitation that hinders the accurate parameterisation of models and prediction of transmission dynamics, particularly in systems where seropositivity may wane or become undetectable over time despite persistent protection.

Impact of transhumant movement on RVF dynamics

The inclusion of transhumant cattle movement and seasonality in our transmission matrices allowed us to assess its impact on RVF dynamics in The Gambia. Our results suggested that a movement ban does not lead to marked improvement in long-term epidemiological outcomes. Specifically, it only serves to reduce the frequency of outbreaks during the wet season, allowing gradual reduction of herd immunity due to the prolonged IEPs until a critical threshold after which a more severe outbreak occurred in the Sahelian eco-region. In this movement ban scenario, we observed that the estimated proportion of infectious cattle in the Sahelian eco-region reached 50% compared with approximately 5% when movement was permitted. The severe pattern of RVFV infection dynamics has not been reported in the past two decades in The Gambia. These findings are consistent with the role of transhumant cattle movements in possibly maintaining herd immunity and mitigating the magnitude of infections, to align with the observed RVF dynamics in The Gambia.

Despite the absence of susceptible transhumant cattle moving to the river, RVFV transmission was predicted to continue year-round in the Gambia river eco-region. Furthermore, results from our network analyses suggest that contact between household herds did not appear to influence RVFV seroprevalence (see Supplementary figures). These results supported evidence from the Comoros archipelago, indicating that interventions focusing solely on movement restrictions were insufficient in effectively controlling RVFV transmission once the virus had been introduced (Tennant et al., 2021). The limited impact of the movement ban in reducing the number of infections in our study highlighted the complex nature of RVF transmission and emphasised the need for more comprehensive preventive strategies beyond movement restrictions.

Implications of the model and future enhancements

While this study provided valuable insight into the dynamics of RVF transmission in The Gambia, several data constraints should be considered. The household survey data, which served as the foundation for our modelling approach, is limited by gaps in spatial data coverage and temporal resolution, and may not fully represent the diversity of livestock management systems across the country. Consequently, several assumptions were made, introducing uncertainty in model predictions, particularly regarding the exact timing and spatial extent of transhumant movements. Future studies would benefit from more comprehensive data collection efforts, especially focusing on finer-resolution data regarding livestock movements across different regions.

Conclusion

This study contributes to our understanding RVFV epidemiology in The Gambia by integrating transhumant cattle movement into a modelling framework that accounted for eco-regional transmission dynamics. Our findings highlighted seasonal transmission patterns, with year-round infections predicted in the Gambia river eco-region reinforcing the role of ecological suitability in sustaining RVFV. While transhumance played a role in shaping RVFV transmission through a high risk of RVFV re-introduction, herd immunity limited large-scale outbreaks in the Sahelian eco-region. The predicted seropositivity decay underscored the importance of considering waning of antibody titers in serological assessments. Future refinements incorporating high-resolution livestock movement data and environmental variability will enhance predictive accuracy. There is need for strategies informed by transhumance practices in managing RVF, particularly in remote livestock production communities in The Gambia. Understanding livestock owners’ perceptions of seasonal movements, and the determinants of these movements could inform the design of socially acceptable measures such as targeted vaccination and the strategic or pre-emptive application of topical acaricides for ruminants during peak vector activity in the wet season or before transhumant movements into the river floodplain.

Data Availability Statement

The codes that support the findings and produce the figures of this manuscript are available on Github: https://github.com/enz-j/Rvfv_dynamics_gambia

Acknowledgements

We thank the communities and livestock-owning households for participating in the study. We are also deeply grateful to the livestock assistants of the Department of Livestock Services in The Gambia for their contributions to the field work.

Additional information

Ethical clearance

The study protocol received ethical approval from the National Agricultural Research Institute of The Gambia [Approval Reference: NARI/DLS/01/(01)] and the University of Glasgow School of Veterinary Medicine Research Ethics Committee (Reference: EA47/21). Approval for interactions with human participant were obtained from The Gambia Government/MRCG@LSHTM Joint Ethics Committee (Project ID/Review Ref: 26750) and the University of Glasgow College of Medical, Veterinary and Life Sciences Ethics Committee (Project No: 200210057).

We adhere to the STROBE Statement–Checklist of items that should be included in reports of cross-sectional studies.

Consent

Community engagement meetings were held to introduce the study, explain sampling procedures, and obtain informed consent. Participants provided written consent or a thumbprint as a substitute in cases of illiteracy. Consent procedures emphasized voluntary participation, anonymity, and participants’ right to withdraw. Collaboration with the local livestock assistants of the Department of Livestock Services (DLS), who are familiar with and trusted by the communities, facilitated recruitment of participants. All DLS livestock assistants received biosafety and ethics training prior to data collection. No human samples were collected for this study.

Author Contributions

E.J. and D.T.H. conceptualized the study; E.J. acquired the funding; E.J. coordinated and administered the field work; E.J. and D.T.H developed the deterministic and stochastic models; D.E. conducted the network analysis; E.J. wrote the original draft with input from D.T.H. and D.E.; D.T.H. and S.C. supervised the project; All authors reviewed, edited and approved the manuscript.

Funding

This work was funded by the Foreign, Commonwealth & Development Office of the UK Government through the Commonwealth Scholarship Commission for a PhD studentship for Essa Jarra at University of Glasgow. The Commonwealth Scholarship Commission was not involved in study design, data collection and interpretation, or the decision to submit the work for publication. The opinions expressed in this publication are entirely those of the authors.

Additional files

Supplementary File 1