Identification of stock units of Mullus barbatus (Linnaeus, 1758) in the Italian seas: development of genomic 2b-RAD markers and biological variation analysis.

This thesis is settled within the STOCKMAPPING project, which represents one of the studies that were developed in the framework of RITMARE Flagship project. The main goals of STOCKMAPPING were the creation of a genomic mapping for stocks of demersal target species and the assembling of a database of population genomic, in order to identify stocks and stocks boundaries. The thesis focuses on three main objectives representing the core for the initial assessment of the methodologies and structure that would be applied to the entire STOCKMAPPING project: individuation of an analytical design to identify and locate stocks and stocks boundaries of Mullus barbatus, application of a multidisciplinary approach to validate biological methods and an initial assessment and improvement for the genotyping by sequencing technique utilized (2b-RAD). The first step is the individuation of an analytical design that has to take in to account the biological characteristics of red mullet and being representative for STOCKMAPPING commitments. In this framework a reduction and selection steps was needed due to budget reduction. Sampling areas were ranked according the individuation of four priorities. To guarantee a multidisciplinary approach the biological data associated to the collected samples were used to investigate differences between sampling areas and GSAs. Genomic techniques were applied to red mullet for the first time so an initial assessment of molecular protocols for DNA extraction and 2b-RAD processing were needed. At the end 192 good quality DNAs have been extracted and eight samples have been processed with 2b-RAD. Utilizing the software Stacks for sequences analyses a great number of SNPs markers among the eight samples have been identified. Several tests have been performed changing the main parameter of the Stacks pipeline in order to identify the most explicative and functional sets of parameters.

Hundertneunundzwanzig Xe NMR for investigations of complex and dynamic systems

In this work, the remarkable versatility and usefulness of applications of Xe-129 NMR experiments is further extended. The application of Xe-129 NMR spectroscopy to very different system is studied, including dynamic and static, solid and liquid, porous and non-porous systems. Using the large non-equilibrium polarization created by hyperpolarization of Xe-129, time-resolved NMR measurements can be used for the online-monitoring of dynamic systems. In the first part of this work, several improvements for medical applications of hyperpolarized Xe-129 are achieved and their feasibility shown experimentally. A large gain in speed and reproducibility of the accumulation process of Xe-129 as ice and an enhancement of the usable polarization in any experiment requiring prior accumulation are achieved. An enhancement of the longitudinal relaxation time of Xe-129 is realized by admixture of a buffer gas during the storage of hyperpolarized Xe-129. Pursuing the efforts of simplifying the accumulation process and enhancing the storage time of hyperpolarized Xe-129 will allow for a wider use of the hyperpolarized gas in (medical) MRI experiments. Concerning the use of hyperpolarized Xe-129 in MRI, the influence of the diffusion coefficient of the gas on parameters of the image contrast is experimentally demonstrated here by admixture of a buffer gas and thus changing the diffusion coefficient. In the second part of this work, a polymer system with unique features is probed by Xe-129 NMR spectroscopy, proving the method to be a valuable tool for the characterization of the anisotropic properties of semicrystalline, syndiotactic polystyrene films. The polymer films contain hollow cavities or channels with sizes in the sub-nanometer range, allowing for adsorption of Xe-129 and subsequent NMR measurements. Despite the use of a ’real-world’ system, the transfer of the anisotropic properties from the material to adsorbed Xe-129 atoms is shown, which was previously only known for fully crystalline materials. The anisotropic behavior towards atomar guests inside the polymer films is proven here for the first time for one of the phases. For the polymer phase containing nanochannels, the dominance of interactions between Xe-129 atoms in the channels compared to interactions between Xe atoms and the channel walls are proven by measurements of a powder sample of the polymer material and experiments including the rotation of the films in the external magnetic field as well as temperature-dependent measurements. The characterization of ’real-world’ systems showing very high degrees of anisotropy by Xe-129 are deemed to be very valuable in future applications. In the last part of this work, a new method for the online monitoring of chemical reactions has been proposed and its feasibility and validity are experimentally proven. The chemical shift dependence of dissolved Xe-129 on the composition of a reaction mixture is used for the online monitoring of free-radical miniemulsion polymerization reactions. Xe-129 NMR spectroscopy provides an excellent method for the online monitoring of polymerization reactions, due to the simplicity of the Xe-129 NMR spectra and the simple relationship between the Xe-129 chemical shift and the reaction conversion. The results of the time-resolved Xe-129 NMR measurements are compared to those from calorimetric measurements, showing a good qualitative agreement. The applicability of the new method to reactions other than polymerization reactions is investigated by the online monitoring of an enzymatic reaction in a miniemulsion. The successful combination of the large sensitivity of Xe-129, the NMR signal enhancements due to hyperpolarization, and the solubility of Xe-129 gives access to the large new field of investigations of chemical reaction kinetics in dynamic and complex systems like miniemulsions.

