Functional error modeling for uncertainty quantification in hydrogeology

Approximate models (proxies) can be employed to reduce the computational costs of estimating uncertainty. The price to pay is that the approximations introduced by the proxy model can lead to a biased estimation. To avoid this problem and ensure a reliable uncertainty quantification, we propose to combine functional data analysis and machine learning to build error models that allow us to obtain an accurate prediction of the exact response without solving the exact model for all realizations. We build the relationship between proxy and exact model on a learning set of geostatistical realizations for which both exact and approximate solvers are run. Functional principal components analysis (FPCA) is used to investigate the variability in the two sets of curves and reduce the dimensionality of the problem while maximizing the retained information. Once obtained, the error model can be used to predict the exact response of any realization on the basis of the sole proxy response. This methodology is purpose-oriented as the error model is constructed directly for the quantity of interest, rather than for the state of the system. Also, the dimensionality reduction performed by FPCA allows a diagnostic of the quality of the error model to assess the informativeness of the learning set and the fidelity of the proxy to the exact model. The possibility of obtaining a prediction of the exact response for any newly generated realization suggests that the methodology can be effectively used beyond the context of uncertainty quantification, in particular for Bayesian inference and optimization.

GOCE: precise orbit determination for the entire mission

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) was the first Earth explorer core mission of the European Space Agency. It was launched on March 17, 2009 into a Sun-synchronous dusk-dawn orbit and re-entered into the Earth’s atmosphere on November 11, 2013. The satellite altitude was between 255 and 225 km for the measurement phases. The European GOCE Gravity consortium is responsible for the Level 1b to Level 2 data processing in the frame of the GOCE High-level processing facility (HPF). The Precise Science Orbit (PSO) is one Level 2 product, which was produced under the responsibility of the Astronomical Institute of the University of Bern within the HPF. This PSO product has been continuously delivered during the entire mission. Regular checks guaranteed a high consistency and quality of the orbits. A correlation between solar activity, GPS data availability and quality of the orbits was found. The accuracy of the kinematic orbit primarily suffers from this. Improvements in modeling the range corrections at the retro-reflector array for the SLR measurements were made and implemented in the independent SLR validation for the GOCE PSO products. The satellite laser ranging (SLR) validation finally states an orbit accuracy of 2.42 cm for the kinematic and 1.84 cm for the reduced-dynamic orbits over the entire mission. The common-mode accelerations from the GOCE gradiometer were not used for the official PSO product, but in addition to the operational HPF work a study was performed to investigate to which extent common-mode accelerations improve the reduced-dynamic orbit determination results. The accelerometer data may be used to derive realistic constraints for the empirical accelerations estimated for the reduced-dynamic orbit determination, which already improves the orbit quality. On top of that the accelerometer data may further improve the orbit quality if realistic constraints and state-of-the-art background models such as gravity field and ocean tide models are used for the reduced-dynamic orbit determination.

High-Spin Molecules: Synthesis, X-ray Characterization, and Magnetic Behavior of Two New Cyano-Bridged Ni II 9 Mo V 6 and Ni II 9 W V 6 Clusters with a S = 12 Ground State

The preparations, X-ray structures, and magnetic characterizations are presented for two new pentadecanuclear cluster compounds:  [NiII{NiII(MeOH)3}8(μ-CN)30{MV(CN)3}6]·xMeOH·yH2O (MV = MoV (1) with x = 17, y = 1; MV = WV (2) with x = 15, y = 0). Both compounds crystallize in the monoclinic space group C2/c, with cell dimensions of a = 28.4957(18) Å, b = 19.2583(10) Å, c = 32.4279(17) Å, β = 113.155(6)°, and Z = 4 for 1 and a = 28.5278(16) Å, b = 19.2008(18) Å, c = 32.4072(17) Å, β = 113.727(6)°, and Z = 4 for 2. The structures of 1 and 2 consist of neutral cluster complexes comprising 15 metal ions, 9 NiII and 6 MV, all linked by μ-cyano ligands. Magnetic susceptibilities and magnetization measurements of compounds 1 and 2 in the crystalline and dissolved state indicate that these clusters have a S = 12 ground state, originating from intracluster ferromagnetic exchange interactions between the μ-cyano-bridged metal ions of the type NiII−NC−MV. Indeed, these data show clearly that the cluster molecules stay intact in solution. Ac magnetic susceptibility measurements reveal that the cluster compounds exhibit magnetic susceptibility relaxation phenomena at low temperatures since, with nonzero dc fields, χ‘ ‘M has a nonzero value that is frequency dependent. However, there appears no out-of-phase (χ‘ ‘M) signal in zero dc field down to 1.8 K, which excludes the expected signature for a single molecule magnet. This finding is confirmed with the small uniaxial magnetic anisotropy value for D of 0.015 cm-1, deduced from the high-field, high-frequency EPR measurement, which distinctly reveals a positive sign in D. Obviously, the overall magnetic anisotropy of the compounds is too low, and this may be a consequence of a small single ion magnetic anisotropy combined with the highly symmetric arrangement of the metal ions in the cluster molecule.

