Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level set method”

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 847 records · Page 47

Monte Carlo Control Loops for Cosmic Shear Cosmology with DES Year 1

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, Point Spread Function (PSF), an Approximate Bayesian Computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and the fast estimation of the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $\Omega_m$, $\sigma_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$, defined as $S_8D_{\rm{IA}} = \sigma_8(\Omega_m/0.3)^{0.5}D_{\rm{IA}}$, where $D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find $S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of roughly 60\% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation of Material Models For Automotive Carbon Fiber Composite Structures Via Physical And Crash Testing (VMM Composites Project)

As automotive fuel economy requirements increase, the push for reducing overall vehicle weight will likely include the consideration of materials that have not previously been part of mainstream vehicle design and manufacturing, including carbon fiber composites. Vehicle manufacturers currently rely on computer-aided engineering (CAE) methods as part of the design and development process, so going forward, the ability to accurately and predictably model carbon fiber composites will be necessary. If composites are to be used for structural components, this need applies to both, crash and quasi-static modeling. This final report covers the results of a five-year, $6.89M, 50% cost-shared research project between Department of Energy (DOE) and the US Advanced Materials Partnership (USAMP) under Cooperative Agreement DE-EE-0005661 known as “Validation of Material Models for Automotive Carbon Fiber Composite Structures Via Physical and Crash Testing (VMM).” The objective of the VMM Composites Project was to validate and assess the ability of physics-based material models to predict crash performance of automotive primary load-carrying carbon fiber composite structures. Simulation material models that were evaluated included micro-mechanics based meso-scale models developed by the University of Michigan (UM) and micro-plane models by Northwestern University (NWU) under previous collaborations with the DOE and Automotive Composites Consortium/USAMP, as well as five commercial crash codes: LS-DYNA, RADIOSS, VPS/PAM-CRASH, Abaqus, and GENOA-MCQ. CAE predictions obtained from seven organizations were compared with experimental results from quasi-static testing and dynamic crash testing of a thermoset carbon fiber composite front-bumper and crush-can (FBCC) system gathered under multiple loading conditions. This FBCC design was developed to demonstrate progressive crush, virtual simulation, tooling, fabrication, assembly, non-destructive evaluation and crash testing advances in order to assess the correlation of the predicted results to the physical tests. The FBCC was developed to meet a goal of 30-35% mass reduction while aiming for equivalent energy absorption as a steel component for which baseline experimental results were obtained from testing in the same crash modes. The project also evaluated crash performance of thermoplastic composite structures fabricated from commercial prepreg materials and low cost carbon fiber sourced from Oak Ridge National Laboratory. The VMM Project determined that no set of predictions from a CAE supplier were found to be universally accurate among all the six crash modes evaluated. In general, crash modes that were most dependent on the properties of the prepreg were more accurate than those that were dependent on the behavior of the joints. The project found that current CAE modeling methods or best practices for carbon fiber composites have not achieved standardization, and accuracy of CAE is highly reliant on the experience of its users. Coupon tests alone are not sufficient to develop an accurate material model, but it is necessary to bridge the gap between the coupon data and performance of the actual structure with a series of subcomponent level tests. Much of the unreliability of the predictions can be attributed to shortcomings in our ability to mathematically link the effects of manufacturing and material variability into the material models. This is a subject of ongoing research in the industry. The final report is organized by key technical tasks to describe how the validation project developed, modeled and compared crash data obtained on the composite FBCC to the multiple sets of CAE predictions. Highlights of the report include a discussion of the quantitative comparison between predictions and experimental data, as well as an in-depth discussion of remaining technological gaps that exist in the industry, which are intended to spur innovations and improvements in CAE technology.

36 MATERIALS SCIENCE↗

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗

Methods for evaluation and treatment of data shift

This is a code repository for a set of tools for detecting and mitigating data shifts in machine learning. The goal of the tools is to provide capabilities for determining when new data sets differ from training data sets and for adapting existing models to new data or correcting data shifts (via domain adaptation). The components will be written in Python, a high-level programming language that takes advantage of the Python ecosystem of high-quality open-source packages for machine learning and signal processing.

Parikh, Nidhi↗

SU-E-T-427: Cell Surviving Fractions Derived From Tumor-Volume Variation During Radiotherapy for Non-Small Cell Lung Cancer: Comparison with Predictive Assays

