Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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 1,351 records · Page 75

Action-based dynamical modeling for the Milky Way disk

We present RoadMapping, a full-likelihood dynamical modeling machinery that aims to recover the Milky Way’s (MW) gravitational potential from large samples of stars in the Galactic disk. RoadMapping models the observed positions and velocities of stars with a parameterized, three-integral distribution function (DF) in a parameterized axisymmetric potential. We investigate through differential test cases with idealized mock data how the breakdown of model assumptions and data properties affect constraints on the potential and DF. Our key results are: (i) If the MW’s true potential is not included in the assumed model potential family, we can—in the axisymmetric case—still find a robust estimate for the potential, with only ≲10% difference in surface density within |z|⩽1.1 kpc inside the observed volume. (ii) Modest systematic differences between the true and model DF are inconsequential. E.g., when binning stars to define sub-populations with simple DFs, binning errors do not affect the modeling as long as the DF parameters of neighboring bins differ by <20%. In addition, RoadMapping ensures unbiased potential estimates for either (iii) small misjudgements of the spatial selection function (i.e., ≲15% at the survey volume’s edge), (iv) if distances are known to within 10%, or (v) if proper motion uncertainties are known within 10% or are smaller than δμ≲1 mas yr{sup −1}. Challenges are the rapidly increasing computational costs for large sample sizes. Overall, RoadMapping is well suited to making precise new measurements of the MW’s potential with data from the upcoming Gaia releases.

79 ASTRONOMY AND ASTROPHYSICS↗

Approximating Photo- z PDFs for Large Surveys

Modern galaxy surveys produce redshift probability density functions (PDFs) in addition to traditional photometric redshift (photo-z) point estimates. However, the storage of photo-z PDFs may present a challenge with increasingly large catalogs, as we face a trade-off between the accuracy of subsequent science measurements and the limitation of finite storage resources. This work presents qp, a Python package for manipulating parameterizations of one-dimensional PDFs, as suitable for photo-z PDF compression. We use qp to investigate the performance of three simple PDF storage formats (quantiles, samples, and step functions) as a function of the number of stored parameters on two realistic mock data sets, representative of upcoming surveys with different data qualities. We propose some best practices for choosing a photo-z PDF approximation scheme and demonstrate the approach on a science case using performance metrics on both ensembles of individual photo-z PDFs and an estimator of the overall redshift distribution function. We show that both the properties of the set of PDFs we wish to approximate and the fidelity metric(s) chosen affect the optimal parameterization. Additionally, we find that quantiles and samples outperform step functions, and we encourage further consideration of these formats for PDF approximation.

79 ASTRONOMY AND ASTROPHYSICS↗

EQUATION OF STATE FOR NUCLEONIC AND HYPERONIC NEUTRON STARS WITH MASS AND RADIUS CONSTRAINTS

We obtain a new equation of state for the nucleonic and hyperonic inner core of neutron stars that fulfils the 2 M {sub ⊙} observations as well as the recent determinations of stellar radii below 13 km. The nucleonic equation of state is obtained from a new parameterization of the FSU2 relativistic mean-field functional that satisfies these latest astrophysical constraints and, at the same time, reproduces the properties of nuclear matter and finite nuclei while fulfilling the restrictions on high-density matter deduced from heavy-ion collisions. On the one hand, the equation of state of neutron star matter is softened around saturation density, which increases the compactness of canonical neutron stars leading to stellar radii below 13 km. On the other hand, the equation of state is stiff enough at higher densities to fulfil the 2 M {sub ⊙} limit. By a slight modification of the parameterization, we also find that the constraints of 2 M {sub ⊙} neutron stars with radii around 13 km are satisfied when hyperons are considered. The inclusion of the high magnetic fields present in magnetars further stiffens the equation of state. Hyperonic magnetars with magnetic fields in the surface of ∼10{sup 15} G and with values of ∼10{sup 18} G in the interior can reach maximum masses of 2 M {sub ⊙} with radii in the 12–13 km range.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring Alignments between Galaxies and the Cosmic Web at z ∼ 2–3 Using IGM Tomography