Boundary condition identification of a bridge pier using experimental data and FEM modal updating

Over the past twenty years, new technologies have required an increasing use of mathematical models in order to understand better the structural behavior: finite element method is the one mostly used. However, the reliability of this method applied to different situations has to be tried each time. Since it is not possible to completely model the reality, different hypothesis must be done: these are the main problems of FE modeling. The following work deals with this problem and tries to figure out a way to identify some of the unknown main parameters of a structure. This main research focuses on a particular path of study and development, but the same concepts can be applied to other objects of research. The main purpose of this work is the identification of unknown boundary conditions of a bridge pier using the data acquired experimentally with field tests and a FEM modal updating process. This work doesn’t want to be new, neither innovative. A lot of work has been done during the past years on this main problem and many solutions have been shown and published. This thesis just want to rework some of the main aspects of the structural optimization process, using a real structure as fitting model.

Development of Dynamic Constraint Models for a Model Based Transient Calibration Process

Model based calibration has gained popularity in recent years as a method to optimize increasingly complex engine systems. However virtually all model based techniques are applied to steady state calibration. Transient calibration is by and large an emerging technology. An important piece of any transient calibration process is the ability to constrain the optimizer to treat the problem as a dynamic one and not as a quasi-static process. The optimized air-handling parameters corresponding to any instant of time must be achievable in a transient sense; this in turn depends on the trajectory of the same parameters over previous time instances. In this work dynamic constraint models have been proposed to translate commanded to actually achieved air-handling parameters. These models enable the optimization to be realistic in a transient sense. The air handling system has been treated as a linear second order system with PD control. Parameters for this second order system have been extracted from real transient data. The model has been shown to be the best choice relative to a list of appropriate candidates such as neural networks and first order models. The selected second order model was used in conjunction with transient emission models to predict emissions over the FTP cycle. It has been shown that emission predictions based on air-handing parameters predicted by the dynamic constraint model do not differ significantly from corresponding emissions based on measured air-handling parameters.

Cardiac transplantation with hearts from donors after circulatory declaration of death: haemodynamic and biochemical parameters at procurement predict recovery following cardioplegic storage in a rat model

OBJECTIVES: Donation after circulatory declaration of death (DCDD) could significantly improve the number of cardiac grafts for transplantation. Graft evaluation is particularly important in the setting of DCDD given that conditions of cardio-circulatory arrest and warm ischaemia differ, leading to variable tissue injury. The aim of this study was to identify, at the time of heart procurement, means to predict contractile recovery following cardioplegic storage and reperfusion using an isolated rat heart model. Identification of reliable approaches to evaluate cardiac grafts is key in the development of protocols for heart transplantation with DCDD. METHODS: Hearts isolated from anaesthetized male Wistar rats (n = 34) were exposed to various perfusion protocols. To simulate DCDD conditions, rats were exsanguinated and maintained at 37°C for 15-25 min (warm ischaemia). Isolated hearts were perfused with modified Krebs-Henseleit buffer for 10 min (unloaded), arrested with cardioplegia, stored for 3 h at 4°C and then reperfused for 120 min (unloaded for 60 min, then loaded for 60 min). Left ventricular (LV) function was assessed using an intraventricular micro-tip pressure catheter. Statistical significance was determined using the non-parametric Spearman rho correlation analysis. RESULTS: After 120 min of reperfusion, recovery of LV work measured as developed pressure (DP)-heart rate (HR) product ranged from 0 to 15 ± 6.1 mmHg beats min(-1) 10(-3) following warm ischaemia of 15-25 min. Several haemodynamic parameters measured during early, unloaded perfusion at the time of heart procurement, including HR and the peak systolic pressure-HR product, correlated significantly with contractile recovery after cardioplegic storage and 120 min of reperfusion (P < 0.001). Coronary flow, oxygen consumption and lactate dehydrogenase release also correlated significantly with contractile recovery following cardioplegic storage and 120 min of reperfusion (P < 0.05). CONCLUSIONS: Haemodynamic and biochemical parameters measured at the time of organ procurement could serve as predictive indicators of contractile recovery. We believe that evaluation of graft suitability is feasible prior to transplantation with DCDD, and may, consequently, increase donor heart availability.

