Probabilistic simulation of big climate data for robust quantification of changes in compound hazard events.

Understanding changes in extreme compound hazard events is important for climate mitigation and policy. By definition, such events are rare so robust quantification of their future changes is challenging. An approach is presented, for probabilistic modelling and simulation of climate model data, which is invariant to the event definition since it models the underlying weather variables. The approach is based on the idea of a 'moving window' in conjunction with Generalised Additive Models (GAMs) and Bayesian inference. As such, it is robust to the data size and completely parallelizable, while it fully quantifies uncertainty allowing also for comprehensive model checking. Lastly, Gaussian anamorphosis is used to capture dependency across weather variables. The approach results in probabilistic simulations to enable extrapolation beyond the original data range and thus robust quantification of future changes of rare events. We illustrate by application to daily temperature, humidity and precipitation from a regional climate model.

A Bayesian modelling framework to quantify multiple sources of spatial variation for disease mapping.

Spatial connectivity is an important consideration when modelling infectious disease data across a geographical region. Connectivity can arise for many reasons, including shared characteristics between regions and human or vector movement. Bayesian hierarchical models include structured random effects to account for spatial connectivity. However, conventional approaches require the spatial structure to be fully defined prior to model fitting. By applying penalized smoothing splines to coordinates, we create two-dimensional smooth surfaces describing the spatial structure of the data while making minimal assumptions about the structure. The result is a non-stationary surface which is setting specific. These surfaces can be incorporated into a hierarchical modelling framework and interpreted similarly to traditional random effects. Through simulation studies, we show that the splines can be applied to any symmetric continuous connectivity measure, including measures of human movement, and that the models can be extended to explore multiple sources of spatial structure in the data. Using Bayesian inference and simulation, the relative contribution of each spatial structure can be computed and used to generate hypotheses about the drivers of disease. These models were found to perform at least as well as existing modelling frameworks, while allowing for future extensions and multiple sources of spatial connectivity.

The impact of climate suitability, urbanisation, and connectivity on the expansion of dengue in 21st century Brazil.

Dengue is hyperendemic in Brazil, with outbreaks affecting all regions. Previous studies identified geographical barriers to dengue transmission in Brazil, beyond which certain areas, such as South Brazil and the Amazon rainforest, were relatively protected from outbreaks. Recent data shows these barriers are being eroded. In this study, we explore the drivers of this expansion and identify the current limits to the dengue transmission zone. We used a spatio-temporal additive model to explore the associations between dengue outbreaks and temperature suitability, urbanisation, and connectivity to the Brazilian urban network. The model was applied to a binary outbreak indicator, assuming the official threshold value of 300 cases per 100,000 residents, for Brazil's municipalities between 2001 and 2020. We found a nonlinear relationship between higher levels of connectivity to the Brazilian urban network and the odds of an outbreak, with lower odds in metropoles compared to regional capitals. The number of months per year with suitable temperature conditions for Aedes mosquitoes was positively associated with the dengue outbreak occurrence. Temperature suitability explained most interannual and spatial variation in South Brazil, confirming this geographical barrier is influenced by lower seasonal temperatures. Municipalities that had experienced an outbreak previously had double the odds of subsequent outbreaks. We identified geographical barriers to dengue transmission in South Brazil, western Amazon, and along the northern coast of Brazil. Although a southern barrier still exists, it has shifted south, and the Amazon no longer has a clear boundary. Few areas of Brazil remain protected from dengue outbreaks. Communities living on the edge of previous barriers are particularly susceptible to future outbreaks as they lack immunity. Control strategies should target regions at risk of future outbreaks as well as those currently within the dengue transmission zone.

Spatial connectivity in mosquito-borne disease models: a systematic review of methods and assumptions.

Spatial connectivity plays an important role in mosquito-borne disease transmission. Connectivity can arise for many reasons, including shared environments, vector ecology and human movement. This systematic review synthesizes the spatial methods used to model mosquito-borne diseases, their spatial connectivity assumptions and the data used to inform spatial model components. We identified 248 papers eligible for inclusion. Most used statistical models (84.2%), although mechanistic are increasingly used. We identified 17 spatial models which used one of four methods (spatial covariates, local regression, random effects/fields and movement matrices). Over 80% of studies assumed that connectivity was distance-based despite this approach ignoring distant connections and potentially oversimplifying the process of transmission. Studies were more likely to assume connectivity was driven by human movement if the disease was transmitted by an Aedes mosquito. Connectivity arising from human movement was more commonly assumed in studies using a mechanistic model, likely influenced by a lack of statistical models able to account for these connections. Although models have been increasing in complexity, it is important to select the most appropriate, parsimonious model available based on the research question, disease transmission process, the spatial scale and availability of data, and the way spatial connectivity is assumed to occur.

