Search NASA⌕ Search

SEARCH · Search NASA

Results for “reporting and analysis of anomalies”

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.

340 records · Page 19

CFD simulation of coaxial injectors

The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Landrum, D. Brian↗

Analysis of Base Metal Microstructures and Mechanical Properties of a Single-Piece, Spin-Formed Forward Pressure Vessel Bulkhead

In 2012, a pathfinder forward pressure vessel bulkhead (FPVBH) for the Orion Multi-Purpose Crew Vehicle (MPCV) was fabricated using aluminum (Al) alloy 2219. The demonstration article was offered as a replacement for the baseline configuration, which at that time was a multi-piece welded construction using aluminum-lithium (Al-Li) alloy 2195. In 2014, the NASA Engineering and Safety Center (NESC) funded a Phase I effort to explore spin-forming as a manufacturing method to produce a complex-shaped, single-piece FPVBH for the crew module (CM) [1]. The primary focus of the NESC effort was to expand spin-forming technology to the fabrication of an Al-Li 2195 FPVBH. As a result, limited testing and analysis was performed on the original Al 2219 article. In this study, a sufficiently large plate of Al 2219 was not available for the fabrication of a FPVBH. Therefore, the forming blank comprised two plates (from different lots) butt welded together using a single friction stir weld. The welded forming blank was then spin formed into the FPVBH configuration and heat treated to the T62 temper. Mechanical property testing of the spin-formed 2219-T62 FPVBH revealed that the tensile properties were comparable to the Metallic Material Properties Development and Standardization (MMPDS) A-basis allowable for Al 2219-T6 wrought products. However, post-spin forming application of the standard heat treatment produced large, recrystallized grains in the material. The response varied within the two plates; one plate lot exhibited nearly complete recrystallization, while the other lot contained bands of recrystallization. Although no reduction in tensile properties was observed due to these microstructural variations, a more detailed investigation was recommended to better understand the effect of spin forming on the recrystallization mechanism and the impact on other material properties. The current study examined the microstructure and mechanical properties of the fully processed Al 2219-T62 FPVBH with emphasis on comparing differences between the two plates. Testing was limited to the cone region of the FPVBH. The two specific tasks in this study were: a) Characterization of the microstructure and mechanical properties in the two plates used in the FPVBH. b) Isolation of the impact of each step in the thermo-mechanical processing (TMP) sequence of the FPVBH on recrystallization and grain growth. Mechanical property testing sampled regions containing both large, recrystallized and small, unrecrystallized grains. Tensile strengths and elongations were equivalent in both regions and exceeded MMPDS A-basis allowables. Fracture toughness was equivalent in both regions and fracture morphology was typical of aluminum alloys. The fatigue precrack region of the fracture toughness specimens appeared more faceted in the recrystallized plate. This suggests that property testing that is more sensitive to grain size, such as fatigue crack growth rate and stress corrosion cracking, may be warranted. Alcoa reviewed the production records for the two plate lots and confirmed that composition of both plates was within the allowable range for Al 2219 and that plate processing steps were all performed within production limits. However, the Fe content was notably higher in the fully recrystallized plate, resulting in a greater volume fraction of large Fe-bearing inclusions and Al-Cu constituent particles. Alcoa did not provide details regarding plate processing variables for the two plate lots. It was surmised that deformation by hot rolling was the most accurate TMP simulation of commercial spin forming within a laboratory setting. Both hot rolling and spin forming comprise combinations of shear and compressive deformation. The through-thickness strain caused by hot rolling is symmetrical about the mid-plane, as the rollers contact both sides of the material. Shear strains are maximum at the inner and outer surfaces and trend toward zero at the mid-thickness. In contrast, spin forming utilizes a single roller on the outer surface, producing a gradient in shear stress through the thickness. Shear strains in spin-formed material are highest at the outer surface due to direct contact with the roller, and trend toward zero at the inner surface. Consequently, there will undoubtedly be some discrepancies in the mechanical response of hot-rolled and spin-formed products in the through-thickness direction. Interrupted TMP simulations were performed on remnant plate from the FPVBH forming blank in order to isolate the effects of a post-weld anneal, spin-forming deformation/thermal cycling, and solution heat treatment (SHT) on development of the recrystallized microstructures. Analysis confirmed that large, recrystallized grains formed exclusively during SHT and only after a certain deformation level was exceeded. Remnant plate from both lots that received thermal processing only did not exhibit recrystallization but exhibited fully recrystallized microstructures when processed by hot rolling. This indicates that (1) the level of deformation exceeded the critical level for recrystallization in both plates and (2) that the level of deformation during hot rolling was greater than that which accumulated during spin forming of the FPVBH. Comparison of the hot-rolled plate with the spin-formed material showed that the recrystallized grain size was larger in the FPVBH, providing further evidence that deformation levels were lower in the FPVBH. The most plausible explanation for the differing degrees of recrystallization in the two plates is particle stimulated nucleation (PSN) at the Fe-bearing inclusions and large Al-Cu constituent particles. In the plate with higher solute content, the greater volume fraction of these particles reduces the deformation level required to promote recrystallization. During spin forming the critical deformation level was likely exceeded for the higher Fe content plate resulting in a fully recrystallized microstructure. The lower Fe content plate developed a microstructure that exhibited bands of recrystallization due to the combination of a non-uniform distribution of particles and deformation. The critical deformation required for recrystallization was likely only exceeded during spin forming in regions of higher particle content. While there were differences in Fe content between these lots of Al 2219 plate, the composition of both was within defined alloy limits and Alcoa reported no anomalies during plate production. The observed differences in recrystallization after spin forming may reflect the inherent variability between these two lots of plate. The spin-forming vendor, Spincraft, reported that similar variations in recrystallized microstructures have been observed during examination of multiple spin-forming trials of a different 2xxx series aluminum alloy. Hence, tighter chemical specification limits on impurity levels in Al 2219 plate may help ensure uniform, predictable microstructures. Finally, the tensile and fracture toughness properties were not affected by the microstructural variations noted.

