Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling variability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

Spitzer Variability Properties of Low-gravity L Dwarfs

We present Spitzer Space Telescope variability monitoring observations of three low-gravity L dwarfs with previous detections of variability in the near-IR: 2MASS J0045+16, 2MASS J0501–00, and 2MASS J1425–36. We detect significant periodic variability in two of our targets, 2MASS J0045+16 and 2MASS J0501–00. We do not detect variability in 2MASS J1425–36. Combining our new rotation periods with rotational velocities, we calculate inclination angles of 22° ± 1°, 60{sub -8}{sup +13°}, and 52{sub -13}{sup +19°} for 2MASS J0045+16, 2MASS J0501–00, and 2MASS J1425–36, respectively. Our three new objects are consistent with the tentative relations between inclination, amplitude, and color anomaly previously reported. Objects with the highest variability amplitudes are inclined equator on, while the maximum observed amplitude decreases as the inclination angle decreases. We also find a correlation between the inclination angle and (J–K){sub 2MASS} color anomaly for the sample of objects with measured inclinations. Compiling the entire sample of brown dwarfs with Spitzer variability detections, we find no enhancement in amplitude for young, early-L dwarfs compared to the field dwarf population. We find a possible enhancement in amplitude of low-gravity late-L dwarfs at 4.5 μm. We do not find a correlation between amplitude ratio and spectral type for field dwarfs or for the young population. Finally, we compile the rotation periods of a large sample of brown dwarfs with ages 1 Myr–1 Gyr and compare the rotation rates predicted by evolutionary models assuming angular momentum conservation. We find that the rotation rates of the current sample of brown dwarfs fall within the expected range set by evolutionary models and breakup limits.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of Weathering on Mechanical, Antimicrobial Properties and Biodegradation Process of Silver Loaded TPE Compounds

The incorporation of antimicrobial metals such as silver is an alternative to protect the material against microbial attack. However, loaded polymer can lose its antimicrobial properties after some time of use, and the additive may even leak out into the environment becoming harmful to non-target organisms. This study aims to evaluate the mechanical properties and antimicrobial activity of silver containing thermoplastic elastomer (TPE) samples exposed to weathering and the influence of additive incorporation in material biodegradation in the soil. For this purpose, silver ions (Ag{sup +}-bentonite, Ag{sup +}-phosphate) and silver nanoparticles (AgNp-silica) based additives were blended in a formulation of SEBS, polypropylene and mineral oil. The test samples were exposed to natural ageing over nine months, and were then evaluated according to their mechanical properties, antimicrobial activity, and degree of crystallinity and surface characteristics. The biodegradation process before and after natural ageing was evaluated through the generation of carbon dioxide. The results show that the action of natural ageing reduced the mechanical properties of loaded and unloaded TPE, and modified the degree of crystallinity and the chemical characteristic of the TPE surface. The presence or type of additive did not influence material resistance after being exposed to weathering. A decrease in antimicrobial activity in samples after natural ageing was observed. At a variable level and according to the chemical content, generation of carbon dioxide from TPE samples was greater in aged samples than in unexposed ones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of Gaussian Process Modeling to Analysis of Functional Unreliability

