Search NASA⌕ Search

SEARCH · Search NASA

Results for “Observation Data”

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,711 records · Page 95

The JOCR program

Observational data on large-scale plasma structures in comets is of value in (1) analyzing the interaction between the solar wind and comets, (2) using comets as solar wind probes, and (3) using comets as an astrophysical plasma laboratory with the CO+ plasma serving as tracers of the magnetic field. To obtain this data, the Joint Observatory for Cometary Research was built on South Baldy in the Magdalena Mountains west of Socorro, New Mexico, at an altitude of 10,615 feet. The telescope is a 14" f/2 Schmidt camera which records an 80 deg x 10 deg field onto 4 x 5 inch plates. The camera is designed for fixed focus operation, i.e., within design limits, the focus is not a function of temperature. A vacuum platen system allows the use of film for color photography. Sample observations of comets Kohoutek, West, Kobayashi-Berger-Milon and Bradfield are interpreted.

Brandt, J. C.↗

Information-incorporated gene network construction with FDR control

Abstract Motivation Large-scale gene expression studies allow gene network construction to uncover associations among genes. To study direct associations among genes, partial correlation-based networks are preferred over marginal correlations. However, FDR control for partial correlation-based network construction is not well-studied. In addition, currently available partial correlation-based methods cannot take existing biological knowledge to help network construction while controlling FDR. Results In this paper, we propose a method called Partial Correlation Graph with Information Incorporation (PCGII). PCGII estimates partial correlations between each pair of genes by regularized node-wise regression that can incorporate prior knowledge while controlling the effects of all other genes. It handles high-dimensional data where the number of genes can be much larger than the sample size and controls FDR at the same time. We compare PCGII with several existing approaches through extensive simulation studies and demonstrate that PCGII has better FDR control and higher power. We apply PCGII to a plant gene expression dataset where it recovers confirmed regulatory relationships and a hub node, as well as several direct associations that shed light on potential functional relationships in the system. We also introduce a method to supplement observed data with a pseudogene to apply PCGII when no prior information is available, which also allows checking FDR control and power for real data analysis. Availability and implementation R package is freely available for download at https://cran.r-project.org/package=PCGII.

59 BASIC BIOLOGICAL SCIENCES↗

Estimation of Surface Sensible Heat Flux due to Precipitation over CONUS and Its Impact on Urban Extreme Precipitation Modeling

The surface sensible heat flux induced by precipitation (Q P ) is a consequence of the temperature difference between the surface and the rain droplets. Despite its seemingly negligible nature, Q P is frequently omitted from both meteorological and climatological models. Nevertheless, it is important to acknowledge the numerous occasions in which the instantaneous values of Q P can be significant, particularly during extreme precipitation events. This study undertakes a comprehensive assessment of Q P across the contiguous United States (CONUS) utilizing high-resolution reanalysis, observational data, and numerical modeling to examine the influence of Q P on precipitation and the surface energy budget. The findings indicate that the spatial distribution of Q P climatology is analogous to that of precipitation, with magnitudes ranging from 2 to 3 W m -2 predominantly over the Midwest and Southeast regions. A seasonal analysis of QP reveals that the highest values occurring during the June–August (JJA) period, averaging 3.18 W m -2 . Peak Q P values of approximately 4 W m -2 are observed during JJA over the Great Plains region. We hypothesize that the Q P during an extreme precipitation event would be nonnegligible and have a significant impact on the local weather. To test this conjecture, we perform high-resolution simulations with and without Q P during an extreme precipitation event over the Chicago Metropolitan Area (CMA). The results show that the Q P may be a dominant factor compared to other components of surface heat flux during the zenith of precipitation hours. Also, Q P has the potential to not only diminish precipitation but also alter and reconfigure the remaining surface energy budget components.

54 ENVIRONMENTAL SCIENCES↗

Frequency shifts of ionospheric nfH resonances.

New observational data are analyzed to interpret mechanisms responsible for large positive frequency shifts of Alouette II plasma resonances (corresponding to the first and second electron cyclotron harmonics) relative to frequency values expected from model magnetic field calculations. It is shown that the frequency shifts can be consistently explained by plasma wave dispersion effects combined with sounder transmitter frequency deviation (positive offset of several kilohertz) and a negative offset (several tens of gammas) in the geomagnetic field relative to the model field. Plasma wave dispersion effects are observed on the electron cyclotron second harmonic resonance when it is in the vicinity of the resonance observed near the upper hybrid frequency. The observations suggest that an oblique echo model may be required for interpretation of the electron cyclotron second harmonic resonance.

