Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “missing information”

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 883 records · Page 49

Review of Quantitative Monitoring Methodologies for Emissions Verification and Accounting for Carbon Dioxide Capture and Storage for California’s Greenhouse Gas Cap-and-Trade and Low-Carbon Fuel Standard Programs

The Cap-and-Trade and Low Carbon Fuel Standard (LCFS) programs being administered by the California Air Resources Board (CARB) include Carbon Dioxide Capture and Storage (CCS) as a potential means to reduce greenhouse gas (GHG) emissions. However, there is currently no universal standard approach that quantifies GHG emissions reductions for CCS and that is suitable for the quantitative needs of the Cap-and-Trade and LCFS programs. CCS involves emissions related to the capture (e.g., arising from increased energy needed to separate carbon dioxide (CO 2 ) from a flue gas and compress it for transport), transport (e.g., by pipeline), and storage of CO 2 (e.g., due to leakage to the atmosphere from geologic CO 2 storage sites). In this project, we reviewed and compared monitoring, verification, and accounting (MVA) protocols for CCS from around the world by focusing on protocols specific to the geologic storage part of CCS. In addition to presenting the review of these protocols, we highlight in this report those storage-related MVA protocols that we believe are particularly appropriate for CCS in California. We find that none of the existing protocols is completely appropriate for California, but various elements of all of them could be adopted and/or augmented to develop a rigorous, defensible, and practical surface leakage MVA protocol for California. The key features of a suitable surface leakage MVA plan for California are that it: (1) informs and validates the leakage risk assessment, (2) specifies use of the most effective monitoring strategies while still being flexible enough to accommodate special or site-specific conditions, (3) quantifies stored CO 2 , and (4) offers defensible estimates of uncertainty in monitored properties. California’s surface leakage MVA protocol needs to be applicable to the main CO 2 storage opportunities (in California and in other states with entities participating in California’s Cap-and-Trade or LCFS programs), specifically CO 2 -enhanced oil recovery (CO 2 -EOR), CO 2 injection into depleted gas reservoirs (with or without CO 2 -enhanced gas recovery (CO 2 -EGR)), as well as deep saline storage. Regarding the elements of an effective surface leakage MVA protocol, our recommendations for California are that: (1) both CO 2 and methane (CH 4 ) surface leakage should be monitored, especially for enhanced recovery scenarios, (2) emissions from all sources not directly related to injection and geologic storage (e.g., from capture, or pipeline transport) should be monitored and reported under a plan separate from the surface leakage MVA plan that is included as another component of the quantification methodology (QM), (3) the primary objective of the surface leakage MVA plan should be to quantify surface leakage of CO 2 and CH 4 and its uncertainty, with consideration of best-practices and state-of-the-art approaches to monitoring including attribution assessment, (4) effort should be made to monitor CO 2 storage and migration in the subsurface to anticipate future surface leakage monitoring needs, (5) detailed descriptions of specific monitoring technologies and approaches should be provided in the MVA plan, (6) the main purpose of the CO 2 injection project (CO 2 -EOR, CO 2 -EGR, or pure geologic carbon sequestration (GCS)) needs to be stated up front, (7) approaches to dealing with missing data and quantifying uncertainty need to be described, and (8) post-injection monitoring should go on for a period consistent with or longer than that prescribed by the U.S. EPA.

54 ENVIRONMENTAL SCIENCES↗

Digital Particle Image Velocimetry (DPIV) Used for Space-Time Correlations in Nozzle Flow

An optical measurement technique known as Digital Particle Image Velocimetry (DPIV) was used previously to characterize the first- and second-order statistical properties of both cold and hot jet flows from externally mixed nozzles in NASA Glenn Research Center's Nozzle Acoustic Test Rig. In this technique, an electronic camera records particles entrained in a flow as a laser light sheet is pulsed at two instances in time. Correlation processing of the recorded particle image pairs yields the two-component velocity field across the imaged plane of the flow. The information acquired using DPIV is being used to improve our understanding of the decay of turbulence in jet flows-a critical element for understanding the acoustic properties of the flow. Recently, two independent DPIV systems were installed in Glenn's Small Hot Jet Acoustic Rig, enabling multiplane correlations in time and space. The data were collected over a range of different Mach numbers and temperature ratios. DPIV system 1 was fixed to a large traverse rig, and DPIV system 2 was mounted on a small traverse system mounted on the large traverse frame. The light sheets from the two DPIV systems were aligned to lie in the same axial plane, with DPIV system 2 being independently traversed downstream along the flow direction. For each measurement condition, the DPIV systems were started at a fully overlapping orientation. A polarization separation technique was used to avoid cross-talk between the two systems. Then, the DPIV systems fields were shifted axially apart, in successively increasing steps. The downstream DPIV system 2 was triggered at a short time delay after the upstream DPIV system 1, where the time delay was proportional to the convective flow velocity in the shear layer of the jet flow and the axial separation of the two DPIV systems. The acquired data were processed to obtain the instantaneous velocity vector maps over a range of time delays and spatial separations. The velocity fields from the different DPIV systems were then cross-correlated to determine the degree of correlation remaining in the flow as the downstream convection distance was increased. The new data provide Lagrangian measurements of the convective turbulent structures in the shear layer of an exhaust nozzle. These measurements, obtained in both cold and hot flows, will be used to validate and correct models for space-time velocity correlations-long a missing key to predicting jet noise.

Wernet, Mark P.↗

Dynamics of A Vibration Isolation System Including Inertia of the Human Body

Using an exercise device in a spacecraft is liable to transmit an unacceptable amount of vibration to that vehicle. This is commonly mitigated by a Vibration Isolation System (VIS), whose dynamics must be analyzed to confirm that the oscillatory forces on the spacecraft remain within allowed range, both from a structural and microgravity perspective (see, e.g., [1]). When modeling a VIS for countermeasures devices, one common approach is to record forces and moments applied on the floor while exercising, and then drive the VIS simulation by applying these recorded loads to the exercise platform part of the VIS mechanical model. This approach misses the fact that when exercising on a moving platform, the force and moment on it will differ from that on the stationary floor due to inertial effects involving the human body. For example, standing up on a platform as it gives under the subject’s feet reduces the foot force on it, and such inertial effects are especially complex for rotational motion. In principle, one could model both the motion of the human body and dynamics of the VIS mechanism in a single combined simulation, e.g., employing a tool such as the commonly used biomechanical simulation OpenSim. Here, the joints of the human body would be driven kinematically along prescribed exercise trajectories while the dynamics engine computed the response of the VIS degrees of freedom. However, mechanism designers and biomechanics experts have their own established tools, making it very desirable to have a way of decoupling the biomechanics from the VIS modeling, simulation, and analyses. We have derived a set of equations that rigorously accomplishes this goal, and have implemented them as an interface function that provides an alternative driving mechanism for an existing force-based VIS analysis simulation. When enabled, the simulated human/VIS system dynamics is now driven by this function, instead of the recorded force methodology described above. The function is designed to accept input from a data file containing the required time-stamped human motion and inertia terms corresponding to the specific exercise in question. This data file is generated by an OpenSim plugin written for that purpose. The existing VIS analytical simulation is developed using NASA’s Trick Simulation Environment [2], as well as its MBDyn multibody dynamics [3] package. The presentation will provide a detailed overview of the mathematical formulation, assumptions, plugin implementation, software interfaces, and results for a sample set of representative exercises. The results from this work aim to better inform VIS design efforts, as well as countermeasure device/protocol designs with respect to exercise type and frequency effects on vehicle structural and microgravity restrictions.

