Aerosol-cloud-precipitation interactions in marine stratocumulus clouds

Marine stratocumulus clouds are generally optically thick and shallow, exerting a net cooling influence on climate. Changes in atmospheric aerosol levels alter cloud microphysics (e.g., droplet size) and cloud macrophysics (e.g., liquid water path, cloud thickness), thereby affecting cloud albedo and Earth’s radiative balance. To understand the aerosol-cloud-precipitation interactions and to explore the dynamical effects, three-dimensional large-eddy simulations (LES) with detailed bin-resolved microphysics are performed to explore the diurnal variation of marine stratocumulus clouds under different aerosol levels and environmental conditions. It is shown that the marine stratocumulus cloud albedo is sensitive to aerosol perturbation under clean background conditions, and to environmental conditions such as large-scale divergence rate and free tropospheric humidity.

Based on the in-situ Eastern Pacific Emitted Aerosol Cloud Experiment (E-PEACE) during Jul. and Aug. 2011, and A-Train satellite observation of 589 individual ship tracks during Jun. 2006-Dec. 2009, an analysis of cloud albedo responses in ship tracks is presented. It is found that the albedo response in ship tracks depends on the mesoscale cloud structure, the free tropospheric humidity, and cloud top height. Under closed cell structure (i.e., cloud cells ringed by a perimeter of clear air), with sufficiently dry air above cloud tops and/or higher cloud top heights, the cloud albedo can become lower in ship tracks. Based on the satellite data, nearly 25% of ship tracks exhibited a decreased albedo. The cloud macrophysical responses are crucial in determining both the strength and the sign of the cloud albedo response to aerosols.

To understand the aerosol indirect effects on global marine warm clouds, multisensory satellite observations, including CloudSat, MODIS, CALIPSO, AMSR-E, ECMWF, CERES, and NCEP, have been applied to study the sensitivity of cloud properties to aerosol levels and to large scale environmental conditions. With an estimate of anthropogenic aerosol fraction, the global aerosol indirect radiative forcing has been assessed.

As the coupling among aerosol, cloud, precipitation, and meteorological conditions in the marine boundary layer is complex, the integration of LES modeling, in-situ aircraft measurements, and global multisensory satellite data analyses improves our understanding of this complex system.

Studies of ambient organic and inorganic aerosol in Southern California

The negative impacts of ambient aerosol particles, or particulate matter (PM), on human health and climate are well recognized. However, owing to the complexity of aerosol particle formation and chemical evolution, emissions control strategies remain difficult to develop in a cost effective manner. In this work, three studies are presented to address several key issues currently stymieing California's efforts to continue improving its air quality.

Gas-phase organic mass (GPOM) and CO emission factors are used in conjunction with measured enhancements in oxygenated organic aerosol (OOA) relative to CO to quantify the significant lack of closure between expected and observed organic aerosol concentrations attributable to fossil-fuel emissions. Two possible conclusions emerge from the analysis to yield consistency with the ambient organic data: (1) vehicular emissions are not a dominant source of anthropogenic fossil SOA in the Los Angeles Basin, or (2) the ambient SOA mass yields used to determine the SOA formation potential of vehicular emissions are substantially higher than those derived from laboratory chamber studies. Additional laboratory chamber studies confirm that, owing to vapor-phase wall loss, the SOA mass yields currently used in virtually all 3D chemical transport models are biased low by as much as a factor of 4. Furthermore, predictions from the Statistical Oxidation Model suggest that this bias could be as high as a factor of 8 if the influence of the chamber walls could be removed entirely.

Once vapor-phase wall loss has been accounted for in a new suite of laboratory chamber experiments, the SOA parameterizations within atmospheric chemical transport models should also be updated. To address the numerical challenges of implementing the next generation of SOA models in atmospheric chemical transport models, a novel mathematical framework, termed the Moment Method, is designed and presented. Assessment of the Moment Method strengths and weaknesses provide valuable insight that can guide future development of SOA modules for atmospheric CTMs.

Finally, regional inorganic aerosol formation and evolution is investigated via detailed comparison of predictions from the Community Multiscale Air Quality (CMAQ version 4.7.1) model against a suite of airborne and ground-based meteorological measurements, gas- and aerosol-phase inorganic measurements, and black carbon (BC) measurements over Southern California during the CalNex field campaign in May/June 2010. Results suggests that continuing to target sulfur emissions with the hopes of reducing ambient PM concentrations may not the most effective strategy for Southern California. Instead, targeting dairy emissions is likely to be an effective strategy for substantially reducing ammonium nitrate concentrations in the eastern part of the Los Angeles Basin.