This paper applies Gaussian Process (GP) modeling to analysis of the functional unreliability of a “passive system.” GPs have been used widely in many ways [1]. The present application uses a GP for emulation of a system simulation code. Such an emulator can be applied in several distinct ways, discussed below. All applications illustrated in this paper have precedents in the literature; the present paper is an application of GP technology to a problem that was originally analyzed [2] using neural networks (NN), and later [3, 4] by a method called “Alternating Conditional Expectations” (ACE). This exercise enables a multifaceted comparison of both the processes and the results. Given knowledge of the range of possible values of key system variables, one could, in principle, quantify functional unreliability by sampling from their joint probability distribution, and performing a system simulation for each sample to determine whether the function succeeded for that particular setting of the variables. Using previously available system simulation codes, such an approach is generally impractical for a plant-scale problem. It has long been recognized, however, that a well-trained code emulator or surrogate could be used in a sampling process to quantify certain performance metrics, even for plant-scale problems. “Response surfaces” were used for this many years ago. But response surfaces are at their best for smoothly varying functions; in regions of parameter space where key system performance metrics may behave in complex ways, or even exhibit discontinuities, response surfaces are not the best available tool. This consideration was one of several that drove the work in [2]. In the present paper, (1) the original quantification of functional unreliability using NN [2], and later ACE [3], is reprised using GP; (2) additional information provided by the GP about uncertainty in the limit surface, generally unavailable in other representations, is discussed; (3) a simple forensic exercise is performed, analogous to the inverse problem of code calibration, but with an accident management spin: given an observation about containment pressure, what can we say about the system variables? References 1. For an introduction to GPs, see (for example) Gaussian Processes for Machine Learning, C. E. Rasmussen and C. K. I. Williams (MIT, 2006). 2. Reliability Quantification of Advanced Reactor Passive Safety Systems, J. J. Vandenkieboom, PhD Thesis (University of Michigan, 1996). 3. Z. Cui, J. C. Lee, J. J. Vandenkieboom, and R. W. Youngblood, “Unreliability Quantification of a Containment Cooling System through ACE and ANN Algorithms,” Trans. Am. Nucl. Soc. 85, 178 (2001). 4. Risk and Safety Analysis of Nuclear Systems, J. C. Lee and N. J. McCormick (Wiley, 2011). See especially §11.2.4.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A SYSTEMATIC SEARCH FOR PERIODICALLY VARYING QUASARS IN PAN-STARRS1: AN EXTENDED BASELINE TEST IN MEDIUM DEEP SURVEY FIELD MD09

We present a systematic search for periodically varying quasars and supermassive black hole binary (SMBHB) candidates in the Pan-STARRS1 (PS1) Medium Deep Survey’s MD09 field. From a color-selected sample of 670 quasars extracted from a multi-band deep-stack catalog of point sources, we locally select variable quasars and look for coherent periods with the Lomb–Scargle periodogram. Three candidates from our sample demonstrate strong variability for more than ∼3 cycles, and their PS1 light curves are well fitted to sinusoidal functions. We test the persistence of the candidates’ apparent periodic variations detected during the 4.2 years of the PS1 survey with archival photometric data from the SDSS Stripe 82 survey or new monitoring with the Large Monolithic Imager at the Discovery Channel Telescope. None of the three periodic candidates (including PSO J334.2028+1.4075) remain persistent over the extended baseline of 7–14 years, corresponding to a detection rate of <1 in 670 quasars in a search area of ≈5 deg{sup 2}. Even though SMBHBs should be a common product of the hierarchal growth of galaxies, and periodic variability in SMBHBs has been theoretically predicted, a systematic search for such signatures in a large optical survey is strongly limited by its temporal baseline and the “red noise” associated with normal quasar variability. We show that follow-up long-term monitoring (≳5 cycles) is crucial to our search for these systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Midmorning Point Sampling May Not Accurately Represent Nitrous Oxide Emissions Following Fertilizer Applications

A common approach for measuring N 2 O emissions is to collect midmorning or early evening gas samples from experiments utilizing the chamber-based flux methodology. However, due to high spatial and temporal variability, N 2 O estimates based on midmorning or early evening sampling may not provide accurate estimates of total emissions. This study determined the impact of sampling collection timing on the precision and accuracy of N 2 O emissions estimates. Nitrous oxide emissions, air and soil temperatures, and soil moisture were measured for 21 d following the application of 224 kg urea-N ha –1 on 20 Sept. 2017, 11 Oct. 2017, and 1 May 2018, at six time intervals (0130–0230, 0530–0630, 0930–1030, 1330–1430, 1730–1830, and 2130–2230 h) over a 24-h period. Based on multiple daily measurements, point samples collected between 0930 and 1030 h (midmorning) were inconsistent in their ability to predict N 2 O emissions. However, samples collected between 2130 and 2230 h (early evening) were similar to average emissions. The number of randomly collected point samples to be within 20% of the mean 80% of the time over a 21-d period ranged from 13 samples for fertilizer applied on 20 Sept. 2017 to 48 samples for fertilizer applied on 11 Oct. 2017. This research indicates that management and climatic variability affect N 2 O emissions, and that accurate sampling protocols vary across management and climates. To reduce uncertainty, N 2 O sampling protocol should be tested under conditions likely to occur and where possible, near-continuous measurement systems should be adopted.

54 ENVIRONMENTAL SCIENCES↗

Self-supervised physics-informed generative networks for phase retrieval from a single X-ray hologram

