Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plume Surface Interaction”

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.

286 records · Page 16

Thrust Stand Measurements Using Alternative Propellants in the Microwave Assisted Discharge Inductive Plasma Accelerator

Storable propellants (for example water, ammonia, and hydrazine) are attractive for deep space propulsion due to their naturally high density at ambient interplanetary conditions, which obviates the need for a cryogenic/venting system. Water in particular is attractive due to its ease of handling and availability both terrestrially and extra-terrestrially. While many storable propellants are reactive and corrosive, a propulsion scheme where the propellant is insulated from vulnerable (e.g. metallic) sections of the assembly would be well-suited to process these otherwise incompatible propellants. Pulsed inductive plasma thrusters meet this criterion because they can be operated without direct propellant-electrode interaction. During operation of these devices, electrical energy is capacitively stored and then discharged through an inductive coil creating a time-varying current in the coil that interacts with a plasma covering the face of the coil to induce a plasma current. Propellant is accelerated and expelled at a high exhaust velocity (O(10-100 km/s)) by the Lorentz body force arising from the interaction of the magnetic field and the induced plasma current. While this class of thruster mitigates the life-limiting issues associated with electrode erosion, many pulsed inductive plasma thrusters require high pulse energies to inductively ionize propellant. The Microwave Assisted Discharge Inductive Plasma Accelerator (MAD-IPA) is a pulsed inductive plasma thruster that addressees this issue by partially ionizing propellant inside a conical inductive coil before the main current pulse via an electron cyclotron resonance (ECR) discharge. The ECR plasma is produced using microwaves and a static magnetic field from a set of permanent magnets arranged to create a thin resonance region along the inner surface of the coil, restricting plasma formation, and in turn current sheet formation, to a region where the magnetic coupling between the plasma and the theta-pinch coil is high. The use of a conical theta-pinch coil also serves to provide neutral propellant containment and plasma plume focusing that is improved relative to the more common planar geometry of the Pulsed Inductive Thruster (PIT). In this paper, we describe thrust stand measurements performed to evaluate the specific impulse and thrust efficiency of the MAD-IPA for a variety of propellants. Propellants tested include both widely-used, non-reactive noble gases like argon, and rarely-used propellants such as water, hydrazine and ammonia. Dependencies of impulse data on propellant species are discussed in the context of the current sheet formation and electromagnetic plasma acceleration processes.

Hallock, Ashley K.↗

Thrust Stand Measurements of the Microwave Assisted Discharge Inductive Plasma Accelerator

Pulsed inductive plasma thrusters [1-3] are spacecraft propulsion devices in which electrical energy is capacitively stored and then discharged through an inductive coil. This type of pulsed thruster is electrodeless, with a time-varying current in the coil interacting with a plasma covering the face of the coil to induce a plasma current. Propellant is accelerated and expelled at a high exhaust velocity (O(10-100 km/s)) by the Lorentz body force arising from the interaction of the magnetic field and the induced plasma current. While this class of thruster mitigates the life-limiting issues associated with electrode erosion, pulsed inductive plasma thrusters require high pulse energies to inductively ionize propellant. The Microwave Assisted Dis- charge Inductive Plasma Accelerator (MAD-IPA), shown in Fig. 1, is a pulsed inductive plasma thruster that addressees this issue by partially ionizing propellant inside a conical inductive coil before the main current pulse via an electron cyclotron resonance (ECR) discharge. The ECR plasma is produced using microwaves and a static magnetic field from a set of permanent magnets arranged to create a thin resonance region along the inner surface of the coil, restricting plasma formation, and in turn current sheet formation, to a region where the magnetic coupling between the plasma and the theta-pinch coil is high. The use of a conical theta-pinch coil also serves to provide neutral propellant containment and plasma plume focusing that is improved relative to the more common planar geometry of the Pulsed Inductive Thruster (PIT) [1, 2]. In this paper, we describe thrust stand measurements performed to characterize the performance (specific impulse, thrust efficiency) of the MAD-IPA thruster. Impulse data are obtained at various pulse energies, mass flow rates and inductive coil geometries. Dependencies on these experimental parameters are discussed in the context of the current sheet formation and electromagnetic plasma acceleration processes.