Countermeasures↗

Oak Ridge National Laboratory Technical Input for the Nuclear Regulatory Commission Review of the 2017 Edition of ASME Section III, Division 5, ‘High Temperature Reactors’

To assist the Nuclear Regulatory Commission in its decision making on endorsement of the American Society for Mechanical Engineers Boiler and Pressure Vessel Code Section III, Division 5 (2017 Edition) for development of advanced non-light water reactors, the following Division 5 portions were reviewed: Article HBB-2000 Material; Article HCB-2000 Material; Article HGB-2000 Material; Mandatory Appendix HBB-I-14 Tables and Figures; and, Nonmandatory Appendix HBB-U Guidelines for Restricted Material Specifications to Improve Performance in Certain Service Applications. In addition to the 2017 Edition, the same parts of the 2019 Edition have also been reviewed as indicated in various sections of the report. This review was conducted by a collaboration of national laboratory and private sector participants with significant industrial experience, including some heavy lifting and deep diving from Clarus Consulting, LLC., all intended to achieve an objective, independent, and practical perspective. The report provides recommendations, descriptions of the evaluation methods, and the source references for the data used. To build confidence required for endorsement of the Code, this review was conducted as a verification and validation of the above Code contents. The objective of verification is to ensure that the Code is free of error – direct or implied; contains the information needed for its use, including proper coverage of the Code-specified materials for the intended application, and completeness and adequacy of references to other portions of the Code. The objective of validation is to authenticate that the Code tabulations and graphs represent design inputs consistent with what are determined using rules and methods specified by the Code. The authentication process used data that were assembled and/or generated independent of Code development, while the methods of analysis followed Code-specified methods where appropriate. The designated portions for this review cover the five alloys codified for high temperature reactor applications in Division 5, i.e. 316 SS, 304 SS, 800H, 2¼Cr-1Mo, and 9Cr-1Mo-V, regarding their general requirements, permitted specifications and design stress intensity values for pressure-retaining applications, deterioration in service, fatigue acceptance test, permissible weld materials, tensile and yield strength, expected minimum stress-to-rupture values (including for Alloy 718), weld stress rupture factors, permissible materials for bolting use, and restricted specifications in certain service applications. Additionally, stress intensity values for bolting materials including 316 SS, 304 SS and alloy 718 were reviewed. Analysis and discussion are also provided on contents outside of these designated Code portions where it was deemed relevant and necessary to develop a technically sound understanding of issues relating to the designated portions. Due to unavailability of sufficient test data on welds during the review period, the weld stress rupture factors in Tables HBB-I-10.14A to E, which cover a total of ten tables for the five alloys welded with twenty-eight different weld metals (some with similar properties), have been deferred to a future review effort. The review identified mainly two types of issues. The first type includes instances where the Code is found factually incomplete or incorrect, such as obsolete materials specifications listings, missing tabulation of stresses for bolting. Changes to the Code are recommended in these cases. The second type of issue includes instances where the Code tabulations and graphs are found to be less conservative than the review analysis results. In these cases, recommendations are made for further review and consideration where the difference in conservatism exceeds 10%, which is our threshold for questioning technical adequacy, meriting a risk assessment by the Nuclear Regulatory Commission and/or reactor designers. It is noted that this effort has been executed using all available data and established methods of analysis, including methods and criteria specified and used by the Code. As such, the findings that are presented in quantitative detail, in a format for convenient comparison with the Code, and with identification of where further review is recommended, should provide a sound technical basis for decisions about quantifying the implications of the reduced design margins and technical adequacy/inadequacy to form a basis for conditioning specific Code tabulation values on endorsement. Recommendations for specific changes to the Code, however, entail design conservatism considerations beyond the scope of this review effort, and are not made in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Triaxial Probe Magnetic Data Analysis

The Triaxial Magnetic Moment Analysis software uses measured magnetic field test data to compute dipole and quadrupole moment information from a hardware element. It is used to support JPL projects needing magnetic control and an understanding of the spacecraft-generated magnetic fields. Evaluation of the magnetic moment of an object consists of three steps: acquisition, conditioning, and analysis. This version of existing software was extensively rewritten for easier data acquisition, data analysis, and report presentation, including immediate feedback to the test operator during data acquisition. While prior JPL computer codes provided the same data content, this program has a better graphic display including original data overlaid with reconstructed results to show goodness of fit accuracy and better appearance of the report graphic page. Data are acquired using three magnetometers and two rotations of the device under test. A clean acquisition user interface presents required numeric data and graphic summaries, and the analysis module yields the best fit (least squares) for the magnetic dipole and/or quadrupole moment of a device. The acquisition module allows the user to record multiple data sets, selecting the best data to analyze, and is repeated three times for each of the z-axial and y-axial rotations. In this update, the y-axial rotation starting position has been changed to an option, allowing either the x- or z-axis to point towards the magnetometer. The code has been rewritten to use three simultaneous axes of magnetic data (three probes), now using two "rotations" of the device under test rather than the previous three rotations, thus reducing handling activities on the device under test. The present version of the software gathers data in one-degree increments, which permits much better accuracy of the fit ted data than the coarser data acquisition of the prior software. The data-conditioning module provides a clean data set for the analysis module. For multiple measurements at a given degree, the first measurement is used. For omitted measurements, the missing field is estimated by linear interpolation between the two nearest measurements. The analysis module was rewritten for the dual rotation, triaxial probe measurement process and now has better moment estimation accuracy, based on the finer one degree of data acquisition resolution. The magnetic moments thus computed are used as an input to summarize the total spacecraft field.

Shultz, Kimberly↗

