A multicentre, retrospective case-control study assessing the role of trabecular bone score (TBS) in menopausal Caucasian women with low areal bone mineral density (BMDa): Analysing the odds of vertebral fracture.

INTRODUCTION: The trabecular bone score (TBS) is a new parameter that is determined from grey level analysis of DXA images. It relies on the mean thickness and volume fraction of trabecular bone microarchitecture. This was a preliminary case-control study to evaluate the potential diagnostic value of TBS, both alone and combined with bone mineral density (BMDa), in the assessment of vertebral fracture. METHODS: Out of a subject pool of 441 Caucasian, postmenopausal women between the ages of 50 and 80 years, we identified 42 women with osteoporosis-related vertebral fractures, and compared them with 126 age-matched women without any fractures (1 case: 3 controls). Primary outcomes were BMDa and TBS. Inter-group comparisons were undertaken using Student's t-tests and Wilcoxon signed ranks tests for parametric and non-parametric data, respectively. Odds ratios for vertebral fracture were calculated for each incremental one standard deviation decrease in BMDa and TBS, and areas under the receiver operating curve (AUC) calculated and sensitivity analysis were conducted to compare BMDa alone, TBS alone, and the combination of BMDa and TBS. Subgroup analyses were performed specifically for women with osteopenia, and for women with T-score-defined osteoporosis. RESULTS: Across all subjects (n=42, 126) weight and body mass index were greater and BMDa and TBS both less in women with fractures. The odds of vertebral fracture were 3.20 (95% CI, 2.01-5.08) for each incremental decrease in TBS, 1.95 (1.34-2.84) for BMDa, and 3.62 (2.32-5.65) for BMDa + TBS combined. The AUC was greater for TBS than for BMDa (0.746 vs. 0.662, p=0.011). At iso-specificity (61.9%) or iso-sensitivity (61.9%) for both BMDa and TBS, TBS + BMDa sensitivity or specificity was 19.1% or 16.7% greater than for either BMDa or TBS alone. Among subjects with osteoporosis (n=11, 40) both BMDa (p=0.0008) and TBS (p=0.0001) were lower in subjects with fractures, and both OR and AUC (p=0.013) for BMDa + TBS were greater than for BMDa alone (OR=4.04 [2.35-6.92] vs. 2.43 [1.49-3.95]; AUC=0.835 [0.755-0.897] vs. 0.718 [0.627-0.797], p=0.013). Among subjects with osteopenia, TBS was lower in women with fractures (p=0.0296), but BMDa was not (p=0.75). Similarly, the OR for TBS was statistically greater than 1.00 (2.82, 1.27-6.26), but not for BMDa (1.12, 0.56-2.22), as was the AUC (p=0.035), but there was no statistical difference in specificity (p=0.357) or sensitivity (p=0.678). CONCLUSIONS: The trabecular bone score warrants further study as to whether it has any clinical application in osteoporosis detection and the evaluation of fracture risk.

Kernel-based Mapping of Meteorological Fields in Complex Orography

Due to the advances in sensor networks and remote sensing technologies, the acquisition and storage rates of meteorological and climatological data increases every day and ask for novel and efficient processing algorithms. A fundamental problem of data analysis and modeling is the spatial prediction of meteorological variables in complex orography, which serves among others to extended climatological analyses, for the assimilation of data into numerical weather prediction models, for preparing inputs to hydrological models and for real time monitoring and short-term forecasting of weather.In this thesis, a new framework for spatial estimation is proposed by taking advantage of a class of algorithms emerging from the statistical learning theory. Nonparametric kernel-based methods for nonlinear data classification, regression and target detection, known as support vector machines (SVM), are adapted for mapping of meteorological variables in complex orography.With the advent of high resolution digital elevation models, the field of spatial prediction met new horizons. In fact, by exploiting image processing tools along with physical heuristics, an incredible number of terrain features which account for the topographic conditions at multiple spatial scales can be extracted. Such features are highly relevant for the mapping of meteorological variables because they control a considerable part of the spatial variability of meteorological fields in the complex Alpine orography. For instance, patterns of orographic rainfall, wind speed and cold air pools are known to be correlated with particular terrain forms, e.g. convex/concave surfaces and upwind sides of mountain slopes.Kernel-based methods are employed to learn the nonlinear statistical dependence which links the multidimensional space of geographical and topographic explanatory variables to the variable of interest, that is the wind speed as measured at the weather stations or the occurrence of orographic rainfall patterns as extracted from sequences of radar images. Compared to low dimensional models integrating only the geographical coordinates, the proposed framework opens a way to regionalize meteorological variables which are multidimensional in nature and rarely show spatial auto-correlation in the original space making the use of classical geostatistics tangled.The challenges which are explored during the thesis are manifolds. First, the complexity of models is optimized to impose appropriate smoothness properties and reduce the impact of noisy measurements. Secondly, a multiple kernel extension of SVM is considered to select the multiscale features which explain most of the spatial variability of wind speed. Then, SVM target detection methods are implemented to describe the orographic conditions which cause persistent and stationary rainfall patterns. Finally, the optimal splitting of the data is studied to estimate realistic performances and confidence intervals characterizing the uncertainty of predictions.The resulting maps of average wind speeds find applications within renewable resources assessment and opens a route to decrease the temporal scale of analysis to meet hydrological requirements. Furthermore, the maps depicting the susceptibility to orographic rainfall enhancement can be used to improve current radar-based quantitative precipitation estimation and forecasting systems and to generate stochastic ensembles of precipitation fields conditioned upon the orography.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.

Validation of and comparison between a semidistributed rainfall-runoff hydrological model (PREVAH) and a spatially distributed snow-evolution model (SnowModel) for snow cover prediction in mountain ecosystems

Snow cover is an important control in mountain environments and a shift of the snow-free period triggered by climate warming can strongly impact ecosystem dynamics. Changing snow patterns can have severe effects on alpine plant distribution and diversity. It thus becomes urgent to provide spatially explicit assessments of snow cover changes that can be incorporated into correlative or empirical species distribution models (SDMs). Here, we provide for the first time a with a lower overestimation comparison of two physically based snow distribution models (PREVAH and SnowModel) to produce snow cover maps (SCMs) at a fine spatial resolution in a mountain landscape in Austria. SCMs have been evaluated with SPOT-HRVIR images and predictions of snow water equivalent from the two models with ground measurements. Finally, SCMs of the two models have been compared under a climate warming scenario for the end of the century. The predictive performances of PREVAH and SnowModel were similar when validated with the SPOT images. However, the tendency to overestimate snow cover was slightly lower with SnowModel during the accumulation period, whereas it was lower with PREVAH during the melting period. The rate of true positives during the melting period was two times higher on average with SnowModel with a lower overestimation of snow water equivalent. Our results allow for recommending the use of SnowModel in SDMs because it better captures persisting snow patches at the end of the snow season, which is important when modelling the response of species to long-lasting snow cover and evaluating whether they might survive under climate change.