Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time scale separation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72

Engineering Thermal Transport across Layered Graphene–MoS 2 Superlattices

Layering two-dimensional van der Waals materials provides a high degree of control over atomic placement, which could enable tailoring of vibrational spectra and heat flow at the sub-nanometer scale. In this work, using spatially resolved ultrafast thermoreflectance and spectroscopy, we uncover the design rules governing cross-plane heat transport in superlattices assembled from monolayers of graphene (G) and MoS 2 (M). Using a combinatorial experimental approach, we probe nine different stacking sequences, G, GG, MG, GGG, GMG, GGMG, GMGG, GMMG, and GMGMG, and identify the effects of vibrational mismatch, interlayer adhesion, and junction asymmetry on thermal transport. Pure G sequences display evidence of quasi-ballistic transport, whereas adding even a single M layer strongly disrupts heat conduction. The experimental data are described well by molecular dynamics simulations, which include thermal expansion, accounting for the effect of finite temperature on the interlayer spacing. The simulations show that an increase of ~2.4% in the layer separation of GMGMG, relative to its value at 300 K, can lead to a doubling of the thermal resistance. Using these design rules, we experimentally demonstrate a five-layer GMGMG superlattice “thermal metamaterial” with an ultralow effective cross-plane thermal conductivity comparable to that of air.

2D materials↗

Drop Size Distribution - Based Separation of Stratiform and Convective Rain

For applications in hydrology and meteorology, it is often desirable to separate regions of stratiform and convective rain from meteorological radar observations, both from ground-based polarimetric radars and from space-based dual frequency radars. In a previous study by Bringi et al. (2009), dual frequency profiler and dual polarization radar (C-POL) observations in Darwin, Australia, had shown that stratiform and convective rain could be separated in the log10(Nw) versus Do domain, where Do is the mean volume diameter and Nw is the scaling parameter which is proportional to the ratio of water content to the mass weighted mean diameter. Note, Nw and Do are two of the main drop size distribution (DSD) parameters. In a later study, Thurai et al (2010) confirmed that both the dual-frequency profiler based stratiform-convective rain separation and the C-POL radar based separation were consistent with each other. In this paper, we test this separation method using DSD measurements from a ground based 2D video disdrometer (2DVD), along with simultaneous observations from a collocated, vertically-pointing, X-band profiling radar (XPR). The measurements were made in Huntsville, Alabama. One-minute DSDs from 2DVD are used as input to an appropriate gamma fitting procedure to determine Nw and Do. The fitted parameters - after averaging over 3-minutes - are plotted against each other and compared with a predefined separation line. An index is used to determine how far the points lie from the separation line (as described in Thurai et al. 2010). Negative index values indicate stratiform rain and positive index indicate convective rain, and, moreover, points which lie somewhat close to the separation line are considered 'mixed' or 'transition' type precipitation. The XPR observations are used to evaluate/test the 2DVD data-based classification. A 'bright-band' detection algorithm was used to classify each vertical reflectivity profile as either stratiform or convective, depending on whether or not a clearly-defined melting layer is present at an expected height, and if present, maximum reflectivity within the melting layer as well as the corresponding height are determined. We will present results of quantitative comparisons between the XPR observations-based classifications and the simultaneous 2DVD data-based classifications. Time series comparisons will be presented for thirteen events in Huntsville.

Thurai, Merhala↗

Engineered Multifunctional Surfaces for Fluid Handling

