Bayesian Hierarchical Modeling Framework Research Articles

Interactions between runoff volume, timing, and annual temperatures shape migration phenology of a threatened adfluvial sucker

AbstractMigratory fish species have been particularly impacted by changes to the hydrologic and climatic cues to which migration and spawning behaviours have been adapted across generations. While conservation and recovery programs increasingly implement flow management actions to promote successful migration and spawning, uncertainty regarding how spawning migration phenology responds to changing environmental conditions can limit the ability to effectively target such recovery actions. Here, we use a Bayesian hierarchical modelling framework to analyse spawning migration phenology of individually tagged June suckers (Chasmistes liorus) – a federally threatened, long‐lived, iteroparous, adfluvial species endemic to Utah Lake (Utah, USA) and its tributaries. We then examine how annual hydrologic and thermal conditions relate to different components of annual migration phenology, including peak migration date, in‐stream residence time, and among‐individual variation in migration timing. Peak migration date occurred earlier in years with warmer spring air temperatures (a proxy for water temperatures), though this effect interacted with peak runoff timing. Both residence time and among‐individual variation in migration timing were greater in years with larger spring discharge and later peak flows. Residence time was also longer in warmer years. These results highlight how natural and anthropogenic changes to river flow and thermal regimes are likely to impact June sucker migration timing and duration, and our approach can be applied to other migratory species to identify the external drivers of the different components of migration phenology.

Improved hierarchical Bayesian modeling framework with arbitrary polynomial chaos for probabilistic model updating

Bayesian finite element model updating techniques have found widespread application in the structural health monitoring. The conventional Bayesian modeling framework (CBMF) identifies the posterior distribution of structural parameters using data from a single experiment. However, structural parameters exhibit variability due to changes in environmental or experimental conditions, an aspect overlooked by CBMF. In recent years, the hierarchical Bayesian modeling framework (HBMF) has gradually attracted attention. This framework incorporates an additional layer of hyperparameters and utilizes multiple sets of experimental data to describe the variability in structural parameters caused by changes in experimental conditions. Unfortunately, the increase in hyperparameters and data size leads the HBMF to face the dilemma of exceedingly high computational costs. To address this issue, this study proposes an improved hierarchical Bayesian modeling framework (IHBMF). To enhance the computational efficiency of IHBMF, the importance sampling method is employed to simplify the computation of the likelihood function for hyperparameters. Subsequently, a surrogate model based on arbitrary polynomial chaos expansion is introduced to further reduce the computational cost of IHBMF. The framework's accuracy and efficiency were validated through analysis of a simply supported beam and a steel pedestrian bridge. The results demonstrate that IHBMF not only delivers a more accurate quantification of structural parameter uncertainty than CBMF but also exhibits higher computational efficiency compared to HBMF, showcasing its superior capability in structural health monitoring applications.

Bayesian Forecasting of Mortality Rates for Small Areas Using Spatiotemporal Models.

Estimation and prediction of subnational mortality rates for small areas are essential planning tools for studying health inequalities. Standard methods do not perform well when data are noisy, a typical behavior of subnational datasets. Thus, reliable estimates are difficult to obtain. I present a Bayesian hierarchical model framework for prediction of mortality rates at a small or subnational level. By combining ideas from demography and epidemiology, the classical mortality modeling framework is extended to include an additional spatial component capturing regional heterogeneity. Information is pooled across neighboring regions and smoothed over time and age. To make predictions more robust and address the issue of model selection, a Bayesian version of stacking is considered using leave-future-out validation. I apply this method to forecast mortality rates for 96 regions in Bavaria, Germany, disaggregated by age and sex. Uncertainty surrounding the forecasts is provided in terms of prediction intervals. Using posterior predictive checks, I show that the models capture the essential features and are suitable to forecast the data at hand. On held-out data, my predictions outperform those of standard models lacking a regional component.

Drought-related wildfire accounts for one-third of the forest wildfires in subtropical China

