Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Zero Correction”

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 613 records · Page 34

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium E h values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9.3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-1 and Table ES-2, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measurement of the Dispersion in the CeC Beamlines

The coherent electron cooling requires to match the relativistic factors of the hadron and electron beam within 10 -4 accuracy. Initial matching is based on the energy measurement of the electron beam using solenoid with accuracy about 10 -3 (the hadron relativistic factor can be found from the revolution frequency). Further matching is done with observation of the recombination rate with fine tuning based on the cooling rate. The last two steps require change of the electron energy which is achieved by change of the linac voltage. The interaction of the hadron and electron beams is affected by the transverse overlap of two beams. To minimize the effect, we need zero dispersion in the common section. The merge of the beams a dogleg transfer line is used containing two dipoles and three quadrupoles. The zeroing of dispersion is performed by settings of the quadrupoles. We have developed two ways to measure the dispersion. In the first mostly commonly used we change both dipoles by the same value creating the mismatch between the electron energy and bending magnet. In the compensated dogleg the position in the second dipole should not change, and the trajectory angle will change by the same amount as deflecting angle but in opposite direction. Therefore, the trajectory in the common section will not change. The value of the dispersion can be calculated from the beam displacement divided by the relative chang of the dipoles’ current. The second method is based on the change of the electron beam energy using linac. This method requires electron orbit alignment with linac axis to avoid beam steering. The results are shown in Fig.1 and Fig. 2. The first figure shows the dispersion prior the correction, and the second figure after the correction.

43 PARTICLE ACCELERATORS↗

Trends in the Vertical Distribution of Ozone: A Comparison of Two Analyses of Ozonesonde Data

We present the results of two independent analyses of ozonesonde measurements of the vertical profile of ozone. For most of the ozonesonde stations we use data that were recently reprocessed and reevaluated to improve their quality and internal consistency. The two analyses give similar results for trends in ozone. We attribute differences in results primarily to differences in data selection criteria and in utilization of data correction factors, rather than in statistical trend models. We find significant decreases in stratospheric ozone at all stations in middle and high latitudes of the northern hemisphere from 1970 to 1996, with the largest decreases located between 12 and 21 km, and trends of -3 to -10 %/decade near 17 km. The decreases are largest at the Canadian and the most northerly Japanese station, and are smallest at the European stations, and at Wallops Island, U.S.A. The mean mid-latitude trend is largest, -7 %/decade, from 12 to 17.5 km for 1970-96. For 1980-96, the decrease is more negative by 1-2 %/decade, with a maximum trend of -9 %/decade in the lowermost stratosphere. The trends vary seasonally from about 12 to 17.5 km, with largest ozone decreases in winter and spring. Trends in tropospheric ozone are highly variable and depend on region. There are decreases or zero trends at the Canadian stations for 1970-96, and decreases of -2 to -8 %/decade for the mid-troposphere for 1980-96; the three European stations show increases for 1970-96, but trends are close to zero for two stations for 1980-96 and positive for one; there are increases in ozone for the three Japanese stations for 1970-96, but trends are either positive or zero for 1980-96; the U.S. stations show zero or slightly negative trends in tropospheric ozone after 1980. It is not possible to define reliably a mean tropospheric ozone trend for northern mid-latitudes, given the small number of stations and the large variability in trends. The integrated column trends derived from the sonde data are consistent with trends derived from both surface based and satellite measurements of the ozone column.

Loogan, J. A.↗

Cosmological dynamics of string theory axion strings

The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu-Goldstone boson of global U⁢(1) PQ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we show that string theory axions behave fundamentally differently from field theory axions in the early Universe. Field theory axions may form axion strings if the PQ phase transition takes place after inflation. In contrast, we show that string theory axions do not generically form axion strings. In special inflationary paradigms, such as D-brane inflation, string theory axion strings may form; however, their tension is parametrically larger than that of field theory axion strings. We then show that such QCD axion strings overproduce the DM abundance for all allowed QCD axion masses and are thus ruled out, except in scenarios with large warping. A loop-hole to this conclusion arises in the axiverse, where an axion string could be composed of multiple different axion mass eigenstates; a heavier eigenstate could collapse the network earlier, allowing for the QCD axion to produce the correct DM abundance and also generating observable gravitational wave signals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electronic Structure of InAs and InSb Surfaces: Density Functional Theory and Angle-Resolved Photoemission Spectroscopy