The influence of freestream turbulence, freestream nuclei populations and a drag-reducing polymer on cavitation inception on two axisymmetric bodies

The influence upon the basic viscous flow about two axisymmetric bodies of (i) freestream turbulence level and (ii) the injection of small amounts of a drag-reducing polymer (Polyox WSR 301) into the test model boundary layer was investigated by the schlieren flow visualization technique. The changes in the type and occurrence of cavitation inception caused by the subsequent modifications in the viscous flow were studied. A nuclei counter using the holographic technique was built to monitor freestream nuclei populations and a few preliminary tests investigating the consequences of different populations on cavitation inception were carried out.

Both test models were observed to have a laminar separation over their respective test Reynolds number ranges. The separation on one test model was found to be insensitive to freestream turbulence levels of up to 3.75 percent. The second model was found to be very susceptible having its critical velocity reduced from 30 feet per second at a 0.04 percent turbulence level to 10 feet per second at a 3.75 percent turbulence level. Cavitation tests on both models at the lowest turbulence level showed the value of the incipient cavitation number and the type of cavitation were controlled by the presence of the laminar separation. Cavitation tests on the second model at 0.65 percent turbulence level showed no change in the inception index, but the appearance of the developed cavitation was altered.

The presence of Polyox in the boundary layer resulted in a cavitation suppression comparable to that found by other investigators. The elimination of the normally occurring laminar separation on these bodies by a polymer-induced instability in the laminar boundary layer was found to be responsible for the suppression of inception.

Freestream nuclei populations at test conditions were measured and it was found that if there were many freestream gas bubbles the normally present laminar separation was elminated and travelling bubble type cavitation occurred - the value of the inception index then depended upon the nuclei population. In cases where the laminar separation was present it was found that the value of the inception index was insensitive to the free stream nuclei populations.

A quantitative investigation of the solar modulation of cosmic-ray protons and helium nuclei

The differential energy spectra of cosmic-ray protons and He nuclei have been measured at energies up to 315 MeV/nucleon using balloon- and satellite-borne instruments. These spectra are presented for solar quiet times for the years 1966 through 1970. The data analysis is verified by extensive accelerator calibrations of the detector systems and by calculations and measurements of the production of secondary protons in the atmosphere.

The spectra of protons and He nuclei in this energy range are dominated by the solar modulation of the local interstellar spectra. The transport equation governing this process includes as parameters the solar-wind velocity, V, and a diffusion coefficient, K(r,R), which is assumed to be a scalar function of heliocentric radius, r, and magnetic rigidity, R. The interstellar spectra, j_D, enter as boundary conditions on the solutions to the transport equation. Solutions to the transport equation have been calculated for a broad range of assumed values for K(r,R) and j_D and have been compared with the measured spectra.

It is found that the solutions may be characterized in terms of a dimensionless parameter, ψ(r,R) = ^∞∫_r V dr'/(K(r',R). The amount of modulation is roughly proportional to ψ. At high energies or far from the Sun, where the modulation is weak, the solution is determined primarily by the value of ψ (and the interstellar spectrum) and is not sensitive to the radial dependence of the diffusion coefficient. At low energies and for small r, where the effects of adiabatic deceleration are found to be large, the spectra are largely determined by the radial dependence of the diffusion coefficient and are not very sensitive to the magnitude of ψ or to the interstellar spectra. This lack of sensitivity to j_D implies that the shape of the spectra at Earth cannot be used to determine the interstellar intensities at low energies.

Values of ψ determined from electron data were used to calculate the spectra of protons and He nuclei near Earth. Interstellar spectra of the form j_D α (W - 0.25m)^-2.65 for both protons and He nuclei were found to yield the best fits to the measured spectra for these values of ψ, where W is the total energy and m is the rest energy. A simple model for the diffusion coefficient was used in which the radial and rigidity dependence are separable and K is independent of radius inside a modulation region which has a boundary at a distance D. Good agreement was found between the measured and calculated spectra for the years 1965 through 1968, using typical boundary distances of 2.7 and 6.1 A.U. The proton spectra observed in 1969 and 1970 were flatter than in previous years. This flattening could be explained in part by an increase in D, but also seemed to require that a noticeable fraction of the observed protons at energies as high at 50 to 100 MeV be attributed to quiet-time solar emission. The turnup in the spectra at low energies observed in all years was also attributed to solar emission. The diffusion coefficient used to fit the 1965 spectra is in reasonable agreement with that determined from the power spectra of the interplanetary magnetic field (Jokipii and Coleman, 1968). We find a factor of roughly 3 increase in ψ from 1965 to 1970, corresponding to the roughly order of magnitude decrease in the proton intensity at 250 MeV. The change in ψ might be attributed to a decrease in the diffusion coefficient, or, if the diffusion coefficient is essentially unchanged over that period (Mathews et al., 1971), might be attributed to an increase in the boundary distance, D.

