Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse sampling”

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 109 records · Page 6

Phase-weighted slant stacking for surface wave dispersion measurement

SUMMARY Surface wave retrieval from ambient noise records using seismic interferometry techniques has been widely used for multiscale shear wave velocity (Vs) imaging. One key step during Vs imaging is the generation of dispersion spectra and the extraction of a reliable dispersion curve from the retrieved surface waves. However, the sparse array geometry usually affects the ability for high-frequency (>1 Hz) seismic signals’ acquisition. Dispersion measurements are degraded by array response due to sparse sampling and often present smeared dispersion spectra with sidelobe artefacts. Previous studies usually focus on interferograms’ domain (e.g. cross-correlation function) and attempt to enhance coherent signals before dispersion measurement. We propose an alternative technique to explicitly deblur dispersion spectra through use of a phase-weighted slant-stacking algorithm. Numerical examples demonstrate the strength of the proposed technique to attenuate array responses as well as incoherent noise. Three different field examples prove the flexibility and superiority of the proposed technique: the first data set consists of ambient noise records acquired using a nodal seismometer array; the second data set utilizes distributed acoustic sensing (DAS) and a marine fibre-optic cable to acquire a similar ambient noise data set; the last data set is a vibrator-based active-source surface wave data. The enhanced dispersion measurements provide cleaner and higher-resolution spectra without distortions which will assist both human interpreters as well as ML algorithms in efficiently picking curves for subsequent Vs inversion.

58 GEOSCIENCES↗

Cosmological cross-correlations and nearest neighbour distributions

ABSTRACT Cross-correlations between data sets are used in many different contexts in cosmological analyses. Recently, k-nearest neighbour cumulative distribution functions (kNN-CDF) were shown to be sensitive probes of cosmological (auto) clustering. In this paper, we extend the framework of NN measurements to describe joint distributions of, and correlations between, two data sets. We describe the measurement of joint kNN-CDFs, and show that these measurements are sensitive to all possible connected N-point functions that can be defined in terms of the two data sets. We describe how the cross-correlations can be isolated by combining measurements of the joint kNN-CDFs and those measured from individual data sets. We demonstrate the application of these measurements in the context of Gaussian density fields, as well as for fully non-linear cosmological data sets. Using a Fisher analysis, we show that measurements of the halo-matter cross-correlations, as measured through NN measurements are more sensitive to the underlying cosmological parameters, compared to traditional two-point cross-correlation measurements over the same range of scales. Finally, we demonstrate how the NN cross-correlations can robustly detect cross-correlations between sparse samples – the same regime where the two-point cross-correlation measurements are dominated by noise.

79 ASTRONOMY AND ASTROPHYSICS↗

Distributed Stochastic Optimization of a Neural Representation Network for Time-Space Tomography Reconstruction

4D time-space reconstruction of dynamic events or deforming objects using X-ray computed tomography (CT) is an important inverse problem in non-destructive evaluation. Conventional back-projection based reconstruction methods assume that the object remains static for the duration of several tens or hundreds of X-ray projection measurement images (reconstruction of consecutive limited-angle CT scans). However, this is an unrealistic assumption for many in-situ experiments that causes spurious artifacts and inaccurate morphological reconstructions of the object. To solve this problem, we propose to perform a 4D time-space reconstruction using a distributed implicit neural representation (DINR) network that is trained using a novel distributed stochastic training algorithm. Our DINR network learns to reconstruct the object at its output by iterative optimization of its network parameters such that the measured projection images best match the output of the CT forward measurement model. Here, we use a forward measurement model that is a function of the DINR outputs at a sparsely sampled set of continuous valued 4D object coordinates. Unlike previous neural representation architectures that forward and back propagate through dense voxel grids that sample the object's entire time-space coordinates, we only propagate through the DINR at a small subset of object coordinates in each iteration resulting in an order-of-magnitude reduction in memory and compute for training. DINR leverages distributed computation across several compute nodes and GPUs to produce high-fidelity 4D time-space reconstructions. We use both simulated parallel-beam and experimental cone-beam X-ray CT datasets to demonstrate the superior performance of our approach.

