Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrometeor”

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 289 records · Page 16

Estimates of Random Error in Satellite Rainfall Averages

Satellite rain estimates are most accurate when obtained with microwave instruments on low earth-orbiting satellites. Estimation of daily or monthly total areal rainfall, typically of interest to hydrologists and climate researchers, is made difficult, however, by the relatively poor coverage generally available from such satellites. Intermittent coverage by the satellites leads to random "sampling error" in the satellite products. The inexact information about hydrometeors inferred from microwave data also leads to random "retrieval errors" in the rain estimates. In this talk we will review approaches to quantitative estimation of the sampling error in area/time averages of satellite rain retrievals using ground-based observations, and methods of estimating rms random error, both sampling and retrieval, in averages using satellite measurements themselves.

Bell, Thomas L.↗

CRYSTAL/FACE

Droplet Measurement Technologies (DMT), under funding from NASA, participated in the CRYSTAL/FACE field campaign in July, 2002 with measurements of cirrus cloud hydrometeors in the size range from 0.5 to 1600 microns. The measurements were made with the DMT Cloud, Aerosol and Precipitation Spectrometer (CAPS) that was flown on NASA's WB57F. With the exception of the first research flight when the data system failed two hours into the mission, the measurement system performed almost flawlessly during the thirteen flights. The measurements from the CAPS have been essential for interpretation of cirrus cloud properties and their impact on climate. The CAPS data set has been used extensively by the CRYSTAL/FACE investigators and as of the date of this report, have been included in five published research articles, 10 conference presentations and six other journal articles currently in preparation.

Baumgardner, Darrel↗

Microphysical Characterization of Ice in Hurricane Erin for Wideband Passive Microwave Comparisons

In order to better understand the characteristics of frozen hydrometeors in hurricane systems, computed brightness temperatures were compared with observations of Hurricane Erin from CAMEX-4 radiometers on the ER-2 aircraft. These wideband observations ranged from 10.7 to 183 plus or minus 10 GHz. In order to initialize the cloud profiles used in the radiative transfer calculations, radar reflectivities from an ER-2 instrument were converted to mass contents at specified height levels. Dropsonde data from the ER-2 and modeled Hurricane Andrew data were used to generate profiles of temperature, pressure, and relative humidity and to provide surface wind speed and ocean temperature. Four different ice (size and density) parameterizations were used to compute brightness temperature values. The comparisons between the calculated and observed brightness temperatures show that in general the Grenfell and Warren particle size distribution does the best job of matching the observations for most of the high frequencies. These high frequencies (greater than or equal to 85 GHz) are particularly sensitive to frozen particles. The Sehkon and Srivastava particle size distribution produced brightness temperatures that were too cold in most cases. On the other hand, the fluffy snow particle parameterizations with and without a melting layer produced brightness temperatures that were much warmer than the observations at these high frequencies. All low frequency calculations track the observations well except in the eye of Hurricane Erin where the surface wind speed and/or near surface temperatures were not properly adjusted to warm the brightness temperature calculations to values closer to the observations. Overall, this work shows that selecting a proper ice particle parameterization can be difficult and additional in situ measurements are needed to define appropriate parameterizations.

Skofronick-Jackson, G.↗

A Study of the Response of Deep Tropical Clouds to Mesoscale Processes: Sensitivities to Microphysics, Radiation, and Surface Fluxes - Part 2

The Goddard Cumulus Ensemble (GCE) model is used to examine the sensitivities of surface fluxes, explicit radiation, and ice microphysical processes on multi-day simulations of deep tropical convection over the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). The simulations incorporate large-scale advective temperature and moisture forcing, as well as large-scale momentum, that are updated every time step on a periodic lateral boundary grid. This study shows that when surface fluxes are eliminated, the mean atmosphere is much cooler and drier, convection and CAPE are much weaker, precipitation is less, and cloud coverage in stratiform regions much greater. Surface fluxes using the TOGA COARE flux algorithm are weaker than with the aerodynamic formulation, but closer to the observed fluxes. In addition, similar trends noted above for the case without surface fluxes are produced for the TOGA flux case, albeit to a much lesser extent. The elimination of explicit shortwave and longwave radiation is found to have only minimal effects on the mean thermodynamics, convection, and precipitation. However explicit radiation does have a significant impact on cloud temperatures and structure above 200 mb and on the overall mean vertical circulation. The removal of ice processes produces major changes in the structure of the cloud. Much of the liquid water is transported aloft and into anvils above the melting layer (600 mb), leaving narrow, but intense bands of rainfall in convective regions. The elimination of melting processes leads to greater hydrometeor mass below the melting layer, and produces a much warmer and moister boundary layer, leading to a greater mean CAPE. Finally, the elimination of the graupel species has only a small impact on mean total precipitation, thermodynamics, and dynamics of the simulation, but does produce much greater snow mass just above the melting layer. Some of these results differ from previous CRM studies of tropical systems, which is likely due to the type of simulated system, total time integration, and model setup.