Benson, R. F.↗

A Viewing Angle-Kinetic Luminosity Scheme for BL Lacertae Objects

We propose a unified classification for BL Lacertae objects (BLs), focusing on the synchrotron peak frequency of the spectral energy distribution. The unification scheme is based on the angle theta that describes the orientation of the relativistic jet and on the electron kinetic luminosity Lamba(sub kin) of the jet. We assume that Lamba(sub kin) scales with the size of the jet (r) in a self-similar fashion (Lamba(sub kin) is proportional to r(exp 2)), as supported by observational data. The jets are self-similar in geometry and have the same pressure and median magnetic field at the inlet, independent of size. The self-similarity is broken for the highest energy electrons, which radiate mainly at high frequencies, since for large sources they suffer more severe radiative energy losses over a given fraction of the jet length. We calculate the optically thin synchrotron spectrum using an accelerating inner jet model based on simple relativistic gas dynamics and show that it can fit the observed infrared-to-X-ray spectrum of PKS 2155-304. We couple the accelerating jet model to the unification scheme and compare the results to complete samples of BLs. The negative apparent evolution of X-ray-selected BLs is explained as a result of positive evolution of the jet electron kinetic luminosity Lamba(sub kin). We review observational arguments in favor of the existence of scaled-down accretion disks and broad emission-line regions in BLs. The proposed unification scheme can explain the lack of observed broad emission lines in X-ray-selected BLs as well as the existence of those lines preferentially in luminous radio-selected BLs. Finally, we review observational arguments that suggest the extension of this unification scheme to all blazars.

Georganopoulos, Markos↗

Use of an Existing Airborne Radon Data Base in the Verification of the NASA/AEAP Core Model

The primary objective of this project was to apply the tropospheric atmospheric radon (Rn222) measurements to the development and verification of the global 3-D atmospheric chemical transport model under development by NASA's Atmospheric Effects of Aviation Project (AEAP). The AEAP project had two principal components: (1) a modeling effort, whose goal was to create, test and apply an elaborate three-dimensional atmospheric chemical transport model (the NASA/AEAP Core model to an evaluation of the possible short and long-term effects of aircraft emissions on atmospheric chemistry and climate--and (2) a measurement effort, whose goal was to obtain a focused set of atmospheric measurements that would provide some of the observational data used in the modeling effort. My activity in this project was confined to the first of these components. Both atmospheric transport and atmospheric chemical reactions (as well the input and removal of chemical species) are accounted for in the NASA/AEAP Core model. Thus, for example, in assessing the effect of aircraft effluents on the chemistry of a given region of the upper troposphere, the model must keep track not only of the chemical reactions of the effluent species emitted by aircraft flying in this region, but also of the transport into the region of these (and other) species from other, remote sources--for example, via the vertical convection of boundary layer air to the upper troposphere. Radon, because of its known surface source and known radioactive half-life, and freedom from chemical production or loss, and from removal from the atmosphere by physical scavenging, is a recognized and valuable tool for testing the transport components of global transport and circulation models.

Kritz, Mark A.↗

Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population

With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.

Andre Horstmann↗

PHOTOSPHERIC PROPERTIES OF WARM EUV LOOPS AND HOT X-RAY LOOPS

We investigate the photospheric properties (vector magnetic fields and horizontal velocity) of a well-developed active region, NOAA AR 10978, using the Hinode Solar Optical Telescope specifically to determine what gives rise to the temperature difference between ''warm loops'' (1-2 MK), which are coronal loops observed in EUV wavelengths, and ''hot loops'' (>3 MK), coronal loops observed in X-rays. We found that outside sunspots, the magnetic filling factor in the solar network varies with location and is anti-correlated with the horizontal random velocity. If we accept that the observed magnetic features consist of unresolved magnetic flux tubes, this anti-correlation can be explained by the ensemble average of flux-tube motion driven by small-scale random flows. The observed data are consistent with a flux tube width of ∼77 km and horizontal flow at ∼2.6 km s{sup –1} with a spatial scale of ∼120 km. We also found that outside sunspots, there is no significant difference between warm and hot loops either in the magnetic properties (except for the inclination) or in the horizontal random velocity at their footpoints, which are identified with the Hinode X-Ray Telescope and the Transition Region and Coronal Explorer. The energy flux injected into the coronal loops by the observed photospheric motion of the magnetic fields is estimated to be 2 × 10{sup 6} erg s{sup –1} cm{sup –2}, which is the same for both warm and hot loops. This suggests that coronal properties (e.g., loop length) play a more important role in giving rise to temperature differences of active-region coronal loops than photospheric parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Spherical harmonic analysis of a synoptic climatology generated with a global general circulation model

