Search NASA⌕ Search

SEARCH · Search NASA

Results for “Resolution refinement”

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 631 records · Page 35

A New Glacial Varve Chronology Along the Southern Laurentide Ice Sheet That Spans the Younger Dryas–Holocene Boundary

Glacial varves can detail ice-margin positions and provide a proxy for meltwater discharge at resolutions comparable to those of the Greenland ice core archives, and thus they can be critical paleorecords for assessing the response of both ancient and modern ice sheets to climate change. Here we provide an ~1500 yr varve chronology straddling the Younger Dryas (YD)–Holocene boundary (11.65 cal. kyr B.P.), the first such chronology in North America. The varves are from glacial Lake Agassiz (central North America). The chronology is pinned on accelerator mass spectrometry radiocarbon-dated terrestrial macrofossils at the base of a widespread red-clay bed deposited during flooding from the Lake Superior basin. We illustrate the utility of this record by examining ice-margin retreat and melting through the late Younger Dryas and across the Holocene boundary. The ice margin receded at a constant rate, not only during the late YD, but for at least 300 yr after the onset of the Holocene. In contrast, varve thicknesses increased at the boundary, and a moraine formed over a 50 yr period, perhaps in response to the warming climate. Our expectation is that this time series will continue to be developed, expanded, and refined because it promises to be the longest and most geographically extensive glacial varve data set in North America.

Glacial varves↗

High-Resolution, Lightweight, and Low-cost X-Ray Optics for the Lynx Observatory

We describe an approach to build an x-ray mirror assembly that can meet Lynx’s requirements of high-angular resolution, large effective area, light weight, short production schedule, and low-production cost. Adopting a modular hierarchy, the assembly is composed of 37,492 mirror segments, each of which measures ∼100 mm × 100 mm × 0.5 mm. These segments are integrated into 611 modules, which are individually tested and qualified to meet both science performance and spaceflight environment requirements before they in turn are integrated into 12 metashells. The 12 metashells are then integrated to form the mirror assembly. This approach combines the latest precision polishing technology and the monocrystalline silicon material to fabricate the thin and lightweight mirror segments. Because of the use of commercially available equipment and material and because of its highly modular and hierarchical building-up process, this approach is highly amenable to automation and mass production to maximize production throughput and to minimize production schedule and cost. As of fall 2018, the basic elements of this approach, including substrate fabrication, coating, alignment, and bonding, have been validated by the successful building and testing of single-pair mirror modules. In the next few years, the many steps of the approach will be refined and perfected by repeatedly building and testing mirror modules containing progressively more mirror segments to fully meet science performance, spaceflight environments, as well as programmatic requirements of the Lynx mission and other proposed missions, such as AXIS.

Zhang, William W.↗

Fast event-based electron counting for small-molecule structure determination by MicroED

Electron counting helped realize the resolution revolution in single-particle cryoEM and is now accelerating the determination of MicroED structures. Its advantages are best demonstrated by new direct electron detectors capable of fast (kilohertz) event-based electron counting (EBEC). This strategy minimizes the inaccuracies introduced by coincidence loss (CL) and promises rapid determination of accurate structures. We used the Direct Electron Apollo camera to leverage EBEC technology for MicroED data collection. Given its ability to count single electrons, the Apollo collects high-quality MicroED data from organic small-molecule crystals illuminated with incident electron beam flux densities as low as 0.01–0.045 e − /Å 2 /s. Under even the lowest flux density (0.01 e − /Å 2 /s) condition, fast EBEC data produced ab initio structures of a salen ligand (268 Da) and biotin (244 Da). Each structure was determined from a 100° wedge of data collected from a single crystal in as few as 50 s, with a delivered fluence of only ∼0.5 e − /Å 2 . Fast EBEC data collected with a fluence of 2.25 or 3.33 e − /Å 2 also facilitated a 1.5 Å structure of thiostrepton (1665 Da). While refinement of these structures appeared unaffected by CL, a CL adjustment applied to EBEC data further improved the distribution of intensities measured from the salen ligand and biotin crystals. However, CL adjustment only marginally improved the refinement of their corresponding structures, signaling the already high counting accuracy of detectors with counting rates in the kilohertz range. Overall, by delivering low-dose structure-worthy data, fast EBEC collection strategies open new possibilities for high-throughput MicroED.