Spatial co-localisation of extreme weather events: a clear and present danger.

Extreme weather events have become a dominant feature of the narrative surrounding changes in global climate with large impacts on ecosystem stability, functioning and resilience; however, understanding of their risk of co-occurrence at the regional scale is lacking. Based on the UK Met Office's long-term temperature and rainfall records, we present the first evidence demonstrating significant increases in the magnitude, direction of change and spatial co-localisation of extreme weather events since 1961. Combining this new understanding with land-use data sets allowed us to assess the likely consequences on future agricultural production and conservation priority areas. All land-uses are impacted by the increasing risk of at least one extreme event and conservation areas were identified as the hotspots of risk for the co-occurrence of multiple event types. Our findings provide a basis to regionally guide land-use optimisation, land management practices and regulatory actions preserving ecosystem services against multiple climate threats.

A modelling approach for correcting reporting delays in disease surveillance data.

One difficulty for real-time tracking of epidemics is related to reporting delay. The reporting delay may be due to laboratory confirmation, logistical problems, infrastructure difficulties, and so on. The ability to correct the available information as quickly as possible is crucial, in terms of decision making such as issuing warnings to the public and local authorities. A Bayesian hierarchical modelling approach is proposed as a flexible way of correcting the reporting delays and to quantify the associated uncertainty. Implementation of the model is fast due to the use of the integrated nested Laplace approximation. The approach is illustrated on dengue fever incidence data in Rio de Janeiro, and severe acute respiratory infection data in the state of Paraná, Brazil.

Light pollution is associated with earlier tree budburst across the United Kingdom.

The ecological impact of night-time lighting is of concern because of its well-demonstrated effects on animal behaviour. However, the potential of light pollution to change plant phenology and its corresponding knock-on effects on associated herbivores are less clear. Here, we test if artificial lighting can advance the timing of budburst in trees. We took a UK-wide 13 year dataset of spatially referenced budburst data from four deciduous tree species and matched it with both satellite imagery of night-time lighting and average spring temperature. We find that budburst occurs up to 7.5 days earlier in brighter areas, with the relationship being more pronounced for later-budding species. Excluding large urban areas from the analysis showed an even more pronounced advance of budburst, confirming that the urban 'heat-island' effect is not the sole cause of earlier urban budburst. Similarly, the advance in budburst across all sites is too large to be explained by increases in temperature alone. This dramatic advance of budburst illustrates the need for further experimental investigation into the impact of artificial night-time lighting on plant phenology and subsequent species interactions. As light pollution is a growing global phenomenon, the findings of this study are likely to be applicable to a wide range of species interactions across the world.

Transmission parameters estimated for Salmonella typhimurium in swine using susceptible-infectious-resistant models and a Bayesian approach.

BACKGROUND: Transmission models can aid understanding of disease dynamics and are useful in testing the efficiency of control measures. The aim of this study was to formulate an appropriate stochastic Susceptible-Infectious-Resistant/Carrier (SIR) model for Salmonella Typhimurium in pigs and thus estimate the transmission parameters between states. RESULTS: The transmission parameters were estimated using data from a longitudinal study of three Danish farrow-to-finish pig herds known to be infected. A Bayesian model framework was proposed, which comprised Binomial components for the transition from susceptible to infectious and from infectious to carrier; and a Poisson component for carrier to infectious. Cohort random effects were incorporated into these models to allow for unobserved cohort-specific variables as well as unobserved sources of transmission, thus enabling a more realistic estimation of the transmission parameters. In the case of the transition from susceptible to infectious, the cohort random effects were also time varying. The number of infectious pigs not detected by the parallel testing was treated as unknown, and the probability of non-detection was estimated using information about the sensitivity and specificity of the bacteriological and serological tests. The estimate of the transmission rate from susceptible to infectious was 0.33 [0.06, 1.52], from infectious to carrier was 0.18 [0.14, 0.23] and from carrier to infectious was 0.01 [0.0001, 0.04]. The estimate for the basic reproduction ration (R0) was 1.91 [0.78, 5.24]. The probability of non-detection was estimated to be 0.18 [0.12, 0.25]. CONCLUSIONS: The proposed framework for stochastic SIR models was successfully implemented to estimate transmission rate parameters for Salmonella Typhimurium in swine field data. R0 was 1.91, implying that there was dissemination of the infection within pigs of the same cohort. There was significant temporal-cohort variability, especially at the susceptible to infectious stage. The model adequately fitted the data, allowing for both observed and unobserved sources of uncertainty (cohort effects, diagnostic test sensitivity), so leading to more reliable estimates of transmission parameters.