Spherical harmonic analysis was used to analyze the observed climatological (C) fields of temperature at 850 mb, geopotential height at 500 mb, and sea level pressure. The spherical harmonic method was also applied to the corresponding "model climatological" fields (M) generated by a general circulation model, the "GISS climate model." The climate model was initialized with observed data for the first of December 1976 at 00. GMT and allowed to generate five years of meteorological history. Monthly means of the above fields for the five years were computed and subjected to spherical harmonic analysis. It was found from the comparison of the spectral components of both sets, M and C, that the climate model generated reasonable 500 mb geopotential heights. The model temperature field at 850 mb exhibited a generally correct structure. However, the meridional temperature gradient was overestimated and overheating of the continents was observed in summer.

Christidis, Z. D.↗

Fluid Stratification Separate Effects Analysis, Testing and Benchmarking

High Temperature Gas Reactors (HTGRs) are positioned to disrupt local and global markets via their unique ability to produce carbon-free process heat, high efficiency power generation, and passively safe operational features. However, significant impediments still exist to delay deployment of this particular technology, including a lack of experimental data, verified code application, and lack of consensus with regards to severe accident progression. In particular, air ingress accidents represent a particular challenge to designers and engineers, as they represent low probability, but highly complex, accident scenarios. Including phenomena such as molecular diffusion, free convection, and complex heat and mass transfer paths, experimental and traceable data is essential to maturing the state of the industry. Therefore, this work presents an experimental investigation of the transition to natural convection in HTGR applications using the Stratified Flow Separate Effects Test Facility, housed at Oregon State University. In particular, this work will present data that challenges the assumption that molecular diffusion is a significant factor in this severe accident in the reference facility of the General Atomic 600 MW th Gas Turbine-Modular Helium Reactor (GT-MHR). Rather, experimentally produced data shows a statistically suggestive retardant effect of cross duct orientation, indicating that ingress mechanics may be fundamentally altered via facility geometry. Experimentally, this is achieved using a simplified cross duct that may be positioned in one of two ways so as to provide either horizontal or vertical access to the lower plenum area. Onset of natural convection (ONC) is measured using an oxygen sensor probe, immersed in the helium working fluid, so as to provide direct indication of air presence in the upper plenum. This document will present the doctoral dissertation produced by this work, along with its analysis, experimental data, observations, and conclusions, in addition to select documents that will highlight critical changes during the development process. Those documents will include the quality assurance plan, all subsequent change requests, and the engineering transfer from Harris Thermal in order to establish the most complete quality record possible.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of integrated assessment model hindcast experiments: a case study of the GCAM 3.0 land use module

Abstract. Hindcasting experiments (conducting a model forecast for a time period in which observational data are available) are being undertaken increasingly often by the integrated assessment model (IAM) community, across many scales of models. When they are undertaken, the results are often evaluated using global aggregates or otherwise highly aggregated skill scores that mask deficiencies. We select a set of deviation-based measures that can be applied on different spatial scales (regional versus global) to make evaluating the large number of variable–region combinations in IAMs more tractable. We also identify performance benchmarks for these measures, based on the statistics of the observational dataset, that allow a model to be evaluated in absolute terms rather than relative to the performance of other models at similar tasks. An ideal evaluation method for hindcast experiments in IAMs would feature both absolute measures for evaluation of a single experiment for a single model and relative measures to compare the results of multiple experiments for a single model or the same experiment repeated across multiple models, such as in community intercomparison studies. The performance benchmarks highlight the use of this scheme for model evaluation in absolute terms, providing information about the reasons a model may perform poorly on a given measure and therefore identifying opportunities for improvement. To demonstrate the use of and types of results possible with the evaluation method, the measures are applied to the results of a past hindcast experiment focusing on land allocation in the Global Change Assessment Model (GCAM) version 3.0. The question of how to more holistically evaluate models as complex as IAMs is an area for future research. We find quantitative evidence that global aggregates alone are not sufficient for evaluating IAMs that require global supply to equal global demand at each time period, such as GCAM. The results of this work indicate it is unlikely that a single evaluation measure for all variables in an IAM exists, and therefore sector-by-sector evaluation may be necessary.

58 GEOSCIENCES↗

Search for heavy resonances decaying into two Higgs bosons or into a Higgs boson and a W or Z boson in proton-proton collisions at 13 TeV