Johnson, Daniel↗

Classification of Tropical Oceanic Precipitation using High Altitude Aircraft: Microwave and Electric Field Measurements

During the 1998 and 2001 hurricane seasons of the western Atlantic Ocean and Gulf of Mexico, the Advanced Microwave Precipitation Radiometer (AMPR), the ER-2 Doppler (EDOP) radar, and the Lightning Instrument Package (LIP) were flown aboard the National Aeronautics and Space Administration ER-2 high altitude aircraft as part of the Third Convection and Moisture Experiment (CAMEX-3) and the Fourth Convection and Moisture Experiment (CAMEX-4). Several hurricanes, tropical storms, and other precipitation systems were sampled during these experiments. An oceanic rainfall screening technique has been developed using AMPR passive microwave observations of these systems collected at frequencies of 10.7, 19.35,37.1, and 85.5 GHz. This technique combines the information content of the four AMPR frequencies regarding the gross vertical structure of hydrometeors into an intuitive and easily executable precipitation mapping format. The results have been verified using vertical profiles of EDOP reflectivity and lower altitude horizontal reflectivity scans collected by the National Oceanic and Atmospheric Administration WP-3D Orion radar. Matching the rainfall classification results with coincident electric field information collected by the LIP readily identifies convective rain regions within the precipitation fields. This technique shows promise as a real-time research and analysis tool for monitoring vertical updraft strength and convective intensity from airborne platforms such as remotely operated or uninhabited aerial vehicles. The technique is analyzed and discussed for a wide variety of precipitation types using the 26 August 1998 observations of Hurricane Bonnie near landfall.

Hood, Robbie E.↗

Multivariate Statistical Inference of Lightning Occurrence, and Using Lightning Observations

Two classes of multivariate statistical inference using TRMM Lightning Imaging Sensor, Precipitation Radar, and Microwave Imager observation are studied, using nonlinear classification neural networks as inferential tools. The very large and globally representative data sample provided by TRMM allows both training and validation (without overfitting) of neural networks with many degrees of freedom. In the first study, the flashing / or flashing condition of storm complexes is diagnosed using radar, passive microwave and/or environmental observations as neural network inputs. The diagnostic skill of these simple lightning/no-lightning classifiers can be quite high, over land (above 80% Probability of Detection; below 20% False Alarm Rate). In the second, passive microwave and lightning observations are used to diagnose radar reflectivity vertical structure. A priori diagnosis of hydrometeor vertical structure is highly important for improved rainfall retrieval from either orbital radars (e.g., the future Global Precipitation Mission "mothership") or radiometers (e.g., operational SSM/I and future Global Precipitation Mission passive microwave constellation platforms), we explore the incremental benefit to such diagnosis provided by lightning observations.

Boccippio, Dennis↗

A Model for Estimation of Rain Rate on Tropical Land from TRMM Microwave Imager Radiometer Observations

