Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Tumor necrosis factor (TNF) signaling, but not TWEAK (TNF-like weak inducer of apoptosis)-triggered cIAP1 (cellular inhibitor of apoptosis protein 1) degradation, requires cIAP1 RING dimerization and E2 binding.

Cellular inhibitor of apoptosis (cIAP) proteins, cIAP1 and cIAP2, are important regulators of tumor necrosis factor (TNF) superfamily (SF) signaling and are amplified in a number of tumor types. They are targeted by IAP antagonist compounds that are undergoing clinical trials. IAP antagonist compounds trigger cIAP autoubiquitylation and degradation. The TNFSF member TWEAK induces lysosomal degradation of TRAF2 and cIAPs, leading to elevated NIK levels and activation of non-canonical NF-kappaB. To investigate the role of the ubiquitin ligase RING domain of cIAP1 in these pathways, we used cIAP-deleted cells reconstituted with cIAP1 point mutants designed to interfere with the ability of the RING to dimerize or to interact with E2 enzymes. We show that RING dimerization and E2 binding are required for IAP antagonists to induce cIAP1 degradation and protect cells from TNF-induced cell death. The RING functions of cIAP1 are required for full TNF-induced activation of NF-kappaB, however, delayed activation of NF-kappaB still occurs in cIAP1 and -2 double knock-out cells. The RING functions of cIAP1 are also required to prevent constitutive activation of non-canonical NF-kappaB by targeting NIK for proteasomal degradation. However, in cIAP double knock-out cells TWEAK was still able to increase NIK levels demonstrating that NIK can be regulated by cIAP-independent pathways. Finally we show that, unlike IAP antagonists, TWEAK was able to induce degradation of cIAP1 RING mutants. These results emphasize the critical importance of the RING of cIAP1 in many signaling scenarios, but also demonstrate that in some pathways RING functions are not required.

Which prior knowledge? Quantification of in vivo brain 13C MR spectra following 13C glucose infusion using AMARES.

The recent developments in high magnetic field 13C magnetic resonance spectroscopy with improved localization and shimming techniques have led to important gains in sensitivity and spectral resolution of 13C in vivo spectra in the rodent brain, enabling the separation of several 13C isotopomers of glutamate and glutamine. In this context, the assumptions used in spectral quantification might have a significant impact on the determination of the 13C concentrations and the related metabolic fluxes. In this study, the time domain spectral quantification algorithm AMARES (advanced method for accurate, robust and efficient spectral fitting) was applied to 13 C magnetic resonance spectroscopy spectra acquired in the rat brain at 9.4 T, following infusion of [1,6-(13)C2 ] glucose. Using both Monte Carlo simulations and in vivo data, the goal of this work was: (1) to validate the quantification of in vivo 13C isotopomers using AMARES; (2) to assess the impact of the prior knowledge on the quantification of in vivo 13C isotopomers using AMARES; (3) to compare AMARES and LCModel (linear combination of model spectra) for the quantification of in vivo 13C spectra. AMARES led to accurate and reliable 13C spectral quantification similar to those obtained using LCModel, when the frequency shifts, J-coupling constants and phase patterns of the different 13C isotopomers were included as prior knowledge in the analysis.

Sparse regularization for fiber ODF reconstruction: from the suboptimality of ℓ2 and ℓ1 priors to ℓ0.

Diffusion MRI is a well established imaging modality providing a powerful way to probe the structure of the white matter non-invasively. Despite its potential, the intrinsic long scan times of these sequences have hampered their use in clinical practice. For this reason, a large variety of methods have been recently proposed to shorten the acquisition times. Among them, spherical deconvolution approaches have gained a lot of interest for their ability to reliably recover the intra-voxel fiber configuration with a relatively small number of data samples. To overcome the intrinsic instabilities of deconvolution, these methods use regularization schemes generally based on the assumption that the fiber orientation distribution (FOD) to be recovered in each voxel is sparse. The well known Constrained Spherical Deconvolution (CSD) approach resorts to Tikhonov regularization, based on an ℓ(2)-norm prior, which promotes a weak version of sparsity. Also, in the last few years compressed sensing has been advocated to further accelerate the acquisitions and ℓ(1)-norm minimization is generally employed as a means to promote sparsity in the recovered FODs. In this paper, we provide evidence that the use of an ℓ(1)-norm prior to regularize this class of problems is somewhat inconsistent with the fact that the fiber compartments all sum up to unity. To overcome this ℓ(1) inconsistency while simultaneously exploiting sparsity more optimally than through an ℓ(2) prior, we reformulate the reconstruction problem as a constrained formulation between a data term and a sparsity prior consisting in an explicit bound on the ℓ(0)norm of the FOD, i.e. on the number of fibers. The method has been tested both on synthetic and real data. Experimental results show that the proposed ℓ(0) formulation significantly reduces modeling errors compared to the state-of-the-art ℓ(2) and ℓ(1) regularization approaches.

A Modulated Closed Form solution for Quantitative Susceptibility Mapping - A thorough evaluation and comparison to iterative methods based on edge prior knowledge.