The threat of wildfires to human safety and ecosystem services is widely known to increase due to climate change characterized by drought. In subtropical China, seasonal drought during the non-monsoon period contributes to increase of fire risk. However, both the relationship between drought and wildfire and the extent of drought effect are poorly understood. Based on satellite active fire detections and monthly drought metrics, we quantified the impact of drought on forest wildfire frequency in subtropical China during 2001–2020, using a Bayesian hierarchical modelling framework. We observed excess wildfire, defined here as above normal wildfire activity level, is significantly correlated to drought across subtropical China, and drought-wildfire relationship showed threshold behavior that the probability of excess wildfire increased significantly when the vapor pressure deficit (VPD) crossed 0.24 standard deviations of its long-term average. From 2001 to 2020, excess wildfire caused by drought accounted for approximately 31 % of the total number of forest fire points during the fire season (November to May). The drought-related excess wildfire was concentrated in 2007–2010 and 2014, and their spatial distribution showed a dipole between the western and eastern subtropical China in individual years. Moreover, we evaluated the spatiotemporal sensitivity of wildfire to drought by Bayesian model coefficients. There was higher sensitivity than average level of the entire region in 54 % of subtropical China, especially in northwestern Yunnan. Our results indicate that drought plays a dominant role on wildfire regime at the regional scale in subtropical China. As the drying trend intensifies in the future, wildfires will show different patterns due to the large differences in the sensitivity of wildfire to drought in subtropical China.

A hierarchical Bayesian modeling framework for identification of Non-Gaussian processes

Non-Gaussian processes are frequently encountered in engineering problems, posing a challenge when it comes to identification. The main challenge in the identification arises from the fact that a non-Gaussian process can be treated as a collection of infinite dimensional non-Gaussian variables. The application of the hierarchical Bayesian modeling (HBM) framework is constrained due to the inherent complexity of dimensionality and non-Gaussian characteristics associated with these variables. To tackle the issue of dimensionality, the improved orthogonal series expansion (iOSE) representing a non-Gaussian process by time functions with non-Gaussian coefficients, which are readily obtained from discretizing the process at some specific time points, is introduced within the HBM framework. In particular, the iOSE is embedded to convert the identification of a non-Gaussian process into a finite number of non-Gaussian coefficients. Regarding their non-Gaussian characteristics, polynomial chaos expansion (PCE) is used to quantify the uncertainty of the non-Gaussian coefficients with parameters in PCE treated as hyper parameters to be estimated by the HBM framework. The proposed framework is applicable to the identification of both stationary and nonstationary non-Gaussian processes. The effectiveness of non-Gaussian process quantification by the proposed framework is demonstrated using simulated data of a non-stationary extreme value process. Theapplicability of this approach for non-Gaussian process identification is validated by accurately identifying a stochastic load in a structural dynamic problem. Furthermore, it is successfully applied to the reconstruction of random mode shapes of a building arising from different environmental conditions.

How effective were Australian Quarantine Stations in mitigating transmission aboard ships during the influenza pandemic of 1918-19?

The influenza pandemic of 1918-19 was the most devastating pandemic of the 20th century. It killed an estimated 50-100 million people worldwide. In late 1918, when the severity of the disease was apparent, the Australian Quarantine Service was established. Vessels returning from overseas and inter-state were intercepted, and people were examined for signs of illness and quarantined. Some of these vessels carried the infection throughout their voyage and cases were prevalent by the time the ship arrived at a Quarantine Station. We study four outbreaks that took place on board the Medic, Boonah, Devon, and Manuka in late 1918. These ships had returned from overseas and some of them were carrying troops that served in the First World War. By analysing these outbreaks under a stochastic Bayesian hierarchical modeling framework, we estimate the transmission rates among crew and passengers aboard these ships. Furthermore, we ask whether the removal of infectious, convalescent, and healthy individuals after arriving at a Quarantine Station in Australia was an effective public health response.

Susceptible windows of prenatal and postnatal fine particulate matter exposures and attention-deficit hyperactivity disorder symptoms in early childhood