single-piece↗

Process Anomaly Detection for Sparsely Labeled Events in Nuclear Power Plants

An essential aspect of online monitoring, subtle anomaly detection increases the detection lead time for equipment failure and enables a nuclear power plant (NPP) to mitigate unexpected partial or full outages, resulting in significant cost saving to the plant. Once an anomaly is detected by plant staff, its cause and severity are investigated. Because the vast majority of anomalies require some level of investigation, including some that require time-consuming examination, before they are passed over to the engineering organization for further analysis, plants are often equipped with tools to assist the staff in performing anomaly detection. Those tools operate as a black box and are often based on statistical methods that establish sensor correlations using preconfigured mathematical models and flag correlation deviations as anomalies. Due to the number of anomalies detected at a given NPP on a daily basis, a significant number of flagged anomalies usually await examination for days or weeks. A primary cause of this backlog is that the methods used by the tools generate many false positives. Though this is usually attributed to oversensitive model settings due to very narrow normal operation bands, it can also be associated with the model development being inadequate for the process being monitored, or with missing model inputs that could have explained misclassified positives. The performance of anomaly detection tools impacts their plant acceptance and utilization, especially when the effort to address false positives generated by the tool depletes the value or cost saved by using that tool. Thus, means to advance anomaly detection performance have been investigated by the Department of Energy’s Light Water Reactor Sustainability program. Previous and ongoing efforts have targeted unsupervised machine-learning (ML) methods, which do not require the labeling of any data fed into the ML model. By contrast, in supervised anomaly detection methods, every data point is labeled as either a normal or abnormal process condition, and the model is trained to replicate the classification process. Supervised methods usually outperform unsupervised methods, due to the added value in differentiating normal from anomalous states of the monitored process. An NPP’s corrective action program requires it to track and document, via a dedicated report, the resolution of any issues that occur within the plant. Once created, each report is reviewed by a plant screening committee, and several classifications and decisions are made. Recently, a collaborating NPP developed an artificial intelligence and ML-based classifier to categorize a condition report (CR) into classes that can serve to label the data as normal or anomalous. Applying CRs as labels represents a semi-supervised use case. Semi-supervised ML assumes that labels exist for some data points (i.e., labeled anomalies, in this case) but not for the rest. In this effort, semi-supervised ML methods were used to fuse data from CRs with anomaly detection methods in order to test the hypothesis that partially labeled anomalies would improve the accuracy of the anomaly detection methods. Specifically, two methods were used. The first is the deep Semi-supervised Anomaly Detection (deep SAD) method, which can handle labels ranging from fully unsupervised to fully supervised cases. The second is a newly designed ML method developed specifically for this effort and referred to as the high-order feature (HOF)-based method. To evaluate these two methods in controlled environments, synthetic data generators were developed and used. The first datasets used a spring-mass-damper (SMD) system simulator commonly found in mechanical engineering references. This was used to create two use cases: a one- and a three-mass system. Anomalies were introduced by changing the spring and damper coefficients while the system was actuated by random forces. The second datasets used the commercial Dymola-Modelica software to build a simplified nuclear reactor model. Anomalies were added in the form of corrupted sensor readings and/or control commands. The deep SAD method was tested using the SMD system, while the HOF method was tested using both datasets. Application of the deep SAD semi-supervised ML method demonstrated that labels can generate increased confidence in detecting true anomalies. This helped increase the number of true positives and decrease the number of false negatives—something that would aid in addressing the backlog of possible anomalies. Application of the HOF method demonstrated that labels can aid in down selecting from a candidate set of features to a more optimal subset in order to better differentiate between normal and anomalous conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Utah FORGE 2439: A Multi-Component Approach to Characterizing In-Situ Stress

Core-based in-situ stress estimation, Triaxial Ultrasonic Velocity (labTUV) data, and Deformation Rate Analysis (DRA) data for Utah FORGE well 16A(78)-32 using triaxial ultrasonic velocity and deformation rate analysis. Report documenting a multi-component approach to characterizing in-situ stress at the U.S. DOE FORGE EGS site: laboratory, modeling and field measurement. Core-based methods for in-situ stress estimation were applied using samples from 5 intervals within the Utah FORGE 16A(78)-32 well. At three of these locations, Triaxial Ultrasonic Velocity (labTUV) tests were performed, resulting in experimentally-determined relationships between wave velocities and stresses. Non-monotonic increase in the velocity-stress relationships are inferred provide evidence of stress history and are therefore used to estimate in-situ stress magnitudes. Additionally, Deformation Rate Analysis (DRA) tests were run on core plugs from various orientations at each of the 5 sampling locations. These, too, provide evidence of stress history based on stress-strain behavior. A novel Weight of Evidence (WoE) method was developed as a means of synthesizing in-situ stress evidence from these two types of tests. Results indicate the minimum horizontal stress gradient ranges from 0.58 psi/ft to 0.69 psi/ft, with 4 of the 5 values between 0.66 psi/ft and 0.69 psi/ft. The vertical stress gradient ranges from 1.05 psi/ft to 1.12 psi/ft, with 4 of the 5 zones given results between 1.09 psi/ft and 1.12 psi/ft. The maximum horizontal stress gradient ranges from 0.98 psi/ft to 1.34 psi/ft, with 4 of the 5 zones falling between 0.98 psi/ft and 1.24 psi/ft. The stress regime thus appears to be on the edge between normal faulting and strike-slip faulting, potentially flipping back and forth between the two regimes due to variability of rock properties, structures such as faults, and/or thermal anomalies.

15 GEOTHERMAL ENERGY↗