Multiscale and Machine Learning Modeling for Process-informed Microstructure Prediction in Additively Manufactured Materials Using MALAMUTE

Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy, aims to develop and qualify additively-manufactured materials for nuclear applications. The key challenges to these efforts are the microstructural variabilities observed on the AM products and their impact on the properties and performance of the material in extreme environments. AMMT is using a combination of high-through-put experimental and modeling techniques to accelerate the qualification efforts. Conventionally, in-situ and ex-situ characterizations and testing are performed to correlate different aspects of the AM process to the final product and its performance. However, adopting a trial-and-error approach to experimentally evaluate the vast range of process parameters required to capture the microstructural variabilities is cost-prohibitive. Modeling and simulation provide a comparatively inexpensive way to understand and correlate the microstructural evolution to the processing conditions. The modeling and simulation work-packages within the AMMT program aims to use physics-based and machine learning modeling capabilities to develop a digital twin for AM that can correlate the process conditions to the final product and establish a process-structure-property-performance (PSPP) correlation for AM materials. The melting and subsequent solidification that occurs during the AM process is a complex phenomenon that requires multiscale multiphysics analysis. Idaho National Laboratory’s (INL) Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) software, provides an ideal platform for developing the multiphysics multiscale model to explore the intricacies of the microstructural evolution during the AM processes within a single framework. Furthermore, given that such full-fidelity simulations can be computationally intensive, reduced order models are necessary to explore the PSPP space for AM materials in an efficient, reliable, and cost-effective way. This work package focuses on understanding the role of process variabilities on the various microstructural characteristics of the AM materials. Microstructures unique to AM materials, such as compositional micro-heterogeneity and dislocation cells, are of particular interest here since they can influence the creep properties and radiation performance. In fiscal year (FY) 24, we significantly advanced upon our work in the last fiscal year, both on physics-based and ML models. The alloy solidification model available in MOOSE has been extended to incorporate the thermodynamic properties and free energy relevant to 316SS. The model demonstrates the Cr segregation that occurs during solidifcation. It is demonstrated that rate of solidification and solute segregation is primarily influence by the cooling rate dictating the level of freezing. This work captures the microstructural variabilities at the subgrain level that are often missing in the part-scale models. With an aim to connect the microstructural evolution model to realistic process conditions, a reduced order model is developed for predicting the thermal conditions around meltpool from high-fidelity process simulations. Furthermore, machine learning approach is used to accelerate the temperature prediction during the AM process. In the following years, MALAMUTE will be used to connect different aspects of the models and quantitatively predict the microstructural evolution. The developed ML-based surrogate model will consider the process conditions as the input to predict the microstructural features in a cost-effective way. The generated microstructures can be used by other work packages under AMMT to evaluate the properties and environmental response of the material at the mesoscale. Thus, this work help identify the key microstructural features at the subgrain level that are significant in property/performance prediction of the AM products. This work will provide inputs to the large-scale process variability models to reevaluate and validate assumptions/simplifications made in the part-scale models. Furthermore, through active learning this work will help identify the data need from both modeling and experimental sides for development of a robust digital twin for AM.

36 MATERIALS SCIENCE↗

Advances in Cross-Cutting Ideas for Computational Climate Science

This report presents results from the DOE-sponsored workshop titled, ``Advancing X-Cutting Ideas for Computational Climate Science Workshop,'' known as AXICCS, held on September 12--13, 2016 in Rockville, MD. The workshop brought together experts in climate science, computational climate science, computer science, and mathematics to discuss interesting but unsolved science questions regarding climate modeling and simulation, promoted collaboration among the diverse scientists in attendance, and brainstormed about possible tools and capabilities that could be developed to help address them. Emerged from discussions at the workshop were several research opportunities that the group felt could advance climate science significantly. These include (1) process-resolving models to provide insight into important processes and features of interest and inform the development of advanced physical parameterizations, (2) a community effort to develop and provide integrated model credibility, (3) including, organizing, and managing increasingly connected model components that increase model fidelity yet complexity, and (4) treating Earth system models as one interconnected organism without numerical or data based boundaries that limit interactions. The group also identified several cross-cutting advances in mathematics, computer science, and computational science that would be needed to enable one or more of these big ideas. It is critical to address the need for organized, verified, and optimized software, which enables the models to grow and continue to provide solutions in which the community can have confidence. Effectively utilizing the newest computer hardware enables simulation efficiency and the ability to handle output from increasingly complex and detailed models. This will be accomplished through hierarchical multiscale algorithms in tandem with new strategies for data handling, analysis, and storage. These big ideas and cross-cutting technologies for enabling breakthrough climate simulation advancements also need the "glue" of outreach and learning across the scientific domains to be successful. The workshop identified several strategies to allow productive, continuous engagement across those who have a broad knowledge of the various angles of the problem. Specific ideas to foster education and tools to make material progress were discussed. Examples include follow-on cross-cutting meetings that enable unstructured discussions of the types this workshop fostered. A concerted effort to recruit undergraduate and graduate students from all relevant domains and provide them experience, training, and networking across their immediate expertise is needed. This will broaden and expand their exposure to the future needs and solutions, and provide a pipeline of scientists with a diversity of knowledge and know-how. Providing real-world experience with subject matter experts from multiple angles may also motivate the students to attack these problems and even come up with the missing solutions.

54 ENVIRONMENTAL SCIENCES↗

Advances in Cross-Cutting Ideas for Computational Climate Science

This report presents results from the DOE-sponsored workshop titled, Advancing X-Cutting Ideas for Computational Climate Science Workshop,'' known as AXICCS, held on September 12--13, 2016 in Rockville, MD. The workshop brought together experts in climate science, computational climate science, computer science, and mathematics to discuss interesting but unsolved science questions regarding climate modeling and simulation, promoted collaboration among the diverse scientists in attendance, and brainstormed about possible tools and capabilities that could be developed to help address them. Emerged from discussions at the workshop were several research opportunities that the group felt could advance climate science significantly. These include (1) process-resolving models to provide insight into important processes and features of interest and inform the development of advanced physical parameterizations, (2) a community effort to develop and provide integrated model credibility, (3) including, organizing, and managing increasingly connected model components that increase model fidelity yet complexity, and (4) treating Earth system models as one interconnected organism without numerical or data based boundaries that limit interactions. The group also identified several cross-cutting advances in mathematics, computer science, and computational science that would be needed to enable one or more of these big ideas. It is critical to address the need for organized, verified, and optimized software, which enables the models to grow and continue to provide solutions in which the community can have confidence. Effectively utilizing the newest computer hardware enables simulation efficiency and the ability to handle output from increasingly complex and detailed models. This will be accomplished through hierarchical multiscale algorithms in tandem with new strategies for data handling, analysis, and storage. These big ideas and cross-cutting technologies for enabling breakthrough climate simulation advancements also need the "glue" of outreach and learning across the scientific domains to be successful. The workshop identified several strategies to allow productive, continuous engagement across those who have a broad knowledge of the various angles of the problem. Specific ideas to foster education and tools to make material progress were discussed. Examples include follow-on cross-cutting meetings that enable unstructured discussions of the types this workshop fostered. A concerted effort to recruit undergraduate and graduate students from all relevant domains and provide them experience, training, and networking across their immediate expertise is needed. This will broaden and expand their exposure to the future needs and solutions, and provide a pipeline of scientists with a diversity of knowledge and know-how. Providing real-world experience with subject matter experts from multiple angles may also motivate the students to attack these problems and even come up with the missing solutions.

