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.

Statistical modeling of bowing control applied to violin sound synthesis

Excitation-continuous music instrument control patterns are often not explicitly represented in current sound synthesis techniques when applied to automatic performance. Both physical model-based and sample-based synthesis paradigmswould benefit from a flexible and accurate instrument control model, enabling the improvement of naturalness and realism. Wepresent a framework for modeling bowing control parameters inviolin performance. Nearly non-intrusive sensing techniques allow for accurate acquisition of relevant timbre-related bowing control parameter signals.We model the temporal contour of bow velocity, bow pressing force, and bow-bridge distance as sequences of short Bézier cubic curve segments. Considering different articulations, dynamics, and performance contexts, a number of note classes are defined. Contours of bowing parameters in a performance database are analyzed at note-level by following a predefined grammar that dictates characteristics of curve segment sequences for each of the classes in consideration. As a result, contour analysis of bowing parameters of each note yields an optimal representation vector that is sufficient for reconstructing original contours with significant fidelity. From the resulting representation vectors, we construct a statistical model based on Gaussian mixtures suitable for both the analysis and synthesis of bowing parameter contours. By using the estimated models, synthetic contours can be generated through a bow planning algorithm able to reproduce possible constraints caused by the finite length of the bow. Rendered contours are successfully used in two preliminary synthesis frameworks: digital waveguide-based bowed stringphysical modeling and sample-based spectral-domain synthesis.

Interobserver and intraobserver variability of apparent diffusion coefficient measurements in treated malignant hepatic lesions: measurements in the whole lesion versus in the area with the most restricted diffusion

Purpose: The increase of apparent diffusion coefficient (ADC) in treated hepatic malignancies compared to pre-therapeutic values has been interpreted as treatment success; however, the variability of ADC measurements remains unknown. Furthermore, ADC has been usually measured in the whole lesion, while measurements should be probably centered on the area with the most restricted diffusion (MRDA) as it represents potential tumoral residue. Our objective was to compare the inter/intraobserver variability of ADC measurements in the whole lesion and in MRDA. Material and methods: Forty patients previously treated with chemoembolization or radiofrequency were evaluated (20 on 1.5T and 20 on 3.0T). After consensual agreement on the best ADC image, two readers measured the ADC values using separate regions of interest that included the whole lesion and the whole MRDA without exceeding their borders. The same measurements were repeated two weeks later. Spearman test and the Bland-Altman method were used. Results: Interobserver correlation in ADC measurements in the whole lesion and MRDA was as follows: 0.962 and 0.884. Intraobserver correlation was, respectively, 0.992 and 0.979. Interobserver limits of variability (mm2/sec*10-3) were between -0.25/+0.28 in the whole lesion and between -0.51/+0.46 in MRDA. Intraobserver limits of variability were, respectively: -0.25/+0.24 and -0.43/+0.47. Conclusion: We observed a good inter/intraobserver correlation in ADC measurements. Nevertheless, a limited variability does exist, and it should be considered when interpreting ADC values of hepatic malignancies.

Pavement Management Performance Modeling: Evaluating the Existing PCI Equations, RB14-013 2014

The work described in this report documents the activities performed for the evaluation, development, and enhancement of the Iowa Department of Transportation (DOT) pavement condition information as part of their pavement management system operation. The study covers all of the Iowa DOT’s interstate and primary National Highway System (NHS) and non-NHS system. A new pavement condition rating system that provides a consistent, unified approach in rating pavements in Iowa is being proposed. The proposed 100-scale system is based on five individual indices derived from specific distress data and pavement properties, and an overall pavement condition index, PCI-2, that combines individual indices using weighting factors. The different indices cover cracking, ride, rutting, faulting, and friction. The Cracking Index is formed by combining cracking data (transverse, longitudinal, wheel-path, and alligator cracking indices). Ride, rutting, and faulting indices utilize the International Roughness Index (IRI), rut depth, and fault height, respectively.

Transcranial Doppler Pulsatility Index: What it is and What it Isn't.

