A diverse portfolio of marine protected areas can better advance global conservation and equity.

Marine protected areas (MPAs) are widely used for ocean conservation, yet the relative impacts of various types of MPAs are poorly understood. We estimated impacts on fish biomass from no-take and multiple-use (fished) MPAs, employing a rigorous matched counterfactual design with a global dataset of >14,000 surveys in and around 216 MPAs. Both no-take and multiple-use MPAs generated positive conservation outcomes relative to no protection (58.2% and 12.6% fish biomass increases, respectively), with smaller estimated differences between the two MPA types when controlling for additional confounding factors (8.3% increase). Relative performance depended on context and management: no-take MPAs performed better in areas of high human pressure but similar to multiple-use in remote locations. Multiple-use MPA performance was low in high-pressure areas but improved significantly with better management, producing similar outcomes to no-take MPAs when adequately staffed and appropriate use regulations were applied. For priority conservation areas where no-take restrictions are not possible or ethical, our findings show that a portfolio of well-designed and well-managed multiple-use MPAs represents a viable and potentially equitable pathway to advance local and global conservation.

A cross-sectional analysis of medical conditions and environmental factors associated with fractional exhaled nitric oxide (FeNO) in women and children from the ISA birth cohort, Costa Rica.

BACKGROUND: Fractional exhaled nitric oxide (FeNO) is a marker of airway inflammation. Elevated FeNO has been associated with environmental exposures, however, studies from tropical countries are limited. Using data from the Infants' Environmental Health Study (ISA) birth cohort, we evaluated medical conditions and environmental exposures' association with elevated FeNO. METHODS: We performed a cross-sectional analysis of 277 women and 293 8-year old children who participated in the 8-year post-partum visit in 2019. We measured FeNO and collected information on medical conditions and environmental exposures including smoke from waste burning, work in banana plantations, and home pesticide use. We defined elevated FeNO as >25 ppb for women and >20 ppb for children. To evaluate factors associated with elevated FeNO, we used logistic regression models adjusted for obesity in women and unadjusted in children. RESULTS: Overall elevated FeNO was common (20% of women, 13% of children). Rhinitis diagnosis was significantly associated with elevated FeNO in both women (odds ratio (OR): 3.67 95% Confidence Interval (CI): 1.81,7.35) and children (OR: 8.18 95%CI: 3.15, 21.22); wheeze was associated with elevated FeNO in women (OR: 4.50 95% CI: 2.25, 8.99). Environmental exposures were associated with elevated FeNO, but not significantly. Waste burning was associated with elevated FeNO in both women (OR: 1.58 95%CI 0.68, 4.15) and children (OR: 2.49 95%CI:0.82, 10.79). Para-occupational pesticide exposures were associated with elevated FeNO in women and children. For women, having a partner working in agriculture was associated with elevated FeNO (OR: 1.61 95%CI:0.77, 3.58) and for children, maternal work in agriculture was associated with elevated FeNO. (OR 2.08 95%CI 0.86, 4.67) CONCLUSION: Rhinitis and wheeze were associated with elevated FeNO in this rural, agricultural population. Smoke from waste burning as well as para-occupational pesticide exposure may contribute to elevated FeNO in rural communities.

A nonparametric test of group distributional differences for hierarchically clustered functional data.

Biological sex and gender are critical variables in biomedical research, but are complicated by the presence of sex-specific natural hormone cycles, such as the estrous cycle in female rodents, typically divided into phases. A common feature of these cycles are fluctuating hormone levels that induce sex differences in many behaviors controlled by the electrophysiology of neurons, such as neuronal membrane potential in response to electrical stimulus, typically summarized using a priori defined metrics. In this paper, we propose a method to test for differences in the electrophysiological properties across estrous cycle phase without first defining a metric of interest. We do this by modeling membrane potential data in the frequency domain as realizations of a bivariate process, also depending on the electrical stimulus, by adopting existing methods for longitudinal functional data. We are then able to extract the main features of the bivariate signals through a set of basis function coefficients. We use these coefficients for testing, adapting methods for multivariate data to account for an induced hierarchical structure that is a product of the experimental design. We illustrate the performance of the proposed approach in simulations and then apply the method to experimental data.

Bayesian nonparametric quantile process regression and estimation of marginal quantile effects.