Resolving diffusion in clay minerals at different time scales: Combination of experimental and modeling approaches

Since no single experimental or modeling technique provides data that allow a description of transport processes in clays and clay minerals at all relevant scales, several complementary approaches have to be combined to understand and explain the interplay between transport relevant phenomena. In this paper molecular dynamics simulations (MD) were used to investigate the mobility of water in the interlayer of montmorillonite (Mt), and to estimate the influence of mineral surfaces and interlayer ions on the water diffusion. Random Walk (RW) simulations based on a simplified representation of pore space in Mt were used to estimate and understand the effect of the arrangement of Mt particles on the meso- to macroscopic diffusivity of water. These theoretical calculations were complemented with quasielastic neutron scattering (QENS) measurements of aqueous diffusion in Mt with two pseudo-layers of water performed at four significantly different energy resolutions (i.e. observation times). The size of the interlayer and the size of Mt particles are two characteristic dimensions which determine the time dependent behavior of water diffusion in Mt. MD simulations show that at very short time scales water dynamics has the characteristic features of an oscillatory motion in the cage formed by neighbors in the first coordination shell. At longer time scales, the interaction of water with the surface determines the water dynamics, and the effect of confinement on the overall water mobility within the interlayer becomes evident. At time scales corresponding to an average water displacement equivalent to the average size of Mt particles, the effects of tortuosity are observed in the meso- to macroscopic pore scale simulations. Consistent with the picture obtained in the simulations, the QENS data can be described using a (local) 3D diffusion at short observation times, whereas at sufficiently long observation times a 2D diffusive motion is clearly observed. The effects of tortuosity measured in macroscopic tracer diffusion experiments are in qualitative agreement with RW simulations. By using experimental data to calibrate molecular and mesoscopic theoretical models, a consistent description of water mobility in clay minerals from the molecular to the macroscopic scale can be achieved. In turn, simulations help in choosing optimal conditions for the experimental measurements and the data interpretation. (C) 2014 Elsevier B.V. All rights reserved.

Reevaluation of the MASTER-2009 MLI and H-10 Debris Modeling

The current paper is an excerpt from the doctoral thesis ”Multi-Layer Insulation as Contribution to Orbital Debris”written at the Institute of Aerospace Systems of the Technische Universit ̈at of Braunschweig. The Multi-Layer In-sulation (MLI) population included in ESA’s MASTER-2009 (M eteoroid and Space-Debris Terrestrial Environment Reference) software is based on models for two mechanisms: One model simulates the release of MLI debris during fragmentation events while another estimates the continuo us release of larger MLI pieces due to aging related deterioration of the material. The aim of the thesis was to revise the MLI models from the base up followed by a re-validation of the simulated MLI debris population. The validation is based on comparison to measurement data of the GEO and GTO debris environment obtained by the Astronomical Institute of the University of Bern (AIUB) using ESA’s Space Debris Telescope (ESASDT), the 1-m Zeiss telescope located at the Optical Ground Station (OGS) at the Teide Observatory at Tenerife, Spain. The re-validation led to the conclusion that MLI may cover a much smaller portion of the observed objects than previously published. Further investigation of the resulting discrepancy revealed that the contribution of altogether nine known Ariane H-10 upper stage explosion events which occurred between 1984 and 2002 has very likely been underestimated in past simulations.

Doppler Orbit Determination of Deep Space Probes by the Bernese GNSS Software: First Results of the Combined Orbit Determination from DSN and Inter-Satellite Ka-Band Data from the Grail Mission