Secondary Organic Aerosol Composition Studies Using Mass Spectrometry

Trace volatile organic compounds emitted by biogenic and anthropogenic sources into the atmosphere can undergo extensive photooxidation to form species with lower volatility. By equilibrium partitioning or reactive uptake, these compounds can nucleate into new aerosol particles or deposit onto already-existing particles to form secondary organic aerosol (SOA). SOA and other atmospheric particulate matter have measurable effects on global climate and public health, making understanding SOA formation a needed field of scientific inquiry. SOA formation can be done in a laboratory setting, using an environmental chamber; under these controlled conditions it is possible to generate SOA from a single parent compound and study the chemical composition of the gas and particle phases. By studying the SOA composition, it is possible to gain understanding of the chemical reactions that occur in the gas phase and particle phase, and identify potential heterogeneous processes that occur at the surface of SOA particles. In this thesis, mass spectrometric methods are used to identify qualitatively and qualitatively the chemical components of SOA derived from the photooxidation of important anthropogenic volatile organic compounds that are associated with gasoline and diesel fuels and industrial activity (C12 alkanes, toluene, and o-, m-, and p-cresols). The conditions under which SOA was generated in each system were varied to explore the effect of NO_x and inorganic seed composition on SOA chemical composition. The structure of the parent alkane was varied to investigate the effect on the functionalization and fragmentation of the resulting oxidation products. Relative humidity was varied in the alkane system as well to measure the effect of increased particle-phase water on condensed-phase reactions. In all systems, oligomeric species, resulting potentially from particle-phase and heterogeneous processes, were identified. Imines produced by reactions between (NH₄)₂SO₄ seed and carbonyl compounds were identified in all systems. Multigenerational photochemistry producing low- and extremely low-volatility organic compounds (LVOC and ELVOC) was reflected strongly in the particle-phase composition as well.

Investigation of Fundamental Processes Governing Secondary Organic Aerosol Formation in Laboratory Chambers

Our understanding of the processes and mechanisms by which secondary organic aerosol (SOA) is formed is derived from laboratory chamber studies. In the atmosphere, SOA formation is primarily driven by progressive photooxidation of SOA precursors, coupled with their gas-particle partitioning. In the chamber environment, SOA-forming vapors undergo multiple chemical and physical processes that involve production and removal via gas-phase reactions; partitioning onto suspended particles vs. particles deposited on the chamber wall; and direct deposition on the chamber wall. The main focus of this dissertation is to characterize the interactions of organic vapors with suspended particles and the chamber wall and explore how these intertwined processes in laboratory chambers govern SOA formation and evolution.

A Functional Group Oxidation Model (FGOM) that represents SOA formation and evolution in terms of the competition between functionalization and fragmentation, the extent of oxygen atom addition, and the change of volatility, is developed. The FGOM contains a set of parameters that are to be determined by fitting of the model to laboratory chamber data. The sensitivity of the model prediction to variation of the adjustable parameters allows one to assess the relative importance of various pathways involved in SOA formation.

A critical aspect of the environmental chamber is the presence of the wall, which can induce deposition of SOA-forming vapors and promote heterogeneous reactions. An experimental protocol and model framework are first developed to constrain the vapor-wall interactions. By optimal fitting the model predictions to the observed wall-induced decay profiles of 25 oxidized organic compounds, the dominant parameter governing the extent of wall deposition of a compound is identified, i.e., wall accommodation coefficient. By correlating this parameter with the molecular properties of a compound via its volatility, the wall-induced deposition rate of an organic compound can be predicted based on its carbon and oxygen numbers in the molecule.

Heterogeneous transformation of δ-hydroxycarbonyl, a major first-generation product from long-chain alkane photochemistry, is observed on the surface of particles and walls. The uniqueness of this reaction scheme is the production of substituted dihydrofuran, which is highly reactive towards ozone, OH, and NO3, thereby opening a reaction pathway that is not usually accessible to alkanes. A spectrum of highly-oxygenated products with carboxylic acid, ester, and ether functional groups is produced from the substituted dihydrofuran chemistry, thereby affecting the average oxidation state of the alkane-derived SOA.