The aim of this study is to perform a thorough comparison of quantitative susceptibility mapping (QSM) techniques and their dependence on the assumptions made. The compared methodologies were: two iterative single orientation methodologies minimizing the l2, l1TV norm of the prior knowledge of the edges of the object, one over-determined multiple orientation method (COSMOS) and anewly proposed modulated closed-form solution (MCF). The performance of these methods was compared using a numerical phantom and in-vivo high resolution (0.65mm isotropic) brain data acquired at 7T using a new coil combination method. For all QSM methods, the relevant regularization and prior-knowledge parameters were systematically changed in order to evaluate the optimal reconstruction in the presence and absence of a ground truth. Additionally, the QSM contrast was compared to conventional gradient recalled echo (GRE) magnitude and R2* maps obtained from the same dataset. The QSM reconstruction results of the single orientation methods show comparable performance. The MCF method has the highest correlation (corrMCF=0.95, r(2)MCF =0.97) with the state of the art method (COSMOS) with additional advantage of extreme fast computation time. The l-curve method gave the visually most satisfactory balance between reduction of streaking artifacts and over-regularization with the latter being overemphasized when the using the COSMOS susceptibility maps as ground-truth. R2* and susceptibility maps, when calculated from the same datasets, although based on distinct features of the data, have a comparable ability to distinguish deep gray matter structures.

A forensic perspective of the AFL investigation into peptides: an antidoping investigation case study.

BACKGROUND: The World Anti-Doping Agency (WADA) is introducing enhancements to doping investigations in its 2015 Code, which include improved sharing of information between antidoping organisations (including sporting bodies) and enhanced accountability of athlete support staff. These additions will improve the control of links between sports doping and organised crime. In February 2013 the Australian Crime Commission released a report that linked several professional sporting codes, professional athletes with links to organised crime, performance enhancing drugs and illicit substances. Following this report the Australian Football League (AFL) partnered the Australian national antidoping organisation to investigate peptide use in Australian football. METHODS: This review compared the model proposed by Marclay, a hypothetical model for anti-doping investigations that proposed a forensic intelligence and analysis approach, to use the forensic capabilities of the AFL investigation to test the model's relevance to an actual case. RESULTS: The investigation uncovered the use of peptides used to enhance athlete performance. The AFL investigation found a high risk of doping where athlete support staff existed in teams with weak corporate governance controls. A further finding included the need for the investigation to provide a timely response in professional team sports that were sensitive to the competition timing. In the case of the AFL the team was sanctioned prior to the finals as an interim outcome for allowing the risk of use of performance-enhancing substances. Doping violation charges are still being considered. DISCUSSION: Antidoping strategies should include the investigation of corporate officers in team doping circumstances, the mandatory recording of all athlete substance use during competition and training phases, the wider sharing of forensic intelligence with non-sporting bodies particularly law enforcement and collaboration between antidoping and sporting organisations in doping investigations. CONCLUSIONS: The AFL investigation illustrated the importance of the 2015 WADA Code changes and highlighted the need for a systematic use of broad forensic intelligence activities in the investigation of doping violations.

Comparison of quantiferon-TB gold with tuberculin skin test for detecting latent tuberculosis infection prior to liver transplantation.

Screening for latent tuberculosis infection (LTBI) is recommended prior to organ transplantation. The Quantiferon-TB Gold assay (QFT-G) may be more accurate than the tuberculin skin test (TST) in the detection of LTBI. We prospectively compared the results of QFT-G to TST in patients with chronic liver disease awaiting transplantation. Patients were screened for LTBI with both the QFT-G test and a TST. Concordance between test results and predictors of a discordant result were determined. Of the 153 evaluable patients, 37 (24.2%) had a positive TST and 34 (22.2%) had a positive QFT-G. Overall agreement between tests was 85.1% (kappa= 0.60, p < 0.0001). Discordant test results were seen in 12 TST positive/QFT-G negative patients and in 9 TST negative/QFT-G positive patients. Prior BCG vaccination was not associated with discordant test results. Twelve patients (7.8%), all with a negative TST, had an indeterminate result of the QFT-G and this was more likely in patients with a low lymphocyte count (p = 0.01) and a high MELD score (p = 0.001). In patients awaiting liver transplantation, both the TST and QFT-G were comparable for the diagnosis of LTBI with reasonable concordance between tests. Indeterminate QFT-G result was more likely in those with more advanced liver disease.

The robustness of the weak selection approximation for the evolution of altruism against strong selection.

The weak selection approximation of population genetics has made possible the analysis of social evolution under a considerable variety of biological scenarios. Despite its extensive usage, the accuracy of weak selection in predicting the emergence of altruism under limited dispersal when selection intensity increases remains unclear. Here, we derive the condition for the spread of an altruistic mutant in the infinite island model of dispersal under a Moran reproductive process and arbitrary strength of selection. The simplicity of the model allows us to compare weak and strong selection regimes analytically. Our results demonstrate that the weak selection approximation is robust to moderate increases in selection intensity and therefore provides a good approximation to understand the invasion of altruism in spatially structured population. In particular, we find that the weak selection approximation is excellent even if selection is very strong, when either migration is much stronger than selection or when patches are large. Importantly, we emphasize that the weak selection approximation provides the ideal condition for the invasion of altruism, and increasing selection intensity will impede the emergence of altruism. We discuss that this should also hold for more complicated life cycles and for culturally transmitted altruism. Using the weak selection approximation is therefore unlikely to miss out on any demographic scenario that lead to the evolution of altruism under limited dispersal.