The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, in this work, the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces is studied using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). Large-scale first principles simulations are enabled by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, a “bulk unfolding” scheme is implemented by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, a good agreement is found between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states, which may be detected by ARPES with low photon energies. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. The combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, for example, topological qubits utilizing the Majorana zero modes.

42 ENGINEERING↗

Solar-Cycle Changes in GONG P-Mode Widths and Amplitudes 1995-1998

We search for a solar cycle variation in mode widths and amplitudes derived from 3-month GONG time series. The variation of mode width and amplitude observed in GONG data are the combined effects of fill factor, temporal variation, and measurement uncertainties. The largest variation is caused by the fill factor resulting in modes with increased width and reduced amplitude when fill is lower. We assume that the solar cycle variation is the only other systematic variation beside the temporal window function effect. We correct all currently available data sets for the fill factor and simultaneously derive the solar cycle variation. We find an increase of about 3% on average in mode width from the previous minimum to Oct. 1998 and a decrease of about 7% and 6% in mode amplitude and mode area (width x amplitude). We find no l dependence of the solar-cycle changes. As a function of frequency, these changes show a maximum between 2.7 and 3.3 mHz with about 47% higher than average values for mode width and about 29% and 36% higher ones for mode amplitude and area. We estimate the significance of these rather small changes by a pre-whitening method and find that the results are significant at or above the 99.9% level with mode area showing the highest level of significance and mode width the lowest. The variation in background amplitude is most likely not significant and is consistent with a zero change.

Komm R. W.↗

Toward Trustworthy Autonomous Science: A Two-Year Community Roadmap

One year ago, the AISLE roadmap argued that autonomous laboratories operated as isolated islands and proposed a grassroots network organized around five critical dimensions. The field has since moved faster than that roadmap anticipated: multi-agent systems have produced experimentally validated hypotheses, self-driving laboratories have grown more interoperable and orchestrated, reasoning-trained and domain foundation models have raised the capability ceiling, and the Genesis Mission has placed autonomous experimentation at the center of U.S. federal science strategy, with industry emerging as a primary actor. Progress has met a sobering counter-current, including a corrected flagship discovery result, benchmarks showing that agents which rival experts on closed-ended questions still complete only a fraction of open-ended research, and fabricated citations surfacing at leading venues. We read this as the defining tension of the field: producing a candidate discovery is no longer the hard part, but verifying it is, and this asymmetry now limits autonomous science more than raw model capability. Accordingly, we update the roadmap around seven dimensions, revisiting the original five and elevating two former cross-cutting concerns, trust, verification, and reproducibility, and safety, security, and governance, to first-class status. We assess the original milestones (M1 through M14) as achieved, partially achieved, reframed, or open, add four new milestones (M15 through M18) for the elevated dimensions, and scope the path forward to a two-year horizon, with the first year concentrating on interfaces, protocol adoption, and the scaffolding of verification, and the second targeting federation, zero-trust coordination, and governance. Throughout, we position the grassroots network as the interoperability fabric that lets national programs, international initiatives, and commercial platforms connect rather than re-silo.

99 GENERAL AND MISCELLANEOUS↗

Rough Surfaces: Is the Dark Stuff Just Shadow?: "Who knows what evil lurks in the hearts of men? The shadow knows!"

Remote observations of the surfaces of airless planetary objects are fundamental to inferring the physical structure and compositional makeup of the surface material. A number of forward models have been developed to reproduce the photometric behavior of these surfaces, based on specific, assumed structural properties such as macroscopic roughness and associated shadowing. Most work of this type is applied to geometric albedos, which are affected by complicated effects near zero phase angle that represent only a tiny fraction of the net energy reflected by the object. Other applications include parameter fits to resolved portions of some planetary surface as viewed over a range of geometries. The spherical albedo of the entire object (when it can be determined) captures the net energy balance of the particle more robustly than the geometric albedo. In most treatments involving spherical albedos, spherical albedos and particle phase functions are often treated as if they are independent, neglecting the effects of roughness. In this paper we take a different approach. We note that whatever function captures the phase angle dependence of the brightness of a realistic rough, shadowed, flat surface element relative to that of a smooth granular surface of the same material, it is manifested directly in both the integral phase function and the spherical albedo of the object. We suggest that, where broad phase angle coverage is possible, spherical albedos may be easily corrected for the effects of shadowing using observed (or assumed) phase functions, and then modeled more robustly using smooth-surface regolith radiative transfer models without further imposed (forward-modeled) shadowing corrections. Our approach attributes observed "power law" phase functions of various slope (and "linear" ranges of magnitude-vs.-phase angle) to shadowing, as have others, and goes on to suggest that regolith-model-based inferences of composition based on shadow-uncorrected spherical albedos overestimate the amount of absorbing material contained in the regolith.