Designs incorporating variations in capillary geometry and hydrophilic and/or antibacterial surface properties have been developed that are capable of passive gas/liquid separation and passive water flow. These designs can incorporate capillary grooves and/or surfaces arranged to create linear and circumferential capillary geometry at the micro and macro scale, radial fin configurations, micro holes and patterns, and combinations of the above. The antibacterial property of this design inhibits the growth of bacteria or the development of biofilm. The hydrophilic property reduces the water contact angle with a treated substrate such that water spreads into a thin layer atop the treated surface. These antibacterial and hydrophilic properties applied to a thermally conductive surface, combined with capillary geometry, create a novel heat exchanger capable of condensing water from a humid, two-phase water and gas flow onto the treated heat exchanger surfaces, and passively separating the condensed water from the gas flow in a reduced gravity application. The overall process to generate the antibacterial and hydrophilic properties includes multiple steps to generate the two different surface properties, and can be divided into two major steps. Step 1 uses a magnetron-based sputtering technique to implant the silver atoms into the base material. A layer of silver is built up on top of the base material. Completion of this step provides the antibacterial property. Step 2 uses a cold-plasma technique to generate the hydrophilic surface property on top of the silver layer generated in Step 1. Completion of this step provides the hydrophilic property in addition to the antibacterial property. Thermally conductive materials are fabricated and then treated to create the antibacterial and hydrophilic surface properties. The individual parts are assembled to create a condensing heat exchanger with antibacterial and hydrophilic surface properties and capillary geometry, which is capable of passive phase separation in a reduced gravity application. The plasma processes for creating antibacterial and hydrophilic surface properties are suitable for applications where water is present on an exposed surface for an extended time, such that bacteria or biofilms could form, and where there is a need to manage the water on the surface. The processes are also suitable for applications where only the hydrophilic property is needed. In particular, the processes are applicable to condensing heat exchangers (CHXs), which benefit from the antibacterial properties as well as the hydrophilic properties. Water condensing onto the control surfaces of the CHX will provide the moist conditions necessary for the growth of bacteria and the formation of biofilms. The antibacterial properties of the base layer (silver) will mitigate and prevent the growth of bacteria and formation of biofilms that would otherwise reduce the CHX performance. In addition, the hydrophilic properties reduce the water contact angle and prevent water droplets from bridging between control surfaces. Overall, the hydrophilic properties reduce the pressure drop across the CHX.

Thomas, Chris↗

Disturbed, diffuse, or just missing? A global study of the H I content of Hickson compact groups

Hickson compact groups (HCGs) are dense configurations of four to ten galaxies, whose H I morphology appears to follow an evolutionary sequence of three phases, with gas initially confined to galaxies, then significant amounts spread throughout the intra-group medium, and finally with almost no gas remaining in the galaxies themselves. It has also been suggested that several groups may harbour a diffuse H I component that is resolved out by interferometric observations. The H I deficiency of HCGs is expected to increase as the H I morphological phase progresses along the evolutionary sequence. If this is the case, H I deficiency would be a rough proxy for the age and evolutionary state of a HCG. We aim to test this hypothesis for the first time using a large sample of HCGs and to investigate the evidence for diffuse H I in HCGs. We performed a uniform reduction of all publicly available VLA H I observations (38 HCGs) with a purpose-built pipeline that also maximises the reproducibility of this study. The resulting H I data cubes were then analysed with the latest software tools to perform a manual separation of emission features into those belonging to galaxies and those extending into the intra-group medium. We thereby classified the H I morphological phase of each group as well as quantified their H I deficiency compared to galaxies in isolation. We find little evidence that H I deficiency can be used as a proxy for the evolutionary phase of a compact group in either of the first two phases, with the distribution of H I deficiency being consistent in both. However, for the final phase, the distribution clearly shifts to high H I deficiencies, with more than 90% of the expected H I content typically missing. Across all HCGs studied, we identify a few cases where there is strong evidence for a diffuse gas component in the intra-group medium, which might be detectable with improved observations. We also classify a new sub-phase where groups contain a lone H I -bearing galaxy, but are otherwise devoid of gas. The new morphological phase we have identified is likely the result of an evolved, gas-poor group acquiring a new, gas-rich member. The large spread of H I deficiencies in the first two morphological phases suggests that there is a broad range of initial H I content in HCGs, which is perhaps influenced by large-scale environment, and that the timescale for morphological changes is, in general, considerably shorter than the timescale for the destruction or consumption of neutral gas in these systems.

79 ASTRONOMY AND ASTROPHYSICS↗

A Transportable Gravity Gradiometer Based on Atom Interferometry

