Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “current distribution mapping”

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 325 records · Page 18

E-Beam Irradiation Beamline at Jefferson Lab for Wastewater Treatment

Electron beam (e-beam) irradiation is an advanced method for removing the emerging organic contaminants in wastewater, such as 1,4-dioxane, and per- and polyfluoroalkyl sub?stances (PFAS), which have the potential to be regulated in the near future but are difficult to degrade by conventional treatment methods. This dissertation was based on a collabo?rating study between Jefferson Lab (Jlab) and Hampton Roads Sanitation District (HRSD) to evaluate the impacts of e-beam irradiation on the degradation of 1,4-dioxane and PFAS. It presents the detailed design of an e-beam irradiation beamline, the beamline simulation using general particle tracer (GPT), the installation, and the successful commissioning. The beamline is an extended section to the Upgraded Injector Test Facility (UITF) at Jlab, which is a continuous wave (CW), superconducting radio frequency (SRF) linear accelerator with an energy up to 10 MeV. To let the sample volume receive a significant and uniform dose distribution, the sample target was designed by investigating the dose distribution simulated with a Monte-Carlo code, FLUKA. The optimized beam energy was 8 MeV and the beam current for irradiation was ? 100 nA. With the use of a solenoid and a raster, the electron beam was expanded to the optimized size, 1? ? 15 mm, at the sample target. This dis?sertation also presents the irradiation experiments for 64 1,4-dioxane samples and 48 PFAS samples in different water matrices, including real wastewater, which were prepared and an?alyzed by HRSD. The absorbed doses of samples were obtained by designing a new method that combines the FLUKA simulation and dose mapping experiments. The irradiation re?sults showed that around 95% of 1,4-dioxane was removed at a low dose of ? 2 kGy, and PFAS was possibly removed at a lower dose compared to similar research. This work is a first implementation of SRF technology for wastewater irradiation. It pro?vides a valuable new research platform for wastewater remediation. This will be in support of future industrial applications of advanced SRF technologies and better understanding of the effect of e-beam irradiation on wastewater or water treatment.

Li, Xi↗

Monitoring Airspace Complexity and Determining Contributing Factors

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

complexity↗

Characteristic structures of the highland boundary on Mars: Evidence against a single mega-impact event?

It has been suggested that an early mega-impact event might explain the fumdamental crustal dichotomy on Mars. Detailed morphological mapping of the structures which characterize the boundary between the cratered highlands and northern plains does not support this idea: the distribution of these features along and especially away from the boundary is more consistent with a larger number of smaller but overlapping impacts. A data base was assembled for the study of the highland boundary through quantitative mapping (including crater counts) based on the 1:2 M controlled photomosaics. Whole and partial craters larger than 10 km in diameter, knobby terrain, detached plateaus and intervening plains-forming units have been identified wherever they occur on Mars between +65 and -45 deg latitude. There is evidence for occurrences of possible old cratered terrain (or its remanents) at high latitudes well north of the current or likely former location of the proposed Borealis Basin rim which are difficult to reconcile with an impact structure of the proposed size and location.

Semeniuk, Ann Marie↗

Openet: Applications of Satellite-Based Evapotranspiration Data for Water Resources Management in the Western United States

Advancing water security in overallocated river basins globally requires consistent and reproducible information on consumptive use of water that can anchor the development of data-driven solutions to the challenge of balancing water supply and demand. OpenET is a fully automated system for field-scale (30 m), satellite-based mapping of evapotranspiration (ET) at daily, monthly and annual timesteps. OpenET currently provides spatially contiguous data throughout the 23 westernmost states in the continental US, and includes both current information as well as multi-year timeseries of ET. The OpenET consortium has implemented an ensemble of satellite-based ET models (ALEXI/DisALEXI, eeMETRIC, PT-JPL, geeSEBAL, SIMS and SSEBop) on Google Earth Engine, which provides a shared computing platform for collaboration on processing of data from Landsat and other satellites, land cover and meteorological inputs, leading to increased consistency and accuracy across the ensemble of models. Earth Engine also facilitates hosting and distribution of data via open data collections and an application programming interface. We provide updates on the OpenET framework, open data services and data access tools, approach to geographic expansion, recent accuracy assessments, and describe how a user-driven design approach has facilitated successful applications of OpenET data for a wide range of water resource management activities. Applications to date include: use of ET data to improve quantification of ET and consumptive use in Oregon, Utah and the Upper Colorado River Basin; streamlining of water use reporting requirements in the California Delta; support for calculation of water budgets for the implementation of the Sustainable Groundwater Management Act in California; and integration into decision support tools for irrigation management. The use cases demonstrate how satellite-derived ET data that are easily accessed and seen as broadly accepted can accelerate adoption of innovative water management practices at scale, and support advances in the sustainability of water supplies. Uptake and use of data by the OpenET science community has also led to advances in our understanding of the impacts of landcover change, irrigation intensification and wildfire events on hydrology and the water security.