36 MATERIALS SCIENCE↗

SU-E-I-20: Dead Time Count Loss Compensation in SPECT/CT: Projection Versus Global Correction

Purpose: To compare projection-based versus global correction that compensate for deadtime count loss in SPECT/CT images. Methods: SPECT/CT images of an IEC phantom (2.3GBq 99mTc) with ∼10% deadtime loss containing the 37mm (uptake 3), 28 and 22mm (uptake 6) spheres were acquired using a 2 detector SPECT/CT system with 64 projections/detector and 15 s/projection. The deadtime, Ti and the true count rate, Ni at each projection, i was calculated using the monitor-source method. Deadtime corrected SPECT were reconstructed twice: (1) with projections that were individually-corrected for deadtime-losses; and (2) with original projections with losses and then correcting the reconstructed SPECT images using a scaling factor equal to the inverse of the average fractional loss for 5 projections/detector. For both cases, the SPECT images were reconstructed using OSEM with attenuation and scatter corrections. The two SPECT datasets were assessed by comparing line profiles in xyplane and z-axis, evaluating the count recoveries, and comparing ROI statistics. Higher deadtime losses (up to 50%) were also simulated to the individually corrected projections by multiplying each projection i by exp(-a*Ni*Ti), where a is a scalar. Additionally, deadtime corrections in phantoms with different geometries and deadtime losses were also explored. The same two correction methods were carried for all these data sets. Results: Averaging the deadtime losses in 5 projections/detector suffices to recover >99% of the loss counts in most clinical cases. The line profiles (xyplane and z-axis) and the statistics in the ROIs drawn in the SPECT images corrected using both methods showed agreement within the statistical noise. The count-loss recoveries in the two methods also agree within >99%. Conclusion: The projection-based and the global correction yield visually indistinguishable SPECT images. The global correction based on sparse sampling of projections losses allows for accurate SPECT deadtime loss correction while keeping the study duration reasonable.

60 APPLIED LIFE SCIENCES↗

SU-E-I-21: Deformation Mapping and Shape Prediction with 3D Tumor Volume Morphing

Purpose: To compare projection-based versus global correction that compensate for deadtime count loss in SPECT/CT images. Methods: SPECT/CT images of an IEC phantom (2.3GBq 99mTc) with ∼10% deadtime loss containing the 37mm (uptake 3), 28 and 22mm (uptake 6) spheres were acquired using a 2 detector SPECT/CT system with 64 projections/detector and 15 s/projection. The deadtime, Ti and the true count rate, Ni at each projection, i was calculated using the monitor-source method. Deadtime corrected SPECT were reconstructed twice: (1) with projections that were individually-corrected for deadtime-losses; and (2) with original projections with losses and then correcting the reconstructed SPECT images using a scaling factor equal to the inverse of the average fractional loss for 5 projections/detector. For both cases, the SPECT images were reconstructed using OSEM with attenuation and scatter corrections. The two SPECT datasets were assessed by comparing line profiles in xyplane and z-axis, evaluating the count recoveries, and comparing ROI statistics. Higher deadtime losses (up to 50%) were also simulated to the individually corrected projections by multiplying each projection i by exp(-a*Ni*Ti), where a is a scalar. Additionally, deadtime corrections in phantoms with different geometries and deadtime losses were also explored. The same two correction methods were carried for all these data sets. Results: Averaging the deadtime losses in 5 projections/detector suffices to recover >99% of the loss counts in most clinical cases. The line profiles (xyplane and z-axis) and the statistics in the ROIs drawn in the SPECT images corrected using both methods showed agreement within the statistical noise. The count-loss recoveries in the two methods also agree within >99%. Conclusion: The projection-based and the global correction yield visually indistinguishable SPECT images. The global correction based on sparse sampling of projections losses allows for accurate SPECT deadtime loss correction while keeping the study duration reasonable.

60 APPLIED LIFE SCIENCES↗

WE-G-207-03: Mask Guided Image Reconstruction (MGIR): A Novel Method for Ultra-Low-Dose 3D and Enhanced 4D Cone-Beam Computer-Tomography