Hallock, Ashley K.↗

Sputter ejection of matter from Io

Direct collisional interaction of magnetospheric particles, particularly 520-eV S ions, with Io, cause sputter removal of matter. It is estimated that direct sputtering of a full-disk S-containing atmosphere with an exobase at a few hundred km, can provide up to 5 x 10 to the tenth S atoms per sq cm-s. Supplies of S and O required to stabilize the torus are estimated to be from 10 to the 10th to 10 to the 12th per sq cm-s. Sputtering rates are calculated for an atmosphere containing a one percent concentration of Na and K, and are found to be large enough to supply the fluxes required to maintain the Na and K clouds. Sputtering is found to remove heavy molecules from the atmosphere, and the rate of direct sputtering of unprotected surfaces is calculated for ejections of S and Na. Atomic species on the surface are ejected at a rate proportional to the surface abundance; and plume sputtering, avalanche cascading, and ionic saltation which lead to spatial and temporal variations in the number of ejected particles are observed.

Haff, P. K.↗

Optical characterization of LDEF contaminant film

Dark brown molecular film deposits were found at numerous locations on the Long Duration Exposure Facility (LDEF) and have been documented in great detail by several investigators. The exact deposition mechanism for these deposits is as yet unknown, although direct and scattered atomic oxygen, and solar radiation interacting with materials outgassing products have all been implicated in the formation process. Specimens of the brown molecular film were taken from below the flange of the experimental tray located at position D10 on the LDEF. The tray was one of two, comprising the same experiment, the other being located on the wake facing side of the LDEF satellite at position B4. Having access to both trays, we were able to directly compare the effect that orientation with respect to the atomic oxygen flux vector had on the formation of the brown molecular film deposits. The film is thickest on surfaces facing toward the exterior, i.e. the tray corner, as can be seen by comparing the lee and wake aspects of the rivets. The patterns appear to be aligned not with the velocity vector but with the corner of the tray suggesting that flux to the surface is due to scattered atomic oxygen rather than direct ram impingement. The role of scattered flux is further supported by more faint plume patterns on the sides of the tray. The angle of these plumes is strongly aligned with the ram direction but the outline of the deposit implies that incident atoms are scattered by collisions with the edges of the opening resulting in a directed, but diffuse, flux of atomic oxygen to the surface. Spectral reflectance measurements in the 2 to 10 micron (4000 to 1000 wavenumbers) spectral range are presented for the film in the 'as deposited' condition and for the free standing film. The material was analyzed by FTIR (Fourier Transform Infrared) microspectroscopy using gold as the reference standard. The 'as deposited' specimen was on an aluminum rivet taken from beneath the tray flange while the free film was obtained by chipping some of the material from the rivet. The transmission spectrum over the 2 to 10 micron range for the free film is presented. This spectrum appears to be essentially the same as that presented by Crutcher et.al. for films formed at vent sites which faced into the ram direction and suggested to originate from urethanes and silicones used on the LDEF. Banks et. al. state that silicones, when exposed to atomic oxygen, release polymeric scission fragments which deposit on surfaces and form a glassy, dark contaminant layer upon further atomic oxygen exposure and solar irradiation.

Blakkolb, Brian K.↗

Testing of Mars Retropropulsion Concepts in the Langley Unitary Plan Wind Tunnel