BACKGROUND: Transcranial Doppler (TCD) pulsatility index (PI) has traditionally been interpreted as a descriptor of distal cerebrovascular resistance (CVR). We sought to evaluate the relationship between PI and CVR in situations, where CVR increases (mild hypocapnia) and decreases (plateau waves of intracranial pressure-ICP). METHODS: Recordings from patients with head-injury undergoing monitoring of arterial blood pressure (ABP), ICP, cerebral perfusion pressure (CPP), and TCD assessed cerebral blood flow velocities (FV) were analyzed. The Gosling pulsatility index (PI) was compared between baseline and ICP plateau waves (n = 20 patients) or short term (30-60 min) hypocapnia (n = 31). In addition, a modeling study was conducted with the "spectral" PI (calculated using fundamental harmonic of FV) resulting in a theoretical formula expressing the dependence of PI on balance of cerebrovascular impedances. RESULTS: PI increased significantly (p < 0.001) while CVR decreased (p < 0.001) during plateau waves. During hypocapnia PI and CVR increased (p < 0.001). The modeling formula explained more than 65% of the variability of Gosling PI and 90% of the variability of the "spectral" PI (R = 0.81 and R = 0.95, respectively). CONCLUSION: TCD pulsatility index can be easily and quickly assessed but is usually misinterpreted as a descriptor of CVR. The mathematical model presents a complex relationship between PI and multiple haemodynamic variables.

Automation of DEM Cutting for Hydrologic/Hydraulic Modeling, TR-631, 2015

Hydrologic analysis is a critical part of transportation design because it helps ensure that hydraulic structures are able to accommodate the flow regimes they are likely to see. This analysis is currently conducted using computer simulations of water flow patterns, and continuing developments in elevation survey techniques result in higher and higher resolution surveys. Current survey techniques now resolve many natural and anthropogenic features that were not practical to map and, thus, require new methods for dealing with depressions and flow discontinuities. A method for depressional analysis is proposed that uses the fact that most anthropogenically constructed embankments are roughly more symmetrical with greater slopes than natural depressions. An enforcement method for draining depressions is then analyzed on those depressions that should be drained. This procedure has been evaluated on a small watershed in central Iowa, Walnut Creek of the South Skunk River, HUC12 # 070801050901, and was found to accurately identify 88 of 92 drained depressions and place enforcements within two pixels, although the method often tries to drain prairie pothole depressions that are bisected by anthropogenic features.

Isozymic variability in a Brazilian collection of annatto (Bixa orellana L.)

The objectives of this work were to optimize the isozyme electrophoresis technique for Bixa orellana, and use isozyme markers for a preliminary survey on the genetic variability in Brazilian annatto germplasm accessions. Collection consisted of seed samples from sixty open pollinated trees, representing two Northern and four Southern geographic provenances. The extraction, electrophoresis, and interpretation of annatto isozymes are described. Three out of the twenty-one identified isozyme loci were polymorphic in the collection. The percentage of polymorphic loci (P = 21.05) and the expected heterozygosity in annatto (H T = 0.064) were low, compared to other tropical woody species. A UPGMA phenogram, constructed with Nei's genetic distances, clearly separated the germplasm provenant from North and Central Brazil. Variability was significantly higher among the accessions from Maranhão. A sharp genetic differentiation was detected between accessions from Maranhão and Pará States, despite their geographical proximity. The distinctive isozyme polymorphism, observed in the accessions from Maranhão, together with reports on local morphological heterogeneity in annatto fruit shape, color, and pubescence, calls for more detailed genetic and taxonomic investigation.

Modeling Merging Behavior at Lane Drops, InTrans Projects 13-466 and 06-277, 2015

In work-zone configurations where lane drops are present, merging of traffic at the taper presents an operational concern. In addition, as flow through the work zone is reduced, the relative traffic safety of the work zone is also reduced. Improving work-zone flow-through merge points depends on the behavior of individual drivers. By better understanding driver behavior, traffic control plans, work zone policies, and countermeasures can be better targeted to reinforce desirable lane closure merging behavior, leading to both improved safety and work-zone capacity. The researchers collected data for two work-zone scenarios that included lane drops with one scenario on the Interstate and the other on an urban arterial roadway. The researchers then modeled and calibrated these scenarios in VISSIM using real-world speeds, travel times, queue lengths, and merging behaviors (percentage of vehicles merging upstream and near the merge point). Once built and calibrated, the researchers modeled strategies for various countermeasures in the two work zones. The models were then used to test and evaluate how various merging strategies affect safety and operations at the merge areas in these two work zones.

