seasonality is modeled by a forcing that reduces the basic reproductive number by a factor α min ranging from 0.1 to 1 (no seasonality) [20] . The forcing is described by a sinusoidal function of 12 months-period that reaches its peak during Winter time and its minimum during Summer time in each Hemisphere, with the two Hemispheres with opposite phases. Given the population and mobility data, infection dynamics parameters, and initial conditions, GLEaM performs the simulation of stochastic realizations of the worldwide unfolding of the epidemic. From these in silico epidemics a variety of information can be gathered, such as prevalence, morbidity, number of secondary cases, number of imported cases, hospitalized patients, amounts of drugs used, and other quantities for each subpopulation with a time resolution of 1 day.