Reduced-rank spatio-temporal modeling of air pollution concentrations in the Multi-Ethnic Study of Atherosclerosis and Air Pollution

There is growing evidence in the epidemiologic literature of the relationship between air pollution and adverse health outcomes. Prediction of individual air pollution exposure in the Environmental Protection Agency (EPA) funded Multi-Ethnic Study of Atheroscelerosis and Air Pollution (MESA Air) study relies on a flexible spatio-temporal prediction model that integrates land-use regression with kriging to account for spatial dependence in pollutant concentrations. Temporal variability is captured using temporal trends estimated via modified singular value decomposition and temporally varying spatial residuals. This model utilizes monitoring data from existing regulatory networks and supplementary MESA Air monitoring data to predict concentrations for individual cohort members. In general, spatio-temporal models are limited in their efficacy for large data sets due to computational intractability. We develop reduced-rank versions of the MESA Air spatio-temporal model. To do so, we apply low-rank kriging to account for spatial variation in the mean process and discuss the limitations of this approach. As an alternative, we represent spatial variation using thin plate regression splines. We compare the performance of the outlined models using EPA and MESA Air monitoring data for predicting concentrations of oxides of nitrogen (NO$_x$)-a pollutant of primary interest in MESA Air-in the Los Angeles metropolitan area via cross-validated $R^2$. Our findings suggest that use of reduced-rank models can improve computational efficiency in certain cases. Low-rank kriging and thin plate regression splines were competitive across the formulations considered, although TPRS appeared to be more robust in some settings.Comment: Published in at http://dx.doi.org/10.1214/14-AOAS786 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Application of Risk within Net Present Value Calculations for Government Projects

In January 2004, President Bush announced a new vision for space exploration. This included retirement of the current Space Shuttle fleet by 2010 and the development of new set of launch vehicles. The President's vision did not include significant increases in the NASA budget, so these development programs need to be cost conscious. Current trade study procedures address factors such as performance, reliability, safety, manufacturing, maintainability, operations, and costs. It would be desirable, however, to have increased insight into the cost factors behind each of the proposed system architectures. This paper reports on a set of component trade studies completed on the upper stage engine for the new launch vehicles. Increased insight into architecture costs was developed by including a Net Present Value (NPV) method and applying a set of associated risks to the base parametric cost data. The use of the NPV method along with the risks was found to add fidelity to the trade study and provide additional information to support the selection of a more robust design architecture

Pragmatic Estimation of a Spatio-Temporal Air Quality Model With Irregular Monitoring Data

Statistical analyses of the health effects of air pollution have increasingly used GIS-based covariates for prediction of ambient air quality in “land-use” regression models. More recently these regression models have accounted for spatial correlation structure in combining monitoring data with land-use covariates. The current paper builds on these concepts to address spatio-temporal prediction of ambient concentrations of particulate matter with aerodynamic diameter less than 2.5 μm (PM2.5) on the basis of a model representing spatially varying seasonal trends and spatial correlation structures. Our hierarchical methodology provides a pragmatic approach that fully exploits regulatory and other supplemental monitoring data which jointly define a complex spatio-temporal monitoring design. We explain the elements of the computational approach, including estimation of smoothed empirical orthogonal functions (SEOFs) as basis functions for temporal trend, spatial (“land use”) regression by Partial Least Squares (PLS), modeling of spatio-temporal correlation structure, and generalized universal kriging prediction of ambient exposure for subjects in the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air) project. Analyses are demonstrated in detail for the South California study area of the MESA Air project using AQS monitoring data from 2000 to 2006 and supplemental MESA Air monitoring data beginning in 2005. Results of application of the modeling and estimation methodology are presented also for five other MESA Air metropolitan study areas across the country with comments on current and future research developments

Predicting Intra-Urban Variation in Air Pollution Concentrations with Complex Spatio-Temporal Interactions

We describe a methodology for assigning individual estimates of long-term average air pollution concentrations that accounts for a complex spatio-temporal correlation structure and can accommodate unbalanced observations. This methodology has been developed as part of the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air), a prospective cohort study funded by the U.S. EPA to investigate the relationship between chronic exposure to air pollution and cardiovascular disease. Our hierarchical model decomposes the space-time field into a “mean” that includes dependence on covariates and spatially varying seasonal and long-term trends and a “residual” that accounts for spatially correlated deviations from the mean model. The model accommodates complex spatio-temporal patterns by characterizing the temporal trend at each location as a linear combination of empirically derived temporal basis functions, and embedding the spatial fields of coefficients for the basis functions in separate linear regression models with spatially correlated residuals (universal kriging). This approach allows us to implement a scalable single-stage estimation procedure that easily accommodates a significant number of missing observations at some monitoring locations. We apply the model to predict long-term average concentrations of oxides of nitrogen (NOx) from 2005-2007 in the Los Angeles area, based on data from 18 EPA Air Quality System regulatory monitors. The cross-validated R2 is 0.67. The MESA Air study is also collecting additional concentration data as part of a supplementary monitoring campaign. We describe the sampling plan and demonstrate in a simulation study that the additional data will contribute to improved predictions of long-term average concentrations

Perspectives on open access high resolution digital elevation models to produce global flood hazard layers

Global flood hazard models have recently become a reality thanks to the release of open access global digital elevation models, the development of simplified and highly efficient flow algorithms, and the steady increase in computational power. In this commentary we argue that although the availability of open access global terrain data has been critical in enabling the development of such models, the relatively poor resolution and precision of these data now limit significantly our ability to estimate flood inundation and risk for the majority of the planet’s surface. The difficulty of deriving an accurate ‘bare-earth’ terrain model due to the interaction of vegetation and urban structures with the satellite-based remote sensors means that global terrain data are often poorest in the areas where people, property (and thus vulnerability) are most concentrated. Furthermore, the current generation of open access global terrain models are over a decade old and many large floodplains, particularly those in developing countries, have undergone significant change in this time. There is therefore a pressing need for a new generation of high resolution and high vertical precision open access global digital elevation models to allow significantly improved global flood hazard models to be developed