geometric albedos↗

Global Assessment of the SMAP Level-4 Soil Moisture Product Using Assimilation Diagnostics

The Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture and related land surface variables from 31 March 2015 to present with approx. 2.5-day latency. The ensemble-based L4_SM algorithm assimilates SMAP brightness temperature (Tb) observations into the Catchment land surface model. This study describes the spatially distributed L4_SM analysis and assesses the observation-minus-forecast (O-F) Tb residuals and the soil moisture and temperature analysis increments. Owing to the climatological rescaling of the Tb observations prior to assimilation, the analysis is essentially unbiased, with global mean values of approx. 0.37 K for the O-F Tb residuals and practically zero for the soil moisture and temperature increments. There are, however, modest regional (absolute) biases in the O-F residuals (under approx. 3 K), the soil moisture increments (under approx. 0.01 cu m/cu m), and the surface soil temperature increments (under approx. 1 K). Typical instantaneous values are approx. 6 K for O-F residuals, approx. 0.01 (approx. 0.003) cu m/cu m for surface (root-zone) soil moisture increments, and approx. 0.6 K for surface soil temperature increments. The O-F diagnostics indicate that the actual errors in the system are overestimated in deserts and densely vegetated regions and underestimated in agricultural regions and transition zones between dry and wet climates. The O-F auto-correlations suggest that the SMAP observations are used efficiently in western North America, the Sahel, and Australia, but not in many forested regions and the high northern latitudes. A case study in Australia demonstrates that assimilating SMAP observations successfully corrects short-term errors in the L4_SM rainfall forcing.

Reichle, Rolf↗

Using Data Assimilation Diagnostics to Assess the SMAP Level-4 Soil Moisture Product

The Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture and related land surface variables from 31 March 2015 to present with approx.2.5-day latency. The ensemble-based L4_SM algorithm assimilates SMAP brightness temperature (Tb) observations into the Catchment land surface model. This study describes the spatially distributed L4_SM analysis and assesses the observation-minus-forecast (O-F) Tb residuals and the soil moisture and temperature analysis increments. Owing to the climatological rescaling of the Tb observations prior to assimilation, the analysis is essentially unbiased, with global mean values of approx. 0.37 K for the O-F Tb residuals and practically zero for the soil moisture and temperature increments. There are, however, modest regional (absolute) biases in the O-F residuals (under approx. 3 K), the soil moisture increments (under approx. 0.01 cu m/cu m), and the surface soil temperature increments (under approx. 1 K). Typical instantaneous values are approx. 6 K for O-F residuals, approx. 0.01 (approx. 0.003) cu m/cu m for surface (root-zone) soil moisture increments, and approx. 0.6 K for surface soil temperature increments. The O-F diagnostics indicate that the actual errors in the system are overestimated in deserts and densely vegetated regions and underestimated in agricultural regions and transition zones between dry and wet climates. The O-F auto-correlations suggest that the SMAP observations are used efficiently in western North America, the Sahel, and Australia, but not in many forested regions and the high northern latitudes. A case study in Australia demonstrates that assimilating SMAP observations successfully corrects short-term errors in the L4_SM rainfall forcing.

SMAP↗

Lessons Learned in Space Life Support System Testing