Purpose: To show that a distribution of cell surviving fractions S{sub 2} in a heterogeneous group of patients can be derived from tumor-volume variation curves during radiotherapy for non-small cell lung cancer. Methods: Our analysis was based on two data sets of tumor-volume variation curves for heterogeneous groups of 17 patients treated for nonsmall cell lung cancer with conventional dose fractionation. The data sets were obtained previously at two independent institutions by using megavoltage (MV) computed tomography (CT). Statistical distributions of cell surviving fractions S{sup 2} and cell clearance half-lives of lethally damaged cells T1/2 have been reconstructed in each patient group by using a version of the two-level cell population tumor response model and a simulated annealing algorithm. The reconstructed statistical distributions of the cell surviving fractions have been compared to the distributions measured using predictive assays in vitro. Results: Non-small cell lung cancer presents certain difficulties for modeling surviving fractions using tumor-volume variation curves because of relatively large fractional hypoxic volume, low gradient of tumor-volume response, and possible uncertainties due to breathing motion. Despite these difficulties, cell surviving fractions S{sub 2} for non-small cell lung cancer derived from tumor-volume variation measured at different institutions have similar probability density functions (PDFs) with mean values of 0.30 and 0.43 and standard deviations of 0.13 and 0.18, respectively. The PDFs for cell surviving fractions S{sup 2} reconstructed from tumor volume variation agree with the PDF measured in vitro. Comparison of the reconstructed cell surviving fractions with patient survival data shows that the patient survival time decreases as the cell surviving fraction increases. Conclusion: The data obtained in this work suggests that the cell surviving fractions S{sub 2} can be reconstructed from the tumor volume variation curves measured during radiotherapy with conventional fractionation. The proposed method can be used for treatment evaluation and adaptation.

60 APPLIED LIFE SCIENCES↗

Case studies in bias reduction and inference for electronic health record data with selection bias and phenotype misclassification

Electronic health records (EHR) are not designed for population‐based research, but they provide easy and quick access to longitudinal health information for a large number of individuals. Many statistical methods have been proposed to account for selection bias, missing data, phenotyping errors, or other problems that arise in EHR data analysis. However, addressing multiple sources of bias simultaneously is challenging. We developed a methodological framework (R package, SAMBA ) for jointly handling both selection bias and phenotype misclassification in the EHR setting that leverages external data sources. These methods assume factors related to selection and misclassification are fully observed, but these factors may be poorly understood and partially observed in practice. As a follow‐up to the methodological work, we demonstrate how to apply these methods for two real‐world case studies, and we evaluate their performance. In both examples, we use individual patient‐level data collected through the University of Michigan Health System and various external population‐based data sources. In case study (a), we explore the impact of these methods on estimated associations between gender and cancer diagnosis. In case study (b), we compare corrected associations between previously identified genetic loci and age‐related macular degeneration with gold standard external summary estimates. These case studies illustrate how to utilize diverse auxiliary information to achieve less biased inference in EHR‐based research.

60 APPLIED LIFE SCIENCES↗

Carbon Ion Radiation Therapy With Concurrent Gemcitabine for Patients With Locally Advanced Pancreatic Cancer

Purpose: To determine, in the setting of locally advanced pancreatic cancer, the maximum tolerated dose of carbon ion radiation therapy (C-ion RT) and gemcitabine dose delivered concurrently and to estimate local effect and survival. Methods and Materials: Eligibility included pathologic confirmation of pancreatic invasive ductal carcinomas and radiographically unresectable disease without metastasis. Concurrent gemcitabine was administered on days 1, 8, and 15, and the dose levels were escalated from 400 to 1000 mg/m{sup 2} under the starting dose level (43.2 GyE) of C-ion RT. The dose levels of C-ion RT were escalated from 43.2 to 55.2 GyE at 12 fractions under the fixed recommended gemcitabine dose determined. Results: Seventy-six patients were enrolled. Among the 72 treated patients, dose-limiting toxicity was observed in 3 patients: grade 3 infection in 1 patient and grade 4 neutropenia in 2 patients. Only 1 patient experienced a late grade 3 gastric ulcer and bleeding 10 months after C-ion RT. The recommended dose of gemcitabine with C-ion RT was found to be 1000 mg/m{sup 2}. The dose of C-ion RT with the full dose of gemcitabine (1000 mg/m{sup 2}) was safely increased to 55.2 GyE. The freedom from local progression rate was 83% at 2 years using the Response Evaluation Criteria in Solid Tumors. The 2-year overall survival rates in all patients and in the high-dose group with stage III (≥45.6 GyE) were 35% and 48%, respectively. Conclusions: Carbon ion RT with concurrent full-dose gemcitabine was well tolerated and effective in patients with unresectable locally advanced pancreatic cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Examination of How Well Long-Range-Corrected Density Functionals Satisfy the Ionization Energy Theorem

