Search NASA⌕ Search

SEARCH · Search NASA

Results for “Error Rate Predictions”

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 433 records · Page 24

Mapping the Gas Turbulence in the Coma Cluster: Predictions for Astro-H

Astro-H will be able for the first time to map gas velocities and detect turbulence in galaxy clusters. One of the best targets for turbulence studies is the Coma cluster, due to its proximity, absence of a cool core, and lack of a central active galactic nucleus. To determine what constraints Astro-H will be able to place on the Coma velocity field, we construct simulated maps of the projected gas velocity and compute the second-order structure function, an analog of the velocity power spectrum. We vary the injection scale, dissipation scale, slope, and normalization of the turbulent power spectrum, and apply measurement errors and finite sampling to the velocity field. We find that even with sparse coverage of the cluster, Astro-H will be able to measure the Mach number and the injection scale of the turbulent power spectrum-the quantities determining the energy flux down the turbulent cascade and the diffusion rate for everything that is advected by the gas (metals, cosmic rays, etc.). Astro-H will not be sensitive to the dissipation scale or the slope of the power spectrum in its inertial range, unless they are outside physically motivated intervals. We give the expected confidence intervals for the injection scale and the normalization of the power spectrum for a number of possible pointing configurations, combining the structure function and velocity dispersion data. Importantly, we also determine that measurement errors on the line shift will bias the velocity structure function upward, and show how to correct this bias.

clusters: intracluster medium↗

Error Sources and Mitigation Strategies for Thermocouples Integrated in Flexible Thermal Protection System Materials

• The flexible thermal protection system (FTPS) on NASA’s Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) vehicle will be instrumented with thermocouples (TCs) to measure the in-depth thermal response during entry into Earth’s atmosphere. Accurate flight temperature measurements are critical for verifying vehicle performance and reducing thermal model uncertainties. • The TC leads are routed from the measurement location to the data acquisition system within the same FTPS layer that they are monitoring the temperature. This approach eliminates the need to put holes in the FTPS layers. However, the insulated TC leads are exposed to high temperatures and large thermal gradients. • The baseline TCs were commercially available 30 AWG Type K TCs with a binder impregnated glass braid insulation. The glass braid is rated to a maximum continuous use temperature of 482°C. LOFTID’s heat pulse will be on the order of minutes and the maximum predicted temperature beneath the outermost FTPS layers is 1350°C.

R. A. Miller↗

Challenges and Solutions for Precision Solar Pointing on the ISS for the TSIS Instrument

The International Space Station offers an attractive platform for external payloads in many regards, but its large, complex structure is challenging for payloads that require precise, unobstructed pointing. As an instrument that measures solar irradiance on a continual basis, the TSIS payload was designed to robustly account for these ISS challenges and is now performing its mission on the ISS with great success.This paper first describes how TSIS was able to measure the ISS background jitter using High-rate Fine Sun Sensor data that was taken during passive solar transits through its field of view. The measurements showed that ISS jitter at the TSIS location on ELC-3 Site 5 is relatively benign with less than 4 arcseconds 1 at frequencies less than 0.5 Hz. This value is lower than previous estimates but agrees well with recent measurements from other external payloads, including OPALS and SAGE III. Because TSIS was designed to be robust to base motion jitter, it has achieved less than 4 arcseconds 1 of pointing error, which provides an extremely precise pointing platform for the solar irradiance measurements. The other challenging ISS characteristic is its large, complex structure that obscures solar viewing at sunrise and sunset from the zenith-facing location of TSIS. Although TSIS was able to predict solar viewing durations during the design phase with relatively high accuracy, the exact timing and causes of obscuration were difficult to predict and understand after TSIS began operating on the ISS. To aid in the understanding of blockage by ISS structure, detailed obscuration maps were created using HFSS measurements that clearly show the exact gimbal angles where blockage occurs and the specific structural features of the ISS that cause the blockage. These maps have helped TSIS to maximize its solar tracking durations each orbit, which has supported the mission in meeting its science objectives.

background jitter↗

A New Eddy Dissipation Rate Formulation for the Terminal Area PBL Prediction System(TAPPS)