A transportable atom interferometer-based gravity gradiometer has been developed at JPL to carry out measurements of Earth's gravity field at ever finer spatial resolutions, and to facilitate high-resolution monitoring of temporal variations in the gravity field from ground- and flight-based platforms. Existing satellite-based gravity missions such as CHAMP and GRACE measure the gravity field via precise monitoring of the motion of the satellites; i.e. the satellites themselves function as test masses. JPL's quantum gravity gradiometer employs a quantum phase measurement technique, similar to that employed in atomic clocks, made possible by recent advances in laser cooling and manipulation of atoms. This measurement technique is based on atomwave interferometry, and individual laser-cooled atoms are used as drag-free test masses. The quantum gravity gradiometer employs two identical atom interferometers as precision accelerometers to measure the difference in gravitational acceleration between two points (Figure 1). By using the same lasers for the manipulation of atoms in both interferometers, the accelerometers have a common reference frame and non-inertial accelerations are effectively rejected as common mode noise in the differential measurement of the gravity gradient. As a result, the dual atom interferometer-based gravity gradiometer allows gravity measurements on a moving platform, while achieving the same long-term stability of the best atomic clocks. In the laboratory-based prototype (Figure 2), the cesium atoms used in each atom interferometer are initially collected and cooled in two separate magneto-optic traps (MOTs). Each MOT, consisting of three orthogonal pairs of counter-propagating laser beams centered on a quadrupole magnetic field, collects up to 10(exp 9) atoms. These atoms are then launched vertically as in an atom fountain by switching off the magnetic field and introducing a slight frequency shift between pairs of lasers to create a moving rest frame for the trapped atoms. While still in this moving-frame molasses, the laser frequencies are further detuned from the atomic resonance (while maintaining this relative frequency shift) to cool the atom cloud's temperature to 2 K or below, corresponding to an rms velocity of less than 2 cm/s. After launch, the cold atoms undergo further state and velocity selection to prepare for atom interferometry. The atom interferometers are then realized using laser-induced stimulated Raman transitions to perform the necessary manipulations of each atom, and the resulting interferometer phase is measured using laser-induced fluorescence for state-normalized detection. More than 20 laser beams with independent controls of frequency, phase, and intensity are required for this measurement sequence. This instrument can facilitate the study of Earth's gravitational field from surface and air vehicles, as well as from space by allowing gravity mapping from a low-cost, single spacecraft mission. In addition, the operation of atom interferometer-based instruments in space offers greater sensitivity than is possible in terrestrial instruments due to the much longer interrogation times available in the microgravity environment. A space-based quantum gravity gradiometer has the potential to achieve sensitivities similar to the GRACE mission at long spatial wavelengths, and will also have resolution similar to GOCE for measurement at shorter length scales.

Yu, Nan↗

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗

Applying the core-satellite species concept: Characteristics of rare and common riverine dissolved organic matter

Introduction: Dissolved organic matter (DOM) composition varies over space and time, with a multitude of factors driving the presence or absence of each compound found in the complex DOM mixture. Compounds ubiquitously present across a wide range of river systems (hereafter termed core compounds) may differ in chemical composition and reactivity from compounds present in only a few settings (hereafter termed satellite compounds). Here, we investigated the spatial patterns in DOM molecular formulae presence (occupancy) in surface water and sediments across 97 river corridors at a continental scale using the “Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems—WHONDRS” research consortium. Methods: We used a novel data-driven approach to identify core and satellite compounds and compared their molecular properties identified with Fourier-transform ion cyclotron resonance mass spectrometry (FT-ICR MS). Results: In this work, we found that core compounds clustered around intermediate hydrogen/carbon and oxygen/carbon ratios across both sediment and surface water samples, whereas the satellite compounds varied widely in their elemental composition. Within surface water samples, core compounds were dominated by lignin-like formulae, whereas protein-like formulae dominated the core pool in sediment samples. In contrast, satellite molecular formulae were more evenly distributed between compound classes in both sediment and water molecules. Core compounds found in both sediment and water exhibited lower molecular mass, lower oxidation state, and a higher degree of aromaticity, and were inferred to be more persistent than global satellite compounds. Higher putative biochemical transformations were found in core than satellite compounds, suggesting that the core pool was more processed. Discussion: The observed differences in chemical properties of core and satellite compounds point to potential differences in their sources and contribution to DOM processing in river corridors. Overall, our work points to the potential of data-driven approaches separating rare and common compounds to reduce some of the complexity inherent in studying riverine DOM.

54 ENVIRONMENTAL SCIENCES↗

Flexible Micropost Arrays for Shear Stress Measurement