X-ray phase contrast imaging significantly improves the visualization of structures with weak or uniform absorption, broadening its applications across a wide range of scientific disciplines. Propagation-based phase contrast is particularly suitable for time- or dose-critical in vivo/in situ/operando (tomography) experiments because it requires only a single intensity measurement. However, the phase information of the wave field is lost during the measurement and must be recovered. Conventional algebraic and iterative methods often rely on specific approximations or boundary conditions that may not be met by many samples or experimental setups. In addition, they require manual tuning of reconstruction parameters by experts, making them less adaptable for complex or variable conditions. Here we present a self-learning approach for solving the inverse problem of phase retrieval in the near-field regime of Fresnel theory using a single intensity measurement (hologram). A physics-informed generative adversarial network is employed to reconstruct both the phase and absorbance of the unpropagated wave field in the sample plane from a single hologram. Unlike most state-of-the-art deep learning approaches for phase retrieval, our approach does not require paired, unpaired, or simulated training data. This significantly broadens the applicability of our approach, as acquiring or generating suitable training data remains a major challenge due to the wide variability in sample types and experimental configurations. The algorithm demonstrates robust and consistent performance across diverse imaging conditions and sample types, delivering quantitative, high-quality reconstructions for both simulated data and experimental datasets acquired at beamline P05 at PETRA III (DESY, Hamburg), operated by Helmholtz-Zentrum Hereon. Furthermore, it enables the simultaneous retrieval of both phase and absorption information.

36 MATERIALS SCIENCE↗

Diel, seasonal, and inter-annual variation in carbon dioxide effluxes from lakes and reservoirs

Abstract Accounting for temporal changes in carbon dioxide (CO 2 ) effluxes from freshwaters remains a challenge for global and regional carbon budgets. Here, we synthesize 171 site-months of flux measurements of CO 2 based on the eddy covariance method from 13 lakes and reservoirs in the Northern Hemisphere, and quantify dynamics at multiple temporal scales. We found pronounced sub-annual variability in CO 2 flux at all sites. By accounting for diel variation, only 11% of site-months were net daily sinks of CO 2 . Annual CO 2 emissions had an average of 25% (range 3%–58%) interannual variation. Similar to studies on streams, nighttime emissions regularly exceeded daytime emissions. Biophysical regulations of CO 2 flux variability were delineated through mutual information analysis. Sample analysis of CO 2 fluxes indicate the importance of continuous measurements. Better characterization of short- and long-term variability is necessary to understand and improve detection of temporal changes of CO 2 fluxes in response to natural and anthropogenic drivers. Our results indicate that existing global lake carbon budgets relying primarily on daytime measurements yield underestimates of net emissions.

54 ENVIRONMENTAL SCIENCES↗

Parameter Uncertainty for Repository Thermal Analysis

This report is one follow-on to a study of reference geologic disposal design concepts (Hardin et al. 2011a). Based on an analysis of maximum temperatures, that study concluded that certain disposal concepts would require extended decay storage prior to emplacement, or the use of small waste packages, or both. The study used nominal values for thermal properties of host geologic media and engineered materials, demonstrating the need for uncertainty analysis to support the conclusions. This report is a first step that identifies the input parameters of the maximum temperature calculation, surveys published data on measured values, uses an analytical approach to determine which parameters are most important, and performs an example sensitivity analysis. Using results from this first step, temperature calculations planned for FY12 can focus on only the important parameters, and can use the uncertainty ranges reported here. The survey of published information on thermal properties of geologic media and engineered materials, is intended to be sufficient for use in generic calculations to evaluate the feasibility of reference disposal concepts. A full compendium of literature data is beyond the scope of this report. The term “uncertainty” is used here to represent both measurement uncertainty and spatial variability, or variability across host geologic units. For the most important parameters (e.g., buffer thermal conductivity) the extent of literature data surveyed samples these different forms of uncertainty and variability. Finally, this report is intended to be one chapter or section of a larger FY12 deliverable summarizing all the work on design concepts and thermal load management for geologic disposal (M3FT-12SN0804032, due 15Aug2012).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Towards long-term standardised carbon and greenhouse gas observations for monitoring Europe’s terrestrial ecosystems: A review