The TAPPS employs the MASS model to produce mesoscale atmospheric simulations in support of the Wake Vortex project at Dallas Fort-Worth International Airport (DFW). A post-processing scheme uses the simulated three-dimensional atmospheric characteristics in the planetary boundary layer (PBL) to calculate the turbulence quantities most important to the dissipation of vortices: turbulent kinetic energy and eddy dissipation rate. TAPPS will ultimately be employed to enhance terminal area productivity by providing weather forecasts for the Aircraft Vortex Spacing System (AVOSS). The post-processing scheme utilizes experimental data and similarity theory to determine the turbulence quantities from the simulated horizontal wind field and stability characteristics of the atmosphere. Characteristic PBL quantities important to these calculations are determined based on formulations from the Blackadar PBL parameterization, which is regularly employed in the MASS model to account for PBL processes in mesoscale simulations. The TAPPS forecasts are verified against high-resolution observations of the horizontal winds at DFW. Statistical assessments of the error in the wind forecasts suggest that TAPPS captures the essential features of the horizontal winds with considerable skill. Additionally, the turbulence quantities produced by the post-processor are shown to compare favorably with corresponding tower observations.

Charney, Joseph J.↗

Characterization of Errors in Satellite-Based HCHO/NO 2 Tropospheric Column Ratios With Respect to Chemistry, Column-to-PBL Translation, Spatial Representation, and Retrieval Uncertainties

The availability of formaldehyde (HCHO) (a proxy for volatile organic compound reactivity) and nitrogen dioxide (NO 2 ) (a proxy for nitrogen oxides) tropospheric columns from ultraviolet–visible (UV–Vis) satellites has motivated many to use their ratios to gain some insights into the near-surface ozone sensitivity. Strong emphasis has been placed on the challenges that come with transforming what is being observed in the tropospheric column to what is actually in the planetary boundary layer (PBL) and near the surface; however, little attention has been paid to other sources of error such as chemistry, spatial representation, and retrieval uncertainties. Here we leverage a wide spectrum of tools and data to quantify those errors carefully. Concerning the chemistry error, a well-characterized box model constrained by more than 500 h of aircraft data from NASA's air quality campaigns is used to simulate the ratio of the chemical loss of HO 2 + RO 2 (LRO x ) to the chemical loss of NO x (LNO x ). Subsequently, we challenge the predictive power of HCHO/NO 2 ratios (FNRs), which are commonly applied in current research, in detecting the underlying ozone regimes by comparing them to LRO x /LNO x . FNRs show a strongly linear (R 2 =0.94) relationship with LRO x /LNO x , but only on the logarithmic scale. Following the baseline (i.e., ln(LRO x /LNO x ) = −1.0 ± 0.2) with the model and mechanism (CB06, r2) used for segregating NO x -sensitive from VOC-sensitive regimes, we observe a broad range of FNR thresholds ranging from 1 to 4. The transitioning ratios strictly follow a Gaussian distribution with a mean and standard deviation of 1.8 and 0.4, respectively. This implies that the FNR has an inherent 20 % standard error (1σ) resulting from not accurately describing the RO x –HO x cycle. We calculate high ozone production rates (PO 3 ) dominated by large HCHO × NO 2 concentration levels, a new proxy for the abundance of ozone precursors. The relationship between PO 3 and HCHO × NO 2 becomes more pronounced when moving towards NO x -sensitive regions due to nonlinear chemistry; our results indicate that there is fruitful information in the HCHO × NO 2 metric that has not been utilized in ozone studies. The vast amount of vertical information on HCHO and NO 2 concentrations from the air quality campaigns enables us to parameterize the vertical shapes of FNRs using a second-order rational function permitting an analytical solution for an altitude adjustment factor to partition the tropospheric columns into the PBL region. We propose a mathematical solution to the spatial representation error based on modeling isotropic semivariograms. Based on summertime-averaged data, the Ozone Monitoring Instrument (OMI) loses 12 % of its spatial information at its native resolution with respect to a high-resolution sensor like the TROPOspheric Monitoring Instrument (TROPOMI) (> 5.5 × 3.5 km 2 ). A pixel with a grid size of 216 km 2 fails at capturing ∼ 65 % of the spatial information in FNRs at a 50 km length scale comparable to the size of a large urban center (e.g., Los Angeles). We ultimately leverage a large suite of in situ and ground-based remote sensing measurements to draw the error distributions of daily TROPOMI and OMI tropospheric NO 2 and HCHO columns. At a 68 % confidence interval (1σ), errors pertaining to daily TROPOMI observations, either HCHO or tropospheric NO 2 columns, should be above 1.2–1.5 × 10 16 molec. cm −2 to attain a 20 %–30 % standard error in the ratio. This level of error is almost non-achievable with the OMI given its large error in HCHO. The satellite column retrieval error is the largest contributor to the total error (40 %–90 %) in the FNRs. Due to a stronger signal in cities, the total relative error (< 50 %) tends to be mild, whereas areas with low vegetation and anthropogenic sources (e.g., the Rocky Mountains) are markedly uncertain (> 100 %). Our study suggests that continuing development in the retrieval algorithm and sensor design and calibration is essential to be able to advance the application of FNRs beyond a qualitative metric.