Increased fuel costs, heightened environmental protection requirements, and noise abatement continue to place drag reduction at the forefront of aerospace research priorities. Unfortunately, shortfalls still exist in the fundamental understanding of boundary-layer airflow over aerodynamic surfaces, especially regarding drag arising from skin friction. For example, there is insufficient availability of instrumentation to adequately characterize complex flows with strong pressure gradients, heat transfer, wall mass flux, three-dimensionality, separation, shock waves, and transient phenomena. One example is the acoustic liner efficacy on aircraft engine nacelle walls. Active measurement of shear stress in boundary layer airflow would enable a better understanding of how aircraft structure and flight dynamics affect skin friction. Current shear stress measurement techniques suffer from reliability, complexity, and airflow disruption, thereby compromising resultant shear stress data. The state-of-the-art for shear stress sensing uses indirect or direct measurement techniques. Indirect measurements (e.g., hot-wire, heat flux gages, oil interferometry, laser Doppler anemometry, small scale pressure drag surfaces, i.e., fences) require intricate knowledge of the studied flow, restrictive instrument arrangements, large surface areas, flow disruption, or seeding material; with smaller, higher bandwidth probes under development. Direct measurements involve strain displacement of a sensor element and require no prior knowledge of the flow. Unfortunately, conventional "floating" recessed components for direct measurements are mm to cm in size. Whispering gallery mode devices and Fiber Bragg Gratings are examples of recent additions to this type of sensor with much smaller (m) sensor components. Direct detection techniques are often single point measurements and difficult to calibrate and implement in wind tunnel experiments. In addition, the wiring, packaging, and installation of delicate micro-electromechanical devices impede the use of most direct shear sensors. Similarly, the cavity required for sensing element displacement is sensitive to particulate obstruction. This work was focused on developing a shear stress sensor for use in subsonic wind tunnel test facilities applicable to an array of test configurations. The non-displacement shear sensors described here have minimal packaging requirements resulting in minimal or no disturbance of boundary layer flow. Compared to previous concepts, device installation could be simple with reduced cost and down-time. The novelty lies in the creation of low profile (nanoscale to 100 μm) micropost arrays that stay within the viscous sub-layer of the airflow. Aerodynamic forces, which are related to the surface shear stress, cause post deflection and optical property changes. Ultimately, a reliable, accurate shear stress sensor that does not disrupt the airflow has the potential to provide high value data for flow physics researchers, aerodynamicists, and aircraft manufacturers leading to greater flight efficiency arising from more in-depth knowledge on how aircraft design impacts near surface properties.

Wohl, Christopher J.↗

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Soil organic matter diagenetic state informs boreal forest ecosystem feedbacks to climate change

The fate of soil organic carbon (SOC) in boreal forests is dependent on the integrative ecosystem response to climate change. For example, boreal forest productivity is often nitrogen (N) limited, and climate warming can enhance N cycling and primary productivity. However, the net effect of this feedback on the SOC reservoir and its longevity with climate change remain unclear due to difficulty in detecting small differences between large and variable carbon (C) fluxes needed to determine net changes in soil reservoirs. The diagenetic state of SOC – resulting from the physicochemical and biological transformations that alter the original biomolecular composition of detrital inputs to soil over time – is useful for tracing the net response of SOC at the timescales relevant to climate change not usually discernible from fluxes and stocks alone. Here, we test two hypotheses using a mesic boreal forest climate transect: (1) the SOC diagenetic state is maintained across this climosequence, and (2) the maintenance of the SOC diagenetic state is a consequence of coupled soil C and N cycling, signifying the role of enhanced N cycling supporting SOC inputs that maintain SOC stocks within the warmer-climate forests. Shifts in nonvascular to vascular plant inputs with climate observed in these and other boreal forests highlighted the need to carefully separate biogeochemical indicators of SOC source from those signifying diagenetic alteration. We thus evaluated and applied lignin biomarkers to assess the diagenetic alteration of SOC in these boreal forest organic soils and directly compared the lignin diagenetic state with that of soil organic nitrogen (SON) assessed through amino acid composition. The lignin diagenetic state remained constant across the climate transect, indicating a balance between the input and removal of lignin in these mesic boreal forests. When combined with previous knowledge of these forest ecosystems, including the diagenetic state of SON and direct measures of C fluxes and stocks, the results indicate a coupled increase in C and N cycling with climate warming that supports forest productivity and maintains SOC stocks. This balance could markedly shift as other factors begin to limit forest productivity (e.g., trace nutrients, water) with further climate change or affect forest nutrient allocation (e.g., forest age or compositional change). Further application of the approach presented here could be used to detect the limits of this and other ecosystem–climate feedbacks, by providing a tractable and parameterizable index of the lignin state across large spatial scales, necessary for ecosystem-scale parameterizations.

54 ENVIRONMENTAL SCIENCES↗

The X-43A Flight Research Program: Lessons Learned on the Road to Mach 10