Few prior studies have explored windows of susceptibility to fine particulate matter (PM2.5) in both the prenatal and postnatal periods and children's attention-deficit/hyperactivity disorder (ADHD) symptoms.We analyzed data from 1416 mother-child pairs from the Spanish INMA (INfancia y Medio Ambiente) Study (2003–2008). Around 5 years of age, teachers reported the number of ADHD symptoms (i.e., inattention, hyperactivity/impulsivity) using the ADHD Diagnostic and Statistical Manual of Mental Disorders. Around 7 years of age, parents completed the Conners' Parent Rating Scales, from which we evaluated the ADHD index, cognitive problems/inattention, hyperactivity, and oppositional subscales, reported as age- and sex-standardized T-scores. Daily residential PM2.5 exposures were estimated using a two-stage random forest model with temporal back-extrapolation and averaged over 1-week periods in the prenatal period and 4-week periods in the postnatal period. We applied distributed lag non-linear models within the Bayesian hierarchical model framework to identify susceptible windows of prenatal or postnatal exposure to PM2.5 (per 5-μg/m3) for ADHD symptoms. Models were adjusted for relevant covariates, and cumulative effects were reported by aggregating risk ratios (RRcum) or effect estimates (βcum) across adjacent susceptible windows.A similar susceptible period of exposure to PM2.5 (1.2–2.9 and 0.9–2.7 years of age, respectively) was identified for hyperactivity/impulsivity symptoms assessed ~5 years (RRcum = 2.72, 95% credible interval [CrI] = 1.98, 3.74) and increased hyperactivity subscale ~7 years (βcum = 3.70, 95% CrI = 2.36, 5.03). We observed a susceptibility period to PM2.5 on risk of hyperactivity/impulsivity symptoms ~5 years in gestational weeks 16–22 (RRcum = 1.36, 95% CrI = 1.22, 1.52). No associations between PM2.5 exposure and other ADHD symptoms were observed.We report consistent evidence of toddlerhood as a susceptible window of PM2.5 exposure for hyperactivity in young children. Although mid-pregnancy was identified as a susceptible period of exposure on hyperactivity symptoms in preschool-aged children, this association was not observed at the time children were school-aged.

A Bayesian hierarchical model framework to quantify uncertainty of tropical cyclone precipitation forecasts

A Bayesian hierarchical model framework to quantify uncertainty of tropical cyclone precipitation forecasts

Virtual sensing based on Hierarchical Bayesian Modeling framework using a Laplace-based Gibbs sampler

This paper presents a Bayesian-based approach to realize virtual strain sensing and quantify the uncertainty in these estimates using sparse output-only measurements. This method relies on output-only measurements to identify the modal coordinates of the structure and realizes the virtual sensing according to the modal superposition method. A robust Hierarchical Bayesian modeling (HBM) framework is developed to fully account for the uncertainties arising from modeling errors and measurement noise. Moreover, the HBM framework has the ability to quantify the uncertainty of parameters in the statistical model by introducing an additional layer of the prior distribution. The Gibbs sampler is adopted to obtain the marginal posterior distributions of the hyper-parameters, with the implementation of Gibbs methods relying on the derived analytical expressions of the full conditional posterior distributions. The Laplace asymptotic approximation simplified the form of the full conditional posterior distributions, thus improving the computational efficiency. Given the hyper-parameters’ marginal posterior distribution, the dynamic parameters’ marginal posterior distribution can be calculated according to the total probability theorem. Finally, one can propagate the uncertainty to make predictions of quantities of interest (QoI). Three case studies, including two numerical examples and a laboratory model test, are carried out to verify the accuracy and efficiency of the proposed algorithm. The results show that the estimates of the parameters are close to the true values. Besides, the prediction results of QoI also show high accuracy, verifying the effectiveness of the proposed HBM method in virtual sensing.

Forecasting Magnitude and Frequency of Seasonal Streamflow Extremes Using a Bayesian Hierarchical Framework