Purpose: Recently, compressed sensing (CS) based iterative reconstruction (IR) method is receiving attentions to reconstruct high quality cone beam computed tomography (CBCT) images using sparsely sampled or noisy projections. The aim of this study is to develop a novel baseline algorithm called Mask Guided Image Reconstruction (MGIR), which can provide superior image quality for both low-dose 3DCBCT and 4DCBCT under single mathematical framework. Methods: In MGIR, the unknown CBCT volume was mathematically modeled as a combination of two regions where anatomical structures are 1) within the priori-defined mask and 2) outside the mask. Then we update each part of images alternatively thorough solving minimization problems based on CS type IR. For low-dose 3DCBCT, the former region is defined as the anatomically complex region where it is focused to preserve edge information while latter region is defined as contrast uniform, and hence aggressively updated to remove noise/artifact. In 4DCBCT, the regions are separated as the common static part and moving part. Then, static volume and moving volumes were updated with global and phase sorted projection respectively, to optimize the image quality of both moving and static part simultaneously. Results: Examination of MGIR algorithm showed that high quality of both low-dose 3DCBCT and 4DCBCT images can be reconstructed without compromising the image resolution and imaging dose or scanning time respectively. For low-dose 3DCBCT, a clinical viable and high resolution head-and-neck image can be obtained while cutting the dose by 83%. In 4DCBCT, excellent quality 4DCBCT images could be reconstructed while requiring no more projection data and imaging dose than a typical clinical 3DCBCT scan. Conclusion: The results shown that the image quality of MGIR was superior compared to other published CS based IR algorithms for both 4DCBCT and low-dose 3DCBCT. This makes our MGIR algorithm potentially useful in various on-line clinical applications. Provisional Patent: UF#15476; WGS Ref. No. U1198.70067US00.

60 APPLIED LIFE SCIENCES↗

SU-G-IeP1-13: Sub-Nyquist Dynamic MRI Via Prior Rank, Intensity and Sparsity Model (PRISM)

Purpose: Accelerated dynamic MRI is important for MRI guided radiotherapy. Inspired by compressive sensing (CS), sub-Nyquist dynamic MRI has been an active research area, i.e., sparse sampling in k-t space for accelerated dynamic MRI. This work is to investigate sub-Nyquist dynamic MRI via a previously developed CS model, namely Prior Rank, Intensity and Sparsity Model (PRISM). Methods: The proposed method utilizes PRISM with rank minimization and incoherent sampling patterns for sub-Nyquist reconstruction. In PRISM, the low-rank background image, which is automatically calculated by rank minimization, is excluded from the L1 minimization step of the CS reconstruction to further sparsify the residual image, thus allowing for higher acceleration rates. Furthermore, the sampling pattern in k-t space is made more incoherent by sampling a different set of k-space points at different temporal frames. Results: Reconstruction results from L1-sparsity method and PRISM method with 30% undersampled data and 15% undersampled data are compared to demonstrate the power of PRISM for dynamic MRI. Conclusion: A sub- Nyquist MRI reconstruction method based on PRISM is developed with improved image quality from the L1-sparsity method.

60 APPLIED LIFE SCIENCES↗

TU-AB-BRA-06: BEST IN PHYSICS (JOINT IMAGING-THERAPY): An MRI Compatible Externally and Internally Deformable Lung Motion Phantom for Multi-Modality IGRT