A search is presented for massive narrow resonances decaying either into two Higgs bosons, or into a Higgs boson and a W or Z boson. The decay channels considered are $ \mathrm{H}\mathrm{H}\to \mathrm{b}\overline{\mathrm{b}}{\tau}^{+}{\tau}^{-} $ and $ \mathrm{V}\mathrm{H}\to \mathrm{q}\overline{\mathrm{q}}{\tau}^{+}{\tau}^{-} $ , where H denotes the Higgs boson, and V denotes the W or Z boson. This analysis is based on a data sample of proton-proton collisions collected at a center-of-mass energy of 13 TeV by the CMS Collaboration, corresponding to an integrated luminosity of 35.9 fb$^{−1}$. For the TeV-scale mass resonances considered, substructure techniques provide ways to differentiate among the hadronization products from vector boson decays to quarks, Higgs boson decays to bottom quarks, and quark- or gluon-induced jets. Reconstruction techniques are used that have been specifically optimized to select events in which the tau lepton pair is highly boosted. The observed data are consistent with standard model expectations and upper limits are set at 95% confidence level on the product of cross section and branching fraction for resonance masses between 0.9 and 4.0 TeV. Exclusion limits are set in the context of bulk radion and graviton models:spin-0 radion resonances are excluded below a mass of 2.7 TeV at 95% confidence level. In the spin-1 heavy vector triplet framework, mass-degenerate W′ and Z′ resonances with dominant couplings to the standard model gauge bosons are excluded below a mass of 2.8 TeV at 95% confidence level. These are the first limits for massive resonances at the TeV scale with these decay channels at $ \sqrt{s}=13 $ TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the Role of Airborne Ice Nucleating Particles in Primary and Secondary Ice Formation Processes in Convective Midlatitude Clouds

Abstract Formation of ice in mixed-phase clouds may proceed via primary nucleation and secondary ice production (SIP). Primary nucleation involves the presence of ice nucleating particles (INPs), while SIP follows primary nucleation, resulting in ice crystal number concentrations N ice exceeding the number of INPs. The ice formation pathways in young congestus clouds are not well understood nor are the environmental conditions favorable for the onset of SIP. Coincident airborne measurements of INPs and N ice collected during the Secondary Production of Ice in Cumulus Experiment (SPICULE) campaign over the U.S. central Great Plains are reported for young congestus clouds. Number concentrations of INPs and N ice at temperatures between −10° and −20°C were used to categorize whether SIP was active in each analyzed cloud pass. Observational data suggested that fragmentation of freezing droplets (FFDs) was occurring and the most likely SIP mechanism for ice enhancement during the early stages of the sampled clouds, where higher cloud-base temperatures and stronger updrafts produced conditions favorable for the onset of FFD. Simple model simulations of a congestus cloud were used to compare predictions of N ice from newly implemented SIP mechanisms (including FFD) alongside those from the existing Hallett–Mossop (HM) mechanism. Although predicted N ice , when all SIP mechanisms were active, was in good agreement with observations, the dominant SIP mechanism in the model was predicted to be HM. The rate of heterogeneous freezing in raindrops was likely inhibiting realistic FFD rates in the simulated cloud. These discrepancies between the observations and simulations underscore the need for extended laboratory studies to formulate improved model representations of SIP in convective clouds.