54 ENVIRONMENTAL SCIENCES↗

Search for the Standard Model Higgs boson produced in association with a W Boson in the isolated-track charged-lepton channel using the Collider Detector at Fermilab

The Higgs boson is the only elementary particle predicted by the Standard Model (SM) that has not yet been observed experimentally. If it exists, it explains the spontaneous electroweak symmetry breaking and the origin of mass for gauge bosons and fermions. We test the validity of the SM by performing a search for the associated production of a Higgs boson and a W boson in the channel where the Higgs boson decays to a bottom-antibottom quark pair and the W boson decays to a charged lepton and a neutrino (the WH channel). We study a dataset of proton-antiproton collisions at a centre-of-mass energy √s = 1.96 TeV provided by the Tevatron accelerator, corresponding to an integrated luminosity of 5.7 fb -1 , and recorded using the Collider Detector at Fermilab (CDF).We select events consistent with the signature of exactly one charged lepton (electron or muon), missing transverse energy due to the undetected neutrino (MET) and two collimated streams of particles (jets), at least one of which is required to be identified as originating from a bottom quark. We improve the discrimination of Higgs signal from backgrounds through the use of an artificial neural network. Using a Bayesian statistical inference approach, we set for each hypothetical Higgs boson mass in the range 100-150 GeV/c 2 with 5 GeV/c 2 increments a 95% credibility level (CL) upper limit on the ratio between the Higgs production cross section times branching fraction and the SM prediction. Our main original contributions are the addition of a novel charged lepton reconstruction algorithm with looser requirements (ISOTRK) with respect the electron or muon tight criteria (TIGHT), as well as the introduction of a novel trigger-combination method that allows to maximize the event yield while avoiding trigger correlations and that is used for the ISOTRK category. The ISOTRK candidate is a high-transverse-momentum good-quality track isolated from other activity in the tracking system and not required to match a calorimeter cluster, as for a tight electron candidate, or an energy deposit in the muon detector, as for a tight muon candidate. The ISOTRK category recovers real charged leptons that otherwise would be lost in the non-instrumented regions of the detector. This allows the reconstruction of more W boson candidates, which in turn increases the number of reconstructed WH signal candidate events, and therefore improves the sensitivity of the WH search. For the TIGHT charged lepton categories, we employ charged-lepton-dedicated triggers to improve the rate of WH signal acceptance during data taking. Since there is no ISOTRK-dedicated trigger at CDF, for the ISOTRK charged lepton category we employ three MET-plus-jets-based triggers. For each trigger we first identify the jet selection where the trigger efficiency is flat with respect to jet information (transverse energy and direction of motion in the transverse plane for the two jets in the event) and then we parametrize the trigger efficiency as a function of trigger MET. On an event-by-event basis, for each trigger we compute a trigger efficiency as a function of trigger parametrization, trigger MET, jet information, trigger prescale and information about whether the trigger is defined or not. For the ISOTRK category we combine the three triggers using a novel method, which allows the combination of any number of triggers in order to maximize the event yield while avoiding trigger correlations. On an event-by-event basis, only the trigger with the largest efficiency is used. By avoiding a logical 'OR' between triggers, the loss in the yield of events accepted by the trigger combination is compensated by a smaller and easier-to-compute corresponding systematic uncertainty. The addition of the ISOTRK charged lepton category to the TIGHT category produces an increase of 33% in the WH signal yield and a decrease of 15.5% to 19.0% in the median expected 95% CL cross-section upper limits across the entire studied Higgs mass interval. The improvement in analysis sensitivity is smaller than the improvement in signal yield because the ISOTRK category has a smaller signal over background ratio than the TIGHT category, due to the looser ISOTRK reconstruction criteria. The observed (median expected) 95% CL SM Higgs upper limits on cross section times branching ratio vary between 2.39 x SM (2.73 x SM) for a Higgs mass of 100 GeV/c 2 to 31.1 x SM (31.2 x SM) for a Higgs mass of 150 GeV/c 2 , while the value for a 115 GeV/c 2 Higgs boson is that of 5.08 x SM (3.79 x SM). The novel trigger combination method is already in use by several CDF analyses. It is applicable to any analysis that uses triggers based on MET and jets, such as supersymmetry searches at the ATLAS and CMS experiments at the Large Hadron Collider. In its most general form, the method can be used by any analysis that combines any number of different triggers.

43 PARTICLE ACCELERATORS↗

Forecasting Multi-Step-Ahead Street-Scale Nuisance Flooding using a seq2seq LSTM Surrogate Model for Real-Time Application in a Coastal-Urban City

In coastal-urban cities facing an elevated risk of nuisance flooding (by rain and tide) due to increased heavy rainfall, sea level rise, urbanization, and aging drainage systems, real-time flood forecasting at the street-scale can provide useful information to transportation decision-makers. Physics-Based Models (PBMs) that offer high accuracy come with high computational runtimes and costs that limit their application for real-time flood forecasting. To address this challenge, Machine Learning (ML) surrogate models trained from PBMs have been proposed to provide street-scale flood forecasts. Previous related studies have focused on using Long Short-Term Memory (LSTM) architectures to model hourly flood depth on streets. While LSTM models can capture input sequences effectively, they fall short in accurately preserving output sequences, limiting their suitability for multi-step-ahead forecasts. The seq2seq LSTM architecture offers a key advantage here by capturing the full sequence of input–output, making it potentially more suitable for multi-step-ahead flood forecasts compared to traditional LSTM models. However, seq2seq LSTM has not been tested for street-scale flood forecasting, particularly for rapidly fluctuating nuisance flooding events which require special attention to its temporal sequences. Hence, in this study, we applied the seq2seq LSTM model to explore multi-step-ahead street-scale nuisance flooding and compared its results to the traditional LSTM model as a benchmark model. LSTM and seq2seq LSTM surrogate models were applied to 22 flood-prone streets in Norfolk, Virginia, as a case study with a 4-hr (short-term) and 8-hr (long-term) lead time. The models were trained with environmental (rainfall and tide) and topographic (elevation, Topographic Wetness Index, and Depth-To-Water) features along with PBM-derived water depths for different storm events. The results demonstrated satisfactory performance of both LSTM and seq2seq LSTM surrogate models throughout the forecast period compared to the PBM. However, the seq2seq LSTM showed lower Mean Absolute Error (MAE)/ Root Mean Square Error (RMSE) and higher Nash–Sutcliffe Efficiency (NSE)/ correlation than the LSTM across most lead times, particularly for long-term forecasting due to its supremacy in handling both input–output sequences together, which is missing in the traditional LSTM. For example, in the long-term, the average RMSE ranges were 0.0268–0.0373 m for LSTM and 0.0226–0.0319 m for seq2seq LSTM, while in the short-term, they were 0.0263–0.0293 m and 0.0261–0.0283 m, respectively. Additionally, while both models exhibited similar performance in distinguishing flooded and non-flooded streets for flood depth ≥ 0.1 m, the seq2seq LSTM model demonstrated superior performance for higher flood depths (such as ≥ 0.2 m and ≥ 0.3 m). Once trained, inference took only 0.09 to 0.11 s (short-term) and 0.30 to 0.35 s (long-term) per storm event for the 22 streets, making the application highly suitable for real-time decision-making during nuisance flood events.