Flexible estimation of multiple conditional quantiles is of interest in numerous applications, such as studying the effect of pregnancy-related factors on low and high birth weight. We propose a Bayesian nonparametric method to simultaneously estimate noncrossing, nonlinear quantile curves. We expand the conditional distribution function of the response in I-spline basis functions where the covariate-dependent coefficients are modeled using neural networks. By leveraging the approximation power of splines and neural networks, our model can approximate any continuous quantile function. Compared to existing models, our model estimates all rather than a finite subset of quantiles, scales well to high dimensions, and accounts for estimation uncertainty. While the model is arbitrarily flexible, interpretable marginal quantile effects are estimated using accumulative local effect plots and variable importance measures. A simulation study shows that our model can better recover quantiles of the response distribution when the data are sparse, and an analysis of birth weight data is presented.

Spectral adjustment for spatial confounding.

Adjusting for an unmeasured confounder is generally an intractable problem, but in the spatial setting it may be possible under certain conditions. We derive necessary conditions on the coherence between the exposure and the unmeasured confounder that ensure the effect of exposure is estimable. We specify our model and assumptions in the spectral domain to allow for different degrees of confounding at different spatial resolutions. One assumption that ensures identifiability is that confounding present at global scales dissipates at local scales. We show that this assumption in the spectral domain is equivalent to adjusting for global-scale confounding in the spatial domain by adding a spatially smoothed version of the exposure to the mean of the response variable. Within this general framework, we propose a sequence of confounder adjustment methods that range from parametric adjustments based on the Matérn coherence function to more robust semiparametric methods that use smoothing splines. These ideas are applied to areal and geostatistical data for both simulated and real datasets.

A discrete-time survival model for porcine epidemic diarrhoea virus.

Since the arrival of porcine epidemic diarrhea virus (PEDV) in the United States in 2013, elimination and control programmes have had partial success. The dynamics of its spread are hard to quantify, though previous work has shown that local transmission and the transfer of pigs within production systems are most associated with the spread of PEDV. Our work relies on the history of PEDV infections in a region of the southeastern United States. This infection data is complemented by farm-level features and extensive industry data on the movement of both pigs and vehicles. We implement a discrete-time survival model and evaluate different approaches to modelling the local-transmission and network effects. We find strong evidence in that the local-transmission and pig-movement effects are associated with the spread of PEDV, even while controlling for seasonality, farm-level features and the possible spread of disease by vehicles. Our fully Bayesian model permits full uncertainty quantification of these effects. Our farm-level out-of-sample predictions have a receiver-operating characteristic area under the curve (AUC) of 0.779 and a precision-recall AUC of 0.097. The quantification of these effects in a comprehensive model allows stakeholders to make more informed decisions about disease prevention efforts.

Estimating intervention effects on infectious disease control: The effect of community mobility reduction on Coronavirus spread.

Understanding the effects of interventions, such as restrictions on community and large group gatherings, is critical to controlling the spread of COVID-19. Susceptible-Infectious-Recovered (SIR) models are traditionally used to forecast the infection rates but do not provide insights into the causal effects of interventions. We propose a spatiotemporal model that estimates the causal effect of changes in community mobility (intervention) on infection rates. Using an approximation to the SIR model and incorporating spatiotemporal dependence, the proposed model estimates a direct and indirect (spillover) effect of intervention. Under an interference and treatment ignorability assumption, this model is able to estimate causal intervention effects, and additionally allows for spatial interference between locations. Reductions in community mobility were measured by cell phone movement data. The results suggest that the reductions in mobility decrease Coronavirus cases 4 to 7 weeks after the intervention.

Pyrimethanil and chlorpyrifos air concentrations and pregnant women's urinary metabolites in the Infants' Environmental Health Study (ISA), Costa Rica.