EBEC↗

The Wide-Field Imaging Interferometry Testbed: Progress and Plans

We describe the technique of wide field mosaic imaging for optical/IR interferometers and present early experimental results from a laboratory instrument designed to validate, experiment with, and refine the technique. A conventional single-detector stellar interferometer operating with narrow bandwidth at center wavelength lambda is limited in its field of view to the primary beam of the individual telescope apertures, or approx. lambda/D(sub tel) radians, where is the telescope diameter. Such a field is too small for many applications; often one wishes to image extended sources. We are developing and testing a technique analogous to the mosaic method employed in millimeter and radio astronomy, but applicable to optical/IR Michelson interferometers, in which beam combination is done in the pupil plane. An N(sub pix) x N(sub pix) detector array placed in the image plane of the interferometer is used to record simultaneously the fringe patterns from many contiguous telescope fields, effectively multiplying the field size by N(sub pix)/2, where the factor 2 allows for Nyquist sampling. This mosaic imaging technique will be especially valuable for far IR and submillimeter interferometric space observatories such as the Space Infrared Interferometric Telescope (SPIRIT) and the Submillimeter Probe of the Evolution of Cosmic Structure (SPECS). SPIRIT and SPECS will be designed to provide sensitive, high angular resolution observations of fields several arcminutes in diameter, and views of the universe complementary to those provided by HST, NGST, and ALMA.

Rinehart, S. A.↗

Hydrodynamic simulations of white dwarf–white dwarf mergers and the origin of R Coronae Borealis stars

We study the properties of double white dwarf (DWD) mergers by performing hydrodynamic simulations using the new and improved adaptive mesh refinement code OCTO-TIGER. We follow the orbital evolution of DWD systems of mass ratio $q=0.7$ for tens of orbits until and after the merger to investigate them as a possible origin for R Coronae Borealis (RCB) type stars. We reproduce previous results, finding that during the merger, the helium WD donor star is tidally disrupted within 20–80 min since the beginning of the simulation onto the accretor carbon–oxygen WD, creating a high temperature shell around the accretor. We investigate the possible helium burning in this shell and the merged object’s general structure. Specifically, we are interested in the amount of oxygen-16 dredged-up from the accretor to the hot shell and the amount of oxygen-18 produced. This is critical as the discovery of very low oxygen-16 to oxygen-18 ratios in RCB stars pointed out the merger scenario as a favourable explanation for their origin. A small amount of hydrogen in the donor may help keep the oxygen-16 to oxygen-18 ratios within observational bounds, even if moderate dredge-up from the accretor occurs. In addition, we perform a resolution study to reconcile the difference found in the amount of oxygen-16 dredge-up between smoothed-particle hydrodynamics and grid-based simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Experiments with the Mesoscale Atmospheric Simulation System (MASS) using the synthetic relative humidity

