Robust Uncertainty Quantification and Scalable Computation for Computer Models with Massive Output

Uncertainty quantification (UQ) is both an old and new concept. The current novelty lies in the interactions and synthesis of mathematical models, computer experiments, statistics, field/real experiments, and probability theory, with a particular emphasize on the large-scale simulations by computer models. The challenges not only come from the complication of scientific questions, but also from the size of the information. It is the focus in this thesis to provide statistical models that are scalable to massive data produced in computer experiments and real experiments, through fast and robust statistical inference.

Chapter 2 provides a practical approach for simultaneously emulating/approximating massive number of functions, with the application on hazard quantification of Soufri\`{e}re Hills volcano in Montserrate island. Chapter 3 discusses another problem with massive data, in which the number of observations of a function is large. An exact algorithm that is linear in time is developed for the problem of interpolation of Methylation levels. Chapter 4 and Chapter 5 are both about the robust inference of the models. Chapter 4 provides a new criteria robustness parameter estimation criteria and several ways of inference have been shown to satisfy such criteria. Chapter 5 develops a new prior that satisfies some more criteria and is thus proposed to use in practice.

Sea-bed photographs (benthos) from the AWI-Hausgarten area along OFOS profile PS62/191-1

Although the use of deep-sea imagery considerably increased during the last decades, reports on nekton falls to the deep seafloor are very scarce. Whereas there are a few reports describing the finding of whale carcasses in the deep north-eastern and south-eastern Pacific, descriptions of invertebrate or vertebrate food-falls at centimetre to metre scale are extremely rare. After 4 years of extensive work at a deep-sea long-term station in northern polar regions (AWI-"Hausgarten"), including large-scale visual observations with various camera systems covering some 10 000 m2 of seafloor at water depths between 1250 and 5600 m, this paper describes the first observation of a fish carcass at about 1280 m water depth, west off Svålbard. The fish skeleton had a total length of 36 cm and an approximated biomass of 0.5 kg wet weight. On the basis of in situ experiments, we estimated a very short residence time of this particular carcass of about 7 h at the bottom. The fast response of the motile deep-sea scavenger community to such events and the rapid utilisation of this kind of organic carbon supply might partly explain the extreme rarity of such an observation.

Data collection of Calanus finmarchicus reproduction life history traits in the North Atlantic Ocean

Observations of egg production rates (EPR) for female Calanus finmarchicus were compared from different regions of the North Atlantic. The regions were diverse in size and sampling frequency, ranging from a fixed time series station in the Lower St Lawrence Estuary, off Rimouski, where nearly 200 experiments were carried out between May and December from 1994 to 2006, to a large-scale survey in the Northern Norwegian Sea, where about 50 experiments were carried out between April and June from 2002 to 2004. For this analysis the stations were grouped mostly along geographic lines, with only limited attention being paid to oceanographic features. There is some overlap between regions, however, where stations were sometimes kept together when they were sampled on the same cruise. As well some stations other than off Rimouski were occupied more than once during different years and/or in different seasons.

A suite of global, cross-scale topographic variables for environmental and biodiversity modeling, links to files in GeoTIFF format

Topographic variation, the spatial variation in elevation and terrain features, underpins a myriad of patterns and processes in geography and ecology and is key to understanding the variation of life on the planet. The characterization of this variation is scale-dependent, i.e. it varies with the distance over which features are assessed and with the spatial grain (grid cell resolution) of analysis. A fully standardized and global multivariate product of different terrain features has the potential to support many large-scale basic research and analytical applications, however to date, such technique is unavailable. Here we used the digital elevation model products of global 250 m GMTED and near-global 90 m SRTM to derive a suite of topographic variables: elevation, slope, aspect, eastness, northness, roughness, terrain roughness index, topographic position index, vector ruggedness measure, profile and tangential curvature, and 10 geomorphological landform classes. We aggregated each variable to 1, 5, 10, 50 and 100 km spatial grains using several aggregation approaches (median, average, minimum, maximum, standard deviation, percent cover, count, majority, Shannon Index, entropy, uniformity). While a global cross-correlation underlines the high similarity of many variables, a more detailed view in four mountain regions reveals local differences, as well as scale variations in the aggregated variables at different spatial grains. All newly-developed variables are available for download at http://www.earthenv.org and can serve as a basis for standardized hydrological, environmental and biodiversity modeling at a global extent.

UTILIZATION OF GIS AND SPATIAL ANALYSIS TECHNIQUES IN ARTISANAL AND SMALL SCALE MINING TO LOCATE A CENTRALIZED PROCESSING CENTRE

One of the global phenomena with threats to environmental health and safety is artisanal mining. There are ambiguities in the manner in which an ore-processing facility operates which hinders the mining capacity of these miners in Ghana. These problems are reviewed on the basis of current socio-economic, health and safety, environmental, and use of rudimentary technologies which limits fair-trade deals to miners. This research sought to use an established data-driven, geographic information (GIS)-based system employing the spatial analysis approach for locating a centralized processing facility within the Wassa Amenfi-Prestea Mining Area (WAPMA) in the Western region of Ghana. A spatial analysis technique that utilizes ModelBuilder within the ArcGIS geoprocessing environment through suitability modeling will systematically and simultaneously analyze a geographical dataset of selected criteria. The spatial overlay analysis methodology and the multi-criteria decision analysis approach were selected to identify the most preferred locations to site a processing facility. For an optimal site selection, seven major criteria including proximity to settlements, water resources, artisanal mining sites, roads, railways, tectonic zones, and slopes were considered to establish a suitable location for a processing facility. Site characterizations and environmental considerations, incorporating identified constraints such as proximity to large scale mines, forest reserves and state lands to site an appropriate position were selected. The analysis was limited to criteria that were selected and relevant to the area under investigation. Saaty’s analytical hierarchy process was utilized to derive relative importance weights of the criteria and then a weighted linear combination technique was applied to combine the factors for determination of the degree of potential site suitability. The final map output indicates estimated potential sites identified for the establishment of a facility centre. The results obtained provide intuitive areas suitable for consideration

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the Iberia centre fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. These anisotropies have been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

Linking environmental variables with regional- scale variability in ecological structure and standing stock of carbon within UK kelp forests

Kelp forests represent some of the most productive and diverse habitats on Earth. Understanding drivers of ecological patterns at large spatial scales is critical for effective management and conservation of marine habitats. We surveyed kelp forests dominated by Laminaria hyperborea (Gunnerus) Foslie 1884 across 9° latitude and >1000 km of coastline and measured a number of physical parameters at multiple scales to link ecological structure and standing stock of carbon with environmental variables. Kelp density, biomass, morphology and age were generally greater in exposed sites within regions, highlighting the importance of wave exposure in structuring L. hyperborea populations. At the regional scale, wave-exposed kelp canopies in the cooler regions (the north and west of Scotland) were greater in biomass, height and age than in warmer regions (southwest Wales and England). The range and maximal values of estimated standing stock of carbon contained within kelp forests was greater than in historical studies, suggesting that this ecosystem property may have been previously undervalued. Kelp canopy density was positively correlated with large-scale wave fetch and fine-scale water motion, whereas kelp canopy biomass and the standing stock of carbon were positively correlated with large-scale wave fetch and light levels and negatively correlated with temperature. As light availability and summer temperature were important drivers of kelp forest biomass, effective management of human activities that may affect coastal water quality is necessary to maintain ecosystem functioning, while increased temperatures related to anthropogenic climate change may impact the structure of kelp forests and the ecosystem services they provide.

Linking environmental variables with regional- scale variability in ecological structure and standing stock of carbon within UK kelp forests

Kelp forests represent some of the most productive and diverse habitats on Earth. Understanding drivers of ecological patterns at large spatial scales is critical for effective management and conservation of marine habitats. We surveyed kelp forests dominated by Laminaria hyperborea (Gunnerus) Foslie 1884 across 9° latitude and >1000 km of coastline and measured a number of physical parameters at multiple scales to link ecological structure and standing stock of carbon with environmental variables. Kelp density, biomass, morphology and age were generally greater in exposed sites within regions, highlighting the importance of wave exposure in structuring L. hyperborea populations. At the regional scale, wave-exposed kelp canopies in the cooler regions (the north and west of Scotland) were greater in biomass, height and age than in warmer regions (southwest Wales and England). The range and maximal values of estimated standing stock of carbon contained within kelp forests was greater than in historical studies, suggesting that this ecosystem property may have been previously undervalued. Kelp canopy density was positively correlated with large-scale wave fetch and fine-scale water motion, whereas kelp canopy biomass and the standing stock of carbon were positively correlated with large-scale wave fetch and light levels and negatively correlated with temperature. As light availability and summer temperature were important drivers of kelp forest biomass, effective management of human activities that may affect coastal water quality is necessary to maintain ecosystem functioning, while increased temperatures related to anthropogenic climate change may impact the structure of kelp forests and the ecosystem services they provide.

Recent advances in laser-driven neutron sources

Due to the limited number and high cost of large-scale neutron facilities, there has been a growing interest in compact accelerator-driven sources. In this context, several potential schemes of laser-driven neutron sources are being intensively studied employing laser-accelerated electron and ion beams. In addition to the potential of delivering neutron beams with high brilliance, directionality and ultra-short burst duration, a laser-driven neutron source would offer further advantages in terms of cost-effectiveness, compactness and radiation confinement by closed-coupled experiments. Some of the recent advances in this field are discussed,
showing improvements in the directionality and flux of the laser-driven neutron beams.

Modeling long-term variability and change of soil moisture and groundwater level - from catchment to global scale

The water stored in and flowing through the subsurface is fundamental for sustaining human activities and needs, feeding water and its constituents to surface water bodies and supporting the functioning of their ecosystems. Quantifying the changes that affect the subsurface water is crucial for our understanding of its dynamics and changes driven by climate change and other changes in the landscape, such as in land-use and water-use. It is inherently difficult to directly measure soil moisture and groundwater levels over large spatial scales and long times. Models are therefore needed to capture the soil moisture and groundwater level dynamics over such large spatiotemporal scales. This thesis develops a modeling framework that allows for long-term catchment-scale screening of soil moisture and groundwater level changes. The novelty in this development resides in an explicit link drawn between catchment-scale hydroclimatic and soil hydraulics conditions, using observed runoff data as an approximation of soil water flux and accounting for the effects of snow storage-melting dynamics on that flux. Both past and future relative changes can be assessed by use of this modeling framework, with future change projections based on common climate model outputs. By direct model-observation comparison, the thesis shows that the developed modeling framework can reproduce the temporal variability of large-scale changes in soil water storage, as obtained from the GRACE satellite product, for most of 25 large study catchments around the world. Also compared with locally measured soil water content and groundwater level in 10 U.S. catchments, the modeling approach can reasonably well reproduce relative seasonal fluctuations around long-term average values. The developed modeling framework is further used to project soil moisture changes due to expected future climate change for 81 catchments around the world. The future soil moisture changes depend on the considered radiative forcing scenario (RCP) but are overall large for the occurrence frequency of dry and wet events and the inter-annual variability of seasonal soil moisture. These changes tend to be higher for the dry events and the dry season, respectively, than for the corresponding wet quantities, indicating increased drought risk for some parts of the world.

Konzeption, Implementierung und Anwendung eines automatisierten aerothermodynamischen Vorentwurfsprozesses für Axialturbinen

In dieser Arbeit wird ein Prozess für den frühen aerothermodynamischen Entwurf von Axialturbinen konzipiert und durch Kopplung einzelner Computerprogramme im DLR Göttingen realisiert. Speziell für die Erstauslegung von Geometrien und die Vorhersage von globalen Leistungsdaten beliebiger Axialturbinen wurde ein neues Programm erzeugt. Dessen effiziente Anwendung wird mit einer zu diesem Zweck konzipierten grafischen Entwurfsumgebung ausgeführt. Kennzeichnend für den Vorentwurfsprozess in dieser Arbeit ist die Anwendung von ein- und zweidimensionaler Strömungssimulation sowie der hohe Grad an Verknüpfung der verwendeten Programme sowohl auf prozesstechnischer wie auch auf datentechnischer Ebene. Dabei soll dem sehr frühen Entwurf eine deutlich stärkere Rolle zukommen als bisher üblich und im Gegenzug die Entwurfszeit mit höher auflösenden Vorentwurfsprogrammen reduziert werden. Die Anwendung der einzelnen Programme im Rahmen von Subprozessen wird anhand von exemplarischen Turbinenkonfigurationen in der Arbeit ebenso dargestellt, wie die Validierung des gesamten Entwurfsprozesses anhand der Auslegung einer folgend realisierten und erfolgreich operierenden Axialturbine eines Triebwerkssimulators für Flugzeug-Windkanalmodelle (TPS). Neben der Erleichterung von manueller Entwurfstätigkeit durch grafische Benutzerinteraktion kommt in einzelnen Subprozessen eine automatisierte Mehrziel-Optimierung zum Einsatz.

Impacts des dépositions atmosphériques de fer sur les assemblages phytoplanctoniques et la production de diméthylsulfure dans le Pacifique Nord-Est contemporain et dans le contexte de l’acidification des océans

La croissance du phytoplancton est limitée par les faibles concentrations de fer (Fe) dans près de 40% de l’océan mondial. Le Pacifique subarctique Nord-Est représente une de ces zones limitées en fer et désignées High Nutrient - Low Chlorophyll (HNLC). Cet écosystème, dominé par des cellules de petite taille telles les prymnésiophytes, est caractérisé par de très faibles concentrations estivales de chlorophylle a et de fortes concentrations de macronutriments. Il a été maintes fois démontré que les ajouts de fer, sous différentes formes chimiques (habituellement FeSO4), dans les zones HNLC, stimulent la croissance et modifient la structure des communautés planctoniques en favorisant la croissance des cellules de grande taille, notamment les diatomées. Ces effets sur la communauté planctonique ont le potentiel d’influencer les grands mécanismes régulateurs du climat, tels la pompe biologique de carbone et la production de diméthylsulfure (DMS). Les poussières provenant des déserts du nord de la Chine sont reconnues depuis longtemps comme une source sporadique importante de fer pour le Pacifique Nord-Est. Malgré leur importance potentielle, l’influence directe exercée par ces poussières sur l’écosystème planctonique de cette zone HNLC n’a jamais été étudiée. Il s’agit d’une lacune importante puisque le fer associé aux poussières est peu soluble dans l’eau de mer, que la proportion biodisponible n’est pas connue et que les poussières peuvent avoir un effet inhibiteur chez le phytoplancton. Cette thèse propose donc, dans un premier temps, de mesurer pour la première fois l’effet de la fertilisation de la communauté planctonique du Pacifique Nord-Est par un gradient de concentrations de poussières désertiques naturelles. Cette première expérimentation a démontré que le fer contenu dans les poussières asiatiques est biodisponible et qu’une déposition équivalente à celles prenant place au printemps dans le Pacifique Nord-Est peut résulter en une stimulation significative de la prise de nutriments et de la croissance du phytoplancton. Mes travaux ont également montré que l’ajout de 0,5 mg L-1 de poussières peut résulter en la production d’autant de biomasse algale que l’ajout de FeSO4, l’espèce chimique utilisée lors des expériences d’enrichissement en fer à grande échelle. Cependant, les ajouts de FeSO4 favorisent davantage les cellules de petite taille que les ajouts de poussières, observation démontrant que le FeSO4 n’est pas un proxy parfait des poussières asiatiques. Dans un deuxième temps, je me suis intéressée à une source alternative de fer atmosphérique, les cendres volcaniques. Mon intérêt pour cette source de fer a été attisé par les observations d’une floraison spectaculaire dans le Pacifique Nord-Est, ma région d’étude, associée à l’éruption de 2008 du volcan Kasatochi dans les îles Aléoutiennes. Forte de mon expérience sur les poussières, j’ai quantifié l’effet direct de ces cendres volcaniques sur la communauté planctonique du Pacifique Nord-Est. Mes résultats ont montré que le fer contenu dans les cendres volcaniques est également biodisponible pour le phytoplancton. Ils ont également montré que cette source de fer peut être aussi importante que les poussières désertiques dans la régulation de la croissance du phytoplancton dans cette partie de l’océan global à l’échelle millénaire. Dans un troisième temps, j’ai estimé comment l’acidification des océans modulera les réponses des communautés planctoniques aux dépositions naturelles de fer mises en évidence lors de mes expériences précédentes. Pour ce faire, j’ai effectué des enrichissements de poussière dans de l’eau de mer au pH actuel de 8.0 et dans l’eau de mer acidifiée à un pH de 7.8. Mes résultats ont montré une diminution du taux de croissance du phytoplancton dans le milieu acidifié mais pas de changement notable dans la structure de la communauté. Les ajouts de poussières et de cendres, de même que les variations de pH, n’ont pas eu d’effet significatif sur la production de DMS et de son précurseur le diméthylsulfoniopropionate (DMSP), probablement en raison de la courte durée (4 jours) des expériences. L’ensemble des résultats de cette thèse montre que le fer contenu dans diverses sources atmosphériques naturelles est biodisponible pour le phytoplancton du Pacifique Nord-Est et que des taux de déposition réalistes peuvent stimuler la croissance de manière notable dans les premiers jours suivant une tempête désertique ou une éruption volcanique. Finalement, les résultats de mes expériences à stresseurs multiples Fer/acidification suggèrent une certaine résistance des communautés phytoplanctoniques à la diminution du pH prédite d’ici la fin du siècle pour les eaux de surface des océans.

Flexible modeling and estimation for high-dimensional graphs

Thesis (Ph.D.)--University of Washington, 2016-08

From Aari to Zulu: Massively Multilingual Creation of Language Tools using Interlinear Glossed Text

Thesis (Ph.D.)--University of Washington, 2016-08

A large study on the effect of code obfuscation on the quality of java code

Context: Obfuscation is a common technique used to protect software against mali- cious reverse engineering. Obfuscators manipulate the source code to make it harder to analyze and more difficult to understand for the attacker. Although different ob- fuscation algorithms and implementations are available, they have never been directly compared in a large scale study. Aim: This paper aims at evaluating and quantifying the effect of several different obfuscation implementations (both open source and commercial), to help developers and project manager to decide which one could be adopted. Method: In this study we applied 44 obfuscations to 18 subject applications covering a total of 4 millions lines of code. The effectiveness of these source code obfuscations has been measured using 10 code metrics, considering modularity, size and complexity of code. Results: Results show that some of the considered obfuscations are effective in mak- ing code metrics change substantially from original to obfuscated code, although this change (called potency of the obfuscation) is different on different metrics. In the pa- per we recommend which obfuscations to select, given the security requirements of the software to be protected.