For this work, we calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r 2 ) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee–Yang–Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee–Yang–Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An assessment of low-lying excitation energies and triplet instabilities of organic molecules with an ab initio Bethe-Salpeter equation approach and the Tamm-Dancoff approximation

The accurate prediction of singlet and triplet excitation energies is an area of intense research of significant fundamental interest and critical for many applications. Most calculations of singlet and triplet energies use time-dependent density functional theory (TDDFT) in conjunction with an approximate exchange-correlation functional. In this work, we examine and critically assess an alternative method for predicting low-lying neutral excitations with similar computational cost, the ab initio Bethe-Salpeter equation (BSE) approach, and compare results against high-accuracy wavefunction-based methods. We consider singlet and triplet excitations of 27 prototypical organic molecules, including members of Thiel’s set, the acene series, and several aromatic hydrocarbons exhibiting charge-transfer-like excitations. Analogous to its impact in TDDFT, we find that the Tamm-Dancoff approximation (TDA) overcomes triplet instabilities in the BSE approach, improving both triplet and singlet energetics relative to higher level theories. Lastly, we find that BSE-TDA calculations built on effective DFT starting points, such as those utilizing optimally tuned range-separated hybrid functionals, can yield accurate singlet and triplet excitation energies for gas-phase organic molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

Neighborhood sociome factors and pediatric asthma exacerbations: Protective role of tree crown density and importance of pharmacy access in Chicago's south side

Abstract Background Pediatric asthma exacerbations remain a critical public health concern, particularly in historically underserved urban settings. Objective This study investigates sociome factors—the social context of disease—associated with asthma exacerbations among children living in Chicago's South Side, leveraging clinical and publicly available generalizable census tract‐level datasets from agencies including ChiVes, the City of Chicago Data Portal, EPA, Census Bureau, HUD, NOAA, and more. The aim is to uncover novel hypotheses for potential new interventions. Methods A generalized linear model assessed associations with the outcome of asthma exacerbations while accounting for clustering at the patient level. Predictors included all variables from the Sociome Data Commons, including social, environmental, behavioral, economic, housing, and school variables. Results Predictors of decreased risk included patient age (+4.8 years, −22%), tree crown density (+6% coverage, −17%), parks per acre (+0.41, −8%), and labor market engagement (+0.8 points, −9%). Conversely, predictors of increased risk included increased distance to the nearest pharmacy (+0.28 miles, +12%), limited English skills (+2.3%, +10%), higher inequality (+0.08 points, +8%), and visits in the Spring (+11%) and Fall (+20%). Conclusion The results suggest that tree crown density, a novel finding in the context of asthma exacerbations, may play a protective role. Limited access to health care facilities such as pharmacies continues to complicate care. Clinical Implications These findings provide hypotheses for future interventions for long‐standing asthma disparities.

Allergy↗

SU-E-J-262: Segmentation in Therapy: Impact of Display

Purpose: To determine the calibration of monitors used for structure delineation in radiotherapy and to investigate the impact of window/level on the delineations. Method and Materials: hotometers were used to measure ambient light and luminance for 12 monitors used for delineating structures in radiotherapy. Luminance was determined using AAPM TG-18 patterns TG18-LN12-01 to TG18-LN12-18. The DICOM grayscale display function (GSDF) was calculated for each monitor. To investigate the impact of window and level eight observers each outlined two sets of structures (four structures per set) on CT simulation images. In the first set window and level were initialized at random values and the observers adjusted as needed. In the second set window and level were locked at preselected values for all observers. Results: Ambient light variability at different displays was 24% and was higher than radiology reading rooms. Monitor luminance curves were substantially different than the corresponding GSDF. Monitors demonstrated higher contrast than the GSDF in the lower third of driving levels and lower contrast in the upper third of driving levels. Increased contrast at low driving levels may be to compensate for high ambient light which would otherwise reduce contrast. Abrupt changes in luminance and loss of all contrast at both upper and lower driving levels was observed. Interobserver variability was 379HU for window and 217HU for level. Outline variability was lower when the window and level settings were similar and when they were preselected. Structures without sharp boundaries (e.g., partial volume and a wide window value) resulted in substantial interobserver variability even when window/level were preselected. Conclusion: Substantial variation in monitor calibration and ambient light exists in radiotherapy. Further, the substantially different window/level contributes to outline variability. Future work will quantify the contributions of monitor calibration and window/level on outline variability and accuracy.

60 APPLIED LIFE SCIENCES↗

Physical Uncertainty Bounds (PUB)