During an aerospace engineer's undergraduate studies, he or she will attend classes in aerodynamics, thermodynamics, structures, stability and control, dynamics, design, propulsion, and computer science, along with the related courses in mathematics, physics, statistics, and chemistry required to understand the material. Upon graduation, the new engineer will have acquired a basic knowledge of how to build an aerospace vehicle. What only comes through experience, however, is the understanding of the inevitable imperfect process through which an aerospace vehicle is built. This is the adventure of turning a basic concept into functional hardware. Engineers working on a project must often deal with ambiguous situations. They are routinely asked by management to provide risk assessments of a project, yet even after careful analysis uncertainties remain. The project must be accomplished within finite limits of time and money. The question an engineer answers is whether the solution to potential problem is worth the cost and schedule delay, or if the solution might actually be worse than the problem it is meant to solve. Review protocols are established to ensure that an unknown has not been overlooked. But these cannot protect against an unknown unknown. Examples of these situations can be found in the history of the X-43A Hyper-X (Hypersonic Experiment) program. In this NASA project, a supersonic combustion ramjet (scramjet) engine was flight tested on a subscale vehicle. The X-43A Hyper-X Research Vehicle (HXRV) was launched from a B-52B mothership, then boosted to the test speed by a modified Pegasus rocket first stage, called the Hyper-X Launch Vehicle (HXLV). Once at the proper speed and altitude, the X-43A separated from the booster, stabilized itself, and then the engine test began. Although wind-tunnel scramjet engine tests had begun in the late 1950s, before the Hyper-X program there had never been an actual in-flight test of such an engine integrated with an appropriate airframe. Thus, while the scramjet had successfully operated in the artificial airflow of wind tunnels, the concept had yet to be proven in real air. These conditions meant changes in density and temperature, as well as changes in angle of attack and sideslip of a free-flying vehicle. A wind tunnel is limited in its ability to simulate these subtle factures, which have a major impact on almost any vehicle, but especially that of a scramjet's performance. The Hyper-X project was to provide a real-world benchmark of the ground test data. The full scale X-43A engine would be operated in the wind tunnel, and then flown, and the data from its operation would then be compared with projections. If these matched, the wind tunnel data would be considered a reliable design tool for future scramjet. If there were significant differences, the reasons for these would have to be identified. Until such information was available, scramjets would lack the technological maturity to be considered for future space launch or high-speed atmospheric flight vehicles.

Peebles, Curtis↗

Commercial Smallsat Data Acquisition Program On-ramp #2 Airbus U.S. Synthetic Aperture Radar (SAR) Evaluation Report

In 2017, NASA’s Earth Science Division (ESD) launched the Private-Sector Small Constellation Satellite Data Product Pilot, now referred to as the Commercial Smallsat Data Acquisition (CSDA) program. The objective of CSDA is to identify, evaluate, and acquire commercial remote sensing data that support NASA’s Earth science research and application activities. The Pilot successfully concluded in early 2020, when CSDA transitioned into a sustained program with on-ramping opportunities for new vendors as the industry emerges with new candidates and capabilities. In October 2019, a Request for Information (RFI) seeking capability statements from parties interested in providing data from spaceborne platforms was released for the CSDA on-ramp #2 evaluations. To be responsive to the RFI, the commercial satellite constellations had to consist of three or more operating spacecraft actively collecting data in a non-geostationary orbit with full latitudinal coverage and be U.S. companies. Two vendors responded to the RFI and were evaluated by a committee composed of NASA ESD leadership, program managers, and scientists. Both vendors satisfied the RFI requirements and were asked to respond to a Request for Proposal (RFP). After review of the proposals, NASA entered into a Blanket Purchase Agreement (BPA) with Airbus Defense and Space GEO, Inc. (Airbus) U.S. in September 2021 and with BlackSky Geospatial Solutions, Inc. (BlackSky) in November 2021. In this report, CSDA provides an evaluation of the usefulness of data provided by the Airbus U.S. Synthetic Aperture Radar (SAR) satellite constellation, consisting of TerraSAR-X (launched in 2007), TanDEM-X (launched in 2010), and PAZ (launched in 2018), for advancing NASA’s Earth system science research and applications. The evaluation of the BlackSky commercial data will be provided in a separate report. To conduct the Airbus evaluation, NASA’s ESD augmented 13 existing research projects that could potentially benefit from, and had the expertise to evaluate, the commercial data being considered for longer-term purchase. Investigators from NASA’s Research and Analysis Program science focus areas and from NASA’s Applied Sciences Program elements participated in the evaluation. A summary of the research areas evaluated by the Principal Investigator (PI) teams is presented in Figure 3. CSDA also funded a dedicated activity to evaluate the satellite data quality (calibration and geolocation) independently by assessing the accuracy of data from Airbus. Evaluation activities were carried out by the selected PIs from December 7, 2022, to December 7, 2023. Delivery of datasets requested by the researchers began in January 2023. The vendors were evaluated on the accessibility of data, accuracy and completeness of metadata, and promptness and quality of user support services. Datasets purchased during the evaluation have been archived by NASA and will be made available to current and future government-funded researchers in accordance with the End User License Agreement (EULA). This synthesis report distills and integrates the findings of research reports commissioned by NASA for the Airbus evaluation. This report also includes recommendations that inform the way ahead for the program. The scientific results from the evaluations demonstrated that the commercial data from Airbus were able to advance NASA research and applications. However, the PIs encountered limitations that diminished the usefulness of the data due to the amount of effort that was required to access, preprocess, and analyze these data. One significant issue encountered was the limited spatial and temporal coverage of the data in the Airbus archive that could be used to conduct time series analyses or assessments over large spatial scales. Overall, however, the utility and the quality of the evaluated data outweighed the difficulties encountered, and NASA has concluded that the Airbus SAR data would complement NASA’s existing Earth observation capabilities and Airbus U.S. would qualify to participate in the sustained phase of the program.

