Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cluster detection”

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 721 records · Page 40

Ion mobility spectrometry of hydrazine, monomethylhydrazine, and ammonia in air with 5-nonanone reagent gas

Hydrazine (HZ) and monomethylhydrazine (MMH) in air were monitored continuously using a hand-held ion mobility spectrometer equipped with membrane inlet, 63Ni ion source, acetone reagent gas, and ambient temperature drift tube. Response characteristics included detection limit, 6 ppb; linear range, 10-600 ppb; saturated response, >2 ppm; and stable response after 15-30 min. Ammonia interfered in hydrazines detection through a product ion with the same drift time as that for MMH and HZ. Acetone reagent gas was replaced with 5-nonanone to alter drift times of product ions and separate ammonia from MMH and HZ. Patterns in mobility spectra, ion identifications from mass spectra, and fragmentation cross-sections from collisional-induced dissociations suggest that drift times are governed by ion-cluster equilibria in the drift region of the mobility spectrometer. Practical aspects including calibration, stability, and reproducibility are reported from the use of a hand-held mobility spectrometer on the space shuttle Atlantis during mission STS-37.

Flight Experiment↗

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI.

Comparison↗

Small-scale Interstellar Structure Toward the Open Cluster Chi Persei; Program ID: BO29 (Cycle 2)

The primary objective of this program was to obtain FUSE observations of the interstellar H2 absorption toward six early-type stars in the core of the open star cluster Chi Persei. High resolution optical observations of the interstellar Na I absorption toward these stars and others in the core of Chi Per have revealed a rich variety of small-scale diffuse ISM structure in both the distant Perseus and more local Orion spiral arms. At the 1 arc minute angular scales typically probed by the separations of these stars, this structure corresponds to respective physical length scales of 0.6 and 0.15 pc in the Perseus and Orion arm gas. The principal scientific goal of the FUSE observations was to determine if this small-scale Na I structure toward the core of Chi Per was reflected in the abundance and/or excitation of interstellar H2. In September 2001, we obtained FUSE observations of the Chi Per stars BD +56 563, BD +56 571, BD +56 573, BD +56 574, BD +56 575, and BD +56 578 for this program. The data reveal that the Perseus arm gas exhibits significant small-scale variations in the strength of its molecular hydrogen absorption while the Orion arm H2 absorption is relatively uniform. We presented these results at the January 2003 American Astronomical Society meeting in Seattle, WA (points, Meyer, Lauroesch, Nguyen, and Friedman 2002, Bulletin of the American Astronomical Society, 34,1229). Since Na I is not a dominant ion in H I clouds, the variations detected in the Perseus arm Na I line profiles could be due to spatial differences in either the physical structure or the physical conditions in the gas. Although many of the H2 absorption profiles are too saturated or too weak to be definitive in this regard, it is clear that small-scale variations are present in H2 gas with excitation as high as the J=4 rotational level. In contrast, lower spatial resolution H I 21 cm studies of the Perseus arm gas across the face of the h and Chi Per double cluster have shown smoother, continuous variations on the larger scale. In order to investigate if the global velocity structure of H2 toward h and Chi Per more closely follows the chaotic structure seen in interstellar Na I absorption or the smoother, continuous H I structure, we have expanded our FUSE study with a Cycle 4 program (D902, PI: Points) to survey the H2 absorption toward a number of other stars across the face of the double cluster. We have obtained FUSE spectra of 14 additional h and Chi Per stars through this program and are now combining all of the data into a comprehensive study that will be submitted for publication next year.

Meyer, David M.↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

A simulation experiment of photochemical reactions in the mesosphere

An ionospheric simulation experiment has been performed in a large vacuum chamber. The chamber is filled with NO and other gases including N2, O2, CO2, NH3 and H2O in the pressure range of 0.01 torr. A lamp which produces photons at 1236- and 1165-A by means of microwave discharge in krypton is utilized as an ionization source. In addition to 30+ large quantities of the water cluster ions 55+, H3O(+).(H2O)2, 73+, H3O(+).(H2O)3 and 91+, H3O(+).(H2O)4 were observed when nitric oxide and water were present. This closely approximates the condition of the terrestrial D region. After long periods of UV irradiation 74+ and 104+ ions grow in intensity. These ions are tentatively identified as NO(+).N2O and NO(+).NO.N2O. In addition the series 18+, 36+, 54+, and 72+ is detected which can be labeled NH4(+), NH4(+).(H2O), NH4(+).(H2O)2 and NH4(+).(H2O)3. These same species of ions are observed with the introduction of ammonia into the chamber. Presumably both N2O and NH3 are products of the photolysis.

