pith. sign in

arxiv: 2511.16802 · v2 · submitted 2025-11-20 · 🧬 q-bio.PE · math.DS

A model for mosquito-borne epidemic outbreaks with information-dependent protective behaviour

Simone De Reggi , Andrea Pugliese , Mattia Sensi , Cinzia Soresina This is my paper

Pith reviewed 2026-05-17 20:04 UTC · model grok-4.3

classification 🧬 q-bio.PE math.DS
keywords mosquito-borne epidemicsprotective behaviourinformation-dependent responsebasic reproduction numbergeometric singular perturbationrecurrent epidemic waveshost compositionvector transmission
0
0 comments X

The pith

Protective behaviour based on disease information can raise or lower the reproduction number of mosquito-borne epidemics and trigger repeated waves.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper builds a model of mosquito-borne epidemics in which humans adopt protective actions against bites when they receive information about past or current prevalence. Mosquitoes are allowed to feed on both disease-competent and non-competent hosts. The authors first show that the resulting behavioural response can either decrease or increase the basic reproduction number, depending on the strength of the response, the proportion of competent hosts, and the transmission parameters. They then assume that opinion formation occurs on a much faster timescale than infection dynamics and apply a reduction technique to obtain a simpler model of a uniform host population. Analysis of the reduced system reveals that the same behavioural mechanism can produce low-attack-rate outbreaks, speed epidemic fade-out, or instead sustain the disease through recurrent damped waves.

Core claim

The authors establish that behaviour-driven protection may either decrease or increase the basic reproduction number depending on the interaction between behavioural response, host composition, and transmission parameters. Under the assumption that opinion dynamics evolves on a much faster time scale than disease transmission, Geometric Singular Perturbation Theory reduces the original two-group model to one for a homogeneous host population; the reduced system shows that information-induced behavioural changes can facilitate epidemic control or prolong disease persistence, potentially generating recurrent damped epidemic waves and outbreaks with low attack rates.

What carries the argument

Reduction of the two-group model to a homogeneous-host model via Geometric Singular Perturbation Theory applied to fast opinion dynamics.

If this is right

  • Behaviour-driven protection may decrease or increase the basic reproduction number.
  • Protective measures can facilitate epidemic control under suitable host composition and parameters.
  • Behavioural responses may prolong disease persistence and generate recurrent damped epidemic waves.
  • Information-induced changes can produce outbreaks that reach only low attack rates.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Public health messaging that supplies timely prevalence data might be tuned to favour the control-promoting regime rather than the persistence-promoting one.
  • The same fast-behaviour reduction could be tested on other vector-borne diseases that involve mixed host competence.
  • Collecting paired data on behavioural response speed and epidemic wave patterns would directly test whether the timescale separation holds in real populations.

Load-bearing premise

Opinion dynamics evolves on a much faster time scale than disease transmission.

What would settle it

Field or simulation data showing no recurrent damped waves and no low-attack-rate outbreaks despite strong information-dependent protective behaviour would contradict the reduced model's transient predictions.

Figures

Figures reproduced from arXiv: 2511.16802 by Andrea Pugliese, Cinzia Soresina, Mattia Sensi, Simone De Reggi.