Amir H. Souri↗

TEM Analyses of Itokawa Regolith Grains and Lunar Soil Grains to Directly Determine Space Weathering Rates on Airless Bodies

Samples returned from the moon and Asteroid Itokawa by NASA's Apollo Missions and JAXA's Hayabusa Mission, respectively, provide a unique record of their interaction with the space environment. Space weathering effects result from micrometeorite impact activity and interactions with the solar wind. While the effects of solar wind interactions, ion implantation and solar flare particle track accumulation, have been studied extensively, the rate at which these effects accumulate in samples on airless bodies has not been conclusively determined. Results of numerical modeling and experimental simulations do not converge with observations from natural samples. We measured track densities and rim thicknesses of three olivine grains from Itokawa and multiple olivine and anorthite grains from lunar soils of varying exposure ages. Samples were prepared for analysis using a Leica EM UC6 ultramicrotome and an FEI Quanta 3D dual beam focused ion beam scanning electron microscope (FIB-SEM). Transmission electron microscope (TEM) analyses were performed on the JEOL 2500SE 200kV field emission STEM. The solar wind damaged rims on lunar anorthite grains are amorphous, lack inclusions, and are compositionally similar to the host grain. The rim width increases as a smooth function of exposure age until it levels off at approximately 180 nm after approximately 20 My (Fig. 1). While solar wind ion damage can only accumulate while the grain is in a direct line of sight to the Sun, solar flare particles can penetrate to mm-depths. To assess whether the track density accurately predicts surface exposure, we measured the rim width and track density in olivine and anorthite from the surface of rock 64455, which was never buried and has a surface exposure age of 2 My based on isotopic measurements. The rim width from 64455 (60-70nm) plots within error of the well-defined trend for solar wind amorphized rims in Fig. 1. Measured solar flare track densities are accurately reflecting the surface exposure of the grains. Track densities correlate with the amorphous rim thicknesses. While the space-weathered rims of anorthite grains are amorphous, the space-weathered rims on both Itokawa and lunar olivine grains show solar wind damaged rims that are not amorphous. Instead, the rims are nanocrystalline with high dislocation densities and sparse inclusions of nanophase Fe metal. The rim thicknesses on the olivine grains also correlate with track density. The Itokawa olivine grains have track densities that indicate surface exposures of approximately 10(exp 5) years. Longer exposures (up to approximately 10(exp 7) years) do not amorphize the rims, as evidenced by lunar soil olivines with high track densities (approximately 10(exp 11) cm(exp -2)). From the combined data, shown in Fig. 1, it is clear that olivine is damaged (but not amorphized) more rapidly by the solar wind compared to anorthite. The olivine damaged rim forms quickly (in approximately 10(exp 6) y) and saturates at approximately 120nm with longer exposure time. The anorthite damaged rims form more slowly, amorphize, and grow thicker than the olivine rims. This is in agreement with numerical modeling data which predicts that solar wind damaged rims on anorthite will be thicker than olivine. However, the models predict that both olivine and anorthite rims will amorphize and reach equilibrium widths in less than 10(exp 3) y, in contrast to what is observed for natural samples. Laboratory irradiation experiments, which show rapid formation of fully amorphous and blistered surfaces from simulated solar wind exposures are also in contrast to observations of natural samples. These results suggest that there is a flux dependence on the type and extent of irradiation damage that develops in olivine. This flux dependence suggests that great caution be used in extrapolating between high-flux laboratory experiments and the natural case, as demonstrated by. We constrain the space weathering rate through analysis of returned samples. Provided that the track densities and the solar wind damaged rim widths exhibited by the Itokawa grains are typical of the fine-grained regions of Itokawa, then the space weathering rate is on the order of 10(exp 5) y. Space weathering effects in lunar soils saturate within a few My of exposure while those in Itokawa regolith grains formed in approximately 10(exp 5) y. Olivine and anorthite respond differently to solar wind irradiation. The space weathering effects in olivine are particularly difficult to reconcile with laboratory irradiation studies and numerical models. Additional measurements, experiments, and modeling are required to resolve the discrepancies among the observations and calculations involving solar wind amorphization of different minerals on airless bodies.

Berger, Eve L.↗