This paper introduces and motivates the need for a new methodology for determining upper bounds on the uncertainties in simulations of engineered systems due to limited fidelity in the composite continuum-level physics models needed to simulate the systems. We show that traditional uncertainty quantification methods provide, at best, a lower bound on this uncertainty. We propose to obtain bounds on the simulation uncertainties by first determining bounds on the physical quantities or processes relevant to system performance. By bounding these physics processes, as opposed to carrying out statistical analyses of the parameter sets of specific physics models or simply switching out the available physics models, one can obtain upper bounds on the uncertainties in simulated quantities of interest.

36 MATERIALS SCIENCE↗

Evaluating IMRT and VMAT dose accuracy: Practical examples of failure to detect systematic errors when applying a commonly used metric and action levels

Purpose: This study (1) examines a variety of real-world cases where systematic errors were not detected by widely accepted methods for IMRT/VMAT dosimetric accuracy evaluation, and (2) drills-down to identify failure modes and their corresponding means for detection, diagnosis, and mitigation. The primary goal of detailing these case studies is to explore different, more sensitive methods and metrics that could be used more effectively for evaluating accuracy of dose algorithms, delivery systems, and QA devices.Methods: The authors present seven real-world case studies representing a variety of combinations of the treatment planning system (TPS), linac, delivery modality, and systematic error type. These case studies are typical to what might be used as part of an IMRT or VMAT commissioning test suite, varying in complexity. Each case study is analyzed according to TG-119 instructions for gamma passing rates and action levels for per-beam and/or composite plan dosimetric QA. Then, each case study is analyzed in-depth with advanced diagnostic methods (dose profile examination, EPID-based measurements, dose difference pattern analysis, 3D measurement-guided dose reconstruction, and dose grid inspection) and more sensitive metrics (2% local normalization/2 mm DTA and estimated DVH comparisons).Results: For these case studies, the conventional 3%/3 mm gamma passing rates exceeded 99% for IMRT per-beam analyses and ranged from 93.9% to 100% for composite plan dose analysis, well above the TG-119 action levels of 90% and 88%, respectively. However, all cases had systematic errors that were detected only by using advanced diagnostic techniques and more sensitive metrics. The systematic errors caused variable but noteworthy impact, including estimated target dose coverage loss of up to 5.5% and local dose deviations up to 31.5%. Types of errors included TPS model settings, algorithm limitations, and modeling and alignment of QA phantoms in the TPS. Most of the errors were correctable after detection and diagnosis, and the uncorrectable errors provided useful information about system limitations, which is another key element of system commissioning.Conclusions: Many forms of relevant systematic errors can go undetected when the currently prevalent metrics for IMRT/VMAT commissioning are used. If alternative methods and metrics are used instead of (or in addition to) the conventional metrics, these errors are more likely to be detected, and only once they are detected can they be properly diagnosed and rooted out of the system. Removing systematic errors should be a goal not only of commissioning by the end users but also product validation by the manufacturers. For any systematic errors that cannot be removed, detecting and quantifying them is important as it will help the physicist understand the limits of the system and work with the manufacturer on improvements. In summary, IMRT and VMAT commissioning, along with product validation, would benefit from the retirement of the 3%/3 mm passing rates as a primary metric of performance, and the adoption instead of tighter tolerances, more diligent diagnostics, and more thorough analysis.

61 RADIATION PROTECTION AND DOSIMETRY↗

Mass Spectrometric Determination of Uranium and Thorium in High Radiopurity Polymers Using Ultra Low Background Electroformed Copper Crucibles for Dry Ashing

A rapid new method for determining the U and Th mass concentrations in high radiopurity plastics is described, consisting of 1) dry ashing the plastic sample and tracers in low mass crucibles made of ultra low background electroformed copper (ULB EF-Cu) foil cut and folded into boats, 2) dissolving both the ash and the boat in acid, 3) performing a column separation to remove copper, and 4) determining the elements of interest by isotope dilution mass spectrometry. This method was demonstrated on both unfluorinated and fluorinated plastics, demonstrating high tracer recoveries and detection limits to pg/g (i.e., parts per trillion) levels or below, corresponding to μBq/kg of material. Samples of biomedical polyester (Max-Prene® 955) and a fluoropolymer (polyvinylidene fluoride, PVDF) were analyzed in powder raw material forms as well as solids in the form of pellets or injection molded parts. The polyester powder contained 6 pg/g and 2 pg/g for Th and U respectively. These levels correspond to 25 and 25 μBq/kg radioactivity, respectively. Determinations on samples of PVDF powder were typically below 1 pg/g for Th and 2 pg/g for U, corresponding to 4 and 25 μBq/kg radioactivity, respectively. The use of low mass ULB EF-Cu boats for dry ashing successfully overcame the problem of crucible-generated contaminants in the analysis; absolute detection limits, calculated as 3 × standard deviation of the process blanks, were typically 20-100 fg within a sample set. Complete dissolution of the ash and low mass boat provided high tracer recoveries, and provides a convincing method to recover both the tracer and sample isotopes when full equilibration of tracer isotopes with sample isotopes is not possible prior to beginning chemical sample processing on solids.