AbstractWe develop a space‐time Bayesian hierarchical modeling (BHM) framework for two flood risk attributes—seasonal daily maximum flow and the number of events that exceed a threshold during a season (NEETM)—at a suite of gauge locations on a river network. The model uses generalized extreme value (GEV) and Poisson distributions as marginals for these flood attributes with non‐stationary parameters. The rate parameters of the Poisson distribution and location, scale, and shape parameters of the GEV are modeled as linear functions of suitable covariates. Gaussian copulas are applied to capture the spatial dependence. The best covariates are selected using the Watanabe‐Akaike information criterion (WAIC). The modeling framework results in the posterior distribution of the flood attributes at all the gauges and various lead times. We demonstrate the utility of this modeling framework to forecast the flood risk attributes during the summer peak monsoon season (July‐August) at five gauges in the Narmada River basin (NRB) of West‐Central India for several lead times (0–3 months). As potential covariates, we consider climate indices such as El Niño–Southern Oscillation (ENSO), the Indian Ocean Dipole (IOD), and the Pacific Warm Pool Region (PWPR) from antecedent seasons, which have shown strong teleconnections with the Indian monsoon. We also include new indices related to the East Pacific and West Indian Ocean regions depending on the lead times. We show useful long lead skill from this modeling approach which has a strong potential to enable robust risk‐based flood mitigation and adaptation strategies 3 months before flood occurrences.

Understanding the impact of mobility on COVID-19 spread: A hybrid gravity-metapopulation model of COVID-19.

The outbreak of the severe acute respiratory syndrome coronavirus 2 started in Wuhan, China, towards the end of 2019 and spread worldwide. The rapid spread of the disease can be attributed to many factors including its high infectiousness and the high rate of human mobility around the world. Although travel/movement restrictions and other non-pharmaceutical interventions aimed at controlling the disease spread were put in place during the early stages of the pandemic, these interventions did not stop COVID-19 spread. To better understand the impact of human mobility on the spread of COVID-19 between regions, we propose a hybrid gravity-metapopulation model of COVID-19. Our modeling framework has the flexibility of determining mobility between regions based on the distances between the regions or using data from mobile devices. In addition, our model explicitly incorporates time-dependent human mobility into the disease transmission rate, and has the potential to incorporate other factors that affect disease transmission such as facemasks, physical distancing, contact rates, etc. An important feature of this modeling framework is its ability to independently assess the contribution of each factor to disease transmission. Using a Bayesian hierarchical modeling framework, we calibrate our model to the weekly reported cases of COVID-19 in thirteen local health areas in Metro Vancouver, British Columbia (BC), Canada, from July 2020 to January 2021. We consider two main scenarios in our model calibration: using a fixed distance matrix and time-dependent weekly mobility matrices. We found that the distance matrix provides a better fit to the data, whilst the mobility matrices have the ability to explain the variance in transmission between regions. This result shows that the mobility data provides more information in terms of disease transmission than the distances between the regions.

Associations between ambient air pollutants and childhood hand, foot, and mouth disease in Sichuan, China: a spatiotemporal study

Hand, foot, and mouth disease (HFMD) is a major public health concern in the Asia–Pacific region. Previous studies have implied that ambient air pollution may affect the incidence of HFMD, but the results among different regions are inconsistent. We aimed to deepen the understanding of the associations between air pollutants and HFMD by conducting a multicity study. Daily data on childhood HFMD counts and meteorological and ambient air pollution (PM2.5, PM10, NO2, CO, O3, and SO2) concentrations in 21 cities in Sichuan Province from 2015 to 2017 were collected. A spatiotemporal Bayesian hierarchical model framework was established, and then a distributed lag nonlinear models (DLNMs) was constructed to reveal exposure-lag-response relationships between air pollutants and HFMD while controlling for spatiotemporal effects. Furthermore, given the differences in the levels and seasonal trends of air pollutants between the basin region and plateau region, we explored whether these associations varied between different areas (basin and plateau). The associations between air pollutants and HFMD were nonlinear, with different lag responses. Low NO2 concentrations and both low and high PM2.5 and PM10 concentrations were associated with a decreased risk of HFMD. No significant associations between CO, O3, and SO2 and HFMD were found. The associations between air pollutant concentrations and HFMD were different between the basin and plateau regions. Our study revealed associations between PM2.5, PM10, and NO2 concentrations and HFMD, deepening the understanding of the relationships between air pollutants and HFMD. These findings provide evidence to support the formulation of relevant prevention measures and the establishment of an early warning system.

Improving Bayesian radiological profiling of waste drums using Dirichlet priors, Gaussian process priors, and hierarchical modeling