The vapor wall loss correction is applied to several chamber-derived SOA systems generated from both anthropogenic and biogenic sources. Experimental and modeling approaches are employed to constrain the partitioning behavior of SOA-forming vapors onto suspended particles vs. chamber walls. It is demonstrated that deposition of SOA-forming vapors to the chamber wall during photooxidation experiments can lead to substantial and systematic underestimation of SOA. Therefore, it is likely that a lack of proper accounting for vapor wall losses that suppress chamber-derived SOA yields contribute substantially to the underprediction of ambient SOA concentrations in atmospheric models.

Mass measurements of some proton-rich nuclei

The Q values and 0^o cross sections of (He³, n) reactions forming seven proton-rich nuclei have been measured with accuracies varying from 6 to 18 keV. The Q values (in keV) are: Si²⁶ (85), S³⁰ (-573), Ar³⁴ (-759), Ti⁴² (-2865), Cr⁴⁸ (5550), Ni⁵⁶ (4513) and Zn⁶⁰ (818). At least one excited state was found for all but Ti⁴². The first four nuclei complete isotopic spin triplets; the results obtained agree well with charge-symmetry predictions. The last three, all multiples of the α particle, are important in the α and e-process theories of nucleo-synthesis in stars. The energy available for β decay of these three was found by magnetic spectrometer measurements of the (He³, p) Q values of reactions leading to V⁴⁸, Co⁵⁶, and Cu⁶⁰. Many excited states were seen: V⁴⁸ (3), Co⁵⁶ (15), Cu⁶⁰ (23). The first two states of S³⁰ are probably 0⁺ and 2⁺ from (He³, n) angular distribution measurements. Two NaI γ-ray measurements are described: the decay of Ar³⁴ (measured ^Ƭ1/2 = 1.2 ± 0.3s) and the prompt γ-ray spectrum from Fe⁵⁴(He³, nγ)Ni⁵⁶. Possible collective structure in Ni⁵⁶ and Ca⁴⁰, both doubly magic, is discussed.

The (He³, n) neutron energy and yield measurements utilized neutron-induced nuclear reactions in a silicon semiconductor detector. Cross sections for the most important detection processes, Si²⁸ (n, α) Mg²⁵ and Si²⁸ (n, p) Al²⁸, are presented for reactions leading to the first four states of both residual nuclei for neutron energies from 7.3 to 16.4 MeV. Resolution and pulse-height anomalies associated with recoil Mg²⁵ and Al²⁸ ions are discussed. The 0^o cross section for Be⁹ (α, n) C¹², used to provide calibration neutrons, has been measured with a stilbene spectrometer for n_o (5.0 ≤ E_α ≤ 12 MeV), n₁ (4.3 ≤ E_α ≤ 12.0 MeV) and n₂ (6.0 ≤ E_α ≤ 10.1 MeV). Resonances seen in the n_o yield may correspond to nine new levels in C¹³.

A study of T=1, T=3/2, and T=2 states in some light nuclei using (He^3,n) reactions

The (He³, n) reactions on B¹¹, N¹⁵, O¹⁶, and O¹⁸ targets have been studied using a pulsed-beam time-of-flight spectrometer. Special emphasis was placed upon the determination of the excitation energies and properties of states with T = 1 (in Ne¹⁸), T = 3/2 (in N¹³ and F¹⁷) and T = 2 (in Ne²⁰). The identification of the T = 3/2 and T = 2 levels is based on the structure of these states as revealed by intensities and shapes of angular distributions. The reactions are interpreted in terms of double stripping theory. Angular distributions have been compared with plane and distorted wave stripping theories. Results for the four reactions are summarized below:

1) O¹⁶ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV and two previously unreported levels in Ne¹⁸ were observed at E_x = 4.55 ± .015 MeV (Γ = 70 ± 30 keV) and E_x = 5.14 ± .018 MeV (Γ = 100 ± 40 keV).

2) B¹¹ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV. Three T = 3/2 levels in N¹³ have been identified at E_x = 15.068 ± .008 MeV (Γ ˂ 15 keV), E_x = 18.44 ± .04, and E_x 18.98 ± .02 MeV (Γ = 40 ± 20 keV).

3) N¹⁵ (He³, n). The reaction has been studied at incident energies up to 11.88 MeV. T = 3/2 levels in F¹⁷ have been identified at E_x = 11.195 ± .007 MeV (Γ ˂ 20 keV), E_x = 12.540 ± .010 MeV (Γ ˂ 25 keV), and E_x = 13.095 ± .009 MeV (Γ ˂ 25 keV).

4) O¹⁸ (He³, n). The reaction has been studied at incident energies up to 9.0 MeV. The excitation energy of the lowest T = 2 level in Ne²⁰ has been found to be 16.730 ± .006 MeV (Γ ˂ 20 keV).