Figure 1
Figure 1. Figure 1: The split of the human population in protected and non-protected individuals. [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Flow chart for system (5). Straight lines: compartmental movements within each population; dashed lines: infections between populations (mosquitoes infecting humans and viceversa). In the following, we argue in terms of the fractions SP /H, IP /H, SNP /H, INP /H, and IM/M, which, for later convenience, will be still denoted with the same variables that have been using so far. The resulting model reads  … view at source ↗
Figure 3
Figure 3. Figure 3: The behaviour of the function F(p, q) in (15). Note that this is independent of all the parameters of the model but p and q. As for the second term on the LHS of (16), it is clear that p 1 − p + q 2p = p 1 − p(1 − q 2) < 1, for p ∈ (0, 1] and q ∈ [0, 1). Hence l hp 1 − p + q 2p − 1 i < 0, p ∈ (0, 1], q ∈ [0, 1), (18) which implies that the presence of other hosts l always contributes to reducing the value … view at source ↗
Figure 4
Figure 4. Figure 4: Rc as a function of q for different values of p and l, with ρ = 2 and epidemiological parameters as in [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: A simulation for model (5) with l = 0.25, ρ = 2, epidemiological parameters as in [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: A simulation for model (5) with l = 0.25, ρ = 2, epidemiological parameters as in [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: A simulation for model (5) with l = 1; ρ = 5.12, epidemiological parameters as in [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Flow chart for model (22). Then, taking into account (23), we modify (5) as follows:    S ′ P = −βH←MρIM qSP c(p, q) + l + a(J)SNP − w(J)SP , S ′ NP = −βH←MρIM SNP c(p, q) + l − a(J)SNP + w(J)SP , I ′ P = βH←MρIM qSP c(p, q) + l − γIP + a(J)INP − w(J)IP , I ′ NP = βH←MρIM SNP c(p, q) + l − γINP − a(J)INP + w(J)IP , p ′ = a(J) − [a(J) + w(J)] p, I ′ M =… view at source ↗
Figure 9
Figure 9. Figure 9: Comparison between (24) and (45) with model parameters and initial conditions as in Section 5 with l = 0.25, ρ = 2, q = 0 (perfect protection), χ = 104 , n = 1, φ = 20 days, and from bottom to top, w0 = 0.1, 1, 10, 100 (which give Rˆ c ≈ 2.17, 1.85, 1.78, 1.77, respectively). Left: IH (defined as in (40)) as a function of time for model (24) and (45). Right: proportions SP /SH, IP /IH, p and a(J)/[a(J) + w… view at source ↗
Figure 10
Figure 10. Figure 10: Plots for (24) (red for w0 = 1, yellow for w0 = 0.1) and (45) (blue) with model parameters and initial conditions as in Section 5, l = 0.25, ρ = 2, and χ = 104 for, from top to bottom, q = 0, 0.1, 0.2, with n = 1 and φ = 30 days in (25). Left and centre: IH as a function of time. Right: RH as a function of time [PITH_FULL_IMAGE:figures/full_fig_p027_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Plots for (24) (red for w0 = 1, yellow for w0 = 0.1) and (45) (blue) with model parameters and initial conditions as in Section 5, l = 0.25 and ρ = 2, and χ = 104 for, from top to bottom, q = 0, 0.1, 0.2, with n = 2 and φ = 30 days in (25). Left and centre: IH as a function of time. Right: RH as a function of time. 27 [PITH_FULL_IMAGE:figures/full_fig_p027_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Plots for (24) (red for w0 = 1, yellow for w0 = 0.1) and (45) (blue) with model parameters and initial conditions as in Section 5, l = 0.25 and ρ = 2, and χ = 104 for, from top to bottom, q = 0, 0.1, 0.2, with n = 5 and φ = 30 days in (25). Left and : IH as a function of time. Right: RH as a function of time [PITH_FULL_IMAGE:figures/full_fig_p028_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Plots for (24) (red for w0 = 1, yellow for w0 = 0.1) and (45) (blue) with model parameters and initial conditions as in Section 5, χ = 104 , q = 0.5, n = 1 and φ = 30 days in (25). Left: IH as a function of time. Right: RH as a function of time. The dashed lines represent the same simulations for the model without protective behaviour (q = 1). Upper row: l = 0.25 and ρ = 2. Lower row: l = 1 and ρ = 5.12. … view at source ↗
Figure 14
Figure 14. Figure 14: Plots of IH (blue) as a function of time for (45) and EE (green) as a function of SH with model parameters and initial conditions as in Section 5, l = 0.25 and ρ = 2, χ = 104 (upper row), χ = 103 (lower row), q = 0 (left) and q = 0.1 (right), with n = 1 and φ = 30 days in (25) [PITH_FULL_IMAGE:figures/full_fig_p036_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Plots of IH (blue) as a function of time for (45) and EE (green) as a function of SH with model parameters and initial conditions as in Section 5, l = 0.25 and ρ = 2, χ = 104 (upper row), χ = 103 (lower row), q = 0 (left) and q = 0.1 (right), with n = 2 and φ = 30 days in (25). 36 [PITH_FULL_IMAGE:figures/full_fig_p036_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Plots of IH (blue) as a function of time for (45) and EE (green) as a function of SH with model parameters and initial conditions as in Section 5, l = 0.25 and ρ = 2, χ = 104 (upper row), χ = 103 (lower row), q = 0 (left) and q = 0.1 (right), with n = 5 and φ = 30 days in (25) [PITH_FULL_IMAGE:figures/full_fig_p037_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Plots of IH (solid blue) as a function of time for (45), and EE (solid green) and Rˆ e (dashed black) as functions of SH, with model parameters and initial conditions as in Section 5, l = 0.25 and ρ = 2 (left), l = 1 and ρ = 5.12 (right), χ = 104 , q = 0.2, with n = 1 and φ = 30 days in (25). For each plot, the values on the vertical axis on the left (blue colour) are relevant to IH and EE, while those on… view at source ↗
read the original abstract

We investigate a model for a mosquito-borne epidemic in which human hosts may adopt protective behaviour against vector bites in response to information on both past and current disease prevalence. Assuming that mosquitoes can also feed on non-competent hosts (i.e.\ hosts that do not contribute to disease transmission), we first revisit existing results and show that behaviour-driven protection may either decrease or increase the basic reproduction number, depending on the interaction between behavioural response, host composition, and transmission parameters. Assuming that opinion dynamics evolves on a much faster time scale than disease transmission, we then apply Geometric Singular Perturbation Theory to effectively reduce the original two-group model to a model for a homogeneous host population. The reduced system enables a detailed investigation of the impact of information-induced behavioural changes on the transient dynamics of the epidemic, including scenarios in which protective measures lead to outbreaks with low attack rates. Our analysis shows that behavioural responses may either facilitate epidemic control or prolong disease persistence, potentially generating recurrent damped epidemic waves. Numerical simulations are provided to illustrate and support the analytical findings.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The paper develops a compartmental model for mosquito-borne epidemics incorporating information-dependent protective behavior in human hosts and non-competent hosts that mosquitoes can feed upon. It revisits results on the basic reproduction number, showing that protective behavior can either decrease or increase R0 depending on behavioral response strength, host composition, and transmission parameters. Assuming fast opinion dynamics, Geometric Singular Perturbation Theory is applied to reduce the two-group model to an effective homogeneous-host model. The reduced system is then used to analyze transient dynamics, including low attack rates and the potential for recurrent damped epidemic waves, with numerical simulations provided for illustration.

Significance. If the GSPT reduction is valid with uniform attraction and the time-scale separation holds, the work offers a useful framework for examining how information-driven protective behavior modulates epidemic persistence and control in vector-borne diseases. The approach of combining behavioral response modeling with singular perturbation techniques to study transient phenomena such as damped recurrent waves is a methodological strength when the reduction is rigorously justified.

major comments (2)
  1. [GSPT model reduction] In the section applying Geometric Singular Perturbation Theory, the reduction assumes opinion dynamics evolves on a much faster time scale than disease transmission, but no explicit bounds on the separation parameter ε are derived, nor are numerical checks reported confirming that the reduced vector field reproduces key transient quantities (attack rate, damping rate of recurrent waves) for biologically plausible ε values consistent with information spread versus vector-host contact rates. This verification is load-bearing for the central claims about behavioral responses generating recurrent damped waves.
  2. [Basic reproduction number analysis] The claim that behaviour-driven protection may increase R0 (depending on host composition and transmission parameters) is stated in the abstract and introduction as a revisited result, but the manuscript provides no explicit derivation, threshold conditions, or sensitivity analysis in the relevant section to delineate the parameter regimes where the increase occurs versus a decrease.
minor comments (3)
  1. [Model formulation] The notation for the opinion dynamics variable and its coupling to the transmission rate should be introduced with a clear table of symbols in the model formulation section to improve readability.
  2. [Numerical results] Figure captions for the numerical simulations should explicitly list the default parameter values and the range of ε values tested to allow readers to assess the robustness of the illustrated wave patterns.
  3. [Discussion] A brief discussion of the biological plausibility of the assumed functional form for the protective behavior response (e.g., dependence on past vs. current prevalence) would strengthen the interpretation of the results.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful and constructive report. We address each major comment below and outline the revisions we will make to strengthen the manuscript.

read point-by-point responses

  1. Referee: [GSPT model reduction] In the section applying Geometric Singular Perturbation Theory, the reduction assumes opinion dynamics evolves on a much faster time scale than disease transmission, but no explicit bounds on the separation parameter ε are derived, nor are numerical checks reported confirming that the reduced vector field reproduces key transient quantities (attack rate, damping rate of recurrent waves) for biologically plausible ε values consistent with information spread versus vector-host contact rates. This verification is load-bearing for the central claims about behavioral responses generating recurrent damped waves.

    Authors: We agree that the absence of explicit bounds on ε and direct numerical validation for plausible values represents a gap that should be addressed to rigorously support the reduction. In the revised manuscript we will derive approximate bounds on ε by comparing the characteristic time scales of opinion dynamics (governed by information dissemination rates) against vector-host contact rates. We will also add a dedicated numerical subsection comparing the full two-group system and the reduced model for ε ∈ {0.01, 0.05, 0.1}, confirming that attack rates and the damping rates of recurrent waves remain quantitatively close within these ranges, which we argue are consistent with rapid social-media-driven information spread relative to slower epidemiological processes. revision: yes

  2. Referee: [Basic reproduction number analysis] The claim that behaviour-driven protection may increase R0 (depending on host composition and transmission parameters) is stated in the abstract and introduction as a revisited result, but the manuscript provides no explicit derivation, threshold conditions, or sensitivity analysis in the relevant section to delineate the parameter regimes where the increase occurs versus a decrease.

    Authors: The manuscript revisits the basic reproduction number for the two-group model with information-dependent protection in Section 3. Nevertheless, we accept that a clearer, self-contained derivation together with explicit threshold conditions and sensitivity analysis would improve readability and substantiate the claim that protection can increase R0 under certain host compositions and transmission parameters. We will expand the relevant section to include the full next-generation-matrix derivation, the resulting threshold expressions separating the regimes of R0 increase versus decrease, and a brief sensitivity analysis with respect to the behavioral response strength and the proportion of non-competent hosts. revision: yes

Circularity Check

0 steps flagged

No significant circularity; GSPT reduction is an explicit modeling assumption applied to a new system

full rationale

The derivation begins with a standard two-group model incorporating information-dependent protective behavior and non-competent hosts. The central reduction step invokes Geometric Singular Perturbation Theory under the stated assumption that opinion dynamics evolves on a much faster timescale than disease transmission; this is presented as a methodological choice to obtain an effective homogeneous-host model rather than a quantity defined in terms of itself. Subsequent analysis of transient dynamics (attack rates, damped recurrent waves) follows directly from the reduced vector field. No fitted parameters are relabeled as predictions, no load-bearing uniqueness theorems are imported via self-citation, and no ansatz is smuggled through prior work by the same authors. The model remains self-contained against its own equations and external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The model introduces several standard epidemiological parameters (transmission rates, biting rates, recovery rates) plus new behavioral response functions and a separation-of-timescales assumption. No machine-checked proofs or external data are mentioned.

free parameters (2)
  • behavioral response strength and functional form
    The rate and shape of how protection increases with prevalence information is introduced to produce the desired non-monotonic R0 effect and wave behavior.
  • proportion of non-competent hosts
    This fraction is varied to show reversal of protection benefit and is not derived from first principles.
axioms (1)
  • domain assumption Opinion dynamics evolves on a much faster time scale than disease transmission
    Invoked to justify the geometric singular perturbation reduction; appears in the abstract as the premise for model simplification.

pith-pipeline@v0.9.0 · 5489 in / 1544 out tokens · 56848 ms · 2026-05-17T20:04:41.344658+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Kinetic models of opinion-driven epidemic dynamics modulated by graphons

    math.AP 2026-04 unverdicted novelty 7.0

    Kinetic models of opinion-driven epidemics on graphons prove L1 and Sobolev convergence, introduce a structure-preserving scheme, and show that a time-dependent reproduction number analog can generate waves without ex...

Reference graph

Works this paper leans on

76 extracted references · 76 canonical work pages · cited by 1 Pith paper

  1. [1]

    Thermostatted kinetic equations as models for complex systems in physics and life sciences.Phys

    M. Aguiar, V . Anam, K. B. Blyuss, C. D. S. Estadilla, B. V . Guerrero, D. Knopoff, B. W. Kooi, A. K. Srivastav, V . Steindorf, and N. Stollenwerk. Mathematical models for dengue fever epidemiology: A 10-year systematic review.Phys. Life Rev.40 (2022), pp. 65–92.DOI: 10.1016/j.plrev. 2022.02.001

    work page doi:10.1016/j.plrev 2022

  2. [2]

    Alphey, L

    N. Alphey, L. Alphey, and M. B. Bonsall. A model framework to estimate impact and cost of genetics-based sterile insect methods for dengue vector control.PloS one6.10 (2011), e25384.DOI: 10.1371/journal.pone.0025384

    work page doi:10.1371/journal.pone.0025384 2011

  3. [3]

    A. Andò, D. Breda, and G. Gava. How fast is the linear chain trick? A rigorous analysis in the context of behavioral epidemiology.Math. Biosci. Eng.17.5 (2020), pp. 5059–5084.ISSN: 1551-0018.DOI: 10.3934/mbe.2020273

    work page doi:10.3934/mbe.2020273 2020

  4. [4]

    A. Andò, S. De Reggi, F. Scarabel, R. Vermiglio, and J. Wu. Behavior-induced oscillations in epidemic outbreaks with distributed memory: beyond the linear chain trick using numerical methods. Submitted(2025)

    work page 2025

  5. [5]

    M. D. Asfaw, B. Buonomo, and S. M. Kassa. Impact of human behavior on ITNs control strategies to prevent the spread of vector borne diseases.Atti Accad. Peloritana dei Pericolanti Cl. Sci. Fis. Mat. Nat.96.S3 (2018), p. 2.DOI:10.1478/AAPP.96S3A2

    work page doi:10.1478/aapp.96s3a2 2018

  6. [6]

    Bacaër.A Short History of Mathematical Population Dynamics

    N. Bacaër.A Short History of Mathematical Population Dynamics. Springer, London, 2011.DOI: 10.1007/978-0-85729-115-8

    work page doi:10.1007/978-0-85729-115-8 2011

  7. [7]

    C. T Bauch. Imitation dynamics predict vaccinating behaviour.Proc. R. Soc. Lond. B Biol. Sci. 272.1573 (2005), pp. 1669–1675.DOI:10.1098/rspb.2005.3153. 42

    work page doi:10.1098/rspb.2005.3153 2005

  8. [8]

    Bedson, L

    J. Bedson, L. A. Skrip, D. Pedi, S. Abramowitz, S. Carter, M. F. Jalloh, S. Funk, N. Gobat, T. Giles-Vernick, G. Chowell, J. R. de Almeida, R. Elessawi, S. V . Scarpino, R. A. Hammond S., Briand, J. M. Epstein, L. Hébert-Dufresne, and B. M. Althouse. A review and agenda for integrated disease models including social and behavioural factors.Nat. Hum. Behav...

    work page doi:10.1038/s41562-021-01136-2 2021

  9. [9]

    Bolzoni, A

    L. Bolzoni, A. Pugliese, and R. Rosà. The role of heterogeneity on the invasion probability of mosquito-borne diseases in multi-host models.J. Theor. Biol.377 (2015), pp. 25–35.DOI: 10. 1016/j.jtbi.2015.03.027

    work page 2015

  10. [10]

    F. Brauer. A final size relation for epidemic models of vector-transmitted diseases.Infect. Dis. Model.2.1 (2017), pp. 12–20.DOI:10.1016/j.idm.2016.12.001

    work page doi:10.1016/j.idm.2016.12.001 2017

  11. [11]

    F. Brauer. A singular perturbation approach to epidemics of vector-transmitted diseases.Infect. Dis. Model.4 (2019), pp. 115–123.DOI:10.1016/j.idm.2019.04.004

    work page doi:10.1016/j.idm.2019.04.004 2019

  12. [12]

    I. M. Bulai, M. Sensi, and S. Sottile. A geometric analysis of the SIRS compartmental model with fast information and misinformation spreading.Chaos Solit. Fractals185 (2024), p. 115104.DOI: 10.1016/j.chaos.2024.115104

    work page doi:10.1016/j.chaos.2024.115104 2024

  13. [13]

    B. Buonomo. Modeling human response to bed–net promotion campaigns and its impact on malaria transmission. In:Ansari, A. (eds) Adv. Appl. Math.Springer Proc. Math. Stat. Springer, Cham, 2014, pp. 23–30.DOI:https://doi.org/10.1007/978-3-319-06923-4_3

    work page doi:10.1007/978-3-319-06923-4_3 2014

  14. [14]

    Buonomo, E

    B. Buonomo, E. Messina, C. Panico, and A. Vecchio. An integral renewal equation approach to behavioural epidemic models with information index.J. Math. Biol.90.1 (2025), p. 8.DOI: 10.1007/s00285-024-02172-y

    work page doi:10.1007/s00285-024-02172-y 2025

  15. [15]

    Capasso and G

    V . Capasso and G. Serio. A generalization of the Kermack-McKendrick deterministic epidemic model.Math. Biosci.42 (1978), pp. 43–61.DOI:10.1016/0025-5564(78)90006-8

    work page doi:10.1016/0025-5564(78)90006-8 1978

  16. [16]

    Frontiers in Public HealthV olume 10 - 2022(2022) https://doi.org/10.3389/fpubh

    B. Chala and F. Hamde. Emerging and re-emerging vector-borne infectious diseases and the challenges for control: a review.Front. Public Health9 (2021), p. 715759.DOI: 10.3389/fpubh. 2021.715759

    work page doi:10.3389/fpubh 2021

  17. [17]

    Brito da Cruz and H

    A.M.C. Brito da Cruz and H. S. Rodrigues. Personal protective strategies for dengue disease: Simulations in two coexisting virus serotypes scenarios.Math. Comput. Simul.188 (2021), pp. 254– 267.DOI:10.1016/j.matcom.2021.04.002

    work page doi:10.1016/j.matcom.2021.04.002 2021

  18. [18]

    Della Marca, A

    R. Della Marca, A. d’Onofrio, M. Sensi, and S. Sottile. A geometric analysis of the impact of large but finite switching rates on vaccination evolutionary games.Nonlinear Anal. Real World Appl.75 (2024), p. 103986.DOI:10.1016/j.nonrwa.2023.103986

    work page doi:10.1016/j.nonrwa.2023.103986 2024

  19. [19]

    Demers, S

    J. Demers, S. Bewick, J. Calabrese, and W. F. Fagan. Dynamic modelling of personal protection control strategies for vector-borne disease limits the role of diversity amplification.J. R. Soc. Interface.15.145 (2018), p. 20180166.DOI:10.1098/rsif.2018.0166

    work page doi:10.1098/rsif.2018.0166 2018

  20. [20]

    Diekmann and M

    O. Diekmann and M. Gyllenberg. Equations with infinite delay: blending the abstract and the concrete.J. Differ. Equ.252.2 (2012), pp. 819–851.DOI:10.1016/j.jde.2011.09.038

    work page doi:10.1016/j.jde.2011.09.038 2012

  21. [21]

    Diekmann, H

    O. Diekmann, H. Heesterbeek, and T. Britton.Mathematical Tools for Understanding Infectious Disease Dynamics. Princeton Ser. Theor. Comput. Biol. Princeton University Press, Princeton, 2013. DOI:10.1515/9781400845620

    work page doi:10.1515/9781400845620 2013

  22. [22]

    Diekmann, J

    O. Diekmann, J. A. P. Heesterbeek, and J. A. J. Metz. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations.J. Math. Biol.28 (1990), pp. 365–382.DOI:10.1007/BF00178324. 43

    work page doi:10.1007/bf00178324 1990

  23. [23]

    Diekmann, J

    O. Diekmann, J. A. P. Heesterbeek, and M. G. Roberts. The construction of next-generation matrices for compartmental epidemic models.J. R. Soc. Interface.7.47 (2010), pp. 873–885.DOI: 10.1098/rsif.2009.0386

    work page doi:10.1098/rsif.2009.0386 2010

  24. [24]

    M., Gils, S

    O. Diekmann, S. A. Van Gils, S. M. V . Lunel, and H.-O. Walther.Delay Equations: Functional-, complex-, and Nonlinear Analysis. V ol. 110. Applied Mathematical Sciences. Springer, New York, 1995.DOI:10.1007/978-1-4612-4206-2

    work page doi:10.1007/978-1-4612-4206-2 1995

  25. [25]

    Accessed on November 17, 2025.URL: https : / / www

    European Centre for Disease Prevention and Control (ECDC).Dengue. Accessed on November 17, 2025.URL: https : / / www . ecdc . europa . eu / en / all - topics - z / dengue / surveillance- and- disease- data/autochthonous- transmission- dengue- virus-eueea-previous-years

    work page 2025

  26. [26]

    Accessed on November 19, 2025.URL: https://www.ecdc.europa.eu/en/mosquito-borne- diseases

    European Centre for Disease Prevention and Control (ECDC).Mosquito-borne diseases. Accessed on November 19, 2025.URL: https://www.ecdc.europa.eu/en/mosquito-borne- diseases

    work page 2025

  27. [27]

    d’Onofrio and P

    A. d’Onofrio and P. Manfredi. Information-related changes in contact patterns may trigger oscilla- tions in the endemic prevalence of infectious diseases.J. Theor. Biol.256 (2009), pp. 473–478.DOI: 10.1016/j.jtbi.2008.10.005

    work page doi:10.1016/j.jtbi.2008.10.005 2009

  28. [28]

    d’Onofrio and P

    A. d’Onofrio and P. Manfredi. Behavioral SIR models with incidence-based social-distancing. Chaos, Solit. Fractals159 (2022), p. 112072.DOI:10.1016/j.chaos.2022.112072

    work page doi:10.1016/j.chaos.2022.112072 2022

  29. [29]

    Van den Driessche and J

    P. Van den Driessche and J. Watmough. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission.Math. Biosci.180.1-2 (2002), pp. 29–48.DOI: 10.1016/S0025-5564(02)00108-6

    work page doi:10.1016/s0025-5564(02)00108-6 2002

  30. [30]

    Dye and G

    C. Dye and G. Hasibeder. Population dynamics of mosquito-borne disease: effects of flies which bite some people more frequently than others.Trans. R. Soc. Trop. Med. Hyg.80.1 (1986), pp. 69–77. DOI:10.1016/0035-9203(86)90199-9

    work page doi:10.1016/0035-9203(86)90199-9 1986

  31. [31]

    d’Onofrio, P

    A. d’Onofrio, P. Manfredi, and M. Iannelli. Dynamics of partially mitigated multi-phasic epidemics at low susceptible depletion: phases of COVID-19 control in Italy as case study.Math. Biosci.340 (2021), p. 108671.ISSN: 0025-5564.DOI:10.1016/j.mbs.2021.108671

    work page doi:10.1016/j.mbs.2021.108671 2021

  32. [32]

    d’Onofrio, P

    A. d’Onofrio, P. Manfredi, and E. Salinelli. Vaccinating behaviour, information, and the dynamics of SIR vaccine preventable diseases.Theor. Popul. Biol.71.3 (2007), pp. 301–317.DOI: 10.1016/ j.tpb.2007.01.001

    work page 2007

  33. [33]

    Esteva and C

    L. Esteva and C. Vargas. Analysis of a dengue disease transmission model.Math. Biosci.150.2 (1998), pp. 131–151.DOI:10.1016/S0025-5564(98)10003-2

    work page doi:10.1016/s0025-5564(98)10003-2 1998

  34. [34]

    Fenichel

    N. Fenichel. Geometric singular perturbation theory for ordinary differential equations.J. Differ. Equ.31 (1979), pp. 53–98.DOI:10.1016/0022-0396(79)90152-9

    work page doi:10.1016/0022-0396(79)90152-9 1979

  35. [35]

    S. Funk, M. Salathé, and V . Jansen. Modelling the Influence of Human Behaviour on the Spread of Infectious Diseases: A Review.Journal of the Royal Society, Interface / the Royal Society7 (May 2010), pp. 1247–56.DOI:10.1098/rsif.2010.0142

    work page doi:10.1098/rsif.2010.0142 2010

  36. [36]

    U. J. Giménez-Mujica, J. Velázquez-Castro, and A. Anzo-Hernández. Final size of the epidemic for metapopulation vector-borne diseases.J. Math. Anal. Appl.526.1 (2023), p. 127200.DOI: 10.1016/j.jmaa.2023.127200

    work page doi:10.1016/j.jmaa.2023.127200 2023

  37. [37]

    Hasibeder and C

    G. Hasibeder and C. Dye. Population dynamics of mosquito-borne disease: persistence in a com- pletely heterogeneous environment.Theor. Popul. Biol.33.1 (1988), pp. 31–53.DOI: 10.1016/ 0040-5809(88)90003-2

    work page 1988

  38. [38]

    J. A. P. Heesterbeek. A brief history of R0 and a recipe for its calculation.Acta Biotheor.50.3 (2002), pp. 189–204.DOI:10.1023/a:1016599411804. 44

    work page doi:10.1023/a:1016599411804 2002

  39. [39]

    G. Hek. Geometric singular perturbation theory in biological practice.J. Math. Biol.60.3 (2010), pp. 347–386.DOI:10.1007/s00285-009-0266-7

    work page doi:10.1007/s00285-009-0266-7 2010

  40. [40]

    Hu and L

    L. Hu and L. Nie. Stability and Hopf Bifurcation Analysis of a Multi-Delay Vector-Borne Disease Model with Presence Awareness and Media Effect.Fractal Fract.7.12 (2023), p. 831.DOI: 10. 3390/fractalfract7120831

    work page 2023

  41. [41]

    Jardón-Kojakhmetov, C

    H. Jardón-Kojakhmetov, C. Kuehn, A. Pugliese, and M. Sensi. A geometric analysis of the SIR, SIRS and SIRWS epidemiological models.Nonlinear Anal. Real World Appl.58 (2021), p. 103220. DOI:10.1016/j.nonrwa.2020.103220

    work page doi:10.1016/j.nonrwa.2020.103220 2021

  42. [42]

    Jardón-Kojakhmetov, C

    H. Jardón-Kojakhmetov, C. Kuehn, A. Pugliese, and M. Sensi. A geometric analysis of the SIRS epidemiological model on a homogeneous network.J. Math. Biol.83.4 (2021), p. 37.DOI: 10. 1007/s00285-021-01664-5

    work page 2021

  43. [43]

    Kaklamanos, A

    P. Kaklamanos, A. Pugliese, M. Sensi, and S. Sottile. A geometric analysis of the SIRS model with secondary infections.SIAM J. Appl. Math.84.2 (2024), pp. 661–686.DOI: 10 . 1137 / 23M1565632

    work page 2024

  44. [44]

    G. F. Killeen and T. A. Smith. Exploring the contributions of bed nets, cattle, insecticides and excitorepellency to malaria control: a deterministic model of mosquito host-seeking behaviour and mortality.Trans. R. Soc. Trop. Med. Hyg.101.9 (2007), pp. 867–880.DOI: 10.1016/j.trstmh. 2007.04.022

    work page doi:10.1016/j.trstmh 2007

  45. [45]

    An introduction to population protocols

    C. Kuehn.Multiple Time Scale Dynamics. V ol. 191. Springer, Cham, 2015.DOI: 10.1007/978- 3-319-12316-5

    work page doi:10.1007/978- 2015

  46. [46]

    Laxmi, C. N. Ngonghala, and S. Bhattacharyya. An evolutionary game model of individual choices and bed net use: elucidating key aspect in malaria elimination strategies.R. Soc. Open Sci.9.11 (2022), p. 220685.DOI:10.1098/rsos.220685

    work page doi:10.1098/rsos.220685 2022

  47. [47]

    Macdonald.The Epidemiology and Control of Malaria.London: Oxford Univ

    G. Macdonald.The Epidemiology and Control of Malaria.London: Oxford Univ. Pr., 1957

    work page 1957

  48. [48]

    MacDonald.Biological delay systems: linear stability theory

    N. MacDonald.Biological delay systems: linear stability theory. Cambridge Studies in Mathematical Biology. Cambridge Univeristy Press, Cambridge, 1989

    work page 1989

  49. [49]

    Magal, O

    P. Magal, O. Seydi, and G. Webb. Final size of an epidemic for a two-group SIR model.SIAM J. Appl. Math.76.5 (2016), pp. 2042–2059.DOI:10.1137/16M1065392

    work page doi:10.1137/16m1065392 2016

  50. [50]

    Reversibility and asymmetri c conflict in event structures

    P. Magal, O. Seydi, and G. Webb. Final size of a multi-group SIR epidemic model: Irreducible and non-irreducible modes of transmission.Math. Biosci.301 (2018), pp. 59–67.DOI: 10.1016/j. mbs.2018.03.020

    work page doi:10.1016/j 2018

  51. [51]

    Manfredi and A

    P. Manfredi and A. d’Onofrio.Modeling the Interplay Between Human Behavior and the Spread of Infectious Diseases. Springer New York, NY, 2013.DOI:10.1007/978-1-4614-5474-8

    work page doi:10.1007/978-1-4614-5474-8 2013

  52. [52]

    McCallum, N

    H. McCallum, N. Barlow, and J. Hone. How should pathogen transmission be modelled?Trends Ecol. Evol.16.6 (2001), pp. 295–300.DOI:10.1016/S0169-5347(01)02144-9

    work page doi:10.1016/s0169-5347(01)02144-9 2001

  53. [53]

    Journal of The Royal Society Interface 14(135), 20170326 (2017)

    E. Miller, J. Dushoff, and A. Huppert. The risk of incomplete personal protection coverage in vector-borne disease.J. R. Soc. Interface.13.115 (2016), p. 20150666.DOI: 10.1098/rsif. 2015.0666

    work page doi:10.1098/rsif 2016

  54. [54]

    Miller and A

    E. Miller and A. Huppert. The effects of host diversity on vector-borne disease: the conditions under which diversity will amplify or dilute the disease risk.PLoS One8.11 (2013), e80279.DOI: 10.1371/journal.pone.0080279

    work page doi:10.1371/journal.pone.0080279 2013

  55. [55]

    A. K. Misra, A. Sharma, and L. Jia. A mathematical model for control of vector borne diseases through media campaigns.Discrete Contin. Dyn. Syst. - B18.7 (2013).DOI: 10.3934/dcdsb. 2013.18.1909. 45

    work page doi:10.3934/dcdsb 2013

  56. [56]

    S. J. Moore, C. R . Davies, N. Hill, and M. M. Cameron. Are mosquitoes diverted from repellent- using individuals to non-users? Results of a field study in Bolivia.Trop. Med. Int. Health.12.4 (2007), pp. 532–539.DOI:10.1111/j.1365-3156.2006.01811.x

    work page doi:10.1111/j.1365-3156.2006.01811.x 2007

  57. [57]

    E. A. Newton and P. Reiter. A model of the transmission of dengue fever with an evaluation of the impact of ultra-low volume (ULV) insecticide applications on dengue epidemics.Am. J. Trop. Med. Hyg.47.6 (1992), pp. 709–720.DOI:10.4269/ajtmh.1992.47.709

    work page doi:10.4269/ajtmh.1992.47.709 1992

  58. [58]

    Nishiura

    H. Nishiura. Mathematical and statistical analyses of the spread of dengue.Dengue Bull.30 (2006)

    work page 2006

  59. [59]

    S. T. Ogunlade, M. T. Meehan, A. I. Adekunle, and E. S. McBryde. A systematic review of mathematical models of dengue transmission and vector control: 2010–2020.Viruses15.1 (2023), p. 254.DOI:10.3390/v15010254

    work page doi:10.3390/v15010254 2010

  60. [60]

    Pellis, P

    L. Pellis, P. J. Birrell, J. Blake, C. E. Overton, F. Scarabel, H. B. Stage, E. Brooks-Pollock, L. Danon, I. Hall, T. A. House, et al. Estimation of reproduction numbers in real time: conceptual and statistical challenges.J. R. Stat. Soc. Ser. A. Stat. Soc.185.Supplement_1 (2022), S112–S130.DOI: 10.1111/rssa.12955

    work page doi:10.1111/rssa.12955 2022

  61. [61]

    Poletti, B

    P. Poletti, B. Caprile, M. Ajelli, A. Pugliese, and S. Merler. Spontaneous behavioural changes in response to epidemics.J. Theor. Biol.260.1 (2009), pp. 31–40.DOI: 10.1016/j.jtbi.2009. 04.029

    work page doi:10.1016/j.jtbi.2009 2009

  62. [62]

    Pugliese and S

    A. Pugliese and S. De Reggi. A personal overview of epidemic models for mosquito-borne infections. Submitted(2025)

    work page 2025

  63. [63]

    Roosa and N

    K. Roosa and N. H. Fefferman. A general modeling framework for exploring the impact of individual concern and personal protection on vector-borne disease dynamics.Parasit. Vectors.15.1 (2022), p. 361.DOI:10.1186/s13071-022-05481-7

    work page doi:10.1186/s13071-022-05481-7 2022

  64. [64]

    Rosà and A

    R. Rosà and A. Pugliese. Effects of tick population dynamics and host densities on the persistence of tick-borne infections.Math. Biosci.208.1 (2007), pp. 216–240.DOI: 10.1016/j.mbs.2006. 10.002

    work page doi:10.1016/j.mbs.2006 2007

  65. [65]

    Ross.The Prevention of Malaria.London: Churchill, 1911

    R. Ross.The Prevention of Malaria.London: Churchill, 1911

    work page 1911

  66. [66]

    Gambette and K

    S. Schecter. Geometric singular perturbation theory analysis of an epidemic model with spontaneous human behavioral change.J. Math. Biol.82.6 (2021), p. 54.DOI: 10.1007/s00285- 021- 01605-2

    work page doi:10.1007/s00285- 2021

  67. [67]

    D. L. Smith, K. E. Battle, S. I. Hay, C. M. Barker, T. W. Scott, and F. E. McKenzie. Ross, Macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens.PLoS pathogens8.4 (2012), e1002588.DOI:10.1371/journal.ppat.1002588

    work page doi:10.1371/journal.ppat.1002588 2012

  68. [68]

    Thongsripong, J

    P. Thongsripong, J. M. Hyman, D. D. Kapan, and S. N. Bennett. Human–mosquito contact: a missing link in our understanding of mosquito-borne disease transmission dynamics.Ann. Entomol. Soc. Am.114.4 (2021), pp. 397–414.DOI:10.1093/aesa/saab011

    work page doi:10.1093/aesa/saab011 2021

  69. [69]

    Valerio, F

    L. Valerio, F. Marini, G. Bongiorno, L. Facchinelli, M. Pombi, B. Caputo, M. Maroli, and A. Della Torre. Host-feeding patterns of Aedes albopictus (Diptera: Culicidae) in urban and rural contexts within Rome province, Italy.Vector-Borne and Zoonotic Dis.10.3 (2010), pp. 291–294. DOI:10.1089/vbz.2009.0007

    work page doi:10.1089/vbz.2009.0007 2010

  70. [70]

    Z. Wang, C. T. Bauch, S. Bhattacharyya, A. d’Onofrio, P. Manfredi, M. Perc, N. Perra, M. Salathé, and D. Zhao. Statistical physics of vaccination.Physics Reports664 (2016), pp. 1–113.ISSN: 0370-1573.DOI:10.1016/j.physrep.2016.10.006

    work page doi:10.1016/j.physrep.2016.10.006 2016

  71. [71]

    Accessed on November 17, 2025.URL: https://www.who.int/news- room/fact- sheets/detail/vector- borne- diseases

    World Health Organization (WHO).Vector-borne diseases. Accessed on November 17, 2025.URL: https://www.who.int/news- room/fact- sheets/detail/vector- borne- diseases. 46

    work page 2025

  72. [72]

    Accessed on October 16, 2025.URL: https://www

    World Health Organization (WHO).Dengue. Accessed on October 16, 2025.URL: https://www. who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue

    work page 2025

  73. [73]

    M. J. Wonham, M. A. Lewis, J. Rencławowicz, and P. Van den Driessche. Transmission assumptions generate conflicting predictions in host–vector disease models: a case study in West Nile virus.Ecol. Lett.9.6 (2006), pp. 706–725.DOI:10.1111/j.1461-0248.2006.00912.x

    work page doi:10.1111/j.1461-0248.2006.00912.x 2006

  74. [74]

    L. Yakob. How do biting disease vectors behaviourally respond to host availability?Parasit. Vectors. 9.1 (2016), p. 468.DOI:10.1186/s13071-016-1762-4

    work page doi:10.1186/s13071-016-1762-4 2016

  75. [75]

    Zhang, F

    X. Zhang, F. Scarabel, K. Murty, and J. Wu. Renewal equations for delayed population behaviour adaptation coupled with disease transmission dynamics: A mechanism for multiple waves of emerging infections.Math. Biosci.365 (2023), p. 109068.DOI: 10 . 1016 / j . mbs . 2023 . 109068

    work page 2023

  76. [76]

    Zhong, M

    H. Zhong, M. Lima, A. Pugliese, M. Sabbatino, C. Soresina, and T. Tekeli.On the effectiveness of odor-baited traps on mosquito-borne infections. 2025. arXiv:2509.23397. 47

    work page arXiv 2025