Applications↗

King Ghidorah Supercluster: Mapping the light and dark matter in a new supercluster at z = 0.55 using the subaru hyper suprime-cam

ABSTRACT This paper reports our discovery of the most massive supercluster, termed the King Ghidorah Supercluster (KGSc), at z = 0.50–0.64 in the Third Public Data Release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP PDR3) over 690 deg2, as well as an initial result for a galaxy and dark matter mapping. The primary structure of the KGSc comprises triple broad weak-lensing (WL) peaks over 70 comoving Mpc. Such extensive WL detection at z > 0.5 can only currently be achieved using the wide-field high-quality images produced by the HSC-SSP. The structure is also contiguous with multiple large-scale structures across a ∼400 comoving Mpc scale. The entire field has a notable overdensity (δ = 14.7 ± 4.5) of red-sequence clusters. Additionally, large-scale underdensities can be found in the foreground along the line of sight. We confirmed the overdensities in stellar mass and dark matter distributions to be tightly coupled and estimated the total mass of the main structure to be 1 × 1016 solar masses, according to the mock data analyses based on large-volume cosmological simulations. Further, upcoming wide-field multi-object spectrographs such as the Subaru Prime Focus Spectrograph may aid in providing additional insights into distant superclusters beyond the 100 Mpc scale.

Shimakawa, Rhythm (ORCID:0000000344422750)↗

Global Bathymetric Prediction For Ocean Modeling and Marine Geophysics

We proposed to construct a complete bathymetric map of the oceans at a 3-10 km resolution by combining all of the available depth soundings collected over the past 30 years with high resolution marine gravity information provided by the Geosat, ERS-1/2, and Topex/Poseidon altimeters. Detailed bathymetry is essential for understanding physical oceanography and marine geophysics. Currents and tides are controlled by the overall shapes of the ocean basins as well as the smaller sharp ocean ridges and seamounts. Because erosion rates are low in the deep oceans, detailed bathymetry reveals the mantle convection patterns, the plate boundaries, the cooling/subsidence of the oceanic lithosphere, the oceanic plateaus, and the distribution of off-ridge volcanoes. We proposed to: (1) Accumulate all available depth soundings collected over the past 30 years; (2) Use the short wavelength (< 160 km) satellite gravity information to interpolate between sparse ship soundings; (3) Improve the resolution of the marine gravity field using enhanced estimates along repeat altimeter profiles together with the dense altimeter measurements; (4) Refine/improve bathymetric predictions using the improved resolution gravity field and also by investigating computer-intensive methods for bathymetric prediction such as inverse theory; and (5) Produce a 'Globe of the Earth' similar to the globe of Venus prepared by the NASA Magellan investigation. This will also include the best available digital land data.

Sandwell, David T.↗

Conformal maps of viscous electron flow in the Gurzhi crossover

Herein we investigate the impact of geometric constriction on the viscous flow of electron liquid through quantum point contacts. We provide analysis on the electric potential distribution given the setup of a slit configuration and use the method of conformal mapping to obtain analytical results. The potential profile can be tested and contrasted experimentally with the scanning tunneling potentiometry technique. We discuss intricate physics that underlies the Gurzhi effect, i.e., the enhancement of conductivity in the viscous flow, and compare results for different boundary conditions. In addition, we calculate the temperature dependence of the momentum relaxation time as a result of impurity assisted quasiballistic interference effects and discuss various correlational corrections that lead to the violation of Matthiessen's rule in the hydrodynamic regime. We caution that spatially inhomogeneous profiles of current in the Gurzhi crossover between Ohmic and Stokes flows might also appear in the nonhydrodynamic limit where nonlocality plays an important role. This conclusion is corroborated by calculation of dispersive conductivity in the weakly impure limit.

2-dimensional systems↗

California coast nearshore processes study