Patnaude, Ryan J. [Department of Atmospheric Scien↗

The Attraction Indian Buffet Distribution

We propose the attraction Indian buffet distribution (AIBD), a distribution for binary feature matrices influenced by pairwise similarity information. Binary feature matrices are used in Bayesian models to uncover latent variables (i.e., features) that explain observed data. The Indian buffet process (IBP) is a popular exchangeable prior distribution for latent feature matrices. In the presence of additional information, however, the exchangeability assumption is not reasonable or desirable. The AIBD can incorporate pairwise similarity information, yet it preserves many properties of the IBP, including the distribution of the total number of features. Thus, much of the interpretation and intuition that one has for the IBP directly carries over to the AIBD. A temperature parameter controls the degree to which the similarity information affects feature-sharing between observations. Unlike other nonexchangeable distributions for feature allocations, the probability mass function of the AIBD has a tractable normalizing constant, making posterior inference on hyperparameters straight-forward using standard MCMC methods. A novel posterior sampling algorithm is proposed for the IBP and the AIBD. We demonstrate the feasibility of the AIBD as a prior distribution in feature allocation models and compare the performance of competing methods in simulations and an application.

97 MATHEMATICS AND COMPUTING↗

The rate of binary black hole mergers inferred from advanced LIGO observations surrounding GW150914

A transient gravitational-wave signal, GW150914, was identified in the twin Advanced LIGO detectors on 2015 September 2015 at 09:50:45 UTC. To assess the implications of this discovery, the detectors remained in operation with unchanged configurations over a period of 39 days around the time of the signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of the 16 days of simultaneous two-detector observational data is estimated to have a false-alarm rate (FAR) of <4.9×10{sup −6} yr{sup −1}, yielding a p-value for GW150914 of <2×10{sup −7}. Parameter estimation follow-up on this trigger identifies its source as a binary black hole (BBH) merger with component masses (m{sub 1},m{sub 2})=(36{sub −4}{sup +5},29{sub −4}{sup +4}) M{sub ⊙} at redshift z=0.09{sub −0.04}{sup +0.03} (median and 90% credible range). Here, we report on the constraints these observations place on the rate of BBH coalescences. Considering only GW150914, assuming that all BBHs in the universe have the same masses and spins as this event, imposing a search FAR threshold of 1 per 100 years, and assuming that the BBH merger rate is constant in the comoving frame, we infer a 90% credible range of merger rates between 2--53 Gpc{sup −3} yr{sup −1} (comoving frame). Incorporating all search triggers that pass a much lower threshold while accounting for the uncertainty in the astrophysical origin of each trigger, we estimate a higher rate, ranging from 13--600 Gpc{sup −3} yr{sup −1} depending on assumptions about the BBH mass distribution. All together, our various rate estimates fall in the conservative range 2--600 Gpc{sup −3} yr{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of Antarctic crust and upper mantle using MAGSAT and other geophysical data

Antarctica is the perfect proving ground for testing data reduction procedures which are alternatives to the standard procedure employing polynomial fitting. Unreduced data (observed values minus corefield model theo each pass) were averaged in 3 degree bins. The resulting map was then high pass filtered using a finite Fourier transform filter so that spectral peaks between 4200 km and 5280 km were diminished by 1/3, peaks corresponding to wavelengths greater than or equal to 5280 km were diminished by 1/2, and the d.c. component was set equal to zero. The surface subtracted in this way differs from the polynomials used in the standard method in that it is static over the 5 north data window. The map produced is presented and the effects of the static model are compared to those of the dynamic mode.

Bentley, C. R.↗

Empirical and modeled synoptic cloud climatology of the Arctic Ocean

A set of cloud cover data were developed for the Arctic during the climatically important spring/early summer transition months. Parallel with the determination of mean monthly cloud conditions, data for different synoptic pressure patterns were also composited as a means of evaluating the role of synoptic variability on Arctic cloud regimes. In order to carry out this analysis, a synoptic classification scheme was developed for the Arctic using an objective typing procedure. A second major objective was to analyze model output of pressure fields and cloud parameters from a control run of the Goddard Institue for Space Studies climate model for the same area and to intercompare the synoptic climatatology of the model with that based on the observational data.

Barry, R. G.↗

An experimental study of counter-rotating cores in elliptical galaxies

Recent observational studies (Franx and Illingworth 1987; Jedrzejewski and Schechter 1988; Bender 1988; Illingworth and Franx 1989) have shown that some elliptical galaxies have a small region near the center that rotates in the opposite direction from the outer parts of the galaxy. Often the rotation in the central part is much faster than that in the outer part. A few other galaxies show a small region near the center that rotates in the same direction as the rest of the galaxy, but much faster. Either way, the part near the center that shows a strange pattern of rotation (the 'core') has been interpreted as a distinct dynamical subsystem. Very briefly, the observational data seem to be that anomalies show up in rotation curves near the centers of some elliptical galaxies and that galaxies with these strange rotational properties do not show a photometric signature: there are no noticeable bumps in the brightness profile and no unusual shapes of isophotal contours that would suggest an excess of matter concentrated near the center. No strong color variations have been reported. The puzzle is to learn what we can about elliptical galaxies in general, and about galaxies with strange central regions in particular, from these observational facts. The authors' approach is experimental. They make a guess about the form of the dynamically distinct subsystem, and then build a galaxy model to test experimental consequences such as the amount of matter required to produce observable effects and the length of time over which these effects would remain observable. They sidestep questions about how the galaxy might have gotten to be that way in the first place. That gives them more freedom to explore a variety of suggestions about what kind of dynamical system might give rise to the observed rotational patterns.

Miller, R. H.↗