Navigation of deep space probes is most commonly operated using the spacecraft Doppler tracking technique. Orbital parameters are determined from a series of repeated measurements of the frequency shift of a microwave carrier over a given integration time. Currently, both ESA and NASA operate antennas at several sites around the world to ensure the tracking of deep space probes. Just a small number of software packages are nowadays used to process Doppler observations. The Astronomical Institute of the University of Bern (AIUB) has recently started the development of Doppler data processing capabilities within the Bernese GNSS Software. This software has been extensively used for Precise Orbit Determination of Earth orbiting satellites using GPS data collected by on-board receivers and for subsequent determination of the Earth gravity field. In this paper, we present the currently achieved status of the Doppler data modeling and orbit determination capabilities in the Bernese GNSS Software using GRAIL data. In particular we will focus on the implemented orbit determination procedure used for the combined analysis of Doppler and intersatellite Ka-band data. We show that even at this earlier stage of the development we can achieve an accuracy of few mHz on two-way S-band Doppler observation and of 2 µm/s on KBRR data from the GRAIL primary mission phase.

Optical Light Curve Observations to Determine Attitude States of Space Debris

The currently proposed space debris remediation measures include the active removal of large objects and “just in time” collision avoidance by deviating the objects using, e.g., ground-based lasers. Both techniques require precise knowledge of the attitude state and state changes of the target objects. In the former case, to devise methods to grapple the target by a tug spacecraft, in the latter, to precisely propagate the orbits of potential collision partners as disturbing forces like air drag and solar radiation pressure depend on the attitude of the objects. Non-resolving optical observations of the magnitude variations, so-called light curves, are a promising technique to determine rotation or tumbling rates and the orientations of the actual rotation axis of objects, as well as their temporal changes. The 1-meter telescope ZIMLAT of the Astronomical Institute of the University of Bern has been used to collect light curves of MEO and GEO objects for a considerable period of time. Recently, light curves of Low Earth Orbit (LEO) targets were acquired as well. We present different observation methods, including active tracking using a CCD subframe readout technique, and the use of a high-speed scientific CMOS camera. Technical challenges when tracking objects with poor orbit redictions, as well as different data reduction methods are addressed. Results from a survey of abandoned rocket upper stages in LEO, examples of abandoned payloads and observations of high area-to-mass ratio debris will be resented. Eventually, first results of the analysis of these light curves are provided.

Space debris attitude simulation - ιOTA (In-Orbit Tumbling Analysis)

Today, there is little knowledge on the attitude state of decommissioned intact objects in Earth orbit. Observational means have advanced in the past years, but are still limited with respect to an accurate estimate of motion vector orientations and magnitude. Especially for the preparation of Active Debris Removal (ADR) missions as planned by ESA’s Clean Space initiative or contingency scenarios for ESA spacecraft like ENVISAT, such knowledge is needed. ESA's “Debris Attitude Motion Measurements and Modelling” project (ESA Contract No. 40000112447), led by the Astronomical Institute of the University of Bern (AIUB), addresses this problem. The goal of the project is to achieve a good understanding of the attitude evolution and the considerable internal and external effects which occur. To characterize the attitude state of selected targets in LEO and GTO, multiple observation methods are combined. Optical observations are carried out by AIUB, Satellite Laser Ranging (SLR) is performed by the Space Research Institute of the Austrian Academy of Sciences (IWF) and radar measurements and signal level determination are provided by the Fraunhofer Institute for High Frequency Physics and Radar Techniques (FHR). The In-Orbit Tumbling Analysis tool (ιOTA) is a prototype software, currently in development by Hyperschall Technologie Göttingen GmbH (HTG) within the framework of the project. ιOTA will be a highly modular software tool to perform short-(days), medium-(months) and long-term (years) propagation of the orbit and attitude motion (six degrees-of-freedom) of spacecraft in Earth orbit. The simulation takes into account all relevant acting forces and torques, including aerodynamic drag, solar radiation pressure, gravitational influences of Earth, Sun and Moon, eddy current damping, impulse and momentum transfer from space debris or micro meteoroid impact, as well as the optional definition of particular spacecraft specific influences like tank sloshing, reaction wheel behaviour, magnetic torquer activity and thruster firing. The purpose of ιOTA is to provide high accuracy short-term simulations to support observers and potential ADR missions, as well as medium-and long-term simulations to study the significance of the particular internal and external influences on the attitude, especially damping factors and momentum transfer. The simulation will also enable the investigation of the altitude dependency of the particular external influences. ιOTA's post-processing modules will generate synthetic measurements for observers and for software validation. The validation of the software will be done by cross-calibration with observations and measurements acquired by the project partners.

