Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “WINDOWS”

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

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

At least 1,549 records · Page 86

An extension to artifact-free projection overlaps

Purpose: In multipinhole single photon emission computed tomography, the overlapping of projections has been used to increase sensitivity. Avoiding artifacts in the reconstructed image associated with projection overlaps (multiplexing) is a critical issue. In our previous report, two types of artifact-free projection overlaps, i.e., projection overlaps that do not lead to artifacts in the reconstructed image, were formally defined and proved, and were validated via simulations. In this work, a new proposition is introduced to extend the previously defined type-II artifact-free projection overlaps so that a broader range of artifact-free overlaps is accommodated. One practical purpose of the new extension is to design a baffle window multipinhole system with artifact-free projection overlaps. Methods: First, the extended type-II artifact-free overlap was theoretically defined and proved. The new proposition accommodates the situation where the extended type-II artifact-free projection overlaps can be produced with incorrectly reconstructed portions in the reconstructed image. Next, to validate the theory, the extended-type-II artifact-free overlaps were employed in designing the multiplexing multipinhole spiral orbit imaging systems with a baffle window. Numerical validations were performed via simulations, where the corresponding 1-pinhole nonmultiplexing reconstruction results were used as the benchmark for artifact-free reconstructions. The mean square error (MSE) was the metric used for comparisons of noise-free reconstructed images. Noisy reconstructions were also performed as part of the validations. Results: Simulation results show that for noise-free reconstructions, the MSEs of the reconstructed images of the artifact-free multiplexing systems are very similar to those of the corresponding 1-pinhole systems. No artifacts were observed in the reconstructed images. Therefore, the testing results for artifact-free multiplexing systems designed using the extended type-II artifact-free overlaps numerically validated the developed theory. Conclusions: First, the extension itself is of theoretical importance because it broadens the selection range for optimizing multiplexing multipinhole designs. Second, the extension has an immediate application: using a baffle window to design a special spiral orbit multipinhole imaging system with projection overlaps in the orbit axial direction. Such an artifact-free baffle window design makes it possible for us to image any axial portion of interest of a long object with projection overlaps to increase sensitivity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-P-49: Evaluation of Image Quality and Radiation Dose of Various Unenhanced Head CT Protocols

Purpose: To evaluate the diagnostic value of various unenhanced head CT protocols and predicate acceptable radiation dose level for head CT exam. Methods: Our retrospective analysis included 3 groups, 20 patients per group, who underwent clinical routine unenhanced adult head CT examination. All exams were performed axially with 120 kVp. Three protocols, 380 mAs without iterative reconstruction and automAs, 340 mAs with iterative reconstruction without automAs, 340 mAs with iterative reconstruction and automAs, were applied on each group patients respectively. The images were reconstructed with H30, J30 for brain window and H60, J70 for bone window. Images acquired with three protocols were randomized and blindly reviewed by three radiologists. A 5 point scale was used to rate each exam The percentage of exam score above 3 and average scores of each protocol were calculated for each reviewer and tissue types. Results: For protocols without automAs, the average scores of bone window with iterative reconstruction were higher than those without iterative reconstruction for each reviewer although the radiation dose was 10 percentage lower. 100 percentage exams were scored 3 or higher and the average scores were above 4 for both brain and bone reconstructions. The CTDIvols are 64.4 and 57.8 mGy of 380 and 340 mAs, respectively. With automAs, the radiation dose varied with head size, resulting in 47.5 mGy average CTDIvol between 39.5 and 56.5 mGy. 93 and 98 percentage exams were scored great than 3 for brain and bone windows, respectively. The diagnostic confidence level and image quality of exams with AutomAs were less than those without AutomAs for each reviewer. Conclusion: According to these results, the mAs was reduced to 300 with automAs OFF for head CT exam. The radiation dose was 20 percentage lower than the original protocol and the CTDIvol was reduced to 51.2 mGy.

60 APPLIED LIFE SCIENCES↗

Penalized maximum likelihood reconstruction for x-ray differential phase-contrast tomography