BACKGROUND: Only few studies have compared environmental pesticide air concentrations with specific urinary metabolites to evaluate pathways of exposure. Therefore, we compared pyrimethanil and chlorpyrifos concentrations in air with urinary 4-hydroxypyrimethanil (OHP, metabolite of pyrimethanil) and 3,5,6-trichloro-2-pyridinol (TCPy, metabolite of chlorpyrifos) among pregnant women from the Infant's Environmental Health Study (ISA) in Matina County, Costa Rica. METHODS: During pregnancy, we obtained repeat urinary samples from 448 women enrolled in the ISA study. We extrapolated pyrimethanil and chlorpyrifos concentrations measured with passive air samplers (PAS) (n = 48, from 12 schools), across space and time using a Bayesian spatiotemporal model. We subsequently compared these concentrationswith urinary OHP and TCPy in 915 samples from 448 women, usingseparatemixed models andconsidering several covariables. RESULTS: A 10% increase in air pyrimethanil (ng/m3) was associated with a 5.7% (95% confidence interval (CI 4.6, 6.8) increase in OHP (µg/L). Women living further from banana plantations had lower OHP: -0.7% (95% CI -1.2, -0.3) for each 10% increase in distance (meters) as well as women who ate rice and beans ≥15 times a week -23% (95% CI -38, -4). In addition, each 1 ng/m3 increase in chlorpyrifos in air was associated with a 1.5% (95% CI 0.2, 2.8) increase in TCPy (µg/L), and women working in agriculture tended to have increased TCPy (21%, 95% CI -2, 49). CONCLUSION: The Bayesian spatiotemporal models were useful to estimate pyrimethanil and chlorpyrifos air concentrations across space and time. Our results suggest inhalation of pyrimethanil and chlorpyrifos is a pathway of environmental exposure. PAS seems a useful technique to monitor environmental current-use pesticide exposures. For future studies, we recommend increasing the number of locations of environmental air measurements, obtaining all air and urine measurements during the same month, and, ideally, including dermal exposure estimates as well.

Soil Properties and Moisture Synergistically Influence Nontuberculous Mycobacterial Prevalence in Natural Environments of Hawai'i.

Nontuberculous mycobacteria (NTM) are opportunistic pathogens that cause chronic pulmonary disease (PD). NTM infections are thought to be acquired from the environment; however, the basal environmental factors that drive and sustain NTM prevalence are not well understood. The highest prevalence of NTM PD cases in the United States is reported from Hawai'i, which is unique in its climate and soil composition, providing an opportunity to investigate the environmental drivers of NTM prevalence. We used microbiological sampling and spatial logistic regression complemented with fine-scale soil mineralogy to model the probability of NTM presence across the natural landscape of Hawai'i. Over 7 years, we collected and microbiologically cultured 771 samples from 422 geographic sites in natural areas across the Hawaiian Islands for the presence of NTM. NTM were detected in 210 of these samples (27%), with Mycobacterium abscessus being the most frequently isolated species. The probability of NTM presence was highest in expansive soils (those that swell with water) with a high water balance (>1-m difference between rainfall and evapotranspiration) and rich in Fe-oxides/hydroxides. We observed a positive association between NTM presence and iron in wet soils, supporting past studies, but no such association in dry soils. High soil-water balance may facilitate underground movement of NTM into the aquifer system, potentially compounded by expansive capabilities allowing crack formation under drought conditions, representing further possible avenues for aquifer infiltration. These results suggest both precipitation and soil properties are mechanisms by which surface NTM may reach the human water supply. IMPORTANCE Nontuberculous mycobacteria (NTM) are ubiquitous in the environment, being found commonly in soils and natural bodies of freshwater. However, little is known about the environmental niches of NTM and how they relate to NTM prevalence in homes and other human-dominated areas. To characterize NTM environmental associations, we collected and cultured 771 samples from 422 geographic sites in natural areas across Hawai'i, the U.S. state with the highest prevalence of NTM pulmonary disease. We show that the environmental niches of NTM are most associated with highly expansive, moist soils containing high levels of iron oxides/hydroxides. Understanding the factors associated with NTM presence in the natural environment will be crucial for identifying potential mechanisms and risk factors associated with NTM infiltration into water supplies, which are ultimately piped into homes where most exposure risk is thought to occur.

Geostatistical modeling of positive-definite matrices: An application to diffusion tensor imaging.