Genetic variability of the Brazilian hair sheep breeds

The objectives of this work were to investigate the genetic structure of the Brazilian hair sheep breeds and to determine the origin of the Santa Inês breed. Molecular similarity was determined using Randomly Amplified Polymorphic DNA - Polymerase Chain Reaction markers in 238 individuals from five naturalized sheep breeds: Santa Inês (48 animals), Rabo Largo (48), Somali (48), Morada Nova (48) and Bergamasca (46), collected in Goiás, Sergipe, Bahia, and Ceará States as well as in the Federal District. Fifty-four loci were selected from 19 primers, after a pilot test using 140 primers. Qualitative analyses indicate diagnostic markers for all breeds. All breeds were significantly different from each other. Interbreed differences were explained by 14.92% of the total variation. Santa Inês clustered with Bergamasca (97% bootstrap) and with Rabo Largo, composing the third member of the group (81% bootstrap) while Morada Nova and Somali breeds clustered separately. Each breed should be considered as a separate management and conservation unit, and special care should be taken with Rabo Largo, Morada Nova and Somali breeds, represented by small herds in Brazil.

Statistical modeling of violin bowing parameter contours

We present a framework for modeling right-hand gestures in bowed-string instrument playing, applied to violin. Nearly non-intrusive sensing techniques allow for accurate acquisition of relevant timbre-related bowing gesture parameter cues. We model the temporal contour of bow transversal velocity, bow pressing force, and bow-bridge distance as sequences of short segments, in particular B´ezier cubic curve segments. Considering different articulations, dynamics, andcontexts, a number of note classes is defined. Gesture parameter contours of a performance database are analyzed at note-level by following a predefined grammar that dictatescharacteristics of curve segment sequences for each of the classes into consideration. Based on dynamic programming, gesture parameter contour analysis provides an optimal curve parameter vector for each note. The informationpresent in such parameter vector is enough for reconstructing original gesture parameter contours with significant fidelity. From the resulting representation vectors, weconstruct a statistical model based on Gaussian mixtures, suitable for both analysis and synthesis of bowing gesture parameter contours. We show the potential of the modelby synthesizing bowing gesture parameter contours from an annotated input score. Finally, we point out promising applicationsand developments.

Qualitative modeling identifies IL-11 as a novel regulator in maintaining self-renewal in human pluripotent stem cells.

Pluripotency in human embryonic stem cells (hESCs) and induced pluripotent stem cells (iPSCs) is regulated by three transcription factors-OCT3/4, SOX2, and NANOG. To fully exploit the therapeutic potential of these cells it is essential to have a good mechanistic understanding of the maintenance of self-renewal and pluripotency. In this study, we demonstrate a powerful systems biology approach in which we first expand literature-based network encompassing the core regulators of pluripotency by assessing the behavior of genes targeted by perturbation experiments. We focused our attention on highly regulated genes encoding cell surface and secreted proteins as these can be more easily manipulated by the use of inhibitors or recombinant proteins. Qualitative modeling based on combining boolean networks and in silico perturbation experiments were employed to identify novel pluripotency-regulating genes. We validated Interleukin-11 (IL-11) and demonstrate that this cytokine is a novel pluripotency-associated factor capable of supporting self-renewal in the absence of exogenously added bFGF in culture. To date, the various protocols for hESCs maintenance require supplementation with bFGF to activate the Activin/Nodal branch of the TGFβ signaling pathway. Additional evidence supporting our findings is that IL-11 belongs to the same protein family as LIF, which is known to be necessary for maintaining pluripotency in mouse but not in human ESCs. These cytokines operate through the same gp130 receptor which interacts with Janus kinases. Our finding might explain why mESCs are in a more naïve cell state compared to hESCs and how to convert primed hESCs back to the naïve state. Taken together, our integrative modeling approach has identified novel genes as putative candidates to be incorporated into the expansion of the current gene regulatory network responsible for inducing and maintaining pluripotency.