Fugono, N.↗

High resolution far-infrared observations of the evolved H II region M16

A balloon-borne telescope was used to map NGC 6611, and three far-IR sources embedded in an extended ridge of emission were detected. It is found that Source I is an unresolved far-IR source, embedded in a molecular cloud, which despite the association of an H2O maser with its source does not exhibit radio continuum emissions. Far-IR sources II and III are associated with ionized gas-molecular cloud interfaces, with the far-IR radiation arising from dust at the boundary heated by the OB cluster. These observations are consistent with models of the later stages of evolution of an H II region.

Mcbreen, B.↗

Lightning on Venus

On the night side of Venus, the plasma wave instrument on the Pioneer-Venus Orbiter frequently detects strong and impulsive low-frequency noise bursts when the local magnetic field is strong and steady and when the field is oriented to point down to the ionosphere. The signals have characteristics of lightning whistlers, and an attempt was made to identify the sources by tracing rays along the B-field from the Orbiter down toward the surface. An extensive data set strongly indicates a clustering of lightning sources near the Beta and Phoebe Regios, with additional significant clustering near the Atla Regio at the eastern edge of Aphrodite Terra. These results suggest that there are localized lightning sources at or near the planetary surface.

Scarf, F. L.↗

A kinematic search for supernova remnants in giant extragalactic H II regions

We have obtained velocity fields of the Giant H II complexes NGC 5471 in M101, NGC 2363 in NGC 2366, and the largest H II region in NGC 2403 from H-alpha observations using the TAURUS imaging Fabry-Perot interferometer. We have detected five H-alpha sources with velocity profiles which are broad when compared with the surrounding H II region. Region B in NGC 5471 has been previously determined to contain a supernova remnant by the presence of nonthermal radio continuum radiation and enhanced (O I) and (S II) emission (Skillman 1985) and broad H-alpha emission (Chu & Kennicutt 1986). Two broad H-alpha sources in NGC 2363 coincide with regions where strong splitting has been found in the (O III) line (Roy et al. 1991). Two more broad H-alpha sources have been identified in the largest H II region in NGC 2403. Very Large Array (VLA) radio continuum observations with a resolution of 2 sec at lambda(6) and lambda(20) cm of all 3 H II complexes are presented. In addition, high resolution (subarcsecond) VLA images of NGC 5471 were made at lambda(2) and lambda(6) cm. The presence of a nonthermal source in region NGC 5471 B was confirmed while region NGC 5471 A appears to be dominated by thermal emission. The nonthermal spectral index in NGC 2363 A indicates the existence of none or more supernova remnants at the position of a large velocity width source detected in H-alpha emission. No similar nonthermal sources were detected in NGC 2403 #1. Supernovae explosions and stellar winds are considered as causes for these large velocity width sources (LVWS). If the emission from the LVWSs is attributed to single supernova remnants, they are unusually luminous in both nonthermal radio continuum and H-alpha emision. The very large H-alpha luminosities could be a result of high velocity gas being ionized by the neighboring stellar cluster.

Yang, Hui↗

Beamed and Unbeamed X-Ray Emission in FR1 Radio Galaxies

The research exploited ROSAT's sensitivity, together with its spatial and spectral resolution, to separate X-ray emission components in the sources. Prior to ROSAT, the dominant X-ray emission mechanism in radio galaxies as a class was unclear, with correlations between the X-ray and radio emission used on one hand to argue for a nuclear origin for the X-rays, and on the other hand for a thermal origin. Our observations (normally between 10 and 25 ks in length) routinely detected the target sources, and demonstrated that both resolved (thermal) and unresolved X-ray emission are typically present. Highlights of our work included two of the first detections of high-power radio galaxies at high redshift, 3C 280 and 3C 220.1. When combined with the work of two other groups, we find that of the 38 radio galaxies at z > 0.6 in the 3CRR sample, 12 were observed in ROSAT pointed observations and 9 were detected with the four most significant detections exhibiting source extent, including 3C 280 and 3C 220.1. Moreover, we discovered extended emission around five 3CRR quasars at redshift greater than about 0.4, one of which is at z > 0.6. Unification predicts that the X-ray environments of powerful radio galaxies and quasars should be similar, and our results show that powerful radio sources are finding some of the highest-redshift X-ray clusters known to date, pointing to deep gravitational potential wells early in the Universe.