Batuhan Osmanoglu↗

Amorphous Phases on the Surface of Mars

Both primary (volcanic/impact glasses) and secondary (opal/silica, allophane, hisingerite, npOx, S-bearing) amorphous phases appear to be major components of martian surface materials based on orbital and in-situ measurements. A key observation is that whereas regional/global scale amorphous components include altered glass and npOx, local scale amorphous phases include hydrated silica/opal. This suggests widespread alteration at low water-to-rock ratios, perhaps due to snow/ice melt with variable pH, and localized alteration at high water-to-rock ratios. Orbital and in-situ measurements of the regional/global amorphous component on Mars suggests that it is made up of at least three phases: npOx, amorphous silicate (likely altered glass), and an amorphous S-bearing phase. Fundamental questions regarding the composition and the formation of the regional/global amorphous component(s) still remain: Do the phases form locally or have they been homogenized through aeolian activity and derived from the global dust? Is the parent glass volcanic, impact, or both? Are the phases separate or intimately mixed (e.g., as in palagonite)? When did the amorphous phases form? To address the question of source (local and/or global), we need to look for variations in the different phases within the amorphous component through continued modeling of the chemical composition of the amorphous phases in samples from Gale using CheMin and APXS data. If we find variations (e.g., a lack of or enrichment in amorphous silicate in some samples), this may imply a local source for some phases. Furthermore, the chemical composition of the weathering products may give insight into the formation mechanisms of the parent glass (e.g., impact glasses contain higher Al and lower Si [30], so we might expect allophane as a weathering product of impact glass). To address the question of whether these phases are separate or intimately mixed, we need to do laboratory studies of naturally altered samples made up of mixed phases (e.g., palagonite) and synthetic single phases to determine their short-range order structures and calculate their XRD patterns to use in models of CheMin data. Finally, to address the timing of the alteration, we need to study rocks on the martian surface of different ages that may contain glass (volcanic or impact) with MSL and future rovers to better understand how glass alters on the martian surface, if that alteration mechanism is universal, and if alteration spans across long periods of time or if there is a time past which unaltered glass remains.

Rampe, E. B.↗

Ambulation Increases Decompression Sickness in Spacewalk Simulations