Purpose: MRI has become an attractive tool for tumor motion management. Current MR-compatible phantoms are only capable of reproducing translational motion. This study describes the construction and validation of a more realistic, MRI-compatible lung phantom that is deformable internally as well as externally. We demonstrate a radiotherapy application of this phantom by validating the geometric accuracy of the open-source deformable image registration software NiftyReg (UCL, UK). Methods: The outer shell of a commercially-available dynamic breathing torso phantom was filled with natural latex foam with eleven water tubes. A rigid foam cut-out served as the diaphragm. A high-precision programmable, in-house, MRI-compatible motion platform was used to drive the diaphragm. The phantom was imaged on a 3T scanner (Philips, Ingenia). Twenty seven tumor traces previously recorded from lung cancer patients were programmed into the phantom and 2D+t image sequences were acquired using a sparse-sampling sequence k-t BLAST (accn=3, resolution=0.66×0.66×5mm3; acquisition-time=110ms/slice). The geometric fidelity of the MRI-derived trajectories was validated against those obtained via fluoroscopy using the on board kV imager on a Truebeam linac. NiftyReg was used to perform frame by frame deformable image registration. The location of each marker predicted by using NiftyReg was compared with the values calculated by intensity-based segmentation on each frame. Results: In all cases, MR trajectories were within 1 mm of corresponding fluoroscopy trajectories. RMSE between centroid positions obtained from segmentation with those obtained by NiftyReg varies from 0.1 to 0.21 mm in the SI direction and 0.08 to 0.13 mm in the LR direction showing the high accuracy of deformable registration. Conclusion: We have successfully designed and demonstrated a phantom that can accurately reproduce deformable motion under a variety of imaging modalities including MRI, CT and x-ray fluodoscopy, making it an invaluable research tool for validating novel motion management strategies. This work was partially supported through research funding from National Institutes of Health (R01CA169102).

60 APPLIED LIFE SCIENCES↗

Characteristics of thermospheric infrasound predicted using ray tracing and weakly non-linear waveform analyses

Infrasonic signals refracted by thermal gradients in the rarefied upper atmosphere are modeled using a combination of ray tracing and weak shock theory to develop an understanding of thermospheric infrasound signals produced by energetic, transient sources. Canonical arrival structures in the form of u-wave signatures are identified for returns refracted at lower altitudes within the thermosphere, and possible multi-pathing produced by effective sound speed inflections are investigated to elucidate more complex arrival structures, which are found to be spatially localized. Furthermore, variability in the source characteristics is investigated and it is found that whereas some waveform phase information is lost due to finite amplitude effects, arrival characteristics are strongly dependent on the peak overpressure near the source. Variability in the propagation path is considered using archived atmospheric specifications and implies that despite uncertainties related to the dynamic and sparsely sampled nature of the atmosphere, thermospheric signatures might be useful in estimating the yield for explosive sources. Last, thermospheric arrivals from a failed rocket launch, as well as those from several large chemical explosions, are analyzed and it is found that qualitative trends match those predicted, and analyses here provide additional insight into such signatures.

42 ENGINEERING↗

Assimilation of Time-Averaged Pseudoproxies for Climate Reconstruction

The efficacy of a novel ensemble data assimilation (DA) technique is examined in the climate field reconstruction (CFR) of surface temperature. A minimalistic, computationally inexpensive DA technique is employed that requires only a static ensemble of climatologically plausible states. Pseudoproxy experiments are performed with both general circulation model (GCM) and Twentieth Century Reanalysis (20CR) data by reconstructing surface temperature fields from a sparse network of noisy pseudoproxies. The DA approach is compared to a conventional CFR approach based on principal component analysis (PCA) for experiments on global domains. DA outperforms PCA in reconstructing global-mean temperature in all experiments and is more consistent across experiments, with a range of time series correlations of 0.69–0.94 compared to 0.19–0.87 for the PCA method. DA improvements are even more evident in spatial reconstruction skill, especially in sparsely sampled pseudoproxy regions and for 20CR experiments. It is hypothesized that DA improves spatial reconstructions because it relies on coherent, spatially local temperature patterns, which remain robust even when glacial states are used to reconstruct nonglacial states and vice versa. These local relationships, as utilized by DA, appear to be more robust than the orthogonal patterns of variability utilized by PCA. Comparing results for GCM and 20CR data indicates that pseudoproxy experiments that rely solely on GCM data may give a false impression of reconstruction skill.

54 ENVIRONMENTAL SCIENCES↗

Designing and Testing Broadly-Protective Filoviral Vaccines Optimized for Cytotoxic T-Lymphocyte Epitope Coverage