Research infrastructures play a key role in launching a new generation of integrated long-term, geographically distributed observation programmes designed to monitor climate change, better understand its impacts on global ecosystems, and evaluate possible mitigation and adaptation strategies. The pan-European Integrated Carbon Observation System combines carbon and greenhouse gas (GHG; CO 2 , CH 4 , N 2 O, H 2 O) observations within the atmosphere, terrestrial ecosystems and oceans. High-precision measurements are obtained using standardised methodologies, are centrally processed and openly available in a traceable and verifiable fashion in combination with detailed metadata. The Integrated Carbon Observation System ecosystem station network aims to sample climate and land-cover variability across Europe. In addition to GHG flux measurements, a large set of complementary data (including management practices, vegetation and soil characteristics) is collected to support the interpretation, spatial upscaling and modelling of observed ecosystem carbon and GHG dynamics. The applied sampling design was developed and formulated in protocols by the scientific community, representing a trade-off between an ideal dataset and practical feasibility. Finally, the use of open-access, high-quality and multi-level data products by different user communities is crucial for the Integrated Carbon Observation System in order to achieve its scientific potential and societal value.

54 ENVIRONMENTAL SCIENCES↗

Effect of Biomass Water Dynamics in Cosmic-Ray Neutron Sensor Observations: A Long-Term Analysis of Maize–Soybean Rotation in Nebraska

Precise soil water content (SWC) measurement is crucial for effective water resource management. This study utilizes the Cosmic-Ray Neutron Sensor (CRNS) for area-averaged SWC measurements, emphasizing the need to consider all hydrogen sources, including time-variable plant biomass and water content. Near Mead, Nebraska, three field sites (CSP1, CSP2, and CSP3) growing a maize–soybean rotation were monitored for 5 (CSP1 and CSP2) and 13 (CSP3) years. Data collection included destructive biomass water equivalent (BWE) biweekly sampling, epithermal neutron counts, atmospheric meteorological variables, and point-scale SWC from a sparse time domain reflectometry (TDR) network (four locations and five depths). In 2023, dense gravimetric SWC surveys were collected eight (CSP1 and CSP2) and nine (CSP3) times over the growing season (April to October). The N0 parameter exhibited a linear relationship with BWE, suggesting that a straightforward vegetation correction factor may be suitable (fb). Results from the 2023 gravimetric surveys and long-term TDR data indicated a neutron count rate reduction of about 1% for every 1 kg m−2 (or mm of water) increase in BWE. This reduction factor aligns with existing shorter-term row crop studies but nearly doubles the value previously reported for forests. This long-term study contributes insights into the vegetation correction factor for CRNS, helping resolve a long-standing issue within the CRNS community.

Chemistry↗

Adversarial sampling of unknown and high-dimensional conditional distributions

Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximations need to be made when the distribution is unknown. In this paper the sampling method, as well as the inference of the underlying distribution, are both handled with a data-driven method known as generative adversarial networks (GAN), which trains two competing neural networks to produce a network that can effectively generate samples from the training set distribution. In practice, it is often necessary to draw samples from conditional distributions. When the conditional variables are continuous, only one (if any) data point corresponding to a particular value of a conditioning variable may be available, which is not sufficient to estimate the conditional distribution. This work handles this problem using an a priori estimation of the conditional moments of a PDF. Herein, two approaches, stochastic estimation, and an external neural network are compared for computing these moments; however, any preferred method can be used. The algorithm is demonstrated in the case of the deconvolution of a filtered turbulent flow field. It is shown that all the versions of the proposed algorithm effectively sample the target conditional distribution with minimal impact on the quality of the samples compared to state-of-the-art methods. Additionally, the procedure can be used as a metric for the diversity of samples generated by a conditional GAN (cGAN) conditioned with continuous variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Physical properties, internal structure, and the three‐dimensional petrography of CI chondrites