The earlier problems can be found and corrected, the easier and cheaper it is to fix them. Doing less testing saves cost and time but doing too little testing increases the risk of operational failures causing large costs and delays. Integrated test is necessary to determine if the subsystems work together and the overall architecture performs as intended. This report reviews the testing lessons learned from the NASA Systems Engineering Handbook, a National Research Council report, and five reviews of International Space Station (ISS) lessons learned. The five reviews all mention two important points. First, that testing should be performed on the final integrated system, one as close as possible to the intended flight system. Second, “test as you fly,” while operating as planned in an environment as close as possible to the expected flight environment. Other lessons are the need for extensive preflight ground testing, the need to establish and defend an adequate budget, the problems using protoflight hardware on ISS, and the benefit of having ISS as a zero gravity test bed. The major ISS life support systems, carbon dioxide, water recycling, and oxygen recovery, were protoflight systems with little testing before launch to ISS. The failure rates these systems have been much greater than predicted and this has caused dissatisfaction with the protoflight approach. The more costly traditional approach is building qualification and test units in addition to flight units. The test units are used to test, analyze, and fix failure modes. Other work shows that there is an optimum cost-effective intuitive appeal of a human ecosystem in space.

Life support↗

Lessons Learned in Space Life Support System Testing

The earlier problems can be found and corrected, the easier and cheaper it is to fix them. Doing less testing saves cost and time but doing too little testing increases the risk of operational failures causing large costs and delays. Integrated test is necessary to determine if the subsystems work together and the overall architecture performs as intended. This report reviews the testing lessons learned from the NASA Systems Engineering Handbook, a National Research Council report, and five reviews of International Space Station (ISS) lessons learned. The five reviews all mention two important points. First, that testing should be performed on the final integrated system, one as close as possible to the intended flight system. Second, “test as you fly,” while operating as planned in an environment as close as possible to the expected flight environment. Other lessons are the need for extensive preflight ground testing, the need to establish and defend an adequate budget, the problems using protoflight hardware on ISS, and the benefit of having ISS as a zero gravity test bed. The major ISS life support systems, carbon dioxide removal, water recycling, and oxygen recovery, were protoflight systems with little testing before launch to ISS. The failure rates of these systems have been much greater than predicted and this has caused dissatisfaction with the protoflight approach. The more costly traditional approach builds qualification and test units in addition to flight units. The test units are used to find, analyze, and fix failure modes. Other work shows that there is an optimum cost-effective amount of testing when redundant systems must have a specified reliability and confidence.

Harry W. Jones↗

A Zero-Emission Process for Direct Reduction of Iron by Hydrogen Plasma in a Rotary Kiln Reactor

This project’s goal was to demonstrate a hydrogen plasma (H-plasma)-rotary kiln process for reducing iron ore to iron as part of the steel manufacturing process. The H-plasma provides a greater thermodynamic driving force for reducing iron ores than thermal processes such as the DRI process, enabling lower reaction temperatures. We estimated that our process technology can reduce energy consumption by 45% compared to the blast furnace process and ~15% compared to the DRI process. Steel manufacturing produces about 1.8 tons of CO2/ton of steel with iron ore reduction accounting for about one-third of the CO2 produced in the overall manufacturing process. We estimated our process has the potential to reduce GHG emissions from ironmaking by 35% with today’s grid and by up to 88% with a future low-carbon grid while being cost competitive with the current blast furnace route. We demonstrated reduction of hematite and magnetite rich materials at temperatures from 600 to 800°C. We achieved 90-95% metallization on 100 gr samples in batch reduction experiments in the H-plasma rotary kiln furnace at 600-650°C. Attempts to perform tests in a continuous operation mode identified problems with the ore feed mechanism. We identified solutions but there was not time nor budget to correct these for this project

36 MATERIALS SCIENCE↗

A Zero-Emission Process for Direct Reduction of Iron by Hydrogen Plasma in a Rotary Kiln Reactor

This project’s goal was to demonstrate a hydrogen plasma (H-plasma)-rotary kiln process for reducing iron ore to iron as part of the steel manufacturing process. The H-plasma provides a greater thermodynamic driving force for reducing iron ores than thermal processes such as the DRI process, enabling lower reaction temperatures. We estimated that our process technology can reduce energy consumption by 45% compared to the blast furnace process and ~15% compared to the DRI process. Steel manufacturing produces about 1.8 tons of CO2/ton of steel with iron ore reduction accounting for about one-third of the CO2 produced in the overall manufacturing process. We estimated our process has the potential to reduce GHG emissions from ironmaking by 35% with today’s grid and by up to 88% with a future low-carbon grid while being cost competitive with the current blast furnace route. We demonstrated reduction of hematite and magnetite rich materials at temperatures from 600 to 800°C. We achieved 90-95% metallization on 100 gr samples in batch reduction experiments in the H-plasma rotary kiln furnace at 600-650°C. Attempts to perform tests in a continuous operation mode identified problems with the ore feed mechanism. We identified solutions but there was not time nor budget to correct these for this project

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A field measure of the shade fraction