Early Spring in Europe: A Result of More Dominant North-Atlantic Southwesterlies?

Abstract A 1999 study reports an advancement of spring in Europe by 0.2 days per year in the 30 years since 1960. Our analysis indicates that this trend results directly from a change in the late-winter surface winds over the eastern North Atlantic: the southwesterly direction became more dominant, and the speed of these southwesterlies increased slightly. Splitting the 52-year NCEP reanalysis dataset into the First Half, FH (1948-1973)), and the Second Half, SH (1974-1999), we analyze the wind direction for the February mean at three sites at 45N: site A at 30W, site B at 20W, and site C at 10W. The incidence (number of years) of the southwesterlies in SH Vs. (FH) at these sites respectively increased in SH as follows: 24(18), 19(12), 14(l 1); whereas the incidence of northeasterlies decreased: 0(2), 1(2), and 1(6). When the February mean wind is southwesterly, the monthly mean sensible heat flux from the ocean at these sites takes zero or slightly negative values, that is, the surface air is warmer than the ocean. Analyzing the scenario in the warm late winter 1990, we observe that the sensible heat flux from the ocean surface in February 1990 shows a "tongue" of negative values extending southwest from southern England to 7N. This indicates that the source of the maritime air advected into Europe lies to the south of the "tongue." Streamline analysis suggests that the Southwestern or southcentral North Atlantic is the source. For February 1990, we find strong, ascending motions over Europe at 700 mb, up to -0.4 Pa/s as monthly averages. Associated with the unstable low-levels of the troposphere are positive rain and cloud anomalies. Thus, positive in situ feedback over land in late winter (when shortwave absorption is not significant) apparently further enhances the surface temperature through an increase in the greenhouse effect due to increased water vapor and cloudiness.

Otterman, J.↗

Cybersecurity for Grid Connected eXtreme Fast Charging (XFC) Station (CyberX) (Final Scientific/Technical Report)

This report summarizes the activities conducted under the DOE VTO funded project DE- EE0008451, where ABB Inc. (ABB), in collaboration with Idaho National Laboratory (INL), APS Global (APS), and XOS Trucks (XOS) pursued the development of a cyber-resilient extreme fast charging (XFC) management system. This project entitled Cybersecurity for Grid Connected eXtreme Fast Charging (XFC) Station (CyberX) focuses on a resilient architecture for smart charging EV Supply Equipment (EVSE) device control and Coordinated Anomaly Detection System (CADS) features that can be added at the charging site depot level to increase cybersecurity. The project was split into two budget periods focused first on developing the threat model and resilient control concepts and second on testing, improving, and validating those developed resilient control algorithms and features with a focus on key vulnerabilities identified during the threat assessment portion of the project. During the first budget period of the CyberX project, the ABB led team focused on activities to identify, model, and quantitatively prioritize high-impact attack scenarios with potential cyber-physical effects while also modeling and developing concepts for a resilient control system that could securely address integration of DERs and other resources with EV charging. Development of the security focused XFC management system (XMS) was accomplished first by offline simulation using a developed XFC station or depot with 480V input level and simulating measurement inputs to monitoring and control systems in concept development. A representative distribution grid model was developed supporting an EV charging site model with BESS and 6 general EV charging models. These EV charging models allowed multiple configurations of charging level, multiple connected protection and measurement devices, and simulation function to show general compromise of EV, BESS, and protection features based on parallel threat analysis. During the second budget period, the EV site and supporting systems model was developed in more detail and converted from offline model to real-time to real-time with EV charging hardware in the loop (HIL). The resilient control architecture developed as concept in the first part of the project was further tested and validated for integration of local energy resources and XFC charging station site equipment while maintaining cybersecure operating principles. The proposed resilient architecture for smart charging and cybersecurity features consists of two main concepts developed and tested within the project. The first concept is an XFC management system (XMS) consisting of a hardware gateway, software platform, and Supervisory Control and Data Acquisition (SCADA) or Distribution Management System integration components. The second concept is a Coordinated Anomaly Detection System (CADS) which forms a primarily software-related subsystem of the total CyberX solution focused on monitoring system measurements, estimation of measurement states, and predicting current at the utility point of interaction based on machine learning for anomaly detection.

33 ADVANCED PROPULSION SYSTEMS↗

Electrically-Active Convection and Tropical Cyclogenesis in the Atlantic and East Pacific

