Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud boundaries”

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,081 records · Page 60

Individual Aerosol Particles from Biomass Burning in Southern Africa Compositions and Aging of Inorganic Particles: Compositions and Aging of Inorganic Particles - 2

Individual aerosol particles collected over southern Africa during the SAFARI 2000 field study were studied using transmission electron microscopy and field-emission scanning electron microscopy. The sizes, shapes, compositions, mixing states, surface coatings, and relative abundances of aerosol particles from biomass burning, in boundary layer hazes, and in the free troposphere were compared, with emphasis on aging and reactions of inorganic smoke particles. Potassium salts and organic particles were the predominant species in the smoke, and most were internally mixed. More KCl particles occur in young smoke, whereas more K2SO4 and KNO3 particles were present in aged smoke. This change indicates that with the aging of the smoke, KCl particles from the fires were converted to K2SO4 and KNO3 through reactions with sulfur- and nitrogen- bearing species from biomass burning as well as other sources. More soot was present in smoke from flaming grass fires than bush and wood fires, probably due to the predominance of flaming combustion in grass fires. The high abundance of organic particles and soluble salts can affect the hygroscopic properties of biomass-burning aerosols and therefore influence their role as cloud condensation nuclei. Particles from biomass burning were important constituents of the regional hazes.

Li, Jia↗

Performance of the GEOS-3/Terra Data Assimilation System in the Northern Stratospheric Winter 1999/2000

As part of NASA's support for the Terra satellite, which became operational in January 2000, the Data Assimilation Office introduced a new version of the GEOS data assimilation system (DAS) in November 1999. This system, GEOS-3/Terra, differs from its predecessor in several ways, notably through an increase in horizontal resolution (from 2-by-2.5 degrees to 1-by-1 degree), a slightly lower upper boundary (0.1 instead of 0.01hPa) with fewer levels (48 as opposed to 70), and substantial changes to the tropospheric physics package. This paper will address the performance of the GEOS-3/Terra DAS in the stratosphere. it focusses on the analyses (produced four times daily) and the five-day forecasts (produced twice daily). These were important for the meteorological support of the SAGE-3 Ozone Loss and Validation Experiment, based in Kiruna, Northern Sweden, in the winter of 1999/2000. It is shown that the analyses of basic meteorological fields (temperature, geopotential height, and horizontal wind) are in good agreement with those from other centers. The analyses captured the cold polar vortex which persisted through most of the winter. It is shown that forecasts (up to five days) tend to have a warm bias, which is important for the prediction of polar stratospheric clouds, which are triggered by temperatures of 195K (or lower). The importance of accurate upper tropospheric forecasts in predicting the stratospheric flow is highlighted in the context of the evolution of the shape of the stratospheric polar vortex. A prominent blocking high in the Atlantic region in January was an important factor determining the shape of the distorted lower stratospheric vortex; the predictive skill of these features was strongly coupled in the GEOS-3/Terra system.

Pawson, S.↗

Validation of AIRS V6 Surface Temperature over Greenland with GCN and NOAA Stations

This work compares the temporal and spatial characteristics of the AIRSAMSU (Atmospheric Infrared Sounder Advanced Microwave Sounding Unit A) Version 6 and MODIS (Moderate resolution Imaging Spectroradiometer) Collection 5 derived surface temperatures over Greenland. To estimate uncertainties in space-based surface temperature measurements, we re-projected the MODIS Ice Surface Temperature (IST) to 0.5 by 0.5 degree spatial resolution. We also re-gridded AIRS Skin Temperature (Ts) into the same grid but classified with different cloud conditions and surface types. These co-located data sets make intercomparison between the two instruments relatively straightforward. Using this approach, the spatial comparison between the monthly mean AIRS Ts and MODIS IST is in good agreement with RMS 2K for May 2012. This approach also allows the detection of any long-term calibration drift and the careful examination of calibration consistency in the MODIS and AIRS temperature data record. The temporal correlations between temperature data are also compared with those from in-situ measurements from GC-Net (GCN) and NOAA stations. The coherent time series of surface temperature evident in the correlation between AIRS Ts and GCN temperatures suggest that at monthly time scales both observations capture the same climate signal over Greenland. It is also suggested that AIRS surface air temperature (Ta) can be used to estimate the boundary layer inversion.