We report the rational design and in vivo testing of mosaic proteins for a polyvalent pan-filoviral vaccine using a computational strategy designed for the Human Immunodeficiency Virus type 1 (HIV-1) but also appropriate for Hepatitis C virus (HCV) and potentially other diverse viruses. Mosaics are sets of artificial recombinant proteins that are based on natural proteins. The recombinants are computationally selected using a genetic algorithm to optimize the coverage of potential cytotoxic T lymphocyte (CTL) epitopes. Because evolutionary history differs markedly between HIV-1 and filoviruses, we devised an adapted computational technique that is effective for sparsely sampled taxa; our first significant result is that the mosaic technique is effective in creating high-quality mosaic filovirus proteins. The resulting coverage of potential epitopes across filovirus species is superior to coverage by any natural variants, including current vaccine strains with demonstrated cross-reactivity. The mosaic cocktails are also robust: mosaics substantially outperformed natural strains when computationally tested against poorly sampled species and more variable genes. Furthermore, in a computational comparison of cross-reactive potential a design constructed prior to the Bundibugyo outbreak performed nearly as well against all species as an updated design that included Bundibugyo. These points suggest that the mosaic designs would be more resilient than natural-variant vaccines against future Ebola outbreaks dominated by novel viral variants. We demonstrate in vivo immunogenicity and protection against a heterologous challenge in a mouse model. This design work delineates the likely requirements and limitations on broadly-protective filoviral CTL vaccines.

59 BASIC BIOLOGICAL SCIENCES↗

An approach for estimating the uncertainty in ParaDiS predictions

This report outlines an approach for computing the uncertainties in the predictions of computationally expensive models. While general, we use ParaDiS , a dislocation dynamics simulator originating in Lawrence Livermore National Laboratory, as the target application. ParaDiS is a mesoscale model, and uses submodels constructed/upscaled from microscale (molecular statics and dynamics) simulations. ParaDiS outputs, in turn, are upscaled and used in continuum (macroscale) simulations, e.g., those performed by ALE3D. This report addresses how one may quantify the uncertainties introduced by upscaling (both from microscale to mesoscale, and mesoscale to continuum), and the dependence of uncertainties in ParaDiS predictions on those of the inputs. This dependence is established via sensitivity analysis, and we address how this may be performed with a minimum of ParaDiS runs, given its immense computational cost. This includes constructing a smaller version of the model, sparse sampling of the parameter space, and exploiting the asymptotic nature of the time-evolution of the model outputs. The report concludes with a discussion of the computational resources required to perform this uncertainty quantification study.

97 MATHEMATICS AND COMPUTING↗

Climate data, analysis and models for the study of natural variability and anthropogenic change