It has been hypothesized that deep, intense convective-scale hot towers may aid the process of tropical cyclogenesis and intensification through dynamic and thermodynamic feedbacks on the larger meso-to-synoptic scale circulation. In this study, we make use of NCEP Reanalysis data and Tropical Rainfall Measurement Mission (TRMM) lightning imaging sensor (LIS), precipitation radar (PR), and microwave imager (TMI) data to investigate the role that widespread and/or intense lightning-producing convection (i.e., electrically-hot towers) present in African easterly waves (AEWs) may play in tropical cyclogenesis over the Atlantic, Caribbean, and East Pacific regions. NCEP Reanalysis 700 hPa meridional winds for the months of June to November for the years 2001-2009 were analyzed for the domain of 5 deg. N-20 deg. N and 130 deg. W-20 deg. E in order to partition individual AEWs into northerly, southerly, trough, and ridge phases. Subsequently, information from National Hurricane Center (NHC) storm reports was used to divide the waves into developing and non-developing waves. In addition, information from the NHC reports was used to further divide the developing waves into those waves that spawned storms that only developed to tropical storm strength and those that spawned storms that reached hurricane strength. The developing waves were also divided by the region in which they developed. To assess the evolution of convection associated with the AEWs as they propagated across our analysis domain, the full 130 deg. W-20 deg. E domain was divided into five longitude bands, and waves were analyzed for each band. To help determine the gross nature of the smaller convective scale, composites were created of all developing and non-developing waves as a function of AEW wave phase over the full analysis domain and each longitude band by compositing TRMM PR, TMI, LIS, and IR brightness temperature data extracted from the NASA global-merged infrared brightness temperature dataset. Finally, similar composites were created using various NCEP variables to assess the nature of the larger scale environment and circulation. Results suggest a clear distinction between developing and non-developing waves as developing waves near their region of development in terms of the intensity of convection (indicated by lightning flash rate), coverage of cold cloudiness (indicated by the percentage of a 2.5 deg. by 2.5 deg. box covered by IR brightness temperatures less than 240 deg. K and 210 deg. K), and large-scale variables, such as midlevel moisture and upper-level upward motion. For example, waves that developed in the East Pacific longitude band (i.e., 130 deg. W 95 deg. W) were observed in that band to have a flash rate of 56.4 flashes per day, a coverage by brightness temperatures less than 240 deg. K equal to 15.9%, a coverage by brightness temperatures less than 210 deg. K equal to 2.2%, a 700-hPa specific humidity anomaly of 0.4 g per kilogram, and a 300-hPa omega value of -0.04 Pascals per second in the trough phase compared to the non-developing wave trough values of 22.1 flashes per day, a coverage by brightness temperatures less than 240 deg. K equal to 8.1%, a coverage by brightness temperatures less than 210 deg. K equal to 0.9%, a 700-hPa specific humidity anomaly of -0.3 g per kilogram, and a 300-hPa omega value of -0.01 Pascals per second. Further analysis is being conducted to determine if the aforementioned behavior is observed for developing waves farther from their region of development and to determine any significant differences between waves that spawned storms that reached tropical storm strength and those that spawned storms that reached hurricane strength.

Leppert, Kenneth D., II↗

Supporting Crew Autonomy in Deep Space Exploration: Preliminary Onboard Capability Requirements and Proposed Research Questions. Technical Report of the Autonomous Crew Operations Technical Interchange Meeting

Communication delays are a critical challenge posed by long duration deep space exploration. Space missions historically have relied on an ever-present Mission Control Center (MCC) to direct operations in near real-time. As unanticipated anomalies that defeat fault detection and resolution systems do arise, the lack of real-time communication will significantly weaken what the MCC support represents: a reliable safety net for the flight crew through its deep and diverse areas of expertise and investigative resources. As a consequence, future space vehicles and habitats need to be equipped with capabilities to support the flight crew to operate with little or no ground support. Considerations must be given to vehicle and mission designs that will fortify the traditionally ground-centered safety net and forge new support systems, when communication delays exist. In August 2018, NASA’s Human Research Program, through its Human Factors and Behavioral Performance Element, convened a Technical Interchange Meeting (TIM) on Autonomous Crew Operations at NASA Ames Research Center. The goal of the meeting was to gather input from NASA centers, industry, academia, and branches of the Department of Defense (DoD) to address how intelligent technologies can be applied to augment onboard capabilities to support crew anomaly response. The TIM featured 24 presentations by 29 speakers and hosted a total of 59 attendees, including 43 from 5 NASA centers (Ames, Johnson, Langley, Marshall, and Jet Propulsion Lab) and 4 from the DoD (3 from Army Research Lab and 1 from Naval Postgraduate School), with remaining attendees from academia (e.g., UC Davis, CMU) and industry (e.g., IBM, Siemens). Discussions were centered around three themes: standards and guidelines, lessons learned in analog environments, and technologies. To help provide a framework for discussion, a concept matrix describing anomaly response processes was created prior to the TIM (Figure 1, page 6). The matrix captures the steps involved (monitoring and detection, diagnosis, solution development and evaluation, solution implementation and verification, resolution documentation) as well as the resources and capabilities required to support these steps (data, knowledge, analysis, synthesis, resource management). A wallpaper size printout of the matrix was utilized at the TIM to solicit attendee inputs along the three themes; the activity garnered 108 submissions of ideas. Overall, what emerged from TIM discussions was a picture of mismatch between crew anomaly response needs and support that can be provided by existing intelligent technologies. The needs are broad, spanning multiple steps and processes/resources, with many of which lacking support from existing technologies, such as knowledge management throughout the steps of problem solving (especially in resolution documentation) and manpower management. The solutions provided by existing intelligent technologies are specific to the steps/processes that they are designed to support and constrained to solving only problems similar to those that have occurred before. What is lacking from technologies is typically made up by humans, specifically their complex critical thinking, creative problem solving, and domain expertise. In the end, the TIM highlighted the pressing need to support responses to onboard anomalies during autonomous crew operations, particularly those that have eluded the system tests, inspection, and other assurance processes. Such anomalies can potentially threaten crew and vehicle safety, as well as significantly impact overall operations with additional workload. These fairly rare events are difficult to anticipate and prepare for, given the state-of-the-art in intelligent technologies. This is true even for anomalies that stem from “unknown knowns”—cases in which there is sufficient external information to characterize the problem but the overall pattern fails to be recognized by the problem solver, or in which the internal knowledge needed to solve a problem is held tacitly and potentially accessible by the problem solver but not articulated. It follows that the ability to tackle anomalies lies not only with the availability of relevant information and knowledge but also their accessibility in times of need. To that end, we propose research questions along the following three broad themes: • How intelligent technologies can help make relevant knowledge and information available? • How intelligent technologies can help make relevant knowledge and information accessible? • How intelligent technologies can help support the crew operating as a team in anomaly response processes?

autonomous crew operations↗

Evaluation of a Surface Spill Operation to Return Adult Steelhead Overshoots Downstream of McNary Dam