Worrall, Diana M.↗

High-resolution imaging of globular cluster cores

An approach based on the maximum entropy method aimed at increasing angular resolution to study globular cluster cores is presented. To perform the image restoration the Gull-Skilling (1989) MEMSYS-3 code for maximum entropy reconstruction of arbitrary data sets was used. This software was recently applied to restoration of ESO images of the R136 object in the core of the 30 Doradus nebula. It was demonstrated that the software made it possible to restore an image at subpixel spatial scales which facilitates the detection of very high-resolution structure in the restored image.

Weir, N.↗

A reprocessing model for the ultraviolet and optical light from 4U 1820-30

We show that the recently discovered optical and ultraviolet light from the X-ray burst source 4U 1820-30 in the globular cluster NGC 6624 is due to reprocessing of the X-rays in the outer regions of an optically thick, geometrically thin accretion disk. We suggest that observation of orbital modulation of the reprocessed light, due to the variable contribution made as the heated face of the companion turns toward and away from the observer, would provide constraints on the inclination of the binary orbit, and we suggest that detection in the reprocessed flux of the 'red noise' already observed in the X-rays would provide useful constraints on the geometry and physics of the accretion disk.

Arons, Jonathan↗

Clusters of DNA induced by ionizing radiation: formation of short DNA fragments. I. Theoretical modeling

We have developed a general theoretical model for the interaction of ionizing radiation with chromatin. Chromatin is modeled as a 30-nm-diameter solenoidal fiber comprised of 20 turns of nucleosomes, 6 nucleosomes per turn. Charged-particle tracks are modeled by partitioning the energy deposition between primary track core, resulting from glancing collisions with 100 eV or less per event, and delta rays due to knock-on collisions involving energy transfers >100 eV. A Monte Carlo simulation incorporates damages due to the following molecular mechanisms: (1) ionization of water molecules leading to the formation of OH, H, eaq, etc.; (2) OH attack on sugar molecules leading to strand breaks: (3) OH attack on bases; (4) direct ionization of the sugar molecules leading to strand breaks; (5) direct ionization of the bases. Our calculations predict significant clustering of damage both locally, over regions up to 40 bp and over regions extending to several kilobase pairs. A characteristic feature of the regional damage predicted by our model is the production of short fragments of DNA associated with multiple nearby strand breaks. The shapes of the spectra of DNA fragment lengths depend on the symmetries or approximate symmetries of the chromatin structure. Such fragments have subsequently been detected experimentally and are reported in an accompanying paper (B. Rydberg, Radiat, Res. 145, 200-209, 1996) after exposure to both high- and low-LET radiation. The overall measured yields agree well quantitatively with the theoretical predictions. Our theoretical results predict the existence of a strong peak at about 85 bp, which represents the revolution period about the nucleosome. Other peaks at multiples of about 1,000 bp correspond to the periodicity of the particular solenoid model of chromatin used in these calculations. Theoretical results in combination with experimental data on fragmentation spectra may help determine the consensus or average structure of the chromatin fibers in mammalian DNA.

Non-NASA Center↗

Changes induced on the surfaces of small Pd clusters by the thermal desorption of CO

The stability and adsorption/desorption properties of supported Pd crystallites less than 5 nm in size were studied by Auger electron spectroscopy and repeated flash thermal desorption of CO. The Pd particles were grown epitaxially on heat-treated, UHV-cleaved mica at a substrate temperature of 300 C and a Pd impingement flux of 10 to the 13th atoms/sq cm s. Auger analysis allowed in situ measurement of relative particle dispersion and contamination, while FTD monitored the CO desorption properties. The results show that significant changes in the adsorption properties can be detected. Changes in the Pd Auger signal and the desorption spectrum during the first few thermal cycles are due to particle coalescence and facetting and the rate of this change is dependent on the temperature and duration of the desorption. Significant reductions in the amplitude of the desorptions peak occur during successive CO desorptions which are attributed to increases of surface carbon, induced by the desorption of CO. The contamination process could be reversed by heat treatment in oxygen or hydrogen