Osteoporotic hip fractures: bisphosphonates sales and observed turning point in trend. A population-based retrospective study.

The aim is to examine the temporal trends of hip fracture incidence in Portugal by sex and age groups, and explore the relation with anti-osteoporotic medication. From the National Hospital Discharge Database, we selected from 1st January 2000 to 31st December 2008, 77,083 hospital admissions (77.4% women) caused by osteoporotic hip fractures (low energy, patients over 49years-age), with diagnosis codes 820.x of ICD 9-CM. The 2001 Portuguese population was used as standard to calculate direct age-standardized incidence rates (ASIR) (100,000 inhabitants). Generalized additive and linear models were used to evaluate and quantify temporal trends of age specific rates (AR), by sex. We identified 2003 as a turning point in the trend of ASIR of hip fractures in women. After 2003, the ASIR in women decreased on average by 10.3 cases/100,000 inhabitants, 95% CI (-15.7 to -4.8), per 100,000 anti-osteoporotic medication packages sold. For women aged 65-69 and 75-79 we identified the same turning point. However, for women aged over 80, the year 2004 marked a change in the trend, from an increase to a decrease. Among the population aged 70-74 a linear decrease of incidence rate (95% CI) was observed in both sexes, higher for women: -28.0% (-36.2 to -19.5) change vs -18.8%, (-32.6 to -2.3). The abrupt turning point in the trend of ASIR of hip fractures in women is compatible with an intervention, such as a medication. The trends were different according to gender and age group, but compatible with the pattern of bisphosphonates sales.

Assessing risk profiles for Salmonella serotypes in breeding pig operations in Portugal using a Bayesian hierarchical model.

BACKGROUND: The EU Regulation No 2160/2003 imposes a reduction in the prevalence of Salmonella in pigs. The efficiency of control programmes for Salmonella in pigs, reported among the EU Member States, varies and definitive eradication seems very difficult. Control measures currently recommended for Salmonella are not serotype-specific. Is it possible that the risk factors for different Salmonella serotypes are different? The aim of this study was to investigate potential risk factors for two groups of Salmonella sp serotypes using pen faecal samples from breeding pig holdings representative of the Portuguese pig sector. METHODS: The data used come from the Baseline Survey for the Prevalence of Salmonella in breeding pigs in Portugal. A total of 1670 pen faecal samples from 167 herds were tested, and 170 samples were positive for Salmonella. The presence of Salmonella in each sample (outcome variable) was classified in three categories: i) no Salmonella, ii) Salmonella Typhimurium or S. Typhimurium-like strains with the antigenic formula: 1,4,5,12:i:-, , and iii) other serotypes. Along with the sample collection, a questionnaire concerning herd management and potential risk factors was utilised. The data have a "natural" hierarchical structure so a categorical multilevel analysis of the dataset was carried out using a Bayesian hierarchical model. The model was estimated using Markov Chain Monte Carlo methods, implemented in the software WinBUGS. RESULTS: The significant associations found (when compared to category "no Salmonella"), for category "serotype Typhimurium or S. Typhimurium-like strains with the antigenic formula: 1,4,5,12:i:-" were: age of breeding sows, size of the herd, number of pigs/pen and source of semen. For the category "other serotypes" the significant associations found were: control of rodents, region of the country, source of semen, breeding sector room and source of feed. CONCLUSIONS: The risk factors significantly associated with Salmonella shedding from the category "serotype Typhimurium or serotype 1,4,5,12:i:-" were more related to animal factors, whereas those associated with "other serotypes" were more related to environmental factors. Our findings suggest that different control measures could be used to control different Salmonella serotypes in breeding pigs.