The author has identified the following significant results. Large scale sediment plumes from intermittent streams and rivers form detectable seasonal patterns on ERTS-1 imagery. The ocean current systems, as plotted from three California coast ERTS mosaics, were identified. Offshore patterns of sediment in areas such as the Santa Barbara Channel are traceable. These patterns extend offshore to heretofore unanticipated ranges as shown on the ERTS-1 imagery. Flying spot scanner enhancements of NASA tapes resulted in details of subtle and often invisible (to the eye) nearshore features. The suspended sediments off San Francisco and in Monterey Bay are emphasized in detail. These are areas of extremely changeable offshore sediment transport patterns. Computer generated contouring of radiance levels resulted in maps that can be used in determining surface and nearsurface suspended sediment distribution. Tentative calibrations of ERTS-1 spectral brightness against sediment load have been made using shipboard measurements. Information from the combined enhancement and interpretation techniques is applicable to operational coastal engineering programs.

Pirie, D. M.↗

Integrated health monitoring and controls for rocket engines

Current research in intelligent control systems at the Lewis Research Center is described in the context of a functional framework. The framework is applicable to a variety of reusable space propulsion systems for existing and future launch vehicles. It provides a 'road map' technology development to enable enhanced engine performance with increased reliability, durability, and maintainability. The framework hierarchy consists of a mission coordination level, a propulsion system coordination level, and an engine control level. Each level is described in the context of the Space Shuttle Main Engine. The concept of integrating diagnostics with control is discussed within the context of the functional framework. A distributed real time simulation testbed is used to realize and evaluate the functionalities in closed loop.

Merrill, W. C.↗

Integrated health monitoring and controls for rocket engines

Current research in intelligent control systems at the Lewis Research Center is described in the context of a functional framework. The framework is applicable to a variety of reusable space propulsion systems for existing and future launch vehicles. It provides a 'road map' technology development to enable enhanced engine performance with increased reliability, durability, and maintainability. The framework hierarchy consists of a mission coordination level, a propulsion system coordination level, and an engine control level. Each level is described in the context of the Space Shuttle Main Engine. The concept of integrating diagnostics with control is discussed within the context of the functional framework. A distributed real time simulation testbed is used to realize and evaluate the functionalities in closed loop.

Merrill, W. C.↗

Estimates of the Lightning NOx Profile in the Vicinity of the North Alabama Lightning Mapping Array

The NASA Marshall Space Flight Center Lightning Nitrogen Oxides Model (LNOM) is applied to August 2006 North Alabama Lightning Mapping Array (NALMA) data to estimate the (unmixed and otherwise environmentally unmodified) vertical source profile of lightning nitrogen oxides, NOx = NO + NO2. Data from the National Lightning Detection Network (Trademark) (NLDN) is also employed. This is part of a larger effort aimed at building a more realistic lightning NOx emissions inventory for use by the U.S. Environmental Protection Agency (EPA) Community Multiscale Air Quality (CMAQ) modeling system. Overall, special attention is given to several important lightning variables including: the frequency and geographical distribution of lightning in the vicinity of the NALMA network, lightning type (ground or cloud flash), lightning channel length, channel altitude, channel peak current, and the number of strokes per flash. Laboratory spark chamber results from the literature are used to convert 1-meter channel segments (that are located at a particular known altitude; i.e., air density) to NOx concentration. The resulting lightning NOx source profiles are discussed.

Koshak, William J.↗

The PISCES 2 parallel programming environment

PISCES 2 is a programming environment for scientific and engineering computations on MIMD parallel computers. It is currently implemented on a flexible FLEX/32 at NASA Langley, a 20 processor machine with both shared and local memories. The environment provides an extended Fortran for applications programming, a configuration environment for setting up a run on the parallel machine, and a run-time environment for monitoring and controlling program execution. This paper describes the overall design of the system and its implementation on the FLEX/32. Emphasis is placed on several novel aspects of the design: the use of a carefully defined virtual machine, programmer control of the mapping of virtual machine to actual hardware, forces for medium-granularity parallelism, and windows for parallel distribution of data. Some preliminary measurements of storage use are included.

Pratt, Terrence W.↗

Applicability of Micro X-Ray Fluorescence Spectroscopy to Astromaterials Curation and Research