Hydrogen atom in space with a compactified extra dimension and potential defined by Gauss' law

We investigate the consequences of one extra spatial dimension for the stability and energy spectrum of the non-relativistic hydrogen atom with a potential defined by Gauss' law, i.e. proportional to 1 /| x | 2 . The additional spatial dimension is considered to be either infinite or curled-up in a circle of radius R. In both cases, the energy spectrum is bounded from below for charges smaller than the same critical value and unbounded from below otherwise. As a consequence of compactification, negative energy eigenstates appear: if R is smaller than a quarter of the Bohr radius, the corresponding Hamiltonian possesses an infinite number of bound states with minimal energy extending at least to the ground state of the hydrogen atom.

Keep at bay! - Abnormal personal space regulation as marker of paranoia in schizophrenia.

BACKGROUND During threat, interpersonal distance is deliberately increased. Personal space regulation is related to amygdala function and altered in schizophrenia, but it remains unknown whether it is particularly associated with paranoid threat. METHODS We compared performance in two tests on personal space between 64 patients with schizophrenia spectrum disorders and 24 matched controls. Patients were stratified in those with paranoid threat, neutral affect or paranoid experience of power. In the stop-distance paradigm, participants indicated the minimum tolerable interpersonal distance. In the fixed-distance paradigm, they indicated the level of comfort at fixed interpersonal distances. RESULTS Paranoid threat increased interpersonal distance two-fold in the stop-distance paradigm, and reduced comfort ratings in the fixed-distance paradigm. In contrast, patients experiencing paranoid power had high comfort ratings at any distance. Patients with neutral affect did not differ from controls in the stop-distance paradigm. Differences between groups remained when controlling for gender and positive symptom severity. Among schizophrenia patients, the stop-distance paradigm detected paranoid threat with 93% sensitivity and 83% specificity. CONCLUSIONS Personal space regulation is not generally altered in schizophrenia. However, state paranoid experience has distinct contributions to personal space regulation. Subjects experiencing current paranoid threat share increased safety-seeking behavior.

Direct visualization of the outer membrane of mycobacteria and corynebacteria in their native state.

The cell envelope of mycobacteria, which include the causative agents of tuberculosis and leprosy, is crucial for their success as pathogens. Despite a continued strong emphasis on identifying the multiple chemical components of this envelope, it has proven difficult to combine its components into a comprehensive structural model, primarily because the available ultrastructural data rely on conventional electron microscopy embedding and sectioning, which are known to induce artifacts. The existence of an outer membrane bilayer has long been postulated but has never been directly observed by electron microscopy of ultrathin sections. Here we have used cryo-electron microscopy of vitreous sections (CEMOVIS) to perform a detailed ultrastructural analysis of three species belonging to the Corynebacterineae suborder, namely, Mycobacterium bovis BCG, Mycobacterium smegmatis, and Corynebacterium glutamicum, in their native state. We provide new information that accurately describes the different layers of the mycobacterial cell envelope and challenges current models of the organization of its components. We show a direct visualization of an outer membrane, analogous to that found in gram-negative bacteria, in the three bacterial species examined. Furthermore, we demonstrate that mycolic acids, the hallmark of mycobacteria and related genera, are essential for the formation of this outer membrane. In addition, a granular layer and a low-density zone typifying the periplasmic space of gram-positive bacteria are apparent in CEMOVIS images of mycobacteria and corynebacteria. Based on our observations, a model of the organization of the lipids in the outer membrane is proposed. The architecture we describe should serve as a reference for future studies to relate the structure of the mycobacterial cell envelope to its function.

Modeling the Histidine–Phenylalanine Interaction: The NH···π Hydrogen Bond of Imidazole·Benzene