Static and Dynamic Error Correction of a Computer-Aided Mechanical Navigation Linkage for Arthroscopic Hip Surgery

While beneficially decreasing the necessary incision size, arthroscopic hip surgery increases the surgical complexity due to loss of joint visibility. To ease such difficulty, a computer-aided mechanical navigation system was developed to present the location of the surgical tool relative to the patient¿s hip joint. A preliminary study reduced the position error of the tracking linkage with limited static testing trials. In this study, a correction method, including a rotational correction factor and a length correction function, was developed through more in-depth static testing. The developed correction method was then applied to additional static and dynamic testing trials to evaluate its effectiveness. For static testing, the position error decreased from an average of 0.384 inches to 0.153 inches, with an error reduction of 60.5%. Three parameters utilized to quantify error reduction of dynamic testing did not show consistent results. The vertex coordinates achieved 29.4% of error reduction, yet with large variation in the upper vertex. The triangular area error was reduced by 5.37%, however inconsistent among all five dynamic trials. Error of vertex angles increased, indicating a shape torsion using the developed correction method. While the established correction method effectively and consistently reduced position error in static testing, it did not present consistent results in dynamic trials. More dynamic paramters should be explored to quantify error reduction of dynamic testing, and more in-depth dynamic testing methodology should be conducted to further improve the accuracy of the computer-aided nagivation system.

Identification of SNPs in the cystic fibrosis interactome influencing pulmonary progression in cystic fibrosis

There is growing evidence that the great phenotypic variability in patients with cystic fibrosis (CF) not only depends on the genotype, but apart from a combination of environmental and stochastic factors predominantly also on modifier gene effects. It has been proposed that genes interacting with CF transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC) are potential modifiers. Therefore, we assessed the impact of single-nucleotide polymorphisms (SNPs) of several of these interacters on CF disease outcome. SNPs that potentially alter gene function were genotyped in 95 well-characterized p.Phe508del homozygous CF patients. Linear mixed-effect model analysis was used to assess the relationship between sequence variants and the repeated measurements of lung function parameters. In total, we genotyped 72 SNPs in 10 genes. Twenty-five SNPs were used for statistical analysis, where we found strong associations for one SNP in PPP2R4 with the lung clearance index (P ≤ 0.01), the specific effective airway resistance (P ≤ 0.005) and the forced expiratory volume in 1 s (P ≤ 0.005). In addition, we identified one SNP in SNAP23 to be significantly associated with three lung function parameters as well as one SNP in PPP2R1A and three in KRT19 to show a significant influence on one lung function parameter each. Our findings indicate that direct interacters with CFTR, such as SNAP23, PPP2R4 and PPP2R1A, may modify the residual function of p.Phe508del-CFTR while variants in KRT19 may modulate the amount of p.Phe508del-CFTR at the apical membrane and consequently modify CF disease.

[Feline hemoplasmas in Switzerland: identification of a novel species, diagnosis, prevalence, and clinical importance]

Two feline hemotropic mycoplasma spp. (aka hemoplasma) have previously been recognized. We recently discovered a third novel species in a cat with hemolytic anemia, designated 'Candidatus Mycoplasma turicensis', which is closely related to rodent haemoplasmas. This novel species induced anemia after experimental transmission to two SPF cats. Three quantitative real-time PCR assays were newly designed and applied to an epidemiological study surveying the Swiss pet cat population. Blood samples from 713 healthy and ill cats were analyzed. Up to 104 parameters per cat (detailed questionnaire, case history, laboratory parameters and retroviral infections) were evaluated. 'Candidatus Mycoplasma haemominutum' infection was more prevalent (8.5%) than Mycoplasma haemofelis (0.5%) and 'Candidatus Mycoplasma turicensis' (1%). Hemoplasma infections were associated with male gender, outdoor access, and old age, but not with disease or anemia. Infections were more frequently found in the South and West of Switzerland. Several hemoplasma infected cats, some acutely infected, others co-infected with FIV or FeLV, showed hemolytic anemia indicating that additional factors might play a role in the pathogenesis of the disease.

Detection of regional blood perfusion changes in epileptic seizures with dynamic brain perfusion CT--a pilot study