This study, funded by the U.S. Army Corps of Engineers (USACE), was conducted by the Pacific Northwest National Laboratory to evaluate the efficacy of operating one Top Spill Weir (TSW) at McNary Dam outside the normal TSW operation dates for juvenile salmon passage. Of interest is whether the TSW is an effective downstream passage route for adult steelhead overshoots. Overshoots are fish which, having passed upstream at McNary Dam, must pass downstream to return to their natal stream to spawn. This report covers Fall (2019) and Spring (2020) study periods. The study design arranged the available 24 hours of TSW spill per week into weekly blocks with sub-blocks differentiated into day and night operations with TSW discharge periods of 4- and 8-hours duration. Hydroacoustic techniques were used to sample adult fish passage at the TSW and at turbine units 1 and 10. The experimental design contrasted TSW spill periods of differing duration and at different times of the day. The small number of fish detected passing the TSW, and the smaller number of fish detected passing the turbine units, were best suited to an ad hoc, exploratory approach to evaluating the effect of TSW spill. It is worthwhile to note that the operations data obtained for the fall study period had a greater than expected number of gaps and apparent anomalies that we believe were a result of how the data were aggregated. These problems are not particularly problematic for the present study, because the available data still provide a good indication of whether the TSW was operating at each point in time. TSW flows in both fall and spring data sets were able to be cleaned up using the established relationship between forebay elevation and TSW discharge rate. Additional cleanup of operations data would be needed, however, if more quantitative evaluation of dam-wide flow and passage relationships were needed. This study contrasted TSW spill periods of differing duration and at different times of the day. The small number of fish detected passing the TSW, and the smaller number of fish detected passing the turbine units were best suited to an ad hoc, exploratory approach to evaluating the effect of TSW spill. Spring adult steelhead passage numbers estimated using hydroacoustics were notably lower than during the fall study period, and that was consistent with our analysis of PIT tagged fish likely to be in the vicinity during each study period. Detections of fish in BlueView sampling areas upstream of the TSW and powerhouse were not correlated with detections of fish passing hydroacoustic sampling areas, which suggests that fish approaching the face of the dam can move around the forebay before passing. Other fish detected in the forebay in large numbers, such as shad, were able to be filtered out of steelhead passage counts and did not appear to be influencing hydroacoustic passage rate estimates. A pulse of passage at the initial TSW opening was weakly evident, but trends across 4- and 8-hour operational periods did not show a distinct decline in passage over time as TSW operation continued. Our findings do not indicate a reason to choose one 8-hour period over two 4-hour periods, or vice versa. This suggests that the duration of spill periods can be chosen based on operational or other considerations. Passage rates were consistently higher during the daytime TSW discharge periods, relative to nighttime TSW discharge periods. The experimental design of the current study used start times near dawn for day periods and near dusk for night periods.

59 BASIC BIOLOGICAL SCIENCES↗

Phasor-Measurement-Unit-Based Data Analytics Using Digital Twin and PhasorAnalytics Software

A major objective of this project was to apply GE’s commercial machine learning and data analytics toolsets to large-scale, real-world, anonymized Phasor Measurement Unit (PMU) datasets in order to extract signatures, correlated and/or causal factors, and precursor patterns associated with significant power system phenomena. The project had a particular emphasis on extraction of insights relevant to asset health monitoring, real-time load modeling and cybersecurity monitoring. Additionally, the team was directed to undertake a comprehensive data quality analysis for the provided datasets and encouraged to estimate the ‘machine-learning readiness’ of the datasets by documenting any major obstacles to the application of commercial machine learning algorithms. To accomplish the aforementioned objectives, the project team’s work centered around the identification of key event signatures and application of the identified event signatures for event detection and event classification. The industry-validated, semi-supervised machine learning strategy employed for event signature identification involved several major tasks, including data-preprocessing, generation of an overabundance of features, normal data identification, normality modeling, and event signature identification through a methodical, quantitative ranking of features in order of relevance to each studied event type. Throughout the project, data quality issues and mitigation techniques were investigated. In this report, insights are provided regarding the readiness of the provided synchrophasor datasets for application of machine learning and data analytics. The methodologies employed for this technical strategy are summarized in this report. With regards to data preprocessing and feature generation, the provided Training and Test Datasets were ingested into GE’s big data environment. Subsequently, the team applied bad data cleansing and data imputation scripts, event detection scripts, and application programming interfaces (APIs) to the datasets for convenient data access. The project team completed development and validation of dozens of physics-based, statistics-based and transformation-based feature functions used for the extraction of over 60 synchrophasor features. Using a new parallel feature generation technology developed on this project, over 60 features have been rapidly generated for the full two years’ worth of Training and Test Dataset data associated with both the Eastern and Western interconnects. Even accommodating for temporal down-sampling inherent to the feature extraction procedure, this parallel feature generation activity resulted in a massive feature set with a storage requirement approximately equal to that of the raw training dataset itself. With regards to normal data identification and normality modeling, a normality model was built using the feature data extracted from the Training Dataset and iteratively refined subsequent to incremental adjustments and expansions of the Training Dataset feature data. With respect to event characterization and signature identification, an event signature identification pipeline was developed and used in conjunction with the normality model to identify over 15 event signatures for key event categories within the Training Dataset. The identified event signatures were used to characterize hundreds of key events in terms of relative severity, duration, and location of the event. An investigation was undertaken to identify correlated and causal factors involved in transformer events. A separate investigation into temporal trends in ring-down analysis results was undertaken to determine possible associations between system dynamics and various other factors such as loading, season or year. To validate the identified event signatures, additional work was undertaken to develop signature-based anomaly detection and classification tools suitable for convenient application to the synchrophasor datasets. The anomaly detection and classification tools, suitable for online application, were then applied to the entirety of the Eastern Interconnect Training and Test Datasets. Performance of the event detection and classification tools was evaluated upon receipt of the Test Dataset event logs (i.e., the labels for events contained in the Test Dataset), and promising results were obtained despite several challenges (documented herein) associated with application of supervised or semi-supervised machine learning methods to large-scale, anonymized datasets. Finally, the detection and classification tools were used to detect, classify, and characterize thousands of new events not included in the original event logs provided by the DOE within both the Training and Test Datasets.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Play Fairway Analysis: Structurally Controlled Geothermal Systems in the Eastern Great Basin Extensional Regime, Utah