Musculoskeletal activity has the potential to both improve and compromise decompression safety. Exercise enhances inert gas elimination during oxygen breathing prior to decompression (prebreathe), but it may also promote bubble nuclei formation (nucleation), which can lead to gas phase separation and bubble growth and increase the risk of decompression sickness (DCS). The timing, pattern and intensity of musculoskeletal activity and the level of tissue supersaturation may be critical to the net effect. Understanding the relationships is important to evaluate exercise prebreathe protocols and quantify decompression risk in gravity and microgravity environments. Data gathered during NASA's Prebreathe Reduction Program (PRP) studies combined oxygen prebreathe and exercise followed by low pressure (4.3 psi; altitude equivalent of 30,300 ft [9,235 m]) microgravity simulation to produce two protocols used by astronauts preparing for extravehicular activity. Both the Phase II/CEVIS (cycle ergometer vibration isolation system) and ISLE (in-suit light exercise) trials eliminated ambulation to more closely simulate the microgravity environment. The CEVIS results (35 male, 10 female) serve as control data for this NASA/Duke study to investigate the influence of ambulation exercise on bubble formation and the subsequent risk of DCS. METHODS Four experiments will replicate the CEVIS exercise-enhanced oxygen prebreathe protocol, each with a different exception. The first of these is currently underway. Experiment 1 - Subjects complete controlled ambulation (walking in place with fixed cadence and step height) during both preflight and at 4.3 psi instead of remaining nonambulatory throughout. Experiment 2 - Subjects remain non-ambulatory during the preflight period and ambulatory at 4.3 psi. Experiment 3 - Subjects ambulate during the preflight period and remain non-ambulatory at 4.3 psi. Experiment 4 - The order of heavy and light exercise employed in the CEVIS protocol is reversed, with the light exercise occurring first (subjects remain non-ambulatory throughout). Decompression stress is assessed with non-invasive ultrasound during each of 14 epochs of a 4 hour simulated spacewalk at 4.3 psi; aural Doppler is used to monitor bubbles (Spencer grade 0-IV scale) passing through the pulmonary artery, and two-dimensional echocardiographic imaging is used to look for left ventricular gas emboli (LVGE; the presence of which is a test termination criterion). Venous blood is collected at baseline and twice following repressurization to determine if the decompression stress is correlated with microparticles (cell fragments) accumulation. The plan is to test 25-50 subjects in each experiment. Fisher Exact Tests (one-tailed) are used to compare test and control groups. Trials are suspended when the DCS or grade IV VGE observations reach 70% confidence of DCS risk >15% and grade IV VGE risk >20%. RESULTS Experiment 1 was concluded with 20 complete trials (15 male, 5 female) since the statistical outcome would not change with five additional trials. The observed DCS was significantly greater in Experiment 1 than in CEVIS trials (4/20 [20%] vs. 0/45 [0%], respectively, p=0.007), as was the frequency of peak grade IV VGE (6/21 [29%; including one additional subject that presented grade IV VGE but whose trial was ended before completion when LVGE were observed] vs. 3/45 [7%], respectively, p=0.024). Experiment 3 trials are now underway, with 11 trials completed (10 male, 1 female). Preliminary results indicate no difference in observed DCS between Experiment 3 and CEVIS trials (1/11 [9%] vs. 0/45 [0%], respectively, p=0.196), or between Experiment 3 and Experiment 1 trials (p=0.405). The frequency of peak grade IV VGE in Experiment 3 (2/11 [18%]) did not differ from CEVIS or Experiment 1 trials (p=0.251 and p=0.425, respectively). Microparticle patterns are widely variable and still under analysis. DISCUSSION The results of the Experiment 1 trials support the thesis that decompression stress is increased by ambulation exercise, given the higher incidence of DCS and grade IV VGE when compared to the non-ambulatory PRP CEVIS trials. Experiment 3 trials are incomplete, but suggest that the effect of ambulation during ground level preflight oxygen breathing alone, when subjects are undersaturated with inert gas, may not differ in risk from ambulation at both preflight and spacesuit pressures, the latter when subjects are supersaturated with inert gas. Further trials are needed to confirm the relative effects of ambulation in undersaturated vs. supersaturated states and to determine whether light exercise facilitates the removal of heavy exercise-induced nucleation (Experiment 4).

Pollock, N. W.↗

An Experimental Study of Small-Scale Variability of Raindrop Size Distribution

An experimental study of small-scale variability of raindrop size distributions (DSDs) has been carried out at Wallops Island, Virginia. Three Joss-Waldvogel disdrometers were operated at a distance of 0.65, 1.05, and 1.70 km in a nearly straight line. The main purpose of the study was to examine the variability of DSDs and its integral parameters of liquid water content, rainfall, and reflectivity within a 2-km array: a typical size of Cartesian radar pixel. The composite DSD of rain events showed very good agreement among the disdrometers except where there were noticeable differences in midsize and large drops in a few events. For consideration of partial beam filling where the radar pixel was not completely covered by rain, a single disdrometer reported just over 10% more rainy minutes than the rainy minutes when all three disdrometers reported rainfall. Similarly two out of three disdrometers reported5%more rainy minutes than when all three were reporting rainfall. These percentages were based on a 1-min average, and were less for longer averaging periods. Considering only the minutes when all three disdrometers were reporting rainfall, just over one quarter of the observations showed an increase in the difference in rainfall with distance. This finding was based on a 15-min average and was even less for shorter averaging periods. The probability and cumulative distributions of a gamma-fitted DSD and integral rain parameters between the three disdrometers had a very good agreement and no major variability. This was mainly due to the high percentage of light stratiform rain and to the number of storms that traveled along the track of the disdrometers. At a fixed time step, however, both DSDs and integral rain parameters showed substantial variability. The standard deviation (SD) of rain rate was near 3 mm/h, while the SD of reflectivity exceeded 3 dBZ at the longest separation distance. These standard deviations were at 6-min average and were higher at shorter averaging periods. The correlations decreased with increasing separation distance. For rain rate, the correlations were higher than previous gauge-based studies. This was attributed to the differences in data processing and the difference in rainfall characteristics in different climate regions. It was also considered that the gauge sampling errors could be a factor. In this regard, gauge measurements were simulated employing existing disdrometer dataset. While a difference was noticed in cumulative distribution of rain occurrence between the simulated gauge and disdrometer observations, the correlations in simulated gauge measurements did not differ from the disdrometer measurements.