Purpose: The purpose of this work is to propose a cost function with regularization to iteratively reconstruct attenuation, phase, and scatter images simultaneously from differential phase contrast (DPC) acquisitions, without the need of phase retrieval, and examine its properties. Furthermore this reconstruction method is applied to an acquisition pattern that is suitable for a DPC tomographic system with continuously rotating gantry (sliding window acquisition), overcoming the severe smearing in noniterative reconstruction. Methods: We derive a penalized maximum likelihood reconstruction algorithm to directly reconstruct attenuation, phase, and scatter image from the measured detector values of a DPC acquisition. The proposed penalty comprises, for each of the three images, an independent smoothing prior. Image quality of the proposed reconstruction is compared to images generated with FBP and iterative reconstruction after phase retrieval. Furthermore, the influence between the priors is analyzed. Finally, the proposed reconstruction algorithm is applied to experimental sliding window data acquired at a synchrotron and results are compared to reconstructions based on phase retrieval. Results: The results show that the proposed algorithm significantly increases image quality in comparison to reconstructions based on phase retrieval. No significant mutual influence between the proposed independent priors could be observed. Further it could be illustrated that the iterative reconstruction of a sliding window acquisition results in images with substantially reduced smearing artifacts. Conclusions: Although the proposed cost function is inherently nonconvex, it can be used to reconstruct images with less aliasing artifacts and less streak artifacts than reconstruction methods based on phase retrieval. Furthermore, the proposed method can be used to reconstruct images of sliding window acquisitions with negligible smearing artifacts.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-F-J-200: An Improved Method for Event Selection in Compton Camera Imaging for Particle Therapy

Purpose: The uncertainty in the beam range in particle therapy limits the conformality of the dose distributions. Compton scatter cameras (CC), which measure the prompt gamma rays produced by nuclear interactions in the patient tissue, can reduce this uncertainty by producing 3D images confirming the particle beam range and dose delivery. However, the high intensity and short time windows of the particle beams limit the number of gammas detected. We attempt to address this problem by developing a method for filtering gamma ray scattering events from the background by applying the known gamma ray spectrum. Methods: We used a 4 stage Compton camera to record in list mode the energy deposition and scatter positions of gammas from a Co-60 source. Each CC stage contained a 4×4 array of CdZnTe crystal. To produce images, we used a back-projection algorithm and four filtering Methods: basic, energy windowing, delta energy (ΔE), or delta scattering angle (Δθ). Basic filtering requires events to be physically consistent. Energy windowing requires event energy to fall within a defined range. ΔE filtering selects events with the minimum difference between the measured and a known gamma energy (1.17 and 1.33 MeV for Co-60). Δθ filtering selects events with the minimum difference between the measured scattering angle and the angle corresponding to a known gamma energy. Results: Energy window filtering reduced the FWHM from 197.8 mm for basic filtering to 78.3 mm. ΔE and Δθ filtering achieved the best results, FWHMs of 64.3 and 55.6 mm, respectively. In general, Δθ filtering selected events with scattering angles < 40°, while ΔE filtering selected events with angles > 60°. Conclusion: Filtering CC events improved the quality and resolution of the corresponding images. ΔE and Δθ filtering produced similar results but each favored different events.

60 APPLIED LIFE SCIENCES↗

SU-F-T-526: A Comparative Study On Gating Efficiency of Varian RPM Device and Calypso System