NH···π hydrogen bonds occur frequently between the amino acid side groups in proteins and peptides. Data-mining studies of protein crystals find that ~80% of the T-shaped histidine···aromatic contacts are CH···π, and only ~20% are NH···π interactions. We investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet-cooled imidazole·benzene (Im·Bz) complex as a model for the NH···π interaction between histidine and phenylalanine. Ground- and excited-state dispersion-corrected density functional calculations and correlated methods (SCS-MP2 and SCS-CC2) predict that Im·Bz has a Cs-symmetric T-shaped minimum-energy structure with an NH···π hydrogen bond to the Bz ring; the NH bond is tilted 12° away from the Bz C₆ axis. IR depletion spectra support the T-shaped geometry: The NH stretch vibrational fundamental is red shifted by −73 cm⁻¹ relative to that of bare imidazole at 3518 cm⁻¹, indicating a moderately strong NH···π interaction. While the Sₒ(A1g) → S₁(B₂u) origin of benzene at 38 086 cm⁻¹ is forbidden in the gas phase, Im·Bz exhibits a moderately intense Sₒ → S₁ origin, which appears via the D₆h → Cs symmetry lowering of Bz by its interaction with imidazole. The NH···π ground-state hydrogen bond is strong, De=22.7 kJ/mol (1899 cm⁻¹). The combination of gas-phase UV and IR spectra confirms the theoretical predictions that the optimum Im·Bz geometry is T shaped and NH···π hydrogen bonded. We find no experimental evidence for a CH···π hydrogen-bonded ground-state isomer of Im·Bz. The optimum NH···π geometry of the Im·Bz complex is very different from the majority of the histidine·aromatic contact geometries found in protein database analyses, implying that the CH···π contacts observed in these searches do not arise from favorable binding interactions but merely from protein side-chain folding and crystal-packing constraints. The UV and IR spectra of the imidazole·(benzene)₂ cluster are observed via fragmentation into the Im·Bz+ mass channel. The spectra of Im·Bz and Im·Bz₂ are cleanly separable by IR hole burning. The UV spectrum of Im·Bz₂ exhibits two 000 bands corresponding to the Sₒ → S₁ excitations of the two inequivalent benzenes, which are symmetrically shifted by −86/+88 cm⁻¹ relative to the 000 band of benzene.

3D-modeling of Mercury's solar wind sputtered surface-exosphere environment

The efficiency of sputtered refractory elements by H+ and He++ solar wind ions from Mercury's surface and their contribution to the exosphere are studied for various solar wind conditions. A 3D solar wind-planetary interaction hybrid model is used for the evaluation of precipitation maps of the sputter agents on Mercury's surface. By assuming a global mineralogical surface composition, the related sputter yields are calculated by means of the 2013 SRIM code and are coupled with a 3D exosphere model. Because of Mercury's magnetic field, for quiet and nominal solar wind conditions the plasma can only precipitate around the polar areas, while for extreme solar events (fast solar wind, coronal mass ejections, interplanetary magnetic clouds) the solar wind plasma has access to the entire dayside. In that case the release of particles form the planet's surface can result in an exosphere density increase of more than one order of magnitude. The corresponding escape rates are also about an order of magnitude higher. Moreover, the amount of He++ ions in the precipitating solar plasma flow enhances also the release of sputtered elements from the surface in the exosphere. A comparison of our model results with MESSENGER observations of sputtered Mg and Ca elements in the exosphere shows a reasonable quantitative agreement. (C) 2015 Elsevier Ltd. All rights reserved.

Modeling asteroid collisions and impact processes

As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.

Early Life History and a Modeling Framework for Lobster (Homarus Americanus) Populations in the Gulf of Maine

Beginning in the late 1980s, lobster (Homarus americanus) landings for the state of Maine and the Bay of Fundy increased to levels more than three times their previous 20-year means. Reduced predation may have permitted the expansion of lobsters into previously inhospitable territory, but we argue that in this region the spatial patterns of recruitment and the abundance of lobsters are substantially driven by events governing the earliest life history stages, including the abundance and distribution of planktonic stages and their initial settlement as Young-of-Year (YOY) lobsters. Settlement densities appear to be strongly driven by abundance of the pelagic postlarvae. Postlarvae and YOY show large-scale spatial patterns commensurate with coastal circulation, but also multi-year trends in abundance and abrupt shifts in abundance and spatial patterns that signal strong environmental forcing. The extent of the coastal shelf that defines the initial settlement grounds for lobsters is important to future population modeling. We address one part of this definition by examining patterns of settlement with depth, and discuss a modeling framework for the full life history of lobsters in the Gulf of Maine.