'Shade' has a technical definition peculiar to linear spectral mixture analysis of imaging spectrometer data: it is the reduction in radiance from a surface due to lighting conditions and geometry, and includes topographic shading described by photometric functions as well as shadowing at all scales. 'Shade' is an important constituent of nearly all remotely sensed images, and is one endmember resolved in spectral mixture analysis, where it is represented as a fraction of the measured radiance and a characteristic spectrum. This spectrum is typically the null vector, provided the data have been corrected for atmospheric and instrument effects: i.e., 'shade' is the radiance from an ideal black surface. In topographic shading, irradiance is reduced - typically in proportion to cos(i), where i (incidence angle) is the angle between the sun and the local surface normal vectors. Therefore, the radiance is lowered by a multiplicative factor. Shadowing occurs when i is greater than 90 deg, or when sunlight is blocked by adjacent high terrain; the only irradiance is down-welling skylight and bounce light from adjacent terrain. In spectral mixture analysis, 'shade' is regarded as an additive term. In this regard, it is an accurate description of the proportion of a scene that consists of ideal shadows ('checkerboard mixing'); however, 'shade' represents the multiplicative cos(i) factor as well, as here it should be interpreted as the proportion of shadow that would darken the scene an equivalent amount. In either case, the 'shade' fraction is lessened by adjacency effects, because the scene has a non-zero reflectivity instead of the ideal black surface generally assumed.

Gillespie, Alan R.↗

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗

Properties of the Hermean Regolith: New Optical Reflectance Spectra, Comparison with Lunar Anorthosites and Mineralogical Modelling - 5

We present new optical (0.4-0.65 micron) spectra of Mercury and lunar pure anorthosite locations, obtained quasi-simultaneously with the Nordic Optical Telescope (NOT) in 2002. A comparative study is performed with the model of Lucey et al. between iron-poor, mature, pure anorthosite (less than 90% plagioclase feldspar) Clementine spectra from the lunar farside and a combined 0.4-1.0 micron mercurian spectrum, obtained with the NOT, calculated for standard photometric geometry. Mercury is located at more extreme locations in the Lucey ratio-reflectance diagrams than any known lunar soil, specifically with respect to the extremely iron-poor mature anorthosites. Though quantitative prediction of FeO and TiO2 abundances cannot be made without a more generally applicable model, we find qualitatively that the abundances of both these oxides must be near zero for Mercury. We utilize the theory of Hapke, with realistic photometric parameters, to model laboratory spectra of matured mineral powders and lunar soils, and remotely sensed spectra of lunar anorthosites and Mercury. An important difference between fabricated and natural powders is the high value for the internal scattering parameter necessary to interpret the spectra for the former, and the requirement of rough and non-isotropically scattering surfaces in the modelling of the latter. The mature lunar anorthosite spectra were well modelled with binary mixtures of calcic feldspars and olivines, grain sizes of 25-30 micron and a concentration of submicroscopic metallic iron (SMFe) of 0.12-0.15% in grain coatings. The mercurian spectrum is not possible to interpret from terrestrial mineral powder spectra without introducing an average particle scattering function for the bulk soil that increases in backscattering efficiency with wavelength. The observed spectrum is somewhat better predicted with binary mixture models of feldspars and pyroxenes, that single-component regoliths consisting of either albite or diopside. Correct spectral reflectance values were predicted with a concentration of 0.1 wt% SMFe in coatings of 15-30 micron sized grains. Since reasonable cosmogonical formation scenarios for Mercury, or meteoritic infall, predict iron concentrations at least this high, we draw the conclusion that the average grain size of Mercury is about a factor of two smaller than for average returned lunar soil samples. The 0.6-2.5 micron spectrum of McCord and Clark is used to further limit the possible range of mineralogical composition of Mercury. It is found that an intimately mixed and matured 3 : 1 labradorite-to-enstatite regolith composition best matches both the optical and near-infrared spectra, yielding an abundance of approx. 1.2 wt% FeO and -approx. 0 wt% TiO2.

Warell, J.↗