physical properties and the nature of their breccation, we investigated nine samples of the Ivuna and Orgueil CI chondrites ranging in size from 1 mm to 4 cm in approximate diameter. The combined mass of unique material investigated in this work is 113 g. For our investigations, we use ideal gas pycnometry, 3-D laser scanning, x-ray computed microtomography (μCT), and accompanying digital data extraction techniques. We found that the bulk density of the samples ranged from 1.61 to 2.10 g cm −3 . Larger samples tend to have a lower bulk density. Grain density (ranging from 2.44 to 2.55 g cm −3 ) is significantly less variable than the bulk density in our samples and the quantity of porosity (ranging from 14.6% to 33.8%) is the dominant factor in determining the bulk density of CI chondrite material. Our μCT results show that the visible porosity across all sizes of our CI chondrite samples is in the form of cracks, but these cracks can account for less than two-thirds of the porosity in the CI chondrites. Other porosity is not visible, even at μCT resolutions of 2.7 μm voxel edge −1 and we conclude that it is sub-micron in nature. It is not clear if the cracks seen in our samples are indigenous to the chondrites or are a result of terrestrial processes. We also find that the CI chondrites are excellent examples of the fractal-like nature of brecciation, where clasts can be observed at all scales we imaged. The breccias are composed of sub-equant-shaped and sub-rounded-textured clasts like melt-free impact breccias on other solar system bodies. From our μCT volume and digital data extraction, we determine that the Ivuna CI chondrite breccia is organized: the mostly sub-equant clasts within our ~2 cm chunk of Ivuna have a mean diameter of 1.33 mm and their aligned longest axes define a lineation structure. We speculate that the lineation was imparted after fragmentation of the clasts by slight shear on the parent asteroid which could be the result of seismic-related granular flow or mild non-axial impact-related compaction. These data will help to place returned asteroidal material from asteroids 162173 Ryugu and 101955 Bennu and the CI chondrites into a mutual geological context.

CI chondrite↗

Systematic modification of functionality in disordered elastic networks through free energy surface tailoring

A combined machine learning–physics–based approach is explored for molecular and materials engineering. Specifically, collective variables, akin to those used in enhanced sampled simulations, are constructed using a machine learning model trained on data gathered from a single system. Through the constructed collective variables, it becomes possible to identify critical molecular interactions in the considered system, the modulation of which enables a systematic tailoring of the system’s free energy landscape. To explore the efficacy of the proposed approach, we use it to engineer allosteric regulation and uniaxial strain fluctuations in a complex disordered elastic network. Its successful application in these two cases provides insights regarding how functionality is governed in systems characterized by extensive connectivity and points to its potential for design of complex molecular systems.

42 ENGINEERING↗

Continuous-time probabilistic models for longitudinal electronic health records

Analysis of longitudinal Electronic Health Record (EHR) data is an important goal for precision medicine. Difficulty in applying Machine Learning (ML) methods, either predictive or unsupervised, stems in part from the heterogeneity and irregular sampling of EHR data. Here, we present an unsupervised probabilistic model that captures nonlinear relationships between variables over continuous-time. This method works with arbitrary sampling patterns and captures the joint probability distribution between variable measurements and the time intervals between them. Inference algorithms are derived that can be used to evaluate the likelihood of future using under a trained model. As an example, we consider data from the United States Veterans Health Administration (VHA) in the areas of diabetes and depression. Likelihood ratio maps are produced showing the likelihood of risk for moderate-severe vs minimal depression as measured by the Patient Health Questionnaire-9 (PHQ-9).

59 BASIC BIOLOGICAL SCIENCES↗

Fluxes all of the time? A primer on the temporal representativeness of FLUXNET: FLUXES ALL OF THE TIME?

FLUXNET, the global network of eddy covariance flux towers, provides the largest synthesized data set of CO 2 , H 2 O, and energy fluxes. To achieve the ultimate goal of providing flux information “everywhere and all of the time,” studies have attempted to address the representativeness issue, i.e., whether measurements taken in a set of given locations and measurement periods can be extrapolated to a space- and time-explicit extent (e.g., terrestrial globe, 1982–2013 climatological baseline). This study focuses on the temporal representativeness of FLUXNET and tests whether site-specific measurement periods are sufficient to capture the natural variability of climatological and biological conditions. FLUXNET is unevenly representative across sites in terms of the measurement lengths and potentials of extrapolation in time. Similarity of driver conditions among years generally enables the extrapolation of flux information beyond measurement periods. Yet such extrapolation potentials are further constrained by site-specific variability of driver conditions. Several driver variables such as air temperature, diurnal temperature range, potential evapotranspiration, and normalized difference vegetation index had detectable trends and/or breakpoints within the baseline period, and flux measurements generally covered similar and biased conditions in those drivers. About 38% and 60% of FLUXNET sites adequately sampled the mean conditions and interannual variability of all driver conditions, respectively. For long-record sites (≥15 years) the percentages increased to 59% and 69%, respectively. However, the justification of temporal representativeness should not rely solely on the lengths of measurements. Whenever possible, site-specific consideration (e.g., trend, breakpoint, and interannual variability in drivers) should be taken into account.

58 GEOSCIENCES↗