Angular distributions of the transitions leading to the above higher isospin states are well described by double stripping theory. Analog correspondences are established by comparing the present results with recent studies (t, p) and (He³, p) reactions on the same targets.

Beta and gamma vibrational bands in deformed nuclei

Energies and relative intensities of gamma transitions in ¹⁵²Sm, ¹⁵²Gd, ¹⁵⁴Gd, ¹⁶⁶Er, and ²³²U following radioactive decay have been measured with a Ge(Li) spectrometer. A peak fitting program has been developed to determine gamma ray energies and relative intensities with precision sufficient to give a meaningful test of nuclear models. Several previously unobserved gamma rays were placed in the nuclear level schemes. Particular attention has been paid to transitions from the beta and gamma vibrational bands, since the gamma ray branching ratios are sensitive tests of configuration mixing in the nuclear levels. As the reduced branching ratios depend on the multipolarity of the gamma transitions, experiments were performed to measure multipole mixing ratios for transitions from the gamma vibrational band. In ¹⁵⁴Gd, angular correlation experiments showed that transitions from the gamma band to the ground state band were predominantly electric quadrupole, in agreement with the rotational model. In ²³²U, the internal conversion spectrum has been studied with a Si(Li) spectrometer constructed for electron spectroscopy. The strength of electric monopole transitions and the multipolarity of some gamma transitions have been determined from the measured relative electron intensities.

The results of the experiments have been compared with the rotational model and several microscopic models. Relative B(E2) strengths for transitions from the gamma band in ²³²U and ¹⁶⁶Er are in good agreement with a single parameter band mixing model, with values of z₂= 0.025(10) and 0.046(2), respectively. Neither the beta nor the gamma band transition strengths in ¹⁵²Sm and ¹⁵⁴Gd can be accounted for by a single parameter theory, nor can agreement be found by considering the large mixing found between the beta and gamma bands. The relative B(E2) strength for transitions from the gamma band to the beta band in ²³²U is found to be five times greater than the strength to the ground state band, indicating collective transitions with strength approximately 15 single particle units.

In Vitro Surfactant and Perfluorocarbon Aerosol Deposition in a Neonatal Physical Model of the Upper Conducting Airways

Objective: Aerosol delivery holds potential to release surfactant or perfluorocarbon (PFC) to the lungs of neonates with respiratory distress syndrome with minimal airway manipulation. Nevertheless, lung deposition in neonates tends to be very low due to extremely low lung volumes, narrow airways and high respiratory rates. In the present study, the feasibility of enhancing lung deposition by intracorporeal delivery of aerosols was investigated using a physical model of neonatal conducting airways. Methods: The main characteristics of the surfactant and PFC aerosols produced by a nebulization system, including the distal air pressure and air flow rate, liquid flow rate and mass median aerodynamic diameter (MMAD), were measured at different driving pressures (4-7 bar). Then, a three-dimensional model of the upper conducting airways of a neonate was manufactured by rapid prototyping and a deposition study was conducted. Results: The nebulization system produced relatively large amounts of aerosol ranging between 0.3 +/- 0.0 ml/min for surfactant at a driving pressure of 4 bar, and 2.0 +/- 0.1 ml/min for distilled water (H(2)Od) at 6 bar, with MMADs between 2.61 +/- 0.1 mu m for PFD at 7 bar and 10.18 +/- 0.4 mu m for FC-75 at 6 bar. The deposition study showed that for surfactant and H(2)Od aerosols, the highest percentage of the aerosolized mass (similar to 65%) was collected beyond the third generation of branching in the airway model. The use of this delivery system in combination with continuous positive airway pressure set at 5 cmH(2)O only increased total airway pressure by 1.59 cmH(2)O at the highest driving pressure (7 bar). Conclusion: This aerosol generating system has the potential to deliver relatively large amounts of surfactant and PFC beyond the third generation of branching in a neonatal airway model with minimal alteration of pre-set respiratory support.

Diseño preliminar de un sistema para generación de mezclas aerosol/vapor

El objetivo de este proyecto es establecer, mediante un prediseño, las especificaciones y requisitos que deberá cumplir el sistema para satisfacer las necesidades que tiene el Centro Nacional de Verificación de Maquinaria del Instituto Nacional de Seguridad e Higiene en el Trabajo, detectar cuáles serían las mejores soluciones para su configuración y realizar una preselección de los componentes principales para preparar un informe que permitan cuantificar el esfuerzo necesario para el diseño, adquisición de componentes, montaje y puesta en marcha del equipo. Este informe deberá recoger información de las distintas opciones para generar estas mezclas en condiciones de presión, temperatura y humedad seleccionadas y proponer un diseño preliminar que abarque la etapa de generación aerosol/vapor.