Artificial intelligence time series forecasting for feed-forward laser stabilization

Laser plasma accelerators, typically operating at 1–10 Hz repetition rates, have the ability to produce high-quality electron beams in compact, all-optical-driven configurations, with the electron beams uniquely suited for a wide variety of accelerator-based applications. However, fluctuations and drifts in the laser delivery to the meter-scaled and below plasma target (the electron beam source) will translate into electron beam source variations that can limit their utility for demanding applications like light sources or linear colliders. Commercially available active feedback laser stabilization systems are intrinsically bandwidth limited due to their integration with multi-inch corrective mirror mounts which minimizes their effectiveness. In this manuscript, we present a Neural Network time series forecaster that can predict laser position fluctuations of the laser delivery to the final target well ahead of a future laser shot. The Root-Mean-Square-Error (RMSE) of the prediction accuracy was < 2 μ m for a 1 / e 2 beam radius of 34 μ m . Our feed-forward approach serves as a first-step in circumventing the bandwidth limitations imposed by the currently available stabilization systems since it allows for mirrors to be moved into position ahead of time to offset the predicted future position drift. This will help advance laser plasma accelerator research by providing greater robustness and stability needed for its applications.

Berger, Curtis↗

Ozone Budgets from a Global Chemistry/ Transport Model and Comparison to Observations from POLARIS

The objective of the Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) field mission was to obtain data to better characterize the summertime seasonal decrease of ozone at mid to high latitudes. The decrease in ozone occurs mainly in the lower stratosphere and is expected to result from in situ chemical destruction. Instrumented balloons and aircraft were used in POLARIS, along with satellites, to measure ozone and chemical species which are involved with stratospheric ozone chemistry. In order to close the seasonal ozone budget, however, ozone transport must also be estimated. Comparison to a global chemistry and transport model (CTM) of the stratosphere indicates how well the summertime ozone loss processes are simulated and thus how well we can predict the ozone response to changing amounts of chemical source gases. Moreover, the model gives insight into the possible relative magnitude of transport contributions to the seasonal ozone decline. Initial comparison to the Goddard CTM, which uses transport winds and temperatures from meteorological data assimilation, shows a high ozone bias in the model and an attenuated summertime ozone loss cycle. Comparison of the model chemical partitioning, and ozone catalytic loss rates to those derived from measurements shows fairly close agreement both at ER-2 altitudes (20 km) and higher. This suggests that the model transport is too active in resupplying ozone to the high latitude region, although chemistry failings cannot be completely ruled out. Comparison of ozone and related species will be shown along with a full diagnosis of the model ozone budget and its possible sources of error.

Kawa, S. Randy↗

Ozone Budgets from a Global Chemistry/Transport Model and Comparison to Observations from POLARIS

The objective of the Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) field mission was to obtain data to better characterize the summertime seasonal decrease of ozone at mid to high latitudes. The decrease in ozone occurs mainly in the lower stratosphere and is expected to result from in situ chemical destruction. Instrumented balloons and aircraft were used in POLARIS, along with satellites, to measure ozone and chemical species which are involved with stratospheric ozone chemistry. In order to close the seasonal ozone budget, however, ozone transport must also be estimated. Comparison to a global chemistry and transport model (CTM) of the stratosphere indicates how well the summertime ozone loss processes are simulated and thus how well we can predict the ozone response to changing amounts of chemical source gases. Moreover, the model gives insight into the possible relative magnitude of transport contributions to the seasonal ozone decline. Initial comparison to the Goddard CTM, which uses transport winds and temperatures from meteorological data assimilation, shows a high ozone bias in the model and an attenuated summertime ozone loss cycle. Comparison of the model chemical partitioning and ozone catalytic loss rates to those derived from measurements shows fairly close agreement both at ER-2 altitudes (20 km) and higher. This suggests that the model transport is too active in resupplying ozone to the high latitude region, although chemistry failings cannot be completely ruled out. Comparison of ozone and related species will be shown along with a full diagnosis of the model ozone budget and its possible sources of error.

Kawa, S. Randolph↗

Improved Estimate of Phobos Secular Acceleration from MOLA Observations