A research team with membership from the University of Utah/Energy & Geoscience Institute, the University of Utah/Dept. of Geology & Geophysics, and the Utah Geological Survey, undertook a play fairway analysis (PFA) for geothermal resources in the Eastern Great Basin (EGB) extensional tectonic regime of western Utah. This is a high-priority region for geothermal exploration because active Basin and Range (B&R) extension with volcanism having a N-S strike is superimposed upon pre-existing E-W belts of plutonic rocks and large-scale structural lineaments. Cumulative heat flow along the N-S strike of the state totals approximately 5 GWt above background stable interior. Three electricity producing power plants currently exist with substantial potential for increase. Succinctly, our PFA approach aims to resolve potential sources of heat and permeability in the region, which are the two principal criteria for establishing a geothermal resource. An initial Phase 1 was carried out using only existing geoscientific data in the area. Criteria selected for focusing heat potential include direct heat flow measurements in boreholes, magnetotelluric (MT) low resistivity anomalies, fluid/gas geochemistry, and proximity to recent volcanic eruptions. Permeability is established through geological structures (fault density, critically stressed areas, seismicity, gravity), and MT low resistivity anomalies. In Phase II of this PFA project, additional geological, geophysical and geochemical data were acquired and analysis carried out primarily over promising composite common risk (CCR) areas initially identified in Phase I in order to focus prospectivity and prepare for drilling recommendations. These prospects are near the Twin Peaks rhyolite field, high heat flow areas north of the producing Cove Fort system, and geophysical structure beneath the Crater Knoll area off the northeast flank of the Mineral Mountains. Additional data included MT site fill-in, structural mapping and analysis using high-resolution imagery, gravity and on-ground mapping, Nodal 3C passive seismic collection, and passive 3He surveying. Heat source and permeability potential are again expressed in terms of their individual common risk segment (CRS) maps, with a color scheme using green for most favorable (low risk) and red for least favorable (high risk). Diverse data types are united through the technique of probability kriging, which establishes prospectivity thresholds for each data type and then computes probability of exceeding that threshold over the PFA area. Modified CRS and CCRS maps are compared to those of Phase I to highlight tighter prospectivity focus. In doing so, the prospectivity threshold for heat was increased significantly to narrow the targeting. In the final Phase III of this project where a recommended deep thermal gradient hole was sited, additional geophysical, geochemical and geological field collection and analysis was carried out to refine drill hole targeting. This includes prospect-scale MT, gravity, structure, passive seismic deployment (Cove Fort area), and detailed 3He isotope profiling. It was the recommendation of the DOE Technical Monitoring Team (TMT) that one or more holes be sited in the north-ern Cove Fort area where legacy TG gradient holes showed high cumulative heat flow. These were to be of moderate depth, 2000-3000 feet, to reach the geothermal fluid table expected to start in excess of 1000 feet depth. The drilling organization stipulated by the DOE/GTO was that of the USGS Research Drilling Program (RDP) centered in Las Vegas, NV. A detailed well plan, appended to this report, was developed principally by Dr. Ben Barker consulting to University of Utah, Dr. Steve Pye on the DOE TMT, Mr. Steven Crawford of the USGS-RDP, and the project PI Phil Wannamaker. However, temperature and possible H 2 S at the systems lead to cancellation of the drilling last-minute as this appeared outside the experience base of the USGS-RDP. We hope to have the opportunity to revisit the test drilling and expand the Play Fairway Analysis of this region at some point in the future.

15 GEOTHERMAL ENERGY↗

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗

Debris/ice/TPS assessment and integrated photographic analysis for Shuttle Mission STS-62

A pre-launch debris inspection of the pad and Shuttle vehicle was conducted on 2 March 1994. The detailed walkdown of Launch Pad 39B and MLP-1 also included the primary flight elements OV-102 Columbia (16th flight), ET-62 (LWT 55), and BI-064 SRB's. There were no significant facility or vehicle anomalies. After the launch on March 4th, a debris inspection of Pad 39B was performed. Damage to the pad overall was minimal. On-orbit photographs taken by the flight crew and two films from the ET/ORB umbilical cameras of the External Tank after separation from the Orbiter revealed no major damage or lost flight hardware that would have been a safety of flight concern. Orbiter performance on final approach appeared normal. Infrared imagery of landing gear deployment showed the loss of thermal barrier from the nose gear wheel well. The missing thermal barrier material was not recovered. The Solid Rocket Boosters were inspected at Hanger AF after retrieval. Both frustums had a combined total of 44 MSA-2 debonds over fasteners. Significant amounts of BTA had been applied to closeouts on the RH frustum, forward skirt, and aft skirt. Hypalon paint was blistered/missing over the areas were the BTA had been applied. The underlying BTA was not sooted (IFA STS-62-B-1). Investigation of this condition has concluded there was insufficient heat rates to cause blistering of the Hypalon until late in the ascent phase. A post landing inspection of OV-102 was conducted after the landing at KSC. The Orbiter TPS sustained a total of 97 hits, of which 16 had a major dimension of 1 inch or larger. The Orbiter lower surface had a total of 36 hits, of which 7 had a major dimension of 1 inch or larger. Based on these numbers and comparison to statistics from previous missions of similar configuration, both the total number of debris hits and the number of hits 1 inch or larger was less than average. Six thermal barriers, total size approximately 36 in. x 3 in. x 1.5 in., and one corner tile piece 4 in. x 4 in., were missing from the nose landing gear doors. Runway infrared cameras recorded these objects falling from the Orbiter when the nose landing gear doors were opened on final approach. A search of the area under the flight path failed to recover the missing thermal batteries, which may be submerged in the large body of water south of the runway. The cause for the loss of the material was attributed to the way the substrate was prepared for bonding at the vendor. Orbiter post landing microchemical sample results revealed a variety of residuals in the Orbiter window samples. These were attributed to SRB BSM exhaust, Orbiter TPS, window protective covers and processing, natural landing site products, and paints/primers from various sources. The residual sampling data do not indicate a single source of damaging debris and all of these materials have been documented previously in analyses for post landing sample reports. Data from residual sampling also showed no trends when compared to previous mission data. A total of ten Post Launch Anomalies, including one IFA, were observed during the STS-62 mission assessment.