Historical Prediction Modeling Approach for Estimating Long-Term Concentrations of PM in Cohort Studies Before the 1999 Implementation of Widespread Monitoring

Introduction: Recent cohort studies use exposure prediction models to estimate the association between long-term residential concentrations of PM2.5 and health. Because these prediction models rely on PM2.5 monitoring data, predictions for times before extensive spatial monitoring present a challenge to understanding long-term exposure effects. The Environmental Protection Agency (EPA) Federal Reference Method (FRM) network for PM2.5 was established in 1999. We evaluated a novel statistical approach to produce high quality exposure predictions from 1980-2010 for epidemiological applications. Methods: We developed spatio-temporal prediction models using geographic predictors and annual average PM2.5 data from 1999 through 2010 from the FRM and the Interagency Monitoring of Protected Visual Environments (IMPROVE) networks. The model consists of a spatially-varying long-term mean, a spatially-varying temporal trend, and spatially-varying and temporally-independent spatio-temporal residuals structured using a universal kriging framework. Temporal trends in annual averages of PM2.5 before 1999 were estimated by using a) extrapolation based on PM2.5 data for 1999-2010 in FRM/IMPROVE, b) PM2.5 sulfate data for 1987-2010 in the Clean Air Status and Trends Network, and c) visibility data for 1980-2010 across the Weather-Bureau-Army-Navy network. We validated the resulting models using PM2.5 data collected before 1999 from IMPROVE, California Air Resources Board dichotomous sampler monitoring (CARB dichot), the Southern California Children’s Health Study (CHS), and the Inhalable Particulate Network (IPN). Results: The PM2.5 prediction model performed well across three trend estimation approaches when validated using IMPROVE and CHS data (R2= 0.84–0.91). Model performance using CARB dichot and IPN data was worse than those in IMPROVE most likely due to inconsistent sampling methods and smaller numbers of monitoring sites. Discussion: Our prediction modeling approach will allow health effects estimation associated with long-term exposures to PM2.5 over extended time periods of up to 30 years

A climate-conditioned catastrophe risk model for UK flooding

We present a transparent and validated climate-conditioned catastrophe flood model for the UK, that simulates pluvial, fluvial and coastal flood risks at 1 arcsec spatial resolution (∼ 20–25 m). Hazard layers for 10 different return periods are produced over the whole UK for historic, 2020, 2030, 2050 and 2070 conditions using the UK Climate Projections 2018 (UKCP18) climate simulations. From these, monetary losses are computed for five specific global warming levels above pre-industrial values (0.6, 1.1, 1.8, 2.5 and 3.3 ∘C). The analysis contains a greater level of detail and nuance compared to previous work, and represents our current best understanding of the UK's changing flood risk landscape. Validation against historical national return period flood maps yielded critical success index values of 0.65 and 0.76 for England and Wales, respectively, and maximum water levels for the Carlisle 2005 flood were replicated to a root mean square error (RMSE) of 0.41 m without calibration. This level of skill is similar to local modelling with site-specific data. Expected annual damage in 2020 was GBP 730 million, which compares favourably to the observed value of GBP 714 million reported by the Association of British Insurers. Previous UK flood loss estimates based on government data are ∼ 3× higher, and lie well outside our modelled loss distribution, which is plausibly centred on the observations. We estimate that UK 1 % annual probability flood losses were ∼ 6 % greater for the average climate conditions of 2020 (∼ 1.1 ∘C of warming) compared to those of 1990 (∼ 0.6 ∘C of warming), and this increase can be kept to around ∼ 8 % if all countries' COP26 2030 carbon emission reduction pledges and “net zero” commitments are implemented in full. Implementing only the COP26 pledges increases UK 1 % annual probability flood losses by 23 % above average 1990 values, and potentially 37 % in a “worst case” scenario where carbon reduction targets are missed and climate sensitivity is high.</p

A Flexible Spatio-Temporal Model for Air Pollution: Allowing for Spatio-Temporal Covariates

Given the increasing interest in the association between exposure to air pollution and adverse health outcomes, the development of models that provide accurate spatio-temporal predictions of air pollution concentrations at small spatial scales is of great importance when assessing potential health effects of air pollution. The methodology presented here has been developed as part of the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air), a prospective cohort study funded by the US EPA to investigate the relationship between chronic exposure to air pollution and cardiovascular disease. We present a spatio-temporal framework that models and predicts ambient air pollution by combining data from several different monitoring networks with the output from deterministic air pollution model(s). The model can accommodate arbitrarily missing observations and allows for a complex spatio-temporal correlation structure. We apply the model to predict long-term average concentrations of gaseous oxides of nitrogen (NOx) ─ one of the primary pollutants of interest in the MESA Air study ─ during a ten year period in the Los Angeles area, based on measurements from the EPA Air Quality System and MESA Air monitoring. The measurements are augmented by a spatio-temporal covariate based on the output from a source dispersion model for traffic related air pollution (Caline3QHC) and the model is evaluated using cross-validation. The predictive ability of the model is good with cross-validated R2 of approximately 0.7 at subject sites. The incorporation of a dispersion model output into the overall prediction model was feasible, but the particular implementation of Caline3QHC used here did not improve predictions in a model that also includes road information. However, excluding the road information the inclusion of model output improves predictions and we find some evidence that the source dispersion model can replace road covariates. The model presented in this paper has been implemented in an R package, SpatioTemporal, which will be available on CRAN shortly