AIRS↗

Integrating Cloud-Based Workflows in Continental-Scale Cropland Extent Classification

Accurate information on cropland spatial distribution is required for global-scale assessments and agricultural land use policies. Cloud computing platforms such as Google Earth Engine (GEE) provide unprecedented opportunities for large-scale classifications of Landsat data. We developed a novel method to fuse pixel-based random forest classification of continental-scale Landsat data on GEE and an object-based segmentation approach known as recursive hierarchical segmentation (RHSeg). Using our fusion method, we produced a continental-scale cropland extent map for North America at 30m spatial resolution for the nominal year 2010. The total cropland area for North America was estimated at 275.18 million hectares (Mha). The overall accuracies of the map are>90% across the continent. This map also compares well with the United States Department of Agriculture (USDA) cropland data layer (CDL), Agriculture and Agri-food Canada (AAFC) annual crop inventory (ACI), and the Mexican government agency Servicio de Informacion Agroalimentaria y Pesquera (SIAP)'s agricultural boundaries. Furthermore, our map compared well with sub-country statistics including state-wise and county-wise cropland statistics in regression models resulting in R2 > 0.84. This key contribution paves the way for more detailed products such as crop intensity, crop type, and crop irrigation, and provides a method for creating high-resolution cropland extent maps for other countries where spatial information about croplands are not as prevalent.

Massey, Richard↗

Small scale structure and mixing at the edge of the Antarctic vortex

Small scale correlations and patterns in the chemical tracers measured from the NASA ER-2 aircraft in the 1987 AAOE campaign can be used to investigate the structure of the edge of the polar vortex and the chemically perturbed region within it. Examples of several types of transport processes can be found in the data. Since ClO and O3 have similar vertical gradients and opposite horizontal gradients near the chemically perturbed region, the correlation between ClO and O3 can be used to study the extent of horizontal transport at the edge of the chemically perturbed region. Horizontal transport dominates the correlation for a latitude band up to 4 degrees on each side of the boundary. This implies a transition zone containing a substantial fraction of the mass of the total polar vortex. Similar horizontal transport can be seen in other tracers as well. It has not been possible to distinguish reversible transport from irreversible mixing. One manifestation of the horizontal transport is that the edge of the chemically perturbed region is often layered rather than a vertical curtain. This can be seen from the frequent reversed vertical gradients of NO2, caused by air with high NO2 overlapping layers with lower mixing ratios. Water and NO2 are positively correlated within the chemically perturbed region. This is the opposite sign to the correlation in the unperturbed stratosphere. The extent of the positive correlation is too great to be attributed solely to horizontal mixing. Instead, it is hypothesized that dehydration and descent are closely connected on a small scale, possibly due to radiative cooling of the clouds that also cause ice to fall to lower altitudes.

Murphy, D. M.↗

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗

Origin, Transport, and Vertical Distribution of Atmospheric Polluntants over the Northern Sourth China Sea During the 7-SEAS-Dongsha Experiment