EFCu Electroformed copper↗

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations.

59 BASIC BIOLOGICAL SCIENCES↗

Estimating groundwater dynamics at a Colorado River floodplain site using historical hydrological data and climate information

Long-term prediction of groundwater dynamics is important for assessing water resources and their impacts on biogeochemical cycling. However, estimating future groundwater dynamics is challenging due to the wide range of spatiotemporal scales in hydrological processes and uncertainty in future climate conditions. In this study, we develop a Bayesian model to combine small-scale historical hydrological data with large-scale climate information to estimate groundwater dynamics at a floodplain site in Rifle, Colorado. Although we have only a few years of groundwater elevation measurements, we have 47 years of streamflow data from a gaging station approximately 43 km upstream and long-term climate prediction on the Upper Colorado River Basin. To estimate future daily groundwater dynamics, we first develop a time series model to downscale the monthly streamflow derived from climate information to daily streamflow, and then transform the daily streamflow to groundwater dynamics at the downstream floodplain site. We use Monte Carlo methods to estimate future groundwater dynamics at the site through sampling from the joint posterior probability distribution. The results suggest that although future groundwater levels are expected to be similar to the current levels, the timing of the high groundwater levels is predicted to occur about 1 month earlier. The developed framework is extendable to other sites to estimate future groundwater dynamics given disparate data sets and climate projections. Additionally, the obtained estimates are being used as input to a site-specific watershed reactive transport models to predict how climate-induced changes will influence future biogeochemical cycling relevant to a variety of ecosystem services.

58 GEOSCIENCES↗

Sizing ramping reserve using probabilistic solar forecasts: A data-driven method

Ramping products have been introduced or proposed in several U.S. power markets to mitigate the impact of load and renewable uncertainties on market efficiency and reliability. Current methods often rely on historical data to estimate the requirements of ramping products and fail to take into account the effects of the latest weather conditions and their uncertainties, which could lead to overly conservative or insufficient requirements. This study proposes a k-nearest-neighbor-based method to give weather-informed estimates of ramping needs based on short-term probabilistic solar irradiance forecasts. Forecasts from multiple sites are employed in conjunction with principal component analysis to derive numerical classifiers to characterize system-level weather conditions. In addition, we develop a data-driven method to optimize the model parameters in a rolling-forward manner. By using real-world data from the California Independent System Operator, we design two metrics to evaluate method performance: 1) frequency of shortage and 2) oversupply of ramping product. Our proposed method presents advantages in comparison with the baseline and a set of benchmark methods: without compromising system reliability, it reduces system ramping requirements by up to 25%, therefore improving both system reliability and economics.

14 SOLAR ENERGY↗

Consensus Recommendations for Radiation Therapy Contouring and Treatment of Vulvar Carcinoma

Purpose: The purpose of this study was to develop a radiation therapy (RT) contouring atlas and recommendations for women with postoperative and locally advanced vulvar carcinoma. Methods and Materials: An international committee of 35 expert gynecologic radiation oncologists completed a survey of the treatment of vulvar carcinoma. An initial set of recommendations for contouring was discussed and generated by consensus. Two cases, 1 locally advanced and 1 postoperative, were contoured by 14 physicians. Contours were compared and analyzed using an expectation-maximization algorithm for simultaneous truth and performance level estimation (STAPLE), and a 95% confidence interval contour was developed. The level of agreement among contours was assessed using a kappa statistic. STAPLE contours underwent full committee editing to generate the final atlas consensus contours. Results: Analysis of the 14 contours showed substantial agreement, with kappa statistics of 0.69 and 0.64 for cases 1 and 2, respectively. There was high specificity for both cases (≥99%) and only moderate sensitivity of 71.3% and 64.9% for cases 1 and 2, respectively. Expert review and discussion generated consensus recommendations for contouring target volumes and treatment for postoperative and locally advanced vulvar cancer. Conclusions: These consensus recommendations for contouring and treatment of vulvar cancer identified areas of complexity and controversy. Given the lack of clinical research evidence in vulvar cancer radiation therapy, the committee advocates a conservative and consistent approach using standardized recommendations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