54 ENVIRONMENTAL SCIENCES↗

WE-G-16A-01: Evolution of Radiation Treatment Planning

Welcome and Introduction - Lawrence N. Rothenberg This symposium is one a continuing series of presentations at AAPM Annual Meetings on the historical aspects of medical physics, radiology, and radiation oncology that have been organized by the AAPM History Committee. Information on previous presentations including “Early Developments in Teletherapy” (Indianapolis 2013), “Historical Aspects of Cross-Sectional Imaging” (Charlotte 2012), “Historical Aspects of Brachytherapy” (Vancouver 2011), “50 Years of Women in Medical Physics” (Houston 2008), and “Roentgen's Early Investigations” (Minneapolis 2007) can be found in the Education Section of the AAPM Website. The Austin 2014 History Symposium will be on “Evolution of Radiation Treatment Planning.” Overview - Radhe Mohan Treatment planning is one of the most critical components in the chain of radiation therapy of cancers. Treatment plans of today contain a wide variety of sophisticated information conveying the potential clinical effectiveness of the designed treatment to practitioners. Examples of such information include dose distributions superimposed on three- or even four-dimensional anatomic images; dose volume histograms, dose, dose-volume and dose-response indices for anatomic structures of interest; etc. These data are used for evaluating treatment plans and for making treatment decisions. The current state-of-the-art has evolved from the 1940s era when the dose to the tumor and normal tissues was estimated approximately by manual means. However, the symposium will cover the history of the field from the late-1950's, when computers were first introduced for treatment planning, to the present state involving the use of high performance computing and advanced multi-dimensional anatomic, functional and biological imaging, focusing only on external beam treatment planning. The symposium will start with a general overview of the treatment planning process including imaging, structure delineation, assignment of dose requirements, consideration of uncertainties, selection of beam configurations and shaping of beams, and calculations, optimization and evaluation of dose distributions. This will be followed by three presentations covering the evolution of treatment planning, which parallels the evolution of computers, availability of advanced volumetric imaging and the development of novel technologies such as dynamic multi-leaf collimators and online image guidance. This evolution will be divided over three distinct periods - prior to 1970's, the 2D era; from 1980 to the mid-1990's, the 3D era; and from the mid 1990's to today, the IMRT era. When the World was Flat: The Two-Dimensional Radiation Therapy Era” - Jacob Van Dyk In the 2D era, anatomy was defined with the aid of solder wires, special contouring devices and projection x-rays. Dose distributions were calculated manually from single field, flat surface isodoses on transparencies. Precalculated atlases of generic dose distributions were produced by the International Atomic Energy Agency. Massive time-shared main frames and mini-computers were used to compute doses at individual points or dose distributions in a single plane. Beam shapes were generally rectangular, with wedges, missing tissue compensators and occasional blocks to shield critical structures. Dose calculations were measurement-based or they used primary and scatter calculations based on scatter-air ratio methodologies. Dose distributions were displayed on line printers as alpha-numeric character maps or isodose patterns made with pen plotters. More than Pretty Pictures: 3D Treatment Planning and Conformal Therapy - Benedick A. Fraass The introduction of computed tomography allowed the delineation of anatomy three-dimensionally and, supported partly by contracts from the National Cancer Institute, made possible the introduction and clinical use of 3D treatment planning, leading to development and use of 3D conformal therapy in the 1980's. 3D computer graphics and 3D anatomical structure definitions made possible Beam's Eye View (BEV) displays, making conformal beam shaping and much more sophisticated beam arrangements possible. These conformal plans significantly improved target dose coverage as well as normal tissue sparing. The use of dose volume histograms, gross/clinical/planning target volumes, MRI and PET imaging, multileaf collimators, and computer-controlled treatment delivery made sophisticated planning approaches practical. The significant improvements in dose distributions and analysis achievable with 3D conformal therapy made possible formal dose escalation and normal tissue tolerance clinical studies that set new and improved expectations for improved local control and decreasing complications in many clinical sites. From the Art to the State of the Art: Inverse Planning and IMRT - Thomas R. Bortfeld While the potential of intensity modulation was recognized in the mid- 1980's, intensity-modulated radiotherapy (IMRT) did not become a reality until the mid-1990's. Broad beams of photons could be sub-divided into narrow beamlets whose intensities could be determined using sophisticated optimization algorithms to appropriately balance tumor dose with normal tissue sparing. The development of dynamic multi-leaf collimators (on conventional linear accelerators as well as in helical delivery devices) enabled the efficient delivery of IMRT. The evolution of IMRT planning is continuing in the form of Volumetric Modulated Arc Therapy (VMAT) and through advanced optimization tools, such as multi-criteria optimization, automated IMRT planning, and robust optimization to protect dose distributions against uncertainties. IMRT also facilitates “dose painting” in which different sub-volumes of the target are prescribed different doses. Clearly, these advancements are being made possible by the increasing power and lower cost of computers and developments in other fields such as imaging and operations research. Summary - Radhe Mohan The history does not end here. The advancement of treatment planning is expected to continue, leading to further automation and improvements in conformality and robustness of dose distributions, particularly in the area of particle therapy. Radiobiological modeling will gain emphasis as part of the planning process. Learning Objectives: The scope of changes in technology and the capabilities of radiation treatment planning The impact of these changes in the quality of treatment plans and optimality of dose distributions The impact of development in other fields (imaging, computers, operations research, etc.) on the evolution of radiation treatment planning.

60 APPLIED LIFE SCIENCES↗