Geostatistical modeling for continuous point-referenced data has extensively been applied to neuroimaging because it produces efficient and valid statistical inference. However, diffusion tensor imaging (DTI), a neuroimaging technique characterizing the brain's anatomical structure, produces a positive-definite (p.d.) matrix for each voxel. Currently, only a few geostatistical models for p.d. matrices have been proposed because introducing spatial dependence among p.d. matrices properly is challenging. In this paper, we use the spatial Wishart process, a spatial stochastic process (random field), where each p.d. matrix-variate random variable marginally follows a Wishart distribution, and spatial dependence between random matrices is induced by latent Gaussian processes. This process is valid on an uncountable collection of spatial locations and is almost-surely continuous, leading to a reasonable way of modeling spatial dependence. Motivated by a DTI data set of cocaine users, we propose a spatial matrix-variate regression model based on the spatial Wishart process. A problematic issue is that the spatial Wishart process has no closed-form density function. Hence, we propose an approximation method to obtain a feasible Cholesky decomposition model, which we show to be asymptotically equivalent to the spatial Wishart process model. A local likelihood approximation method is also applied to achieve fast computation. The simulation studies and real data application demonstrate that the Cholesky decomposition process model produces reliable inference and improved performance, compared to other methods.

A spatiotemporal recommendation engine for malaria control.

Malaria is an infectious disease affecting a large population across the world, and interventions need to be efficiently applied to reduce the burden of malaria. We develop a framework to help policy-makers decide how to allocate limited resources in realtime for malaria control. We formalize a policy for the resource allocation as a sequence of decisions, one per intervention decision, that map up-to-date disease related information to a resource allocation. An optimal policy must control the spread of the disease while being interpretable and viewed as equitable to stakeholders. We construct an interpretable class of resource allocation policies that can accommodate allocation of resources residing in a continuous domain and combine a hierarchical Bayesian spatiotemporal model for disease transmission with a policy-search algorithm to estimate an optimal policy for resource allocation within the pre-specified class. The estimated optimal policy under the proposed framework improves the cumulative long-term outcome compared with naive approaches in both simulation experiments and application to malaria interventions in the Democratic Republic of the Congo.

A SPATIAL CAUSAL ANALYSIS OF WILDLAND FIRE-CONTRIBUTED PM2.5 USING NUMERICAL MODEL OUTPUT.

Wildland fire smoke contains hazardous levels of fine particulate matter (PM2.5), a pollutant shown to adversely effect health. Estimating fire attributable PM2.5 concentrations is key to quantifying the impact on air quality and subsequent health burden. This is a challenging problem since only total PM2.5 is measured at monitoring stations and both fire-attributable PM2.5 and PM2.5 from all other sources are correlated in space and time. We propose a framework for estimating fire-contributed PM2.5 and PM2.5 from all other sources using a novel causal inference framework and bias-adjusted chemical model representations of PM2.5 under counterfactual scenarios. The chemical model representation of PM2.5 for this analysis is simulated using Community Multiscale Air Quality Modeling System (CMAQ), run with and without fire emissions across the contiguous U.S. for the 2008-2012 wildfire seasons. The CMAQ output is calibrated with observations from monitoring sites for the same spatial domain and time period. We use a Bayesian model that accounts for spatial variation to estimate the effect of wildland fires on PM2.5 and state assumptions under which the estimate has a valid causal interpretation. Our results include estimates of the contributions of wildfire smoke to PM2.5 for the contiguous U.S. Additionally, we compute the health burden associated with the PM2.5 attributable to wildfire smoke.

Landscape-level variation in Bt crops predict Helicoverpa zea (Lepidoptera: Noctuidae) resistance in cotton agroecosystems.

BACKGROUND: Helicoverpa zea (Boddie) damage to Bt cotton and maize has increased as a result of widespread Bt resistance across the USA Cotton Belt. Our objective was to link Bt crop production patterns to cotton damage through a series of spatial and temporal surveys of commercial fields to understand how Bt crop production relates to greater than expected H. zea damage to Bt cotton. To do this, we assembled longitudinal cotton damage data that spanned the Bt adoption period, collected cotton damage data since Bt resistance has been detected, and estimated local population susceptibility using replicated on-farm studies that included all Bt pyramids marketed in cotton. RESULTS: Significant year effects of H. zea damage frequency in commercial cotton were observed throughout the Bt adoption period, with a recent damage increase after 2012. Landscape-level Bt crop production intensity over time was positively associated with the risk of H. zea damage in two- and three-toxin pyramided Bt cotton. Helicoverpa zea damage also varied across Bt toxin types in spatially replicated on-farm studies. CONCLUSIONS: Landscape-level predictors of H. zea damage in Bt cotton can be used to identify heightened Bt resistance risk areas and serves as a model to understand factors that drive pest resistance evolution to Bt toxins in the southeastern United States. These results provide a framework for more effective insect resistance management strategies to be used in combination with conventional pest management practices that improve Bt trait durability while minimizing the environmental footprint of row crop agriculture. © 2021 Society of Chemical Industry. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Statistical downscaling with spatial misalignment: Application to wildland fire PM2.5 concentration forecasting.