Pattern dystrophy with high intrafamilial variability associated with Y141C mutation in the peripherin/RDS gene and successful treatment of subfoveal CNV related to multifocal pattern type with anti-VEGF (ranibizumab) intravitreal injections.

OBJECTIVE: To identify disease causing mutation in three generations of a Swiss family with pattern dystrophy and high intrafamilial variability of phenotype. To assess the effect of intravitreal ranibizumab injections in the treatment of subfoveal choroidal neovascularization associated with pattern dystrophy in one patient. METHODS: Affected family members were ascertained for phenotypic and genotypic characterization. Ophthalmic evaluations included fundus photography, autofluorescence imaging, optical coherence tomography, and International Society for Clinical Electrophysiology of Vision standard full-field electroretinography. When possible family members had genetic testing. The proband presented with choroidal neovascularization and had intravitreal injections as needed according to visual acuity and optical coherence tomography. RESULTS: Proband had a multifocal type pattern dystrophy, and his choroidal neovascularization regressed after four intravitreal injections. The vision improved from 0.8 to 1.0, and optical coherence tomography showed complete anatomical restoration. A butterfly-shaped pattern was observed in her cousin, whereas a fundus pulverulentus pattern was seen in a second cousin. Aunt had a multifocal atrophic appearance, simulating geographic atrophy in age-related macular degeneration. The Y141C mutation was identified in the peripherin/RDS gene and segregated with disease in the family. CONCLUSION: This is the first report of marked intrafamilial variation of pattern dystrophy because of peripherin/RDS Y141C mutation. Intravitreal ranibizumab injections might be a valuable treatment for associated subfoveal choroidal neovascularization.

Hierarchical modeling gave plausible estimates of associations between metabolic syndrome and components of antiretroviral therapy.

OBJECTIVE: Hierarchical modeling has been proposed as a solution to the multiple exposure problem. We estimate associations between metabolic syndrome and different components of antiretroviral therapy using both conventional and hierarchical models. STUDY DESIGN AND SETTING: We use discrete time survival analysis to estimate the association between metabolic syndrome and cumulative exposure to 16 antiretrovirals from four drug classes. We fit a hierarchical model where the drug class provides a prior model of the association between metabolic syndrome and exposure to each antiretroviral. RESULTS: One thousand two hundred and eighteen patients were followed for a median of 27 months, with 242 cases of metabolic syndrome (20%) at a rate of 7.5 cases per 100 patient years. Metabolic syndrome was more likely to develop in patients exposed to stavudine, but was less likely to develop in those exposed to atazanavir. The estimate for exposure to atazanavir increased from hazard ratio of 0.06 per 6 months' use in the conventional model to 0.37 in the hierarchical model (or from 0.57 to 0.81 when using spline-based covariate adjustment). CONCLUSION: These results are consistent with trials that show the disadvantage of stavudine and advantage of atazanavir relative to other drugs in their respective classes. The hierarchical model gave more plausible results than the equivalent conventional model.

Crustal structure and reflectivity of the Swiss Alps from 3-Dimensional seismic modeling. 2. Penninic Nappes

Surface geological mapping, laboratory measurements of rock properties, and seismic reflection data are integrated through three-dimensional seismic modeling to determine the likely cause of upper crustal reflections and to elucidate the deep structure of the Penninic Alps in eastern Switzerland. Results indicate that the principal upper crustal reflections recorded on the south end of Swiss seismic line NFP20-EAST can be explained by the subsurface geometry of stacked basement nappes. In addition, modeling results provide improvements to structural maps based solely on surface trends and suggest the presence of previously unrecognized rock units in the subsurface. Construction of the initial model is based upon extrapolation of plunging surface. structures; velocities and densities are established by laboratory measurements of corresponding rock units. Iterative modification produces a best fit model that refines the definition of the subsurface geometry of major structures. We conclude that most reflections from the upper 20 km can be ascribed to the presence of sedimentary cover rocks (especially carbonates) and ophiolites juxtaposed against crystalline basement nappes. Thus, in this area, reflections appear to be principally due to first-order lithologic contrasts. This study also demonstrates not only the importance of three-dimensional effects (sideswipe) in interpreting seismic data, but also that these effects can be considered quantitatively through three-dimensional modeling.