During the spring of 2010, comprehensive in situ measurements were made for the first time on a small atoll (Dongsha Island) in the northern South China Sea (SCS), a key region of the 7-SEAS (the Seven South East Asian Studies) program. This paper focuses on characterizing the source origins, transport processes, and vertical distributions of the Asian continental outflows over the region, using measurements including mass concentration, optical properties, hygroscopicity, and vertical distribution of the aerosol particles, as well as the trace gas composition. Cluster analysis of backward trajectories classified 52% of the air masses arriving at ground level of Dongsha Island as having a continental origin, mainly from northern China to the northern SCS, passing the coastal area and being confined in the marine boundary layer (0-0.5 km). Compared to aerosols of oceanic origin, the fine mode continental aerosols have a higher concentration, extinction coefficient, and single-scattering albedo at 550 nm (i.e., 19 vs. 14 microg per cubic meter in PM(sub 2.5); 77 vs. 59 M per meter in beta(sub e); and 0.94 vs. 0.90 in omega, respectively). These aerosols have a higher hygroscopicity (f at 85% RH = 2.1) than those in the upwind inland regions, suggesting that the aerosols transported to the northern SCS were modified by the marine environment. In addition to the near-surface aerosol transport, a significant upper-layer (3-4 km) transport of biomass-burning aerosols was observed. Our results suggest that emissions from both China and Southeast Asia could have a significant impact on the aerosol loading and other aerosol properties over the SCS. Furthermore, the complex vertical distribution of aerosols-coinciding-with-clouds has implications for remote-sensing observations and aerosol-cloud-radiation interactions.

pollutants↗

"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 Dryline on 22 May 2002 during IHOP_2002: Convective-Scale Measurements at the Profiling Site

A detailed analysis of the structure of a double dryline observed over the Oklahoma panhandle during the first International H2O Project (IHOP_2002) convective initiation (CI) mission on 22 May 2002 is presented. A unique and unprecedented set of high temporal and spatial resolution measurements of water vapor mixing ratio, wind, and boundary layer structure parameters were acquired using the National Aeronautics and Space Administration (NASA) scanning Raman lidar (SRL), the Goddard Lidar Observatory for Winds (GLOW), and the Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE), respectively. These measurements are combined with the vertical velocity measurements derived from the National Center for Atmospheric Research (NCAR) Multiple Antenna Profiler Radar (MAPR) and radar structure function from the high-resolution University of Massachusetts frequency-modulated continuous-wave (FMCW) radar to reveal the evolution and structure of the late afternoon double-dryline boundary layer. The eastern dryline advanced and then retreated over the Homestead profiling site in the Oklahoma panhandle, providing conditions ripe for a detailed observation of the small-scale variability within the boundary layer and the dryline. In situ aircraft data, dropsonde and radiosonde data, along with NCAR S-band dual-polarization Doppler radar (S-Pol) measurements, are also used to provide the larger-scale picture of the double-dryline environment. Moisture and temperature jumps of about 3 g kg(sup -1) and 1 -2 K, respectively, were observed across the eastern radar fine line (dryline), more than the moisture jumps (1-2 g kg(sup -1)) observed across the western radar fine line (secondary dryline). Most updraft plumes observed were located on the moist side of the eastern dryline with vertical velocities exceeding 3 m s(sup -1) and variable horizontal widths of 2-5 km, although some were as wide as 7-8 km. These updrafts were up to 1.5 g kg(sup -1) moister than the surrounding environment. Although models suggested deep convection over the Oklahoma panhandle and several cloud lines were observed near the dryline, the dryline itself did not initiate any storms over the intensive observation region (IOR). Possible reasons for this lack of convection are discussed. Strong capping inversion and moisture detrainment between the lifting condensation level and the level of free convection related to an overriding drier air, together with the relatively small near-surface moisture values (less than 10 g kg(sup -1)), were detrimental to CI in this case.

Demoz, Belay↗

Clutter Assessment for an Autonomous Multi-Agent Search Mission