This study is intended to examine the impact of the synthetic relative humidity on the model simulation of mesoscale convective storm environment. The synthetic relative humidity is derived from the National Weather Services surface observations, and non-conventional sources including aircraft, radar, and satellite observations. The latter sources provide the mesoscale data of very high spatial and temporal resolution. The synthetic humidity data is used to complement the National Weather Services rawinsonde observations. It is believed that a realistic representation of initial moisture field in a mesoscale model is critical for the model simulation of thunderstorm development, and the formation of non-convective clouds as well as their effects on the surface energy budget. The impact will be investigated based on a real-data case study using the mesoscale atmospheric simulation system developed by Mesoscale Environmental Simulations Operations, Inc. The mesoscale atmospheric simulation system consists of objective analysis and initialization codes, and the coarse-mesh and fine-mesh dynamic prediction models. Both models are a three dimensional, primitive equation model containing the essential moist physics for simulating and forecasting mesoscale convective processes in the atmosphere. The modeling system is currently implemented at the Applied Meteorology Unit, Kennedy Space Center. Two procedures involving the synthetic relative humidity to define the model initial moisture fields are considered. It is proposed to perform several short-range (approximately 6 hours) comparative coarse-mesh simulation experiments with and without the synthetic data. They are aimed at revealing the model sensitivities should allow us both to refine the specification of the observational requirements, and to develop more accurate and efficient objective analysis schemes. The goal is to advance the MASS (Mesoscal Atmospheric Simulation System) modeling expertise so that the model output can provide reliable guidance for thunderstorm forecasting.

Chang, Chia-Bo↗

Orbit Determination Support for the Microwave Anisotropy Probe (MAP)

THe microwave Anisotropy Probe (MAP) ia the second Medium Class Explorer (MIDEX) mission of the National Aeronautics and Space Administration (NASA). The main goal of the MAP observatory is to measure the temperature fluctuations, known as anisotropy, of the cosmic microwave background (CBG) radiation over the entire sky and to produce a map of the CMB anisotropies with an angular resolution of approximately 3 degrees. MAP was launched from the Cape Canaveral Air Force Station Complex 17 aboard a Delta II 7425-10 expendable launch vehicle at exactly 19:46:46.183 UTC on June 30, 20001. The spacecraft receiver a nominal direct insertion by the Delta expendable launch vehicle into a 185-km circular orbit with a 28.7 deg. inclination. MAP was than maneuvered into a sequence of phasing loops designed to set up a lunar gravity-assisted acceleration of the spacecraft onto a transfer trajectory to a lissajous orbit about the Earth-Sun L2 Lagrange point, about 1.5 million km from Earth. The science mission minimum lifetime is two years of observations at L2 with a desired lifetime of 4 years. The MAP transfer orbit consisted of 3.5 phasing loops.The MAP trajectory schematic all the way through L2 is shown. The first loop had a period of 7 days, the second and third loops were 10 days long, and the last half loop was 5 days. The periselene (i.e., lunar encounter or swingby) took place approximately 30 days after launch. After the periselene, the spacecraft cruised for approximately 60 days before it arrived in the vicinity of the L2 libration point. Two mid-course correction (MCC) maneuvers were performed to refine MAP's post-launch trajectory-one after periselene and one prior to arrival at vicinity of L2. Now that MAP is at its operational L2 lissajous orbit, the MAP satellite is commanded to perform occasional station-keeping (SK) maneuvers in order to maintain its orbit around L2. Because of its complex orbital characteristics, the mission provided a unique challenge to orbit determination (OD) support in many orbital regimes.

Truong, Son H.↗

Comparing the interactions between particulate matter and cloud properties over two populated cities in Texas using WRF-Chem fine-resolution modeling

Accurate modeling of aerosol-cloud interactions is essential for reliable weather and air quality simulations, given their significant impact on precipitation patterns, cloud dynamics, and aerosol distributions. This study employed the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) to examine the impact of enhanced meteorological simulations, achieved through advanced microphysics parameterization supported by data assimilation techniques, on air quality across Texas on August 19 and 20, 2022. We tested four distinct configurations: (1) the Morrison two-moment bulk microphysics scheme, (2) Morrison's with observation nudging, (3) the Spectral Bin Microphysics (SBM), and (4) SBM with observation nudging. While the SBM scheme is known for its detailed representation of aerosol-cloud interactions, our focus was on how improvements in meteorological accuracy translate to more precise air quality simulations. Our findings demonstrated a progressive improvement in simulation accuracy, starting with the Morrison's scheme and further enhanced by adopting the SBM scheme, complemented by incorporating observation nudging. Specifically, the combination of the SBM scheme and the nudging substantially enhanced the model's ability to capture convective precipitation events, as shown by better alignment with NEXRAD radar reflectivity, with R increasing from –0.21 to 0.82, IOA from 0.10 to 0.87, and NMB decreasing from 99% to 34% in Houston. The enhanced meteorology translated into more accurate PM 2.5 concentration simulations, particularly through the more accurate representation of aerosol washout during precipitation events. In Houston, the SBM scheme with nudging improved the model's PM 2.5 simulations significantly, with NMB decreasing from –20% to 5% and IOA improving from 0.43 to 0.61. In San Antonio, improvements were also notable, with NMB improved from –27% to –22%, R increased from 0.48 to 0.82, and IOA increased from 0.66 to 0.86. Furthermore, our results underscore the crucial role of accurate meteorological simulations in refining our understanding of aerosol behaviors in relation to precipitation patterns, directly enhancing the reliability and effectiveness of air quality modeling.