Tokay, Ali↗

Attitude Ground System (AGS) For The Magnetospheric Multi-Scale (MMS) Mission

The Magnetospheric Multiscale (MMS) mission is a Solar-Terrestrial Probe mission consisting of four identically instrumented spin-stabilized spacecraft flying in an adjustable pyramid-like formation around the Earth. The formation of the MMS spacecraft allows for three-dimensional study of the phenomenon of magnetic reconnection, which is the primary objective of the mission. The MMS spacecraft were launched early on March 13, 2015 GMT. Due to the challenging and very constricted attitude and orbit requirements for performing the science, as well as the need to maintain the spacecraft formation, multiple ground functionalities were designed to support the mission. These functionalities were incorporated into a ground system known as the Attitude Ground System (AGS). Various AGS configurations have been used widely to support a variety of three-axis-stabilized and spin-stabilized spacecraft missions within the NASA Goddard Space Flight Center (GSFC). The original MMS operational concept required the AGS to perform highly accurate predictions of the effects of environmental disturbances on the spacecraft orientation and to plan the attitude maneuvers necessary to stay within the science attitude tolerance. The orbit adjustment requirements for formation control drove the need also to perform calibrations that have never been done before in support of NASA GSFC missions. The MMS mission required support analysts to provide fast and accurately calibrated values of the inertia tensor, center of mass, and accelerometer bias for each MMS spacecraft. During early design of the AGS functionalities, a Kalman filter for estimating the attitude, body rates, center of mass, and accelerometer bias, using only star tracker and accelerometer measurements, was heavily analyzed. A set of six distinct filters was evaluated and considered for estimating the spacecraft attitude and body rates using star tracker data only. Four of the six filters are closely related and were compared during support of the Time History of Events and Macroscale Interactions during Substorms (THEMIS) and Space Technology-5 (ST-5) missions. These analyses exposed high dependency and sensitivity on the knowledge of the spacecraft inertia tensor for both body rates and accelerometer bias estimation. The conclusion of the analysis led to the design of an inertia tensor calibration technique using only star tracker data. The second most important result of the analysis was the design of two separate Kalman filters to estimate the spacecraft attitude and body rates and the accelerometer bias instead of a single combined filter. In this paper, the calibration results of the mass properties, as well as the performance of the spacecraft attitude and body rates filters using flight data are presented and compared against the mission requirements.

Estimation↗

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗

Data-informed grid refinement to improve traveltime accuracy in the regional seismic traveltime (RSTT) model

The regional seismic traveltime (RSTT) model predicts traveltimes of regional seismic phases accounting for 3-D structure of the crust and the upper mantle on a global scale. Previous versions of the RSTT model have been implemented using nodes separated by ∼1° spacing across the globe. A regional-scale study using regional Pn and Pg traveltimes across Israel and the Middle East demonstrated that data driven, systematic grid refinement reduces traveltime residuals and enhances resolution of smaller tectonic features in regions having dense ray coverage. High density Pn ray coverage in the western US, Europe, Middle East and East Asia can likewise provide the resolution that allows systematic global grid refinement of the RSTT model. In this study, we use a large number of Pn ray paths originating from events located with an epicentral location uncertainty of 25 km (GT25) or better. We conduct targeted grid refinements at 1.0°, 0.5°, 0.25° and 0.125° on a global scale, producing a refined RSTT model that yields a 21.6 per cent reduction in median event location error in Europe and the Middle East, when compared with the original global RSTT model presented in Begnaud et al. The new model also resolves finer tectonic structures in regions with high Pn ray density.

58 GEOSCIENCES↗