Purpose: In general, the linear accelerator is gated using respiratory signal obtained by way of external sensors to account for the breathing motion during radiotherapy. One of the commonly used gating devices is the Varian RPM device. Calypso system that uses electromagnetic tracking of implanted or surface transponders could also be used for gating. The aim of this study is to compare the gating efficiency of RPM device and the calypso system by phantom studies. Methods: An ArcCheck insert was used as the phantom with a Gafchromic film placed in its holder. The ArcCheck insert was placed on a Motion Sim platform and moved in the longitudinal direction simulating a respiratory motion with a period of 5 seconds and amplitude of ±6mm. The Gafchromic film was exposed to a 2 × 2cm{sup 2} field, i) with the phantom static, ii) phantom moving but ungated iii) gated with gating window of 2mm and 3mm. This was repeated with Calypso system using surface transponders with the same gating window. The Gafchromic films were read with an EPSON 11000 flatbed scanner and analysed with ‘Medphysto’ software. Results: The full width at half maximum (FWHM) as measured with film at the level of the film holder was 1.65cm when the phantom was static. FWHM measured with phantom moving and without gating was 1.16 cm and penumbra was 7 mm (80–20%) on both sides. When the beam was gated with 2 mm gating window the FWHM was 1.8 cm with RPM device and 1.9 cm with Calypso. Similarly, when the beam was gated with 3 mm window, the FWHM was 1.9cm with RPM device and 2cm with Calypso. Conclusion: This work suggests that the gating efficiency of RPM device is better than that of the Calypso with surface transponder, with reference to the latency in gating.

60 APPLIED LIFE SCIENCES↗

SU-G-BRB-07: Developing a QA Procedure for Gated VMAT SABR Treatments Using 10 MV Beam in Flattening-Filter Free Mode

Purpose: To develop a QA procedure for gated VMAT stereotactic ablative radiotherapy (SABR) treatments. Methods: An interface was constructed to attach the translational stage of a Quasar respiratory motion phantom to a pinpoint ion chamber insert and move the ion chamber inside an ArcCheck diode array. The Quasar phantom controller used a patient specific breathing pattern to translate the ion chamber in a superior-inferior direction inside the ArcCheck. An amplitude-based RPM tracking system was specified to turn the beam on during the exhale phase of the breathing pattern. SABR plans were developed using Eclipse for liver PTVs ranging in size from 3-12 cm in diameter using a 2-arc VMAT technique. Dose was measured in the middle of the penumbra region, where the high dose gradient allowed for sensitive detection of any inaccuracies in gated dose delivery. The overall fidelity of the dose distribution was confirmed using ArcCheck. The sensitivity of the gating QA procedure was investigated with respect to the following four parameters: PTV size, duration of exhale, baseline drift, and gating window size. Results: The difference between the measured dose to a point in the penumbra and the Eclipse calculated dose was under 2% for small residual motions. The QA procedure was independent of PTV size and duration of exhale. Baseline drift and gating window size, however, significantly affected the penumbral dose measurement, with differences of up to 30% compared to Eclipse. Conclusion: This study described a highly sensitive QA procedure for gated VMAT SABR treatments. The QA outcome was dependent on the gating window size and baseline drift. Analysis of additional patient breathing patterns will be required to determine a clinically relevant gating window size and an appropriate tolerance level for this procedure.

60 APPLIED LIFE SCIENCES↗

TH-CD-207B-03: How to Quantify Temporal Resolution in X-Ray MDCT Imaging?

Purpose: In modern CT scanners, a quantitative metric to assess temporal response, namely, to quantify the temporal resolution (TR), remains elusive. Rough surrogate metrics, such as half of the gantry rotation time for single source CT, a quarter of the gantry rotation time for dual source CT, or measurements of motion artifact’s size, shape, or intensity have previously been used. In this work, a rigorous framework which quantifies TR and a practical measurement method are developed. Methods: A motion phantom was simulated which consisted of a single rod that is in motion except during a static period at the temporal center of the scan, termed the TR window. If the image of the motion scan has negligible motion artifacts compared to an image from a totally static scan, then the system has a TR no worse than the TR window used. By repeating this comparison with varying TR windows, the TR of the system can be accurately determined. Motion artifacts were also visually assessed and the TR was measured across varying rod motion speeds, directions, and locations. Noiseless fan beam acquisitions were simulated and images were reconstructed with a short-scan image reconstruction algorithm. Results: The size, shape, and intensity of motion artifacts varied when the rod speed, direction, or location changed. TR measured using the proposed method, however, was consistent across rod speeds, directions, and locations. Conclusion: Since motion artifacts vary depending upon the motion speed, direction, and location, they are not suitable for measuring TR. In this work, a CT system with a specified TR is defined as having the ability to produce a static image with negligible motion artifacts, no matter what motion occurs outside of a static window of width TR. This framework allows for practical measurement of temporal resolution in clinical CT imaging systems. Funding support: GE Healthcare; Conflict of Interest: Employee, GE Healthcare.