We present three methodological improvements of our recently proposed approach for Bayesian inference of the radionuclide inventory in radioactive waste drums, from radiological measurements. First we resort to the Dirichlet distribution for the prior distribution of the isotopic vector. The Dirichlet distribution possesses the attractive property that the elements of its vector samples sum up to 1. Second, we demonstrate that such Dirichlet priors can be incorporated within an hierarchical modeling of the prior uncertainty in the isotopic vector, when prior information about isotopic composition is available. Our used Bayesian hierarchical modeling framework makes use of this available information but also acknowledges its uncertainty by letting to a controlled extent the information content of the indirect measurement data (i.e., gamma and neutron counts) shape the actual prior distribution of the isotopic vector. Third, we propose to regularize the Bayesian inversion by using Gaussian process (GP) prior modeling when inferring 1D spatially-distributed mass or, equivalently, activity distributions. As of uncertainty in the efficiencies, we keep using the same stylized drum modeling approach as proposed in our previous work to account for the source distribution uncertainty across the vertical direction of the drum. A series of synthetic tests followed by application to a real waste drum show that combining hierarchical modeling of the prior isotopic composition uncertainty together with GP prior modeling of the vertical Pu profile across the drum works well. We also find that our GP prior can handles both cases with and without spatial correlation. Of course, our GP prior modeling framework only makes sense in the context of spatial inference. Furthermore, the computational times involved by our approach are on the order of a few hours, say about 2, to provide uncertainty estimates for all variables of interest in the considered inverse problem. This warrants further investigations to speed up the inference.

An analytically tractable solution for hierarchical Bayesian model updating with variational inference scheme

The hierarchical Bayesian modelling (HBM) framework has recently been proposed to properly account for the model parameter uncertainty in structural dynamics. This framework postulates a hierarchical prior for the model parameters that depend on the hyper parameters to be identified using the multiple datasets. However, the number of hyper parameters increases linearly with the number of model parameters as well as the number of datasets, making the framework computationally challenging and analytically intractable. To this end, this study employs a variational inference scheme within the framework, deriving explicit expressions for the posterior distributions of the hyper parameters and further the predictive distribution of the model parameters so as to avoid resorting to sampling-based strategy, thus efficiently enhancing the computational efficiency of the HBM framework. In particular, the posterior distribution of the hyper parameters is derived as a Normal-Inverse-Wishart (NIW) distribution, and the predictive distribution of the model parameters are expressed as a normal distribution conditional on the samples of its hyper parameters. Such derivations offer valuable insights into the interpretations of different sources of uncertainties existing in the procedure of model updating. Specifically, it reveals that the ensemble uncertainty of the model parameters consists of both the identification uncertainty obtained from each dataset as well as the test-to-test variability across the overall datasets. Additionally, it suggests that, when the test-to-test variability dominates the uncertainty of the model parameters, the HBM framework with negligible identification uncertainty can reduce to the same expression with a frequentist perspective. Linear models with modal properties data and nonlinear models with time histories data of building systems are utilized to verify and demonstrate the effectiveness of the proposed framework. Results indicate that the proposed formulations provide sufficiently accurate solutions with efficient computational performance, compared to the full sampling approach regarded as the reference of the HBM framework.

Using multi-sourced big data to correlate sleep deprivation and road traffic noise: A US county-level ecological study

BackgroundRoad traffic noise is a serious public health problem globally as it has adverse psychological and physiologic effects (i.e., sleep). Since previous studies mainly focused on individual levels, we aim to examine associations between road traffic noise and sleep deprivation on a large scale; namely, the US at county level. MethodsInformation from a large-scale sleep survey and national traffic noise map, both obtained from government's open data, were utilized and processed with Geographic Information System (GIS) techniques. To examine the associations between traffic noise and sleep deprivation, we used a hierarchical Bayesian spatial modelling framework to simultaneously adjust for multiple socioeconomic factors while accounting for spatial correlation. FindingsWith 62.90% of people not getting enough sleep, a 10 dBA increase in average sound-pressure level (SPL) or Ls10 (SPL of the relatively noisy area) in a county, was associated with a 49% (OR: 1.49; 95% CrIs:1.19–1.86) or 8% (1.08; 1.00–1.16) increase in the odds of a person in a particular county not getting enough sleep. No significant association was observed for Ls90 (SPL of the relatively quiet area). A 10% increase in noise exposure area or population ratio was associated with a 3% (1.03; 1.01–1.06) or 4% (1.04; 1.02–1.06) increase in the odds of a person within a county not getting enough sleep. InterpretationTraffic noise can contribute to variations in sleep deprivation among counties. This study suggests that policymakers could set up different noise-management strategies for relatively quiet and noisy areas and incorporate geospatial noise indicators, such as exposure population or area ratio. Furthermore, urban planners should consider urban sprawl patterns differently in terms of noise-induced sleep problems.