Autoload: a pipeline for expanding the holdings of an Institutional Repository enabled by ResourceSync

Providing local access to locally produced content is a primary goal of the Institutional Repository (IR). Guidelines, requirements, and workflows are among the ways in which institutions attempt to ensure this content is deposited and preserved, but some content is always missed. At Los Alamos National Laboratory, the library implemented a service called LANL Research Online (LARO), to provide public access to a collection of publicly shareable LANL researcher publications authored between 2006 and 2016. LARO exposed the fact that we have full text for only about 10% of eligible publications for this time period, despite a review and release requirement that ought to have resulted in a much higher deposition rate. This discovery motivated a new effort to discover and add more full text content to LARO. Autoload attempts to locate and harvest items that were not deposited locally, but for which archivable copies exist. Here we describe the Autoload pipeline prototype and how it aggregates and utilizes Web services including Crossref, SHERPA/RoMEO, and oaDOI as it attempts to retrieve archivable copies of resources. Autoload employs a bootstrapping mechanism based on the ResourceSync standard, a NISO standard for resource replication and synchronization. We implemented support for ResourceSync atop the LARO Solr index, which exposes metadata contained in the local IR. This allowed us to utilize ResourceSync without modifying our IR. Here, we close with a brief discussion of other uses we envision for our ResourceSync-Solr implementation, and describe how a new effort called Signposting can replace cumbersome screen scraping with a robust autodiscovery path to content which leverages Web protocols.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Lessons from Earth Aerobiology for Venus Astrobiology

Venus’s clouds have often suggested as a possible habitat. The constraints governing putative airborne life in such a habitat in turn inform priorities and strategies for remote and in situ exploration, methods by which resulting biosignatures might be detected, exoplanet habitability assessment, and planetary protection concerns. Lessons drawn from studying Earth’s aerobiosphere can help improve this understanding. There are altitude ranges within Venus’s clouds in which temperature, pressure, particle size, and radiation appear to be within the limits of microbial life on Earth, and life cycles involving S- and Fe-based redox metabolism have been proposed. However, given the lack of a habitable surface reservoir, a long-term stable Venus aerobiosphere would require that the reproduction rate of the airborne microbes be faster than the settling rate of airborne microbes due to gravity, or eventually the population would be depleted; put another way, the mean generation time would need to exceed the mean residence time. This creates a joint constraint of aerosol dynamics, potential nutrient availability and energy influx, and bioenergetic costs such as desiccation and radiation damage. Even at an optimistic estimate of 75% H2SO4, Venus aerosol water activity (aw) is still ~0.02, far below the observed microbial growth limit of ~0.6. Long-term desiccation with brief spurts of repair and growth in response to transient water influx, such as from volcanism, is the most likely model for Earth-like life on Venus – a ‘desert bloom’ scenario. Several high-priority science goals – cloud aerosol composition, internal radiative flux, and aerosol residence time and circulation models – thus will also improve our understanding of Venus in an astrobiology context. In Earth’s troposphere, warm water clouds can carry 103 – 105 cells/mL, some metabolically active. However, Earth’s stratospheric sulfate aerosol layer may be a better analogue: supercooled sulfuric acidaerosols (acid weight fraction 0.6 – 0.85, 0.1 – 1 μm diameter, 0.1 – 1 cm-3) with little water activity, long residence times, high UV radiation, and only sporadic influx from surface particle sources. Though ‘hot spots’ can occur associated with tropospheric mixing, viable cells in stratospheric samples are rare (~102 cells/m3), and primarily inactive forms such as spores. It is not yet clear whether such bioaerosols are associated with sulfate aerosols or simply co-located, and reproduction in situ has not yet been observed. In this model of a sparse, largely dormant Venus ecosystem, a single transect on descent is likely to pass through a low-water, inactive region, missing potential signs of habitability or biosignatures. A targeted strategy would sample through an aerial region with some upwelling from surface sources, andtake multiple transects separated in time and space. This is compatible with other in situ science goals seeking to understand the dynamics and heterogeneities of Venus’s clouds.

Lower atmosphere↗

Collective charge and spin dynamics in low-dimensional itinerant electron systems with spin-orbit coupling

This research project in theoretical condensed-matter physics aims to study the collective charge and spin excitations of itinerant quasi-two-dimensional electronic systems such as doped semiconductor nanostructures and carbon-based materials, in the presence of impurities and spin-orbit coupling. In previous work under DOE support, we had discovered an interesting interplay between many-body effects and Rashba and Dresselhaus spin-orbit interactions: contrary to many situation in spintronics where spin-orbit interactions lead to loss of spin coherence, collective spin modes are robust against such coherence losses, and behave as macroscopic quantum objects with a collective magnetic moment precessing about an effective spin-orbit magnetic field. Experimental data for GaAs and CdMnTe quantum wells, obtained by our collaborators at the University Pierre et Marie Curie in Paris, confirmed our predictions and suggested that the observed effect is universal. However, a full theoretical explanation of the effect is still missing, and a wealth of new experimental data remains to be analyzed. We therefore plan to address the following three interrelated objectives: 1. Theory of collective modes in chiral systems. We will explore the various ways in which spin-orbit interactions affect collective spin-density and spin-flip excitations in semiconductor quantum wells. The goal will be to develop a universal formalism for chiral collective modes, with a particular emphasis on analytic results for dispersion relations, which will be compared with experimental results from inelastic light scattering. We will calculate the linewidths of spin plasmons and spin-flip waves due to impurity scattering and spin Coulomb drag from first principles. 2.Many-body effects in noncollinear magnetism. We will describe electronic many-body effects using time-dependent density-functional theory (TDDFT) in the linear-response regime. Quantum wells with spin-orbit coupling are noncollinear itinerant magnetic systems at the crossover between two and three dimensions, which poses subtle challenges for standard exchange-correlation (xc) functionals such as the local-density approximation (LDA). We will develop new orbital-dependent xc functionals for noncollinear magnetic systems by generalizing the so-called STLS approach. 3. Beyond single semiconductor quantum wells. We will consider quantum well bilayer or multilayer systems with spin-orbit coupling, where the coupling between the quasi-2D electron systems leads to acoustic and optical modes (in phase and out of phase). How are these modes affected if the two layers have the same or different chirality? Plasmon modes have been widely studied in graphene and other topological materials, and we will study collective spin excitations in graphene and other carbon-based materials in the presence of spin-orbit coupling and magnetic fields. The potential impact of the proposed research will be both on a fundamental and a practical level. A universal theory of chiral collective modes will advance our understanding of the interplay between many-body effects and spin-orbit interactions in a wide variety of materials and systems. A generalized STLS approach for noncollinear magnetism will be of wide interest beyond the systems studied here. This work may also point towards new ways of manipulating and controlling spin waves in low-dimensional itinerant electron systems: this could be used to encode, transport and process information, which would be of interest to the field of magnon spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Detection of DoS Attacks Using ARFIMA Modeling of GOOSE Communication in IEC 61850 Substations