Introduction: The Astromaterials Acquisition and Curation Office at NASA’s Johnson Space Center (JSC) curates NASA’s astromaterial sample collections which includes: Apollo samples, Luna samples, Ant-arctic meteorites, cosmic dust particles, microparticle impacts into space-flown materials, Genesis solar wind atoms, Stardust comet Wild-2 particles, Stardust inter-stellar particles, Hayabusa asteroid Itokawa particles, Hayabusa 2 asteroid Ryugu particles, and future OSIRIS-Rex asteroid Bennu particles (landing in Sep-tember, 2023) [1–3]. To enhance JSC’s advanced cu-ration capabilities, we have recently installed a high-performance micro-X-ray fluorescence (µXRF) spec-trometer to assist in sample characterization through rapid, non-destructive, in-situ elemental analyses that do not require the sample preparation protocols (i.e., polishing and carbon-coating) commonly needed for electron beam analyses. With this new instrument, we are capable of detecting all elements down to carbon in a variable-pressure or He-purged chamber for anal-ysis of a wide range of sample types. Here we describe the instrumental set-up, capabilities, and applicability of µXRF analysis to astromaterials curation and re-search. Instrumentation and Methodology: The X-ray fluorescence and computed tomography lab (X-FaCT) lab at JSC is now equipped with a Bruker M4 Tornado Plus µXRF (Fig. 1). This system is an energy-dispersive x-ray spectrometer equipped with two 60 mm2 silicon drift detectors (SDD) that are able to be used simulta-neously for output count rates ~500,000 cps. New light element windows allow detecting and analyzing the entire elemental range from carbon to americium. Two x-ray tubes (micro-focus Rh with polycapillary lenses and W with collimators of 0.5, 1.0, 2.0, and 4.5 mm) with max excitation parameters of 50 kV, 30 W and 50 kV, 40 W, respectively, allow for more flexibility of the analysis of high energy lines. The motorized X-Y-Z stage has a mapping range of 190 x 160 mm and can support samples up to 7 kg (~15.5 lbs) and a height of 120 mm [4]. Analytical modes include elemental analysis (down to ~20 µm spot size) via point, line, or area of bulk materials (rock surfaces, thin sections, thick sections, etc.) as well as coating analysis (determination of thickness and composition) of samples. This system has a variable vacuum chamber (1 mbar to 1 atm) that is also equipped with a He-purge system which accommodates vacuum sensitive samples while still allowing detection of light elements at atmospheric pressure. Utility and Applicability of µXRF in Astro-materials research and exploration science (ARES): Elemental analysis using µXRF is commonly em-ployed for both terrestrial and planetary geological science disciplines [5]. It is especially useful for analy-sis of astromaterials given the limited sample prepara-tion required, which is not feasible for certain materi-als. Here we show select applications of µXRF anal-yses of astromaterials that can, have, and will be done at JSC’s X-FaCT lab. Point analysis: In-situ spot analyses (~20 µm spot size) on a cut slab of Martian meteorite NWA 10922 allowed for the discovery, qualitative elemental analy-sis and determination of different feldspar minerals [6]. These point analyses served as an effective prelim-inary step for subsequent quantitative analyses. Ana-lytical standards can be employed for more accurate quantification of µXRF spot analyses. Area analysis: This analytical mode measures all detectable elements (from C to Am) at each pixel (>5 µm pixel size) in a user-defined area. The results are shown as elemental maps which can be extracted as 16-bit TIFF’s for further data processing. In Fig. 2. we show elemental distribution maps of the high-Ti basalt 73001,531 that have been processed using ImageJ software. From these maps you can accurately and quickly (this map took ~50 mins.) identify mineral components, such as pyroxene, plagioclase, oxides, and phosphates, compositional zoning, and mineral textures. Detection of high-Z phases: µXRF techniques are es-pecially effective at analyzing trace minerals with high-atomic-number (high-Z) elements because the high-energy characteristic X-rays used (relative to SEM EDS) allow for mapping of K lines in elements up to La (typically SEM maps use L X-ray lines for elements >Zn, and these can often have interferences). Thus, µXRF is especially suited for identifying minerals like zircon, baddeleyite, REE-rich phosphates, Fe-rich met-als, oxides, sulfides, and phosphides [4]. In Fig. 3 we show elemental distribution maps for 73001,530 where we are able to correlate the original video image with, Zr, Si, Y and Hf elemental maps together identi-fying the location of a zircon. In this location you would expect lower Si compared to surrounding mate-rial, as well as higher Zr, Y, and Hf content compared to surrounding material, all of which is confirmed by our XRF ele-mental distribution maps (Figure 2.) Conclusions: The new M4 Tornado Plus µXRF within the Astromaterials Acquisition and Curation office at NASA JSC allows for rapid and non-destructive elemental analysis of astromaterials with limited or no sample preparation. µXRF analyses pro-vide crucial compositional knowledge for the prelimi-nary examination and curation of astromaterials. This instrument enhances the advanced curation capabili-ties in the X-FaCT laboratory at JSC by allowing pro-ductive, cohesive, and non-destructive multi-modal x-ray analyses on astromaterial samples, which is neces-sary for the comprehensive curation and study of our current and future astromaterial collections. Addition-ally, µXRF can provide complimentary information to researchers for studies on astromaterials. References: [1] Allen, C. et al., (2011). Chemie De Erde Geochemistry, 71, 1-20. [2] McCubbin, F. M. et al., (2016) 47th LPSC, abstract #2668 [3] Zeigler, R. A. et al., (2017) 48th LPSC, abstract #2772 [4] Bruker User Manual [5] Young et al., (2016) Appl. Geochemistry, 72, 77-87 [6] Mor-ris, R. V. et al., (2023) 54th LPSC.