60 APPLIED LIFE SCIENCES↗

Sci-Fri PM: Radiation Therapy, Planning, Imaging, and Special Techniques - 06: Patient-specific QA Procedure for Gated VMAT SABR Treatments using 10x Beam in Flattening-Filter Free Mode

Purpose: To develop a highly sensitive patient specific QA procedure for gated VMAT stereotactic ablative radiotherapy (SABR) treatments. Methods: A platform was constructed to attach the translational stage of a Quasar respiratory motion phantom to a pinpoint ion chamber insert and move the ion chamber inside the ArcCheck. The Quasar phantom controller uses a patient-specific breathing pattern to translate the ion chamber in a superior-inferior direction inside the ArcCheck. With this system the ion chamber is used to QA the correct phase of the gated delivery and the ArcCheck diodes are used to QA the overall dose distribution. This novel approach requires a single plan delivery for a complete QA of a gated plan. The sensitivity of the gating QA procedure was investigated with respect to the following parameters: PTV size, exhale duration, baseline drift, gating window size. Results: The difference between the measured dose to a point in the penumbra and the Eclipse calculated dose was under 2% for small residual motions. The QA procedure was independent of PTV size and duration of exhale. Baseline drift and gating window size, however, significantly affected the penumbral dose measurement, with differences of up to 30% compared to Eclipse. Conclusion: This study described a highly sensitive QA procedure for gated VMAT SABR treatments. The QA outcome was dependent on the gating window size and baseline drift. Analysis of additional patient breathing patterns is currently undergoing to determine a clinically relevant gating window size and an appropriate tolerance level for this procedure.

60 APPLIED LIFE SCIENCES↗

Aerial 3D Building Reconstruction from Drone Imagery (A3DBR) v1

This toolkit is composed of several modules for extracting buildings geometrical and thermal characteristics from RGB and thermal imagery captured using a drone. - Building 3D reconstruction module: leverage a photogrammetry software to construct a 3D point cloud from RGB drone imagery, which is then used in conjunction with image processing and geometric methods to extract building footprint and building height (i.e., 3D model of the building). - Windows to wall ratio estimation module: leverage deep learning semantic segmentation modeling to detect windows on 2D drone RGB images. The detected windows are then projected onto the extracted building 3D model (using building 3D reconstruction module) and their area is computed to obtain window to wall ratio estimation. - Thermal anomalies detection module: leverage image processing and machine learning algorithm to detect on 2D drone thermal images potential thermal anomalies within building's facades and roofs.

Granderson, Jessica↗

Eco Look-Ahead Control of Battery Electric Vehicles and Roadway Grade Effects

Roadway grade is one of the major variables that affects a vehicle’s energy consumption. This study demonstrates the potential benefits of an eco look-ahead control for battery electric vehicles (BEVs). The proposed eco look-ahead controller was developed for the lead vehicle of a platoon of BEVs. The eco look-ahead control predicts the optimum speed and acceleration levels within a preset speed window to minimize energy consumption considering the instantaneous energy used and the regenerated energy. The developed BEV eco look-ahead control system integrates a BEV energy consumption model and a powertrain model to save fuel while maintaining vehicle speed within a user-specified window. In addition, the study demonstrates the effects of roadway grade on BEV energy consumption. The results show that the energy consumption of BEVs are significantly reduced on downhill roads compared with internal combustion engine vehicles owing to regenerative braking. Specifically, testing showed that BEVs could even produce extra energy on downhill roads. The study also tested the eco look-ahead control on an 18-km section of I-81 and found that the system produced savings of 8.51% and 26.88% in the uphill and downhill directions, respectively. This study demonstrated that regenerative energy in BEVs is a critical factor in energy efficiency and the developed eco look-ahead control significantly improved the energy efficiency for BEVs.