Future Mars human lander missions will require using retrorockets beginning at supersonic conditions, something which has never been done before at Mars. Significant aerosciences challenges exist due to interactions between the retrorocket exhaust plumes, the freestream flow, and the vehicle that alter the aerodynamic behavior. Historical retropropulsion ground and flight data do not provide a basis for calibrating the computational uncertainties for aerodynamic interference forces and moments on proposed Mars descent vehicles, either because of insufficient data or dissimilar vehicles geometries and/or conditions. A new retropropulsion wind tunnel test was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to address some of the deficiencies of previous tests and to provide a basis for computational model comparisons on relevant geometries. The test included improvements in model design and measured data. Two sub-scale model geometries derived from full-scale vehicle concepts were tested: a blunt low lift-to-drag vehicle and a more slender geometry. The first model includes variations in nozzle parameters that were estimated to have first-order effects on the measurements: location, pointing direction, size, and exit-to-throat area ratio. Both models were instrumented with discrete steady and unsteady pressure sensors on the aerodynamic and interior flow path surfaces. The heats shields of both models also were covered with pressure sensitive paint that effectively measures a steady pressure distribution in an otherwise unsteady flow field. This pressure distribution was integrated over the surface geometry in order to estimate the steady aerodynamic force and moment coefficient contribution from the painted surfaces; this approach was used as a replacement for using a flow-through balance. This paper will cover the test objectives, model design, test matrix, and available test results for each model configuration at different freestream Mach numbers, model attitude, and non-dimensional thrust magnitude. A subsequent publication will cover uncertainty quantification of the test data and comparisons to computational fluid dynamics model results.

Powered Descent↗

Development and application of computational aerothermodynamics flowfield computer codes

Research was performed in the area of computational modeling and application of hypersonic, high-enthalpy, thermo-chemical nonequilibrium flow (Aerothermodynamics) problems. A number of computational fluid dynamic (CFD) codes were developed and applied to simulate high altitude rocket-plume, the Aeroassist Flight Experiment (AFE), hypersonic base flow for planetary probes, the single expansion ramp model (SERN) connected with the National Aerospace Plane, hypersonic drag devices, hypersonic ramp flows, ballistic range models, shock tunnel facility nozzles, transient and steady flows in the shock tunnel facility, arc-jet flows, thermochemical nonequilibrium flows around simple and complex bodies, axisymmetric ionized flows of interest to re-entry, unsteady shock induced combustion phenomena, high enthalpy pulsed facility simulations, and unsteady shock boundary layer interactions in shock tunnels. Computational modeling involved developing appropriate numerical schemes for the flows on interest and developing, applying, and validating appropriate thermochemical processes. As part of improving the accuracy of the numerical predictions, adaptive grid algorithms were explored, and a user-friendly, self-adaptive code (SAGE) was developed. Aerothermodynamic flows of interest included energy transfer due to strong radiation, and a significant level of effort was spent in developing computational codes for calculating radiation and radiation modeling. In addition, computational tools were developed and applied to predict the radiative heat flux and spectra that reach the model surface.

Venkatapathy, Ethiraj↗

Water in the Oceanic Lithosphere: Salt Lake Crater Xenoliths, Oahu, Hawaii