This paper presents a method for evaluating the amount of clutter in a region where autonomous vehicles in a multi-agent system must operate based on LIDAR point cloud measurements. The point cloud is used to generate an occupancy grid which is then projected onto a 2D plane of vehicle motion, constituting an image. A series of Gaussian radial basis functions (GRBFs) is created, each centered at an occupied pixel in a 2D image, and summed together to form the clutter field. The clutter field is a representation of the density and permeability of the space at each coordinate. The clutter field is then approximated such that iso-clutter contours are simple geometric objects so that intelligent machine assets can easily query the distance between them and any given point in an environment. In this way, agents are able to determine whether to enter into or steer away from areas of interest. Each vehicle has a clutter threshold representing the clutter value of the space in which it can safely maneuver. The iso-clutter contour corresponding to a vehicle’s clutter threshold is treated as the boundary of an obstacle to be avoided. A simulation is presented where a multi-agent system is tasked with persistent observation of a cluttered area. Each vehicle in the simulation has a different clutter threshold. The vehicles use a potential field-based guidance algorithm, and an allocation of vehicles to specific regions of the space emerges.

multi-agent↗

X-Band Radar for Studies of Tropical Storms from High Altitude UAV Platform

The increased role of unmanned aerial vehicles (UAV) in NASA's suborbital program has created a strong interest in the development of instruments with new capabilities, more compact sizes and reduced weights than the instruments currently operated on manned aircrafts. There is a strong demand and tremendous potential for using high altitude UAV (HUAV) to carry weather radars for measurements of reflectivity and wind fields from tropical storms. Tropical storm genesis frequently occurs in ocean regions that are inaccessible to piloted aircraft due to the long off shore range and the required periods of time to gather significant data. Important factors of interest for the study of hurricane genesis include surface winds, profiled winds, sea surface temperatures, precipitation, and boundary layer conditions. Current satellite precipitation and surface wind sensors have resolutions that are too large and revisit times that are too infrequent to study this problem. Furthermore, none of the spaceborne sensors measure winds within the storm itself. A dual beam X-band Doppler radar, UAV Radar (URAD), is under development at the NASA Goddard Space Flight Center for the study of tropical storms from HUAV platforms, such as a Global Hawk. X-band is the most desirable frequency for airborne weather radars since these can be built in a relatively compact size using off-the-shelf components which cost significantly less than other higher frequency radars. Furthermore, X-band radars provide good sensitivity with tolerable attenuation in storms. The low-cost and light-weight URAD will provide new capabilities for studying hurricane genesis by analyzing the vertical structure of tropical cyclones as well as 3D reflectivity and wind fields in clouds. It will enable us to measure both the 3D precipitation structure and surface winds by using two antenna beams: fixed nadir and conical scanning each produced by its associated subsystem. The nadir subsystem is a magnetron based radar modified from a marine radar transceiver. It is capable of measuring vertical reflectivity and velocity profile while being a lower-cost, smaller size, and lighter weight version of the NASA ER-2 Doppler Radar (EDOP), which has flown during many NASA field campaigns and has provided valuable scientific information on hurricanes and weather phenomena. Unfortunately, EDOP is too large and heavy for most UAV platforms, but the experience gained with this instrument provided us with the heritage to build a new low-cost, light-weight, smaller system that will be capable of flying on UAVs. The scanning subsystem uses a TWT transmitter and provides measurements of 3D reflectivity/wind fields in-clouds. Conical scanning of the radar beam at a 35 deg. incidence angle will also provide information of surface wind speed and direction derived from the surface return over a single 360 deg. sweep. URAD data system will be Linux based with the capability of autonomous operation. It will utilize cutting edge digital receiver and FPGA technologies to carry out the data acquisition and processing tasks. High speed navigation data from the aircraft will also be captured and saved along with radar data for 3D measurement field reconstruction and aircraft motion correction. There is a tremendous potential for UAVs to carry down-looking weather radars for measurements of reflectivity, horizontal and vertical winds from tropical storms. With operation from HUAV platforms, the dual beam X-band radar under development promises to provide greatly needed information for tropical storm research.

Rodriquez, Shannon↗

Assessment of Aerosol Distributions from GEOS-5 Using the CALIPSO Feature Mask