Integration of Information and Communication Technology (ICT) in modern smart grids (SGs) offers many advantages including the use of renewables and an effective way to protect, control and monitor the energy transmission and distribution. To reach an optimal operation of future energy systems, availability, integrity and confidentiality of data should be guaranteed. Research on the cyber-physical security of electrical substations based on IEC 61850 is still at an early stage. In the present work, we first model the network traffic data in electrical substations, then, we present a statistical Anomaly Detection (AD) method to detect Denial of Service (DoS) attacks against the Generic Object Oriented Substation Event (GOOSE) network communication. According to interpretations on the self-similarity and the Long-Range Dependency (LRD) of the data, an Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) model was shown to describe well the GOOSE communication in the substation process network. Based on this ARFIMA-model and in view of cyber-physical security, an effective model-based AD method is developed and analyzed. Two variants of the statistical AD considering statistical hypothesis testing based on the Generalized Likelihood Ratio Test (GLRT) and the cumulative sum (CUSUM) are presented to detect flooding attacks that might affect the availability of the data. Our work presents a novel AD method, with two different variants, tailored to the specific features of the GOOSE traffic in IEC 61850 substations. The statistical AD is capable of detecting anomalies at unknown change times under the realistic assumption of unknown model parameters. The performance of both variants of the AD method is validated and assessed using data collected from a simulation case study. We perform several Monte-Carlo simulations under different noise variances. The detection delay is provided for each detector and it represents the number of discrete time samples after which an anomaly is detected. In fact, our statistical AD method with both variants (CUSUM and GLRT) has around half the false positive rate and a smaller detection delay when compared with two of the closest works found in the literature. Our AD approach based on the GLRT detector has the smallest false positive rate among all considered approaches. Whereas, our AD approach based on the CUSUM test has the lowest false negative rate thus the best detection rate. Depending on the requirements as well as the costs of false alarms or missed anomalies, both variants of our statistical detection method can be used and are further analyzed using composite detection metrics.

IEC 61850 electrical substations↗

Retaining Systems Engineering Model Meaning Through Transformation: Demo 2

Digital engineering strategies typically assume that digital engineering models interoperate seamlessly across the multiple different engineering modeling software applications involved, such as model- based systems engineering (MBSE), mechanical computer-aided design (MCAD), electrical computer-aided design (ECAD), and other engineering modeling applications. The presumption is that the data schema in these modeling software applications are structured in the familiar flat- tabular schema like any other software application. Engineering domain-specific applications (e.g., systems, mechanical, electrical, simulation) are typically designed to solve domain-specific problems, necessarily excluding explicit representations of non-domain information to help the engineer focus on the domain problems (system definition, design, simulation). Such exclusions become problematic in inter-domain information exchange. The obvious assumptions of one domain might not be so obvious to experts in another domain. Ambiguity in domain-specific language can erode the ability to enable different domain modeling applications to interoperate, unless the underlying language is understood and used as the basis for translation from one application to another. The engineering modeling software application industry has struggled for decades to enable these applications to interoperate. Industry standards have been developed, but they have not unified the industry. Why is this? The authors assert that the industry has relied on traditional database integration methods. The basic issue prohibiting successful application integration then is that traditional database-driven integration does not consider the distinct languages of each domain. An engineering models meaning is expressed through the underlying language of that engineering domain. In essence, traditional integration methods do not retain the semantic context (meaning) of the model. The basis of this research stems from the widely held assumption that systems engineering models are (or can be) structured according to the underlying semantic ontology of the model. This assumption can be imagined from two thoughts. 1) Digital systems engineering models are often represented using graph theory (the graph of a complex systems model can contain millions of nodes and edges). When examining the nodes one at a time and following the outbound edges of each node one by one, one can end up with rudimentary statements about the model (i.e., node A relates to node B), as in a semantic graph. 2) Likewise, from the study of natural languages, a sentence can be structured into unambiguous triples of subject-predicate-object within formal and highly expressive semantic ontologies. The rudimentary statements about a systems model discerned with graph theory closely mimic the triples used in the ontologies that try to structure natural languages. In other words, a systems models semantic graph can be (or is) structured into an ontology. Additionally, it is well established in industry that through natural language processing (NLP), which provides the means to create language structures, that computers can interpret ontological graphs. Therefore, the authors hypothesized that if the integrity of the underlying semantic structure of a systems model is retained, the contextual meaning of the model is retained. By structuring system models into the triples of the underlying ontology during the transformation from one MBSE application to another, the authors have provided a proof of the concept that the meaning of a system model can be retained during transformation. The authors assert that this is the missing ingredient in effective systems model-to-model interoperability. ACKNOWLEDGEMENTS The authors would like to thank the FY19 Model Interoperability team members who provided a solid foundation for the FY20 team to leverage: John McCloud, for the work he did to guide us toward the right use of technology that will appropriately discover and manipulate ontologies. Carlos Tafoya, for the work he did to develop an application programming interface (API)/Adapter that would export ontology-based data from GENESYS. Peter Chandler, for the work he did to architect our overall integration solution, with an eye toward the future that would influence a large-scale federated production-level systems engineering digital model ecosystem.

42 ENGINEERING↗

Final Technical Report