Water can be present in nominally anhydrous minerals of peridotites in the form of hydrogen bonded to structural oxygen. Such water in the oceanic upper mantle could have a significant effect on its physical and chemical properties. However, the water content of the MORB source has been inferred indirectly from the compositions of basalts. Direct determinations on abyssal peridotites are scarce because they have been heavily hydrothermally altered. Here we present the first water analyses of minerals from spinel peridotite xenoliths of Salt Lake Crater, Oahu, Hawaii, which are exceptionally fresh. These peridotites are thought to represent fragments of the Pacific oceanic lithosphere that was refertilized by alkalic Hawaiian melts. A few have unradiogenic Os and radiogenic Hf isotopes and may be fragments of an ancient (~2 Ga) depleted and recycled lithosphere. Water contents in olivine (Ol), orthopyroxene (Opx), and clinopyroxene (Cpx) were determined by FTIR spectrometry. Preliminary H_{2}O contents show ranges of 8-10 ppm for Ol, 151-277 ppm for Opx, and 337-603 ppm for Cpx. Reconstructed bulk rock H_{2}O contents range from 88-131 ppm overlapping estimates for the MORB source. Water contents between Ol minerals of the same xenolith are heterogeneous and individual OH infrared bands vary within a mineral with lower 3230 cm^{-1} and higher 3650-3400 cm^{-1} band heights from core to edge. This observation suggests disturbance of the hydrogen in Ol likely occurring during xenolith entrainment to the surface. Pyroxene water contents are higher than most water contents in pyroxenes from continental peridotite xenoliths and higher than those of abyssal peridotites. Cpx water contents decrease with increasing degree of depletion (e.g. increasing Fo in Ol and Cr# in spinel) consistent with an incompatible behavior of water. However Cpx water contents also show a positive correlation with LREE/HREE ratios and LREE concentrations consistent with refertilization. Opx water contents increase with increasing degree of depletion and decrease with LREE/HREE ratios which is inconsistent with the incompatible behavior of H. Calculated water contents of melts in equilibrium with Cpx or Opx range from 1.4 to 3.8 wt % which is higher than that of all Hawaiian lavas. Calculated melts in equilibrium with Cpx and Opx have variable but mostly high H_{2}O/Ce ratios (194 to 1146) consistent with those of rejuvenated stage lavas from Niihau and the South Arch volcanic field, but unlike the drier shield building stage tholeiites. Whether the high water contents recorded in Salt Lake Crater xenoliths were acquired before and/or during interaction of the oceanic lithosphere with the Hawaiian plume will be discussed.

Peslier, Anne H.↗

Surface composition of solid-rocket exhausted aluminum oxide particles

Particulate samples of aluminum oxide were collected on Teflon filters from the exhaust plume of the Space Shuttle (STS-61A, October 30, 1985) over the altitude interval 4.6-7.6 km immediately after launch. These particles were analyzed using SEM, energy-dispersive X-ray analysis, electron spectroscopy for chemical analysis, X-ray fluorescent spectroscopy, and conventional wet-chemical techniques. The samples were 0.6-1.0 percent surface-chlorided (chlorided meaning predominantly aluminum chlorides and oxychlorides, possibly including other adsorbed forms of chloride) by weight. This level of chloriding is about one-third of the amount determined previously from laboratory-prepared alumina and surface site samples of solid-rocket-produced alumina (SRPA) after both had been exposed to moist HCl vapor at temperatures down to ambient. This level is equivalent to previous laboratory results with samples exposed to moist HCl at temperatures above the boiling point of water. It is suggested that the present lower chloriding levels, determined for samples from a 'dry' Shuttle exhaust cloud, underscore the importance of a liquid water/hydrochloric acid phase in governing the extent of surface chloriding of SRPA. The reduced chloriding is not trivial with respect to potential physical/chemical modification of the SRPA particle surfaces and their corresponding interaction with the atmosphere.

Cofer, Wesley R., III↗

Optical imaging and high spatial resolution 21 cm H I observations of the peculiar galaxy NGC 2782 (Arp 215)

We have used the Very Large Array (VLA) B and C Arrays to make 21 cm H I observations of the peculiar galaxy NGC 2782 (Arp 215). These observations are complementary to previously published D Array VLA data, which revealed the presence of a long (5 min to 54 kpc) H I plume near the western side of this galaxy. We have also obtained BVRI H alpha images of the main body of this galaxy using the McDonald Observatory 30 inch telescope. The optical images of this galaxy show a strong stellar tail extending to the east, opposite the H I plume. Within the disk of NGC 2782, unsharp masking of the optical images at all 4 broadband wavelengths reveals three bright 'ripples', separated by approximately 15 sec. The light profiles across these ripples are symmetric, without a sharp outer edge. H alpha is strong at the starburst nucleus and along the northern and western sections of the inner ripple. The new higher resolution H I data show that the atomic gas is very clumpy. We have identified ten H I clumps in the long western plume, with H I masses of approximately 10(exp 8) solar mass, similar to those of dwarf galaxies, and column densities of approximately 10(exp 21) cm(exp -2) over surface areas of approximately 10 kpc(exp 2). No CO (1-0) emission has been detected from this plume, suggesting that it is material stripped from the outer edge of a disk galaxy. The H alpha peaks, in contrast, are generally not coincident with H I peaks. No H I is seen at the tip of the eastern extension. The H I distribution near this structure is ring-like rather than tail-like as in the optical data. We have detected redshifted H I absorption toward the central continuum source, indicating gas infall into the nuclear region. Using a restricted 3-body dynamical model, we have successfully reproduced the basic properties of NGC 2782 with an off-center collision between two galaxies, where a lower mass disk companion (M(sub 2)/M(sub 1) approximately 0.25) collides almost head-on with a larger disk galaxy. In this scenario, the companion becomes the eastern extension, and the long western tail is material pulled out from the larger galaxy. The lack of H I at the tip of the eastern tail may be due to stripping during the collision. Alternative formation scenarios include a merger or interaction between two disk galaxies, where the ripples may be due to a disturbed spiral density wave or galactic winds due to the central starburst.

Smith, Beverly J.↗

Axisymmetric Implementation for 3D-Based DSMC Codes

The primary objective in developing NASA s DSMC Analysis Code (DAC) was to provide a high fidelity modeling tool for 3D rarefied flows such as vacuum plume impingement and hypersonic re-entry flows [1]. The initial implementation has been expanded over time to offer other capabilities including a novel axisymmetric implementation. Because of the inherently 3D nature of DAC, this axisymmetric implementation uses a 3D Cartesian domain and 3D surfaces. Molecules are moved in all three dimensions but their movements are limited by physical walls to a small wedge centered on the plane of symmetry (Figure 1). Unfortunately, far from the axis of symmetry, the cell size in the direction perpendicular to the plane of symmetry (the Z-direction) may become large compared to the flow mean free path. This frequently results in inaccuracies in these regions of the domain. A new axisymmetric implementation is presented which aims to solve this issue by using Bird s approach for the molecular movement while preserving the 3D nature of the DAC software [2]. First, the computational domain is similar to that previously used such that a wedge must still be used to define the inflow surface and solid walls within the domain. As before molecules are created inside the inflow wedge triangles but they are now rotated back to the symmetry plane. During the move step, molecules are moved in 3D but instead of interacting with the wedge walls, the molecules are rotated back to the plane of symmetry at the end of the move step. This new implementation was tested for multiple flows over axisymmetric shapes, including a sphere, a cone, a double cone and a hollow cylinder. Comparisons to previous DSMC solutions and experiments, when available, are made.

Stewart, Benedicte↗

Characterization of Asian Dust Properties Near Source Region During ACE-Asia

Asian dust typically originates in desert areas far from polluted urban regions. During transport, dust layers can interact with anthropogenic sulfate and soot aerosols from heavily polluted urban areas. Added to the complex effects of clouds and natural marine aerosols, dust particles reaching the marine environment can have drastically different properties than those from the source. Thus, understanding the unique temporal and spatial variations of Asian aerosols is of special importance in regional-to-global climate issues such as radiative forcing, the hydrological cycle, and primary biological productivity in the mid-Pacific Ocean. During ACE-Asia campaign, we have acquired ground- based (temporal) and satellite (spatial) measurements to infer aerosol physical/optical/radiative properties, column precipitable water amount, and surface reflectivity over this region. The inclusion of flux measurements permits the determination of aerosol radiative flux in addition to measurements of loading and optical depth. At the time of the Terra/MODIS, SeaWiFS, TOMS and other satellite overpasses, these ground-based observations can provide valuable data to compare with satellite retrievals over land. In this paper, we will demonstrate new capability of the Deep Blue algorithm to track the evolution of the Asian dust storm from sources to sinks. Although there are large areas often covered by clouds in the dust season in East Asia, this algorithm is able to distinguish heavy dust from clouds over the entire regions. Examination of the retrieved daily maps of dust plumes over East Asia clearly identifies the sources contributing to the dust loading in the atmosphe~~. We have compared the satellite retrieved aerosol optical thickness to the ground-based measurements and obtained a reasonable agreement between these two. Our results also indicate that there is a large difference in the retrieved value of spectral single scattering albedo of windblown dust between different sources in East Asia.

Tsay, Si-Chee↗

MISR Decadal Observations of Mineral Dust: Property Characterization and Climate Applications

The Multi-angle Imaging SpectroRadiometer (MISR) provides a unique, independent source of data for studying dust emission and transport. MISR's multiple view angles allow the retrieval of aerosol properties over bright surfaces, and such retrievals have been shown to be sensitive to the non-sphericity of dust aerosols over both land and water. MISR stereographic views of thick aerosol plumes allow height and instantaneous wind derivations at spatial resolutions of better than 1.1 km horizontally and ~200m vertically. We will discuss the radiometric and stereo-retrieval capabilities of MISR specifically for dust, and demonstrate the use of MISR data in conjunction with other available satellite observations for dust property characterization and climate studies.First, we will discuss MISR non-spherical (dust) fraction product over the global oceans. We will show that over the Atlantic Ocean, changes in the MISR-derived non-spherical AOD fraction illustrate the evolution of dust during transport. Next, we will present a MISR satellite perspective on dust climatology in major dust source regions with a particular emphasis on the West Africa and Middle East and discuss MISR's unique strengths as well as current product biases. Finally, we will discuss MISR dust plume product and climatological applications.

dust sources↗

"Coulombic Viscosity" In Granular Materials: Planetary and Astrophysical Implications

The term "Coulombic viscosity" is introduced here to define an empirically observed phenomenon from experiments conducted in both microgravity, and in ground-based 1-g conditions. In the latter case, a sand attrition device was employed to test the longevity of aeolian materials by creating two intersecting grain-circulation paths or cells that would lead to most of the grain energy being expended on grain-to-grain collisions (simulating dune systems). In the areas in the device where gravitationally-driven grain-slurries recycled the sand, the slurries moved with a boundary-layer impeded motion down the chamber walls. Excessive electrostatic charging of the grains during these experiments was prevented by the use of an a.c. corona (created by a Tesla coil) through which the grains passed on every cycle. This created both positive and negative ions which neutralized the triboelectrically-generated grain charges. When the corona was switched on, the velocity of the wall-attached slurries increased by a factor of two as approximately determined by direct observation. What appeared to be a freely-flowing slurry of grains impeded only by intergranular mechanical friction, had obviously been significantly retarded in its motion by electrostatic forces between the grains; with the charging reduced, the grains were able to move past one another without a flow "viscosity" imposed by the Coulombic intergranular forces. A similar phenomenon was observed during microgravity experiments aboard Space Shuttle in USML-1 & USML-2 spacelabs where freely-suspended clouds of sand were being investigated for their potential to for-m aggregates. In this environment, the grains were also charged electrostatically (by natural processes prior to flight), but were free from the intervention of gravity in their interactions. The grains were dispersed into dense clouds by bursts of air turbulence and allowed to form aggregates as the ballistic and turbulent motions damped out. During this very brief (30-60 sec) damping period, motion of the grains was observed to be retarded by the electrostatic interactions. The fact that the grains almost instantly formed aggregates was evidence that their ballistic motions had been constrained and redirected by the dipole-dipole interactions that led to filamentary aggregate development. Undoubtedly, the "Coulombic viscosity" of the cloud assisted in damping grain motion so rapidly. The electrostatically-induced grain-cloud viscosity or drag exerted on grain motion, is a complex function of three major parameters: charge magnitude, charge sign, and mean intergranular distance. The above experiments illustrate one particular type of granular behavior. The discussion here will therefore be restricted to drag relationships: (a) between grains that are naturally charged triboelectrically and thus exhibit dipole-dipole attractions between one another even if there are slight net charges present (which can be overwhelmed by dipole coupling at short distances), and (b) between grains that are densely spaced where the intergranular distance varies between zero and some value (usually tens or hundreds of grain diameters) that permits each grain to detect the dipole moment of another grain -- the distance is not so great that other grains appears as neutral electrical "singularities. I. Aeolian transport: During motion of grains in a saltation cloud (on Earth, Mars, or Venus), triboelectric charging must occur as a result of multiple grain contacts, and by friction with the entraining air. A situation might develop that is similar to the one described above in the attrition device: grain motion becoming significantly retarded (reduced flux) as grains find it increasingly difficult to either separate from the surface, or to pass one another without Coulombic retarding forces. A "Coulombic drag" will exist at flux initiation and increase with time to work in direct opposition to the aerodynamic drag that drives the grain motion. It is predicted that this will lead to an increase with time of both the aerodynamic and bed-dilatancy thresholds (3). Because of Paschen discharge effects in the martian atmosphere, the electrostatic charging in a saltation cloud may be partially abated, but this will lead to greater grain mobility, more charging, and thus to a charge-discharge steady state mediated by mechanical interactions. II. Dry colluvial systems: Sand avalanches on dunes, dry debris flows, talus flows, avalanches, and pyroclastic surges are examples of gravity-driven, dense granular flows where rock/grain fragmentation and grain-to-grain interactions cause triboelectrification (sometimes augmented by other electrical charging processes), and where the grain densities of the systems are such that strong dipole-dipole interactions between grains might be expected to be present. Because it is expected that the Coulombic forces between grains will cause a sluggishness or enhanced granular-flow viscosity, the motion of a grain mass will be retarded or damped so that this will assist, ultimately, in terminating the flow. The greatest Coulombic viscosity will be created in the most highly charged systems, which will also be the most energetic. Thus, grain flows have some tendency to be self-limiting by internal energy partitioning; gravitational potential is converted to Coulombic potential, which manifests itself as a drag force between the grains. III. Volcanic eruption plumes and impact ejecta curtains: The violence of these systems leads to powerful electrical charging of particulates. Lightning storms emanating from volcanic plumes are a testimony to the levels of charging. As pyroclastic grains interact forcefully and frequently within eruption plumes, it is reasonable to predict that the internal turbulent motions of the plume will be significantly damped by the Coulombic viscosity exerted by grain charges. Additional information is contained in the original.

Marshall, J. R.↗

The Long-term Transport and Radiative Impacts of the 2017 British Columbia Pyrocumulonimbus Smoke Aerosols in the Stratosphere

Interactions of meteorology with wildfires in British Columbia, Canada during August 2017 led to three major pyrocumulonimbus (pyroCb) events that resulted in the injection of large amounts of smoke aerosols and other combustion products at the local upper troposphere and lower stratosphere (UTLS). These plumes of UTLS smoke with elevated values of aerosol extinction and backscatter compared to the background state were readily tracked by multiple satellite-based instruments as they spread across the Northern Hemisphere (NH). The plumes were observed in the lower stratosphere for about 8-10 months following the fire injections, with a stratospheric aerosol e-folding time of about 5 months. To investigate the radiative impacts of these events on the Earth system, we performed a number of simulations with the Goddard Earth Observing System (GEOS) atmospheric general circulation model (AGCM). Observations from multiple remote-sensing instruments were used to calibrate the injection parameters (location, amount, composition and heights) and optical properties of the smoke aerosols in the model. The resulting simulations of three-dimensional smoke transport were evaluated for a year from the day of injections using daily observations from OMPS-LP (Ozone Mapping Profiler Suite Limb Profiler). The model simulated rate of ascent, hemispheric spread and residence time (or e-folding time) of the smoke aerosols in the stratosphere are in close agreement with OMPS-LP observations. We found that both aerosol self-lofting and the large-scale atmospheric motion play important roles in lifting the smoke plumes from near the tropopause altitudes (~12 km) to about 22-23 km into the atmosphere. Further, our estimations of the radiative impacts of the pyroCb-emitted smoke aerosols showed that the smoke caused an additional warming of the atmosphere by about 0.6-1 W/m2 (zonal mean) that persisted for about 2-3 months after the injections in regions north of 40oN. The surface experienced a comparable magnitude of cooling. The atmospheric warming is mainly located in the stratosphere, coincident with the location of the smoke plumes, leading to an increase in zonal mean shortwave (SW) heating rates of 0.02-0.04 K/day during September 2017.

pyrocumulonimbus (pyroCb) events↗

Overview of ACE-Asia Spring 2001 Investigations on Aerosol Radiative Effects and Related Aerosol Properties

A primary, ACE-Asia objective was to quantify the interactions between aerosols and radiation in the Asia-Pacific region. Toward this end, radiometric and related aerosol measurements were made from ocean, land, air and space platforms. Models that predict aerosol fields guided the measurements and are helping integrate and interpret results. Companion overview's survey these measurement and modeling components. Here we illustrate how these components were combined to determine aerosol radiative. impacts and their relation to aerosol properties. Because clouds can obscure or change aerosol direct radiative effects, aircraft and ship sorties to measure these effects depended on predicting and finding cloud-free areas and times with interesting aerosols present. Pre-experiment satellite cloud climatologies, pre-flight aerosol and cloud forecasts, and in-flight guidance from satellite imagery all helped achieve this. Assessments of aerosol regional radiative impacts benefit from the spatiotemporal coverage of satellites, provided satellite-retrieved aerosol properties are accurate. Therefore, ACE-Asia included satellite retrieval tests, as part of many comparisons to judge the consistency (closure) among, diverse measurements. Early results include: (1) Solar spectrally resolved and broadband irradiances and optical depth measurements from the C-130 aircraft and at Kosan, Korea yielded aerosol radiative forcing efficiencies, permitting comparisons between efficiencies of ACE-Asia and INDOEX aerosols, and between dust and "pollution" aerosols. Detailed results will be presented in separate papers. (2) Based on measurements of wavelength dependent aerosol optical depth (AOD) and single scattering albedo the estimated 24-h a average aerosol radiative forcing efficiency at the surface for photosynthetically active radiation (400 - 700 nm) in Yulin, China is approx. 30 W sq m per AOD(500 nm). (3) The R/V Brown cruise from Honolulu to Sea of Japan sampled an aerosol optical depth gradient, with AOD(500 nm) extremes from 0.1 to 1.1. On the Pacific transit from Honolulu to Hachijo AOD(500 nm) averaged 0.2, including increases to 0.4 after several storms, suggesting the strong impact of wind-generated seasalt. The AOD maximum, found in the Sea of Japan, was influenced by dust and anthropogenic sources. (4) In Beijing, single scattering albedo retrieved from AERONET sun-sky radiometry yielded midvisible SSA=0.88 with strong wavelength dependence, suggesting a significant black carbon component. SSA retrieved during dust episodes was approx. 0.90 and variable but wavelength neutral reflecting the presence of urban haze with the dust. Downwind at Anmyon Island SSA was considerably higher, approx. 0.94, but wavelength neutral for dust episodes and spectrally dependent during non dust periods. (5) Satellite retrievals show major aerosol features moving from Asia over the Pacific; however, determining seasonal-average aerosol effects is hampered by sampling frequency and large-scale cloud systems that obscure key parts of aerosol patterns. Preliminary calculations using, satellite-retrieved AOD fields and initial ACE-Asia aerosol properties (including sulfates, soot, and dust) yield clear-sky aerosol radiative effects in the seasonal-average ACE-Asia plume exceeding those of manmade greenhouse gases. Quantifying all-sky direct aerosol radiative effects is complicated by the need to define the height of absorbing aerosols with respect to cloud decks.

Russell, Philip B.↗

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