A-train sensors such as MODIS, MISR, and CALIPSO are used to determine aerosol properties, and in the process a means of estimating aerosol type (e.g. smoke vs. dust). Correct classification of aerosol type is important for climate assessment, air quality applications, and for comparisons and analysis with aerosol transport models. The Aerosols-Clouds-Ecosystems (ACE) satellite mission proposed in the NRC Decadal Survey describes a next generation aerosol and cloud suite similar to the current A-train, including a lidar. The future ACE lidar must be able to determine aerosol type effectively in conjunction with modeling activities to achieve ACE objectives. Here we examine the current capabilities of CALIPSO and the NASA Goddard Earth Observing System general circulation model and data assimilation system (GEOS-5), to place future ACE needs in context. The CALIPSO level 2 feature mask includes vertical profiles of aerosol layers classified by type. GEOS-5 provides global 3D aerosol mass for sulfate, sea salt, dust, and black and organic carbon. A GEOS aerosol scene classification algorithm has been developed to provide estimates of aerosol mixtures and extinction profiles along the CALIPSO orbit track. In previous work, initial comparisons between GEOS-5 derived aerosol mixtures and CALIPSO derived aerosol types were presented for July 2007. In general, the results showed that model and lidar derived aerosol types did not agree well in the boundary layer. Agreement was poor over Europe, where CALIPSO indicated the presence of dust and pollution mixtures yet GEOS-5 was dominated by pollution with little dust. Over the ocean in the tropics, the model appeared to contain less sea salt than detected by CALIPSO, yet at high latitudes the situation was reserved. Agreement between CALIPSO and GEOS-5, aerosol types improved above the boundary layer, primarily in dust and smoke dominated regions. At higher altitudes (> 5 km), the model contained aerosol layers not detected by CALIPSO. Here we present new results for a full year study using the new Version 3 CALIPSO data and most recent GEOS-5 model results.

Welton, Ellsworth↗

Evaluating Precipitation Features and Rainfall Characteristics in a Multi-scale Modeling Framework

Cloud and precipitation systems over the tropics and subtropics are simulated with a multi-scale modeling framework (MMF) and compared against the TRMM radar precipitation features (RPFs) product. A methodology, in close analogy to the TRMM RPFs, is developed to analyze simulated cloud precipitating structures from the embedded two-dimensional cloud-resolving models (CRMs) within an MMF. Despite the two-dimensionality of the CRMs, the simulated RPFs population distribution, and horizontal and vertical structure are in good agreement with TRMM observations. However, some deficits are also found in the model simulations. The model tends to overestimate mean convective precipitation rates for RPFs with a size less than 100 km, contributing to the excessive precipitation biases in the warm pool and western Pacific, western and northern India Ocean, and eastern Pacific commonly found in most MMFs. For large features with a size greater than 150 km, both convective and stratiform rain rates are underestimated. The distribution of maximum radar echo top heights as a function of RPF size is well simulated except the model tends to underestimate the occurrence frequency of maximum heights greater than 15 km. The maximum echo top heights for convective cells embedded within large RPFs with a size greater than 150 km are also underestimated. The cyclic lateral boundary with a limited model domain generates artificial occurrences for RPFs with a size close to the model domain size, producing a significant contribution to the total rainfall due to their sizes. This cyclic lateral boundary effect can be easily identified and quantified in both probability and cumulative distribution functions of RPFs. The geophysical distribution of the population of the largest RPFs in the control experiment shows they are mainly located in the Subtropics but also partially contribute to the common MMF biases of excessive precipitation in the Tropics. Sensitivity experiments using CRMs with different domain sizes and different grid spacings show larger domains (higher resolution) tend to shift the RPFs distribution to large (small) sizes. The cyclic lateral boundary biases increase as CRM domain size decreases. The impacts of model horizontal and vertical resolution on simulated convective systems are also investigated.

Chern, Jiun-Dar↗

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↗