54 ENVIRONMENTAL SCIENCES↗

Remote Sensing of In-Flight Icing Conditions: Operational, Meteorological, and Technological Considerations

Remote-sensing systems that map aircraft icing conditions in the flight path from airports or aircraft would allow icing to be avoided and exited. Icing remote-sensing system development requires consideration of the operational environment, the meteorological environment, and the technology available. Operationally, pilots need unambiguous cockpit icing displays for risk management decision-making. Human factors, aircraft integration, integration of remotely sensed icing information into the weather system infrastructures, and avoid-and-exit issues need resolution. Cost, maintenance, power, weight, and space concern manufacturers, operators, and regulators. An icing remote-sensing system detects cloud and precipitation liquid water, drop size, and temperature. An algorithm is needed to convert these conditions into icing potential estimates for cockpit display. Specification development requires that magnitudes of cloud microphysical conditions and their spatial and temporal variability be understood at multiple scales. The core of an icing remote-sensing system is the technology that senses icing microphysical conditions. Radar and microwave radiometers penetrate clouds and can estimate liquid water and drop size. Retrieval development is needed; differential attenuation and neural network assessment of multiple-band radar returns are most promising to date. Airport-based radar or radiometers are the most viable near-term technologies. A radiometer that profiles cloud liquid water, and experimental techniques to use radiometers horizontally, are promising. The most critical operational research needs are to assess cockpit and aircraft system integration, develop avoid-and-exit protocols, assess human factors, and integrate remote-sensing information into weather and air traffic control infrastructures. Improved spatial characterization of cloud and precipitation liquid-water content, drop-size spectra, and temperature are needed, as well as an algorithm to convert sensed conditions into a measure of icing potential. Technology development also requires refinement of inversion techniques. These goals can be accomplished with collaboration among federal agencies including NASA, the FAA, the National Center for Atmospheric Research, NOAA, and the Department of Defense. This report reviews operational, meteorological, and technological considerations in developing the capability to remotely map in-flight icing conditions from the ground and from the air.

Ryerson, Charles C.↗

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗

Casual Rerouting of AERONET Sun/Sky Photometers: Toward a New Network of Ground Measurements Dedicated to the Monitoring of Surface Properties?