Many galaxy formation models predict alignments between galaxy spin and the cosmic web (i.e., directions of filaments and sheets), leading to an intrinsic alignment between galaxies that creates a systematic error in weak-lensing measurements. These effects are often predicted to be stronger at high redshifts (z ≳ 1) that are inaccessible to massive galaxy surveys on foreseeable instrumentation, but IGM tomography of the Lyα forest from closely spaced quasars and galaxies is starting to measure the z ∼ 2–3 cosmic web with requisite fidelity. Using mock surveys from hydrodynamical simulations, we examine the utility of this technique, in conjunction with coeval galaxy samples, to measure alignment between galaxies and the cosmic web at z ∼ 2.5. We show that IGM tomography surveys with ≲5 h {sup −1} Mpc sightline spacing can accurately recover the eigenvectors of the tidal tensor, which we use to define the directions of the cosmic web. For galaxy spins and shapes, we use a model parameterized by the alignment strength, Δ〈cos θ〉, with respect to the tidal tensor eigenvectors from the underlying density field, and also consider observational effects such as errors in the galaxy position angle, inclination, and redshift. Measurements using the upcoming ∼1 deg{sup 2} CLAMATO tomographic survey and 600 coeval zCOSMOS-Deep galaxies should place 3σ limits on extreme alignment models with Δ〈cos θ〉∼0.1, but much larger surveys encompassing >10,000 galaxies, such as Subaru PFS, will be required to constrain models with Δ〈cos θ〉∼0.03. These measurements will constrain models of galaxy–cosmic web alignment and test tidal torque theory at z ∼ 2, improving our understanding of the physics of intrinsic alignments.

79 ASTRONOMY AND ASTROPHYSICS↗

Limb Darkening and Planetary Transits: Testing Center-to-limb Intensity Variations and Limb-darkening Directly from Model Stellar Atmospheres

The transit method, employed by Microvariability and Oscillation of Stars ( MOST ), Kepler , and various ground-based surveys has enabled the characterization of extrasolar planets to unprecedented precision. These results are precise enough to begin to measure planet atmosphere composition, planetary oblateness, starspots, and other phenomena at the level of a few hundred parts per million. However, these results depend on our understanding of stellar limb darkening, that is, the intensity distribution across the stellar disk that is sequentially blocked as the planet transits. Typically, stellar limb darkening is assumed to be a simple parameterization with two coefficients that are derived from stellar atmosphere models or fit directly. In this work, we revisit this assumption and compute synthetic planetary-transit light curves directly from model stellar atmosphere center-to-limb intensity variations (CLIVs) using the plane-parallel Atlas and spherically symmetric SAtlas codes. We compare these light curves to those constructed using best-fit limb-darkening parameterizations. We find that adopting parametric stellar limb-darkening laws leads to systematic differences from the more geometrically realistic model stellar atmosphere CLIV of about 50–100 ppm at the transit center and up to 300 ppm at ingress/egress. While these errors are small, they are systematic, and they appear to limit the precision necessary to measure secondary effects. Our results may also have a significant impact on transit spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Toward Rapid Transient Identification and Characterization of Kilonovae

With the increasing sensitivity of advanced gravitational-wave (GW) detectors, the first joint detection of an electromagnetic and GW signal from a compact binary merger will hopefully happen within this decade. However, current GW likelihood sky areas span ∼100--1000 deg{sup 2}, and thus it is a challenging task to identify which, if any, transient corresponds to the GW event. In this study, we make a comparison between recent kilonova/macronova light-curve models for the purpose of assessing potential light-curve templates for counterpart identification. We show that recent analytical and parameterized models for these counterparts result in qualitative agreement with more complicated radiative transfer simulations. Our analysis suggests that with improved light-curve models with smaller uncertainties it will become possible to extract information about ejecta properties and binary parameters directly from the light-curve measurement. Even tighter constraints are obtained in cases for which GW and kilonova parameter estimation results are combined. It will therefore be important to make comparisons and potentially combine parameter estimation with the kilonova and GW results. However, to be prepared for upcoming detections, more realistic kilonova models are needed. These will require numerical relativity with more detailed microphysics, better radiative transfer simulations, and a better understanding of the underlying nuclear physics.

79 ASTRONOMY AND ASTROPHYSICS↗

The Twins Embedding of Type Ia Supernovae. I. The Diversity of Spectra at Maximum Light