During the past decades, considerable theoretical efforts have been devoted to studying the electronic and geometric structures and related properties of surfaces. Such efforts are particularly important for systems like the actinides for which experimental work is relatively difficult to perform due to material problems and toxicity. The actinides are characterized by a gradual filling of the 5f-electron shell with the degree of localization increasing with the atomic number Z along the last series of the periodic table. The open shell of the 5f electrons determines the atomic, molecular, and solid state properties of the actinide elements and their compounds and understanding the quantum mechanics of the 5f electrons is the defining issue in the chemistry and physics of actinide elements. These elements are also characterized by the increasing prominence of relativistic effects and their studies can, in fact, help us understand the role of relativity throughout the periodic table. However, the electronic and geometric structures of the actinides, specifically the trans-uranium actinides and the roles of the 5f electrons in chemical bonding are still not well understood. This is crucial not only for our understanding of the actinides but also for the fact that the actinides constitute 'the missing link' between the d transition elements and the lanthanides. The 5f orbitals have properties intermediate between those of localized 4f and delocalized 3d orbitals. Thus, a proper understanding of the actinides will help us understand the behavior of the lanthanides and transition metals as well. In fact, there is an urgent need for continued extensive and detailed theoretical research in this area to provide significant and deep understandings of the electronic and geometric structures of the actinides. In this work, we have performed electronic structure studies for plutonium (Pu), americium (Am), and curium (Cm) surfaces, and molecular adsorptions on Pu and Am surfaces. In particular, the region at the boundary of Pu and Am, is widely believed to be the crossover region between d-like itinerant and f-like localized behavior The eventual goal is a complete understanding of the surface chemistry and physics processes of all actinide surfaces, defining the chemistry and physics of such heavy elements. Among the actinides, plutonium, with five 5f electrons in the solid state, is arguably the most complex, fascinating, and enigmatic element known to mankind and has attracted extraordinary scientific and technological interests because of its unique properties, generating a significant body of research in diverse areas, including superconductivity. Pu has, at least, six stable allotropes between room temperature and melting at atmospheric pressure, indicating that the valence electrons can hybridize into a number of complex bonding arrangements. Central and critical questions relate to the electronic structure, localization of the 5f electrons and the magnetism of Pu. For the light-actinides, from Th to Pu, the 5f electrons are believed to be delocalized, hybridizing with the 6d and 7s electrons. For the heavier actinides, Am and beyond, the 5f electrons are localized with the 5f orbitals progressively lower in energy relative to the 6d configuration. Hence, Pu is in a position where the 5f electronic behavior changes from itinerant to localized. As far as magnetism is concerned, a majority of the theoretical calculations continues to claim the existence of magnetism while almost all the experimental results do not find any support for such claims. The second element of interest to us, namely americium, occupies a central position in the actinide series with respect to the involvement of 5f electrons in metallic bonding. It is widely believed that the 5f electrons in Am are localized and that Am undergoes a series of crystallographic phase changes with pressure. Fully-relativistic all electron surface studies of the different phases of Am, initially for the dhcp and the fcc surfaces, can and have provided us with valuable information about chemical bonding in Am and the transitions from f-electron delocalization to f-electron localization in trans-uranium compounds. In particular, a comparative study of the electronic structures of the Pu and Am surfaces using the techniques of all-electron modern density functional theory and beyond can provide significant information about the role of 5f electrons in bond formation as also the localization of the 5f electrons, matters of considerable controversies. The change from metallic 5f bonding into local-moment nonbonding configurations that takes place between Pu and Am is rather unique in the periodic table and is at the very heart of our understanding of electronic structure. We believe that, considering the narrow bandwidth of surface states, any transition from itinerant to localized behavior first takes place at the actinide surfaces with possible reconstructions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessment of the Short-Term Radiometric Stability between Terra MODIS and Landsat 7 ETM+ Sensors

The Landsat 7 (L7) Enhanced Thematic Mapper (ETM+) sensor was launched on April 15th, 1999 and has been in operation for over nine years. It has six reflective solar spectral bands located in the visible and shortwave infrared part of the electromagnetic spectrum (0.5 - 2.5 micron) at a spatial resolution of 30 m. The on-board calibrators are used to monitor the on-orbit sensor system changes. The ETM+ performs solar calibrations using on-board Full Aperture Solar Calibrator (FASC) and the Partial Aperture Solar Calibrator (PASC). The Internal Calibrator Lamp (IC) lamps, a blackbody and shutter optics constitute the on-orbit calibration mechanism for ETM+. On 31 May 2003, a malfunction of the scan-line corrector (SLC) mirror assembly resulted in the loss of approximately 22% of the normal scene area. The missing data affects most of the image with scan gaps varying in width from one pixel or less near the centre of the image to 14 pixels along the east and west edges of the image, creating a wedge-shaped pattern. However, the SLC failure has no impacts on the radiometric performance of the valid pixels. On December 18, 1999, the Moderate Resolution Imaging Spectroradiometer (MODIS) Proto-Flight Model (PFM) was launched on-board the NASA's EOS Terra spacecraft. Terra MODIS has 36 spectral bands with wavelengths ranging from 0.41 to 14.5 micron and collects data over a wide field of view angle (+/-55 deg) at three nadir spatial resolutions of 250 m, 500 in 1 km for bands 1 to 2, 3 to 7, and 8 to 36, respectively. It has 20 reflective solar bands (RSB) with spectral wavelengths from 0.41 to 2.1 micron. The RSB radiometric calibration is performed by using on-board solar diffuser (SD), solar diffuser stability monitor (SDSM), space-view (SV), and spectro-radiometric calibration assembly (SRCA). Through the SV port, periodic lunar observations are used to track radiometric response changes at different angles of incidence (AOI) of the scan mirror. As a part of the AM Constellation satellites, Terra MODIS flies approximately 30 minutes behind L7 ETM+ in the same orbit. The orbit of L7 is repetitive, circular, sunsynchronous, and near polar at a nominal altitude of 705 km (438 miles) at the Equator. The spacecraft crosses the Equator from north to south on a descending node between 10:00 AM and 10:15 AM. Circling the Earth at 7.5 km/sec, each orbit takes nearly 99 minutes. The spacecraft completes just over 14 orbits per day, covering the entire Earth between 81 degrees north and south latitude every 16 days. The longest continuous imaging swath that L7 sensor can collect is for a 14-minute subinterval contact period which is equivalent to 35 full WRS-2 scenes. On the other hand, Terra can provide the entire corresponding orbit with wider swath at any given ETM+ collection without contact time limitation. There are six spectral matching band pairs between MODIS (bands 3, 4, 1, 2, 6, 7) and ETM+ (bands 1, 2, 3, 4, 5, 7) sensor. MODIS has narrower spectral responses than ETM+ in all the bands. A short-term radiometric stability was evaluated using continuous ETM+ scenes within the contact period and the corresponding half orbit MODIS scenes. The near simultaneous earth observations (SNO) were limited by the smaller swath size of ETM+ (187 km) as compared to MODIS (2330 km). Two sets of continuous granules for MODIS and ETM+ were selected and mosaiced based on pixel geolocation information for non cloudy pixels over the North American continent. The Top-of- Atmosphere (TOA) reflectances were computed for the spectrally matching bands between ETM+ and MODIS over the regions of interest (ROI). The matching pixel pairs were aggregated from a finer to a coarser pixel resolution and the TOA reflectance values covering a wide dynamic range of the sensors were compared and analyzed. Considering the uncertainties of the absolute calibration of the both sensors, radiometric stability was verified for the band pairs. The Railroad Valley Playa, Nada (RVPN) was included in the path of this continuous orbit, which served as a verification point between the shortterm and the long-term trending results from previous studies. This work focuses on monitoring the short-term on-orbit stability of MODIS and the ETM+ RSB. It also provides an assessment of the absolute calibration differences between the two sensors over their wide dynamic ranges.

Choi, Taeyoung↗