Katnik, Gregory N.↗

Field Validation of MVA Technology for Offshore CCS: Novel Ultra-High-Resolution 3D Marine Seismic Technology (P-Cable) (Final Report)

The objectives of the proposed study were to deploy and validate a specific monitoring technology, high-resolution 3D marine seismic (HR3D), appropriate for large-demonstration and commercial-scale offshore CCS sites. The project accomplished successful acquisition two HR3D seismic surveys. The first HR3D dataset was over the offshore injection site of the Tomakomai, Japan integrated pilot CCS project, which at the time of survey acquisition was actively injecting CO 2 . The first survey also represented a successful international collaboration between the DOE NETL program and Japan’s national CCS program and was the first successful acquisition and use of HR3D over an active CO 2 injection site (Meckel, Feng et al. 2019). The Tomakomai HR3D survey successfully tested a novel 4-streamer HR3D system array in which, for the first time, no cross-cable (aka “P-Cable”) was utilized and only four GeoEel streamers were used instead of the standard 12-streamer configuration. Consequently, this was not, strictly speaking, a deployment of the “P-Cable” system of (Planke and Berndt 2004) but rather a modified version, thereof, and it is the first known demonstration of the modified system configuration. One very positive outcome from the Japanese collaboration earlier in the project was the ability to learn from the Japanese how they used tail buoys with GPS to determine the position of the seismic source and receivers in time and space. Based on that experience, GCCC designed and built six GPS receivers that could be used to position the streamer receivers and the seismic source via tail buoys. A fundamental advance that was made on the original design, was the ability to directly power the tail buoy GPS units and transfer data through the streamers (i.e., vs. the batteries used at Tomakomai). The bulkiness of the GPS batteries caused drag and episodic surging of the buoys, which affected data quality by lifting up the tail end of the streamers so the receivers were not at the same depth. The units were tested onshore for accuracy and functionality, and the design was subsequently and successfully tested in marine acquisition mode during the SLP survey acquisition. The marine acquisition test and survey satisfied Subtasks 2.2.2, Novel Positioning Technology Selection and Subtask 2.2.3, Novel Positioning Technology Deployment. Results of the novel positioning technology selection (Subtask 2.2.2) were considered successful and will be incorporated in future HR3D seismic acquisition projects to reduce costs, improve deployment safety at sea, and integrate both seismic and data recording via a single data transfer through the streamers to the recording system. The project also established a permitting process through NETL NEPA compliance, which included an Environmental Assessment in a marine setting and is required for conducting these types of surveys using Federal funding. The permitting process charted a “boilerplate,” which can allow future surveys related to other funded projects to move forward more expeditiously. Future improvements that could be considered are more robust seals on the GPS module and stronger materials (especially joints) on tail buoy fabrication. These would increase fixed costs, but would be advisable and probably more economic long-term if multiple HR3D surveys are planned. Project Accomplishments include: • Pre-survey Sensitivity Study • Marine geochemistry methods and data analysis • Successful HR3D seismic dataset acquired @ Tomakomai active CO 2 injection marine site • Developed advanced seismic processing techniques • No NRMS anomalies detected in overburden; Demonstration of containment • Repeatability study • Second survey collected @ San Luis Pass, TX • 4D application using positioning techniques developed in the project for monitoring were successful

3D seismic GPS positioning↗

Soil Moisture Active Passive (SMAP) Project Assessment Report for Version 4 of the L4_SM Data Product

This report provides an assessment of Version 4 of the SMAP Level 4 Surface and Root Zone Soil Moisture (L4_SM) product, released on 14 June 2018. The assessment includes comparisons of L4_SM soil moisture and temperature estimates with in situ measurements from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product, including observation-minus-forecast (O-F) brightness temperature residuals and soil moisture analysis increments.Together, the core validation site comparisons and the statistics of the assimilation diagnostics areconsidered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to upscaling errors from the point-scale to the grid-cell scale of the data product.The Version 4 L4_SM product benefits from an improved land surface modeling system and from retrospective surface meteorological forcing data that are as consistent as possible with the present-day datain terms of their climatology. Specifically, the model changes include revised parameters and parameterizations for (i) the surface energy balance, (ii) recharge from below of the model's surface excess reservoir, and (iii) the snow depletion curve. Updated ancillary inputs include improved datasets for landcover, topography, and vegetation height. The Version 4 algorithm further includes a revised approach to precipitation corrections that improves the precipitation climatology in Africa and the high-latitudes. Moreover, for system calibration the model is forced retrospectively with MERRA-2 reanalysis data, which are more consistent with the near-real time GEOS forward processing (FP) data used during the SMAP period than the retrospective GEOS data that were available for previous L4_SM versions. An analysis of the time-average surface and root zone soil moisture shows that the global pattern ofarid and humid regions is captured by the Version 4 L4_SM estimates. Owing to the changes in the landsurface modeling system, surface soil moisture is typically drier by several volumetric percent in Version 4 compared to Version 3, whereas root zone soil moisture is wetter in Version 4 in some regions and drierin others. Because of these climatological differences, the Version 3 and Version 4 products should not be combined into a single dataset for use in applications.Results from the core validation site comparisons indicate that Version 4 of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the RMSE after removal of the long-term mean difference (ubRMSE). The overall ubRMSE of the 3-hourly L4_SM dataat the 9 km scale is 0.039 m3 m-3 for surface soil moisture and 0.029 m3 m-3 for root zone soil moisture,below the 0.04 m3 m-3 requirement. The L4_SM estimates are an improvement over estimates from a model-only Nature Run version 7.2 (NRv7.2), which demonstrates the beneficial impact of the SMAP brightness temperature data. Overall, L4_SM surface and root zone soil moisture estimates are more skillful than NRv7.2 estimates, with statistically significant improvements at the 5% level for surface soil moisture R and anomaly R values. Results from comparisons of the L4_SM product to i