Engineering↗

Theoretical analysis of the role of complex transition dipole phase in XUV transient-absorption probing of charge migration

We theoretically investigate the role of complex dipole phase in the attosecond probing of charge migration. The iodobromoacetylene ion (ICCBr + ) is considered as an example, in which one can probe charge migration by accessing both the iodine and bromine ends of the molecule with different spectral windows of an extreme-ultraviolet (XUV) pulse. The analytical expression for transient absorption shows that the site-specific information of charge migration is encoded in the complex phase of cross dipole products for XUV transitions between the I-4 d and Br-3 d spectral windows. Ab-initio quantum chemistry calculations on ICCBr + reveal that there is a constant π phase difference between the I-4 d and Br-3 d transient-absorption spectral windows, irrespective of the fine-structure energy splittings. Transient absorption spectra are simulated with a multistate model including the complex dipole phase, and the results correctly reconstruct the charge-migration dynamics via the quantum beats in the two element spectral windows, exhibiting out-of-phase oscillations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Disordered spontaneously buckled optical gratings for improved lighting applications

Technologies for directing daylight deeper within a building space are highly sought after for energy efficiency applications in order to offset artificial lighting costs and to improve workplace productivity via the use of natural light. Vertical window coatings that can perform this task by redistributing sunlight deeper into a space are especially attractive as they are significantly more straightforward to incorporate into a wide variety of architectures as well as to retrofit into existing facades as compared to roof-based skylights or bulky horizontal daylight shelf type options. The potential energy savings are even greater when one takes into account the fact that such technologies would mitigate harsh glare, allowing window shades to be open for longer portions of the day. However, low-cost and readily scalable techniques are essential for widespread adoption of these window coating technologies. Here, we describe a potentially low-cost method to create a window coating that could enhance daylight penetration, requiring only a thin film of polymeric material deposited on an elastomeric substrate. The resulting structure is a disordered, spontaneously buckled optical grating that spreads incident light without noticeable chromatic dispersion due to its stochastic patterning. The described method has the potential to improve energy efficiency and while maintaining acceptable optical clarity.

De Vetter, Brent M.↗

Chemical and bonding analysis of liquids using liquid cell electron microscopy

Liquid cell transmission electron microscopy (TEM) has become an essential tool for studying the structure and properties of both hard and soft condensed-matter samples, as well as liquids themselves. Liquid cell sample holders, often consisting of two thin window layers separating the liquid sample from the high vacuum of the microscope column, have been designed to control in situ conditions, including temperature, voltage/current, or flow through the window region. While high-resolution and time-resolved TEM imaging probes the structure, shape, and dynamics of liquid cell samples, information about the chemical composition and spatially resolved bonding is often difficult to obtain due to the liquid thickness, the window layers, the holder configuration, or beam-induced radiolysis. Here, we review different approaches to quantitative liquid cell electron microscopy, including recent developments to perform energy-dispersive x-ray and electron energy-loss spectroscopy experiments on samples in a liquid environment or the liquid itself. We also cover graphene liquid cells and other ultrathin window layer holders.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning-Based Identification of the Interface Regions for Coupling Local and Nonlocal Models

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. However, the coupling process can be challenging, requiring expertise to identify the interface between local and nonlocal regions. Here, this study introduces a machine learning-based approach to automatically detect the regions in which the local and nonlocal models should be used. The method uses loading functions evaluated at grid points to decide the model selection at those points. Training of the networks is based on datasets provided by classes of loading functions for which reference coupling configurations are computed using accurate coupled solutions, where accuracy is measured in terms of the relative error between the solution to the coupling approach and the solution to the nonlocal model. We study two approaches that vary in data structure. The first, the full-domain input data approach, uses the entire load vector and outputs a complete label vector, performing a global classification. The second, a window-based approach, processes loads into windows and addresses the problem as a node-wise classification where each window's central point is classified individually. The classification problems are solved via deep learning algorithms based on convolutional neural networks. The performance of these approaches is studied on one-dimensional numerical examples using F1-scores and accuracy metrics. Notably, the windowing approach achieves an accuracy of 0.96 and an F1-score of 0.97, highlighting its potential to automate coupling processes effectively and enhance computational efficiency in material science applications.