Under age five children survival times in Nigeria: a Bayesian spatial modeling approach

BackgroundNigeria is among the top five countries in the world with the highest under-five mortality rates. In addition to the general leading causes of under-five mortality, studies have shown that disparity in sociocultural values and practices across ethnic groups in Nigeria influence child survival, thus there is a need for scientific validation. This study quantified the survival probabilities and the impact of socioeconomic and demographic factors, proximate and biological determinants, and environmental factors on the risk of under-five mortality in Nigeria.MethodsThe Kaplan-Meier survival curve, Nelson Aalen hazard curve, and components survival probabilities were estimated. The Exponential, Gamma, Log-normal, Weibull, and Cox hazard models in a Bayesian mixed effect hierarchical hazard modeling framework with spatial components were considered, and the Deviance and Watanabe Akaike information criteria were used to select the best model for inference. A 5% level of significance was assumed throughout this work. The 2018 Nigeria Demographic and Health Survey dataset was used, and the outcome variable was the time between birth and death or birth and the date of interview for children who were alive on the day of the interview.ResultsFindings show that the probability of a child dying within the first two months is 0.04, and the probability of a boy child dying before attaining age five is 0.106, while a girl child is 0.094 probability. Gender, maternal education, household wealth status, source of water and toilet facility, residence, mass media, frequency of antenatal and postnatal visits, marital status, place of delivery, multiple births, who decide healthcare use, use of bednet are significant risk factors of child mortality in Nigeria. The mortality risk is high among the maternal age group below 24 and above 44years, and birth weight below 2.5Kg and above 4.5Kg. The under-five mortality risk is severe in Kebbi, Kaduna, Jigawa, Adamawa, Gombe, Kano, Kogi, Nasarawa, Plateau, and Sokoto states in Nigeria.ConclusionThis study accentuates the need for special attention for the first two months after childbirth as it is the age group with the highest expected mortality. A practicable way to minimize death in the early life of children is to improve maternal healthcare service, promote maternal education, encourage delivery in healthcare facilities, positive parental attitude to support multiple births, poverty alleviation programs for the less privileged, and a prioritized intervention to Northern Nigeria.

Are declining auction prices a result of bidder behavior or auctioneer strategy? Evidence from cattle auctions

AbstractWe test the ability of three previously advanced mechanisms to generate the declining price anomaly frequently observed in sequential auctions by employing a Bayesian hierarchical model framework to simulate heterogeneous bidders in sequential cattle auctions. One hypothesis is that declining prices are due to convergence in the buyers' bidding functions as bidders learn from each other. The second posits that as bidders fulfill their demand by winning earlier lots, the decrease in aggregate demand causes prices to fall. The third is the auction designer's strategic arrangement of lots. Results show the most plausible mechanism to be a strategic ordering of the lots.

Bayesian Hierarchical Modelling for Uncertainty Quantification in Operational Thermal Resistance of LED Systems

Remaining useful life (RUL) prediction is central to prognostics and reliability assessment of light-emitting diode (LED) systems. Their unknown long-term service life remaining when subject to specific operating conditions is affected by various sources of uncertainty stemming from production of individual system components, application of the whole system, measurement and operation. To enhance the reliability of model-based predictions, it is essential to account for all of these uncertainties in a systematic manner. This paper proposes a Bayesian hierarchical modelling framework for inverse uncertainty quantification (UQ) in LED operation under thermal loading. The main focus is on the LED systems’ operational thermal resistances, which are subject to system and application variability. Posterior inference is based on a Markov chain Monte Carlo (MCMC) sampling scheme using the Metropolis–Hastings (MH) algorithm. Performance of the method is investigated for simulated data, which allow to focus on different UQ aspects in applications. Findings from an application scenario in which the impact of disregarded uncertainty on RUL prediction is discussed highlight the need for a comprehensive UQ to allow for reliable predictions.