The Earth Surface Mineral Dust Source Investigation: An Earth Science Imaging Spectroscopy Mission

The Earth Surface Mineral Dust Source Investigation, EMIT, is planned to operate from the International Space Station starting no earlier than the fall of 2021. EMIT will use visible to short wavelength infrared imaging spectroscopy to determine the mineral composition of the arid land dust source regions of the Earth to advance our knowledge of the radiative forcing effect of these aerosols. Mineral dust emitted into the atmosphere under high wind conditions is an element of the Earth system with many impacts to the Earth’s energy balance, atmosphere, surface, and oceans. The Earth’s mineral dust cycle with source, transport, and deposition phases are studied with advanced Earth System Models. Because the chemical composition, optical and surface properties of soil particles vary strongly with the mineral composition of the source, these models require knowledge of surface soil mineral dust source composition knowledge to accurately understand dust impacts on the Earth system now and in the future. At present, compositional knowledge of the Earth’s mineral dust source regions from existing data sets is uncertain as a result of limited measurements. EMIT will use spectroscopically derived surface mineral composition to update the prescribe boundary conditions for state-of-the-art Earth System Models. The EMIT-initialized models will be used to investigate the impact of direct radiative forcing in the Earth system that depends strongly on the composition of the mineral dust aerosols emitted into the atmosphere. These new measurements and related products will be used to address the EMIT science objectives and made available to the science community for additional investigations. An overview of the EMIT science, development, and mission is presented in this paper.

Green, Robert O.↗

Detection of Water Hazards for Autonomous Robotic Vehicles

Four methods of detection of bodies of water are under development as means to enable autonomous robotic ground vehicles to avoid water hazards when traversing off-road terrain. The methods involve processing of digitized outputs of optoelectronic sensors aboard the vehicles. It is planned to implement these methods in hardware and software that would operate in conjunction with the hardware and software for navigation and for avoidance of solid terrain obstacles and hazards. The first method, intended for use during the day, is based on the observation that, under most off-road conditions, reflections of sky from water are easily discriminated from the adjacent terrain by their color and brightness, regardless of the weather and of the state of surface waves on the water. Accordingly, this method involves collection of color imagery by a video camera and processing of the image data by an algorithm that classifies each pixel as soil, water, or vegetation according to its color and brightness values (see figure). Among the issues that arise is the fact that in the presence of reflections of objects on the opposite shore, it is difficult to distinguish water by color and brightness alone. Another issue is that once a body of water has been identified by means of color and brightness, its boundary must be mapped for use in navigation. Techniques for addressing these issues are under investigation. The second method, which is not limited by time of day, is based on the observation that ladar returns from bodies of water are usually too weak to be detected. In this method, ladar scans of the terrain are analyzed for returns and the absence thereof. In appropriate regions, the presence of water can be inferred from the absence of returns. Under some conditions in which reflections from the bottom are detectable, ladar returns could, in principle, be used to determine depth. The third method involves the recognition of bodies of water as dark areas in short-wavelength infrared (SWIR) images. This method is based on the fact, well known among experts in remote sensing, that water bodies of any appreciable depth appear very dark in near-infrared, overhead imagery. Even under a thick layer of marine fog, SWIR illumination is present. Hence, this method may work even in the presence of clouds, though it is unlikely to work at night. Snow and ice also exhibit very strong absorption at wavelengths greater than about 1.4 m. Hence, the wavelength range of about 1.5 to 1.6 m might be useable in this method for recognizing water, snow, and ice. One notable drawback of this method is that useful look-ahead distance could be limited by surface reflections. The fourth method, intended for use at night, involves the contrast between water and terrain in thermal-infrared (medium-wavelength infrared) imagery. This method is based on the fact that at night, water is usually warmer than the adjacent terrain. Look-ahead distance could be limited in this method because, for reasons not yet fully understood, water appears to darken in the thermal infrared with increasing distance.