We report on new observations of the orbital position of Phobos, and use them to obtain a new and improved estimate of the rate of secular acceleration in longitude due to tidal dissipation within Mars. Phobos is the inner-most natural satellite of Mars, and one of the few natural satellites in the solar system with orbital period shorter than the rotation period of its primary. As a result, any departure from a perfect elastic response by Mars in the tides raised on it by Phobos will cause a transfer of angular momentum from the orbit of Phobos to the spin of Mars. Since its discovery in 1877, Phobos has completed over 145,500 orbits, and has one of the best studied orbits in the solar system, with over 6000 earth-based astrometric observations, and over 300 spacecraft observations. As early as 1945, Sharpless noted that there is a secular acceleration in mean longitude, with rate (1.88 + 0.25) 10(exp -3) degrees per square year. In preparation for the 1989 Russian spacecraft mission to Phobos, considerable work was done compiling past observations, and refining the orbital model. All of the published estimates from that era are in good agreement. A typical solution (Jacobson et al., 1989) yields (1.249 + 0.018) 10(exp -3) degrees per square year. The MOLA instrument on MGS is a laser altimeter, and was designed to measure the topography of Mars. However, it has also been used to make observations of the position of Phobos. In 1998, a direct range measurement was made, which indicated that Phobos was slightly ahead of the predicted position. The MOLA detector views the surface of Mars in a narrow field of view, at 1064 nanometer wavelength, and can detect shadows cast by Phobos on the surface of Mars. We have found 15 such serendipitous shadow transit events over the interval from xx to xx, and all of them show Phobos to be ahead of schedule, and getting progressively farther ahead of the predicted position. In contrast, the cross-track positions are quite close to the predicted values. Assuming that the along-track discrepancy is small enough that we can linearize the corrections, we model the mean orbital longitude as a quadratic function of time, and solve for corrections to the constant, linear, and quadratic terms. The time span of the recent observations is insufficient to properly resolve this issue alone, but when the 127 years of prior observations are added, we find a solution which reduces misfit to the new observations considerably, and makes no significant change to the fit to earlier observations. Our estimate for the secular acceleration term is (1.367 + 0.006) degrees per square year. The corresponding rate of energy dissipation is 3.34 MW. From a geophysical perspective, a more interesting parameter than the secular acceleration itself is the tidal lag angle, or tidal quality factor Q, for Mars. Unfortunately, the limiting error source in that determination is remaining uncertainty in the tidal Love numbers at harmonic degrees 2,3, and even 4. Until those parameters are better constrained, improvement in the orbital model of Phobos will not provide corresponding benefits for understanding the interior of Mars.

Bills, Bruce↗

Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements.

Baurle, R. A.↗

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗

Reflectivity and other physicochemical properties of Mn-substituted goethites and hematites

Goethite-bearing samples with values of Mn(s) (Mn/(Mn+Fe) mole fraction) up to 0.206 were synthesized by precipitation from alkaline solution. Samples with Mn(s) less than or equal 0.061 were single-phase Mn-goethites: samples with higher Mn(s) values contained another Mn-bearing phase (probably jacobsite). Mn-hematites were prepared by dehydroxylation of corresponding Mn-goethites at 500 C. Orthorhombic a and b unit cell dimensions of Mn-goethites changed in a linear manner with Mn(s), but not at rates predicted by the Vegrad law. Hexagonal unit cell dimensions of Mn-hematites did not vary with Mn(s). Moessbauer parameters isomer shift (IS), quadrupole splitting (QS), and hyperfine field (B(sub hf)) were measured at 293 and 15 K. For all single-phase Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061), magnetic splitting was observed at both temperatures. At 293 K, small but systematic decreases in B(sub hf) were observed with increasing Mn substitution; IS and QS were not dependent on Mn(s). Mn substitution strongly lowered the Morin transition temperature of hematite. At 15 K, the Morin transition was not present for Mn(s) greater than 0.020(4). The saturation magnetization of Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061) was the expected zero (within error) for antiferromagnetic goethite and for hematites obtained from dehydroxylation of goethites. Mn-geothites with Mn(s) greater than 0.061 were magnetic because of the presence of strongly magnetic jacobsite. For reflectivity spectra, bands resulting from MN(3+) were centered near 454 and 596 nm for Mn-goethites and near 545 and 700 nm for Mn-hematites. There is evidence for a approximately 700 nm band in spectral data for Martian bright regions, but association of it with Mn(3+) is not a unique interpretation. Comparison of laboratory and Martian spectral data implies that Mn(s) less than 0.032 for the Mn(3+) content of Martian hematites.

Vempati, R. K.↗

Dispersion and Lifetime of the SO2 Cloud from the August 2008 Kasatochi Eruption