E W O'Neal↗

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Remote Sensing of Global Fire Patterns, Aerosol Optical Thickness, and Carbon Monoxide During April 1994

Fires play a crucial role in several ecosystems. They are routinely used to burn forests in order to accommodate the needs of the expanding population, clear land for agricultural purposes, eliminate weeds and pests, regenerate nutrients in grazing and crop lands and produce energy for cooking and heating purposes. Most of the fires on earth are related to biomass burning in the tropics, although they are not confined to these latitudes. The boreal and tundra regions also experience fires on a yearly basis. The current study examines global fire patterns, Aerosol Optical Thickness (AOT) and carbon monoxide concentrations during April 9-19, 1994. Recently, global Advanced Very High Resolution Radiometer (AVHRR) data at nadir ground spatial resolution of 1 km are made available through the NASA/NOAA Pathfinder project. These data from April 9-19, 1994 are used to map fires over the earth. In summary, our analysis shows that fires from biomass burning appear to be the dominant factor for increased tropospheric CO concentrations as measured by the MAPS. The vertical transport of CO by convective activities, along with horizontal transport due to the prevailing winds, are responsible for the observed spatial distribution of CO.

Christopher, Sundar A.↗

Estimates of the Lightning NOx Profile in the Vicinity of the North Alabama Lightning Mapping Array

The NASA Marshall Space Flight Center Lightning Nitrogen Oxides Model (LNOM) is applied to August 2006 North Alabama Lightning Mapping Array (LMA) data to estimate the raw (i.e., unmixed and otherwise environmentally unmodified) vertical profile of lightning nitrogen oxides, NOx = NO + NO 2 . This is part of a larger effort aimed at building a more realistic lightning NOx emissions inventory for use by the U.S. Environmental Protection Agency (EPA) Community Multiscale Air Quality (CMAQ) modeling system. Data from the National Lightning Detection Network TM (NLDN) is also employed. Overall, special attention is given to several important lightning variables including: the frequency and geographical distribution of lightning in the vicinity of the LMA network, lightning type (ground or cloud flash), lightning channel length, channel altitude, channel peak current, and the number of strokes per flash. Laboratory spark chamber results from the literature are used to convert 1-meter channel segments (that are located at a particular known altitude; i.e., air density) to NOx concentration. The resulting raw NOx profiles are discussed.

Koshak, William J.↗

Use of Maps, GTE and UARS-MLS Data in Understanding Tropospheric Processes Critical to Model Development

Two of the major problems in using meteorological models to explain observed tropospheric trace constituent distributions and thereby to understand the global budgets of the tracers are to properly define the vertical layered structure in the free atmosphere, and to understand the contribution of advection processes in generating horizontal inhomogeneities at all scales. We proposed to tackle these problems through the examination of an extensive collection of trace constituent data from research and commercial aircraft in conjunction with meteorological data from the European Center for Medium-Range Weather Forecasts. The physical mechanisms responsible for these advection and layering processes were explored and their implications for theories and models assessed. In addition, we calculated examples of how thin layers (not currently resolved by models) affect the radiative heating/cooling rates. We developed an improved algorithm for trace constituent layer detection, and used it to analyze over 100,000 km of ozone and humidity vertical profiles collected by instruments piggybacked on commercial aircraft. We have also performed radiative heating/cooling calculations on some observed layer structures, and have demonstrated that the sharp edges on humidity layers can have strong local effects of self-stabilization (dry layers) or destabilization (wet layers), and have suggested that clear-air turbulence might result from the latter case through convective instability. We also systematically analyzed the mesoscale variabilities, not only of the trace constituents, but also of the aircraft-measured dynamical variables.

Newell, Reginald E.↗