BACKGROUND AND PURPOSE: Perfusion CT (P-CT) is used for acute stroke management, not, however, for evaluating epilepsy. To test the hypothesis that P-CT may identify patients with increased regional cerebral blood flow during subtle status epilepticus (SSE), we compared P-CT in SSE to different postictal conditions. METHODS: Fifteen patients (mean age 47 years, range 21-74) underwent P-CT immediately after evaluation in our emergency room. Asymmetry indices between affected and unaffected hemispheres were calculated for regional cerebral blood volume (rCBV), regional cerebral blood flow (rCBF), and mean transit time (MTT). Regional perfusion changes were compared to EEG findings. RESULTS: Three patients in subtle status epilepticus (group 1) had increased regional perfusion with electro-clinical correlate. Six patients showed postictal slowing on EEG corresponding to an area of regional hypoperfusion (group 2). CT and EEG were normal in six patients with a first epileptic seizure (group 3). Cluster analysis of asymmetry indices separated SSE from the other two groups in all three parameters, while rCBF helped to distinguish between chronic focal epilepsies and single events. CONCLUSION: Preliminary results indicate that P-CT may help to identify patients with SSE during emergency workup. This technique provides important information to neurologists or emergency physicians in the difficult clinical differential diagnosis of altered mental status due to subtle status epilepticus.

Effects of an alveolar recruitment maneuver on cardiovascular and respiratory parameters during total intravenous anesthesia in ponies

OBJECTIVE: To evaluate pulmonary and cardiovascular effects of a recruitment maneuver (RM) combined with positive end-expiratory pressure (PEEP) during total intravenous anesthesia in ponies. ANIMALS: 6 healthy adult Shetland ponies. PROCEDURE: After premedication with detomidine (10 microg/kg, IV), anesthesia was induced with climazolam (0.06 mg/kg, IV) and ketamine (2.2 mg/kg, IV) and maintained with a constant rate infusion of detomidine (0.024 mg/kg/h), climazolam (0.036 mg/kg/h), and ketamine (2.4 mg/kg/h). The RM was preceded by an incremental PEEP titration and followed by a decremental PEEP titration, both at a constant airway pressure difference (deltaP) of 20 cm H2O. The RM consisted of a stepwise increase in deltaP by 25, 30, and 35 cm H2O obtained by increasing peak inspiratory pressure (PIP) to 45, 50, and 55 cm H2O, while maintaining PEEP at 20 cm H2O. Hemodynamic and pulmonary variables were analyzed at every step of the PEEP titration-RM. RESULTS: During the PEEP titration-RM, there was a significant increase in PaO 2 (+12%), dynamic compliance (+ 62%), and heart rate (+17%) and a decrease in shunt (-19%) and mean arterial blood pressure (-21%) was recorded. Cardiac output remained stable. CONCLUSIONS AND CLINICAL RELEVANCE: Although baseline oxygenation was high, Pa(O2) and dynamic compliance further increased during the RM. Despite the use of high PIP and PEEP and a high tidal volume, limited cardiovascular compromise was detected. A PEEP titration-RM may be used to improve oxygenation in anesthetized ponies. During stable hemodynamic conditions, PEEP titration-RM can be performed with acceptable adverse cardiovascular effects.

Predicting asthma control and exacerbations: chronic asthma as a complex dynamic model

PURPOSE OF REVIEW: Predicting asthma episodes is notoriously difficult but has potentially significant consequences for the individual, as well as for healthcare services. The purpose of this review is to describe recent insights into the prediction of acute asthma episodes in relation to classical clinical, functional or inflammatory variables, as well as present a new concept for evaluating asthma as a dynamically regulated homeokinetic system. RECENT FINDINGS: Risk prediction for asthma episodes or relapse has been attempted using clinical scoring systems, considerations of environmental factors and lung function, as well as inflammatory and immunological markers in induced sputum or exhaled air, and these are summarized here. We have recently proposed that newer mathematical methods derived from statistical physics may be used to understand the complexity of asthma as a homeokinetic, dynamic system consisting of a network comprising multiple components, and also to assess the risk for future asthma episodes based on fluctuation analysis of long time series of lung function. SUMMARY: Apart from the classical analysis of risk factor and functional parameters, this new approach may be used to assess asthma control and treatment effects in the individual as well as in future research trials.

Theoretical considerations to the verification of dynamic multileaf collimated fields with a SLIC-type portal image detector