Hemispherical dispersion of the SO2 cloud from the August 2008 Kasatochi eruption is analyzed using satellite data from the Ozone Monitoring Instrument (OMI) and the Goddard Trajectory Model (GTM). The operational OMI retrievals underestimate the total SO2 mass by 20-30% on 8-11 August, as compared with more accurate offline Extended Iterative Spectral Fit (EISF) retrievals, but the error decreases with time due to plume dispersion and a drop in peak SO2 column densities. The GTM runs were initialized with and compared to the operational OMI SO2 data during early plume dispersion to constrain SO2 plume heights and eruption times. The most probable SO2 heights during initial dispersion are estimated to be 10-12 km, in agreement with direct height retrievals using EISF algorithm and IR measurements. Using these height constraints a forward GTM run was initialized on 11 August to compare with the month-long Kasatochi SO2 cloud dispersion patterns. Predicted volcanic cloud locations generally agree with OMI observations, although some discrepancies were observed. Operational OMI SO2 burdens were refined using GTM-predicted mass-weighted probability density height distributions. The total refined SO2 mass was integrated over the Northern Hemisphere to place empirical constraints on the SO2 chemical decay rate. The resulting lower limit of the Kasatochi SO2 e-folding time is approx.8-9 days. Extrapolation of the exponential decay back in time yields an initial erupted SO2 mass of approx.2.2 Tg on 8 August, twice as much as the measured mass on that day.

Krotkov, N. A.↗

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements.

Williams, Christopher A.↗

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗

GEMS

The prime scientific objectives of the Gravity and Extreme Magnetism SMEX, GEMS, are to determine the effects of the spin of black holes, the configurations of the magnetic fields of magnetars, and the structure of the supernova shocks which accelerate cosmic rays. In the cases of both stellar black holes and supermassive black holes, sensitivity to 1% polarization is needed to make diagnostic measurements of the net polarizations predicted for probable disk and corona models. GEMS can reach this goal for several Seyferts and quasars and measure the polarizations of representatives of a variety of other classes of X-ray sources, such as rotation-powered and accretion-powered pulsars. GEMS uses foil mirrors to maximize the collecting area achievable within the SMEX constraints. The polarimeters at the mirror foci are Time Projection Chambers which use the photoelectic effect to measure the polarization of the incident photon. We have built laboratory models with good efficiency and modulation in the 2-10 keV range. An attached small student experiment would add 0.5 keV sensitivity for bright soft sources. The instrument has a point spread function which allows measurement of structures in the brighter nearby supernova remnants. GEMS' Orbital Sciences spacecraft will rotate at a rate of 0.1 revolutions per minute during observations, so that systematic errors due to the detector can be detected and corrected. A program of 35 sources can be observed in 9 months. GEMS is designed for a two year lifetime which will allow a General Observer program that would more than double the number of sources measured. For subsets of black holes, neutron stars and supernova remnants, GEMS will measure the polarization of several sources, solving important questions while establishing the sensitivity required for future missions.

Swank, Jean↗

The GREEN ‘omics of Nutrient Feedbacks to Soil Warming

The GREEN ‘omics of Nutrient Feedbacks in Soil project advanced the DOE Biological and Environmental Research (BER) mission by developing and applying isotope-enabled ’omics tools to understand how soil microbes regulate carbon and nutrient cycling. Guided by the Growth Rate, growth Efficiency, and stoichiometry of Essential Nutrients (GREEN ’omics) framework, the project aimed to build a predictive, systems-level understanding of microbial traits that control ecosystem biogeochemistry. In a collaboration among Northern Arizona University (lead), West Virginia University, Lawrence Livermore National Laboratory, and Pacific Northwest National Laboratory, we combined quantitative stable isotope probing (qSIP), Chip-SIP, NanoSIMS, and genome-resolved metagenomics across long-term experiments in Arctic, boreal, temperate, and tropical ecosystems. The project produced three key outcomes: 1) We showed that community-weighted temperature sensitivities of bacterial growth (Q10) can predict ecosystem-scale soil respiration responses across diverse soils. 2) We provided the first in situ evidence for density-dependent population dynamics in soil bacteria and demonstrated that nutrient additions intensify competition, concentrating carbon use into fewer taxa. 3) We improved and extended isotope-enabled ’omics methods by quantifying qSIP measurement error to guide experimental design and coupling SIP with genome-resolved metagenomics to reveal cross-kingdom interactions among bacteria, fungi, and viruses. Together, these results show that a small number of microbial traits and taxa exert disproportionate control over soil carbon and nutrient cycling, providing critical data and methods to improve representation of microbial processes in Earth system models.

54 ENVIRONMENTAL SCIENCES↗