In this work, we study the spectral diversity of Type Ia supernovae (SNe Ia) at maximum light using high signal-to-noise spectrophotometry of 173 SNe Ia from the Nearby Supernova Factory. We decompose the diversity of these spectra into different extrinsic and intrinsic components, and we construct a nonlinear parameterization of the intrinsic diversity of SNe Ia that preserves pairings of "twin"SNe Ia. We call this parameterization the "Twins Embedding."Our methodology naturally handles highly nonlinear variability in spectra, such as changes in the photosphere expansion velocity, and uses the full spectrum rather than being limited to specific spectral line strengths, ratios, or velocities. We find that the time evolution of SNe Ia near maximum light is remarkably similar, with 84.6% of the variance in common to all SNe Ia. After correcting for brightness and color, the intrinsic variability of SNe Ia is mostly restricted to specific spectral lines, and we find intrinsic dispersions as low as ~0.02 mag between 6600 and 7200 Å. With a nonlinear three-dimensional model plus one dimension for color, we can explain 89.2% of the intrinsic diversity in our sample of SNe Ia, which includes several different kinds of "peculiar"SNe Ia. A linear model requires seven dimensions to explain a comparable fraction of the intrinsic diversity. We show how a wide range of previously established indicators of diversity in SNe Ia can be recovered from the Twins Embedding. In a companion article, we discuss how these results can be applied to the standardization of SNe Ia for cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Electromagnetic Signatures from the Tidal Tail of a Black Hole—Neutron Star Merger

Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of ~3 with orientation while the peak in the optical bands varies by ~3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the V-band peak magnitude to ~-14 despite the lanthanide-rich composition.

79 ASTRONOMY AND ASTROPHYSICS↗

A Probabilistic Autoencoder for Type Ia Supernova Spectral Time Series

We construct a physically parameterized probabilistic autoencoder (PAE) to learn the intrinsic diversity of Type Ia supernovae (SNe Ia) from a sparse set of spectral time series. The PAE is a two-stage generative model, composed of an autoencoder that is interpreted probabilistically after training using a normalizing flow. We demonstrate that the PAE learns a low-dimensional latent space that captures the nonlinear range of features that exists within the population and can accurately model the spectral evolution of SNe Ia across the full range of wavelength and observation times directly from the data. By introducing a correlation penalty term and multistage training setup alongside our physically parameterized network, we show that intrinsic and extrinsic modes of variability can be separated during training, removing the need for the additional models to perform magnitude standardization. We then use our PAE in a number of downstream tasks on SNe Ia for increasingly precise cosmological analyses, including the automatic detection of SN outliers, the generation of samples consistent with the data distribution, and solving the inverse problem in the presence of noisy and incomplete data to constrain cosmological distance measurements. We find that the optimal number of intrinsic model parameters appears to be three, in line with previous studies, and show that we can standardize our test sample of SNe Ia with an rms of 0.091 ± 0.010 mag, which corresponds to 0.074 ± 0.010 mag if peculiar velocity contributions are removed.

79 ASTRONOMY AND ASTROPHYSICS↗

New Covariant Density Functionals of Nuclear Matter for Compact Star Simulations

We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of the not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well.

79 ASTRONOMY AND ASTROPHYSICS↗

A Full Accounting of the Visible Mass in SDSS MaNGA Disk Galaxies

We present a study of the ratio of visible mass to total mass in spiral galaxies to better understand the relative amount of dark matter present in galaxies of different masses and evolutionary stages. Using the velocities of the Hα emission line measured in spectroscopic observations from the Sloan Digital Sky Survey (SDSS) MaNGA Data Release 17 (DR 17), we evaluate the rotational velocity of over 5500 disk galaxies at their 90% elliptical Petrosian radii, R 90 . We compare this to the velocity expected from the total visible mass, which we compute from the stellar, H i, H 2 , and heavy metals and dust masses. H 2 mass measurements are available for only a small subset of galaxies observed in SDSS MaNGA DR17, so we derive a parameterization of the H 2 mass as a function of absolute magnitude in the r band using galaxies observed as part of SDSS DR7. With these parameterizations, we calculate the fraction of visible mass within R 90 that corresponds to the observed velocity. Based on statistically analyzing the likelihood of this fraction, we conclude that the null hypothesis (no dark matter) cannot be excluded at a confidence level better than 95% within the visible extent of the disk galaxies. We also find that when all mass components are included, the ratio of visible to total mass within the visible extent of star-forming disk galaxies increases with galaxy luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗

Real-time Recovery Efficiencies and Performance of the Palomar Transient Factory's Transient Discovery Pipeline

In this paper, we present the transient source detection efficiencies of the Palomar Transient Factory (PTF), parameterizing the number of transients that PTF found versus the number of similar transients that occurred over the same period in the survey search area but were missed. PTF was an optical sky survey carried out with the Palomar 48 inch telescope over 2009–2012, observing more than 8000 square degrees of sky with cadences of between one and five days, locating around 50,000 non-moving transient sources, and spectroscopically confirming around 1900 supernovae. We assess the effectiveness with which PTF detected transient sources, by inserting $\simeq 7$ million artificial point sources into real PTF data. We then study the efficiency with which the PTF real-time pipeline recovered these sources as a function of the source magnitude, host galaxy surface brightness, and various observing conditions (using proxies for seeing, sky brightness, and transparency). The product of this study is a multi-dimensional recovery efficiency grid appropriate for the range of observing conditions that PTF experienced and that can then be used for studies of the rates, environments, and luminosity functions of different transient types using detailed Monte Carlo simulations. Finally, we illustrate the technique using the observationally well-understood class of type Ia supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

A Compound Poisson Generator Approach to Point-source Inference in Astrophysics

Abstract The identification and description of point sources is one of the oldest problems in astronomy, yet even today the correct statistical treatment for point sources remains one of the field’s hardest problems. For dim or crowded sources, likelihood-based inference methods are required to estimate the uncertainty on the characteristics of the source population. In this work, a new parametric likelihood is constructed for this problem using compound Poisson generator (CPG) functionals that incorporate instrumental effects from first principles. We demonstrate that the CPG approach exhibits a number of advantages over non-Poissonian template fitting (NPTF)—an existing method—in a series of test scenarios in the context of X-ray astronomy. These demonstrations show that the effect of the point-spread function, effective area, and choice of point-source spatial distribution cannot, generally, be factorized as they are in NPTF, while the new CPG construction is validated in these scenarios. Separately, an examination of the diffuse-flux emission limit is used to show that most simple choices of priors on the standard parameterization of the population model can result in unexpected biases: when a model comprising both a point-source population and diffuse component is applied to this limit, nearly all observed flux will be assigned to either the population or to the diffuse component. A new parameterization is presented for these priors that properly estimates the uncertainties in this limit. In this choice of priors, CPG correctly identifies that the fraction of flux assigned to the population model cannot be constrained by the data.

79 ASTRONOMY AND ASTROPHYSICS↗

Sensitivity Studies of Dust Ice Nuclei Effect on Cirrus Clouds with the Community Atmosphere Model CAM5

In this study the effect of dust aerosol on upper tropospheric cirrus clouds through heterogeneous ice nucleation is investigated in the Community Atmospheric Model version 5 (CAM5) with two ice nucleation parameterizations. Both parameterizations consider homogeneous and heterogeneous nucleation and the competition between the two mechanisms in cirrus clouds, but differ significantly in the number concentration of heterogeneous ice nuclei (IN) from dust. Heterogeneous nucleation on dust aerosol reduces the occurrence frequency of homogeneous nucleation and thus the ice crystal number concentration in the Northern Hemisphere (NH) cirrus clouds compared to simulations with pure homogeneous nucleation. Global and annual mean shortwave and longwave cloud forcing are reduced by up to 2.0 ± 0.1 W m-2 (1σ uncertainty) and 2.4 ± 0.1 W m-2, respectively due to the presence of dust IN, with the net cloud forcing change of -0.40 ± 0.20 W m-2. Comparison of model simulations with in situ aircraft data obtained in NH mid-latitudes suggests that homogeneous ice nucleation may play an important role in the ice nucleation at these regions with temperatures of 205–230 K. However, simulations overestimate observed ice crystal number concentrations in the tropical tropopause regions with temperatures of 190–205 K, and overestimate the frequency of occurrence of high ice crystal number concentration (> 200 L-1) and underestimate the frequency of low ice crystal number concentration (< 30 L-1) at NH mid-latitudes. These results highlight the importance of quantifying the number concentrations and properties of heterogeneous IN (including dust aerosol) in the upper troposphere from the global perspective.