The verification possibilities of dynamically collimated treatment beams with a scanning liquid ionization chamber electronic portal image device (SLIC-EPID) are investigated. The ion concentration in the liquid of a SLIC-EPID and therefore the read-out signal is determined by two parameters of a differential equation describing the creation and recombination of the ions. Due to the form of this equation, the portal image detector describes a nonlinear dynamic system with memory. In this work, the parameters of the differential equation were experimentally determined for the particular chamber in use and for an incident open 6 MV photon beam. The mathematical description of the ion concentration was then used to predict portal images of intensity-modulated photon beams produced by a dynamic delivery technique, the sliding window approach. Due to the nature of the differential equation, a mathematical condition for 'reliable leaf motion verification' in the sliding window technique can be formulated. It is shown that the time constants for both formation and decay of the equilibrium concentration in the chamber is in the order of seconds. In order to guarantee reliable leaf motion verification, these time constants impose a constraint on the rapidity of the image-read out for a given maximum leaf speed. For a leaf speed of 2 cm s(-1), a minimum image acquisition frequency of about 2 Hz is required. Current SLIC-EPID systems are usually too slow since they need about a second to acquire a portal image. However, if the condition is fulfilled, the memory property of the system can be used to reconstruct the leaf motion. It is shown that a simple edge detecting algorithm can be employed to determine the leaf positions. The method is also very robust against image noise.

Application and assessment of a robust elastic motion correction algorithm to dynamic MRI

The purpose of this study was to assess the performance of a new motion correction algorithm. Twenty-five dynamic MR mammography (MRM) data sets and 25 contrast-enhanced three-dimensional peripheral MR angiographic (MRA) data sets which were affected by patient motion of varying severeness were selected retrospectively from routine examinations. Anonymized data were registered by a new experimental elastic motion correction algorithm. The algorithm works by computing a similarity measure for the two volumes that takes into account expected signal changes due to the presence of a contrast agent while penalizing other signal changes caused by patient motion. A conjugate gradient method is used to find the best possible set of motion parameters that maximizes the similarity measures across the entire volume. Images before and after correction were visually evaluated and scored by experienced radiologists with respect to reduction of motion, improvement of image quality, disappearance of existing lesions or creation of artifactual lesions. It was found that the correction improves image quality (76% for MRM and 96% for MRA) and diagnosability (60% for MRM and 96% for MRA).

Shear band evolution in Zr/Hf-based bulk metallic glasses under static and dynamic indentations