This paper presents an innovative method for observing vegetation health at a very high spatial resolution (~5 × 5 cm) and low cost by upgrading an existing Aerosol RObotic NETwork (AERONET) ground station dedicated to the observation of aerosols in the atmosphere. This study evaluates the capability of a sun/sky photometer to perform additional surface reflectance observations. The ground station of Toulouse, France, which belongs to the AERONET sun/sky photometer network, is used for this feasibility study. The experiment was conducted for a 5-year period (between 2016 and 2020). The sun/sky photometer was mounted on a metallic structure at a height of 2.5 m, and the acquisition software was adapted to add a periodical (every hour) ground-observation scenario with the sun/sky photometer observing the surface instead of being inactive. Evaluation is performed by using a classical metric characterizing the vegetation health: the normalized difference vegetation index (NDVI), using as reference the satellite NDVI derived from a Sentinel-2 (S2) sensor at 10 × 10 m resolution. Comparison for the 5-year period showed good agreement between the S2 and sun/sky photometer NDVIs (i.e., bias = 0.004, RMSD = 0.082, and R = 0.882 for a mean value of S2A NDVI around 0.6). Discrepancies could have been due to spatial-representativeness issues (of the ground measurement compared to S2), the differences between spectral bands, and the quality of the atmospheric correction applied on S2 data (accuracy of the sun/sky photometer instrument was better than 0.1%). However, the accuracy of the atmospheric correction applied on S2 data in this station appeared to be of good quality, and no dependence on the presence of aerosols was observed. This first analysis of the potential of the CIMEL CE318 sun/sky photometer to monitor the surface is encouraging. Further analyses need to be carried out to estimate the potential in different AERONET stations. The occasional rerouting of AERONET stations could lead to a complementary network of surface reflectance observations. This would require an update of the software, and eventual adaptations of the measurement platforms to the station environments. The additional cost, based on the existing AERONET network, would be quite limited. These new surface measurements would be interesting for measurements of vegetation health (monitoring of NDVI, and also of other vegetation indices such as the leaf area and chlorophyll indices), for validation and calibration exercise purposes, and possibly to refine various scientific algorithms (i.e., algorithms dedicated to cloud detection or the AERONET aerosol retrieval algorithm itself). CIMEL is ready to include the ground scenario used in this study in all new sun/sky photometers.

Dominique Carrer↗

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗

Time-Dependent Response Versus Scan Angle for MODIS Reflective Solar Bands