Over the tropical land regions observations of the 85 GHz brightness temperature (T(sub 85v)) made by the TRMM Microwave Imager (TMI) radiometer when analyzed with the help of rain rate (R(sub pR)) deduced from the TRMM Precipitation Radar (PR) indicate that there are two maxima in rain rate. One strong maximum occurs when T(sub 85) has a value of about 220 K and the other weaker one when T(sub 85v) is much colder approx. 150 K. Together with the help of earlier studies based on airborne Doppler Radar observations and radiative transfer theoretical simulations, we infer the maximum near 220 K is a result of relatively weak scattering due to super cooled rain drops and water coated ice hydrometeors associated with a developing thunderstorm (Cb) that has a strong updraft. The other maximum is associated with strong scattering due to ice particles that are formed when the updraft collapses and the rain from the Cb is transit2oning from convective type to stratiform type. Incorporating these ideas and with a view to improve the estimation of rain rate from existing operational method applicable to the tropical land areas, we have developed a rain retrieval model. This model utilizes two parameters, that have a horizontal scale of approx. 20km, deduced from the TMI measurements at 19, 21 and 37 GHz (T(sub 19v), T(sub 21v), T(sub 37v). The third parameter in the model, namely the horizontal gradient of brightness temperature within the 20 km scale, is deduced from TMI measurements at 85 GHz. Utilizing these parameters our retrieval model is formulated to yield instantaneous rain rate on a scale of 20 km and seasonal average on a mesoscale that agree well with that of the PR.

Prabhakara, C.↗

A Convective Vorticity Vector Associated With Tropical Convection: A 2D Cloud-Resolving Modeling Study

Although dry/moist potential vorticity is a useful physical quantity for meteorological analysis, it cannot be applied to the analysis of 2D simulations. A convective vorticity vector (CVV) is introduced in this study to analyze 2D cloud-resolving simulation data associated with 2D tropical convection. The cloud model is forced by the vertical velocity, zonal wind, horizontal advection, and sea surface temperature obtained from the TOGA COARE, and is integrated for a selected 10-day period. The CVV has zonal and vertical components in the 2D x-z frame. Analysis of zonally-averaged and mass-integrated quantities shows that the correlation coefficient between the vertical component of the CVV and the sum of the cloud hydrometeor mixing ratios is 0.81, whereas the correlation coefficient between the zonal component and the sum of the mixing ratios is only 0.18. This indicates that the vertical component of the CVV is closely associated with tropical convection. The tendency equation for the vertical component of the CVV is derived and the zonally-averaged and mass-integrated tendency budgets are analyzed. The tendency of the vertical component of the CVV is determined by the interaction between the vorticity and the zonal gradient of cloud heating. The results demonstrate that the vertical component of the CVV is a cloud-linked parameter and can be used to study tropical convection.

Gao, Shou-Ting↗

Graupel and Snow Determination in Convection Using Dual-Wavelength Radar

Convective storms contain significant amounts of snow and graupel that each have different scattering properties for radar and radiometer remote sensing. TRMM TMI algorithms have relied heavily on radiative transfer calculations based on cloud model predicted microphysics information. In general, without in situ aircraft measurements, it is very difficult to separate snow and graupel using remote measurements. At the 94 GHZ frequency, snow is weakly attenuating while graupel is strongly attenuating. With a non-attenuating frequency such as X-band, one can distinguish Mie scattering from larger hydrometeors as well as attenuation from graupel. The presence of supercooled water will also cause attenuation of the 94 GHz radar signal, but this is mainly in convective updrafts. We will present preliminary results aimed at distinguishing snow from graupel using dual-wavelength (9.6 GHz, 94 GHz) measurements from aircraft overflights of convective storms in Florida. Measurements were made with the 9.6 GHz ER-2 Doppler Radar (EDOP) and the 94 GHz Cloud Radar System (CRS) over a number of different precipitation systems. Analysis of the vertical profiles at the two wavelengths provides information on snow and graupel regions, as well as supercooled water. The presentation will include examples from different cases and discussion of the results.

Heymsfield, Gerald↗

Use of Dual-wavelength Radar for Snow Parameter Estimates

Use of dual-wavelength radar, with properly chosen wavelengths, will significantly lessen the ambiguities in the retrieval of microphysical properties of hydrometeors. In this paper, a dual-wavelength algorithm is described to estimate the characteristic parameters of the snow size distributions. An analysis of the computational results, made at X and Ka bands (T-39 airborne radar) and at S and X bands (CP-2 ground-based radar), indicates that valid estimates of the median volume diameter of snow particles, D(sub 0), should be possible if one of the two wavelengths of the radar operates in the non-Rayleigh scattering region. However, the accuracy may be affected to some extent if the shape factors of the Gamma function used for describing the particle distribution are chosen far from the true values or if cloud water attenuation is significant. To examine the validity and accuracy of the dual-wavelength radar algorithms, the algorithms are applied to the data taken from the Convective and Precipitation-Electrification Experiment (CaPE) in 1991, in which the dual-wavelength airborne radar was coordinated with in situ aircraft particle observations and ground-based radar measurements. Having carefully co-registered the data obtained from the different platforms, the airborne radar-derived size distributions are then compared with the in-situ measurements and ground-based radar. Good agreement is found for these comparisons despite the uncertainties resulting from mismatches of the sample volumes among the different sensors as well as spatial and temporal offsets.

Liao, Liang↗

On Study of Air/Space-borne Dual-Wavelength Radar for Estimates of Rain Profiles

In this study, a framework is discussed to apply air/space-borne dual-wavelength radar for the estimation of characteristic parameters of hydrometeors. The focus of our study is on the Global Precipitation Measurements (GPM) precipitation radar, a dual-wavelength radar that operates at Ku (13.8 GHz) and Ka (35 GHz) bands. As the droplet size distributions (DSD) of rain are expressed as the Gamma function, a procedure is described to derive the median volume diameter (D(sub 0)) and particle number concentration (N(sub T)) of rain. The correspondences of an important quantity of dual-wavelength radar, defined as deferential frequency ratio (DFR), to the D(sub 0) in the melting region are given as a function of the distance from the 0 C isotherm. A self-consistent iterative algorithm that shows a promising to account for rain attenuation of radar and infer the DSD without use of surface reference technique (SRT) is examined by applying it to the apparent radar reflectivity profiles simulated from the DSD model and then comparing the estimates with the model (true) results. For light to moderate rain the self-consistent rain profiling approach converges to unique and correct solutions only if the same shape factors of Gamma functions are used both to generate and retrieve the rain profiles, but does not converges to the true solutions if the DSD form is not chosen correctly. To further examine the dual-wavelength techniques, the self-consistent algorithm, along with forward and backward rain profiling algorithms, is then applied to the measurements taken from the 2nd generation Precipitation Radar (PR-2) built by Jet Propulsion Laboratory. It is found that rain profiles estimated from the forward and backward approaches are not sensitive to shape factor of DSD Gamma distribution, but the self-consistent method is.

Liao, Liang↗

Scattering Properties and Brightness Temperatures Associated with Solid Precipitation

In the past few years, early solid precipitation detection and retrieval algorithms have been developed and shown to be applicable for snowing clouds and blizzards. NOAA has an operational snow versus rain classifier based on AMSU-B observations. Solid precipitation retrieval algorithms reported in the literature over the past two years include those that rely on neural nets, statistics, or physical relationships. All of the algorithms require the use of millimeter-wave radiometer observations. The millimeter-wave frequencies are especially sensitive to the scattering and emission properties of frozen particles due to the ice particle refractive index. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. This investigation probes the sensitivity of scattering and absorption coefficients, and hence computed brightness temperatures, resulting from variations in solid precipitation cloud profiles. The first study compares the single scattering, absorption, and asymmetry parameters associated with snow particles in clouds. Several methodologies are used to convert the physical characteristics (e.g., shape, size distributions, ice-air-water ratios) of ice particles to electromagnetic properties (e.g., absorption, scattering, and asymmetry factors). These methodologies include: conversion to solid ice particles, homogeneous dielectric mixing, or discrete dipole approximation. Changes in the conversion methodology can produce computed brightness temperature differences greater than 50 Kelvin.

Skofronick-Jackson, Gail M.↗

Simulated convective systems using a cloud resolving model: Impact of large-scale temperature and moisture forcing using observations and GEOS-3 reanalysis

The GCE (Goddard Cumulus Ensemble) model, which has been developed and improved at NASA Goddard Space Flight Center over the past two decades, is considered as one of the finer and state-of-the-art CRMs (Cloud Resolving Models) in the research community. As the chosen CRM for a NASA Interdisciplinary Science (IDS) Project, GCE has recently been successfully upgraded into an MPI (Message Passing Interface) version with which great improvement has been achieved in computational efficiency, scalability, and portability. By basically using the large-scale temperature and moisture advective forcing, as well as the temperature, water vapor and wind fields obtained from TRMM (Tropical Rainfall Measuring Mission) field experiments such as SCSMEX (South China Sea Monsoon Experiment) and KWAJEX (Kwajalein Experiment), our recent 2-D and 3-D GCE simulations were able to capture detailed convective systems typical of the targeted (simulated) regions. The GEOS-3 [Goddard EOS (Earth Observing System) Version-3] reanalysis data have also been proposed and successfully implemented for usage in the proposed/performed GCE long-term simulations (i.e., aiming at producing massive simulated cloud data -- Cloud Library) in compensating the scarcity of real field experimental data in both time and space (location). Preliminary 2-D or 3-D pilot results using GEOS-3 data have generally showed good qualitative agreement (yet some quantitative difference) with the respective numerical results using the SCSMEX observations. The first objective of this paper is to ensure the GEOS-3 data quality by comparing the model results obtained from several pairs of simulations using the real observations and GEOS-3 reanalysis data. The different large-scale advective forcing obtained from these two kinds of resources (i.e., sounding observations and GEOS-3 reanalysis) has been considered as a major critical factor in producing various model results. The second objective of this paper is therefore to investigate and present such an impact of large-scale forcing on various modeled quantities (such as hydrometeors, rainfall, and etc.). A third objective is to validate the overall GCE 3-D model performance by comparing the numerical results with sounding observations, as well as available satellite retrievals.

Shie, C.-L.↗

Where is Tropopause?

Much of the earth science that is being proposed for the Stratospheric Observatory for Infrared Astronomy (SOFIA) Upper-Deck Research Facility (SURF) deals with issues related to the tropopause, which will be near SOFIA'S flight level at mid-latitudes. Interpreting in situ or remote aerosol, hydrometeor, and trace gas measurements will require accurate knowledge of the tropopause location. Examples of such measurements are presented, and a brief discussion is given on the Microwave Temperature Profiler (MTP), which the earth science community has used in the past to determine the tropopause height.

in situ↗

Classification of Tropical Oceanic Precipitation using High-Altitude Aircraft Microwave and Electric Field Measurements

During the 1998 and 2001 hurricane seasons of the western Atlantic Ocean and Gulf of Mexico, the Advanced Microwave Precipitation Radiometer (AMPR), the ER-2 Doppler (EDOP) radar, and the Lightning Instrument Package (LIP) were flown aboard the NASA ER-2 high-altitude aircraft as part of the Third Convection and Moisture Experiment (CAMEX-3) and the Fourth Convection and Moisture Experiment (CAMEX-4). Several hurricanes, tropical storms, and other precipitation systems were sampled during these experiments. An oceanic rainfall screening technique has been developed using AMPR passive microwave observations of these systems collected at frequencies of 10.7, 19.35, 37.1, and 85.5 GHz. This technique combines the information content of the four AMPR frequencies regarding the gross vertical structure of hydrometeors into an intuitive and easily executable precipitation mapping format. The results have been verified using vertical profiles of EDOP reflectivity and lower-altitude horizontal reflectivity scans collected by the NOAA WP3D Orion radar. Matching the rainfall classification results with coincident electric field information collected by the LIP readily identifies convective rain regions within the precipitation fields. This technique shows promise as a real-time research and analysis tool for monitoring vertical updraft strength and convective intensity from airborne platforms such as remotely operated or uninhabited aerial vehicles. The technique is analyzed and discussed for a wide variety of precipitation types using the 26 August 1998 observations of Hurricane Bonnie near landfall.

Hood, Robbie E.↗

Use of High-resolution WRF Simulations to Forecast Lightning Threat

Recent observational studies have confirmed the existence of a robust statistical relationship between lightning flash rates and the amount of large precipitating ice hydrometeors in storms. This relationship is exploited, in conjunction with the capabilities of recent forecast models such as WRF, to forecast the threat of lightning from convective storms using the output fields from the model forecasts. The simulated vertical flux of graupel at -15C is used in this study as a proxy for charge separation processes and their associated lightning risk. Six-h simulations are conducted for a number of case studies for which three-dimensional lightning validation data from the North Alabama Lightning Mapping Array are available. Experiments indicate that initialization of the WRF model on a 2 km grid using Eta boundary conditions, Doppler radar radial velocity and reflectivity fields, and METAR and ACARS data yield the most realistic simulations. An array of subjective and objective statistical metrics are employed to document the utility of the WRF forecasts. The simulation results are also compared to other more traditional means of forecasting convective storms, such as those based on inspection of the convective available potential energy field.

McCaul, William E.↗

Variations in Upper-Tropospheric Humidity and Convective Processes as Seen from SSM/T-2

Water vapor feedback, particularly involving water vapor in the upper troposphere (UTH), is widely regarded as the process with the most potential to amplify anthropogenic greenhouse forcing. Yet, our ability to quantify observationally water vapor variations in the current climate and the relationships to convective processes remains rather crude. Remote sensing from polar orbiting instruments has played a major role in documenting UTH variability, supplementing highly undersampled and poorly calibrated rawinsonde measurements. Most of our observational understanding of UTH has come from the 6.7 micrometer channel measurements which are subject to cloud contamination uncertainties. In this work we examine UTH variations present in the Special Sensor Microwave Temperature Sounder 2 (SSM/T-2) sensors flying aboard Defense Military Satellite Program (DMSP) polar orbiting satellites during the period 1993 through 2001. We employ data from the the 183.3 +/- 1 GHz channel which is far less sensitive to cirrus than IR methods. Our primary focus is on obtaining more reliable statistics of interannual behavior; i.e. How close to constant RH are interannual variations in T2 UTH over the tropics? How do temperature and moisture variations contribute regionally? The 1997/1998 strong ENS0 warming event and adjacent cool periods provide a strong signal to study, albeit a perturbation of natural climate variability. Modeling the 183.3 GHz channel using reanalysis temperature data, but with climatological water vapor, allows us to infer the separate contribution by water vapor in the observations. In addition, frozen hydrometeors produced by deep convection are also captured in the 150 GHz oxygen channel, providing an opportUnity to relate the incidence of deep convection to water vapor variability. Our results indicate a much larger variation of 183.3 GHz brightness temperatures would be observed were it not for water vapor variations positively correlated with tropical SSTs. Comparisons are made with previous studies using both IR and microwave observations to characterize UTH response to tropical SSTs.

Robertson, Franklin R.↗

Percentage Contributions from Atmospheric and Surface Features to Computed Brightness Temperatures

Over the past few years, there has become an increasing interest in the use of millimeter-wave (mm-wave) and sub-millimeter-wave (submm-wave) radiometer observations to investigate the properties of ice particles in clouds. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. One methodology used since the 1960's to discern the relationship between the physical state observed and the brightness temperature (TB) is through the temperature weighting function profile. In this research, the temperature weighting function concept is exploited to analyze the sensitivity of various characteristics of the cloud profile, such as relative humidity, ice water path, liquid water path, and surface emissivity. In our numerical analysis, we compute the contribution (in Kelvin) from each of these cloud and surface characteristics, so that the sum of these various parts equals the computed TB. Furthermore, the percentage contribution from each of these characteristics is assessed. There is some intermingling/contamination of the contributions from various components due to the integrated nature of passive observations and the absorption and scattering between the vertical layers, but all in all the knowledge gained is useful. This investigation probes the sensitivity over several cloud classifications, such as cirrus, blizzards, light snow, anvil clouds, and heavy rain. The focus is on mm-wave and submm-wave frequencies, however discussions of the effects of cloud variations to frequencies as low as 10 GHz and up to 874 GHz will also be presented. The results show that nearly 60% of the TB value at 89 GHz comes from the earth's surface for even the heaviest blizzard snow rates. On the other hand, a significant percentage of the TB value comes from the snow in the cloud for 166, and 183 plus or minus 7 GHz for the heavy and medium snow rates. For submm-wave channels, there is no contribution from the surface because these channels cannot probe through clouds, nor normal water vapor amounts in clear air regions. This work is extremely valuable in physically-based retrieval algorithm development research.

Gail Skofronick-Jackson↗