Gridded Temperature Under prior/current support, we completed and published (Jones et al., 2012) the fourth major update to our global land dataset of near-surface air temperatures, CRUTEM4. This is one of the most widely used records of the climate system, having been updated, maintained and further developed with DoE support since the 1980s. We have continued to update the CRUTEM4 (Jones et al., 2012) database that is combined with marine data to produce HadCRUT4 (Morice et al., 2012). The emphasis in our use of station temperature data is to access as many land series that have been homogenized by National Meteorological Services (NMSs, including NCDC/NOAA, Asheville, NC). Unlike the three US groups monitoring surface temperatures in a similar way, we do not infill areas that have no or missing data. We can only infill such regions in CRUTEM4 by accessing more station temperature series. During early 2014, we have begun the extensive task of updating as many of these series as possible using data provided by some NMSs and also through a number of research projects and programs around the world. All the station data used in CRUTEM4 have been available since 2009, but in Osborn and Jones (2014) we have made this more usable using a Google Earth interface (http://www.cru.uea.ac.uk/cru/data/crutem/ge/ ). We have recently completed the update of our infilled land multi-variable dataset (CRU TS 3.10, Harris et al., 2014). This additionally produces complete land fields (except for the Antarctic) for temperature, precipitation, diurnal temperature range, vapour pressure and sunshine/cloud. Using this dataset we have calculated sc-PDSI (self-calibrating Palmer Drought Severity Index) data and compared with other PDSI datasets (Trenberth et al., 2014). Also using CRU TS 3.10 and Reanalysis datasets, we showed no overall increase in global temperature variability despite changing regional patterns (Huntingford et al., 2013). Harris et al. (2014) is an update of an earlier dataset (Mitchell and Jones, 2005) which also had earlier DoE support. The earlier dataset has been cited over 1700 times according to ResearcherID on 31/July/2014 and the recent paper has already been cited 22 times. Analyses of Temperature Data Using the ERA-Interim estimate of the absolute surface air temperature of the Earth (instead of in the more normal form of anomalies) we compared the result against estimates we produced in 1999 with earlier DoE support. The two estimates are surprisingly close (differing by a couple of tenths of a degree Celsius), with the average temperature of the world (for 1981-2010) being very close to 14°C (Jones and Harpham, 2013). We have assessed ERA-Interim against station temperatures from manned and automatic weather station measurements across the Antarctic (Jones and Lister, 2014). Agreement is generally excellent across the Antarctic Peninsula and the sparsely sampled western parts of Antarctica. Differences tend to occur over eastern Antarctica where ERA-Interim is biased warm (up to 6°C) in the interior of the continent and biased cool (up to 6°C) for some of the coastal locations. Opportunities presented themselves during 2012 for collaborative work with a couple of Chinese groups. Three papers develop new temperature series for China as a whole and also for the eastern third of China (Wang et al., 2014, Cao et al., 2013 and Zhao et al., 2014). A dataset of ~400 daily Chinese temperature stations has been added to the CRU datasets. The latter paper finds that urban effects are generally about 10% of the long-term warming trend across eastern China. A fourth paper (Wang et al., 2013) illustrates issues with comparisons between reanalyses and surface temperatures across China, a method that has been widely used by some to suggest urban heating effects are much larger in the region. ERA-Interim can be used but NCEP/NCAR comparisons are very dependent on the period analysed. Earlier a new temperature dataset of homogenized records was developed for China (Li et al., 2009). Urbanization has also been addressed for London (Jones and Lister, 2009) where two rural sites have not warmed more than a city centre site since 1900. Additionally, in Ethymiadis and Jones (2010) we show that land air temperatures agree with marine data around coastal areas, further illustrating that urbanization is not a major component of large-scale surface air temperature change. Early instrumental data (before the development of modern thermometer screens) have always been suspected of being biased warm in summer, due to possible direct exposure to the sun. Two studies (Böhm et al., 2010 and Brunet et al., 2010) show this for the Greater Alpine Region (GAR) and for mainland Spain respectively. The issue is important before about 1870 in the GAR and before about 1900 in Spain. After correction for the problems, summer temperature estimates before these dates are cooler by about 0.4°C. In Jones and Wigley (2010), we discussed the importance of the biases in global temperature estimation. Exposure and to a lesser extent urbanization are the most important biases for the land areas, but both are dwarfed by the necessary adjustments for bucket SST measurements before about 1950. Individual station homogeneity is only important at the local scale. This was additionally illustrated by Hawkins and Jones (2013) where we replicated the temperature record developed by Guy Stewart Callendar in papers in 1938 and 1961. Analyses of Daily Climate Data Work here indicates that ERA-Interim (at least in Europe, Cornes and Jones, 2013, discussed in more detail in this proposal) can be used to monitor extremes (using the ETCCDI software – see Zhang et al., 2011). Additionally, also as a result of Chinese collaboration, a new method of daily temperature homogenization has been developed (Li et al., 2014). In Cornes and Jones (2011) we assessed storm activity in the northeast Atlantic region using daily gridded data. Even though the grid resolution is coarse (5° by 5° lat/long) the changes in storm activity are similar to those developed from the pressure triangle approach with station data. Analyses of humidity and pressure data In Simmons et al. (2010) we showed a reduction in relative humidity over low-latitude and mid-latitude land areas for the 10 years to 2008, based on monthly anomalies of surface air temperature and humidity from ECMWF reanalyses (ERA-40 and ERA-Interim) and our earlier land-only dataset (CRUTEM3) and synoptic humidity observations (HadCRUH). Updates of this station-based humidity dataset (now called HadISDH) extend the record, showing continued reductions (Willett et al., 2013). Analyses of Proxy Temperature Data In Vinther et al. (2010), relationships between the seasonal stable isotope data from Greenland Ice Cores and Greenland and Icelandic instrumental temperatures were investigated for the past 150-200 years. The winter season stable isotope data are found to be influenced by the North Atlantic Oscillation (NAO) and very closely related to SW Greenland temperatures. The summer season stable isotope data display higher correlations with Icelandic summer temperatures and North Atlantic SST conditions than with local SW Greenland temperatures. In Jones et al. (2014) we use these winter isotope reconstructions to show the expected inverse correlation (due to the NAO) with winter-season documentary reconstructions from the Netherlands and Sweden over the last 800 years. Finally, in this section Jones et al. (2013) shows the agreement between tree-ring width measurements from Northern Sweden and Finland and an assessment of the link to explosive volcanic eruptions. An instrumental record for the region in the early 19th century indicates that the summer of 1816 was only slightly below normal, explaining why this year has normal growth for both ring width and density. GCM/RCM/Reanalysis Evaluation In this section we have intercompared daily temperature extremes across Europe in Cornes and Jones (2013) using station data, E-OBS and ERA-Interim. We have additionally considered the impact of the urban issue on the global scale using the results of the Compo et al. (2011) Reanalyses, 20CR. These only make use of SST and station pressure data. Across the world’s land areas, they indicate similar warming since 1900 to that which has occurred (Compo et al., 2013), again illustrating that urbanization is not the cause of the long-term warming. Changes in HadCRUH global land surface specific humidity and CRUTEM3 surface temperatures from 1973 to 1999 were compared to the CMIP3 archive of climate model simulations with 20th Century forcings (Willett et al., 2010). The models reproduce the magnitude of observed interannual variance over all large regions. Observed and modelled trends and temperature-humidity relationships are comparable with the exception of the extra-tropical Southern Hemisphere where observations exhibit no trend but models exhibit moistening.

54 ENVIRONMENTAL SCIENCES↗

Modular hybrid modeling to increase efficiency, explore structural uncertainty, and allow multidimensional complexity scaling in land surface models

Land surface models (LSMs) are indispensable tools for predicting hydrologic extremes, as well as a particularly uncertain component of Earth system models that has stubbornly resisted the convergence in projections over several successive generations of model intercomparisons. This uncertainty in LSMs is poorly quantified and poorly attributed to specific processes, which has hampered efforts to focus research in reducing uncertainty. This has resulted from sparse sampling of the possible uncertainty space—which is high-dimensional and has contributions from parametric, structural, initial, and boundary condition uncertainties—as an artifact of CMIP-type ensembles of opportunity and limitations inherent in observational benchmarks. A new approach is needed to understand and reduce this uncertainty, based around individual LSMs that can represent the breadth of assumptions represented in current CMIP-type efforts, while at the same time exploring that uncertainty in a systematic way, confronting multiple types of observations, and where justified, replacing process representations with ML-driven emulators. We propose an approach of modular hybrid modeling to address these challenges.

58 GEOSCIENCES↗

Parallel Solver Framework for Mixed-Integer PDE-Constrained Optimization

ROL-PEBBL is a C++, MPI-based parallel code for mixed-integer PDE-constrained optimization (MIPDECO). In these problems we wish to optimize (control, design, etc.) physical systems, which must obey the laws of physics, when some of the decision variables must take integer values. ROL-PEBBL combines a code to efficiently search over integer choices (PEBBL = Parallel Enumeration Branch-and-Bound Library) and a code for efficient nonlinear optimization, including PDE-constrained optimization (ROL = Rapid Optimization Library). In this report, we summarize the design of ROL-PEBBL and initial applications/results. For an artificial source-inversion problem, finding sources of pollution on a grid from sparse samples, ROL-PEBBLs solution for the nest grid gave the best optimization guarantee for any general solver that gives both a solution and a quality guarantee.

97 MATHEMATICS AND COMPUTING↗

stochprop Documentation (Release 1.0)

Simulations of infrasonic propagation in the atmosphere typically utilize a single atmospheric specification describing the acoustic sound speed, ambient winds, and density as a function of altitude. Due to the dynamic and sparsely sampled nature of the atmosphere, there is a notable amount of uncertainty in the atmospheric state at a given location and time so that a more robust analysis of infrasonic propagation requires inclusion of this uncertainty. This Python library, stochprop, has been implemented using methods developed jointly by infrasound scientists at Los Alamos National Laboratory (LANL) and the University of Mississippi’s National Center for Physical Acoustics (NCPA). This software library includes methods to quantify variability in the atmospheric state, identify typical seasonal variability in the atmospheric state and generate suites of representative atmospheric states during a given season, as well as perform uncertainty analysis on a specified atmospheric state given some level of uncertainty. These methods have been designed to interface between propagation modeling capabilities such as InfraGA/GeoAc and NCPAprop and signal analysis methods in the LANL InfraPy tool.

58 GEOSCIENCES↗

Pragmatic Uncertainty Quantification and Propagation in Inverse Estimation of Structural Dynamics Parameters given Material Property Uncertainties and Limited Sensor Data

In this report we demonstrate some relatively simple and inexpensive methods to effectively account for various sources of epistemic lack-of-knowledge type uncertainty in inverse problems. The demonstration problem involves inverse estimation of six parameters of a bolted joint that attaches a kettlebell shaped object to a thick plate. The parameters are efficiently inverted in a modal-based model calibration using gradient-based optimization. Two material properties of the kettlebell are treated as uncertain to within given epistemic uncertainty bounds. We apply and test interval and sparse-sample probabilistic approaches to account for uncertainty in the estimated parameters (and various scalar functionals of the parameters as generic quantities of interest, QOIs) due to uncertainties in the material properties. We also investigate the error effects of limited numbers of vibration sensors (accelerometers) on the kettlebell and plate, and therefore abbreviated excitation/response information in the parameter inversions. We propose and demonstrate a Leave-K-Sensors-Out “cross-prediction” UQ approach to estimate related uncertainties on the parameters and QOI functionals. We indicate how uncertainties from material properties and limited sensors are treated in a combined manner. The economical combined UQ approach involves just three to five samples (i.e. three to five inverse simulations), with no added complication or error/uncertainty from use of surrogate models for affordability. Finally, we describe a related economical UQ approach for handling potential parameter solution non-uniqueness and numerical optimization related precision uncertainties in the estimated parameter values. Indicated further research is identified.

36 MATERIALS SCIENCE↗

Shed Light in the DaRk LineagES of the Fungal Tree of Life—STRES

The polyphyletic group of black fungi within the Ascomycota (Arthoniomycetes, Dothideomycetes, and Eurotiomycetes) is ubiquitous in natural and anthropogenic habitats. Partly because of their dark, melanin-based pigmentation, black fungi are resistant to stresses including UV- and ionizing-radiation, heat and desiccation, toxic metals, and organic pollutants. Consequently, they are amongst the most stunning extremophiles and poly-extreme-tolerant organisms on Earth. Even though ca. 60 black fungal genomes have been sequenced to date, [mostly in the family Herpotrichiellaceae (Eurotiomycetes)], the class Dothideomycetes that hosts the largest majority of extremophiles has only been sparsely sampled. By sequencing up to 92 species that will become reference genomes, the “Shed light in The daRk lineagES of the fungal tree of life” (STRES) project will cover a broad collection of black fungal diversity spread throughout the Fungal Tree of Life. Interestingly, the STRES project will focus on mostly unsampled genera that display different ecologies and life-styles (e.g., ant- and lichen-associated fungi, rock-inhabiting fungi, etc.). With a resequencing strategy of 10- to 15-fold depth coverage of up to ~550 strains, numerous new reference genomes will be established. To identify metabolites and functional processes, these new genomic resources will be enriched with metabolomics analyses coupled with transcriptomics experiments on selected species under various stress conditions (salinity, dryness, UV radiation, oligotrophy). The data acquired will serve as a reference and foundation for establishing an encyclopedic database for fungal metagenomics as well as the biology, evolution, and ecology of the fungi in extreme environments.

59 BASIC BIOLOGICAL SCIENCES↗