The Moderate Resolution Imaging Spectroradiometer (MODIS) instruments currently operate onboard the National Aeronautics and Space Administration (NASA's) Terra and Aqua spacecraft, launched on December 18, 1999 and May 4, 2002, respectively. MODIS has 36 spectral bands, among which 20 are reflective solar bands (RSBs) covering a spectral range from 0.412 to 2.13 μm. The RSBs are calibrated on orbit using a solar diffuser (SD) and an SD stability monitor and with additional measurements from lunar observations via a space view (SV) port. Selected pseudo-invariant desert sites are also used to track the RSB on-orbit gain change, particularly for short-wavelength bands. MODIS views the Earth surface, SV, and the onboard calibrators using a two-sided scan mirror. The response versus scan angle (RVS) of the scan mirror was characterized prior to launch, and its changes are tracked using observations made at different angles of incidence from onboard SD, lunar, and Earth view (EV) measurements. These observations show that the optical properties of the scan mirror have experienced large wavelength-dependent degradation in both the visible and near infrared spectral regions. Algorithms have been developed to track the on-orbit RVS change using the calibrators and the selected desert sites. These algorithms have been applied to both Terra and Aqua MODIS Level 1B (L1B) to improve the EV data accuracy since L1B Collection 4, refined in Collection 5, and further improved in the latest Collection 6 (C6). In C6, two approaches have been used to derive the time-dependent RVS for MODIS RSB. The first approach relies on data collected from sensor onboard calibrators and mirror side ratios from EV observations. The second approach uses onboard calibrators and EV response trending from selected desert sites. This approach is mainly used for the bands with much larger changes in their time-dependent RVS, such as the Terra MODIS bands 1-4, 8, and 9 and the Aqua MODIS bands 8- and 9. In this paper, the algorithms of these approaches are described, their performance is demonstrated, and their impact on L1B products is discussed. In general, the shorter wavelength bands have experienced a larger on-orbit RVS change, which, in general, are mirror side and detector dependent. The on-orbit RVS change due to the degradation of band 8 can be as large as 35 percent for Terra MODIS and 20 percent for Aqua MODIS. Vital to maintaining the accuracy of the MODIS L1B products is an accurate characterization of the on-orbit RVS change. The derived time-independent RVS, implemented in C6, makes an important improvement to the quality of the MODIS L1B products.

MODIS↗

Shallow Lunar Seismic Activity and the Current Stress State of the Moon

A vast, global network of more than 3200 lobate thrust fault scarps has been revealed in high resolution Lunar Reconnaissance Orbiter Camera (LROC) images. The fault scarps are very young, less than 50 Ma, based on their small scale and crisp appearance, crosscutting relations with small-diameter impact craters, and rates of infilling of associated small, shallow graben and may be actively forming today. The population of young thrust fault scarps provides a window into the recent stress state of the Moon and offers insight into the origin of global lunar stresses. The distribution of orientations of the fault scarps is non-random, inconsistent with isotropic stresses from late-stage global contraction as the sole source of stress. Modeling shows that tidal stresses contribute significantly to the current stress state of the lunar crust. Tidal stresses (orbital recession and diurnal tides) superimposed on stresses from global contraction result in non-isotropic compressional stress and may produce thrust faults consistent with lobate scarp orientations. At any particular point on the lunar surface, peak compressive stress will be reached at a certain time in the diurnal cycle. Coseismic slip events on currently active thrust faults are expected to be triggered when peak stresses are reached. Analysis of the timing of the 28 the shallow moonquakes recorded by the Apollo seismic network shows that 19 indeed occur when the Moon is closer to apogee, while only 9 shallow events occur when the Moon is closer to perigee. Here we report efforts to refine the model for the current stress state of the Moon by investigating the contribution of polar wander. Progress on relocating the epicentral locations of the shallow moonquakes using an algorithm designed for sparse networks is also reported.

structure↗

A Study of the Non-Thermal X-ray Emission of Shell-Type Supernova Remnants

The term of the forth year of the award is the period from March 15, 2003 to March 14, 2004. During this year, Dr. Thomas Pannuti, who had been performing most of the analyses, made a transition to a new position at Caltech. As of September 2003, Dr. Michael Stage began performing most of the analyses. Dr. Stage has begun constructing a detailed catalog of the spatial and spectral properties of young supernova remnants as described in sections 1 and 4 of the proposal. Specifically, he has focused on the analysis of Chandra ACIS data. The exquisite spatial resolution of the Chandra telescope and the modest spectral resolution of the CCDs make these data ideal. Dr. Stage has developed a standard set of procedures to reduce the data, generate telescope and detector response libraries and handle instrumental and celestial background subtraction. He has also adapted some existing code to automate spectral extraction and fitting. This code enables us to analyze the spectra of tens of thousands of small subregions of supernova remnants. Dr. Stage has recently applied this process to ACIS data for the supernova remnants Cas A and Kepler. The results include maps of the fitted parameters, such as individual line intensities and centroids (i.e. Doppler shifts), the electron temperature, and the absorption column density. These maps are more accurate than simple energy-cut images because it is difficult (and sometimes impossible) to cleanly separate the line emission from the underlying continuum (especially at low energies). It is now possible to identify the locations at which emission from each element is produced. Some of the elements are seen to be layered. The electron temperature maps demonstrate that the outer edges of Cas A and Kepler are very hot with weak lines. This emission is almost certainly synchrotron dominated. Therefore, not only are the maps interesting in and of themselves, but they also provide an unbiased means of easily identifying features in the remnants that can be studied in more detail using, for example, synchrotron models. Over the next year, Dr. Stage will refine the analysis procedures, expand the analyses to include several other remnants, present the results at major scientific conferences and publish a catalog of the results. Since this work is producing impressive results, we have applied for very long Chandra observations of Kepler and Tycho to obtain enough counts to take lull advantage of the technique. A i Ms observation of Cas A is already scheduled. Over the past year, I completed a joint, spectra,l analysis of some X-ray, radio, and gamma-ray data for the supernova remnant SN 1006. The results of this analysis show that the synchrotron and, hence, electron spectrum of the remnant is curved. The amount of curvature in the electron spectrum is quantitatively consistent with predictions of the amount of curvature in the proton spectrum of the remnant. A paper describing this work is nearly complete. When Dr. Pannuti left, he was analyzing X-ray data for the young supernova remnant G266.2-1.2. Dr. Pannuti found the first evidence of thermal X-ray emission from this remnant. Furthermore, like SN 1006, the cutoff frequency varies with azimuth along the bright northeastern rim. This work will be published during the coming year.

Allen, Glenn E.↗

On the Universality of the Kolmogorov Constant in Numerical Simulations of Turbulence

Motivated by a recent survey of experimental data, we examine data on the Kolmogorov spectrum constant in numerical simulations of isotropic turbulence, using results both from previous studies and from new direct numerical simulations over a range of Reynolds numbers (up to 240 on the Taylor scale) at grid resolutions up to 512(exp 3). It is noted that in addition to k(exp -5/3) scaling, identification of a true inertial range requires spectral isotropy in the same wavenumber range. We found that a plateau in the compensated three-dimensional energy spectrum at k(eta) approx. = 0.1 - -0.2, commonly used to infer the Kolmogorov constant from the compensated three-dimensional energy spectrum, actually does not represent proper inertial range behavior. Rather, a proper, if still approximate, inertial range emerges at k(eta) approx. = 0.02 - 0.05 when R(sub lambda) increases beyond 140. The new simulations indicate proportionality constants C(sub 1) and C in the one- and three-dimensional energy spectra respectively about 0.60 and 1.62. If the turbulence were perfectly isotropic then use of isotropy relations in wavenumber space (C(sub 1) = 18/55 C) would imply that C(sub 1) approx. = 0.53 for C = 1.62, in excellent agreement with experiments. However the one- and three-dimensional estimates are not fully consistent, because of departures (due to numerical and statistical limitations) from isotropy of the computed spectra at low wavenumbers. The inertial scaling of structure functions in physical space is briefly addressed. Since DNS is still restricted to moderate Reynolds numbers, an accurate evaluation of the Kolmogorov constant is very difficult. We focus on providing new insights on the interpretation of Kolmogorov 1941 similarity in the DNS literature and do not consider issues pertaining to the refined similarity hypotheses of Kolmogorov (K62).

Yeung, P. K.↗

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗

Polar Dunes Resolved by the Mars Orbiter Laser Altimeter Gridded Topography and Pulse Widths

The Mars Orbiter Laser Altimeter (MOLA) polar data have been refined to the extent that many features poorly imaged by Viking Orbiters are now resolved in densely gridded altimetry. Individual linear polar dunes with spacings of 0.5 km or more can be seen as well as sparsely distributed and partially mantled dunes. The refined altimetry will enable measurements of the extent and possibly volume of the north polar ergs. MOLA pulse widths have been recalibrated using inflight data, and a robust algorithm applied to solve for the surface optical impulse response. It shows the surface root-mean-square (RMS) roughness at the 75-m-diameter MOLA footprint scale, together with a geological map. While the roughness is of vital interest for landing site safety studies, a variety of geomorphological studies may also be performed. Pulse widths corrected for regional slope clearly delineate the extent of the polar dunes. The MOLA PEDR profile data have now been re-released in their entirety (Version L). The final Mission Experiment Gridded Data Records (MEGDR's) are now provided at up to 128 pixels per degree globally. Densities as high as 512 pixels per degree are available in a polar stereographic projection. A large computational effort has been expended in improving the accuracy of the MOLA altimetry themselves, both in improved orbital modeling and in after-the-fact adjustment of tracks to improve their registration at crossovers. The current release adopts the IAU2000 rotation model and cartographic frame recommended by the Mars Cartography Working Group. Adoption of the current standard will allow registration of images and profiles globally with an uncertainty of less than 100 m. The MOLA detector is still operational and is currently collecting radiometric data at 1064 nm. Seasonal images of the reflectivity of the polar caps can be generated with a resolution of about 300 m per pixel.

Gregory A Neumann↗