97 MATHEMATICS AND COMPUTING↗

Regional Characteristics of Observable Foreshocks

Abstract Measures of foreshock occurrence are systematically examined using earthquake catalogs for eight regions (Italy, southern California, northern California, Costa Rica, Onshore Japan, Alaska, Turkey, and Greece) after imposing a magnitude ≥3.0 completeness level. Foreshocks are identified using three approaches: a magnitude-dependent space + fixed-time windowing method, a nearest-neighbor clustering method, and a modified magnitude-dependent space + variable-time windowing method. The method with fixed-time windows systematically yields higher counts of foreshocks than the other two clustering methods. We find similar counts of foreshocks across the three methods when the magnitude aperture is equalized by including only earthquakes in the magnitude range M*−2≤ M< M*, in which M* is the mainshock magnitude. For most of the catalogs (excluding Italy and southern California), the measured b-values of the foreshocks of all region-specific mainshocks are lower by 0.1–0.2 than b-values of respective aftershocks. Allowing for variable-time windows results in relatively high probabilities of having at least one foreshock in Italy (∼43%–56%), compared to other regional catalogs. Foreshock probabilities decrease to 14%–41% for regions such as Turkey, Greece, and Costa Rica. Similar trends are found when requiring at least five foreshocks in a sequence to be considered. Estimates of foreshock probabilities for each mainshock are method dependent; however, consistent regional trends exist regardless of method, with regions such as Italy and southern California producing more observable foreshocks than Turkey and Greece. Some regions with relatively high background seismicity have comparatively low probabilities of detectable foreshock activity when using methods that account for variable background, possibly due to depletion of near-failure fault conditions by background activity.

Geochemistry & Geophysics↗

High Sodium Simulant Testing To Support SB8 Sludge Preparation

Scoping studies were completed for high sodium simulant SRAT/SME cycles to determine any impact to CPC processing. Two SRAT/SME cycles were performed with simulant having sodium supernate concentration of 1.9M at 130% and 100% of the Koopman Minimum Acid requirement. Both of these failed to meet DWPF processing objectives related to nitrite destruction and hydrogen generation. Another set of SRAT/SME cycles were performed with simulant having a sodium supernate concentration of 1.6M at 130%, 125%, 110%, and 100% of the Koopman Minimum Acid requirement. Only the run at 110% met DWPF processing objectives. Neither simulant had a stoichiometric factor window of 30% between nitrite destruction and excessive hydrogen generation. Based on the 2M-110 results it was anticipated that the 2.5M stoichiometric window for processing would likely be smaller than from 110-130%, since it appeared that it would be necessary to increase the KMA factor by at least 10% above the minimum calculated requirement to achieve nitrite destruction due to the high oxalate content. The 2.5M-130 run exceeded the DWPF hydrogen limits in both the SRAT and SME cycle. Therefore, testing of this wash endpoint was halted. This wash endpoint with this minimum acid requirement and mercury-noble metal concentration profile appears to be something DWPF should not process due to an overly narrow window of stoichiometry. The 2M case was potentially processable in DWPF, but modifications would likely be needed in DWPF such as occasionally accepting SRAT batches with undestroyed nitrite for further acid addition and reprocessing, running near the bottom of the as yet ill-defined window of allowable stoichiometric factors, potentially extending the SRAT cycle to burn off unreacted formic acid before transferring to the SME cycle, and eliminating formic acid additions in the frit slurry.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Systematic Evaluation of Traditional and Non-Traditional Data Streams for Integrated Global Biosurveillance (Final Report)