Doering, D. L.↗

Current star formation in S0 galaxies: NGC 4710

Elliptical (E) and lenticular (S0) galaxies lack the substantial interstellar medium (ISM) found in the star-forming spiral galaxies. However, significant numbers of E and S0 galaxies are known to contain detectable amounts of interstellar matter (e.g., Jura 1988). Thus, it is worth investigating whether these galaxies are currently able to form stars from their ISM, or whether they should be consigned to the dustbin of inert objects (Thronson and Bally 1987). The results strongly imply that current star formation is responsible for NGC 4710's far infrared and radio continuum properties. If this is indeed the case, then one expects this star formation to be fueled by molecular gas, which is presumably dominated by H2 and can be traced by the CO-12 J=1 to 0 line. Both Kenney and Young (1988) and Sage and Wrobel (1989) have detected such an emission line from NGC 4710, and infer the presence of more than 10(exp 8) solar mass of H2. The origin of the molecular gas in NGC 4710 remains a mystery. The galaxy is very deficient in HI (Kenney and Young, in preparation), suggesting that it originally was a spiral galaxy from which the outer, mainly atomic, gas was stripped by the ram pressure of the Virgo Cluster's intracluster medium, leaving only a central interstellar medium (ISM) rich in molecular gas. Alternatively, the CO may have originated via stellar mass loss with subsequent cooling, cooling flows, or capture from a gas-rich companion. Information on the morphology and kinematics of the CO can be compared with that of the galaxy's other gases and stars to distinguish among these various possible origins for the molecular gas. Major axis CO mapping with single dishes indicate an unresolved source. Thus, a millimeter array is currently being used to image NGC 4710 in CO to provide the needed morphological and kinematical data.

Wrobel, J. M.↗

Software Testing and Verification in Climate Model Development

Over the past 30 years most climate models have grown from relatively simple representations of a few atmospheric processes to a complex multi-disciplinary system. Computer infrastructure over that period has gone from punch card mainframes to modem parallel clusters. Model implementations have become complex, brittle, and increasingly difficult to extend and maintain. Existing verification processes for model implementations rely almost exclusively upon some combination of detailed analysis of output from full climate simulations and system-level regression tests. In additional to being quite costly in terms of developer time and computing resources, these testing methodologies are limited in terms of the types of defects that can be detected, isolated and diagnosed. Mitigating these weaknesses of coarse-grained testing with finer-grained "unit" tests has been perceived as cumbersome and counter-productive. In the commercial software sector, recent advances in tools and methodology have led to a renaissance for systematic fine-grained testing. We discuss the availability of analogous tools for scientific software and examine benefits that similar testing methodologies could bring to climate modeling software. We describe the unique challenges faced when testing complex numerical algorithms and suggest techniques to minimize and/or eliminate the difficulties.

Clune, Thomas L.↗

A Survey of Variable Extragalactic Sources with XTE's All Sky Monitor (ASM)

The original goal of the project was the near real-time detection of AGN utilizing the SSC 3 of the ASM on XTE which does a deep integration on one 100 square degree region of the sky. While the SSC never performed sufficiently well to allow the success of this goal, the work on the project has led to the development of a new analysis method for coded aperture systems which has now been applied to ASM data for mapping regions near clusters of galaxies such as the Perseus Cluster and the Coma Cluster. Publications are in preparation that describe both the new method and the results from mapping clusters of galaxies.

Jernigan, Garrett↗

A Chandra X-ray Study of Cygnus A: The Nucleus - 2