Can GPS and GRACE data be used to separate past and present-day surface loading in a data-driven approach?

SUMMARY Glacial isostatic adjustment (GIA) and the hydrological cycle are both associated with mass changes and vertical land motion (VLM), which are observed by GRACE and GPS, respectively. Hydrology-related VLM results from the instantaneous response of the elastic solid Earth to surface loading by freshwater, whereas GIA-related VLM reveals the long-term response of the viscoelastic Earth mantle to past ice loading history. Thus, observations of mass changes and VLM are interrelated, making GIA and hydrology difficult to quantify and study independently. In this work, we investigate the feasibility of separating these processes based on GRACE and GPS observations, in a fully data-driven and physically consistent approach. We take advantage of the differences in the spatio-temporal characteristics of the GIA and hydrology fields to estimate the respective contributions of each component using a Bayesian hierarchical modelling framework. A closed-loop synthetic test confirms that our method successfully solves this source separation problem. However, there are significant challenges when applying the same approach with actual observations and the answer to the main question of this study is more nuanced. In particular, in regions where GPS station coverage is sparse, the lack of informative data becomes a limiting factor.

Synergies between environmental degradation and climate variation on malaria re-emergence in southern Venezuela: a spatiotemporal modelling study

Environmental degradation facilitates the emergence of vector-borne diseases, such as malaria, through changes in the ecological landscape that increase human-vector contacts and that expand vector habitats. However, the modifying effects of environmental degradation on climate-disease relationships have not been well explored. Here, we investigate the rapid re-emergence of malaria in a transmission hotspot in southern Venezuela and explore the synergistic effects of environmental degradation, specifically gold-mining activity, and climate variation. In this spatiotemporal modelling study of the 46 parishes of the state of Bolívar, southeast Venezuela, we used data from the Venezuelan Ministry of Health including population data and monthly cases of Plasmodium falciparum malaria and Plasmodium vivax malaria between 1996 and 2016. We estimated mean precipitation and temperature using the ERA5-Land dataset and used monthly anomalies in sea-surface temperature as an indicator of El Niño events between 1996 and 2016. The location of suspected mining sites in Bolívar in 2009, 2017, and 2018 were sourced from the Amazon Geo-Referenced Socio-Environmental Information Network. We estimated measures of cumulative forest loss and urban development by km2 using annual land cover maps from the European Space Agency Climate Change Initiative between 1996 and 2016. We modelled monthly cases of P falciparum and P vivax malaria using a Bayesian hierarchical mixed model framework. We quantified the variation explained by mining activity before exploring the modifying effects of environmental degradation on climate-malaria relationships. We observed a 27% reduction in the additional unexplained spatial variation in incidence of P falciparum malaria and a 23% reduction in P vivax malaria when mining was included in our models. The effect of temperature on malaria was greater in high mining areas than low mining areas, and the P falciparum malaria effect size at temperatures of 26·5°C (2·4 cases per 1000 people [95% CI 1·78-3·06]) was twice as high as the effect in low mining areas (1 case per 1000 people [0·68-1·49]). We show that mining activity in southern Venezuela is associated with hotspots of malaria transmission. Increased temperatures exacerbated malaria transmission in mining areas, highlighting the need to consider how environmental degradation modulates climate effect on disease risk, which is especially important in areas subjected to rapidly rising temperatures and land-use change globally. Our findings have implications for the progress towards malaria elimination in the Latin American region. Our findings are also important for effectively targeting timely treatment programmes and vector-control activities in mining areas with high rates of malaria transmission. Biotechnology and Biological Sciences Research Council, Royal Society, US National Institutes of Health, and Global Challenges Research Fund. For the Spanish translation of the abstract see Supplementary Materials section.