54 ENVIRONMENTAL SCIENCES↗

The AeroCom evaluation and intercomparison of organic aerosol in global models

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a –1 (range 34–144 Tg a &minus;1 ) and the median SOA source strength (natural and anthropogenic) is 19 Tg a –1 (range 13–121 Tg a &minus;1 ). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a –1 (range 16–121 Tg a &minus;1 ), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a –1 ; range 13–20 Tg a –1 , with one model at 37 Tg a &minus;1 ). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a –1 (range 28–209 Tg a &minus;1 ), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to –0.62 (–0.51) based on the comparison against OC (OA) urban data of all models at the surface, –0.15 (+0.51) when compared with remote measurements, and –0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. As a result, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

54 ENVIRONMENTAL SCIENCES↗

Understanding sources of organic aerosol during CalNex-2010 using the CMAQ-VBS

Abstract. Community Multiscale Air Quality (CMAQ) model simulations utilizing the traditional organic aerosol (OA) treatment (CMAQ-AE6) and a volatility basis set (VBS) treatment for OA (CMAQ-VBS) were evaluated against measurements collected at routine monitoring networks (Chemical Speciation Network (CSN) and Interagency Monitoring of Protected Visual Environments (IMPROVE)) and those collected during the 2010 California at the Nexus of Air Quality and Climate Change (CalNex) field campaign to examine important sources of OA in southern California. Traditionally, CMAQ treats primary organic aerosol (POA) as nonvolatile and uses a two-product framework to represent secondary organic aerosol (SOA) formation. CMAQ-VBS instead treats POA as semivolatile and lumps OA using volatility bins spaced an order of magnitude apart. The CMAQ-VBS approach underpredicted organic carbon (OC) at IMPROVE and CSN sites to a greater degree than CMAQ-AE6 due to the semivolatile POA treatment. However, comparisons to aerosol mass spectrometer (AMS) measurements collected at Pasadena, CA, indicated that CMAQ-VBS better represented the diurnal profile and primary/secondary split of OA. CMAQ-VBS SOA underpredicted the average measured AMS oxygenated organic aerosol (OOA, a surrogate for SOA) concentration by a factor of 5.2, representing a considerable improvement to CMAQ-AE6 SOA predictions (factor of 24 lower than AMS). We use two new methods, one based on species ratios (SOA/ΔCO and SOA/Ox) and another on a simplified SOA parameterization, to apportion the SOA underprediction for CMAQ-VBS to slow photochemical oxidation (estimated as 1.5 × lower than observed at Pasadena using −log(NOx : NOy)), low intrinsic SOA formation efficiency (low by 1.6 to 2 × for Pasadena), and low emissions or excessive dispersion for the Pasadena site (estimated to be 1.6 to 2.3 × too low/excessive). The first and third factors are common to CMAQ-AE6, while the intrinsic SOA formation efficiency for that model is estimated to be too low by about 7 × . From source-apportioned model results, we found most of the CMAQ-VBS modeled POA at the Pasadena CalNex site was attributable to meat cooking emissions (48 %, consistent with a substantial fraction of cooking OA in the observations). This is compared to 18 % from gasoline vehicle emissions, 13 % from biomass burning (in the form of residential wood combustion), and 8 % from diesel vehicle emissions. All "other" inventoried emission sources (e.g., industrial, point, and area sources) comprised the final 13 %. The CMAQ-VBS semivolatile POA treatment underpredicted AMS hydrocarbon-like OA (HOA) + cooking-influenced OA (CIOA) at Pasadena by a factor of 1.8 compared to a factor of 1.4 overprediction of POA in CMAQ-AE6, but it did capture the AMS diurnal profile of HOA and CIOA well, with the exception of the midday peak. Overall, the CMAQ-VBS with its semivolatile treatment of POA, SOA from intermediate volatility organic compounds (IVOCs), and aging of SOA improves SOA model performance (though SOA formation efficiency is still 1.6–2 × too low). However, continued efforts are needed to better understand assumptions in the parameterization (e.g., SOA aging) and provide additional certainty to how best to apply existing emission inventories in a framework that treats POA as semivolatile, which currently degrades existing model performance at routine monitoring networks. The VBS and other approaches (e.g., AE6) require additional work to appropriately incorporate IVOC emissions and subsequent SOA formation.