Bulk metallic glasses (BMGs) exhibit superior mechanical properties as compared with other conventional materials and have been proposed for numerous engineering and technological applications. Zr/Hf-based BMGs or tungsten reinforced BMG composites are considered as a potential replacement for depleted uranium armor-piercing projectiles because of their ability to form localized shear bands during impact, which has been known to be the dominant plastic deformation mechanism in BMGs. However, in conventional tensile, compressive and bending tests, limited ductility has been observed because of fracture initiation immediately following the shear band formation. To fully investigate shear band characteristics, indentation tests that can confine the deformation in a limited region have been pursued. In this thesis, a detailed investigation of thermal stability and mechanical deformation behavior of Zr/Hf-based BMGs is conducted. First, systematic studies had been implemented to understand the influence of relative compositions of Zr and Hf on thermal stability and mechanical property evolution. Second, shear band evolution under indentations were investigated experimentally and theoretically. Three kinds of indentation studies were conducted on BMGs in the current study. (a) Nano-indentation to determine the mechanical properties as a function of Hf/Zr content. (b) Static Vickers indentation on bonded split specimens to investigate the shear band evolution characteristics beneath the indention. (c) Dynamic Vickers indentation on bonded split specimens to investigate the influence of strain rate. It was found in the present work that gradually replacing Zr by Hf remarkably increases the density and improves the mechanical properties. However, a slight decrease in glass forming ability with increasing Hf content has also been identified through thermodynamic analysis although all the materials in the current study were still found to be amorphous. Many indentation studies have revealed only a few shear bands surrounding the indent on the top surface of the specimen. This small number of shear bands cannot account for the large plastic deformation beneath the indentations. Therefore, a bonded interface technique has been used to observe the slip-steps due to shear band evolution. Vickers indentations were performed along the interface of the bonded split specimen at increasing loads. At small indentation loads, the plastic deformation was primarily accommodated by semi-circular primary shear bands surrounding the indentation. At higher loads, secondary and tertiary shear bands were formed inside this plastic zone. A modified expanding cavity model was then used to predict the plastic zone size characterized by the shear bands and to identify the stress components responsible for the evolution of the various types of shear bands. The applicability of various hardness—yield-strength ( H −σγ ) relationships currently available in the literature for bulk metallic glasses (BMGs) is also investigated. Experimental data generated on ZrHf-based BMGs in the current study and those available elsewhere on other BMG compositions were used to validate the models. A modified expanding-cavity model, employed in earlier work, was extended to propose a new H −σγ relationship. Unlike previous models, the proposed model takes into account not only the indenter geometry and the material properties, but also the pressure sensitivity index of the BMGs. The influence of various model parameters is systematically analyzed. It is shown that there is a good correlation between the model predictions and the experimental data for a wide range of BMG compositions. Under dynamic Vickers indentation, a decrease in indentation hardness at high loading rate was observed compared to static indentation hardness. It was observed that at equivalent loads, dynamic indentations produced more severe deformation features on the loading surface than static indentations. Different from static indentation, two sets of widely spaced semi-circular shear bands with two different curvatures were observed. The observed shear band pattern and the strain rate softening in indentation hardness were rationalized based on the variations in the normal stress on the slip plane, the strain rate of shear and the temperature rise associated with the indentation deformation. Finally, a coupled thermo-mechanical model is proposed that utilizes a momentum diffusion mechanism for the growth and evolution of the final spacing of shear bands. The influence of strain rate, confinement pressure and critical shear displacement on the shear band spacing, temperature rise within the shear band, and the associated variation in flow stress have been captured and analyzed. Consistent with the known pressure sensitive behavior of BMGs, the current model clearly captures the influence of the normal stress in the formation of shear bands. The normal stress not only reduces the time to reach critical shear displacement but also causes a significant temperature rise during the shear band formation. Based on this observation, the variation of shear band spacing in a typical dynamic indentation test has been rationalized. The temperature rise within a shear band can be in excess of 2000K at high strain rate and high confinement pressure conditions. The associated drop in viscosity and flow stress may explain the observed decrease in fracture strength and indentation hardness. The above investigations provide valuable insight into the deformation behavior of BMGs under static and dynamic loading conditions. The shear band patterns observed in the above indentation studies can be helpful to understand and model the deformation features under complex loading scenarios such as the interaction of a penetrator with armor. Future work encompasses (1) extending and modifying the coupled thermo-mechanical model to account for the temperature rise in quasistatic deformation; and (2) expanding this model to account for the microstructural variation-crystallization and free volume migration associated with the deformation.

Dynamic modeling, simulation and control design of a parafoil-payload system for ship launched aerial delivery system (SLADS)

The objective of this research was to develop a high-fidelity dynamic model of a parafoilpayload system with respect to its application for the Ship Launched Aerial Delivery System (SLADS). SLADS is a concept in which cargo can be transfered from ship to shore using a parafoil-payload system. It is accomplished in two phases: An initial towing phase when the glider follows the towing vessel in a passive lift mode and an autonomous gliding phase when the system is guided to the desired point. While many previous researchers have analyzed the parafoil-payload system when it is released from another airborne vehicle, limited work has been done in the area of towing up the system from ground or sea. One of the main contributions of this research was the development of a nonlinear dynamic model of a towed parafoil-payload system. After performing an extensive literature review of the existing methods of modeling a parafoil-payload system, a five degree-of-freedom model was developed. The inertial and geometric properties of the system were investigated to predict accurate results in the simulation environment. Since extensive research has been done in determining the aerodynamic characteristics of a paraglider, an existing aerodynamic model was chosen to incorporate the effects of air flow around the flexible paraglider wing. During the towing phase, it is essential that the parafoil-payload system follow the line of the towing vessel path to prevent an unstable flight condition called ‘lockout’. A detailed study of the causes of lockout, its mathematical representation and the flight conditions and the parameters related to lockout, constitute another contribution of this work. A linearized model of the parafoil-payload system was developed and used to analyze the stability of the system about equilibrium conditions. The relationship between the control surface inputs and the stability was investigated. In addition to stability of flight, one more important objective of SLADS is to tow up the parafoil-payload system as fast as possible. The tension in the tow cable is directly proportional to the rate of ascent of the parafoil-payload system. Lockout instability is more favorable when tow tensions are large. Thus there is a tradeoff between susceptibility to lockout and rapid deployment. Control strategies were also developed for optimal tow up and to maintain stability in the event of disturbances.