Living in a closely connected and highly mobile world presents many new mechanisms for rapid disease spread and in recent years, global disease surveillance has become a high priority. In addition, much like the contribution of non-traditional medicine to curing diseases, non-traditional data streams are being considered of value in disease surveillance. Los Alamos National Laboratory (LANL) was tasked by the Defense Threat Reduction Agency (DTRA)-Joint Science and Technology Office (JSTO) to determine the relevance of data streams for an integrated global biosurveillance system through the use of defined metrics and methodologies. The long-term, broad objective of this project was to provide information and analysis that can be leveraged from existing and developing national and international disease surveillance systems and methodologies, to create a global disease monitoring system, ultimately providing decision makers with timely information to prepare for and mitigate the spread of disease. Specifically, this project evaluated data streams that are currently being used in surveillance systems and data streams that had the potential for being used in surveillance to enable early disease detection. LANL’s focus was on infectious disease surveillance and the term “biosurveillance” will be used to refer to this scope hereafter. Data stream evaluation was conducted using two different methods: Multi criteria decision analysis (MCDA) using a tool called Logical Decisions ® that assigns utility scores to data streams based on weighted metrics and assigned values specific to data stream categories; and a Surveillance Window concept developed at LANL that assigns a window or windows of time specific to a disease within which information coming from various data streams can be determined to have utility. A cross method analysis was performed between the surveillance window based evaluation and the MCDA-based evaluation to identify data stream categories that showed high utility for both methods. In addition, algorithms that can integrate useful data streams to facilitate early disease detection were also examined. This project provided an understanding of data streams relevant to early warning, early detection and monitoring of disease outbreaks. The report describes the results of LANL’s evaluation of 16 data stream categories using the two methods, and recommendations for next steps. The project evolved significantly over its course, and significant improvements in our approaches were made compared to the initial proposed processes. A key outcome of this effort was that LANL was able to bring together several diverse entities involved in disease surveillance and lay the foundation for new collaborations that straddled military and civilian health surveillance. The robust evaluation framework developed by LANL has generated significant interest and there is interest in adopting this framework for various applications both within disease surveillance as well as other public health initiatives. Finally, LANL developed the Biosurveillance Resource Directory (BRD), a relational database that underwent pilot testing by members of the human, plant and animal disease surveillance community. The BRD is intended to be a global resource to facilitate rapid information access.

60 APPLIED LIFE SCIENCES↗

Carbon Data Assimilation with a Coupled Ensemble Kalman Filter (Final Technical Report)

We proposed (and accomplished) the development of an Ensemble Kalman Filter (EnKF) approach for the estimation of surface carbon fluxes as if they were parameters, augmenting the model with them. Our system is quite different from previous approaches, such as carbon flux inversions, 4D-­Var, and EnKF with approximate background error covariance (Peters et al., 2008). We showed (using observing system simulation experiments, OSSEs) that these differences lead to a more accurate estimation of the evolving surface carbon fluxes at model grid-scale resolution. The main properties of the LETKF-­C are: a) The carbon cycle LETKF is coupled with the simultaneous assimilation of the standard atmospheric variables, so that the ensemble wind transport of the CO2 provides an estimation of the carbon transport uncertainty. b) The use of an assimilation window (6hr) much shorter than the months-long windows used in other methods. This avoids the inevitable “blurring” of the signal that takes place in long windows due to turbulent mixing since the CO2 does not have time to mix before the next window. In this development we introduced new, advanced techniques that have since been adopted by the EnKF community (Kang, 2009, Kang et al., 2011, Kang et al. 2012). These advances include “variable localization” that reduces sampling errors in the estimation of the forecast error covariance, more advanced adaptive multiplicative and additive inflations, and vertical localization based on the time scale of the processes. The main result has been obtained using the LETKF-­C with all these advances, and assimilating simulated atmospheric CO2 observations from different observing systems (surface flask observations of CO2 but no surface carbon fluxes observations, total column CO2 from GoSAT/OCO-­2, and upper troposphere AIRS retrievals). After a spin-­up of about one month, the LETKF-­C succeeded in reconstructing the true evolving surface fluxes of carbon at a model grid resolution. When applied to the CAM3.5 model, the LETKF gave very promising results as well, although only one month is available.

54 ENVIRONMENTAL SCIENCES↗