54 ENVIRONMENTAL SCIENCES↗

On the representation of aerosol activation and its influence on model-derived estimates of the aerosol indirect effect

Abstract. Interactions between aerosol particles and clouds contribute a great deal of uncertainty to the scientific community's understanding of anthropogenic climate forcing. Aerosol particles serve as the nucleation sites for cloud droplets, establishing a direct linkage between anthropogenic particulate emissions and clouds in the climate system. To resolve this linkage, the community has developed parameterizations of aerosol activation which can be used in global climate models to interactively predict cloud droplet number concentrations (CDNCs). However, different activation schemes can exhibit different sensitivities to aerosol perturbations in different meteorological or pollution regimes. To assess the impact these different sensitivities have on climate forcing, we have coupled three different core activation schemes and variants with the CESM-MARC (two-Moment, Multi-Modal, Mixing-state-resolving Aerosol model for Research of Climate (MARC) coupled with the National Center for Atmospheric Research's (NCAR) Community Earth System Model (CESM; version 1.2)). Although the model produces a reasonable present-day CDNC climatology when compared with observations regardless of the scheme used, ΔCDNCs between the present and preindustrial era regionally increase by over 100 % in zonal mean when using the most sensitive parameterization. These differences in activation sensitivity may lead to a different evolution of the model meteorology, and ultimately to a spread of over 0.8 W m −2 in global average shortwave indirect effect (AIE) diagnosed from the model, a range which is as large as the inter-model spread from the AeroCom intercomparison. Model-derived AIE strongly scales with the simulated preindustrial CDNC burden, and those models with the greatest preindustrial CDNC tend to have the smallest AIE, regardless of their ΔCDNC. This suggests that present-day evaluations of aerosol-climate models may not provide useful constraints on the magnitude of the AIE, which will arise from differences in model estimates of the preindustrial aerosol and cloud climatology.

54 ENVIRONMENTAL SCIENCES↗

A global model–measurement evaluation of particle light scattering coefficients at elevated relative humidity

The uptake of water by atmospheric aerosols has a pronounced effect on particle light scattering properties, which in turn are strongly dependent on the ambient relative humidity (RH). Earth system models need to account for the aerosol water uptake and its influence on light scattering in order to properly capture the overall radiative effects of aerosols. Here we present a comprehensive model–measurement evaluation of the particle light scattering enhancement factor f(RH), defined as the particle light scattering coefficient at elevated RH (here set to 85 %) divided by its dry value. The comparison uses simulations from 10 Earth system models and a global dataset of surface-based in situ measurements. In general, we find a large diversity in the magnitude of predicted f(RH) amongst the different models, which can not be explained by the site types. Based on our evaluation of sea salt scattering enhancement and simulated organic mass fraction, there is a strong indication that differences in the model parameterizations of hygroscopicity and model chemistry are driving at least some of the observed diversity in simulated f(RH). Additionally, a key point is that defining dry conditions is difficult from an observational point of view and, depending on the aerosol, may influence the measured f(RH). The definition of dry also impacts our model evaluation, because several models exhibit significant water uptake between RH = 0 % and 40 %. The multisite average ratio between model outputs and measurements is 1.64 when RH = 0 % is assumed as the model dry RH and 1.16 when RH = 40 % is the model dry RH value. The overestimation by the models is believed to originate from the hygroscopicity parameterizations at the lower RH range which may not implement all phenomena taking place (i.e., not fully dried particles and hysteresis effects). This will be particularly relevant when a location is dominated by a deliquescent aerosol such as sea salt. Our results emphasize the need to consider the measurement conditions in such comparisons and recognize that measurements referred to as dry may not be dry in model terms. Recommendations for future model–measurement evaluation and model improvements are provided.

54 ENVIRONMENTAL SCIENCES↗