Reichle, Rolf H.↗

Global Geochemical Variation on the Lunar Surface: A Three-Element Approach

We present a method for displaying the relative abundances of three important elements (Th, Fe, and Ti) on the same map projection of the lunar surface. Using Th-, Fe-, and Ti-elemental abundances from orbital geochemical data and assigning each element a primary color, a false-color map of the lunar surface was created. This approach is similar to the ternary diagram approach presented by Davis and Spudis with some important differences, discussed later. For the present maps, Th abundances were measured by the Lunar Prospector (LP) Gamma-Ray Spectrometer(GRS).The new LPGRS low-altitude dataset was used in this analysis. Iron and Ti weight percentages were based on Clementine spectral reflectance data smoothed to the LP low altitude footprint. This method of presentation was designed to aid in the location and recognition of three principal lunar compositions: ferroan anorthosite (FAN), mare basalts (MB), and the Mg suite/ KREEP-rich rocks on the lunar surface, with special emphasis on the highlands and specific impact basins. In addition to the recognition of these endmember rock compositions, this method is an attempt to examine the relationship between elemental compositions that do not conform readily to previously accepted or observed endmember rocks in various specific regions of interest, including eastern highlands regions centered on 150 deg longitude, and a northern highlands Th-rich region observed. The LP low-altitude data has full width at half-maximum spatial resolution of about 40 km. The Clementine spectral reflectance datasets were adapted using an equal-area, gaussian smoothing routine to this footprint. In addition, these datasets, reported in weight percent of FeO and of Ti02, were adjusted to Fe and Ti weight percentages. Each dataset was then assigned one of the three primary colors: blue for Th, red for Fe, and green for Ti. For each element, the data range was normalized to represent the ratio of each point to the maximum in the dataset. (To view the color table, go to http://cass.jsc.nasa.gov/meetings/moon99/pdf/8033.pdf.) The full range of lunar longitudes is represented, but due to the lack of coverage of the Clementine data for latitudes > 70 deg and <-70 deg, the data for these regions is excluded. The differences between this approach and the ternary diagram approach of Davis and Spudis eliminate some of the uncertainty and ambiguity of the ternary diagram approach. Rather than using a ratio of Th to Ti normalized to CI chondritic ratios, and a ternary diagram with ternary apexes located at specific endmember compositional values, elemental compositions were used independently, eliminating the errors resulting from dividing numbers that can have high uncertainties, especially at low concentration. The three elements used in this method of presentation were chosen for several reasons. One reason for the inclusion of Th in this study is that it is an accurate indicator of KREEP. Iron and Ti concentrations are both low in highland regolith, causing any small fluctuations in Th to stand out very well. In addition, Fe and Ti are good compositional indicators of different mare basalts. Mixed with red for Fe, the green for Ti produces a yellow signal in high-Ti basalts. While universally high in Fe relative to the surrounding highlands, mare basalts have a diverse range of Ti values, making Ti concentration a valuable asset to the classification and identification of different basalt types. Finally, an important constraint in element selection is the availability of the global data, both from LP and Clementine results. Data for Th, Fe, and Ti are among the highest quality of existing lunar remote-sensing data. In addition, LP data for Fe and Ti will become available, enabling these data to be incorporated into the analysis. Using upper-limit values for end member rock compositions calculated from Korotev et al., attempts were made to locate the different endmember compositions of terranes on this diagram. Most strikingly, ferroan anorthosite (Th < and = 0.37 micro g/g; Fe (wt%)< and =2.29; Ti (wt%) < and = 0.22), which should appear as an almost black, reddish color, does not appear on the diagram at any noticeable frequency. Based on this analysis, the suggestion of extensive FAN regions on the lunar surface is not strong, especially at the presently accepted values for Fe and Th. However, to make sure this effect is not due to systematic errors, a thorough investigation of the precision, accuracy, and uncertainties of the Fe, Ti, and Th abundances needs to be carried out, especially at low concentrations. A particular region of interest is an area of high Th concentrations relative to Fe and Ti content north and east of Humboldtianum Crater. First observed by Lawrence et al., this region does not coincide with any visible impact structure and comprises one of the closest approximations to pure blue (high Th, very low Ti and Fe) on the lunar surface. Such an elemental composition does not lend itself readily to classification, and presents something of an anomaly. More detailed analysis of this region is needed to understand its structure and origin. There seems to be a longitudinal asymmetry in the Th concentrations of the highlands regolith. High-Th, low-Ti, and Fe regions are located between 135 deg and 180 deg longitude and between -30 deg and +30 deg latitude. While the Th levels are not high enough to attract attention in a single elemental display, the variation in the abundance of Th relative to Fe and Ti abundances can be clearly seen. The composition that these data suggest is not well represented in the sample return suite. In addition, these regions were largely missed by the Apollo orbital ground tracks, which only covered the outer edge of the areas of interest. The LP orbital Th data represent the first information about the Th concentrations in these regions of the highlands. Additional information contained in original.

Thomsen, D. R.↗