We report Chandra Advanced CCD Imaging Spectrometer and quasi-simultaneous Rossi X-Ray Timing Explorer (RXTE) observations of the nearby, powerful radio galaxy Cygnus A, with the present paper focusing on the properties of the active nucleus. In the Chandra observation, the hard (less than a few keV) X-ray emission is spatially unresolved with a size is approximately 1" (1.5 kpc, H(sub 0) = 50 km/s/Mpc) and coincides with the radio and near-infrared nuclei. In contrast, the soft (less than 2 keV) emission exhibits a bipolar nebulosity that aligns with the optical bipolar continuum and emission-line structures and approximately with the radio jet. In particular, the soft X-ray emission corresponds very well with the [O III] (lambda)5007 and H(alpha) + [N II] lambda(lambda)6548, 6583 nebulosity imaged with Hubble Space Telescope. At the location of the nucleus, there is only weak soft X-ray emission, an effect that may be intrinsic or result from a dust lane that crosses the nucleus perpendicular to the source axis. The spectra of the various X-ray components have been obtained by simultaneous fits to the six detectors. The compact nucleus is detected to 100 keV and is well described by a heavily absorbed power-law spectrum with Gamma(sub h) = 1.52(sup + 0.12, sub -0.12) (similar to other 0.12 narrow-line radio galaxies) and equivalent hydrogen column N(sub H)(nuc) = 2.0(sup +0.1, sub -0.1) x 10(exp 23)/sq cm. This 0.2 column is compatible with the dust obscuration to the near-infrared source for a normal gas-to-dust ratio. The soft (less than 2 keV) emission from the nucleus may be described by a power-law spectrum with the same index (i.e., Gamma(sub l) = Gamma(sub h), although direct fits suggest a slightly larger value for Gamma(sub l). Narrow emission lines from highly ionized neon and silicon, as well as a "neutral" Fe K(alpha) line, are detected in the nucleus and its vicinity (r approximately less than 2 kpc). The equivalent width (EW) of the Fe K(alpha) line (182(sup +40, sub -54) eV) is in good agreement with theoretical predictions for the EW versus N(sub H)(nuc) relationship in various geometries. An Fe K edge is also seen. The RXTE observations indicate a temperature of kT = 6.9(sup +0., sub -1.0) keV for the cluster gas (discussed in Paper III of this series) and cluster emission lines of Fe K(alpha) and Fe K(beta) and/or Ni K(alpha). We consider the possibility that the extended soft X-ray emission is electron-scattered nuclear radiation. Given that 1% of the unabsorbed 2 - 10 keV nuclear radiation would have to be scattered, the necessary gas column [N(sub H)(Scattering) approx. = 3.5 x 10(exp 22)/sq cm] would absorb the X-rays rather than scatter them if the gas is cold. Thus, the scattering plasma must be highly ionized. If this ionization is achieved through photoionization by the nucleus, the ionization parameter zeta greater than 1 ergs cm/s and the electron density n(sub e) approx. = 6 cc given the observed distance of the soft X-ray emission from the nucleus. The electron column density inferred from the X-ray observations is much too low to account for the extended optical scattered light, strongly suggesting that the polarized optical light is scattered by dust. The presence of highly ionized Ne lines in the soft X-ray spectrum requires 20 ergs cm/s approximately less than zeta approximately less than 300 ergs cm/s these lines may originate closer to the nucleus than the extended soft continuum or in a lower density gas. A collisionally ionized thermal model of the extended soft X-rays cannot be ruled out but is unattractive in view of the low metal abundance required (Z = 0.03 Z(mass)). The hard X-ray to far-infrared ratio for the nucleus of Cygnus A is similar to that seen in Seyfert 1 and unobscured radio galaxies. By means of the correlation between hard X-ray luminosity and nuclear optical absolute magnitude for these classes of object, we estimate M(sub B) = -22.4 for Cygnus A, near the .borderline between Seyfert galaxies and QSOs.

Young, Andrew J.↗

The application of Landsat data to mapping avalanche hazards

Two test areas, representing a variety of avalanche hazards, were selected in the San Juan Mountains of Colorado. Midwinter Landsat digital data were analyzed using a clustering technique, and the results compared to 1:24,000 scale maps of avalanche hazards derived from air photo interpretation and field surveys. Confined avalanches were readily identified because of the high contrast between the snow covered avalanche track and the adjacent forested slopes. Unconfined avalanches could not be identified without supplementary topographic data. Spatial characteristics were of primary importance in delineating avalanche tracks. Spatial resolution was the limiting factor in avalanche detection. Landsat data should prove useful for rapid reconnaissance mapping of avalanche hazards, particularly in the absence of other data sources.

Waterman, S.↗