Fine particulate matter, PM2.5, has been documented to have adverse health effects and wildland fires are a major contributor to PM2.5 air pollution in the US. Forecasters use numerical models to predict PM2.5 concentrations to warn the public of impending health risk. Statistical methods are needed to calibrate the numerical model forecast using monitor data to reduce bias and quantify uncertainty. Typical model calibration techniques do not allow for errors due to misalignment of geographic locations. We propose a spatiotemporal downscaling methodology that uses image registration techniques to identify the spatial misalignment and accounts for and corrects the bias produced by such warping. Our model is fitted in a Bayesian framework to provide uncertainty quantification of the misalignment and other sources of error. We apply this method to different simulated data sets and show enhanced performance of the method in presence of spatial misalignment. Finally, we apply the method to a large fire in Washington state and show that the proposed method provides more realistic uncertainty quantification than standard methods.

Multivariate Spatial Prediction of Air Pollutant Concentrations with INLA.

Estimates of daily air pollution concentrations with complete spatial and temporal coverage are important for supporting epidemiologic studies and health impact assessments. While numerous approaches have been developed for modeling air pollution, they typically only consider each pollutant separately. We describe a spatial multipollutant data fusion model that combines monitoring measurements and chemical transport model simulations that leverages dependence between pollutants to improve spatial prediction. For the contiguous United States, we created a data product of daily concentration for 12 pollutants (CO, NOx, NO2, SO2, O3, PM10, and PM2.5 species EC, OC, NO3, NH4, SO4) during the period 2005 to 2014. Out-of-sample prediction showed good performance, particularly for daily PM2.5 species EC (R2 = 0.64), OC (R2 = 0.75), NH4 (R2 = 0.84), NO3 (R2 = 0.73), and SO4 (R2 = 0.80). By employing the integrated nested Laplace approximation (INLA) for Bayesian inference, our approach also provides model-based prediction error estimates. The daily data product at 12km spatial resolution will be publicly available immediately upon publication. To our knowledge this is the first publicly available data product for major PM2.5 species and several gases at this spatial and temporal resolution.

A Review of Spatial Causal Inference Methods for Environmental and Epidemiological Applications.

The scientific rigor and computational methods of causal inference have had great impacts on many disciplines but have only recently begun to take hold in spatial applications. Spatial causal inference poses analytic challenges due to complex correlation structures and interference between the treatment at one location and the outcomes at others. In this paper, we review the current literature on spatial causal inference and identify areas of future work. We first discuss methods that exploit spatial structure to account for unmeasured confounding variables. We then discuss causal analysis in the presence of spatial interference including several common assumptions used to reduce the complexity of the interference patterns under consideration. These methods are extended to the spatiotemporal case where we compare and contrast the potential outcomes framework with Granger causality and to geostatistical analyses involving spatial random fields of treatments and responses. The methods are introduced in the context of observational environmental and epidemiological studies and are compared using both a simulation study and analysis of the effect of ambient air pollution on COVID-19 mortality rate. Code to implement many of the methods using the popular Bayesian software OpenBUGS is provided.

Flexible and Fast Spatial Return Level Estimation Via a Spatially Fused Penalty.

Spatial extremes are common for climate data as the observations are usually referenced by geographic locations and dependent when they are nearby. An important goal of extremes modeling is to estimate the T-year return level. Among the methods suitable for modeling spatial extremes, perhaps the simplest and fastest approach is the spatial generalized extreme value (GEV) distribution and the spatial generalized Pareto distribution (GPD) that assume marginal independence and only account for dependence through the parameters. Despite the simplicity, simulations have shown that return level estimation using the spatial GEV and spatial GPD still provides satisfactory results compared to max-stable processes, which are asymptotically justified models capable of representing spatial dependence among extremes. However, the linear functions used to model the spatially varying coefficients are restrictive and may be violated. We propose a flexible and fast approach based on the spatial GEV and spatial GPD by introducing fused lasso and fused ridge penalty for parameter regularization. This enables improved return level estimation for large spatial extremes compared to the existing methods.