Uncertainty Quantification of Classifying Mechanical Breach in the Crash and Burn Problem

The prediction of mechanical breach in the Crash and Burn problem is important for designing the full system to minimize the probability of loss of containment. Due to uncertainties in the impact angle, material properties, and material models, it is not possible to define an exact critical impact speed at which the system breaches. Furthermore, the cost of running a large scale sampling study to determine the empirical probability of breach is prohibitive. In this work, surrogate models from machine learning, namely logistic regression and artificial neural networks, are introduced to predict binary classification of pass versus breach from a limited set of samples. The structure and calibration of these classifiers is discussed, and a set of metrics for describing the performance of the classifiers is introduced. The classifiers are used on the UUR version of the Crash model to attempt prediction of failure probability and perform variable sensitivity analysis. Where a single sample of the computational model can take hundreds of CPU-hours, training and evaluating a classifier can take seconds or less, thus giving high predictive power in a relatively short time.

42 ENGINEERING↗

Supervised Gamma Process Poisson Factorization

This thesis develops the supervised gamma process Poisson factorization (S- GPPF) framework, a novel supervised topic model for joint modeling of count matrices and document labels. S-GPPF is fully generative and nonparametric: document labels and count matrices are modeled under a uni ed probabilistic framework and the number of latent topics is controlled automatically via a gamma process prior. The framework provides for multi-class classification of documents using a generative max-margin classifier. Several recent data augmentation techniques are leveraged to provide for exact inference using a Gibbs sampling scheme. The first portion of this thesis reviews supervised topic modeling and several key mathematical devices used in the formulation of S-GPPF. The thesis then introduces the S-GPPF generative model and derives the conditional posterior distributions of the latent variables for posterior inference via Gibbs sampling. The S-GPPF is shown to exhibit state-of-the-art performance for joint topic modeling and document classification on a dataset of conference abstracts, beating out competing supervised topic models. The unique properties of S-GPPF along with its competitive performance make it a novel contribution to supervised topic modeling.

97 MATHEMATICS AND COMPUTING↗

Land use regression models for the oxidative potential of fine particles (PM2.5) in five European areas

Highlights: • OP moderately correlated (r ~.6) with PM{sub 2.5} mass, NO{sub 2}, and elemental constituents. • OP{sup AA} and OP{sup GSH} LUR models were developed and evaluated for five areas. • Three reasonably performing (R{sup 2}: .45–.58) OP{sup AA} LUR models in hold-out validation. • OP{sup GSH} models were not robust to new coefficients derived against sub-samples of sites. Oxidative potential (OP) of particulate matter (PM) is proposed as a biologically-relevant exposure metric for studies of air pollution and health. We aimed to evaluate the spatial variability of the OP of measured PM{sub 2.5} using ascorbate (AA) and (reduced) glutathione (GSH), and develop land use regression (LUR) models to explain this spatial variability. We estimated annual average values (m{sup −3}) of OP{sup AA} and OP{sup GSH} for five areas (Basel, CH; Catalonia, ES; London-Oxford, UK (no OP{sup GSH}); the Netherlands; and Turin, IT) using PM{sub 2.5} filters. OP{sup AA} and OP{sup GSH} LUR models were developed using all monitoring sites, separately for each area and combined-areas. The same variables were then used in repeated sub-sampling of monitoring sites to test sensitivity of variable selection; new variables were offered where variables were excluded (p > .1). On average, measurements of OP{sup AA} and OP{sup GSH} were moderately correlated (maximum Pearson's maximum Pearson's R = = .7) with PM{sub 2.5} and other metrics (PM{sub 2.5}absorbance, NO{sub 2}, Cu, Fe). HOV (hold-out validation) R{sup 2} for OP{sup AA} models was .21, .58, .45, .53, and .13 for Basel, Catalonia, London-Oxford, the Netherlands and Turin respectively. For OP{sup GSH}, the only model achieving at least moderate performance was for the Netherlands (R{sup 2} = .31). Combined models for OP{sup AA} and OP{sup GSH} were largely explained by study area with weak local predictors of intra-area contrasts; we therefore do not endorse them for use in epidemiologic studies. Given the moderate correlation of OP{sup AA} with other pollutants, the three reasonably performing LUR models for OP{sup AA} could be used independently of other pollutant metrics in epidemiological studies.

54 ENVIRONMENTAL SCIENCES↗