Matthes, Larry↗

Validation and Development of the GPCP Experimental One-Degree Daily (1DD) Global Precipitation Product

The One-Degree Daily (1DD) precipitation dataset has been developed for the Global Precipitation Climatology Project (GPCP) and is currently in beta test preparatory to release as an official GPCP product. The 1DD provides a globally-complete, observation-only estimate of precipitation on a daily 1 deg. x 1 deg. grid for the period 1997 through early 2000 (by the time of the conference). In the latitude band 40N-40S the 1DD uses the Threshold-Matched Precipitation Index (TMPI), a GPI-like IR product with the pixel-level T(sub b) threshold and (single) conditional rain rate determined locally for each month by the frequency of precipitation in the GPROF SSM/I product and by, the precipitation amount in the GPCP monthly satellite-gauge (SG) combination. Outside 40N-40S the 1DD uses a scaled TOVS precipitation estimate that has month-by-month adjustments based on the TMPI and the SG. Early validation results are encouraging. The 1DD shows relatively large scatter about the daily validation values in individual grid boxes, as expected for a technique that depends on cloud-sensing schemes such as the TMPI and TOVS. On the other hand, the time series of 1DD shows good correlation with validation in individual boxes. For example, the 1997-1998 time series of 1DD and Oklahoma Mesonet values in a grid box in northeastern Oklahoma have the correlation coefficient = 0.73. Looking more carefully at these two time series, the number of raining days for the 1DD is within 7% of the Mesonet value, while the distribution of daily rain values is very similar. Other tests indicate that area- or time-averaging improve the error characteristics, making the data set highly attractive to users interested in stream flow, short-term regional climatology, and model comparisons. The second generation of the 1DD product is currently under development; it is designed to directly incorporate TRMM and other high-quality precipitation estimates. These data are generally sparse because they are observed by low-orbit satellites, so a fair amount of work must be devoted to analyzing the effect of data boundaries. This work is laying, the groundwork for effective use of the NASA Global Precipitation Mission, which will have full Global coverage by low-orbit passive microwave satellites every three hours.

Huffman, George J.↗

Evolution of dynamo-generated magnetic fields in accretion disks around compact and young stars

Geometrically thin, optically thick, turbulent accretion disks are believed to surround many stars. Some of them are the compact components of close binaries, while the others are throught to be T Tauri stars. These accretion disks must be magnetized objects because the accreted matter, whether it comes from the companion star (binaries) or from a collapsing molecular cloud core (single young stars), carries an embedded magnetic field. In addition, most accretion disks are hot and turbulent, thus meeting the condition for the MHD turbulent dynamo to maintain and amplify any seed field magnetic field. In fact, for a disk's magnetic field to persist long enough in comparison with the disk viscous time it must be contemporaneously regenerated because the characteristic diffusion time of a magnetic field is typically much shorter than a disk's viscous time. This is true for most thin accretion disks. Consequently, studying magentic fields in thin disks is usually synonymous with studying magnetic dynamos, a fact that is not commonly recognized in the literature. Progress in studying the structure of many accretion disks was achieved mainly because most disks can be regarded as two-dimensional flows in which vertical and radial structures are largely decoupled. By analogy, in a thin disk, one may expect that vertical and radial structures of the magnetic field are decoupled because the magnetic field diffuses more rapidly to the vertical boundary of the disk than along the radius. Thus, an asymptotic method, called an adiabatic approximation, can be applied to accretion disk dynamo. We can represent the solution to the dynamo equation in the form B = Q(r)b(r,z), where Q(r) describes the field distribution along the radius, while the field distribution across the disk is included in the vector function b, which parametrically depends on r and is normalized by the condition max (b(z)) = 1. The field distribution across the disk is established rapidly, while the radial distribution Q(r) evolves on a considerably longer timescale. It is this evolution that is the subject of this paper.

Stepinski, Tomasz F.↗